E ﬀ ect of Wind Transport of Turbulent Inhomogeneities on Estimation of the Turbulence Energy Dissipation Rate from Measurements by a Conically Scanning Coherent Doppler Lidar

: A method for estimation of the turbulent energy dissipation rate from measurements by a conically scanning pulsed coherent Doppler lidar (PCDL), with allowance for the wind transport of turbulent velocity ﬂuctuations, has been developed. The method has been tested in comparative atmospheric experiments with a Stream Line PCDL (Halo Photonics, Brockamin, Worcester, United Kingdom) and a sonic anemometer. It has been demonstrated that the method provides unbiased estimates of the dissipation rate at arbitrarily large ratios of the mean wind velocity to the linear scanning speed.


Introduction
In the atmospheric boundary layer (ABL), the wind air flow is always turbulent and can be represented as a set of random vortices of various scales, with cascade energy transfer from the largest vortices to small vortices, up to the dissipation of energy into heat.In the inertial subrange of the scales of these eddies, the local spatial structure of the wind flow depends only on the dissipation rate of the kinetic energy of turbulence.Information on the dissipation rate is important for studying the spatial structure of turbulence and the dynamics of ABL and for constructing a mathematical model of ABL designed to solve a wide range of practical problems (weather forecasting, air transport safety, diffusion of atmospheric impurities, etc.).
The advent of pulsed coherent Doppler lidars (PCDLs), providing high spatial and temporal resolution of wind velocity measurements-in particular, all-fiber PCDLs such as Stream Line (Halo Photonics, Brockamin, Worcester, United Kingdom) [1] and WindCube 200s (Leosphere, France) [2]-has led to great potential for the investigation of wind turbulence in the atmospheric boundary layer [3][4][5][6][7][8][9][10][11][12][13][14].In Smalikho and Banakh's work [14], a method was proposed for the determination of such turbulence parameters as the turbulence kinetic energy, turbulence energy dissipation rate, integral scale of turbulence and momentum fluxes from lidar measurements with conical scanning (see Figure 1 or the paper [15]).In the case of conical scanning, a probing beam revolves around the vertical axis at an angle ϕ to the horizontal, with a constant angular rate ω s , and for every measurement height h, the probing volume moves along the circle of the radius R = h/ tan ϕ with a certain linear speed.It is shown in the authors' work [14] that estimates of the turbulence energy dissipation rate obtained from Stream Line lidar data using this method are in a good agreement with the dissipation rate estimates from the data of sonic anemometers.
dissipation rate obtained from Stream Line lidar data using this method are in a good agreement with the dissipation rate estimates from the data of sonic anemometers.An important feature of determination of the wind vector and wind turbulence parameters from lidar data is the need to consider the spatial averaging of the measured radial velocity over the probing volume.Just the consideration of spatial averaging provides for a correct estimation of the dissipation rate from lidar data [14].Even though the probing volume of the Stream Line lidar is small, its longitudinal dimension is z  = 30 m [16].The method proposed in the authors' work [14] assumes that the mean wind velocity is much smaller than the linear speed of motion of the probing volume along the circle of the scanning cone base at the probing height h .However, at small measurement heights and in strong wind, this condition is not always fulfilled.In this paper, the method from the authors' work [14] is generalized to the case of an arbitrary ratio between the mean wind velocity and the linear scanning speed at the probing height.The method was tested in atmospheric experiments, and the results of these tests are presented.

Basic Equations for Estimation of the Turbulence Energy Dissipation Rate from Measurements by PCDL with Conical Scanning
In the case of conical scanning, after primary processing of PCDL echo signals, we obtain an array of estimates of the radial velocity is the scan number.On the assumption that the wind is a stationary process (within one hour) and horizontally statistically homogeneous (within the circle of the scanning cone base), the mean wind vector for the height sin hR   by sine-wave least-square fitting [17].Here, z V is the vertical component of the wind velocity vector { , , } V are its horizontal components, and the angular brackets are for the ensemble averaging.Then, the array of random components of the lidar estimates of radial velocities is calculated for the same height h as ( ; ) ( ; ) ( ) ( ) where ( ) {sin ,cos cos ,cos sin }        S is the unit vector along the optical axis of the probing beam and is the average over N scans.This array is used to calculate the An important feature of determination of the wind vector and wind turbulence parameters from lidar data is the need to consider the spatial averaging of the measured radial velocity over the probing volume.Just the consideration of spatial averaging provides for a correct estimation of the dissipation rate from lidar data [14].Even though the probing volume of the Stream Line lidar is small, its longitudinal dimension is ∆z = 30 m [16].
The method proposed in the authors' work [14] assumes that the mean wind velocity is much smaller than the linear speed of motion of the probing volume along the circle of the scanning cone base at the probing height h.However, at small measurement heights and in strong wind, this condition is not always fulfilled.In this paper, the method from the authors' work [14] is generalized to the case of an arbitrary ratio between the mean wind velocity and the linear scanning speed at the probing height.The method was tested in atmospheric experiments, and the results of these tests are presented.

Basic Equations for Estimation of the Turbulence Energy Dissipation Rate from Measurements by PCDL with Conical Scanning
In the case of conical scanning, after primary processing of PCDL echo signals, we obtain an array of estimates of the radial velocity V L (θ m ;n) for different distances R from the lidar and azimuth angles θ m = m∆θ, where m = 0, 1, 2, . . ., M − 1, ∆θ is the azimuth angle resolution, (M − 1)∆θ = 360 • , and n = 1, 2, 3, 4, . . . is the scan number.On the assumption that the wind is a stationary process (within one hour) and horizontally statistically homogeneous (within the circle of the scanning cone base), the mean wind vector < V(h) >= < V z >, < V x >, < V y > is estimated from the array V L (θ m ; n) for the height h = R sin ϕ by sine-wave least-square fitting [17].Here, V z is the vertical component of the wind velocity vector V = V z , V x , V y , V x , V y are its horizontal components, and the angular brackets are for the ensemble averaging.Then, the array of random components of the lidar estimates of radial velocities is calculated for the same height h as where S(θ) = sin ϕ, cos ϕ cos θ, cos ϕ sin θ is the unit vector along the optical axis of the probing beam and < f (n where l = 1, 2, 3, 4, . . ., L << M.
During scanning, the center of the probing volume moves at a level h = R sin ϕ along the circle of the radius R = R cos ϕ and passes the distance L c = 2πR for one scan.The linear speed of the probing volume along this circle is V s = L c /T scan = ω s R , where ω s = 2π/T scan is the angular scanning rate, and T scan is the duration of one scan.According to Figure 1, the linear scanning speed V s depends on the height h, elevation angle ϕ, and time T scan as We introduced the parameter µ = U/V s , where U =|< U >| is the mean wind speed, U = V x , V y .It was shown in the authors' work [14] that, under the condition µ << 1, when the transport of turbulent inhomogeneities by the mean wind can be neglected in comparison with the linear scanning speed, the azimuth structure function (ensemble average) of radial velocity measured by a conically scanning PCDL can be represented in the form In Equation ( 4), ∆y = ∆θR = ∆θh/ tan ϕ is the transverse dimension of the probing volume (separation between the centers of neighboring probing volumes on the circle of the scanning cone base at the height h), D ⊥ (l∆y) = ε 2/3 A ⊥ (l∆y) is the transverse spatial structure function of the radial velocity averaged over the probing volume, ε is the turbulence energy dissipation rate, and σ e is the instrumental error of the radial velocity estimate.Equation ( 4) is valid for distances within the inertial subrange l∆y < L V , where L V is the integral scale of turbulence, and under the condition l∆θ ≤ π/20 = 9 • .
According to Equation (20) in [14], the function A ⊥ (l∆y) ≡ A(l∆y) can be represented in the form where is the longitudinal transfer function of the low-pass filter, is the transverse transfer function of the low-pass filter, ∆p = cσ p /2, ∆R = cT W /2, c is the speed of light, σ p = τ p / 2 √ ln 2 , τ p is the probing pulse duration, T W is the time window width, and sin c(x) = sin(x)/x [18].In our measurements by the Stream Line lidar, ∆p = 15.3 m and ∆R = 18 m.Equation (21) in the authors' work [14] is used to estimate the turbulence energy dissipation rate under the condition µ = 0 (µ << 1).Let the method for estimation of the dissipation rate by Equation (21) in [14] be called Method 1 and the estimate of the dissipation rate by Method 1 be designated as ε 1 .The method for estimation of the dissipation rate with regard to the transport of turbulent inhomogeneities by the mean wind at µ 0 is referred to as Method 2 and the corresponding estimate of the dissipation rate by Method 2 is designated as ε 2 .
By analogy with [14], we represent the estimate of the dissipation rate in the form where l > 1.The subscripts i = 1 and 2 correspond to Methods 1 and 2, respectively.The function A 1 (l∆y) ≡ A ⊥ (l∆y) is calculated by Equations ( 5)-( 8) while the function A 2 (l∆y) is calculated by the following equation (the derivation can be found in Appendix A): where the integration is performed with respect to the azimuth angle θ, It follows from Equations ( 10)-( 13) that at µ << 1, A 2 (l∆y) ≈ A ⊥ (l∆y) ≡ A 1 (l∆y) and the ratio of lidar estimates of the dissipation rate by Methods 1 and 2, is close to unity.The analysis of Equations ( 10)-( 13), taking into account Equations ( 6)- (8), shows that at fixed µ 0 and ∆y, the ratio ε 1 /ε 2 depends on the elevation angle ϕ.At µ > 1, the ratio ε 1 /ε 2 increases with an increase in µ and can be many times larger than unity.The smaller ∆y is, the larger ratio ε 1 /ε 2 , if the condition µ << 1 is not fulfilled.
Figure 2 shows the ratio of the estimates of the dissipation rate ε 1 /ε 2 as a function of the parameter µ, as calculated by Equations ( 10)-( 13) and Equations ( 6)-( 8) at different elevation angles ϕ and distances ∆y.It follows from Figure 2 that the ratio ε 1 /ε 2 begins to exceed unity markedly at µ, ranging from 0.55 (ϕ = 60 • , ∆y = 3 m) to 1.1 (ϕ = 16 • , ∆y = 12 m), depending on the scanning parameters.An increase in ∆y leads to a decrease in the ratio ε 1 /ε 2 .The larger ϕ is, the faster the increase of ε 1 /ε 2 with an increase in µ.Thus, at ϕ = 60 • , ∆y = 3 m and µ = 2.4, the estimate of the dissipation rate by Method 1 exceeds that obtained by Method 2 by four times.Consequently, calculations by Equations ( 10)-( 13) allow us to determine the limits of applicability of the approach proposed in the authors' work [14] (Method 1) at different ϕ and ∆y.
Equations ( 5)- (13) were used for analysis of the results obtained in experiments on the study of turbulence in the atmospheric boundary layer and were verified in the special experiment.
Equations ( 5)- (13) were used for analysis of the results obtained in experiments on the study of turbulence in the atmospheric boundary layer and were verified in the special experiment.

Experiment 2018
The experiment was conducted on July 6-24, 2018, at the Basic Experimental Observatory (BEO) of the Institute of Atmospheric Optics SB RAS in Tomsk, Russia (56°06′51.41"N,85°06′03.22"E),with the Stream Line PCDL.In this experiment, the scanning was carried out at two alternative elevation angles:  = 1  = 35.3°and  = 2  = 60° [19].The turbulence energy dissipation rate was estimated in this experiment by Method 1 [14] from the azimuth structure function, calculated for horizontal separations V l y L  not exceeding the integral scale of turbulence, which corresponds to the inertial subrange of turbulence.In this case, according to the Kolmogorov-Obukhov hypotheses, estimates of the dissipation rate calculated from lidar measurements at different elevation angles The estimates of the dissipation rate in the authors' work [19] were obtained by Method 1 [14], where the transport of turbulent inhomogeneities by the mean wind is neglected.Therefore, this neglect may be a possible reason for an excess of  = 60° at the height h = 70 m.In the experiment in the authors' work [19], the duration of one scan scan T was 60 s.At this scan duration, the speed of motion of the

Experiment 2018
The experiment was conducted on 6-24 July 2018, at the Basic Experimental Observatory (BEO) of the Institute of Atmospheric Optics SB RAS in Tomsk, Russia (56 • 06 51.41 N, 85 • 06 03.22 E), with the Stream Line PCDL.In this experiment, the scanning was carried out at two alternative elevation angles: ϕ = ϕ 1 = 35.3• and ϕ = ϕ 2 = 60 • [19].The turbulence energy dissipation rate was estimated in this experiment by Method 1 [14] from the azimuth structure function, calculated for horizontal separations l∆y < L V not exceeding the integral scale of turbulence, which corresponds to the inertial subrange of turbulence.In this case, according to the Kolmogorov-Obukhov hypotheses, estimates of the dissipation rate calculated from lidar measurements at different elevation angles ε(ϕ 1 ) and ε(ϕ 2 ) for the same height should coincide.However, this turned out to not quite be the case.In the atmospheric layer at heights of 200-300 m (see Figure 5b in authors' work [18]), the estimates ε(ϕ 1 ) and ε(ϕ 2 ) differed little on average.However, at heights below 150 m, ε(ϕ 2 ) far exceeded ε(ϕ 1 ).
The estimates of the dissipation rate in the authors' work [19] were obtained by Method 1 [14], where the transport of turbulent inhomogeneities by the mean wind is neglected.Therefore, this neglect may be a possible reason for an excess of < ε(ϕ 2 ) > over < ε(ϕ 1 ) > at small heights h, where the speed of motion of the probing volume during the conical scanning at an elevation angle of 60 • can be far smaller than the mean wind velocity.To check this hypothesis, we took the raw data measured by the lidar From the lidar measurements at these parameters, we calculated the mean wind velocity U(ϕ 1 ) and U(ϕ 2 ) by the sine-wave fitting and the dissipation rate ε(ϕ 2 ) and ε(ϕ 1 ) by Methods 1 and 2 using Equations ( 1), ( 2) and ( 5)- (13).We used N = 20 scans to estimate the mean velocity and the azimuth structure function by Equation ( 2).The duration of the measurements at two alternative elevation angles was a little longer than 40 min.Since ∆θ = 3 • , we took l = 3 in Equation ( 5).For the mean velocity, we obtained practically coinciding estimates for U(ϕ 1 ) and U(ϕ 2 ).
Figure 3 depicts the time series of the velocity ratio µ(ϕ i ) = U(ϕ i )/V s (ϕ i ) and turbulence energy dissipation rate ε(ϕ i ) calculated by Methods 1 and 2 for elevation angles of 35.3 • (i = 1) and 60 • (i = 2) at the height h = 70 m.It can be seen that, at an elevation angle of 35.3 • , the parameter µ ranges from 0.5 to 0.8.According to curve 2 in Figure 2, the results of estimation of the dissipation rate by Methods 1 and 2 should not differ widely in this case.At an elevation angle of 60 • , the parameter µ exceeds unity and can achieve 1.9.In this case, according to curve 2 in Figure 2, Method 1 can overestimate the dissipation rate by more than double.Figure 3b confirms this conclusion.If Method 2 is used, the estimates of the dissipation rate at elevation angles of 35.3 • and 60 • coincide on average (see Figure 3c).This fact indicates indirectly the applicability of Method 2 for determination of the turbulence energy dissipation rate from measurements by a conically scanning lidar at an arbitrary ratio between the mean wind velocity and the linear speed of motion of the probing volume.1), ( 2) and ( 5)- (13).We used N = 20 scans to estimate the mean velocity and the azimuth structure function by Equation (2).The duration of the measurements at two alternative elevation angles was a little longer than 40 min.Since   = 3°, we took l = 3 in Equation ( 5).For the mean velocity, we obtained practically coinciding estimates for  exceeds unity and can achieve 1.9.In this case, according to curve 2 in Figure 2, Method 1 can overestimate the dissipation rate by more than double.Figure 3b confirms this conclusion.If Method 2 is used, the estimates of the dissipation rate at elevation angles of 35.3° and 60° coincide on average (see Figure 3c).This fact indicates indirectly the applicability of Method 2 for determination of the turbulence energy dissipation rate from measurements by a conically scanning lidar at an arbitrary ratio between the mean wind velocity and the linear speed of motion of the probing volume.

Experiment 2019
The lidar experiment on the study of turbulence of the statically stable boundary layer was conducted in the coastal zone at the western coast of Lake Baikal near Listvyanka (52 • 50 47 N, 104 • 53 31 E) on 7-24 August of 2019.During this experiment, the Stream Line PCDL was installed at the territory of the Baikal Astrophysical Observatory of the Institute of Solar-Terrestrial Physics SB RAS, a few tens of meters from the building of the Big Solar Vacuum Telescope at a height of 180 m above the Baikal level, with a minimal separation of 340 m from the lake (see Figure 4 in authors' work [15]).A sonic anemometer AMK-03 (Sibanalitpribor, Tomsk, Russia) with a sampling rate of 80 Hz was installed 50 m from the lidar at the top of a mast, at a height of 15 m above the lidar level.The lidar measurements were conducted with conical scanning at the elevation angle ϕ = 60 • .The scanning duration and the azimuth resolution were the same as in experiment 2018 and amounted to T scan = 1 min and ∆θ = 3 • , respectively.
In the experiment, a part of the lidar data was obtained under strong wind in the lower 100 m of the atmospheric layer above the lidar level.In particular, from 20:00 LT of 19 August to 02:30 LT of 20 August, the measurements were conducted at mean wind velocity U at a height of 55 m, ranging from 8 m/s to 12.7 m/s.The linear speed of the probing volume motion at a height of 55 m in the case of scanning at an elevation angle of 60 • is V s = 3.3 m/s.We used these data to check Method 2 in terms of its estimation of the turbulence energy dissipation rate.To obtain estimates of the mean velocity and the dissipation rate ε 1 and ε 2 (respectively, by Methods 1 and 2) by Equation ( 9), in Equations ( 1) and ( 2), we set N = 30 consecutive scans (30 min average).The same 30 min average was applied to obtain estimates of the mean wind velocity and the turbulent energy dissipation rate from measurements of the sonic anemometer.We obtained dissipation rate estimates ε s from temporal structure function D (Uτ) of the longitudinal component of the wind velocity vector measured by the sonic anemometer, where 1 s ≤ τ ≤ 3 s, with the use of Taylor' hypothesis of "frozen" turbulence (see Equation ( 26) in authors' work [14]).Figure 4 demonstrates the results obtained from the lidar measurements.One can see (Figure 3e) that the estimate ε 1 can exceed the estimate ε 2 by an order of magnitude.Using the data from Figure 4a,b, we drew the vertical and temporal profiles of the mean wind velocity U and the velocity ratio µ = U/V s , which are shown as solid curves in Figure 5.The dashed straight lines in this figure correspond to the linear scanning speed V s at different heights h.We can see that, in the 300 m layer adjacent to the ground, the parameter µ decreases with height despite an increase in the velocity U. Above h = 300 m, this parameter takes on values at which the estimates of the dissipation rate ε 1 and ε 2 differ a little from each other.At the height h = 55 m, the parameter µ at its maximum approaches four (see the blue curve in Figure 5d), while ε 1 becomes an order of magnitude larger than ε 2 .Figure 6 shows the time series of the dissipation rate estimates ε 1 (t) and ε 2 (t) obtained from lidar data with Methods 1 and 2, respectively, at heights of 55 m and 100 m, along with ε s (t) calculated from the data measured by the sonic anemometer at a height of 15 m.One can see that if the transport of turbulent inhomogeneities by the mean wind is neglected (Method 1), then the estimates ε 1 (t) at a height of 100 m are close to the estimates ε s (t) at a height of 15 m, while those at a height of 55 m exceed by approximately five times the estimates obtained from the sonic anemometer measurements at a height of 15 m.Although the measurements were conducted under strong wind and intense turbulence, the estimates ε 1 ∼1 m 2 /s 3 obtained from the lidar data at a height of 55 m and ε 1 ∼0.2 m 2 /s 3 at a height of 100 m are obviously overestimated.The use of Method 2 for estimation of the dissipation rate from the lidar data provides the more likely result (Figure 6b).At heights of 55 m and 100 m, the dissipation rate takes on smaller values than those calculated for the height h = 15 m from the sonic anemometer data.They decrease with height and are comparable in terms of the absolute value with the known experimental data (see, for example, [20,21]).
The vertical profiles of the estimates of the turbulent energy dissipation rate obtained from the lidar data with Methods 1 and 2 are shown in Figure 7 as green and blue curves, respectively.It can be seen from Figure 7e that, at a height of 55 m, the estimate 1  is ten times larger than the estimate  The use of Method 2 for estimation of the dissipation rate from the lidar data provides the more likely result (Figure 6b).At heights of 55 m and 100 m, the dissipation rate takes on smaller values than those calculated for the height h = 15 m from the sonic anemometer data.They decrease with height and are comparable in terms of the absolute value with the known experimental data (see, for example, [20,21]).
The vertical profiles of the estimates of the turbulent energy dissipation rate obtained from the lidar data with Methods 1 and 2 are shown in Figure 7 as green and blue curves, respectively.It can be seen from Figure 7e that, at a height of 55 m, the estimate ε 1 is ten times larger than the estimate ε 2 .Red squares indicate the estimates of the dissipation rate from measurements by the sonic anemometer at a height of 15 m.In this experiment, the minimum height of measurement by the lidar is h 0 ≈ 55 m and the height step is ∆h ≈ 15 m.To compare the estimates of the dissipation rate ε 2 (h s ) (Method 2) and ε s obtained from measurements, respectively, with the lidar and sonic anemometer at a height of h s = 15 m, we used the extrapolation of the height profile ε 2 (h k ) (h k = h 0 + k∆h; k = 0, 1, 2, 3) from the layer of 55-100 m to the height h s as follows.First, for each vertical profile ε 2 (h k ) shown in Figure 7 within the layer 55-100 m, we used a linear fitting of lg(ε 2 (h k )) = αh k + β, where α is the slope and β is the intercept, to lg(ε 2 (h k )), by the least squares method.We then calculated values ε 2 (h s ) for the height of the sonic anemometer position.In Figure 7, in the 15-100 m layer, the dashed lines show the profiles ε 2 (h), where h ∈ [h s , h 3 ].Calculations of the value ε 2 (h s )/ε s − 1 ×100% showed that it does not exceed 20%.Thus, if, in the 15-100 m layer, the vertical profile lg(ε 2 (h)) has a height dependence which is close to linear, the lidar estimates of the dissipation rate obtained by Method 2 are in good agreement with the results of measurements with the sonic anemometer (at least, the discrepancy of the results are within the statistical error, examples of which are given below).

Experiment 2020
An experiment focusing on the testing of Method 2 for estimation of the turbulent energy dissipation rate from the data measured by a conically scanning PCDL was carried out in May 2020 at the territory of the BEO.Continuous lidar measurements were performed in two stages: (1) from 10 to 15 May and (2) from 19 to 26 May 2020.The results of estimation of the dissipation rate by Methods 1 and 2 in the experiment were compared with the estimates ε from the data of the sonic anemometer AMK-03 with sampling rate of 80 Hz (Sibanalitpribor, Tomsk, Russia), installed on a 43 m tower with boom length of 2 m at a height of 42 m.The separation between the Stream Line lidar and the tower was 160 m.
The elevation angle ϕ was set equal to 16 • in the measurements, so that the center of probing pulses achieved the height h = 42 m at the distance R = 152 m from the lidar.At a height of 42 m, the radius and the length of the scanning circle were R = Rcosϕ = 146 m and L c = 2πR = 918 m, respectively.The scan duration T scan was set to 5 min.The linear speed of motion of the probing volume along the scanning circle V s = L c /T scan was ≈3 m/s.To estimate one value of the radial velocity, N a = 15,000 laser shots were used for the raw data accumulation, so that the duration of measurement for every azimuth angle was ∆t = N a / f p = 1 s.For one scan, we obtained M = T scan /∆t = 300 estimates of the radial velocity at different azimuth angles (300 rays), with the resolution ∆θ = 360 • /M = 1.2 • and the distance between them along the circle (between centers of probing volumes) ∆y = ∆θR ≈ 3 m.
To obtain estimates of the mean wind velocity and the dissipation rate ε 1 and ε 2 (by Methods 1 and 2, respectively) using Equation ( 9), for the averages in Equations ( 1) and ( 2), we used N = 7 consecutive scans (approximately 35 min average).In Equation ( 9), we took l = 3.The estimates of the dissipation rate ε s from measurements by the sonic anemometer were determined from temporal spectra of the longitudinal (along the mean wind direction) component of the wind vector within the inertial subrange from 0.2 to 2 Hz, where the spectrum was characterized by the "-5/3" power dependence on frequency, with the use of Taylor's "frozen" turbulence hypothesis.In addition, the dissipation rate was estimated with the use of the time series of the wind velocity measured by the sonic anemometer for 35 min.
However, we were able to use the data of the sonic anemometer measured only at a certain interval of wind direction angles, due to the distortions (shading) introduced by the mast.Therefore, in the first part of the experiment (from 10 May to 15 May), only those raw data of the sonic anemometer were suitable for further processing, and these were measured within a period of 17.5 h, starting at 07:40 local time on 13 May 2020.
Figure 8 exemplifies the time series of the mean wind velocity obtained from measurements by the lidar U and by the sonic anemometer U s for 17.5 h starting from 07:40 LT on 13 May 2020.One can see that approximately half of the time, the mean wind velocity U (U s ) exceeded the linear speed of motion of the probing volume V s = 3 m/s.It can be seen that the greatest difference between the estimates of the dissipation rate by Methods 1 and 2 was observed at µ > 2, which is in agreement with the calculated results shown in Figure 2. It follows that the estimates of the dissipation rate obtained by Method 2 from the lidar data and from the sonic anemometer data are in good agreement at all values of the velocity ratio µ given in Figure 8a.The application of Method 1 would lead to a threefold overestimation of the dissipation rate at 17:00 LT.
Using the data of Figure 8, we compared the estimates of the mean wind velocity and the dissipation rate from simultaneous measurements by the sonic anemometer (U s and ε s ) and the lidar (U and ε 2 ).The results of this comparison are shown in Figure 9. On the assumption of statistical independence of the estimates and equal errors in estimation of the mean wind velocity from lidar data U and from sonic anemometer data U s , as well as statistical independence of the estimates and equal errors in estimation of the dissipation rate ε 2 and ε s , we calculated the relative bias B U =< η U > p ×100% and the relative error )/2 × 100% of estimation of the wind velocity along with the relative bias B ε =< η ε > p ×100% and the relative error estimation of the dissipation rate.The operator < . . .> p means an average over all pairs of estimates (U, U s ) or (ε 2 , ε s ), where p is the number of such pairs.To calculate η , we used the data of Figure 9a,b, respectively.Thereby, from the data of Figure 9, we obtained B U = −2%, E U = 4%, B ε = −9%, and E ε = 14%.

       
, we used the data of Figure Using data from Figure 9b, we found that the coefficient of correlation between log(ε s ) and log(ε 2 ) equals 0.987; the results of linear regression are a = 1.0346 and b = 0.0261 (Y i = aX i + b, where X i = lg(ε s ), Y i = lg(ε 2 ), i = 1, 2, 3, . . ., I, I = 205, a and b are the slope and intercept of this linear equation); the coefficient of determination R 2 = 0.98 and p-value as a measure of statistical significance equals 0.03.
In the second part of the experiment (measurements from 19 May to 26 May 2020), the wind direction was continuous for 62 h so that there was no distortion of the wind data due to "shading" by the mast on which the sonic anemometer was installed.Consequently, we were able to compare the results of measurements by the lidar and sonic anemometer for this period of time.
Figure 10 depicts the 62 h time series of the ratio µ of the mean wind velocity to the linear speed of the probing volume and the estimates of the dissipation rate from the measurements by the sonic anemometer and the lidar starting from 16:00 LT on 20 May of 2020.It can be seen that, at µ > 1, Method 1 markedly (up to 2.7 times) overstates the dissipation rate in comparison with ε 2 and ε s .Using the data of Figure 10b, we compared the estimates of the dissipation rate obtained from measurements by the sonic anemometer ε s with the lidar estimates ε 1 and ε 2 .The results of the comparison are shown in Figure 11.It can be seen from Figure 11a that, at strong turbulence, (when ε s ≥ 10 −2 m 2 /s 3 ) and µ > 1, the lidar estimates of the dissipation rate ε 1 , on average, exceed the estimates ε s obtained from measurements by the sonic anemometer.Under the same conditions, the difference between the estimates ε 2 and ε s (Figure 11b) is much smaller.The relative bias of estimates of the dissipation rate calculated by Methods 1 and 2 as judged from the data of Figure 11 is, respectively, B ε = 15% and B ε = −2%.This result allows us to assert that Method 2 gives an unbiased estimate of the dissipation rate at an arbitrary ratio of the mean wind velocity to the linear speed of conical scanning.According to the data of Figure 11, the relative statistical errors of the lidar estimates ε 1 and ε 2 differ less widely, at E ε = 20% and E ε = 16%, respectively.Parameters of the linear regression, including the slope a and the intercept b (equal to 1.134 and 0.382 in the case of Method 1 (Figure 11a; in this case, for calculation of a and b, we used lg(ε 1 ) instead of lg(ε 2 )) and 0.996 and −0.026 in the case of Method 2 (Figure 11b).The coefficient of the correlation between lg(ε s ) and lg(ε 2 ) equals 0.98.Using data from Figure 11b, we determined that the coefficient of determination R 2 = 0.962 and p-value as a measure of statistical significance equals 0.008.
Thus, the comparative analysis of the results of measurements by a conically scanning lidar and a sonic anemometer at the same height, carried out in this section, showed that taking into account the transfer of turbulent inhomogeneities by the mean wind (Method 2) makes it possible to obtain unbiased lidar estimates of the dissipation rate at arbitrary values of parameter µ.
It should be noted that the results presented in Figure 3 (experiment 2018), Figure 4 (experiment 2019) and Figures 8 and 10 (experiment 2020) were obtained from lidar measurements at relatively large signal-to-noise ratios SNR (the ratio of the mean power of the lidar echo signal to the mean power of the detector noise in the receiver passband of 50 MHz).Consequently, the instrumental error in estimating the radial velocity σ e , which substantially depends on SNR and N a , should be small, and this error should insignificantly affect the accuracy of the estimate of the turbulent energy dissipation rate.Therefore, in the experiments of 2018 (N a = 7500) and 2020 (N a = 15,000), the signal-to-noise ratio varied from −10 dB to 0 dB and the error σ e did not exceed 0.1 m/s (σ e was calculated by Equation (23) in authors' work [14] after replacing A(l∆y) with A 2 (l∆y)), and in the 2019 experiment (N a = 7500), it varied from −13 dB to −8 dB and σ e ≤ 0.2 m/s.

Discussion
Pulsed coherent Doppler lidars capable of measuring the radial velocity with high spatial and temporal resolution [1,2,22] are effective technical means of obtaining information on the turbulent energy dissipation rate in the atmospheric boundary layer (in particular, the Stream Line lidar used in this work).The method developed in the authors' work [14] for retrieving the vertical profiles of the dissipation rate from measurements with a conically scanning PCDL has a number of advantages over the previously used approaches [4][5][6][7][8][9][10][11][12].In particular, it allows one to take into account the effect of averaging the radial velocity over the probing volume.However, the application of this method assumes a significant excess of the linear speed of the conical scanning V s over the average wind velocity U (i.e., µ = U/V s << 1), which is not always the case in practice.Indeed, if the dissipation rate is determined from measurements at the height of the atmospheric surface layer at a relatively large elevation angle ϕ and, moreover, with a strong wind, then parameter µ can be many times greater than unity.In this case, the distribution of the radial velocity within the probing volume will change with time, mainly due to the transfer of turbulent inhomogeneities by the mean wind and not due to the movement of the probing volume around the base of the scanning cone.
To obtain an estimate of the dissipation rate (using Equation ( 9)) from the lidar measurements of the azimuth structure function of the radial velocity at any velocity ratio µ, it is necessary to know the equation for the function A 2 (l∆y), which would take into account (1) the motion of the probing volume (on the horizontal plane) during conical scanning, (2) wind transfer of turbulent inhomogeneities and (3) averaging of the radial velocity over the probing volume.In this work, using Taylor's hypothesis of "frozen" turbulence and a number of approximations (see the derivation of the equation in Appendix A), we obtained an equation for the function A 2 (l∆y), taking into account all these three factors (see Equations ( 10)-( 13)).To test this method (Method 2, in which the function A 2 (l∆y) is used in Equation ( 9)), we used the data of lidar experiments, in which any values of the parameter µ were realized, including µ >> 1.
Testing of Method 2 in the 2018 experiment (see Section 3) showed that the reason for the discrepancy in the estimates of the turbulent energy dissipation rate from the measurements by the Stream Line lidar in the atmospheric surface layer at different elevation angles (ϕ 1 = 35.3• and ϕ 2 = 60 • ) is that the method used in the authors' work [19] (Method 1) does not take into account the transfer of turbulent inhomogeneities by the mean wind.The use of Method 2 made it possible to obtain similar results from lidar measurements at angles ϕ 1 and ϕ 2 .
The use (along with the Stream Line lidar) of sonic anemometers in the 2019 (Section 4) and 2020 (Section 5) experiments made it possible to compare the estimates of the dissipation rate from the measurements by the lidar (estimates ε 1 and ε 2 obtained by Methods 1 and 2, respectively) and the sonic anemometer (estimate ε s ).Since, in the 2019 experiment, the height of the sonic anemometer position was 40 m less than the minimum height of the lidar measurement, to compare the estimates of the dissipation rate, we used the extrapolation of the vertical profile ε 2 (h) to the level of the sonic anemometer and also obtained good agreement between the results, including the case at µ >> 1 (see Figure 7).The 2020 experiment, whose data we used to directly compare the dissipation rate estimates from measurements with a lidar and a sonic anemometer at the same height, convincingly proved the applicability of Method 2 for any values of parameter µ.
It is worth using the generalized method at both large wind velocities and small probing heights, up to 200-300 m.At heights above the surface atmospheric layer, the method in the authors' work [14] usually provides quite acceptable (in accuracy) estimates of the dissipation rate.Its application does not require complicated calculations in contrast to the generalized method (compare Equation (20) in [14] with Equations ( 10)-( 13) of this paper).

Figure 1 .
Figure 1.Geometry of measurement by a pulsed coherent Doppler lidar with conical scanning by laser beam.

Figure 1 .
Figure 1.Geometry of measurement by a pulsed coherent Doppler lidar with conical scanning by laser beam.
the average over N scans.This array is used to calculate Remote Sens. 2020, 12, 2802 3 of 20 the azimuth structure function D L (l∆θ) of radial velocity fluctuations (the function averaged over all azimuth angles θ m ) by the following equation:

Figure 2 .
Figure 2. Ratio of lidar estimates of the dissipation rate by Methods 1 and 2 ( 12 /  ) as a function of the ratio of the mean wind velocity to the linear scanning speed (  ) at an elevation angle of 16° (curves 1), 35.3° (curves 2) and 60° (curves 3) and y  = 3 m (solid curves) and 12 m (dashed curves).
for the same height should coincide.However, this turned out to not quite be the case.In the atmospheric layer at heights of 200-300 m (see Figure5bin authors' work[18]differed little on average.However, at heights below 150 m,

1 
the speed of motion of the probing volume during the conical scanning at an elevation angle of 60° can be far smaller than the mean wind velocity.To check this hypothesis, we took the raw data measured by the lidar from 19:40 local time (LT) on July 23, 2018 to 08:20 LT on July 24, 2018, at elevation angles = 35.3°and 2

Figure 2 .
Figure 2. Ratio of lidar estimates of the dissipation rate by Methods 1 and 2 (ε 1 /ε 2 ) as a function of the ratio of the mean wind velocity to the linear scanning speed (µ) at an elevation angle of 16 • (curves 1), 35.3 • (curves 2) and 60 • (curves 3) and ∆y = 3 m (solid curves) and 12 m (dashed curves).
from 19:40 local time (LT) on 23 July 2018 to 08:20 LT on 24 July 2018, at elevation angles ϕ 1 = 35.3• and ϕ 2 = 60 • at the height h = 70 m.In the experiment in the authors' work [19], the duration of one scan T scan was 60 s.At this scan duration, the speed of motion of the probing volume V s (Equation (3)) at a height of 70 m was 10.4 m/s and 4.2 m/s at elevation angles of 35.3 • and 60 • , respectively.With regard to the short time δt ≈ 1 s needed for alternation of the elevation angle between 35.3 • and 60 • , the duration of one cycle T circl = 2(T scan + δt) in the experiment [19] exceeded 2 min by a little.The azimuth resolution was ∆θ = 360 • /M = 3 • , where M = T scan /∆t = 120 is the number of radial velocity estimates for one scan, ∆t = N a / f p = 0.5 s is the time of measurement for estimation of the radial velocity for every azimuth angle, N a = 7500 is the number of laser shots used for the raw data accumulation, and f p = 15 kHz is the pulse repetition frequency for the Stream Line lidar.
a little.The azimuth resolution was   = 360°/M = 3°, where scan / M T t  = 120 is the number of radial velocity estimates for one scan, / ap t N f  = 0.5 s is the time of measurement for estimation of the radial velocity for every azimuth angle, a N = 7500 is the number of laser shots used for the raw data accumulation, and p f = 15 kHz is the pulse repetition frequency for the Stream Line lidar.From the lidar measurements at these parameters, we calculated the mean wind velocity 1 () U  and 2 () U  by the sine-wave fitting and the dissipation rate

Figure 3
Figure 3 depicts the time series of the velocity ratio ( ) ( ) / ( ) i i s i UV      and turbulence energy dissipation rate () i  calculated by Methods 1 and 2 for elevation angles of 35.3° ( i = 1) and 60° ( i = 2) at the height h = 70 m.It can be seen that, at an elevation angle of 35.3°, the parameter  ranges from 0.5 to 0.8.According to curve 2 in Figure 2, the results of estimation of the dissipation rate by Methods 1 and 2 should not differ widely in this case.At an elevation angle of 60°, the parameter 

Figure 3 .
Figure 3.Time series of the ratio of the mean wind velocity to the linear scanning speed  (a) and estimates of the turbulent energy dissipation rate obtained by Methods 1 (b) and 2 (c ) at a height of 70 m.Scanning at an elevation angle of 35.3° (blue curves) and 60° (red curves).Measurements were taken at the BEO from 19:40 LT of 23 July 2018to 08:20 LT of 24 July 2018.

Figure 3 .
Figure 3.Time series of the ratio of the mean wind velocity to the linear scanning speed µ (a) and estimates of the turbulent energy dissipation rate obtained by Methods 1 (b) and 2 (c) at a height of 70 m.Scanning at an elevation angle of 35.3 • (blue curves) and 60 • (red curves).Measurements were taken at the BEO from 19:40 LT of 23 July 2018 to 08:20 LT of 24 July 2018.

Figure 4 .
Figure 4. Height-temporal distributions of the mean wind velocity (a), ratio of the mean wind velocity to the linear scanning speed  (b), estimates of the turbulent energy dissipation rate obtained by Methods 1 (c) and 2 (d) and the ratio of these estimates (e).Measurements were taken on the coast of Lake Baikal from 19:45 LT of 19 August to 02:15 LT of 20 August 2019.

Figure 4 .
Figure 4. Height-temporal distributions of the mean wind velocity (a), ratio of the mean wind velocity to the linear scanning speed µ (b), estimates of the turbulent energy dissipation rate obtained by Methods 1 (c) and 2 (d) and the ratio of these estimates (e).Measurements were taken on the coast of Lake Baikal from 19:45 LT of 19 August to 02:15 LT of 20 August 2019.

22 Figure 5 .
Figure 5. Vertical profiles (a,b) at 21:00 LT (brown curves), 23:00 LT (pale-blue curves) and 01:00 LT (pink curves) and time series (c,d) at a height of 55 m (blue curves), 100 m (red curves) and 300 m (green curves) for the mean wind velocity (a,c) and the ratio of the mean wind velocity to the linear speed of the probing volume motion  (b,d).The data are taken from Figure 4.The dashed lines show the linear scanning speeds at the corresponding heights.

Figure 6
Figure 6 shows the time series of the dissipation rate estimates 1 () t  and 2 () t  obtained from lidar data with Methods 1 and 2, respectively, at heights of 55 m and 100 m, along with () s t  calculated from the data measured by the sonic anemometer at a height of 15 m.One can see that if the transport of turbulent inhomogeneities by the mean wind is neglected (Method 1), then the estimates 1 () t  at a height of 100 m are close to the estimates () s t 

Figure 5 .
Figure 5. Vertical profiles (a,b) at 21:00 LT (brown curves), 23:00 LT (pale-blue curves) and 01:00 LT (pink curves) and time series (c,d) at a height of 55 m (blue curves), 100 m (red curves) and 300 m (green curves) for the mean wind velocity (a,c) and the ratio of the mean wind velocity to the linear speed of the probing volume motion µ (b,d).The data are taken from Figure 4.The dashed lines show the linear scanning speeds at the corresponding heights.

22 Figure 6 .
Figure 6.Time series of the turbulent energy dissipation rate as estimated from the data of the sonic anemometer at a height of 15 m (black curves) and the lidar at a height of 55 m (blue curves) and 100 m (red curves).The lidar estimates of the dissipation rate are obtained with Method 1 (a) and Method 2 (b).The lidar data are taken from Figure 4c,d.

2 )
. Red squares indicate the estimates of the dissipation rate from measurements by the sonic anemometer at a height of 15 m.In this experiment, the minimum height of measurement by the lidar is 0 h  55 m and the height step is h  15 m.To compare the estimates of the dissipation rate 2 () s h  (Method 2) and s  obtained from measurements, respectively, with the lidar and sonic anemometer at a height of s h = 15 m, we used the extrapolation of the height profile 2 () from the layer of 55-100 m to the height s h as follows.First, for each vertical profile 2 () k h shown in Figure7within the layer 55-100 m, we used a linear fitting of  is the slope and  is the intercept, to   of the sonic anemometer position.In Figure7, in the 15-100

Figure 6 .
Figure 6.Time series of the turbulent energy dissipation rate as estimated from the data of the sonic anemometer at a height of 15 m (black curves) and the lidar at a height of 55 m (blue curves) and 100 m (red curves).The lidar estimates of the dissipation rate are obtained with Method 1 (a) and Method 2 (b).The lidar data are taken from Figure 4c,d.
Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 22 m layer, the dashed lines show the profiles 2 () h  does not exceed 20%.Thus, if, in the 15-100 m layer, the vertical profile dependence which is close to linear, the lidar estimates of the dissipation rate obtained by Method 2 are in good agreement with the results of measurements with the sonic anemometer (at least, the discrepancy of the results are within the statistical error, examples of which are given below).

Figure 8 .
Figure 8.Time series of the mean wind velocity and the velocity ratio  (a), the turbulence energy dissipation rate (b) at a height of 42 m as obtained from measurements by the Stream Line lidar (black, green and blue curves) and the sonic anemometer (red curves).The measurements were carried out at the BEO starting from 08:00 LT on May 13 of 2020.For estimation of the dissipation rate, Method 1 (green curve) and Method 2 (blue curve) were used.

Figure 8 .
Figure 8.Time series of the mean wind velocity and the velocity ratio µ (a), the turbulence energy dissipation rate (b) at a height of 42 m as obtained from measurements by the Stream Line lidar (black, green and blue curves) and the sonic anemometer (red curves).The measurements were carried out at the BEO starting from 08:00 LT on 13 May of 2020.For estimation of the dissipation rate, Method 1 (green curve) and Method 2 (blue curve) were used.

1 
measurements by the sonic anemometer s  with the lidar estimates 1  and 2  .The results of the comparison are shown in Figure 11.It can be seen from Figure 11a that, at strong turbulence, (when s   10 −2 m 2 /s ) and 1   , the lidar estimates of the dissipation rate , on average, exceed the estimates s  obtained from measurements by the sonic anemometer.Under the same conditions, the difference between the estimates 2 and s  (Figure11b) is much smaller.

Figure 9 .
Figure 9.Comparison of the estimates of the mean wind velocity (a) and turbulent energy dissipation rate (b) obtained from simultaneous measurements by the lidar (vertical axes) and the sonic anemometer (horizontal axes).The data are taken from Figure 8.In (a) and (b), the brown line is the best fit given by equation i i Y aX b  , where is XU  (a) and lg( ) is X   (b).

Figure 9 .
Figure 9.Comparison of the estimates of the mean wind velocity (a) and turbulent energy dissipation rate (b) obtained from simultaneous measurements by the lidar (vertical axes) and the sonic anemometer (horizontal axes).The data are taken from Figure 8.In (a,b), the brown line is the best fit given by equation Y i = aX i + b, where X i = U s (a) and X i = lg(ε s ) (b).

Figure 10 .
Figure 10.Time series at a height of 42 m for (a) the ratio  of the mean wind velocity to the linear speed of motion of the probing volume, (b) the turbulence energy dissipation rate calculated from lidar data by Method 1 (green curve) and Method 2 (blue curve) and from data of sonic anemometer (red curve) and (c) the ratio of lidar estimates of the dissipation rate calculated by Methods 1 and 2. Measurements at the BEO started from 16:00 LT on May 20 of 2020.

Figure 10 .
Figure 10.Time series at a height of 42 m for (a) the ratio µ of the mean wind velocity to the linear speed of motion of the probing volume, (b) the turbulence energy dissipation rate calculated from lidar data by Method 1 (green curve) and Method 2 (blue curve) and from data of sonic anemometer curve) and (c) the ratio of lidar estimates of the dissipation rate calculated by Methods 1 and 2. Measurements at the BEO started from 16:00 LT on 20 May of 2020.

Figure 10 .
Figure 10.Time series at a height of 42 m for (a) the ratio  of the mean wind velocity to the linear speed of motion of the probing volume, (b) the turbulence energy dissipation rate calculated from lidar data by Method 1 (green curve) and Method 2 (blue curve) and from data of sonic anemometer (red curve) and (c) the ratio of lidar estimates of the dissipation rate calculated by Methods 1 and 2. Measurements at the BEO started from 16:00 LT on May 20 of 2020.

Figure 11 .
Figure 11.Comparison of estimates of the turbulence energy dissipation rate obtained from simultaneous measurements by the lidar (vertical axes) and the sonic anemometer (horizontal axes).Lidar estimates of the dissipation rate were obtained with Methods 1 (a) and 2 (b).The data are taken from Figure 10b.In (a) and (b), the brown line is the best fit given by equation i

Figure 11 .
Figure 11.Comparison of estimates of the turbulence energy dissipation rate obtained from simultaneous measurements by the lidar (vertical axes) and the sonic anemometer (horizontal axes).Lidar estimates of the dissipation rate were obtained with Methods 1 (a) and 2 (b).The data are taken from Figure 10b.In (a,b), the brown line is the best fit given by equation Y i = aX i + b, X i = lg(ε s ).