Research on Optimization Method of Evaporation Duct Prediction Model

: The sea surface roughness parameterization and the universal stability function are key components of the evaporation duct prediction model based on the Monin–Obukhov similarity theory. They determine the model’s performance, which in turn affects the efficiency and accuracy of electromagnetic applications at sea. In this study, we collected layered meteorological and hydrological observation data and preprocessed them to obtain near-surface reference modified refractivity profiles. We then optimized the sea surface roughness parameterization and the universal stability function using particle swarm optimization and simulated annealing algorithms. The results show that the particle swarm optimization algorithm outperforms the simulated annealing algorithm. Compared to the original model, the particle swarm optimization algorithm improved the prediction accuracy of the model by 5.09% under stable conditions and by 9.97% under unstable conditions, demonstrating the feasibility of the proposed method for optimizing the evaporation duct prediction model. Subsequently, we compared the electromagnetic wave propagation path losses under two different evaporation duct heights and modified refractivity profile states, confirming that the modified refractivity profile is more suitable as the accuracy criterion for the evaporation duct prediction model.


Introduction
Due to the influence of oceanic momentum and heat flux, as well as the effects of the interaction between maritime dry and hot air masses, wet and cold air masses, the atmospheric boundary layer over the sea is prone to the formation of atmospheric ducting phenomena, resulting in abnormal refraction of electromagnetic waves [1].Under certain conditions, electromagnetic waves emitted by ship-borne radar and other radio frequency devices can be trapped within the ducting layer, which can achieve "beyond-line-of-sight" propagation with very low transmission losses [2].At the same time, the trapping of electromagnetic waves within the ducting layer also leads to the disappearance of some echo regions, and creates "blind zones" [3].Atmospheric ducting includes elevated duct, surface duct, and evaporation duct.Among them, evaporation duct is a type of surface duct that often occurs at the interface between the sea and the atmosphere.The main reason is that the refractivity can sharply decrease with the height in the near-surface atmosphere due to the evaporation of the sea; it causes the trajectory of electromagnetic wave propagation to bend towards the sea surface and exceeds the curvature of the Earth's surface.Evaporation duct often occurs within a range of approximately 40 m above the sea surface [4].The research shows that the probability of this occurs appears in the South China Sea is often above 80%, and there is permanent evaporation duct in some areas [5].Accurately sensing and predicting atmospheric refractivity is important to scientifically use evaporation duct phenomena.This plays a crucial role in accurately predicting the detection range and blind zones of radar for maritime target detection and achieving beyond-line-of-sight target detection and identification.
Evaporation duct can be detected through direct measurement methods and indirect measurement methods.Direct measurement methods include various tools such as meteorological gradient instruments, microwave refractometers, low-level tethered balloons, or other sounding devices to obtain parameters of the evaporation duct and thereby determine the evaporation duct.This method has high accuracy but is not widely used due to its high cost and rigorous measurement environment.Indirect measurement methods are based on the Monin-Obukhov similarity theory and involving the establishment of a prediction model for the evaporation duct.This method is widely applied due to its real-time capability and simplicity; it is currently the most commonly used method for researching evaporation duct [6].
The output of the evaporation duct prediction model includes the evaporation duct height and the modified refractivity profile.Current research primarily focuses on improving the accuracy of predicting evaporation duct height [2,[7][8][9][10].There are relatively few studies on the modified refractivity profile.Based on the evaporation duct height, it is possible to determine whether electromagnetic waves are trapped within the duct at the emission source and achieving the beyond-line-of-sight propagation.However, to obtain the actual propagation characteristics of electromagnetic waves, one needs the corresponding modified refractivity profile for the given evaporation duct height.Nevertheless, the same evaporation duct height may correspond to multiple modified refractivity profiles.To achieve accurate predictions of the modified refractivity profile and precisely forecast the electromagnetic wave propagation, it is necessary to enhance the model's accuracy in predicting the modified refractivity profile.
The commonly used evaporation duct prediction models mainly include the COARE model [11,12], the BYC model [3], the NPS model [13][14][15], the ECMWF model [16], and the MM5REV model [17].These models use iterative methods to calculate intermediate variables such as Monin-Obukhov length and friction velocity with taking into account atmospheric stability.Although these evaporation duct prediction models are all based on the Monin-Obukhov similarity theory, the accuracy of predicting near-surface refractivity profiles can make a difference in different sea regions and seasons.Tian [18] conducted relevant experiments on the seasonal distribution of evaporation duct in the Arabian Sea, and the results showed that evaporation duct exhibit significant seasonal variations, with lower electromagnetic wave propagation losses in summer than in other seasons.The main reason for these differences lies in the using empirical functions (sea surface roughness parameterization schemes and universal stability functions) in the models [4].Zhang [8] compared the modified refractivity profiles predicted by the BH91, CB05 and G07 models; the results show that the G07 model predictions were closer to the actual measurement and confirm the model's applicability issues.The key to the applicability issues of the models is the empirical function parameter values.Due to the complex nonlinear relationship between the data and the model's empirical function parameters, it is difficult to obtain analytical solutions for these parameters based solely on meteorological and hydrological data.Frederickson [19] modified the universal stability function parameters of the NAVSLaM model under unstable conditions.When using the CASPER-East 2015 dataset, the optimal function parameter values obtained were different from the original function parameters, further confirming that parameter values directly affect the model's prediction accuracy.Therefore, we use global optimization algorithms to optimize both sea surface roughness parameterization schemes and universal stability function parameters in this paper, aiming to improve the accuracy of modified refractivity profiles in the specific sea area.
With the rapid development of technology, the optimization algorithms have been increasingly valued for their convenience.Currently, popular research optimization algorithms include Particle Swarm Optimization (PSO), Genetic Algorithms (GA), Simulated Annealing (SA), and Ant Colony Algorithm (ACA), etc.Among many algorithms, Particle Swarm Optimization algorithm is widely recognized for its advantages, such as having few parameters, fast convergence speed and easy implementation.This algorithm was developed by Kennedy and Eberhart [20], inspired by the social behavior of birds and fish.It is a metaheuristic search algorithm that can effectively solve many complex problems and has significant advantages in solving continuous optimization problems.It is now widely used in areas such as system design, pattern recognition, biological system modeling and signal processing [21].Wang [22] compares the performance of PSO, GA and ACA, the results show that PSO has the advantages of faster convergence speed and higher retrieval accuracy.SA is widely used for large-scale optimization problems, especially when the desired global minimum is hidden among multiple inferior local minima.It was first proposed by Metropolis [23] and later applied to solving combinatorial optimization problems by Kirkpatrick [24], who published a paper in Science.This algorithm is a general probabilistic algorithm that simulates the annealing process of slowly cooling heated metal to reach a stable condition with the minimum energy level [25].Zhao [26] compares the performance of the SA algorithm and GA; the results show that the SA algorithm has faster convergence speed and is more suitable for the field of duct.As PSO and SA algorithm have demonstrated superior performance in the field of evaporation duct, we use these two algorithms to optimize the sea surface roughness parameterization and the universal stability function parameters in this paper.By comparing the performance of the two algorithms, the algorithm with faster convergence speed and better approximation to the optimal solution is selected to improve the prediction accuracy of refractivity profiles in the specific sea area.
The structure of this article is as follows: Section 1 is introduction.Section 2 introduces the theoretical framework and related models.Section 3 presents the optimization algorithms and proposes a model optimization method.Section 4 discusses the collection of measured data, the optimization process, and the optimization results.Section 5 provides a summary of the entire article.Finally, in Section 6, the continuity of universal stability functions and whether the evaporation duct height and modified refractivity profile are suitable as criteria for predicting evaporation duct models are discussed.

Monin Obukhov's Similarity Theory
Due to the small variations in atmospheric horizontal levels over the sea, atmospheric stability states are typically described using stability parameters based on the assumption of atmospheric horizontal uniformity [4]: In Equation ( 1), ζ > 0 denotes that the atmosphere is in the stable condition, ζ < 0 denotes that the atmosphere is in the unstable condition and L is the Monin-Obukhov length.
where κ is the von Kármán constant, κ = 0.4, g is the gravitational acceleration, u * is the friction velocity, T is the atmospheric temperature (K), T * , q * are the scalar parameters of temperature and specific humidity, respectively.
To improve the prediction accuracy of the refractivity profile correction in the specific sea area, in this paper, we separately consider the impact of sea surface roughness parameterization schemes and universal stability functions on the prediction results.Furthermore, optimization of their function parameters is carried out to achieve high-precise prediction of the refractivity profile correction in the specific sea area.

Sea Surface Roughness
The sea surface roughness is an important parameter in the study of the air-sea interaction.However, it is difficult to directly obtain sea surface roughness from the actual environment.Generally, parameterization methods are used to obtain roughness [26].The physical quantities commonly used to describe the roughness of the air-sea interface include momentum roughness z o , scalar roughness z oq and z ot [27].To optimize the coefficients of the sea surface roughness parameterization scheme, we select the following sea surface roughness parameterization schemes in this paper.Among them, there are six types of momentum roughness and one type of scalar roughness.

•
The COARE3.0 sea surface roughness parameterization scheme serves is the foundational parameterization scheme for the COARE model and can be represented as [11] where is the Charnock coefficient, a = 0.011, b = 0.007, c = 0.018, u is wind speed.However, this parameterized scheme does not consider the effect of waves on roughness.• The TY01 sea surface roughness parameterization scheme is a dimensionless roughness parameterization scheme based on a wave slope, in which the sea surface roughness is closely related to wave height [28]: where a = 1200, b = 4.5, h s is the significant wave height, l p is peak wavelength, h s l −1 p is wave slope and v is the kinematic viscosity.• The Oost sea surface roughness parameterization scheme considers both wave age and friction velocity, and its roughness can be expressed as a function of wave age and friction velocity [29]: where 2π , a = 45, b = 4.5, T p is the wave period, c p is the phase speed of dominant waves.
• The COARE 3.5 sea surface roughness parameterization scheme is based on the COARE 3.0 sea surface roughness parameterization scheme which modifies the wind speed dependence of the Charnock parameter; it extends the limited range of the parameterization scheme wind speed to 0~20 m/s [12]: where α c = au − b, a = 0.017, b = 0.005.• The Edson sea surface roughness parameterization scheme subdivides the water depth [30]: where , a = 0.091, b = 2.In deep sea, In shallow sea, c p = gh w , where h w is the depth of the water.If in deep sea, h w = 1, if in shallow sea, h w = 0. • The Porchetta sea surface roughness parameterization scheme finds that the roughness length is related to the angle between the wind direction and the wave direction; the functional form is [31] z o = a cos(bθ)h s u * c p c cos(−dθ) where a = 20, b = 0.45, c = 3.8, d = 0.32.
From the above sea surface roughness parameterization, each equation contains the common term 0.11vu −1 * ; it is the smooth-flow component of the total roughness.In lowwind conditions, the sea surface behaves as a nearly smooth aerodynamic surface; it is irrelevant to roughness [32,33], z o = 0.11vu −1 * .In the above-sea surface roughness parameterization scheme, the scalar roughness representation is the same, and both z oq and z ot can be expressed as where c = 5.5, d = 0.6.

Universal Stability Function
The universal stability function is a bridge that connects both the observed values at a certain height and meteorological parameter profiles near the sea surface.Based on the stable state ζ of the atmospheric stratification, the numerical calculations are performed to correct the meteorological parameter profiles under the neutral atmospheric conditions within the framework of the Monin-Obukhov similarity theory, then obtaining the vertical distribution of physical quantities such as wind speed, temperature, humidity, etc., under the stable or unstable conditions.Currently, the commonly used universal stability functions have different specific mathematical forms and parameter values, which are typically determined based on empirical knowledge [19].To improve the prediction accuracy of the evaporation duct prediction model in the specific sea area, it is necessary to optimize the universal stability function and corresponding parameter values.Therefore, we select eight of the most commonly used universal stability functions and optimize them under the stable and unstable conditions, respectively.
For stable conditions: • BH91 is a summary of the work of Beljaars and Holtslag [13] on surface flux observations and modelling of the MESOGERS-84 experiment in Cabauwin, the Netherlands and the south of France, which offers a new functional form under stable conditions: where w = 5, r = 0.35, t = 10.
• CB05 was studied on the performance of the Monin-Obukhov similarity theory under the stable condition by Cheng and Brutsaert [14], during the Cooperative Atmosphere-Surface Exchange Study-99 (CASES-99) and the following expression was obtained: where a = 6.1, b = 2.5, c = 5.3, d = 1.1.

•
G07 is a dataset collected by Grachev [15] during the SHEBA for a wide range under the stable conditions.It includes strongly stable conditions.New Ψ M and Ψ H functions are derived from this dataset, covering the full range of ζ > 0: where • The COARE 3.0 function can be described as [11]      where dζ = min(50, 0.35ζ), a = 0.667, b = 14.28, c = 8.525.• The COARE 3.5 function is [12] where where a = 6.

Evaporation Duct Prediction Model
In evaporation duct, the accuracy of the calculation results of electromagnetic wave propagation models largely depends on the precise description of the modified refractivity.The modified refractivity is a primary factor that influences the electromagnetic wave propagation characteristics in the atmospheric environment.The modified refractivity can be expressed as where N is the atmospheric refractivity, R is the average radius of earth and z (m) is the altitude.N can be expressed as where T is the atmospheric temperature (K), P is the total atmospheric pressure (hPa), e is the water vapor pressure (hPa).T, P, e can be obtained by the sea surface roughness parameterization scheme and universal stability function in Sections 2.2 and 2.3 [8].Therefore, choosing the appropriate sea surface roughness parameterization scheme and universal stability function can obtain higher accuracy N and thus higher accuracy M.
We select a total of six sea surface roughness parameterization schemes and eight universal stability functions in this paper; there are 48 evaporation duct prediction models after combining them.To achieve high-precision modified refractivity profile predictions in the specific sea area, it is necessary to select the combined models.However, since these functions contain multiple parameters, it is necessary to first optimize the parameters and determine their optimization ranges.Subsequently, particle swarm optimization algorithm and simulated annealing algorithm are used to optimize the selected function parameters for improving the prediction accuracy of the modified refractivity profile in the specific sea area.

The Optimization Algorithm
The optimization problem in this paper can be solved using the following equation, where f optimization (x) is the objective function, and s.t. is the constraint: where i denotes the order of the measurement level, j is the sequence of the sample, I is the total number of measurement levels.x is the measurement of meteorological and hydrological data, such as wind speed and temperature.p is the set of parameters to be optimized.
is the model output.The absolute value of the difference between the five levels of f i,j re f erence (x) and f i,j model (x) is first obtained, and its absolute values are averaged to minimize the rms of all the samples as the objective function.p scale is the corresponding parameter range.
The most widely used optimization algorithms are particle swarm optimization algorithm, simulated annealing algorithm, differential evolution algorithm and its improved variants.Since there are fewer parameters to be optimized in this paper, the classical particle swarm optimization algorithm and simulated annealing algorithm are used to improve the modified refractivity prediction accuracy.
In order to explore the impact of two optimization algorithms on the accuracy of predicting the modified refractivity profile, an evaluation is conducted on the optimization results of both algorithms.The optimization algorithm that approaches the global optimum solution more closely is chosen as the subsequent algorithm for optimizing the modified refractivity profile.

Particle Swarm Optimization Algorithm
Particle swarm algorithm is still a population-based algorithm.It was initially conceived to visually emulate the graceful and unpredictable motions of birds.This algorithm boasts a considerable probability of approximating the global optimal solution, alongside its swift computational speed and superior global search capability compared to conventional optimization algorithms [36].Assuming that the parameters to be optimized are D-dimensional vectors, in the search space, F particles are randomly generated and the position of each particle represents a set of parameters to be optimized, then the position of the f th particle is denoted as The optimal position searched for by the f th particle is P best = p f 1 , p f 2, , . . ., p f D , the optimal position searched for by the population is P gbest = p g1 , p g2 , . . ., p gD .The velocity and position of each particle is updated by the following equation: where k denotes the iteration number.The inertia weight is represented as ω, while c 1 and c 2 are the cognitive and social factors, respectively.r 1 and r 2 are random numbers.The particle count is given by f.The position of the f th particle in a D-dimensional space is described by x f D , with its velocity in the same space represented as v f D .The optimal position discovered by the f th particle in the D-dimensional space is symbolized as p f D .Meanwhile, the swarm's overall best-found position in this space is indicated by p gD .There are a large number of particle swarm optimization algorithms and other improved model of them.Due to the fact that optimization objectives are smaller and the optimization function is simple, we employ the standard PSO algorithm in this paper.In this work, the number of particles is set to 200 and the number of iterations is 100, the inertia weight is 0.8, the cognitive factor is 0.8 and the social factor is 1.5.

Simulated Annealing Algorithm
SA is a simulation of the thermodynamic annealing process using a two-layer loop structure, which essentially is a Monte Carlo algorithm [37].But SA is not a populationbased algorithm.Assuming that the previous state is x(n) and the state changes to x(n + 1) at the next moment, accordingly, the energy of the system changes from E(n) to E(n + 1), and defining the acceptance probability Ps of the system changing from x(n) to x(n + 1) as In this work, the initial temperature is 100, the temperature decay coefficient is 0.95, the number of cycles in the outer layer is 100, and the number of cycles in the inner layer is 200.

Optimization Algorithm for Evaporation Duct Prediction Model
The optimization process of the evaporation duct prediction model, driven by optimization algorithms, comprises four key steps.The first step is data handling.The second step involves model selection followed by the third step where we determine the parame-ters to be optimized and their respective ranges.Lastly, the fourth step encompasses data classification and model optimization.The specific process is shown in Figure 1.

Optimization Algorithm for Evaporation Duct Prediction Model
The optimization process of the evaporation duct prediction model, driven by opti mization algorithms, comprises four key steps.The first step is data handling.The secon step involves model selection followed by the third step where we determine the param eters to be optimized and their respective ranges.Lastly, the fourth step encompasses dat classification and model optimization.The specific process is shown in Figure 1.The main details are described further.First, the collected meteorological and hydro logical data need to be quality controlled, for example, by using the sliding averag method to eliminate the noise problems present in the data.The required data are orga nized into the form of datasets to be passed into the evaporation duct prediction model.
Second, in order to enhance optimization efficiency, we need to choose the combina tion model of the sea surface roughness parameterization scheme and the universal sta bility function under the original parameter conditions, thereby obtaining the combina tion model with the best prediction accuracy in the specific sea area.Figure 2 is mode combination.The main details are described further.First, the collected meteorological and hydrological data need to be quality controlled, for example, by using the sliding average method to eliminate the noise problems present in the data.The required data are organized into the form of datasets to be passed into the evaporation duct prediction model.
Second, in order to enhance optimization efficiency, we need to choose the combination model of the sea surface roughness parameterization scheme and the universal stability function under the original parameter conditions, thereby obtaining the combination model with the best prediction accuracy in the specific sea area.Figure 2 is model combination.
the number of cycles in the outer layer is 100, and the number of cycles in the inner layer is 200.

Optimization Algorithm for Evaporation Duct Prediction Model
The optimization process of the evaporation duct prediction model, driven by optimization algorithms, comprises four key steps.The first step is data handling.The second step involves model selection followed by the third step where we determine the parameters to be optimized and their respective ranges.Lastly, the fourth step encompasses data classification and model optimization.The specific process is shown in Figure 1.The main details are described further.First, the collected meteorological and hydrological data need to be quality controlled, for example, by using the sliding average method to eliminate the noise problems present in the data.The required data are organized into the form of datasets to be passed into the evaporation duct prediction model.
Second, in order to enhance optimization efficiency, we need to choose the combination model of the sea surface roughness parameterization scheme and the universal stability function under the original parameter conditions, thereby obtaining the combination model with the best prediction accuracy in the specific sea area.Figure 2 is model combination.rms is the root mean square deviation of the value between M model and M measurement , which is used to be the fitness of the accuracy of the model prediction results.M model is the modified refractivity profile of the model output, M measurement is the modified refractivity profile by environmentally measured.The smaller the deviation, ∆M rms_selection , the higher the model prediction accuracy: where i denotes the order of the measurement level, I is the total number of measurement levels, j is the sequence of the sample and J signifies the total sample count.For the jth sample at the ith layer, M i,j measurement is the observed modified refractivity and M i,j model is the corresponding predicted value from the model.
Third, after obtaining the optimal combination of models, we need to determine the parameters of the sea surface roughness parameterization scheme (e.g., a, b in TY01) and the parameters of the universal stability function (e.g., a m , d m in Ψ M_stable in the stable condition of G07) and the optimization ranges of the corresponding parameters.If the function comprises multiple coefficients, we select some parameters of this function to ensure an equal number of coefficients to be optimized for each function given in Figure 3. , the higher the model prediction accuracy: , , where i denotes the order of the measurement level, I is the total number of measurement levels, j is the sequence of the sample and J signifies the total sample count.For the jth sample at the ith layer, , i j measurement

M
is the observed modified refractivity and , i j model M is the corresponding predicted value from the model.Third, after obtaining the optimal combination of models, we need to determine the parameters of the sea surface roughness parameterization scheme (e.g., a, b in TY01) and the parameters of the universal stability function (e.g., m a , m d in Ψ _ M stable in the stable condition of G07) and the optimization ranges of the corresponding parameters.If the function comprises multiple coefficients, we select some parameters of this function to ensure an equal number of coefficients to be optimized for each function given in Figure 3.If the function contains multiple parameters, it is necessary to select from these parameters.In the premise of the same variation range for each parameter, we select the parameters that are the most influential to o z and Ψ values, take them as the optimiza- tion parameters through observing the changes range in o z and Ψ by adjusting the pa- rameter values.It should be noted that we need to ensure that the number of parameters to be optimized in both functions is the same when exploring the impact of sea surface roughness parameterization schemes and universal stability functions on prediction results.Subsequently, we use the method of exhaustive search to determine the range of If the function contains multiple parameters, it is necessary to select from these parameters.In the premise of the same variation range for each parameter, we select the parameters that are the most influential to z o and Ψ values, take them as the optimization parameters through observing the changes range in z o and Ψ by adjusting the parameter values.It should be noted that we need to ensure that the number of parameters to be optimized in both functions is the same when exploring the impact of sea surface roughness parameterization schemes and universal stability functions on prediction results.Subsequently, we use the method of exhaustive search to determine the range of parameters to be optimized under the condition that the function's variation trend does not undergo abrupt changes.
Finally, we divide the dataset into stable and unstable condition samples by using Equation (1) and separately optimize both types of the two samples.The sea surface roughness parameterization selected in step two and the original parameters for universal stability functions are used as initial values for the PSO and SA algorithm.The parameter range to be optimized is set as the upper and lower bounds for the optimization algorithm.Parameters are optimized based on fitness Functions ( 29) and ( 30) and iterative iterations are performed to obtain the parameter values corresponding to the minimum fitness value.Unlike the fitness function in step one, since M re f erence may have measurement errors at different heights, we use the fitted and modified refractivity M measurement as the reference value.By comparing and analyzing the accuracy optimization of the refractivity profiles with the two optimization algorithms, the optimization algorithm with the best optimization results is selected to improve the accuracy of predicting refractivity profiles in the specific sea area.

Shipboard Data Acquisition Experiment
The research of the evaporation duct prediction model relies on observation data.Many studies have set up meteorological gradient towers near the coast to collect the necessary meteorological and hydrological data.However, this method can only collect data at a single location and the sensor measurement accuracy is not high.These data cannot represent the spatiotemporal characteristics of the specific sea area.Highprecision meteorological and hydrological observations near the sea surface are crucial for optimizing the evaporation duct prediction model.Due to the lack of high-precision meteorological and hydrological observation data in the specific sea area, this study installed high-precision meteorological and hydrological sensors at five different heights above the sea surface on the research vessel "Xiangyanghong 18" and conducted a 23-day observation of meteorological and hydrological data near the sea surface.High-precision observation data of wind speed, atmospheric temperature, atmospheric humidity, sea surface skin temperature, atmospheric pressure and other parameters were obtained under different spatiotemporal conditions.
To reduce the impact on the measurement accuracy of sensors, it is necessary to ensure that the installation locations of sensors are as open and unobstructed as possible.As shown in Figure 4 The sensor locations as shown in Figure 4.The specifications of the sensors as shown in Table 1.All high-precision observation data are sent to the data collector at the ship's data center.In this paper, we optimize the evaporation duct forecasting model by using meteorological and hydrological observation data that were collected by the research vessel "Xiangyanghong 18" from 1 April to 23 April, as shown in the cruise route map in Figure 5 (cruise number NORC2021-02-1).
Mathematics 2024, 12, x FOR PEER REVIEW

Shipboard Data
The meteorological hydrological sensor used in this study has a data sampling frequency of 1 Hz.A total of 1,987,200 data points were collected from 1 April 2021, 00:00:00 to 23 April 2021, 23:59:59 at the second-level resolution.Before inputting the collected observation data into the evaporation duct prediction model, it is necessary to perform a moving average on the instantaneous observation data to eliminate sensor incidental variations and noise data in the dataset (such as spikes and missing values); thus, we can obtain the observation data that represent the environmental characteristics of the sea area.In this study, a 10 min sliding average window was used for data preprocessing.For

Shipboard Data
The meteorological hydrological sensor used in this study has a data sampling frequency of 1 Hz.A total of 1,987,200 data points were collected from 1 April 2021, 00:00:00 to 23 April 2021, 23:59:59 at the second-level resolution.Before inputting the collected observation data into the evaporation duct prediction model, it is necessary to perform a moving average on the instantaneous observation data to eliminate sensor incidental variations and noise data in the dataset (such as spikes and missing values); thus, we can obtain the observation data that represent the environmental characteristics of the sea area.In this study, a 10 min sliding average window was used for data preprocessing.For true wind observation data, it was first decomposed into north wind and east wind components, then subjected to sliding averaging and synthesis to obtain the required average wind speed and wind direction data.The processed hourly data at whole-hour timestamps were selected as the dataset for further analysis, totaling 552 data points (from 1 April 2021, 00:00:00 to 23 April 2021, 23:00:00).
Figure 6 shows the temporal distribution of data after processing for NORC2021-02-1 cruise.Figure 6a   From Section 2.2, it can be seen that different sea surface roughness parameterization schemes employ various methods to derive sea surface roughness.Some sea surface roughness parameterization schemes (such as Oost, TY01, Porchetta, Edson) utilize wave data, such as significant wave height, wave period, wave direction, etc., to calculate sea surface roughness.In this study, since the shared research cruise did not include sensors for collecting wave data, to ensure data quality, the necessary wave data were obtained from ERA5 reanalysis data [38].From Section 2.2, it can be seen that different sea surface roughness parameterization schemes employ various methods to derive sea surface roughness.Some sea surface roughness parameterization schemes (such as Oost, TY01, Porchetta, Edson) utilize wave data, such as significant wave height, wave period, wave direction, etc., to calculate sea surface roughness.In this study, since the shared research cruise did not include sensors for collecting wave data, to ensure data quality, the necessary wave data were obtained from ERA5 reanalysis data [38].
For layered observations of temperature and humidity, least squares fitting was used to derive the referenced modified refractivity profile in each time instant to act as baselines in the following optimization procedures.
where z is the height from the sea surface, i denotes the order of the measurement level, j is the sample serial number and f 0 , f 1 , f 2 are the fitting coefficients.To eliminate the refractivity noise at different heights in the subsequent optimization, the least-squares fitting process was constrained to ensure that the modified refractivity calculated by the model at the height of the ith layer is the same as the modified refractivity value measured at the corresponding layer.

Evaporation Duct Prediction Model Selection
In accordance with Step 2 in Section 3.4, the model is selected.The sea surface roughness parameterization scheme from Section 2.2 is combined with the universal stability function from Section 2.3 to obtain a new evaporation duct prediction model.The performance of the combined model is evaluated using the NORC2021-02-1 high-precise observation dataset, with Equation (32) as the evaluation criterion.The combined model with the highest prediction accuracy is selected as the model for optimization.
As shown in Figure 7, when using the same universal stability function with different sea surface roughness parameterization schemes, the combination model using the BYC universal stability function has the largest difference in the predicting rms, with an rms difference of 0.0379.The combination model using the G07 universal stability function has the smallest difference in the predicting rms, with an rms difference of 0.004.When using the same sea surface roughness parameterization scheme with different universal stability function combinations, the combination model using the Edson sea surface roughness parameterization scheme has the largest difference in predicting rms, with an rms difference of 0.1072, while the combination model using the TY01 sea surface roughness parameterization scheme has the smallest difference in predicting rms, with an rms difference of 0.073.From the above data, it can be seen that the universal stability function has a greater impact on the model's prediction results than the sea surface roughness parameterization scheme.Among the 48 predicted models for evaporation duct, the combination model using the TY01 sea surface roughness parameterization scheme and the G07 universal stability function in the specific sea area has the highest accuracy, with an rms of 0.4848.Therefore, the combination model using the TY01 sea surface roughness parameterization scheme and the G07 universal stability function is selected as the model to be optimized.

Parameters to Be Optimized
In accordance with Step 3 in Section 3.4, we determine the coefficients to be optim and their respective ranges.For the selected sea surface roughness parameterizatio TY01 and the universal stability function of G07, only variables Ψ _ M stable and Ψ H contain more than two parameters.To ensure that the number of parameters to be o mized in each function is consistent, it is necessary to select parameters containe Ψ _ M stable and Ψ _ H stable .Parameter range [0, 10] is used while keeping the other parame unchanged.Figures 8 and 9 illustrate the influence of various parameters on the valu Ψ _ M stable and Ψ _ H stable .0 0

Parameters to Be Optimized
In accordance with Step 3 in Section 3.4, we determine the coefficients to be optimized their respective ranges.For the selected sea surface roughness parameterization of TY01 and the universal stability function of G07, only variables Ψ M_stable and Ψ H_stable contain more than two parameters.To ensure that the number of parameters to be optimized in each function is consistent, it is necessary to select parameters contained in Ψ M_stable and Ψ H_stable .Parameter range [0, 10] is used while keeping the other parameters unchanged.Figures 8 and 9 illustrate the influence of various parameters on the values of Ψ M_stable and Ψ H_stable .
contain more than two parameters.To ensure that the number of parameters to be optimized in each function is consistent, it is necessary to select parameters contained in Ψ _ M stable and Ψ _ H stable .Parameter range [0, 10] is used while keeping the other parameters unchanged.Figures 8 and 9 illustrate the influence of various parameters on the values of Ψ _ M stable and Ψ _ H stable .
From Figures 8 and 9, it is evident that among the four parameters of Ψ _ M stable , the change in the value of m d has the greatest effect on the value of Ψ _ M stable , followed by the effect of m a , for a consistent range of parameter variations.Among the four parameters of Ψ _ H stable , the change in h c has the greatest influence on the value of Ψ _ H stable , and the influence of h a is the second largest.Therefore, m a and m d are selected as Ψ _ M stable parameters to be optimized, h a and h c are Ψ _ H stable parameters to be optimized.While ensuring the continuity and monotonicity of the Ψ function, the method of enumeration is used to select the range of changes for each parameter.Table 2 lists the optimized parameters and their optimization ranges for the TY01 sea surface roughness parameterization scheme and the G07 universal stability function.From Figures 8 and 9, it is evident that among the four parameters of Ψ M_stable , the change in the value of d m has the greatest effect on the value of Ψ M_stable , followed by the effect of a m , for a consistent range of parameter variations.Among the four parameters of Ψ H_stable , the change in c h has the greatest influence on the value of Ψ H_stable , and the influence of a h is the second largest.Therefore, a m and d m are selected as Ψ M_stable parameters to be optimized, a h and c h are Ψ H_stable parameters to be optimized.
While ensuring the continuity and monotonicity of the Ψ function, the method of enumeration is used to select the range of changes for each parameter.Table 2 lists the optimized parameters and their optimization ranges for the TY01 sea surface roughness parameterization scheme and the G07 universal stability function.

Optimization Results
In accordance with Step 4 in Section 3.4, we divide the observation data into samples corresponding to stable or unstable condition, taking Equation (1) as the criterion for determining atmospheric condition.This division is performed to optimize different functions of the model under different atmospheric conditions in subsequent steps.The TY01 sea surface roughness parameterization scheme and the original parameter values of the G07 universal stability function are used as the initial values for the optimization functions, with the range of optimization parameters set as the upper and lower bounds.The particle swarm optimization algorithm and the simulated annealing algorithm are employed to reference and adjust the refractivity profile M re f erence , the model output modified refractivity profile M model , based on the root mean square deviation ∆M rms_optimization as the fitness evaluation criterion.After 20 experiments, the parameter values corresponding to the minimum fitness value are determined.The optimization results are presented in Table 3. From Table 3, it can be seen that when the number of iterations is the same, both optimization algorithms can approximate the optimal value when optimizing a single function.However, when optimizing multiple functions simultaneously, PSO uses more computational resources because it is a population-based algorithm.Therefore, compared to the SA algorithm, the PSO algorithm is easier to approximate the solution of global optimum.
From Figure 10, we can see that both optimization algorithms converge within 100 iterations, and the particle swarm optimization algorithm converges better than the simulated annealing algorithm overall.The average value of the optimal fitness of the particle swarm optimization algorithm for 20 experiments is 0.3596, and its average running time is 4173.2s.The average value of the optimal fitness of the simulated annealing algorithm for 20 experiments is 0.3599, and its average running time is 4829.7 s.From this, we can conclude that optimization algorithms can be applied to improve the accuracy of predicting the modified refractivity profile.Among them, the PSO algorithm is more effective in improving the accuracy of predicting the modified refractivity profile results.
swarm optimization algorithm for 20 experiments is 0.3596, and its average running time is 4173.2s.The average value of the optimal fitness of the simulated annealing algorithm for 20 experiments is 0.3599, and its average running time is 4829.7 s.From this, we can conclude that optimization algorithms can be applied to improve the accuracy of predicting the modified refractivity profile.Among them, the PSO algorithm is more effective in improving the accuracy of predicting the modified refractivity profile results.When using the above two optimization algorithms to optimize model parameters, the optimization results of samples under unstable conditions are better than those under stable conditions.This result confirms the applicability issues of sea surface roughness parameterization schemes and universal stability functions in the specific sea area.From the optimization results, optimizing parameter Ψ H has the best effect on improving the accuracy of modified refractivity profile predictions, followed by optimizing parameter Ψ M .The impact of optimizing sea surface roughness parameters z o and z ot on the accuracy of modified refractivity profile predictions is relatively small.This conclusion is consistent with Section 3.4; namely, in the specific sea areas, the universal stability function has a greater impact on the accuracy of modified refractivity profile predictions than the sea surface roughness parameterization scheme.Since optimizing parameter Ψ H has a significantly better effect on improving the accuracy of modified refractivity profile predictions than optimizing other parameters, combination optimization including Ψ H is superior to combination optimization that only includes z o and Ψ M .

Summary
We propose a method to improve the accuracy of the prediction model of evaporation duct by using optimization algorithms in this paper.Firstly, high-precise meteorological and hydrological observation data in different heights at the sea surface are collected by installing sensors on a shipborne platform.Wave data are obtained from ERA5 reanalysis data.A least-squares fit is applied to the measured near-sea-surface refractivity profiles based on a fitting formula, resulting in a reference modified refractivity profile near the sea surface.Data collection is used to select sea surface roughness parameterization and universal stability function, yielding a more accurate evaporation duct prediction model for the specific sea area, which is the model to be optimized.The coefficients of the model are determined and coefficient selection is performed for functions with multiple coefficients, ensuring that the number of coefficients to be optimized in each function is the same, the range of values for the coefficients to be optimized is determined.
We use PSO and SA algorithms to optimize the sea surface roughness parameterization and the coefficients of the universal stability function.Comparing with the model's performance before optimization, under stable conditions, the accuracy of NORC2021-02-1 cruise prediction improved by 5.09% and 4.75% when using the PSO and SA algorithms.Under unstable conditions, the prediction accuracy increased by 9.97% with the PSO algorithm and 9.89% with the SA algorithm.These results confirm the feasibility of applying optimization algorithms to improve the accuracy of modified refractivity profile predictions.Using the PSO algorithm is more effective in enhancing the accuracy of modified refractivity profile predictions, indicating that the model optimized with the particle swarm optimization algorithm has better applicability in certain sea areas and seasons.This paper is a preliminary attempt to improve the accuracy of modified refractivity prediction using optimization algorithms, so the classical PSO algorithm and the SA algorithm are used.In subsequent studies, the accuracy of modified reflectivity prediction will be improved by using other optimization algorithms or improved variants of optimization algorithms, such as the differential evolution method, as well as Grew Wolf Optimization [39,40], Biogeography Based Optimization [41] and Social Network Optimization [42], which has been shown to outperform particle swarm algorithms and has been applied to parameter identification.

Discussion
The Monin-Obukhov similarity theory [15] pointed out that Ψ is a continuous and monotonous function, with Ψ(z/L = 0) = 0.In this study, without altering the form of Ψ and its underlying physical laws, model parameters are optimized.Experimental results have confirmed the feasibility of applying the optimization algorithm to the study of evaporation duct prediction models.Most existing research on evaporation duct assumes that the Ψ function is a monotonous function, but there is still not clear consensus on whether the Ψ function in nature is a monotonous function.In the following research, the continuity of the Ψ function will be a key focus.
In addition to obtaining more accurate evaporation duct parameters (such as modified refractivity and evaporation duct height) based on the Monin-Obukhov similarity theory, some studies have utilized machine learning algorithms, such as BP neural networks and LS-SVM, to directly predict the evaporation duct height using known meteorological and hydrological data.While this approach does not adhere to the Monin-Obukhov similarity theory, it can also lead to improved accuracy in predicting the evaporation duct height.We present a statistical analysis of the mean variation of the optimized model's evaporation duct height in this paper.
Because z o , z ot and Ψ M have a minimal impact on the accuracy of the model's predictions, it is not necessary to optimize their coefficients separately.As can be seen from Table 4, the improvement percentage of the evaporation duct height is not a simple linear relationship with the predicted accuracy of the modified refractivity.For example, when optimizing all coefficients, using the particle swarm optimization algorithm yields better results for optimizing the modified refractivity than using the simulated annealing algorithm, but the optimization effect for the evaporation duct height is the opposite.Under unstable conditions of the modified refractivity, the optimization effect is better than under stable conditions, but the optimization effect for the evaporation duct height is similarly opposite.To investigate which parameter of the evaporation duct, either the height or the modified refractivity profile, is more suitable as the accuracy criterion for predicting the evaporation duct prediction model, two groups of samples were selected to observe the propagation path loss of electromagnetic waves.The first group of samples consists of profiles with significant differences in modified refractivity profiles but relatively small differences in evaporation duct heights.The second group of samples consists of profiles with relatively small differences in modified refractivity profiles but significant differences in evaporation duct heights.The parameters set for electromagnetic wave propagation in the sea surface are as follows: the maximum range of electromagnetic wave propagation is set to 150 km, the maximum altitude is set to 100 m, of which the range step is 200 m, the altitude step is 0.1 m, the 3 dB beamwidth is 1 deg, the frequency is 3000 MHz, the dielectric constant is 69.134 and the conductivity is 7.1462 s/m.
The first group of samples shows the height of the evaporation duct and the modified refractivity profile as shown in Figure 11.The antenna height is set to 13 m, which means the antenna is located outside the evaporation duct.We calculate the path loss generated by the two profile sets.the antenna is located outside the evaporation duct.We calculate the path loss generated by the two profile sets.From Figure 12, it can be seen that the path loss is directly proportional to the propagation distance.Within a distance of 0-50 km, the path losses for the M1 and M2 profiles are similar, with no significant differences.However, in the range of 50-100 km, there is a noticeable difference in path loss between the two, with the path loss for the M2 profile significantly lower than that of M1.When the distance exceeds 100 km, the difference in path loss between the two profile groups further increases, with the path loss for M2 profile much lower than that of M1 profile.As shown in the above figure, even when the evaporation duct height is similar, the electromagnetic wave path loss exhibits significant differences due to different refractivity profile corrections.From Figure 12, it can be seen that the path loss is directly proportional to the propagation distance.Within a distance of 0-50 km, the path losses for the M1 and M2 profiles are similar, with no significant differences.However, in the range of 50-100 km, there is a noticeable difference in path loss between the two, with the path loss for the M2 profile significantly lower than that of M1.When the distance exceeds 100 km, the difference in path loss between the two profile groups further increases, with the path loss for M2 profile much lower than that of M1 profile.As shown in the above figure, even when the evaporation duct height is similar, the electromagnetic wave path loss exhibits significant differences due to different refractivity profile corrections.The second group of sample evaporation duct height and modified refractivity profiles is shown in Figure 13.The antenna height is 22 m.The path loss resulting from the two sets of profiles is calculated.The second group of sample evaporation duct height and modified refractivity profiles is shown in Figure 13.The antenna height is 22 m.The path loss resulting from the two sets of profiles is calculated.The second group of sample evaporation duct height and modified refractivity profiles is shown in Figure 13.The antenna height is 22 m.The path loss resulting from the two sets of profiles is calculated.From Figure 14, it can be seen that within a distance of 0-50 km, the path loss difference between profiles M3 and M4 is very small.When the distance exceeds 50 km, the path loss of profile M4 is higher than that of profile M3, but the loss difference is significantly smaller than the first group of profiles.The differences in refractivity profiles are From Figure 14, it can be seen that within a distance of 0-50 km, the path loss difference between profiles M3 and M4 is very small.When the distance exceeds 50 km, the path loss of profile M4 is higher than that of profile M3, but the loss difference is significantly smaller than the first group of profiles.The differences in refractivity profiles are small, even when there is a significant difference in evaporation duct height, the electromagnetic wave path loss does not have too much variation.
Mathematics 2024, 12, x FOR PEER REVIEW 24 of 26 small, even when there is a significant difference in evaporation duct height, the electromagnetic wave path loss does not have too much variation.From the above information, it can be seen that path loss is directly proportional to propagation distance, the influence of evaporation duct height on path loss is much smaller than the influence of the modified refractivity profile on path loss.In other words, the difference in path loss of electromagnetic waves mainly originates from different mod- From the above information, it can be seen that path loss is directly proportional to propagation distance, the influence of evaporation duct height on path loss is much smaller than the influence of the modified refractivity profile on path loss.In other words, the difference in path loss of electromagnetic waves mainly originates from different modified refractivity profiles.In order to obtain a predictive model for evaporation duct with higher accuracy, the modified refractivity profile is a physical quantity that must be considered and it is more suitable as a standard for evaluating the predictive accuracy of the evaporation duct prediction model.Subsequent research should focus on the modified refractivity profile to achieve high-precision prediction of evaporation duct based on predicting the modified refractivity profile.
1, b = 2.5, c = 5.3, d = 1.1.•The BYC model improves the validity of the Monin-Obukhov similarity theory under the low wind speed conditions.The function is[3] 25.  • The COARE 3.0 model, COARE 3.5 model and MM5REV model functions are

Figure 2 .
Figure 2. Model combination.rms is the root mean square deviation of the value between
, high-precision Humidity and Temperature sensors are installed at the Forecastle deck (6 m), Accommodation deck (8.3 m), Compass deck (13.1 m), Mast 1st floor (14.8 m) and Mast 3rd floor (22.3 m) to achieve layered measurements of atmospheric temperature and humidity at different heights near the sea surface.Barometers, Weather stations and Infrared thermometers are installed on the Compass deck to respectively obtain observation data for atmospheric pressure, wind speed and sea surface skin temperature.

Figure 5 .
Figure 5. Cruise trajectories, "∆" marks the start point and " " marks the end point.
represents the wind speed (u) data.Most of the time, the wind speed varies between 0 to 15 m/s, with occasional variations reaching up to 21 m/s. Figure 6b presents relative humidity (RH) data collected at five different heights above the sea surface, namely RH-1, RH-2, RH-3, RH-4 and RH-5 at 6 m, 8.3 m, 13.1 m, 14.8 m and 22.3 m heights, respectively.The RH values range from 43% to 100%, with the majority of the time having RH values above 60%.Figure 6c shows air temperature (T a ) data collected at the same five heights as RH.The five curves represent temperature variations at different heights, ranging from 10 • C to 27 • C. Except for 23 April, the temperature data are mostly below Mathematics 2024, 12, 205 14 of 23 20 • C. Figure 6d displays atmospheric pressure (P) data, which range from 1006 hPa to 1026 hPa. Figure 6e exhibits sea surface temperature (T s ) data, with a range of 12 • C to 26 • C. Finally, Figure 6f presents the Richardson number (R ib ), calculated by the collected data, which indicates atmospheric stability.According to the statistical results of the data, out of 552 hourly samples collected during the NORC2021-02-1 cruise, there are 273 samples when the atmosphere is in a stable condition and 279 samples when it is in an unstable condition, with more samples indicating an unstable atmosphere.face, namely RH-1, RH-2, RH-3, RH-4 and RH-5 at 6 m, 8.3 m, 13.1 m, 14.8 m and 22.3 m heights, respectively.The RH values range from 43% to 100%, with the majority of the time having RH values above 60%.Figure 6c shows air temperature ( a T ) data collected at the same five heights as RH.The five curves represent temperature variations at different heights, ranging from 10 °C to 27 °C.Except for 23 April, the temperature data are mostly below 20 °C. Figure 6d displays atmospheric pressure (P) data, which range from 1006 hPa to 1026 hPa. Figure 6e exhibits sea surface temperature ( s T ) data, with a range of 12 °C to 26 °C.Finally, Figure 6f presents the Richardson number ( ib R ), calculated by the collected data, which indicates atmospheric stability.According to the statistical results of the data, out of 552 hourly samples collected during the NORC2021-02-1 cruise, there are 273 samples when the atmosphere is in a stable condition and 279 samples when it is in an unstable condition, with more samples indicating an unstable atmosphere.

Figure 6 .
Figure 6.NORC2021−02−1 data, (a) is the wind speed, (b) is the relative humidity, (c) is the air temperature, (d) is the atmospheric pressure, (e) is sea surface temperature, (f) is the Richardson number.

Figure 6 .
Figure 6.NORC2021−02−1 data, (a) is the wind speed, (b) is the relative humidity, (c) is the air temperature, (d) is the atmospheric pressure, (e) is sea surface temperature, (f) is the Richardson number.

Figure 10 .
Figure 10.Convergence of PSO algorithm and SA algorithm.Figure 10.Convergence of PSO algorithm and SA algorithm.

Figure 10 .
Figure 10.Convergence of PSO algorithm and SA algorithm.Figure 10.Convergence of PSO algorithm and SA algorithm.

Table 1 .
Sensor specification and installation locations.

Table 1 .
Sensor specification and installation locations.

Table 2 .
Parameters to optimize and corresponding numeric range.

Table 2 .
Parameters to optimize and corresponding numeric range.