3.1. Experiment at the Basic Experimental Observatory of the V.E. Zuev Institute of Atmospheric Optics SB RAS
The lidar experiment at the territory of Basic Experimental Observatory (latitude 56.475448°N, longitude 85.048115°E) of the V.E. Zuev Institute of Atmospheric Optics (IAO) SB RAS in Tomsk was carried out in July 2016. The continuous measurements by the StreamLine lidar with the conical scanning by the probing beam around the vertical axis at the elevation angle
= 35.3° were conducted from 19:00 Local Time (LT) on 20 July 2016 to 14:00 LT on 24 July 2016. The duration of one conical scan was
= 1 min, the time for accumulation of raw lidar data was
= 0.5 s, the number of rays for one complete scan was
= 120, and the azimuth angle resolution was
= 3°. To obtain estimates of wind turbulence parameters, we used the data measured by the lidar for
= 30 scans (30 min). The range gate length
was set equal to 18 m, and the step in height was
≈ 10 m. For the more detailed description of the experiment, see Ref. [
15].
Figure 1 shows the spatiotemporal distributions of the kinetic energy of turbulence, the dissipation rate of turbulent energy, and the integral scale of turbulence obtained for three days: July 21–23. One can clearly see the contrast between the turbulence intensity in daytime and night-time. At night, the turbulence is very weak due to the stability of thermal stratification. In the daytime, intense turbulent mixing occurs. The data for the turbulence kinetic energy and the dissipation rate calculated by Equations (4)–(6) do not contradict the theory of the atmospheric boundary layer [
18,
19,
24], and are in a good agreement with the results of measurements by the sonic anemometer [
15]. The values of the integral scale of turbulence obtained with Equation (8) for the daytime, when the temperature stratification is neutral or unstable, are also in agreement with the known experimental data [
1]. However, in the night-time under stable atmospheric conditions, obtained estimates of
have very large values, which certainly contradicts the theory. Areas of improbable estimates of
are colored in white in
Figure 1d. Correspondingly, the correctness of quantitative estimates of
and
in the night-time also raises doubts. It is worth stressing that during this experiment, the signal-to-noise ratio (SNR) was high enough in order that the error of estimation of turbulent parameters from lidar data was small, even for very weak turbulence.
Next, we consider the results of night-time lidar measurements in more detail. In the period of 20–24 July 2016, LLJ occurred three times: (1) from 22:00 on July 20 to 05:00 on July 21, (2) from 21:00 on July 21 to 06:00 on July 22, and (3) from 22:00 on July 22 to 04:00 on July 23 (hereinafter, local time is used). On the night of July 23–24, the jet stream was not observed.
Figure 2 illustrates the variation of the wind speed and direction with height. Each plot of this figure shows the vertical profiles retrieved from lidar measurements for one scan (1 min) every hour. Taking into account that these measurements were conducted at strongly stable temperature stratification, we can assume that variations of wind profiles for one hour are more likely caused by the nonstationary (mesoscale process) than by small-scale turbulence. Therefore, the application of the high-frequency time filtering by Equation (9) is quite justified.
To reveal why the estimates of the turbulence parameters from night-time lidar data are incorrect, we have analyzed the spatially temporal inhomogeneity of the wind. Deviations of the radial velocity from the results of sine-wave fitting were calculated as functions of the azimuth angle and the distance between the lidar and the probed volume; that is, the distributions
were calculated by Equation (9).
Figure 3 shows an example of the distributions
obtained from lidar measurements at four consecutive conical scans (
) starting from 23:00 LT in the presence of LLJ with the maximum velocity at a height of about 350 m (see
Figure 2e). It can be seen that for the layer adjacent to the ground up to a height of 250 m, the distributions
vary randomly from one scan to another. There are no doubts that these variations are caused by the turbulence. However, above this layer, random wind variations are transformed into a regular periodic structure, which remains practically unchanged during the lifetime of the LLJ, and disappears along with it.
Then, we consider the distribution averaged over time of measurement estimates of radial velocity fluctuations
calculated by Equation (11) with
= 30.
Figure 4 shows two examples of these distributions.
Figure 5 depicts the dependence
on the azimuth angle at different heights obtained from lidar measurements on 23 July 2016 (the data are taken from
Figure 4b).
It follows from the data shown in
Figure 4 and
Figure 5 that temporal averaging smooths out random variations of
, but the regular deviations of lidar estimates of the radial velocity from the sine-wave fitting persist. The regular wind inhomogeneity leads to the overestimation of turbulent energy and the outer scale of turbulence, if Equation (1) is used to obtain the array of fluctuating components of lidar estimates of the radial velocity. To take into account the influence of the regular wind inhomogeneity, the wind turbulence parameters were estimated from the array of lidar estimates of radial velocity fluctuations corrected by Equation (10).
By analogy with Ref. [
15], we compared the azimuth structure functions of the radial velocity
obtained from experimental data with the use of Equation (10) with the structure functions calculated theoretically for the von Kármán model [
21] for the two-dimensional spatial spectrum of wind velocity fluctuations:
where
, and
is described by Equation (27) from Ref. [
15].
This comparison shows that we failed to completely remove the influence of the regular horizontal wind inhomogeneity. However, the value of the parameter:
which was introduced in Ref. [
15] and characterizes the degree of deviation of the experimental azimuth structure functions from the von Kármán model (12) becomes much less than the values of
obtained without the high-frequency spatiotemporal filtering by Equations (9)–(11). On average,
= 0.2, which is approximately twice as large as the average value of
for the layer of intense turbulent mixing [
15].
Figure 6 depicts the time dependence of lidar estimates of the wind turbulence parameters at heights in the vicinity of the maximal wind velocity in LLJ. It is seen from the figure that the use of Equations (9)–(11) in the calculation of fluctuations of the radial velocity measured by the lidar has allowed us to obtain estimates of the integral scale of turbulence under stable conditions in the atmosphere, which do not contradict the theory and the available experimental data [
1,
18,
19]. The obtained results demonstrate that in the central part of LLJ, the integral scale of turbulence varies from 40 m to 150 m. On average, it is about 100 m. The kinetic energy of turbulence does not exceed 0.03 (m/s)
2, and the dissipation rate is less than
m
2/s
3 in the vicinity of the maximal jet stream velocity.
The results depicted in
Figure 6 were obtained from measurements at a SNR from 0.1 to 0.3. At this SNR and the number of probing pulses used for accumulation of raw lidar data
= 7500, the instrumental error
of estimation of the radial velocity by Equation (9) does not exceed 0.05 m/s [
23]. At this instrumental error,
= 120, and
= 30, the condition:
providing for the high accuracy of the estimation of the dissipation rate from lidar data [
15,
23] is fulfilled even at the very weak turbulence when
does not exceed 10
−6 m
2/s
3, and the variance
is, on average, an order of magnitude less than the values of the kinetic energy of turbulence shown in
Figure 6. According to our calculations, for the data in
Figure 6, the relative errors of estimation of the turbulence energy and the dissipation rate vary from 5% to 7%, and those for the integral scale of turbulence range from 9% to 12%.
3.2. Experiment at the Lake Baikal Coast
The lidar experiment at the Lake Baikal Coast was carried out at the western coast near Listvyanka in August 2015. The StreamLine lidar was installed at the territory of the Baikal Astrophysical Observatory of the Institute of Solar-Terrestrial Physics SB RAS (latitude 51.405230°N, longitude 104.673121°E) several tens of meters far from the Big Solar Vacuum Telescope (BSVT). The goal of the experiment was to study atmospheric internal waves (AIW) [
16]. As at the Basic Experimental Observatory, the lidar measurements were conducted with the use of conical scanning by the probing beam about the vertical axis. The following measurement parameters were set: elevation angle
= 60°; angular scanning rate
= 10°/s; number of probing pulses for accumulation of raw data
= 3000; step in the range
= 30 m. For the pulse repetition frequency of the StreamLine lidar,
= 15 kHz, the duration of data measurement for one ray was
= 0.2 s, the resolution in the azimuth angle was
= 2°, the number of rays was
= 180, and the duration of one complete scan was
= 36 s.
The data of measurements from 01:00 to 08:00 LT of 24 August 2015 turned out to be the most suitable for the investigation of wind turbulence. In this period, the LLJ was observed, and the signal-to-noise ratio (SNR) of the lidar echo signal was large. After processing the raw data, we obtained two-dimensional height–time distributions of the signal-to-noise ratio, wind speed and direction, and the vertical component of the wind velocity vector. These distributions are shown in
Figure 7. Each vertical profile was retrieved from measurements for one scan for 36 s. It is seen that the LLJ observed for 7 h was localized at heights of 450–750 m. At a height of about 300 m, the wind direction drastically changed to the opposite one: below this level, the wind was directed from the land, while above it, the wind was directed from the lake. At 07:00 LT, an internal wave occured. It was observed for 30 min. This wave is clearly seen in
Figure 7d for the vertical component of the wind velocity. The localization of the wave is shown by a black rectangle in
Figure 7d. The period of this wave was approximately 6 min. The same wave motion is seen in the data for the average wind speed (
Figure 7b) and direction (
Figure 7c).
In contrast to the lidar experiment at BEO (see the previous section), the analysis of the arrays of estimates of radial velocity fluctuations for the measurement period from 01:00 to 08:00 LT on 24 August 2015 did not reveal a regular inhomogeneity in the distributions , even in the cases of extremely weak turbulence and AIW events. The analysis of the azimuth structure functions of radial velocity that was measured, in particular, at the heights of the central part of the jet stream has shown that these structure functions are well described by the von Kármán model, as in the experiment at BEO. This follows from the calculation of the parameter by Equation (13).
The dissipation rate of the turbulent energy
can be determined from lidar data by Equation (4) only on the assumption of locally homogeneous and isotropic turbulence, and with the azimuth structure functions of radial velocity fluctuations measured within the inertial subrange of turbulent inhomogeneities [
15]. Therefore, in determination of the dissipation rate, there are no strong restrictions to elevation angles
. However, the kinetic energy of turbulence
cannot be calculated from lidar measurements with the use of conical scanning at elevation angles different from 35.3° [
15,
20] because of the anisotropy of turbulence, due to which the variances of the wind vector components
,
, and
differ. For this reason, for the calculation of the integral scale of turbulence
from the Baikal lidar data, which were obtained under the scanning angle
= 60°, we used Equation (7) instead of Equation (8). The relative contribution of the vertical component of the wind velocity vector increased with the increase of the elevation angle, and the variance of the vertical velocity
was always less than the variances of the horizontal components of the wind velocity
. Thus, at
= 60°, the values of
and, consequently,
calculated by Equation (7) were underestimated in comparison with the estimates of these parameters obtained for the same heights from data measured under scanning angle
= 35.3°.
As in the experiments at BEO, to obtain single estimates of the turbulence parameters (, , and ), we used the raw data measured by the lidar for 30 min. In this experiment, the duration of one scan was 36 s. Thus, 50 scans ( = 50) were made for 30 min. The total number of estimates of the radial velocity in the array used for determination of the turbulence parameters was = 180 × 50 = 9000.
Figure 8 depicts the spatiotemporal distributions of the wind turbulence parameters obtained from the Baikal measurements. These data have the time step
= 36 s and the height step
= 26 m. As an example,
Figure 9 and
Figure 10 show the average wind velocity and the time profiles of the turbulence parameters (borrowed from
Figure 8) at the fixed heights, and their vertical profiles at some instants. One can see that turbulence is weak in the layer of 100–200 m and in the top part of the LLJ at heights of 600–700 m. In these locations, the dissipation rate
varied from 10
−6 to
m
2/s
3, and the variance
ranged within 0.005–0.07 (m/s)
2. In the bottom part of the LLJ at heights of 300–500 m from 02:30 to 07:00 LT, the turbulence was much stronger, the dissipation rate took values in the range of 5 × 10
−4 to 2 × 10
−3 m
2/s
3, and the variance changed from 0.1 (m/s)
2 to 0.25 (m/s)
2. The temporal variations of the integral scale of turbulence
at the fixed heights in the layer of 300–500 m are not so large as the temporal variations of
and
at these heights. On average, the integral scale in the layer of 300–500 m is approximately equal to 50 m. With height, the average integral scale increases from 30 m at a height of 100 m to 70–80 m at a height of 500 m, and varies from 80 m to 150 m in the central part of the LLJ at heights of 500–600 m. At the heights where the wind velocity is maximal, the integral scale of turbulence is about 100 m. The same average estimates of the integral scale (~100 m) for the central part of the jet stream at lower heights follow from the lidar measurements at BES at the scanning with the elevation angle
= 35.3°. That is, in the LLJ locations, the integral scale of turbulence is two to three times less than the vertical size of the jet.
The sharp increase of the integral scale
with heights in the layer of 600–700 m, as well as separate peaks of
at heights of 100–300 m, can be caused by errors in estimation of the turbulence parameters. Despite the rather high experimental signal-to-noise ratio (SNR) at which the instrumental error of estimate of the radial velocity
is very small, the value of
may be significantly underestimated [
15,
23] because of the extremely low turbulence observed in the stable atmosphere. As a consequence, the value of
calculated by Equation (7) becomes overestimated in the same proportion.
The relative errors of the estimates of the turbulence energy dissipation rate
, the variance of turbulent variations of the radial velocity
, and the integral scale of turbulence
were calculated based on the closed numerical simulation of lidar measurement of the turbulence parameters by Equation [
22]. The values of
,
,
, and
, which are necessary for the calculation of the errors, were obtained from the raw lidar data of this experiment. In the experiment, the instrumental error
did not exceed 0.12 m/s, at least up to a height of 600 m. At these values, with allowance made for the large number of the degrees of freedom
= 9000, the instrumental error of estimation of the radial velocity exerts practically no influence on the accuracy of determination of the variance of wind velocity fluctuations
in
Figure 8b. According to the calculations based on the simulation, the error
is 6–7%.
The instrumental error more significantly affects the accuracy of estimation of the dissipation rate , especially in the stable boundary layer when the turbulence is very weak. This occurs because the dissipation rate is calculated from the difference of the azimuth structure functions of the radial velocity within the inertial range of turbulence (4), and the number taken as small as possible. It follows from the simulation that at heights of 100–600 m, the relative error of estimation of the dissipation rate may exceed 20% in some cases. This worsens the accuracy of determination of the integral scale of turbulence. However, in the overwhelming majority of the cases, the error of estimation of the dissipation rate is between 5–8%, and the relative error varies within 10–13%.