Research article 21 Mar 2019
Research article  21 Mar 2019
Better turbulence spectra from velocity–azimuth display scanning wind lidar
 ^{1}NTNU, Department of Energy and Process Engineering, Norwegian University of Science and Technology, 7491 Trondheim, Norway
 ^{2}DTU Wind Energy, Technical University of Denmark, 4000 Roskilde, Denmark
 ^{1}NTNU, Department of Energy and Process Engineering, Norwegian University of Science and Technology, 7491 Trondheim, Norway
 ^{2}DTU Wind Energy, Technical University of Denmark, 4000 Roskilde, Denmark
Correspondence: Felix Kelberlau (felix.kelberlau@ntnu.no) and Jakob Mann (jmsq@dtu.dk)
Hide author detailsCorrespondence: Felix Kelberlau (felix.kelberlau@ntnu.no) and Jakob Mann (jmsq@dtu.dk)
Turbulent velocity spectra derived from velocity–azimuth display (VAD) scanning wind lidars deviate from spectra derived from onepoint measurements due to averaging effects and crosscontamination among the velocity components. This work presents two novel methods for minimizing these effects through advanced raw data processing. The squeezing method is based on the assumption of frozen turbulence and introduces a time delay into the raw data processing in order to reduce crosscontamination. The twobeam method uses only certain laser beams in the reconstruction of wind vector components to overcome averaging along the measurement circle. Models are developed for conventional VAD scanning and for both new data processing methods to predict the spectra and identify systematic differences between the methods. Numerical modeling and comparison with measurement data were both used to assess the performance of the methods. We found that the squeezing method reduces crosscontamination by eliminating the resonance effect caused by the longitudinal separation of measurement points and also considerably reduces the averaging along the measurement circle. The twobeam method eliminates this averaging effect completely. The combined use of the squeezing and twobeam methods substantially improves the ability of VAD scanning wind lidars to measure inwind (u) and vertical (w) fluctuations.
Wind speed measurements are an integral element of wind site assessment. Traditionally such measurements have been based on in situ sampling with anemometers attached to tall meteorological masts that reach up to hub height. Such masts are immobile and expensive to erect. It is therefore favorable to implement remotesensing devices, such as conically scanning profiling lidars, that measure wind velocities at adjustable height levels above the ground remotely.
Pulsed and continuouswave wind lidars are the two types of profiling lidars that are currently commercially available. The velocity–azimuth display (VAD) scanning strategy was introduced by Browning and Wexler (1968) and is usually applied for continuouswave profiling lidars like the ZX 300 (previously ZephIR 300) produced by Zephir Ltd. Advanced processing of VADacquired data is the object of investigation here.
Validation studies that compare measurements from meteorological masts and groundbased profiling lidars report good agreement for firstorder statistics, namely the 10 min mean wind velocities and directions (Kindler et al., 2007; Smith et al., 2006; Medley et al., 2015; Kim et al., 2016). The estimation of secondorder statistics of the turbulence in the wind by means of VAD scanning pulsed Doppler lidar was first demonstrated by Eberhard et al. (1989). But such turbulence estimates from VAD scanning lidars deviate from classical measurements with cup or sonic anemometers (Sathe and Mann, 2013; Peña et al., 2009; Canadillas et al., 2010). Sathe et al. (2011) model the secondorder statistics of pulsed and continuouswave profiling lidars. The resulting velocity variances are influenced by the effects that arise from sensing the threedimensional wind field by averaging over spatially distributed volumes. In order to better understand the actual behavior of the lidar in comparison to reference measurements, turbulence spectra of the three wind components u, v and w can provide muchneeded insight. Sathe and Mann (2012) model and analyze turbulence spectra, but only for pulsed lidars that use Doppler beam swing (DBS) scanning. A simplified model for turbulence spectra from VAD scanning wind lidars is presented in Wagner et al. (2009). However, it does not include the effect of crosscontamination and cannot be used to predict the turbulence spectra of real lidars.
The sixbeam method developed by Sathe et al. (2015) is an alternative to VAD scanning that results in more accurate secondorder statistics of turbulence. But its application requires a vertical laser beam and a halfcone opening angle of 45^{∘}, which makes it unusable with commercially available profiling wind lidars.
Newman et al. (2016) propose another method to compensate for the contamination by means of autocorrelation functions derived from collocated mast measurements. This method is, however, only applicable when a meteorological mast is available. In comparing and evaluating the ability of different lidar scanning strategies to measure turbulence, Newman et al. (2016) conclude that crosscontamination of the different velocity components is one of the primary disadvantages of current profiling lidars.
The research presented here demonstrates two methods aimed at overcoming the effects of crosscontamination and averaging along the measurement circle that are inherent in the standard VAD scanning strategy. Both methods are based on modified lineofsight velocity data processing and can be applied to currently available lidars without changes in their hardware. The lineofsight averaging effect remains unresolved.
The first method incorporates Taylor's frozen turbulence hypothesis and introduces a time lag into the wind vector reconstruction process. Bardal and Sætran (2016) measure twopoint correlations of horizontal wind speeds from two meteorological masts that are separated by 79 m in line with the mean wind direction. They find that the crosscorrelation coefficient is around 0.8 when a temporal lag compensates for the time required for the wind to cover the distance between the two measurement points. Without delaying the signal, the crosscorrelation coefficient reaches only half of that value. Applied to VAD scanning lidars, that justifies the assumption that when the processing of lineofsight measurement data is delayed by the time needed to cross the measurement circle, the lidar measurements will be more realistic. This approach is hereafter called “squeezing” and reduces the crosscontamination effect that currently distorts the shape of turbulence spectra acquired with VAD scanning lidars.
The second method is to use only the radial velocities from lines of sight that point into the mean wind direction (downwind) and against it (upwind) to determine the components of the wind that are oriented in line with the mean wind direction (u) and vertical direction (w). This eliminates the averaging along the measurement circle.
The aim of the research presented here is to demonstrate whether one of the two modified data processing algorithms or their combination leads to improved turbulence measurements from standard VAD wind lidars. For each method, we present a numerical model and experimental results. We discuss the effects of the two methods individually and combined.
This research has several practical applications. The reliable elimination of crosscontamination and averaging along the measurement circle would lead to a reduction of the systematic error of wind lidar measurements that is dependent on the prevailing wind conditions and the measurement height. In particular, estimations of the timescale of turbulence could be made with higher certainty, which would support future boundary layer research by means of profiling wind lidars. In addition, estimating the energy content of the wind components at specific wave numbers with higher certainty could also help to better predict the operational wind loads of wind turbines and other structures.
Section 2 summarizes the VAD scanning process and describes, in detail, the averaging and crosscontamination effects it implies for the measurement of turbulence. In Sect. 3 the suggested modified data processing methods are described before they are modeled alongside the conventional processing in Sect. 4. The measurements are described in Sect. 5 before the results are compared with the model predictions in Sect. 6. Section 7 concludes with the most important findings.
2.1 Coordinate system and preliminaries
Figure 1 shows the measurement circle of diameter D_{C} of a VAD scanning lidar and how it is created by the laser beams that are deflected from the zenith by the halfcone opening angle ϕ and rotate around the zenith with continuously changing azimuth angle θ. The beams are focused at a point at distance d_{f} from the lidar, which is located at the origin of a threedimensional lefthanded coordinate system. Five of the laser beams are depicted, four in the cardinal directions and one with an arbitrary azimuth angle. The mean wind direction Θ determined from 10 min intervals is zero when the wind blows from north to south. The wind vector
is composed of the wind components u, v and w that are aligned with the axes of the coordinate system when Θ=0^{∘}. Reynolds decomposition is used for the description of the wind field so that
where u^{′} represents the wind speed fluctuations in all three directions and U is the mean wind velocity vector.
2.2 Taylor's frozen turbulence hypothesis
The frozen turbulence hypothesis published by Taylor (1938) assumes that turbulence is advected by the mean wind velocity U into the mean wind direction Θ. During the transport process the turbulence remains unchanged, i.e., turbulence measured at one point in space gives information about the turbulence found further downwind some time later. That means for a velocity vector field u when U is aligned with the x axis that
The hypothesis is widely used and it is known from experiments that the assumption of frozen turbulence is valid to a high degree for large eddies. For example, Schlipf et al. (2010) measured the inflow velocities of an operating wind turbine at different distances from the rotor plane in order to test the hypothesis of frozen turbulence. They found it to be valid for largescale wind fluctuations with wave numbers $k>\mathrm{1.25}\times {\mathrm{10}}^{\mathrm{1}}$ m^{−1}. Willis and Deardorff (1976) show that the hypothesis lacks validity when
This implies that the validity of the hypothesis depends on the amount of turbulence and that a high degree of validity is expected when the velocity variance is low compared to the mean wind speed.
2.3 VAD measurement principle
Continuouswave wind lidars continuously emit a focused infrared laser beam into the air and detect the small portion of the radiation that is backscattered by particles along the beam path towards the beam's origin. The velocity of the backscattering particles relative to the beam direction is then determined by analyzing the Doppler shift between the frequencies of outgoing and incoming radiation. It is assumed that the backscatterers are lightweight enough to move with the instantaneous wind speed u. The measured radial lineofsight velocities v_{r} are hence equal to the wind velocity projected onto the beam direction. In order to estimate the threedimensional wind vector u, a minimum of three independent lineofsight measurements from different directions must be combined.
When VAD scanning is used, the beam is deflected by a wedge prism by a constant halfcone opening angle ϕ from the zenith and rotated around the zenith with a steadily changing azimuth angle θ. Many radial velocities v_{r} are acquired during one full rotation of the prism. For example in the case of the ZX 300 (previously ZephIR 300), N=49 Doppler spectra are calculated and used to determine the same number of radial velocities. All of them are used to reconstruct one wind vector by applying a leastsquares fit to
where the best fit parameters A, B and C represent the wind data according to
The sign of the radial velocity is usually unknown. We are thus faced with a directional ambiguity of ±180^{∘}, but this does not affect the turbulence analysis here. The wind data v_{hor}, Θ and v_{ver} can be translated into wind vectors u easily.
The wind velocity estimations that result from this processing underlie several effects that distinguish them from onepoint measurements. These effects can be divided into

averaging
 a.
along the lines of sight
 b.
along the measurement circle and
 a.

crosscontamination
 a.
due to longitudinal separation
 b.
due to lateral separation.
 a.
2.4 Averaging effects
2.4.1 Lineofsight averaging
In situ wind speed measurements taken with cup anemometers or ultrasonic anemometers have a small measurement volume that can be considered a point. Lidar measurements, in contrast, sense wind velocities along an extended stretch of the line of sight of the laser beam. In the case of continuouswave lidars, the laser beam leaves the lidar optics with a diameter that corresponds to its effective aperture size a_{0} and is focused onto a focus point. The distance between the lidar optics and the focus point is the focal distance d_{f}. The signal of the backscattered radiation originates from anywhere along the illuminated beam, according to a distribution function that has its maximum at the focus point and is proportional to the intensity of the laser light along the beam (Sonnenschein and Horrigan, 1971).
A definite range gate, such as for pulsed lidars, is therefore not applicable to continuouswave lidars. Instead, the Rayleigh length l_{R} is a measure of the distance between the focus point and the point at which the cross section of the beam has twice the area of the cross section at the focus point. According to Harris et al. (2006), it is given by
where λ is the laser wavelength and a_{0} is the effective aperture diameter. The Rayleigh length is quadratically proportional to the focal distance d_{f} that increases linearly with the selected measurement height level. The degree of lineofsight averaging is thus strongly dependent on the measurement height level and is higher for larger heights. The values of l_{R}, a_{0} and d_{f} for the lidar used in our experiments are given in Table 1.
The intensity of backscattered radiation is a function of the distance s from the focus point along the beam. It is sufficiently well approximated by a Lorentzian function,
where s is the distance from the focus position (Mikkelsen, 2009).
All Doppler spectra that are retrieved during the radial velocity acquisition time are averaged, and the focus point sweeps over a considerable arc of the measurement circle during this time. This arc length l_{A} is
where N is the number of lineofsight measurements v_{r} taken during one rotation. In experimental data, the arc averaging effect is contained in the radial velocities. In the models here, we account for this by averaging along the measurement circle.
The Doppler spectra of each lineofsight measurement resemble the probability density function of the radial wind velocities along the line of sight (Branlard et al., 2013). But by determining one single velocity value for each lineofsight measurement, the turbulence information they contain is filtered out.
The additional temporal averaging along the lines of sight is very low, as one measurement takes only $\frac{\mathrm{1}}{N}$ s. The effect of lineofsight averaging is very strong for high wave numbers but has some effect on long turbulent structures as well. The effect of lineofsight averaging is considered in the numerical models and the discussion in this study. But none of the presented data processing methods can avoid the lineofsight averaging effect.
2.4.2 Measurement circle averaging
As described in Sect. 2.3 lidars use all measurement data of at least one full rotation of the prism to reconstruct one wind vector. The resulting system of equations is overdetermined, and in order to find a solution a quadratic best fit is applied. The more lineofsight velocities that are used to reconstruct a wind vector, the stronger the averaging and thereby the larger the loss of turbulent kinetic energy in the measurement data. The residual of the best fit is a measure of the degree of this form of averaging but is usually not used in the processing.
The diameter D_{C} of the measurement circle is
with h being the measurement height and ϕ the halfcone opening angle. The spatial separation between the points that one reconstructed wind vector is composed of thus linearly increases with measurement height. The larger the cone diameter, the stronger the circle averaging. Turbulence with a length scale below the diameter of the averaging circle is affected the most.
In addition to the spatial separation of the measurement points along the measurement circle, the acquisition time must be considered. The mean wind motion carries the air while it is probed, which might further increase the separation of measurement points in the mean wind direction. The ZephIR 300 measures one full rotation in 1 s, and the distance the air moves within this time is usually small compared to D_{C}. The effect of temporal averaging is therefore often small compared to the spatial averaging. One example for the path of measurements that is averaged over is given in Fig. 4a. The circle diameter represents the spatial averaging, and the shift alongwind with the speed U represents the temporal averaging.
2.5 Crosscontamination
2.5.1 Crosscontamination due to longitudinal separation
Another cause for differences in the shape of turbulence spectra from onepoint measurements and their counterparts from VAD scanning lidars is crosscontamination of different velocity components. VAD scanning lidars combine measurements from spatially separated locations where differing velocities may prevail as if they were collected at one point. This leads to a redistribution of turbulent energy among the velocity components u, v and w. Lidarderived spectra of one of the components can at certain wave numbers show lower energy values than the original wind spectrum of that component but may also show too high values due to a contribution from a different velocity component. To better understand crosscontamination we divide the effect into two different types of separations. First we look into longitudinal separations, i.e., separation along the mean wind direction. Fluctuations at two points separated in line with the wind are highly correlated. If the assumption of frozen turbulence is correct, the coherence would be 1 for all separation lengths and all wave numbers. One example of crosscontamination of correlated fluctuations between two longitudinally separated points is visualized in Fig. 2. The chosen wavelength of the wind fluctuations equals twice the separation distance. This can be called the first resonance wavelength. The resonance wavelengths are given by
The corresponding resonance wave numbers are
where n=1, 2, 3… The resulting values for the first two resonance points are given in Table 1.
The two beam directions in line with and against the mean wind direction can be used to determine u_{lidar} and w_{lidar} by using the formulas on the righthand side of the figure. This example looks at these two lines of sight. The v component can be ignored because transverse fluctuations are not detected by the upstream and downstream beams. The example demonstrates a case with isotropic turbulence, i.e., arbitrary but identical amplitudes for fluctuations in all orientations. Averaging along the lines of sight is ignored here for simplicity. The first column of graphs in the figure isolates the u fluctuations u^{′} and shows the resulting lidarmeasured signal for the two radial velocities in the upwind and downwind directions, i.e., v_{r1} and v_{r2}. When these two signals are combined in the usual way, the reconstructed wind speed components ${u}_{\mathrm{lidar}}^{\prime}$ and ${w}_{\mathrm{lidar}}^{\prime}$ differ strongly from the real inflow conditions u^{′} and w^{′}. The lidar is blind to wind speed fluctuations in the u direction and instead attributes the fluctuations to some extent to the estimation of ${w}_{\mathrm{lidar}}^{\prime}$. The same is done for w^{′} in the second column, and the resulting effect is the reverse. The vertical fluctuations w^{′} are interpreted solely as amplified fluctuations of ${u}_{\mathrm{lidar}}^{\prime}$.
The last column combines the two previous cases and shows the resulting distribution of amplitudes that depends on the halfcone opening angle ϕ. When ϕ<45^{∘} the lidar is more sensitive to vertical variations than to horizontal ones, and the contamination of u^{′} caused by w^{′} is more severe than vice versa.
In a more realistic situation, turbulence is nonisotropic and the amplitude of w^{′} at this first resonance wave number is often considerably lower than the amplitude of u^{′}, which leads to a different distribution of contamination, which can be estimated as follows. We use Eqs. (31) and (33) to define the lidarderived variance in the u direction:
In general, the differences of the lineofsight velocities aligned with the mean wind Δv contain contributions from wind fluctuations in the u and w directions Δv_{u} and Δv_{w}, respectively. Here we look at the resonance case in which Δv_{u}=0 and thus Δv=Δv_{w}. We get
when ϕ=30.6^{∘} as for the lidar we used in this study. The subscript “res” indicates that the equation is only valid for inflow fluctuations at resonance, as in the example given before.
In Sect. 4 we develop a model to predict lidarderived spectra. This model was used to create the plots shown in Fig. 3. Figure 3a shows the modeled spectra of the wind components, u_{wind} and w_{wind}, as solid black and red lines. The parameters of the underlying spectral tensor are given in Table 1. They were chosen to best represent the wind conditions found during the experiment presented in Sect. 6. The model was used to estimate the u component of the wind from two lidar beams that point in the upwind and downwind directions. Also here, we did not include lineofsight averaging to isolate the effect of crosscontamination. The principle of the setup is the same as explained for Fig. 2 but now we see results for all inflow wave numbers and use anisotropic turbulence. The resulting lidarderived spectrum u_{lidar,sum} of the u component of the wind is the sum of the lidar's interpretation of the wind components u_{lidar,u} and u_{lidar,w}. We see that the lidarestimated spectrum of u_{lidar,sum} lies a bit below the target spectrum of u_{wind} for most wave numbers but not at the first and second resonance points that are marked with grey dashed vertical lines. There it exceeds the target spectrum. The reason becomes apparent when we look at the components u_{lidar,u} and u_{lidar,w} that u_{lidar,sum} is composed of. We find that the lidar sees u_{wind} nearly to its full extent for very low wave numbers but when we come close to the resonance points u_{lidar,u} drops to zero. The contribution of the vertical wind u_{lidar,w} shows a mirrored behavior and is amplified according to Eq. (14) since ϕ<45^{∘}.
2.5.2 Crosscontamination due to lateral separation
When the lines of sight under consideration are not longitudinally but laterally separated, they do not face resonance but instead a second form of crosscontamination. The strength of the contamination depends then on the coherence of the turbulence for the given lateral separation. When the fluctuations at the two selected focus points are very coherent i.e., their correlation is close to unity, we can expect that the lidarderived wind speed estimates are correct and no crosscontamination occurs. This can be observed at very low wave numbers at which a high degree of coherence is expected. The other extreme is found at the other end of the spectrum at which small fluctuations measured at both focus points are uncorrelated. The lidarderived spectrum is there a linear combination of the variances of the involved components v and w according to
In the case of fully uncorrelated fluctuations we know that $\mathrm{\Delta}{v}_{v}={v}^{\prime}\mathrm{sin}\mathit{\varphi}$ and $\mathrm{\Delta}{v}_{w}={w}^{\prime}\mathrm{cos}\mathit{\varphi}$ and the variance ${\mathit{\sigma}}_{v,\mathrm{lidar},\mathrm{unc}}^{\mathrm{2}}$ of the lidarderived v velocity is
for the lidar with a halfcone opening angle of ϕ=30.6^{∘}. These two situations and all cases in between are shown in Fig. 3b. The difference from the plots in Fig. 3a is that the two beams that point into and against the v direction are used here to estimate the vspectrum v_{lidar,sum}. The target spectrum of the v component of the wind v_{wind} is given as well as the wspectrum w_{wind} that contaminates the signal. From the v_{lidar,v} and v_{lidar,w} curves it can be seen that at very low wave numbers hardly any contamination occurs but mainly because the wcomponent w_{wind} itself contains a low energy density at low wave numbers. As it increases for higher wave numbers, the contamination also becomes more severe. In this example w_{wind} dominates the lidar spectrum v_{lidar,sum} for all wave numbers above approximately ${k}_{\mathrm{1}}=\mathrm{1.4}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1}. The result is that the lidar overestimates the v variances for all wave numbers. Such an effect is also reported by Wyngaard (1968). Thus, it is essential for accurate turbulence measurements to minimize spatial separation.
VAD scanning along the whole measurement circle is more complex than using only two beams. Examining the two beams aligned with or perpendicular to the mean wind direction is not sufficient to fully understand the effect of crosscontamination. For circle scans, all three wind speed components are involved in contaminating all the beams that do not point in the four cardinal directions. We refer to the model presented in Sect. 4.1 and especially Eqs. (24), (25) and (26) of the spectral weighting functions therein to better understand which components influence another.
The lidar can also be configured to perform a socalled 3 s scan, in which one measurement cycle is built from data from three full rotations. This limits the crosscontamination but comes at the cost of much stronger averaging along the measurement circle, especially in strong wind cases, and a sampling rate that is 3 times slower. The ability to measure turbulence with this approach is so weak that it is not further investigated in this paper. Instead, the next chapter suggests two methods that can be used to reduce both averaging and crosscontamination.
3.1 Squeezed measurement circles
In conventional VAD data processing, each measurement cycle consists of the radial velocities that are acquired during one full rotation of the prism. The data used in the reconstruction of one wind vector thus originates from an air volume with the shape of a cone with a diameter of D_{C} at the height of focus. This results in the abovementioned crosscontamination effects.
One way to eliminate the crosscontamination due to longitudinal separation and mitigate the averaging along the measurement circle lies in making use of Taylor's frozen turbulence hypothesis. As mentioned in Sect. 2.2, the hypothesis assumes that turbulent structures are transported by the mean wind motion without changing. This implies that all turbulent structures that enter the measurement cone at one time are identical after some time t when they leave the cone. The time it takes to cross the measurement circle can be estimated for all azimuth directions θ by
where U is the mean wind velocity calculated from conventional VAD processing.
The basic idea here is to introduce a time lag τ=t into the data processing so that each air package that is involved in the reconstruction of one wind vector is scanned twice: once when it enters and again when it leaves the measurement cone. The composition of the measurement circles is shown in Fig. 4 from a coordinate system that is moving with the mean wind U. In this example D_{C}=92.3 m and U=19.5 m s^{−1}. With conventional VAD data processing, the measurement circle is made up of all N consecutive measurements from one cycle (red segment). By contrast, the lower part of Fig. 4 illustrates the introduction of the time delay τ, in which lineofsight measurements from a total of M=6 different measurement cycles are combined to estimate one wind vector (green segments). In other words, with conventional data processing, a measurement cycle is composed of volumes that are widely spatially distributed. The new proposed method picks measurement data taken from what we term a squeezed measurement circle (SMC).
A restriction that comes with the idea of squeezing is that the circle sample rate f_{S} must be high enough to be able to select measurements that were acquired with a time difference reasonably close to τ. That drastically limits the number of measurement heights that should be selected, especially in strong wind cases. For the measurements analyzed in this paper, the lidar scanned continuously at only one height level, which in general makes sense to measure turbulence effectively.
3.2 Twobeam method
The conventional method of averaging data from all available lines of sight to reconstruct threedimensional wind vectors leads to strong averaging along the measurement circle. The method is known to deliver reliable values for the mean wind speed and direction. The directional information allows it to determine the two beams that lie in the upstream and downstream directions. Within the red and green segments of Fig. 4, small red and green rings indicate these particular beams. These two beams can in a second processing step be used to estimate the u and w components of the wind vectors for turbulence estimations. The resulting values are then not averaged along the measurement circle. This is comparable to the DBS method in cases in which the mean wind blows in line with two of the lines of sight. But an advantage of the twobeam method over the DBS strategy is that the relative angle between the mean wind and the two beams is kept constant in any prevailing wind direction. This is an advantage since beams pointing upwind and downwind are immune to contamination by the crosswind component v.
When the twobeam method is combined with the idea of squeezing, then measurements of the u and w components are taken at virtually one focus point following the flow. Only the lineofsight averaging and some minor longitudinal separation among the different locations along the two beams remain.
That is unfortunately not true when estimating the v component of turbulence. Instead, several problems occur. Intuitively, one would choose a beam direction perpendicular to the mean wind direction in order to estimate the v component of the wind. But the radial velocities in this lineofsight direction are often close to zero, and such estimates from continuouswave lidars are usually not reliable (Mann et al., 2010; Dellwik et al., 2010). The transverse v component must therefore be estimated by either VAD processing or selecting a different third beam direction. In the latter case the results would then be influenced by contamination not only from w but also from the u component. This lies outside the scope of this study. Therefore no v data from measurements are processed with the twobeam method.
Like conventional VAD processing, the SMC method and twobeam method require a wind field that is statistically homogeneous in the horizontal directions to yield correct results.
The mathematics of deducing the lidarmeasured spectrum from the secondorder statistics of turbulence is very convoluted. Therefore, we make the assumption that the measurements are performed much faster than it takes the air to move from one side of the scanning circle to the other; i.e., we assume that $\frac{\mathrm{1}}{{f}_{\mathrm{S}}}\ll \mathit{\tau}$. Effectively, the scanning circle is measured continuously. It is difficult to assess the magnitude of the error committed by the assumption of continuous measurements, but we assume it is negligible.
4.1 VAD and SMC
In order to model spectra obtained from conventionally VADprocessed lidar data, we closely follow the method of Sathe et al. (2011). They use the geometry of the lidar scan and its alongbeam weighting function together with information on the spatial structure of surfacelayer turbulence (Mann, 1994). The focus point of the lidar is at a distance d_{f} away in the direction given by the unit vector
where θ is the azimuth angle and ϕ is the halfcone opening angle. The line of sight or radial wind speed that the lidar is measuring is modeled as
where φ is the spatial weighing function of the continuouswave lidar that we assume to be a Lorentzian function with the Rayleigh length l_{R}. u is the threedimensional velocity field suppressing the time argument since we are assuming Taylor's hypothesis. The integration variable s is the distance along the beam from the focus point. The dot product assures that we obtain the lineofsight velocity. We use x, the coordinate aligned with the mean wind vector, instead of time. e_{1} is the unit vector aligned with x.
The w, u and v components of the velocity are calculated by the first three Fourier coefficients of v_{r} as a function of θ; i.e., w is calculated from
In Sathe et al. (2011) variances are calculated for a conically scanning continuouswave lidar and it is trivial to extend that to spectra. Spectra were in fact calculated in Sathe and Mann (2012) but only for a pulsed system. In Sathe et al. (2011) the variances for a conically scanning continuouswave system, e.g., a ZephIR 300 (Smith et al., 2006; Kindler et al., 2007), were given by
where ^{*} means complex conjugation. The spectral weighting functions α, β and γ are
The spectra measured by the conically scanning lidar will be
and likewise for the u and v components.
where $\mathrm{sinc}x=\frac{\mathrm{sin}x}{x}$ is included in Eq. (28) to account for the finite time of circle scanning before a velocity estimate is obtained. L_{f} is the mean wind speed multiplied with this finite time (see Sathe et al., 2011, for details).
To apply the method of squeezing and model the spectra we obtain from SMC processing, we now substitute Eq. (19) with
Following the exact same steps as in Sathe et al. (2011) but using Eq. (29) instead of Eq. (19) we arrive at Eqs. (21)–(23) but with the complex exponential in Eqs. (24)–(26) exchanged with
4.2 Twobeam method
Only the up and downwind beams to determine the u and w components of the wind vector could introduce less averaging than using the whole circle.
When the mean wind is blowing from the north, the unit vectors in the up and downwind directions are called n^{u} and n^{d}, respectively. Their unit vectors are
and with the opposite sign on the first component for n^{d}.
Parallel to Eq. (19) the lineofsight velocity measured by the upwind beam is assumed to be
The u component estimated by the lidar is normally
where
The correlation function of Δv is
Expanding the product inside the ensemble average (〈〉) and using the definition of the correlation tensor of the velocity field, ${R}_{ij}\left(\mathit{r}\right)\equiv \u2329{u}_{i}\left(\mathit{x}\right){u}_{j}(\mathit{x}+\mathit{r})\u232a$, one obtains
Now we use the relation between the velocity covariance tensor and the spectral velocity tensor
where $\int \mathrm{d}\mathit{k}\equiv {\iiint}_{\mathrm{\infty}}^{\mathrm{\infty}}\mathrm{d}{k}_{\mathrm{1}}\mathrm{d}{k}_{\mathrm{2}}\mathrm{d}{k}_{\mathrm{3}}$, to express the autocovariance function as
By interchanging the order of integration of k and the s's we can cast the expression in terms of the Fourier transform of φ, which in the case of a Lorentzian function is $\widehat{\mathit{\phi}}\left(k\right)=\mathrm{exp}({l}_{\mathrm{R}}k\left\right)$. Thereafter, we Fourier transform R_{Δv} with respect to r to obtain the spectrum F_{Δv}(k_{1}). After that process the first term in Eq. (38) becomes
and upon rearrangement we finally obtain
The derivation of the spectrum obtained from squeezed processing is parallel to the normal spectrum. The only difference lies in the definition of Δv. Now we define it as
Using the exact same steps that led to Eq. (39), we see that the cosine term in that equation has to be substituted with 1 and we get
To obtain the spectrum of u, F_{Δv,(s)} simply has to be divided by $({n}_{\mathrm{1}}^{\mathrm{u}}{n}_{\mathrm{1}}^{\mathrm{d}}{)}^{\mathrm{2}}$ according to Eq. (33).
When obtaining the spectrum of v, we simply exchange the unit vectors of the up and downwind beams n^{u} and n^{d} in all equations by the values of the west and eastbound beams n^{w} and n^{e}. In order to obtain the spectrum of w, Δv defined in Eq. (34) has to be replaced by the sum of both radial velocities v^{u}(x)+v^{d}(x), and F_{Δv,(s)} must eventually be divided by $({n}_{\mathrm{3}}^{\mathrm{u}}+{n}_{\mathrm{3}}^{\mathrm{d}}{)}^{\mathrm{2}}$.
To compare the different methods to calculate spectra from a lidar, Eqs. (39) and (41) have to be evaluated with a model for the spectral tensor. We chose the spectral tensor from Mann (1994) and select the model parameters so that the model spectra resemble the spectra from available sonic measurements. The selected parameters are L=65 m, Γ=4 and αϵ^{⅔}=0.023 ${\mathrm{m}}^{\frac{\mathrm{4}}{\mathrm{3}}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{2}}$. The unfiltered u target model spectrum that we will later compare the model results against is given by
and parallelly for the second and third wind components. The model was tested by comparing the theoretical spectra with results from processing computergenerated wind field turbulence data (Mann, 1998) and was found to predict all four data processing methods, i.e., VAD, SMC, twobeam and squeezed twobeam, accurately for all three wind speed components.
5.1 Test site and instrumentation
The test data were collected at the Danish National Test Center for Large Wind Turbines at Høvsøre. The test site is located in West Jutland, Denmark, 1.7 km east of the North Sea. Apart from the dunes along the coastline, the terrain is nearly flat. The Høvsøre meteorological mast is located to the south of a row of five wind turbines. The reference data were acquired with a Metek USA1 sonic anemometer that is mounted at 80.5 m in height above the ground. It is attached to a 4.3 m long boom pointing north. Mast effects can be observed when the wind is blowing from the south. Turbine wake effects influence the measurement signal when the wind blows from the north. For the data set in this study, the inflow is undisturbed. A detailed description of the test site is given in Peña et al. (2016).
Collocated with the meteorological mast, the lidar measurements were taken by a Qinetiq lidar that was configured to continuously scan at 78 m above the ground. The lidar is comparable to the current ZX 300 (previously ZephIR 300) but the effective aperture size is slightly lower, which results in a longer Rayleigh length and thus greater lineofsight averaging. The lidar was equipped with an optoacoustic modulator that makes it possible to detect the direction of the radial velocities. Lineofsight velocities calculated from the centroid of the Doppler spectra are used in the data processing. The precision of these lidar measurements is not exactly known but is in general better than 1 % (Pedersen et al., 2012).
Measurement data of 32 subsequent 10 min intervals are used. The data were acquired on 20 November 2008 between 10:30 and 15:50 local time. The mean wind velocity measured by the sonic anemometer during this period varied from 14.2 to 22.6 m s^{−1} with an average of 19.5 m s^{−1} and a standard deviation of 2.0 m s^{−1}. The turbulence intensity varied from 4.7 % to 14.0 %, with a mean of 8.8 % and standard deviation of 2.0 %. The wind blew from the northwest and the atmospheric stability was neutral. Table 1 summarizes the most important information about the experimental setup.
5.2 Data processing
The time series of all 10 min intervals derived from all processing methods are used to compute turbulence spectra. The measurement rate for the lidar is 1 Hz. Although it would have been possible in the twobeam processing to calculate measurement values with a rate of 2 Hz by using every newly retrieved radial velocity together with its predecessor, it was decided to use only independent measurements acquired every full second. The sonic anemometer measures with a rate of 20 Hz. These highfrequency data are downsampled by the use of the MATLAB function “resample” to a frequency of 1 Hz. The function includes a lowpass filter to avoid antialiasing. The data rate is thus for all methods 1 Hz. The analyzed frequency range from $\frac{\mathrm{1}}{\mathrm{600}}$ to $\frac{\mathrm{1}}{\mathrm{2}}$ Hz equals the wave number range from roughly $\mathrm{5.4}\times {\mathrm{10}}^{\mathrm{4}}$ to $\mathrm{1.6}\times {\mathrm{10}}^{\mathrm{1}}$ m^{−1}. The spectra are then averaged for all intervals and the results are then binned into 30 logarithmically spaced wave number intervals spread across the wave number axis to avoid high density of values and maintain readability towards higher wave numbers.
The effects of detrending (Hansen and Larsen, 2005) and spike removal (Hojstrup, 1993) on the spectra were both negligible for this data set, so neither was applied here.
6.1 u spectra
Figure 5 shows the spectra of the u fluctuations for all processing methods from measurement data (triangle markers) and the corresponding model predictions (solid lines). We will first discuss the results from processing the whole measurement circle shown in Fig. 5a, followed by the discussion of the results of the twobeam method, shown in Fig. 5b.
6.1.1 Circle processing
To begin with, the model predictions of conventional VAD processing and the new SMC method are compared against each other and with regard to the true u target model spectrum acquired from the spectral tensor according to Eq. (42). The model prediction of the conventionally processed VAD lidar data shows some attenuation of the spectral energy even for very low wave numbers. This can be partly explained by the infinitely long tails of the lineofsight averaging function given in Eq. (8). That means that even very large eddies are slightly weakened by the underlying Lorentzian function. Averaging along the measurement circle might also have some small additional impact on largescale turbulence. Both averaging effects become more and more severe for increasing wave numbers until the measured spectral energy reaches values close to zero at roughly ${k}_{\mathrm{1}}={\mathrm{10}}^{\mathrm{1}}$ m^{−1} and above. The tendency of increasing attenuation with regard to the target spectrum is interrupted around the first resonance frequency that is indicated by a vertical grey dashed line at ${k}_{\mathrm{1}}=\mathrm{3.4}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1}. Here the energy density increases and reaches coincidentally roughly the value of the target spectrum. This behavior is as expected an effect of the crosscontamination with energy from both the w spectrum and to a small extent also from v. A resonance effect at the second resonance frequency is hardly pronounced since the energy is nearly fully consumed by the lineofsight averaging.
The SMC model spectrum predicts a similar shape but without the crosscontamination effect from longitudinal separation. Thus, we find no resonance in the computations. The total variance of the u fluctuations σ^{2}(u^{′}) is lower here since less additional energy from the w component is contained in the u_{SMC} signal. The signal is still contaminated by contributions from other components because the lateral separation cannot be reduced by squeezing. But the averaging along the measurement circle is so strong that for example for wave numbers above around ${k}_{\mathrm{1}}={\mathrm{10}}^{\mathrm{2}}$ m^{−1} less than half of the energy of the target spectrum is expected to be detected by the lidar.
First, when the model is compared with the measurement data, the chosen spectral tensor does not fit the actual wind conditions in the wave number range below ${k}_{\mathrm{1}}={\mathrm{10}}^{\mathrm{2}}$ m^{−1}. The extra energy at low wave numbers compared to the spectral tensor model for this site has been observed before and is related to the inhomogeneous landscape at Høvsøre with its seatoland transition in the main wind direction (Sathe et al., 2015) and mesoscale effects that overlay the expected spectral gap (Larsén et al., 2016). Luckily, this does not severely impede the analysis since the most interesting effects are expected at higher wave numbers and tendencies can still be determined from the relative distances between the markers and lines without matching the absolute values. Next, the comparison of data from sonic measurements and VAD as well as SMCprocessed lidar data shows in the very low wave number range at ${k}_{\mathrm{1}}<\mathrm{3}\times {\mathrm{10}}^{\mathrm{3}}$ m^{−1} that VAD processing and SMC processing produce similar results with a slight tendency towards lower energy densities in the SMCmeasured spectrum that is not found in the model computations. A possible explanation is that the fluctuations of the u and especially the w component in the real wind field are not perfectly correlated, i.e., the frozen turbulence hypothesis that the model assumes is slightly violated. The result is a small contribution of w_{wind} to u_{lidar} that appears to a greater extent in the VADprocessed spectrum. The reason for the difference is that the correlation is closer to unity in the case of SMC processing.
Apart from some exceptions (e.g., at ${k}_{\mathrm{1}}=\mathrm{3}\times {\mathrm{10}}^{\mathrm{3}}$ m^{−1}), a relatively increasing averaging effect towards higher wave numbers is found for the lowest wave numbers as expected. In the wave number range ${k}_{\mathrm{1}}={\mathrm{10}}^{\mathrm{2}}$ to $\mathrm{6}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} the sonic spectrum and the VAD spectrum follow the corresponding modeled spectra nicely through the first resonance point. That shows that the crosscontamination caused by longitudinal separation is present in the measurements and is properly modeled.
The spectrum derived from SMCprocessed data shows a clear tendency towards its modeled spectrum but does not completely reach it. It does not show the resonance effect seen for VAD processing, but the overall energy level is higher than predicted for ${k}_{\mathrm{1}}>{\mathrm{10}}^{\mathrm{2}}$ m^{−1}. It is not possible to determine what causes this deviation. One possible reason is that the model assumes a perfect delay of the measurement timing. In reality this is not possible due to only discrete acquisition times being available. Also the air packages are in reality not always advected with the exact mean wind speed and direction. Both imperfections justify that the behavior of real SMC processing lies in between the modeled SMC and VAD processing.
For ${k}_{\mathrm{1}}>\mathrm{7}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} VAD and SMCprocessed data are nearly identical. As shown in Schlipf et al. (2010), the assumption of frozen turbulence is not valid for high wave numbers. In this region, fluctuations separated by the distances between the relevant focus points are uncorrelated and the squeezing has no effect. The lack of coherence also explains that the values are higher than predicted because the u spectrum is highly contaminated by w and v fluctuations.
6.1.2 Twobeam processing
The plotted model spectrum for the conventional twobeam processing method shows a significantly lower averaging effect compared to whole circle processing methods at all wave numbers except in the very low wave number region, where the methods are expected to perform similarly well.
With the twobeam method it is expected that fluctuations with the highest wave numbers analyzed are to some extent included in the spectrum, while they were close to zero when circle processing was applied. The normal twobeam processing in the model is prone to crosscontamination at both resonance points (vertical dashed lines). This situation is explained in detail in Sect. 2.5. In contrast, the method of squeezing applied to the twobeam processing shows as expected no crosscontamination in the model calculations.
Overall, spectra calculated from the twobeam processed measurement data show good agreement to the model. It is important to keep in mind that, due to the poor fit of the measured spectra of the horizontal wind components and the modeled spectra at low wave numbers, we can compare the relations between the different methods but not absolute values. At low wave numbers, the measured spectra are on average closer to the target spectrum than in the case of circle processing. The slightly lower energy content of squeezed measurements that we observed and explained for circle processing is found here as well. Also, when it comes to deviations from the modeled behavior, like for example the higher energy density at some wave numbers (e.g., ${k}_{\mathrm{1}}=\mathrm{3}\times {\mathrm{10}}^{\mathrm{3}}$ m^{−1}), we find similar tendencies as in circle processing, and the reason is likewise unclear.
The strong crosscontamination at the first resonance frequency is clearly represented in the normal twobeam processing and can be completely avoided by squeezing the two focus points to virtually one point. It is worth mentioning that the squeezing procedure works more like expected when applied to the twobeam method than when applied to the circle processing. This can be explained by the error caused by not having continuous but only discrete delaying times τ available. The relative impact of this error is lower in the case of the twobeam method because then the maximum separation distance D_{C} must be compensated for. In circle processing mode, the shorter separations for which the relative error is larger also contribute to the result.
At ${k}_{\mathrm{1}}>\mathrm{7}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} the two processing methods result in nearly identical values again, and we assume the lack of coherence of short eddies to also be the cause here.
6.2 v spectra
Figure 6 shows the spectra of the v fluctuations for all available data processing methods from both measurement data (triangle markers) and the corresponding model predictions (solid lines). Also here, we first discuss the results from processing the whole measurement circle shown in Fig. 6a, followed by the discussion of the results of the twobeam method shown in Fig. 6b.
6.2.1 Circle processing
The modeled spectra of conventionally VADprocessed lidar measurements predict energy densities that slightly exceed the target spectrum for very long fluctuations with ${k}_{\mathrm{1}}<\mathrm{1.3}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1}. This behavior can be explained by uncorrelated w fluctuations between the eastern and western sides of the measurement circle that contaminate the v signal. This contamination is slightly stronger than averaging that is very weak at low wave numbers.
By contrast, fluctuations shorter than approximately ${k}_{\mathrm{1}}=\mathrm{1.3}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} appear dampened in the spectrum, and fluctuations with higher wave numbers ${k}_{\mathrm{1}}>{\mathrm{10}}^{\mathrm{1}}$ m^{−1} are not even present in the v spectrum due to the strong averaging. Unlike the u spectrum, the v spectrum does not have characteristic behavior around the first resonance wave number. This is not surprising because the lines of sight that are the most important for the detection of v fluctuations lie, according to Eq. (26), orthogonal to the mean wind direction in which turbulence is advected. Thus, no resonance occurs.
When the model spectrum for SMC processing is analyzed, we find a higher variance for all wave numbers above approximately ${k}_{\mathrm{1}}=\mathrm{1.3}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1}. Reduced averaging along the measurement circle is the reason for the higher energy in the SMC spectrum. It is caused by the following: the process of squeezing reduces the longitudinal separation of the focus points ideally to zero while the lateral separation remains unchanged. We know that the lines of site perpendicular to the mean wind direction on both sides of the measurement circle are the most important for the determination of v_{lidar}. Let us assume these are the easterly and westerly beams. The exact east and westbound beams are not affected by the process of squeezing. But for example the northeast and the southeast beams (respectively the northwest and northeast on the other side) see different turbulent structures in conventional VAD processing. With SMC processing, these two beams see the same structure. In the subsequent calculation of the v component all lines of sight are combined and the pairs of radial velocities that lie in line with the mean wind contribute with the average of their amplitudes. This average of amplitudes is lower than the common amplitude measured by the beam pairs under SMC processing. More simply, there is less averaging along the measurement circle when SMC is applied. As a result, the spectrum of SMC shows higher energy densities for all wave numbers at which uncorrelated fluctuations dominate.
Now we compare the measurements with the model. Unfortunately, similar to the u fluctuations, the target spectrum does not represent the sonic measured values properly, especially for low wave numbers. We will therefore concentrate on the tendencies and proportions between the spectra from different methods. While the model predicts the behavior at the lowest wave numbers more or less satisfactorily, we are faced with two outliers at ${k}_{\mathrm{1}}=\mathrm{1.65}\times {\mathrm{10}}^{\mathrm{3}}$ m^{−1} and ${k}_{\mathrm{1}}=\mathrm{2}\times {\mathrm{10}}^{\mathrm{3}}$ m^{−1} for which both the VAD and SMC processing lead to excessive energy estimations. The reason is unclear and not further investigated. At all other wave numbers, the agreement of model and measurements is very satisfactory. In particular, the differences between the two methods are found in the measurements, as predicted. The good agreement between model spectra and measurement spectra at wave numbers above approximately ${k}_{\mathrm{1}}=\mathrm{2}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} might be surprising with regard to the poor agreement of sonic measurements and target spectrum. The reason is that the shape of the lidar v spectra is mainly determined by the crosscontamination from the w component, which, as we describe in Sect. 6.3, agrees better with its model representation.
The identity of VAD and SMCderived measurement spectra that we saw for u fluctuations for ${k}_{\mathrm{1}}>\mathrm{7}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} is found here at ${k}_{\mathrm{1}}>{\mathrm{10}}^{\mathrm{1}}$ m^{−1}. The reason is obvious when we look at the relevant longitudinal separation distances. They are much shorter when processing v fluctuations than u fluctuations, and the assumption of frozen turbulence is more valid for short separation distances. Therefore squeezing can maintain its effect into a somewhat higher wave number region.
6.2.2 Twobeam processing
When the twobeam method is applied, i.e., using only the east and west beams to derive the v component of the wind vector, the method of squeezing has no effect. In comparison with the whole circle processing, the twobeam method is characterized by lower energy estimates at low wave numbers and higher energy estimates at higher wave numbers (see Fig. 6). One reason for the first is assumed to be the lower coherence of v fluctuations separated by the full distance D_{C}. That implies that twobeam processing gets a somewhat lower contribution of v_{wind} to v_{lidar}. A second reason is that there is not crosscontamination from u on v occurring for the twobeam processing. The higher energy content at high wave numbers results from the absence of averaging along the measurement circle.
The model cannot be compared with measurements because the lineofsight velocities of the east and west beams were erroneous. The absolute values we measured are unrealistically biased towards nonzero values. This effect has been previously reported (Mann et al., 2010; Dellwik et al., 2010). We included the model behavior of twobeam processing for the sake of completeness and to show that the availability of reliable measurement data for the east and west beams would be of hardly any use.
6.3 w spectra
Figure 7 shows the spectra of the w fluctuations for all processing methods from both measurement data and the corresponding model predictions. Again, we discuss the results from processing the whole measurement circle first and then the results of the twobeam method.
6.3.1 Circle processing
To begin with, we compare the model predictions of conventional VAD processing and the new SMC method against one another and with regards to the w target spectrum. The results of the actual measurements follow. The model prediction of the conventionally processed VAD lidar data shows some attenuation of the spectral energy even for very low wave numbers. The reason is mainly the infinitely long tails of the lineofsight averaging function and to a lesser extent the averaging along the measurement circle. Both averaging effects become quickly stronger for increasing wave numbers. The spectrum from VAD processed data is expected to drop at the first resonance point marked with a grey dashed vertical line in Fig. 7. This drop is minor due to the overall low energy level present in the spectrum. The spectrum reaches a value near its final minimum with variance values close to zero already at around ${k}_{\mathrm{1}}=\mathrm{5}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1} just after crossing the first resonance point. w fluctuations with higher wave numbers are not detectable with conventional VAD processing. According to the model, the SMC processing improves the situation slightly by removing the longitudinal separation that makes lidar blind to w fluctuations at the resonance points with VAD processing. Squeezing the measurements also helps improve the measurements well above and below the resonance wave number. But still, due to the remaining averaging effects, only a minor fraction of the energy in the vertical wind can be detected with both methods at wave numbers above roughly ${k}_{\mathrm{1}}={\mathrm{10}}^{\mathrm{2}}$ m^{−1}.
The fit between target spectrum and measurement data in the low wave number region is good for the w component. This was not the case for the u and v components. The results of Larsén et al. (2016) show that the spectra for vertical fluctuations are not prone to contributions from the mesoscale spectrum. The measurement data overall support these model predictions and show that the process of squeezing functions well over the entire frequency range in this study. In detail, we only find some mismatch for very low wave numbers at which ${k}_{\mathrm{1}}<{\mathrm{10}}^{\mathrm{3}}$ m^{−1}. The measured spectra lie above the target spectrum here although we expected some attenuation. The discrepancy is caused by the real uwind spectrum being much higher than the underlying target spectrum; see Fig. 5. We already found that largescale u fluctuations are also not perfectly correlated and thus contaminate the measured lidar spectra, which is not considered in the model. At higher wave numbers we find reasonable forecasting of measured w spectra by the model.
6.3.2 Twobeam processing
The modeled twobeam spectra in Fig. 7b lie considerably closer to the target w spectrum for all wave numbers. That can be explained by the absence of circle averaging. The strong influence of resonance visible at the two first resonance wave numbers underlines the importance of squeezing when striving for more realistic spectra from lidar measurements.
At low wave numbers with ${k}_{\mathrm{1}}<{\mathrm{10}}^{\mathrm{2}}$ m^{−1} the measured spectra contain higher energy densities than modeled spectra. A similar but less pronounced effect was found in circle processing only at the lowest wave numbers. The explanation we gave there must therefore be supplemented by mentioning that the assumed decorrelation is stronger for the maximal separations that are involved in the twobeam method. The further comparison of spectra from experiment and model shows that the process of squeezing also leads to the expected effect in the case of using only two beams to determine the w component of the wind vector. As in the case of u fluctuations, this statement must be limited to wave numbers ${k}_{\mathrm{1}}<\mathrm{7}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−1}.
6.4 Extended discussion
The results discussed here are extracted from a single data set that covers one measurement height and a narrow band of mean wind speeds, turbulence conditions and inflow directions at a single location. The reason for working with such a limited data set lies in the fact that very few data are available where a commercial VAD scanning wind lidar, collocated to a meteorological mast, is scanning continuously at one height level, while saving at least the lineofsight velocities. Currently, the only option to save lineofsight velocities acquired by a ZephIR 300 is to stream the data manually to a connected PC. The situation is further complicated by the fact that in the normal “profiling mode” the lidar focuses to a reference height of 38 m periodically for filtering purposes. Therefore, the only known way to focus at one particular altitude continuously is to switch the unit to “turbine mode”. In this way, we acquired some data for the investigation, but their overall quality was lower than the historic data that we eventually selected as the best available data.
In further studies different setups and turbulence conditions should be investigated. Changing the measurement height has the strongest influence on the lidarderived spectra. For example, increasing the measurement height would, first, make the averaging along the measurement circle more severe due to the increased measurement circle diameter. Second, the resonance wave numbers are then shifted towards lower values, which leads to different crosscontamination due to lateral separation. Third, the crosscontamination due to lateral separation becomes even more severe due to the longer separation distances of opposite lineofsight beams. Fourth, a further increase in the focus distance leads to even stronger lineofsight averaging. Fifth, the time lag that is introduced for squeezing must be longer, and the frozen turbulence hypothesis loses some more of its validity. Changing the halfcone opening angle to a smaller value would on the one hand reduce the first three of the aforementioned effects effectively, but on the other hand it would lead to much stronger crosscontamination due to the increased sensitivity to w fluctuations according to Eqs. (14) and (16). Lidar measurements at lower mean wind speeds give the turbulence more time to evolve while crossing the measurement circle, which might lead to a deviation from the predicted spectra at somewhat lower wave numbers than observed in our results. The numerical models will work for all turbulence intensities, and the shape of the spectra is mainly determined by the degree of anisotropy and the turbulence length scale. Atmospheric stability conditions other than neutral would not change the way the lidar measures. But a modified spectral tensor model like the one presented in Chougule et al. (2017) could be used to better compare model values with experimental results.
This paper presents two advanced data processing methods for improving turbulence spectrum estimations with VAD scanning wind lidars, with an aim to reduce crosscontamination and averaging effects. The models of these approaches, developed in Sect. 4, are supported by the comparison with experimental data. Discrepancies can be explained for the most part by the limitations of the frozen turbulence hypothesis that underlies the model calculations yet has slightly reduced validity in real measurements. The fact that the spectra in the experiment do not agree very well with the spectral tensor model is also a cause of differences.
We found that the method of squeezing eliminates the resonance effect caused by the longitudinal separation of combined measurement points successfully. It also considerably reduces the averaging along the measurement circle.
The method of using only two beams for the estimation of the u and w components of the wind vector eliminates the averaging along the measurement circle completely. When it is combined with the method of squeezing, the measurements deviate from the sonic measurements mainly due to lineofsight averaging. This combination of both methods substantially improves the measurability of the w spectrum, which is hardly measurable with current VAD processing.
Accurate measurements of the v spectrum remain difficult, even with the approaches described here. The twobeam method is not applicable to current continuouswave lidars, which in most cases are homodyne. Whether the use of squeezed measurement circles always leads to systematically better results is unclear because the resulting spectra are dominated by contamination from w fluctuations of the wind.
In conventionally processed lidar data, crosscontamination compensates for averaging effects, meaning that in general total variance might be close to target values but for the wrong reasons. For systematically better turbulence measurements from VAD scanning lidars, the findings presented here should be included in raw data processing. Both approaches presented here can be applied to any existing VAD scanning continuouswave profiling lidar unit.
Inquiries about and requests for access to data and source codes used for the analysis in this study should be directed to the authors.
FK contributed the initial idea for SMC, performed the data processing, analyzed the results and wrote the paper. JM suggested the twobeam method, developed the numerical model in Sect. 4 and supplied the measurement data.
The authors declare that they have no conflict of interest.
This research project was supported by Energy and Sensor Systems (ENERSENSE) at the Norwegian University of Science and Technology.
This paper was edited by Marcos Portabella and reviewed by two anonymous referees.
Bardal, L. M. and Sætran, L. R.: Spatial correlation of atmospheric wind at scales relevant for large scale wind turbines, J. Phys. Conf. Ser., 753, 32–33, https://doi.org/10.1088/17426596/753/3/032033, 2016. a
Branlard, E., Pedersen, A. T., Mann, J., Angelou, N., Fischer, A., Mikkelsen, T., Harris, M., Slinger, C., and Montes, B. F.: Retrieving wind statistics from average spectrum of continuouswave lidar, Atmos. Meas. Tech., 6, 1673–1683, https://doi.org/10.5194/amt616732013, 2013. a
Browning, K. A. and Wexler, R.: The Determination of Kinematic Properties of a Wind Field Using Doppler Radar, J. Appl. Meteorol.Clim., 7, 105–113, https://doi.org/10.1175/15200450(1968)007<0105:TDOKPO>2.0.CO;2, 1968. a
Canadillas, B., Bégué, A., and Neumann, T.: Comparison of turbulence spectra derived from LiDAR and sonic measurements at the offshore platform FINO1, 10th German Wind Energy Conference (DEWEK 2010), 17–18 November 2010, Bremen, Germany, 2010. a
Chougule, A., Mann, J., Kelly, M., and Larsen, G.: Modeling Atmospheric Turbulence via Rapid Distortion Theory: Spectral Tensor of Velocity and Buoyancy, J. Atmos. Sci., 74, 949–974, https://doi.org/10.1175/JASD160215.1, 2017. a
Dellwik, E., Mann, J., and Bingöl, F.: Flow tilt angles near forest edges – Part 2: Lidar anemometry, Biogeosciences, 7, 1759–1768, https://doi.org/10.5194/bg717592010, 2010. a, b
Eberhard, W. L., Cupp, R. E., and Healy, K. R.: Doppler Lidar Measurement of Profiles of Turbulence and Momentum Flux, J. Atmos. Ocean. Technol., 6, 809–819, https://doi.org/10.1175/15200426(1989)006<0809:DLMOPO>2.0.CO;2, 1989. a
Hansen, K. S. and Larsen, G. C.: Characterising Turbulence Intensity for Fatigue Load Analysis of Wind Turbines, Wind Eng., 29, 319–329, https://doi.org/10.1260/030952405774857897, 2005. a
Harris, M., Hand, M., and Wright, A.: Lidar for Turbine Control, Tech. Rep. NREL/TP50039154, National Renewable Energy Laboratory, Golden, Colorado, USA, 2006. a
Hojstrup, J.: A statistical data screening procedure, Meas. Sci. Technol., 4, 153–157, 1993. a
Kim, D., Kim, T., Oh, G., Huh, J., and Ko, K.: A comparison of groundbased LiDAR and met mast wind measurements for wind resource assessment over various terrain conditions, J. Wind Eng. Ind. Aerodyn., 158, 109–121, https://doi.org/10.1016/j.jweia.2016.09.011, 2016. a
Kindler, D., Oldroyd, A., Macaskill, A., and Finch, D.: An eight month test campaign of the Qinetiq ZephIR system: Preliminary results, Meteorol. Z., 16, 479–489, https://doi.org/10.1127/09412948/2007/0226, 2007. a, b
Larsén, X. G., Larsen, S. E., and Petersen, E. L.: FullScale Spectrum of BoundaryLayer Winds, Bound.Lay. Meteorol., 159, 349–371, https://doi.org/10.1007/s105460160129x, 2016. a, b
Mann, J.: The spatial structure of neutral atmospheric surfacelayer turbulence, J. Fluid Mech., 273, 141–168, https://doi.org/10.1017/S0022112094001886, 1994. a, b
Mann, J.: Wind Field Simulation, Probab. Eng. Mech., 13, 269–282, https://doi.org/10.1016/S02668920(97)000362, 1998. a
Mann, J., Peña, A., Bingöl, F., Wagner, R., and Courtney, M. S.: Lidar scanning of momentum flux in and above the surface layer, J. Atmos. Ocean. Tech., 27, 959–976, https://doi.org/10.1175/2010JTECHA1389.1, 2010. a, b
Medley, J., Wylie, S., Slinger, C., Harris, M., des Roziers, E. B., Pitter, M., and Mangat, M.: Classification of ZephIR 300 Lidar at the UK Remote Sensor Test Site, EWEA Resource Assessment 2015, Helsinki, Finland, 2015. a
Mikkelsen, T. K.: On mean wind and turbulence profile measurements from groundbased wind lidars: limitations in time and space resolution with continuous wave and pulsed lidar systems, European Wind Energy Conference and Exhibition 2009, 16–19 March 2009, Marseille, France, 6, 4123–4132, 2009. a
Newman, J. F., Klein, P. M., Wharton, S., Sathe, A., Bonin, T. A., Chilson, P. B., and Muschinski, A.: Evaluation of three lidar scanning strategies for turbulence measurements, Atmos. Meas. Tech., 9, 1993–2013, https://doi.org/10.5194/amt919932016, 2016. a, b
Pedersen, A. T., Montes, B. F., Pedersen, J. E., Harris, M., and Mikkelsen, T.: Demonstration of short range wind lidar in a highperformance wind tunnel, in: Poster presented at EWEA 2012 – European Wind Energy Conference & Exhibition, 16–19 March 2012, Copenhagen, Denmark, 2012. a
Peña, A., Hasager, C. B., Gryning, S.E., Courtney, M., Antoniou, I., and Mikkelsen, T.: Offshore Wind Profiling Using Light Detection and Ranging Measurements, Wind Energy, 12, 105–124, https://doi.org/10.1002/we.283, 2009. a
Peña, A., Floors, R., Sathe, A., Gryning, S.E., Wagner, R., Courtney, M. S., Larsen, X. G., Hahmann, A. N., and Hasager, C. B.: Ten Years of BoundaryLayer and WindPower Meteorology at Høvsøre, Denmark, Bound.Lay. Meteorol., 158, 1–26, https://doi.org/10.1007/s1054601500798, 2016. a
Sathe, A. and Mann, J.: Measurement of turbulence spectra using scanning pulsed wind lidars, J. Geophys. Res.Atmos., 117, D01201, https://doi.org/10.1029/2011JD016786, 2012. a, b
Sathe, A. and Mann, J.: A review of turbulence measurements using groundbased wind lidars, Atmos. Meas. Tech., 6, 3147–3167, https://doi.org/10.5194/amt631472013, 2013. a
Sathe, A., Mann, J., Gottschall, J., and Courtney, M. S.: Can wind lidars measure turbulence?, J. Atmos. Ocean. Tech., 28, 853–868, https://doi.org/10.1175/JTECHD1005004.1, 2011. a, b, c, d, e, f
Sathe, A., Mann, J., Vasiljevic, N., and Lea, G.: A sixbeam method to measure turbulence statistics using groundbased wind lidars, Atmos. Meas. Tech., 8, 729–740, https://doi.org/10.5194/amt87292015, 2015. a, b
Schlipf, D., Trabucchi, D., Bischoff, O., Hofsäß, M., Mann, J., Mikkelsen, T., Rettenmeier, A., Trujillo, J. J., and Kühn, M.: Testing of frozen turbulence hypothesis for wind turbine applications with a scanning LIDAR system, ISARS 2010, 28–30 June 2010, Paris, https://doi.org/10.18419/opus3915, 2010. a, b
Smith, D. A., Harris, M., Coffey, A. S., Mikkelsen, T., Jørgensen, H. E., Mann, J., and Danielian, R.: Wind lidar evaluation at the Danish wind test site Høvsøre, Wind Energy, 9, 87–93, https://doi.org/10.1002/we.193, 2006. a, b
Sonnenschein, C. M. and Horrigan, F. A.: SignaltoNoise Relationships for Coaxial Systems that Heterodyne Backscatter from the Atmosphere, Appl. Optics, 10, 1600, https://doi.org/10.1364/AO.10.001600, 1971. a
Taylor, G. I.: The Spectrum of Turbulence, P. Roy. Soc. A, 164, 476–490, https://doi.org/10.1098/rspa.1938.0032, 1938. a
Wagner, R., Mikkelsen, T., and Courtney, M.: Investigation of turbulence measurements with a continuous wave, conically scanning LiDAR, RisøR1682(EN), DTU, Risø, Denmark, 2009. a
Willis, G. E. and Deardorff, J.: On the use of Taylor's translation hypothesis for diffusion in the mixed layer, Q. J. Roy. Meteorol. Soc., 102, 817–822, https://doi.org/10.1002/qj.49710243411, 1976. a
Wyngaard, J. C.: Measurement of smallscale turbulence structure with hot wires, J. Phys. E, 1, 1105–1108, https://doi.org/10.1088/00223735/1/11/310, 1968. a