Optical Turbulence Proﬁle Forecasting and Veriﬁcation in the Offshore Atmospheric Boundary Layer

: A backpropagation neural network (BPNN) approach is proposed for the forecasting and veriﬁcation of optical turbulence proﬁles in the offshore atmospheric boundary layer. To better evaluate the performance of the BPNN approach, the Holloman Spring 1999 thermosonde campaigns (HMNSP99) model for outer scale, and the Hufnagel/Andrew/Phillips (HAP) model for a single parameter are selected here to estimate proﬁles. The results have shown that the agreement between the BPNN approach and the measurement is very close. Additionally, statistical operators are used to quantify the performance of the BPNN approach, and the statistical results also show that the BPNN approach and measured proﬁles are consistent. Furthermore, we focus our attention on the ability of the BPNN approach to rebuild integrated parameters, and calculations show that the BPNN approach is reliable. Therefore, the BPNN approach is reasonable and remarkable for reconstructing the strength of optical turbulence of the offshore atmospheric boundary layer.


Introduction
The atmospheric boundary layer is an important bridge for the exchange of momentum, heat, and water vapor between the earth's surface and the free atmosphere [1]. The strength of optical turbulence is high in the atmospheric boundary layer, which is a distinctive feature that distinguishes the boundary layer from other atmospheric layers. Atmospheric turbulence is the primary reason for laser beam wandering and spreading, and it leads to a serious degradation of image quality and clarity in optical applications. Atmospheric turbulence strength is generally characterized by the structural constant of the refractive index. The fluctuations of the refractive index of the atmosphere C 2 n , also known as optical turbulence, are produced by wind blowing over the earth's surface and temperature gradients between this surface and the air above it [2]. Because of the complexity of atmospheric phenomena, the C 2 n field has a complicated spatio-temporal structure. Problems such as with the intermittent vertical profile make quantitative modeling of optical turbulence a challenging task [3]. Therefore, it is of great significance to accurately estimate the optical turbulence in the atmospheric boundary layer.
Conventionally, C 2 n can be derived using several techniques such as radiosonde, radar, and optical remote sensing, among others [4][5][6][7]. Furthermore, high-frequency wind/temperature measurement instruments are widely utilized, including ultrasonic anemometers and micro-thermometers [8,9]. In 1968, one of the first experiments, the Kansas experiment, was organized outside Liberal, Kansas, USA, to solve the problem in boundary layer research. A sonic anemometer and a platinum wire sensor were installed on a 32 m tower [10]. A second experiment, known as the Minnesota experiment, was performed near Donaldson, Minnesota, USA, in 1973, where a 32 m tower was also installed and a large balloon with turbulence probes was attached to a flying cable. Since the 1970s, many investigations have been performed to determine turbulence height profiles [11]. Coulman analyzed temperature fluctuations in the lowest 200 m of the atmosphere in Culgoora, Australia, in January 1972 with the aid of resistance thermometers and anemometers used at different heights under a captive balloon [12,13]. Bufton, Azouit et al., and McHugh et al. employed two micro-thermal probes placed 1 m apart under a balloon to measure temperature fluctuations up to more than 20 km above the ground [14][15][16][17].
Over the years, several models have been developed that calculate the C 2 n profile. Firstly, the simplest empirical model can be derived from the use of segmented fitting on experimental data, requiring only the input of an altitude variable. For example, the submission laser communication (SLC) model, the Air Force Geophysics Laboratory (AFGL) and Air Force Maui Optical Station (AMOS) model, and the Critical Laser Enhancing Atmospheric Research (CLEAR I) model [18][19][20][21]. Secondly, a parameterization model can be used, based on the basic theory of turbulence to establish the relationship between conventional meteorological parameters and the profiles of C 2 n through the outer scale, such as with the Hufnagel model and Tatarski model [22,23]. Furthermore, another method to forecast the C 2 n profile is given by Felton and Alliss [24], who estimate C 2 n with the use of the weather research and forecast (WRF) numerical weather prediction (NWP) model. However, due to uncertain parameters of this model, the profiles offered by such a model are unsatisfactory. Meanwhile, based on the baroclinic energy of the atmosphere, Kovadlo and Lukin obtained C 2 n profiles of the atmospheric boundary layer using the energy spectrum of the turbulence, ranging from large-scale atmospheric fluctuations to small-scale turbulent fluctuations [25]. Each the conventional method has its own advantages and limitations, and no method is known to be superior [26]. With the development of machine learning theory, neural networks have been gaining interest [27][28][29]. Tracey and Duraisamy employed a neural network to attempt to learn the Spalart-Allmaras turbulence model [30]. This work represents an important step in the use of machine learning algorithms for generalization. Wang and Basu estimated surface-layer optical turbulence at Mauna Loa, Hawaii, and used a Monte Carlo-type approach to optimize the artificial neural network weights and subsequent prediction of C 2 n [31]. Su and Wu et al. developed a new hybrid neural network that combines the backpropagation neural net-work (BPNN) and the adaptive niche-genetic algorithm (AGA) to forecast the near-surface C 2 n [32]. Compared with conventional empirical and parametric model methods, the network method has strong nonlinear mapping capabilities for learning adaptively and is data-driven rather than reliant on empirical formulas. This research proposes a new method for estimating C 2 n profiles with the use of machine learning. The backpropagation (BP) neural network, called BPNN, is designed to forecast C 2 n profiles in the offshore atmospheric boundary layer, with special consideration of input and output variables. The BPNN is validated against in situ measurements with the spatial evolution characteristic of the C 2 n profiles. Compared with the existing model, we concluded that the BPNN approach is reasonable and remarkable. The remainder of this paper is organized as follows: Section 2 introduces three methodologies associated with the investigation of the C 2 n profiles, Section 3 provides the experimental results and discussion, including the use of contrast experiments, and finally, Section 4 concludes the paper and proposes some limitations.

BPNN
The backpropagation neural network (BPNN) was developed by Rumelhart et al. as a solution to the problem of training multi-layer perceptions [33]. It is a feed-forward neural network, and its learning process is shown in Figure 1. First, the original data need to be normalized before the initial weight distributions and thresholds of the network are randomly selected. The input and output parameters, and correction errors of each element in the middle layers, are calculated according to the activation function. Then, the connection weights and thresholds for each layer are adjusted based on acquired errors. After that, learning times are updated to determine whether the errors conform to the allowed range. If so, anti-normalization results are output; otherwise, input, output, and errors are recalculated until the error satisfies the demand. As an artificial intelligence neural network, BPNN has many advantages. First, it is fast, simple, and easy to program. Secondly, it has no parameters to tune apart from the numbers of input. Thirdly, it is a flexible method that does not require prior knowledge about the network, and finally, it does not need any special mention of the features of the function to be learned.

BPNN
The backpropagation neural network (BPNN) was developed by Rumelhart et al. as a solution to the problem of training multi-layer perceptions [33]. It is a feed-forward neural network, and its learning process is shown in Figure 1. First, the original data need to be normalized before the initial weight distributions and thresholds of the network are randomly selected. The input and output parameters, and correction errors of each element in the middle layers, are calculated according to the activation function. Then, the connection weights and thresholds for each layer are adjusted based on acquired errors. After that, learning times are updated to determine whether the errors conform to the allowed range. If so, anti-normalization results are output; otherwise, input, output, and errors are recalculated until the error satisfies the demand. As an artificial intelligence neural network, BPNN has many advantages. First, it is fast, simple, and easy to program. Secondly, it has no parameters to tune apart from the numbers of input. Thirdly, it is a flexible method that does not require prior knowledge about the network, and finally, it does not need any special mention of the features of the function to be learned. To simplify the model, the input threshold is set to 0. Suppose that given N samples, the output of the k sample at the node i of the l − 1 layer is 1 l ik Out − . Then, the input of the k sample at the node j of the l layer is: Out  In , the output of node j of the l layer of the k sample can be calculated: To simplify the model, the input threshold is set to 0. Suppose that given N samples, the output of the k sample at the node i of the l − 1 layer is Out l−1 ik . Then, the input of the k sample at the node j of the l layer is: where Out l−1 ik represents the output of the node i of the l − 1 layer, when the k sample is input. ω l−i ij represents the weight coefficient from the node i of the l − 1 layer to the node j of the l layer. After calculating In l jk , the output of node j of the l layer of the k sample can be calculated: Out l jk = f In l jk (2) Among them, f is used as an activation function and has many manifestations. This article selects the commonly used sigmoid function in BPNN as the activation function, namely: Calculate the error according to the specified expected output, where the commonly used error performance function is square: y jk is the actual output of node j in the output layer, and y jk represents the expected output. The total error is: Define δ l jk = ∂E k /∂In l jk as the partial derivative of the error of the k sample to the j node of the l layer, then the partial derivative of the connection weight is: Subsequent reverse calculation weight correction. First, calculate the weight correction from the i node in the hidden layer (the l layer) to the j node in the output layer (the l + 1 layer). Then, calculate the weight correction from the h node in the input layer (the l − 1 layer) to the i node in the hidden layer. The result is as follows: where Equation (7) represents the weight correction from the hidden layer to the output layer, and Equation (8) represents the weight correction from the input layer to the hidden layer. Finally, modify the weight, where µ is the step size: Conventionally, input data divisions are termed as training, validation, and testing data. The training section uses a specific network structure and a fixed number of hidden nodes to learn the characteristics of the data. The verification section is used to prevent the network from overfitting the data, which will inadvertently reduce the generality of its predictions. After the network is fully configured (based on training and verification), the test section will provide an independent evaluation of network performance.

C 2
n Estimation Model C 2 n estimation models that convert standard meteorological data into C 2 n vertical profiles are established based on the basic theory of turbulence. According to Kolmogorov's theory, the Tatarski C 2 n model is [34]: where a is a constant of 2.8 [35], L 0 is the outer scale, or rather, the largest scale of inertial range turbulence, and M is the potential refractive index gradient defined as: M is related to atmospheric temperature (T), atmospheric pressure (P), and potential temperature (θ). h is the height above ground. From Equations (11)- (13) it can be noted that the vertical temperature gradient can be estimated from sounding data, so choosing the appropriate L 0 is the key to estimating C 2 n . At present, the HMNSP99 model represents the most widely used outer scale model [36].
The HMNSP99 model was developed at the Air Force Research Laboratory (AFRL), Holloman Spring 1998 and Holloman Spring 1999 thermosonde campaigns (New Mexico), and the statistical relationships relating L 0 parameterization to wind shear (S) and temperature lapse rate (dT/dh) [37]. The relationships are expressed as: where S is the vertical shear of the horizontal velocity: where u and v are the north and east horizontal wind components, respectively. A vertical dependence of refractive index structure parameter was described by several authors with one of the most widely used being the Hufnagel-Valley model: where h is the height above the ground (km). W 2 is the square average of the wind speed at a height of 5-20 km. A is the constant. The typical values are W = 21 m/s, A = 1.7 × 10 −4 . Another modification of the Hufnagel model was introduced as the Hufnagel/Andrew/Phillips (HAP) model [38]: where h is the height above ground (m), W is the wind speed (m/s) and C 2 n (h 0 ) is the nominal value of refractive index structure parameter (m −2/3 ) at the height h 0 (m).

Method of Statistical Analysis
In this paper, two statistics, namely, the bias, and root mean square error (RMSE) operators, are applied to evaluate the performance and reliability of the BPNN approach as well as different models in reconstructing optical turbulence strength [39]. The bias contains information on systematic model errors, and the RMSE contains information on the statistical errors plus the systematic errors: where X i is the individual measured values of lgC 2 n , Y i is the corresponding estimated values of lgC 2 n by the neural network and model, and N is the number of samples for a couple (X i , Y i ). Furthermore, utilizing the bias and the RMSE, we retrieve the center root mean square error (CRMSE), which represents the intrinsic uncertainty not affected by the bias, and it provides fundamental information on the systematic bias and statistical uncertainties. Taking into account the actual value of C 2 n , basically in the range of 10 −18 -10 −13 m −2/3 , using lgC 2 n instead of C 2 n makes it convenient for the calculation and visualization of statis-tical results, which then makes the statistical results more comprehensible. This is simply a substituted conversion to represent the error of C 2 n , so it is reliable and valid:

Integrated Astroclimatic Parameters
The integrated value of C 2 n helps us to understand the atmospheric optical quality in terms of seeing ε, as well as the coherence length r 0 , that depends on the local conditions at altitude h. In the adaptive optics (AO) system, the seeing ε is a crucial parameter to identify the most effective window to be used to retrieve r 0 [40]. Meanwhile, the wind speed and C 2 n profiles are employed to calculate τ 0 [41]. All the parameters are given as follows to give the calculation results reference value in astronomical observation, and the point light mode is used here:

Principle of C 2 n and Measurement
The values of C 2 n measured by the micro-thermometer are derived from a pair of horizontally separated micro-temperature probe wires [42]. The refractive index fluctuation mainly results from temperature fluctuation at visible and near-infrared wavelengths [35]. According to Equation (24), the C 2 n is related to the corresponding temperature structure constant (C 2 T ): where T is the air temperature (K) and P is the air pressure (hPa). Kolmogorov defined the C 2 T as the constant of proportionality in the inertial subrange of the temperature structure function D T (r), and the C 2 n can be described as follows [43,44]: where x and r denote the position vector, . . . represents the ensemble average, and l 0 and L 0 are the inner and outer scales of turbulence and have units of m. From November to December in 2017, a field campaign of balloon-borne radiosondes equipped with micro-thermometers and GPS was conducted by the Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences in the offshore area. The balloonborne radiosondes explored the atmosphere twice a day, and all the flights were released at local time 7:00 a.m. and 19:00 p.m. every day. The scene of the launching balloon is shown in Figure 2.
The vertical profiles of C 2 n , and the structure function parameters for temperature (C 2 T ), air pressure, atmospheric temperature, relative humidity, wind speed, and wind direction were obtained by balloon-borne radiosondes. The ascent rate of the balloon was 6 m/s, with a vertical resolution of about 10 m. A micro-thermometer, which is a crucial piece of equipment for detecting the atmospheric temperature structure constant, was erected and whose equivalent noise of temperature fluctuation was about 0.002 • C [45][46][47]. Here, the frequency response range of the micro-thermometer was 0.1-30 Hz, and the corresponding statistical average time was 5 s. and 0 l and 0 L are the inner and outer scales of turbulence and have units of m.
From November to December in 2017, a field campaign of balloon-borne radiosondes equipped with micro-thermometers and GPS was conducted by the Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences in the offshore area. The balloon-borne radiosondes explored the atmosphere twice a day, and all the flights were released at local time 7:00 a.m. and 19:00 p.m. every day. The scene of the launching balloon is shown in Figure 2. The vertical profiles of 2 n C , and the structure function parameters for temperature 2 ( ) T C , air pressure, atmospheric temperature, relative humidity, wind speed, and wind direction were obtained by balloon-borne radiosondes. The ascent rate of the balloon was 6 m/s, with a vertical resolution of about 10 m. A micro-thermometer, which is a crucial piece of equipment for detecting the atmospheric temperature structure constant, was erected and whose equivalent noise of temperature fluctuation was about 0.002 °C [45][46][47]. Here, the frequency response range of the micro-thermometer was 0.1-30 Hz, and the corresponding statistical average time was 5 s.
The monthly average profiles of 2 n C , as well as the profiles of possible influencing factors, are depicted in Figure 3. In this figure we have added the Richardson number, which is used to characterize the influence of temperature stratification and wind speed shear on the development of turbulence [48]. The height reported is the vertical distance The monthly average profiles of C 2 n , as well as the profiles of possible influencing factors, are depicted in Figure 3. In this figure we have added the Richardson number, which is used to characterize the influence of temperature stratification and wind speed shear on the development of turbulence [48]. The height reported is the vertical distance above ground level in all figures of this paper. Furthermore, many studies have proved that the height of the offshore atmospheric boundary layer is below 2 km [49][50][51], so the maximum height of all profiles was fixed at 2 km in all panels. From Figure 3, we can see that the Richardson number profile has peaks when the wind shear and potential temperature gradient profile has large fluctuations at many altitudes. At the same height, the turbulent outer scale profile and the potential refractive index gradient profile also have maximum values, which can cause the appearance of a strong, turbulent, thin layer. that the Richardson number profile has peaks when the wind shear and potential temperature gradient profile has large fluctuations at many altitudes. At the same height, the turbulent outer scale profile and the potential refractive index gradient profile also have maximum values, which can cause the appearance of a strong, turbulent, thin layer.

Setting up of Contrast Experiments
In this study, the training data started from

Setting Up of Contrast Experiments
In this study, the training data started from 15 November 2017, and finished on 14 December 2017 (30 days). The predicted data started from 15 December 2017, and finished on 30 December 2017 (16 days). Using random indices, all the training data were divided into the three portions involved, training data (0.6), validation data (0.25), and testing data (0.15). Seven variables, including height, temperature, relative humidity, pressure, wind speed, potential temperature gradient, and wind shear, are inputs of the BPNN architecture. Following the normalization process of the original data, the network iterated until it reached the maximum iterations (1000), or the performance reached the preset goals. Moreover, the popular Levenberg-Marquardt backpropagation algorithm was employed as the training function, which updated the weight and bias values according to Levenberg-Marquardt optimization [52,53]. The Levenberg-Marquardt backpropagation algorithm is often the fastest backpropagation algorithm and is highly recommended as a first-choice supervised algorithm, although it does require greater memory than other algorithms.
To better evaluate the performance of the BPNN approach, the HMNSP99 model for outer scale and the HAP model for a single parameter were used here to estimate C 2 n profiles. The potential temperature gradient and wind shear, which may be the cause of atmospheric turbulence, were added to the HMNSP99 model, making its estimation of turbulence intensity more credible. The C 2 n profile obtained by the BPNN approach in Figure 4 was predicted by using the trained network and input of basic meteorological data. C 2 n profiles derived from the HMNSP99 and HAP models were evaluated using basic meteorological parameters. In addition, there are C 2 n profiles measured by the radiosonde in the figure. Each figure contains four C 2 n profiles. Each method could obtain 32 profiles, for a total of 32 figures. Due to this large number, only six consecutive comparisons of the C 2 n profile are shown in Figure 4; and as can be seen, the C 2 n vertical distribution decreases sharply in the surface layer, mostly below 100 m. In the boundary layer, the C 2 n is generally in the order of 10 −16 m −2/3 . Both the BPNN approach and HMNSP99 model respect the variability and complexity of atmospheric turbulence better than the HAP model, but comparing the BPNN approach with the HMNSP99 model, it is clearly observed that the BPNN approach outperforms the HMNSP99 model in capturing the spatial evolution of C 2 n . The BPNN approach generally agrees with the measured values both in magnitude and trend. In contrast, the estimated values from the HMNSP99 model have relatively larger drift compared with radiosonde measurements at certain heights. The high variance may be attributed to a wider variety of atmospheric conditions, but overall, the estimated values of C 2 n profiles using the BPNN approach are more consistent with the measured values than the HMNSP99 model and HAP model.

Overall Approach to Performances
To quantify the calculation error of the above BPNN approach, HMNSP99 model and HAP model, the statistical results from individual flights of lgC 2 n between measurements and estimations are demonstrated in Figure 5. As shown in Figure 5, we used the BPNN approach, HMNSP99 model and HAP model to predict 32 C 2 n profiles. Then, the error analysis of the C 2 n profile predicted by each method and the measured C 2 n profile was carried out according to Equations (18)- (20). For each method, 32 error profiles were obtained for each, including the Bias, RMSE, and CRMSE. Then, each error profile was averaged to the height, with an average value obtained. A total of 32 error average values for each method were obtained. It is noted that the bias of the BPNN approach and the HAP model have no similarities except in a few flights. Most bias values calculated from the BPNN approach were between −0.5 and 0.5, but for the HMNSP99 model and HAP model, the maximum value exceeded −1 and 1, respectively. Compared with the results obtained by the HMNSP99 model and HAP model, we can observe that the BPNN approach had the lowest RMSE and CRMSE for 32 individual meteorological balloon flights. Most RMSE and CRMSE values obtained from the BPNN approach were below 1. At the same time, most RMSE and CRMSE values from the HAP model were slightly larger than the BPNN approach, and for the HMNSP99 model, half of the values were more than 1. In general, the BPNN approach had the highest stability and showed more reliable C 2 n profiles with better statistical results than the HMNSP99 and HAP models.
parisons of the 2 n C profile are shown in Figure 4; and as can be seen, the 2 n C vertical distribution decreases sharply in the surface layer, mostly below 100 m. In the boundary layer, the 2 n C is generally in the order of 10 −16 m −2/3 . Both the BPNN approach and HMNSP99 model respect the variability and complexity of atmospheric turbulence better than the HAP model, but comparing the BPNN approach with the HMNSP99 model, it is clearly observed that the BPNN approach outperforms the HMNSP99 model in capturing the spatial evolution of 2 n C . The BPNN approach generally agrees with the measured values both in magnitude and trend. In contrast, the estimated values from the HMNSP99 model have relatively larger drift compared with radiosonde measurements at certain heights. The high variance may be attributed to a wider variety of atmospheric conditions, but overall, the estimated values of 2 n C profiles using the BPNN approach are more consistent with the measured values than the HMNSP99 model and HAP model.

Overall Approach to Performances
To quantify the calculation error of the above BPNN approach, HMNSP99 model and Most RMSE and CRMSE values obtained from the BPNN approach were below 1. At the same time, most RMSE and CRMSE values from the HAP model were slightly larger than the BPNN approach, and for the HMNSP99 model, half of the values were more than 1.
In general, the BPNN approach had the highest stability and showed more reliable 2 n C profiles with better statistical results than the HMNSP99 and HAP models. In addition, the overall statistics comparing the BNPP approach and the other models, derived from all 32 flights performed, are shown in Figure 6, with the vertical profiles of the bias, the RMSE and the CRMSE. According to Equations (18)- (20), N represents the amount of data in the horizontal direction for a couple ( ) at each precise In addition, the overall statistics comparing the BNPP approach and the other models, derived from all 32 flights performed, are shown in Figure 6, with the vertical profiles of the bias, the RMSE and the CRMSE. According to Equations (18)- (20), N represents the amount of data in the horizontal direction for a couple (X i , Y i ) at each precise height. In this section, we use N = 32 at each precise height. In detail, the bias, RMSE and CRMSE were calculated based on the C 2 n profiles with 10 m vertical resolution; then, a five-point moving average was applied to the final profiles. It can be seen from Figure 6 that the absolute values of the bias from the BPNN approach were lower than that from the HMNSP99 model and HAP model at most heights. Most bias values calculated from the BPNN approach were between −0.5 and 0.5 in the offshore atmospheric boundary layer. Note that for the HMNSP99 model there existed a larger bias value from 0 to 1500 m, and the maximum value exceeded −1. The RMSE and CRMSE of the BPNN approach were similar in value, and mainly smaller than 1.25 along the boundary layer. It is notable that the RMSE and CRMSE of the HMNSP99 model were larger than the BPNN approach, more noticeable between 500 and 1000 m. Moreover, almost all the RMSE and CRMSE of the HAP model fell between the BPNN approach and the HMNSP99 model results, which had a range from 1 to 1.25. In conclusion, the BPNN approach was more reliable than the HMNSP99 model and HAP model in reconstructing the strength of optical turbulence through individual and overall statistical analyses.  C . This formula is applicable to visible and near-infrared wavelengths. In order to maintain the consistency of the wavelength, the wavelength used in the calculation of the astroclimatic parameter was 0.55 microns. The results are shown in Figure 7. The BPNN approach provided more reasonable values, which were close to those measured by balloons, and compared with the results obtained by the BPNN approach, it can be seen that integrated astroclimatic parameters calculated from the HMNSP99 model fluctuated greatly. This is unreliable when reconstructing integrated astroclimatic parameters. In this section, we focus attention on the ability of the BPNN approach to rebuild r 0 , ε, and τ 0 . The integrated astroclimatic parameters were obtained via C 2 n values measured by balloons and estimated by the BPNN approach and the HMNSP99 model. According to Equation (24), C 2 n is obtained through C 2 T . This formula is applicable to visible and near-infrared wavelengths. In order to maintain the consistency of the wavelength, the wavelength used in the calculation of the astroclimatic parameter was 0.55 microns. The results are shown in Figure 7. The BPNN approach provided more reasonable values, which were close to those measured by balloons, and compared with the results obtained by the BPNN approach, it can be seen that integrated astroclimatic parameters calculated from the HMNSP99 model fluctuated greatly. This is unreliable when reconstructing integrated astroclimatic parameters.

Conclusions
Optical turbulence has a serious impact on laser beam propagation through the atmosphere and on astronomical observations. Effective predictions of optical turbulence strength are significant for discussing atmospheric motion and transmission characteristics. In this paper, we proposed a simple and efficient technique to predict 2 n C profiles in the offshore atmospheric boundary layer. The proposed BPNN approach was demonstrated to estimate profiles in the offshore atmospheric boundary layer without any prior formula. From the analysis and comparison of the results with the HMNSP99 model and the HAP model, the corresponding conclusions are as follows: (1) The BPNN approach identifies rules from big data without any assumptions, and so may be applied to a larger number of scenarios. Through continuous learning, the BPNN approach starts from basic meteorological parameters and obtains 2 n C profiles, avoiding the uncertainty and inaccuracy of model estimation.
(2) The BPNN approach has the highest computing precision compared with the HMNSP99 and HAP models. Most bias values calculated from the BPNN approach were between −0.5 and 0.5 in the offshore atmospheric boundary layer, and for the HMNSP99 model there existed a larger bias value from 0 to 1500 m, with the maximum value exceeding −1. The RMSE and CRMSE of the BPNN approach and the HMNSP99 model were mainly smaller than 1.25 and 1.75 along the boundary layer, respectively, and as for the HAP model, it had a range from 1 to 1.25.
(3) The profile modes of the offshore atmospheric boundary layer were induced using the BPNN approach. By calculating astroclimatic parameters, it was found that the BPNN approach provided more reasonable values, which are reliable when reconstructing integrated astroclimatic parameters.
Furthermore, the total training data were insufficient, and so the accuracy of the approach must be further improved.

Conclusions
Optical turbulence has a serious impact on laser beam propagation through the atmosphere and on astronomical observations. Effective predictions of optical turbulence strength are significant for discussing atmospheric motion and transmission characteristics. In this paper, we proposed a simple and efficient technique to predict C 2 n profiles in the offshore atmospheric boundary layer. The proposed BPNN approach was demonstrated to estimate profiles in the offshore atmospheric boundary layer without any prior formula. From the analysis and comparison of the results with the HMNSP99 model and the HAP model, the corresponding conclusions are as follows: (1) The BPNN approach identifies rules from big data without any assumptions, and so may be applied to a larger number of scenarios. Through continuous learning, the BPNN approach starts from basic meteorological parameters and obtains C 2 n profiles, avoiding the uncertainty and inaccuracy of model estimation.
(2) The BPNN approach has the highest computing precision compared with the HMNSP99 and HAP models. Most bias values calculated from the BPNN approach were between −0.5 and 0.5 in the offshore atmospheric boundary layer, and for the HMNSP99 model there existed a larger bias value from 0 to 1500 m, with the maximum value exceeding −1. The RMSE and CRMSE of the BPNN approach and the HMNSP99 model were mainly smaller than 1.25 and 1.75 along the boundary layer, respectively, and as for the HAP model, it had a range from 1 to 1.25.
(3) The profile modes of the offshore atmospheric boundary layer were induced using the BPNN approach. By calculating astroclimatic parameters, it was found that the BPNN approach provided more reasonable values, which are reliable when reconstructing integrated astroclimatic parameters.
Furthermore, the total training data were insufficient, and so the accuracy of the approach must be further improved.