## Abstract

The accuracy of the method of azimuth structure function for estimation of the dissipation rate of the kinetic energy of turbulence from an array of radial velocities measured by low-energy micropulse coherent Doppler lidars with conical scanning by a probing beam around the vertical axis has been studied numerically. The applicability of the method in dependence on the turbulence intensity and the signal-to-noise ratio has been determined. The method of azimuth structure function was applied for estimation of the turbulent energy dissipation rate from radial velocities measured by the lidar in the experiments on the coast of Lake Baikal. Two dimensional time–height patterns of the wind turbulence energy dissipation rate were obtained. Part of them were obtained in presence of the atmospheric internal waves (AIWs) and low-level jet streams. It is observed that the wind turbulence in the area occupied by jet streams is very weak. In the process of dissipation of AIWs the wind turbulence strength increases.

© 2017 Optical Society of America

## 1. Introduction

The rate of dissipation of the kinetic energy of turbulence $\epsilon $ is one of the key parameters of the turbulent wind field. In the inertial range of spatial scales of inhomogeneities of the velocity field, the cascade process of turbulent energy transfer from large-scale inhomogeneities to small-scale ones is fully determined by this parameter and obeys the Kolmogorov 2/3 law [1]. The dissipation rate characterizes the intensity of turbulent processes in the atmosphere, and it is important to know how this parameter changes in space and time. Many papers are devoted to development of methods for remote and, in particular, lidar measurement of the dissipation rate (see monograph [2] and review [3] and references therein). There are also more recent publications concerning this issue [4–7]. Thus, in [6] it was proposed to estimate the dissipation rate within the inertial range of turbulence in the atmosphere from the transverse (azimuth) structure function of the radial velocity measured by a pulsed coherent Doppler wind lidar with conical scanning by a probing beam around the vertical axis. It was shown in [6] that estimates of the dissipation rate ${\widehat{\epsilon}}_{L}$ obtained by this method from data of the Halo Photonics Stream Line lidar and values of the dissipation rate ${\widehat{\epsilon}}_{S}$ measured simultaneously by a sonic anemometer at the same height are in good agreement. However, in the general case, the error of estimation of the dissipation rate by this method from wind data obtained with micropulse lidars of the class of Stream Line lidar under various atmospheric conditions was not analyzed yet.

The point is that fiber micropulse coherent Doppler lidars (CDL) similar Stream Line lidar are characterized by the very low signal-to-noise ratio at the background concentration of atmospheric aerosol. In [8], a method was developed for estimation of the wind speed and direction with acceptable accuracy from data measured by such lidars which can contain up to 50% of “bad” [9] estimates of the radial velocity. This method uses the procedure of filtering of good estimates of the radial velocity. There are no similar filtering procedures, which could be used in estimation of wind turbulence parameters. Therefore, for determination of turbulence parameters, only slightly noised data, in which the probability of bad estimate of the radial velocity is close to zero, can be used. As a result, the height limit for reconstruction of profiles of turbulent parameters from lidar data is much lower than the maximal height of reconstruction of wind speed and direction profiles, and it is necessary to know previously which data can be used for estimation of turbulence parameters.

Below, in the next part of this paper, we analyze the accuracy of the method of azimuth structure function (MASF) in dependence on the atmospheric conditions: concentration of aerosol (signal-to-noise ratio) and turbulence intensity (dissipation rate). In third section of the paper, we present the results of lidar studies of wind turbulence in the atmospheric boundary layer in the shore zone of Lake Baikal with this method.

## 2. Theoretical foundation

#### 2.1 Basic equations

During the measurement of the radial velocity by a lidar with the use of conical scanning by the probing beam around the vertical axis, the beam elevation angle $\phi $ is fixed, and the azimuth angle $\theta $ with time $t$ varies from 0° to 360° with the angular scanning rate ${\omega}_{c}$. The raw data of measurements by the Stream Line lidar are arrays of estimates of correlation functions of the complex signal $\widehat{C}(l{T}_{s};{R}_{k},{\theta}_{m},n)$, where $l=0,\text{\hspace{0.17em}}1,\text{\hspace{0.17em}}\mathrm{...},\text{\hspace{0.17em}}6$; ${T}_{s}=1/B$; $B$ = 50 MHz is the receiver bandwidth; ${R}_{k}={R}_{0}+k\text{\hspace{0.17em}}\Delta R$ is the distance from the lidar to the center of the sensing volume; $\Delta R$ is the range step, which is usually taken equal to 18 or 30 m; $k=0,\text{\hspace{0.17em}}1,\text{\hspace{0.17em}}2,\text{\hspace{0.17em}}\mathrm{...},\text{\hspace{0.17em}}K$; ${\theta}_{m}=m\Delta \theta $ is the azimuth angle; $m=1,\text{\hspace{0.17em}}2,\text{\hspace{0.17em}}3,\text{\hspace{0.17em}}\mathrm{...},\text{\hspace{0.17em}}M$; $\Delta \theta $ is the resolution in the azimuth angle, and $n=1,\text{\hspace{0.17em}}2,\text{\hspace{0.17em}}3,\text{\hspace{0.17em}}\mathrm{...},\text{\hspace{0.17em}}N$ is the scan number. For every correlation function $\widehat{C}(l{T}_{s};{R}_{k},{\theta}_{m},n)$, the accumulation of data with the same number ${N}_{a}$ of probing pulses is used. In the case of the Stream Line lidar, the pulse repetition frequency is${f}_{P}$ = 15 kHz. If ${N}_{a}$ = 1500 (the minimal possible value ${N}_{a}$ for this lidar), the time for measurement of the correlation function of the complex lidar signal is ${T}_{C}={N}_{a}/{f}_{P}$ = 0.1 s. Then, of the scanning rate ${\omega}_{c}$ = 10°/s, the resolution in the azimuth angle is $\Delta \theta ={\omega}_{c}{T}_{C}$ = 1°, the number of rays for one full scan is $M$ = 360, and the duration of one scan is ${T}_{\text{scan}}={T}_{C}M$ = 36 s.

For estimating turbulence parameters from wind data, it is important to conduct measurements with the high resolution in the azimuth angle. At the same time, the tangential velocity of motion of the sensing volume ${V}_{t}={\omega}_{c}{R}_{k}\mathrm{cos}\phi $ should exceed significantly the average wind velocity at the height of measurement.

The estimates of the correlation functions $\widehat{C}(l{T}_{s};{R}_{k},{\theta}_{m},n)$ are used to calculate the Doppler spectra ${\widehat{S}}_{D}(f;{R}_{k},{\theta}_{m},n)$ [10], and the positions of their maximum are then used to obtain the array of estimates of the radial velocities ${\widehat{V}}_{r}({R}_{k},{\theta}_{m},n)$. In addition, the array of estimates of the signal-to-noise ratio is determined from the array $\widehat{C}(l{T}_{s};{R}_{k},{\theta}_{m},n)$ as $\text{S}\widehat{N}\text{R}({R}_{k},{\theta}_{m},n)=\widehat{C}(0;{R}_{k},{\theta}_{m},n)/{P}_{\text{N}}-1$, where ${P}_{\text{N}}$ is the noise power. The further procedure of processing of lidar data in order to obtain estimates of the dissipation rate $\widehat{\epsilon}({h}_{k})$ at heights ${h}_{k}={R}_{k}\mathrm{sin}\phi $ consists in the following.

For every height ${h}_{k}$, the vector of average wind velocity $V=\{{V}_{z},{V}_{x},{V}_{y}\}$ is determined from the array of estimates of radial velocities ${\widehat{V}}_{r}({R}_{k},{\theta}_{m},n)$ with the aid of sine-wave fitting [2], and fluctuations of the radial velocity are calculated as

In Eq. (2), ${m}^{\prime}=1,\text{\hspace{0.17em}}2,\text{\hspace{0.17em}}\mathrm{...},\text{\hspace{0.17em}}{M}^{\prime}$ and ${M}^{\prime}\Delta \theta <<\pi /2$ (here $\Delta \theta $ is measured in radians). The condition ${M}^{\prime}\Delta \theta <<\pi /2$ allows the azimuth structure function ${D}_{L}({m}^{\prime}\Delta \theta )$ to be replaced with the transverse [11] structure function of the radial velocity ${D}_{\perp}({m}^{\prime}\Delta {y}_{k})$, where $\Delta {y}_{k}=\Delta \theta {R}_{k}\mathrm{cos}\phi $ is transverse (with respect to the probing beam optical axis) size of the sensing volume. The estimate of the turbulence energy dissipation rate can be found from the difference of lidar estimates of the azimuth structure function $\Delta {\widehat{D}}_{L}({M}^{\prime}\Delta \theta )={\widehat{D}}_{L}({M}^{\prime}\Delta \theta )-{\widehat{D}}_{L}(\Delta \theta )$ [6], 1) if ${M}^{\prime}\Delta y$ does not exceed the upper boundary of the inertial range of turbulence, 2) the probability of bad estimate [9] of the radial velocity ${P}_{b}$ is zero, and 3) ${\widehat{V}}_{r}({R}_{k},{\theta}_{m},n)$ can be represented in the form [2]

The probability of bad estimate of the radial velocity from lidar data depends on atmospheric conditions, mostly, on the concentration of aerosol particles, which manifests itself from the signal-to-noise ratio. Thus, to determine the conditions of applicability of MASF for estimation of the dissipation rate, it is necessary to find the minimal value of the signal-to-noise ratio $\text{SNR}$, at which the probability of bad estimate ${P}_{b}$ in the array of estimates ${\widehat{V}}_{r}({R}_{k},{\theta}_{m},n)$ obtained from measurements micropulse lidar is close to zero (does not exceed an arbitrarily small preset value). Toward this end, we should know the dependence ${P}_{b}(\text{SNR})$. This dependence can be found on the basis of the theory developed in Ref [12], or using numerical simulation of random realizations of lidar estimates of radial velocity.

Using the results of numerical simulation of lidar signals at different values of signal-to-noise ratio $\text{SNR}$(description of the simulation algorithm is given in [2,9]), we found the dependence ${P}_{b}(\text{SNR})$ for the Stream Line lidar with the following parameters: laser wavelength $\lambda $ = 1.5 μm, probing pulse duration ${\tau}_{p}$ = 170 ns, and the width of the time window used to estimate the Doppler spectra ${\tau}_{w}$ = 120 ns. In simulation we set that the noise component of the lidar signal is a white noise with Gaussian probability density function (PDF) and neglected by homogeneity of the radial velocity within the sensing volume. It is valid in the case of weak wind turbulence and absence of a strong wind shear. For spectral accumulation we used ${N}_{a}$ = 1500 laser shots. For each SNR value, we obtained 10^{6} statistically independent radial velocity estimates ${\widehat{V}}_{r}$. From the obtained array of simulated radial velocity estimates, we calculated the PDF of the radial velocity estimate $p({\widehat{V}}_{r})$. Then we determined the probability (or fraction) ${P}_{b}$ of a bad estimate of the radial velocity by minimizing the mean square difference between the PDF $p({\widehat{V}}_{r})$ obtained in the numerical experiment and the model PDF ${p}_{M}({\widehat{V}}_{r})$ described by Eq. (34) in Ref [9].

From the obtained dependence ${P}_{b}(\text{SNR})$at ${N}_{a}$ = 1500 it follows that in the case of $M$ = 360 and $N$ = 20 the probability ${P}_{b}$ is less than ${(MN)}^{-1}$, if $\text{SNR}\ge $0.03 (−15.2 dB). We also obtained dependencies ${P}_{b}(\text{SNR})$ for other values of accumulation number ${N}_{a}$. Analysis of the obtained results has shown that the probability of a bad estimate of the radial velocity is a function of the product of $\text{SNR}$and $\sqrt{{N}_{a}}$, at least for ${N}_{a}\ge $1000. Therefore, if for ${N}_{a}$ = 1500 (${T}_{C}$ = 0.1 s) the condition ${P}_{b}<{(MN)}^{-1}$ ($M$ = 360 and $N$ = 20) is fulfilled, when $\text{SNR}\ge $0.03 (−15.2 dB), then for ${N}_{a}$ = 15000 (${T}_{C}$ = 1 s) this condition is realized at $\text{SNR}\ge 0.03/\sqrt{10}$ (−20.2 dB).

Thus, obtained relations determine the measurement conditions at which the azimuth structure function of the radial velocity can be calculated from the Stream Line lidar array of radial velocity estimates without using a special data filtration procedure.

#### 2.2 Azimuth structure function of the radial velocity

The condition ${M}^{\prime}\Delta {y}_{k}\le {L}_{I}$, where ${L}_{I}$ is the outer scale of the inertial range, allows the turbulence to be considered as isotropic in the calculation of the azimuth structure function of the radial velocity ${D}_{L}({m}^{\prime}\Delta \theta )=\text{\hspace{0.17em}}<{\widehat{D}}_{L}({m}^{\prime}\Delta \theta )>$. Then, assuming horizontal statistical homogeneity and stationarity of the turbulent wind field, with allowance for Eq. (3), for ${D}_{L}({m}^{\prime}\Delta \theta )$ we obtain

where $\overline{D}({m}^{\prime}\Delta \theta )=\text{\hspace{0.17em}}<{[{{\overline{V}}^{\prime}}_{r}({R}_{k},{\theta}_{m}+{m}^{\prime}\Delta \theta ,n)-{{\overline{V}}^{\prime}}_{r}({R}_{k},{\theta}_{m},n)]}^{2}>$ is the azimuth structure function of the radial velocity averaged over the sensing volume and ${{\overline{V}}^{\prime}}_{r}={\overline{V}}_{r}-\text{\hspace{0.17em}}<{V}_{r}>$. For calculation of the structure function $\overline{D}({m}^{\prime}\Delta \theta )$ in Eq. (4) at ${\sigma}_{e}^{2}$ = 0, we used method Monte Carlo based on the algorithm of computer simulation of random realizations of two-dimensional distributions of the radial velocity ${V}_{r}(R+{z}^{\prime},\theta )$ [2,13].Figure 1 shows an example of simulation of random realizations of the dependence on the azimuth angle ${\theta}_{m}$ for the radial velocity ${V}_{r}(R,{\theta}_{m})$ (blue curve) and the lidar estimate of the radial velocity from positions of maxima of the Doppler spectra ${\widehat{V}}_{r}(R,{\theta}_{m})$ (red curve). In simulation we set the average wind speed 10 m/s, the elevation angle $\phi $ = 60°, the distance to the center of the sensing volume $R$ = 344 m (height $h\approx $ 300 m), resolution in the azimuth angle $\Delta \theta $ = 1° ($\Delta {y}_{k}$ = 3 m), the turbulence energy dissipation rate $\epsilon $ = 10^{−2} m^{2}/s^{3}, and the integral scale of correlation of wind velocity fluctuations ${L}_{V}$ = 250 m. Since ${\sigma}_{e}^{2}$ = 0, according to Eq. (3) the lidar estimate of the radial velocity is ${\widehat{V}}_{r}(R,{\theta}_{m})={\overline{V}}_{r}(R,{\theta}_{m})$. One can see that the red curve in Fig. 1 is smoother in comparison with the blue curve, that is, the averaging over the sensing volume takes place. In the case of the Stream Line lidar, the longitudinal dimension of the sensing volume is $\Delta z$ = 30 m.

According to Eq. (4), the difference structure function $\Delta {D}_{L}({m}^{\prime}\Delta \theta )$ = ${D}_{L}({m}^{\prime}\Delta \theta )-{D}_{L}(\Delta \theta )$ is independent of the instrumental error of the estimate of radial velocity ${\sigma}_{e}$ and is a difference structure function averaged over sensing volume $\Delta \overline{D}(y)=\overline{D}({m}^{\prime}\Delta \theta )-\overline{D}(\Delta \theta )$, where $y={m}^{\prime}\Delta {y}_{k}$. Figure 2 demonstrates the effect of averaging over the sensing volume on the difference structure function.

Figure 2 shows the calculated by Monte Carlo method difference structure functions $\Delta D(y)=D({m}^{\prime}\Delta \theta )-D(\Delta \theta )$, where $D({m}^{\prime}\Delta \theta )=\text{\hspace{0.17em}}<{[{{V}^{\prime}}_{r}(R,{\theta}_{m}+{m}^{\prime}\Delta \theta )-{{V}^{\prime}}_{r}(R,{\theta}_{m})]}^{2}>$ is the structure function of the radial velocity (without averaging over the sensing volume), and $\Delta {\overline{D}}_{\mathrm{max}}(y)$, and $\Delta {\overline{D}}_{\text{cnt}}(y)$, where the subscripts “$\mathrm{max}$” and “$\text{cnt}$” indicate that the radial velocity was estimated, respectively, from the position of the spectral maximum and from the centroid of the spectral distribution. It can be seen that, due to the averaging of the radial velocity over the sensing volume, $\Delta {\overline{D}}_{\text{cnt}}(y)$ and $\Delta {\overline{D}}_{\mathrm{max}}(y)$ are less than $\Delta D(y)$ in the considered range of $y$ variation. Therefore, when estimating the dissipation rate from the lidar data, one should necessarily take into account the averaging over the sensing volume.

#### 2.3 Relative error of lidar estimate of the dissipation rate

The relative error of lidar estimate of the dissipation rate ${E}_{\epsilon}=\sqrt{<{(\widehat{\epsilon}/\epsilon -1)}^{2}>}$ was calculated for ${N}_{a}$ = 1500, $N$ = 20, $M$ = 360, $\Delta {y}_{k}$ = 3 m, and ${M}^{\prime}$ = 7 ($y$ = 21 m) in dependence the signal-to-noise ratio $\text{SNR}$ and magnitude of the dissipation rate $\epsilon $.

For the ratio of the estimate of dissipation rate $\widehat{\epsilon}$ to the given (true) value $\epsilon $, the following is valid: $\widehat{\epsilon}/\epsilon ={[\Delta {\widehat{D}}_{L}(y)/\Delta \overline{D}(y)]}^{3/2}$, where $\Delta {\widehat{D}}_{L}(y)={\widehat{D}}_{L}(y)-{\widehat{D}}_{L}(\Delta {y}_{k})$. We assume that ${E}_{\epsilon}<<1$. Then the error ${E}_{\epsilon}$ can be represented in the form

The value of ${\sigma}_{\overline{D}}$ was calculated from the model data of a pulsed lidar with the use of numerical simulation of random realizations ${V}_{r}(R+{z}^{\prime},\theta )$ and estimation of the radial velocity from the positions of maxima of the Doppler spectra. At $N$ = 20, $M$ = 360, $\Delta {y}_{k}$ = 3 m, and ${M}^{\prime}$ = 7, we have obtained ${\sigma}_{\overline{D}}\text{\hspace{0.17em}}$ = 0.04. This value is independent of the dissipation rate $\epsilon $. The model difference transverse structure function averaged over sensing volume $\Delta \overline{D}(y)$ in accordance with the theory [13] is represented in a form

whereTo calculate the instrumental error of estimation of the radial velocity ${\sigma}_{e}$ as a function of $\text{SNR}$, we used the data of numerical simulation of the lidar signal based on the approach [2,8,9,13]. The simulation was performed under assumption of homogenous wind without taking into account the turbulent fluctuations of wind velocity. It follows from the obtained numerical results that the dependence ${\sigma}_{e}(\text{SNR})$ for the Stream Line lidar parameters (${\tau}_{p}$ = 170 ns, ${\tau}_{w}$ = 120 ns and the sampling interval of the complex lidar signal 20 ns) can be approximated well by the simple formula:

where $\Delta v$ = 0.4 m/s is an empirical parameter of fitting.Rough estimates demonstrate that for such small values of the pulse duration and the time window, the turbulent broadening of the Doppler spectrum is negligible as compared to the instrumental broadening, at least in the case of weak and moderate wind turbulence ($\epsilon $< 0.001 m^{2}/s^{3}). Thus, at least in the range of the turbulent intensity determined by the condition $\epsilon $< 0.001 m^{2}/s^{3}, Eq. (12) can be used for estimation of ${\sigma}_{e}$.

Figure 3 shows the calculated dependence of the relative error in the lidar estimate of dissipation rate of the kinetic energy of turbulence ${E}_{\epsilon}\times 100\%$ on the signal-to-noise ratio at the different intensity of wind turbulence. It can be seen that the dissipation rate can be estimated with the acceptable accuracy from measurements by the Stream Line lidar even at the very weak turbulence, when $\epsilon $ = 10^{−6} m^{2}/s^{3}, if $\text{SNR}$ is no less than 0.24.

The data of numerical simulation were used to determine, by Eqs. (6), (7), and (12) the values of $\epsilon $ and $\text{SNR}$, at which ${\sigma}_{I}\le 0.1$. With allowance for ${\sigma}_{\overline{D}}$ = 0.04, the relative error ${E}_{\epsilon}\times 100\%$ in this case should not exceed 16%. Figure 4 shows the dependence of the signal-to-noise ratio on the turbulence energy dissipation rate, at which the relative error ${E}_{\epsilon}\times 100\%$ = 16%.

In Fig. 4, the grey color indicate an area of $\epsilon $ and $\text{SNR}$ values, at which the error of lidar estimate of dissipation rate exceed 16%. The dashed line in the figure shows the level $\text{SNR}$ = 0.03, below which the estimation of the dissipation rate is impossible at any value of $\epsilon $ because of the presence of bad estimates in the measured array of radial velocities. Figure 4 allows us to judge the representativeness of estimation of the dissipation rate $\widehat{\epsilon}$ from the measured signal-to-noise ratio $\text{S}\widehat{N}\text{R}$.

## 3. Experiment

Wind turbulence in the boundary layer of atmosphere was studied on the coast of the Lake Baikal in summer campaigns of 2015 and 2016 with the use of a Stream Line lidar. The lidar was mounted 340 m away from the coast of the lake, on the territory of the Baikal Astrophysical Observatory of Institute of Solar-Terrestrial Physics, Siberian Branch, Russian Academy of Sciences (51°50′47.17”N, 104°53′31.21”E); the lidar altitude above the lake level was 180 m. Using the initial data of these measurements (arrays of radial velocities and signal-to-noise ratios) we retrieved the time–height patterns of the turbulence kinetic energy dissipation rate and calculated the relative error of assessment of the dissipation rate from the lidar data.

#### 3.1 Methodology

Stream Line lidar measurements were carried out under the following parameters: elevation angle $\phi $ = 60°; number of pulses used for data accumulation ${N}_{a}$ = 3000; angular speed of probing beam scanning ${\omega}_{c}$ = 10°/s; measurement time per scanning 36 s; resolution in the azimuth angle $\Delta \theta $ = 2°, and the number of rays per scanning $M$ = 180.

To retrieve the vertical profiles of the wind speed *U*, the wind direction angle ${\theta}_{V}$, and the vertical component of the wind velocity ${V}_{z}$, the method of filtered sine wave fitting was used [2,8]. As shown in [14] for ${N}_{a}$ = 3000, the minimal signal to noise ratio is 0.005 (–23 dB), under which the wind speed estimation error does not exceed 0.1 m/s.

To assess the turbulence energy dissipation rate $\epsilon $, we used the method of transverse (azimuth) structure function of the radial velocity measured within the turbulence inertial interval by the scanning lidar. The time–height patterns of the turbulence energy dissipation rate $\epsilon (h,t)$ were calculated with a height step of 26 m and a time step of 36 s based on the array of radial velocities ${\widehat{V}}_{r}({R}_{k},{\theta}_{m},n)$ measured with the use of data received for $N$ = 25 scannings (measurement time is 15 min) for each vertical profile of $\epsilon $. For each height *h*, the SNR values from the array of $\text{S}\widehat{N}\text{R}({R}_{k},{\theta}_{m},n)$ estimates were averaged over all the azimuth angles ${\theta}_{m}$ and over $N$ = 25 scannings. Then the relative error of the dissipation rate ${E}_{\epsilon}$ was calculated by Eqs. (6) and (7). The error of the radial velocity estimate ${\sigma}_{e}$, which enters into Eq. (7), was determined from the experimental data by two methods: 1) calculation of ${\sigma}_{e}$ by Eq. (12), using measured SNR and 2) estimation of ${\sigma}_{e}$, using Eq. (22) given in [15]. The second method is applicable for arbitrary intensity of the wind turbulence.

#### 3.2 Results

The analysis of the lidar data processing results revealed several AIW events occurred in the presence of jet streams. The measurements on August 23, 2015, showed an AIW with a period of 9 min and longitudinal and vertical components of the velocity of about 1 m/s and 0.3 m/s, respectively, which was propagating through the region of experiment for almost 6 h in the presence of two jet streams. The initial data measured by the lidar from 11:00 to 17:00 on August 23, 2015, satisfied the condition $\text{SNR}\text{\hspace{0.17em}}\text{>}$ 0.05 up to 500 m. We limited to this altitude in the retrieval of the vertical profile of the dissipation rate, while the profiles of the wind velocity components were retrieved up to height of about 900 m.

Figure 5 shows the time–height patterns of the wind speed $U(h,t)$, vertical component of the wind velocity ${V}_{z}(h,t)$, $\epsilon (h,t)$, $\text{SNR}(h,t)$, ${\sigma}_{e}(h,t)$, and ${E}_{\epsilon}(h,t)$. The SNR is below a threshold of 0.005 in the black regions. The values for the error in estimating the radial velocity ${\sigma}_{e}(h,t)$ presented in Fig. 5 were obtained by the second method (using Eq. (22) in [15]), since in many cases the dissipation rate $\epsilon (h,t)$ exceeded 10^{−3} m^{2}/s^{3}. These values were used to determine ${E}_{\epsilon}(h,t)$. SNR varied from 0.05 to 0.2 in the 100–500 m layer during that measurement. The radial velocity error does not exceed 0.16 m/s. The data on the wind speed in Fig. 5 clearly show the presence of two jet streams with speed maxima at altitudes of about 250 m and 750 m, respectively (hereinafter, the altitudes are given relative to the lidar location point). The directions of these jet streams are almost perpendicular to each other (the bottom jet stream was directed from north to south through mountains, and the top jet stream was directed from east to west, i.e., from the lake side) [14]. Near-harmonic oscillations of $U(h,t)$ and ${V}_{z}(h,t)$ are clearly seen in Fig. 5 from 14:20 to 15:00. The data on the turbulence energy dissipation rate in Fig. 5 show that $\epsilon (h,t)$ takes values from a wide range, attaining a maximum of about 0.006 m^{2}/s^{3} before the arising of a jet stream. Relatively strong wind turbulence is also observed under the jet stream. The turbulence is quite weak inside the jet stream, including during an AIW event.

Figures 6 and 7 show the vertical profiles of $\epsilon $, $\text{SNR}$, ${E}_{\epsilon}$, and ${\sigma}_{e}$ measured with an interval of 1 h (the data are taken from Fig. 5).The turbulence energy dissipation rate obviously can be retrieved from lidar measurements with an error of no larger than 30% at $\epsilon ~5\cdot {10}^{-6}$m^{2}/s^{3}. The SNR should be no less than 0.075 in this case. The turbulence spatial structure should include an inertial interval with the upper boundary of no less than a half of the longitudinal size of the sensing volume $\Delta z$ ($\Delta z$ = 30 m for the Stream Line lidar).

In Figs. 6 and 7, the dashed curves show the results for ${\sigma}_{e}$ and ${E}_{\epsilon}$ calculated by Eq. (12) and Eq. (6) (the first method). Using all the data for ${\sigma}_{e}(h,t)={\sigma}_{2}$ given in Fig. 5 (the second method) and the data for ${\sigma}_{e}(h,t)={\sigma}_{1}$ found by the first method, we calculated values $b=\text{\hspace{0.17em}}<\eta >$ and $d=\sqrt{<{\eta}^{2}>}$, where $\eta =2({\sigma}_{1}-{\sigma}_{2})/({\sigma}_{1}+{\sigma}_{2})$ and angular brackets mean spatial-temporal averaging over all obtained data for ${\sigma}_{1}$ and ${\sigma}_{2}$. As a result, we found that $b$ = −0.07 and $d$ = 0.13. It means that the errors of radial velocity estimate (and, correspondingly, the errors of the turbulent energy dissipation rate estimate) calculated by these two methods differ insignificantly.

In August of 2016, in contrast to 2015, when the concentration of aerosol particles was increased due to forest fires, the measurements were conducted at the rather weak lidar echo signal. Long rains washing out the aerosol from the atmosphere were often falling. Therefore, we succeeded in conducting adequate measurements in the absence of rain only for three days. To increase the signal-to-noise ratio at the height of localization of the low-level jet flow, we often had to refocus the probing beam to the closer distance $F$ = 300 m. The elevation angle $\phi $, as in 2015, was taken equal to 60°, the scanning rate was 5°/s (duration of one conical scan of 72 s), and ${N}_{a}$ = 3000 ($\Delta \theta $ = 1°, $M$ = 360). Despite unfavorable weather conditions, we still succeeded in revealing one case of appearance of an atmospheric internal wave from the lidar data.

Figure 8 depicts the two-dimensional distributions of the wind speed $U(h,t)$, wind direction angle ${\theta}_{V}(h,t)$, and vertical component of the wind vector ${V}_{z}(h,t)$ obtained from raw data of the Stream Line lidar as measured on August 6 of 2016. During the measurements, slightly overcast weather was observed approximately to 01:00 Local Time. Nevertheless, a jet flow was observed, in which the height of maximum of the wind velocity decreased from 350 m to 200 m with time. After 5:30 LT, the jet flow practically disappeared, and wind oscillations took place, which especially clearly seen in the figures for the wind direction angle and the vertical component. We believe that these oscillations are associated with appearance of an atmospheric internal wave early in the morning, when the stratification of air temperature is most stable.

Since the measurements were conducted under conditions of weak echo signal, it appeared impossible to obtain two-dimensional distributions of the dissipation rate $\epsilon (h,t)$. Only in the thin layer of 250–350 m, owing to the focusing of the probing beam to 300 m, the signal-to-noise ratio met the requirement of very low probability of bad estimate of the radial velocity. To obtain estimates of the dissipation rate in this layer, we used the array of radial velocities measured by the lidar for 25 scans (measurement duration of 30 min).

Figure 9 shows the dynamics of the wind speed, wind direction angle, vertical component of the wind velocity vector, turbulence energy dissipation rate, and relative error of lidar estimate of the dissipation rate at a height of 300 m (480 m above the surface of Lake Baikal). Every point in Fig. 9(a), Fig. 9(b), Fig. 9(c) corresponds to the measurement for 72 s (one full scan). The error of estimation of the dissipation rate ${E}_{\epsilon}$ was determined as described above in section 2.3. At this height till approximately 07:30 LT, the wind direction was from the shore to the lake. Then, for the relatively short time, the wind direction alternated by approximately 90°.

As can be seen from Fig. 9, since 01:30 LT for five hours at a height of 300 m the vertical component of wind ${V}_{z}$ Fig. 9(c) took only negative values (downward flow). Starting from 04:40 LT, the time profile of ${V}_{z}$ began to oscillate. The period of these oscillations and their amplitude decreased, on average, with time. According to Fig. 9(c), since 05:00 till 07:30 LT the amplitude of oscillations of the vertical velocity was approximately 0.6 m/s. At the same time, the amplitude of oscillations of the horizontal velocity $U$Fig. 9(a) was, at least, two times smaller. This behavior of the internal wave differs significantly from the results of the experiment in 2015 [14], when we observed nearly harmonic oscillations of the three components of the wind velocity vector with the amplitudes, on average, of 1 m/s for the horizontal speed and 0.3 m/s for the vertical one.

We succeeded in obtaining estimates for the turbulence energy dissipation rate with an acceptable accuracy only for time interval, whose total duration is about 30% of the total measurement time (see Fig. 9(e), Fig. 9(f)). The data of Fig. 9 demonstrate that the turbulence inside the jet flow is very weak. It is very similar that the inertial range of turbulence is completely absent, and the flow in the jet is nearly laminar, low frequency wind variations in the jet are mostly caused by mesoscale processes in the atmosphere. We can also suppose that observed increase in the turbulence intensity after 07:30 LT is caused by decay of the internal wave and transfer of its energy to turbulization of the wind flow.

## 4. Conclusions

The paper has proposed the method for calculation of the relative error in lidar estimate of the dissipation rate of the kinetic energy of turbulence. The accuracy of the method of azimuth structure function for estimation of the dissipation rate from the array of radial velocities measured by micropulse coherent Doppler lidars at the conical scanning by the probing beam around the vertical axis has been analyzed numerically. It was shown that the validity of this method depends on the magnitude of the dissipation rate measured and the signal to noise ratio. Obtained in the paper quantitative results allow one to determine beforehand the representativeness of lidar estimation of the turbulence energy dissipation rate from the measured signal to noise ratios. It is especially important when using the low energy micropulse CDL such as Halo Photonics commercial lidar system Stream Line.

The method of azimuth structure function was applied for estimation of the turbulent energy dissipation rate from radial velocities measured by the scanning micropulse CDL Stream Line in the boundary layer of atmosphere on the coast of Lake Baikal. Two dimensional time–height patterns of the wind turbulence energy dissipation rate were obtained and accuracy of lidar estimate of the dissipation rate was determined. Few events of the internal atmospheric waves were registered during the experimental campaigns. All the registered AIWs were observed in the presence of low level jet streams. It is observed that the wind turbulence in the area occupied by jet streams is very weak. Possibly, it is absent at all, and wind velocity variations in the area of jet streams are caused by mesoscale atmospheric processes. In the process of dissipation of AIWs the wind turbulence strength increases. These results are important for understanding the turbulence - wave interaction processes in the atmosphere and evaluation of wind energy variations [7,16].

## Funding

Russian Science Foundation (No. 14-17-00386-П).

## References and links

**1. **A. N. Kolmogorov, “Local structure of turbulence in incompressible viscous fluid at very large Reynolds numbers,” Doklady AN SSSR **30**, 299–303 (1941).

**2. **V. Banakh and I. Smalikho, *Coherent Doppler Wind Lidars in a Turbulent Atmosphere* (Artech House, 2013).

**3. **A. Sathe and J. Mann, “A review of turbulence measurements using ground-based wind lidars,” Atmos. Meas. Tech. **6**(11), 3147–3167 (2013). [CrossRef]

**4. **A. Sathe, R. Banta, L. Pauscher, K. Vogstad, D. Schlipf, and S. Wylie, “Estimating Turbulence Statistics and Parameters from Groundand Nacelle-Based Lidar Measurements:IEA Wind Expert Report,” (2015).

**5. **A. Sathe, J. Mann, N. Vasiljevic, and G. Lea, “A six-beam method to measure turbulence statistics using ground-based wind lidars,” Atmos. Meas. Tech. **8**(2), 729–740 (2015). [CrossRef]

**6. **I. N. Smalikho, V. A. Banakh, A. V. Falits, and Yu. A. Rudi, “Estimation of the turbulent energy dissipation rate from the data measured by the “StreamLine” lidar in the surface layer of atmosphere,” Atmos. Oceanic Opt. **28**(10), 901–905 (2015).

**7. **S. Wu, B. Liu, J. Liu, X. Zhai, C. Feng, G. Wang, H. Zhang, J. Yin, X. Wang, R. Li, and D. Gallacher, “Wind turbine wake visualization and characteristics analysis by Doppler lidar,” Opt. Express **24**(10), A762–A780 (2016). [CrossRef] [PubMed]

**8. **V. A. Banakh and I. N. Smalikho, “Wind sensing in an atmospheric boundary layer by means of micropulse coherent Doppler lidars,” Opt. Spectrosc. **121**(1), 152–159 (2016). [CrossRef]

**9. **R. G. Frehlich and M. J. Yadlowsky, “Performance of mean-frequency estimators for Doppler radar and lidar,” J. Atmos. Ocean. Technol. **11**(5), 1217–1230 (1994). [CrossRef]

**10. **V. A. Banakh, V. A. Brewer, E. L. Pichugina, and I. N. Smalikho, “Measurements of wind velocity and direction with coherent Doppler lidar in conditions of a weak echo signal,” Atmos. Oceanic Opt. **23**(5), 381–388 (2010). [CrossRef]

**11. **N. L. Byzova, V. N. Ivanov, and E. K. Garger, *Turbulence in Atmospheric Boundary Layer* (Gidrometeoizdat, 1989) [in Russian]

**12. **B. J. Rye and R. M. Hardesty, “Discrete spectral peak estimation in incoherent backscatter heterodyne lidar. I: Spectral accumulation and Cramer-Rao lower bound,” IEEE Trans. Geosci. Remote Sens. **31**(1), 16–27 (1993). [CrossRef]

**13. **I. N. Smalikho and V. A. Banakh, “Accuracy of estimation of the turbulent energy dissipation rate from wind measurements with a conically scanning pulsed coherent Doppler lidar,” Part I. Algorithm of data processing, Atmos. Oceanic Opt. **26**(5), 404–410 (2013). [CrossRef]

**14. **V. A. Banakh and I. N. Smalikho, “Lidar observations of atmospheric internal waves in the boundary layer of atmosphere on the coast of Lake Baikal,” Atmos. Meas. Tech. **9**(10), 5239–5248 (2016). [CrossRef]

**15. **I. N. Smalikho and V. A. Banakh, “Measurements of wind turbulence parameters by a conically scanning coherent Doppler lidar in the atmospheric boundary layer,” Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2017-140 (2017). [CrossRef]

**16. **I. N. Smalikho, V. A. Banakh, Y. L. Pichugina, W. A. Brewer, R. M. Banta, J. K. Lundquist, and N. D. Kelley, “Lidar investigation of atmosphere effect on a wind turbine wake,” J. Atmos. Ocean. Technol. **30**(11), 2554–2570 (2013). [CrossRef]