the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Dual-lidar profilers for measuring atmospheric turbulence
Abstract. A dual-lidar system, consisting of WindCube v2.1 profilers with one lidar oriented horizontally at 45° relative to the other, was deployed to estimate along- and cross-wind variances, velocity spectra, and turbulence intensity (TI) using the "variance method". This approach computes second-order statistics directly from the line-of-sight (LOS) velocity variances to infer the full three-dimensional velocity variance components. It is benchmarked against the "traditional method" commonly used in the wind energy sector, which reconstructs instantaneous velocity components from single-lidar LOS measurements before deriving second-order statistics. Both methods are evaluated at a single altitude using a 30-day collocated dataset, with reference measurements from a sonic anemometer classified by atmospheric stability conditions. Two key performance metrics are considered: the mean relative bias error (MRBE) and the relative root mean square error (RRMSE), as defined by DNV (Det Norske Veritas). Spectral analysis of the velocity components shows that, while the traditional method more closely matches the reference spectra at low frequencies, it tends to overestimate spectral energy at higher frequencies. In contrast, the variance method typically underestimates spectral energy in the along-wind component and overestimates it in the crosswind component. Furthermore, linear regression analysis reveals that the variance method captures 90–97 % of the reference variances across all stability regimes, while the traditional method tends to overestimate—especially for cross-wind variance under unstable and neutral conditions (up to 132 %). Overall, the variance method yields lower MRBE and RRMSE values for both along- and cross-wind TI. Specifically, for along-wind TI, MRBE decreases from 10.1 % to 7.5 % and RRMSE from 19.5 % to 11.7 %. For cross-wind TI, MRBE is reduced from up to 17.6 % to 12.3 %, and RRMSE from 23.4 % to 18.2 %.
- Preprint
(5658 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on wes-2025-179', Anonymous Referee #1, 09 Dec 2025
The comment was uploaded in the form of a supplement: https://wes.copernicus.org/preprints/wes-2025-179/wes-2025-179-RC1-supplement.pdfCitation: https://doi.org/
10.5194/wes-2025-179-RC1 -
AC1: 'Reply on RC1', Maxime Thiébaut, 29 Jan 2026
The comment was uploaded in the form of a supplement: https://wes.copernicus.org/preprints/wes-2025-179/wes-2025-179-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Maxime Thiébaut, 29 Jan 2026
-
RC2: 'Comment on wes-2025-179', Anonymous Referee #2, 13 Dec 2025
Thiebaut et al. present an interesting study on turbulence measurements with profiling lidar. The idea to place two lidars with a yaw angle offset and combining them to a lidar with more beams is interesting and innovative. The plots are well prepared and the manuscript is well written. At mutiple points I feel that the description of the results and the methods are a bit unprecise. The comparison to other lidar configurations and techniques on the contrary goes a bit too far in my opinion, because it cannot be justified with the results of the experiment. I thus suggest the manuscript for publication only after major revision. General and specific comments are given below.
General comments
- The database could be better described. A brief statement on wind speed span, median and mean is given, but for example no information about the distribution of wind direction. Is that dataset statistically significant? I think it is, but it is not shown.
- I am not very convinced about the comparison with the 6-beam method. You cannot easily compare the datasets and the lidar parameters are quite different. The comparison of errors and uncertainties on that basis is not sound.
- Some details of the results are not explained in enough detail. For example, the effects for cross-wind variance and how they depend on wind direction, or the differences in the spectra. Why is the low frequency not the same for all methods, why does the variance method have more high frequency noise etc.Specific comments
p.1, l.21: I do not think that you can say that so generally. There are a lot of people who do VAD with pulsed lidars as well. It has advantages, especially for turbulence retrievals, too.
p.2, l.48: Eberhard et al. 1989 requires a full VAD at 35.3° and provides TKE and the covariances, not the single component variances. Later studies by Smalikho, Stephan, Wildmann and Päschke showed that this method is very accurate, if the lidar "intra-beam" volume averaging effects are corrected in the retrieval.
p.3, l.63f: There has recently been a release of a commercial 6-beam lidar https://halo-photonics.com/lidar-systems/beam-6x/, https://halo-photonics.com/lidar-systems/beam-6x-windpower/. It also does not take 15s for an instantaneous measurement any more. I saw further down the manuscript that you discuss this instrument, but i think it would be fair to mention here already.
p.4, l.112: Despite the fact of the foundation being quite impressive, I am not sure how relevant it is for this study. Information about the wind conditions at the site (wind rose, etc.) could be quite interesting instead.
p.7, l.159: Also for the collected dataset, some statistics would be helpful here: wind rose, histograms, of wind, turbulence, stability for example.
p.8, l.185ff: Please provide the thresholds for the despiking
p.10, Eq.13&14: $U$ remains the absolute velocity?
p.11, Eq.19: You switch here from vector notation to Einstein notation (I think), without explaining it. That could be confusing for readers and should be explained.
p.13, l.310: I assume that the "virtual kinematic heat flux" was calculated using the sonic vertical velocity and sonic temperature? Thus not directly the virtual temperature. Could be confusing if you use the same symbol as for the average virtual temperature from the WXT530.
p.15, Tab.2 and p.16, Fig.5: Comparing the methods with completely different datasets is not sound. You would have to make sure that you have the same amount of data for all sorts of wind bins, wind sectors, stability classes, which I assume is not the case here!?
p.16, l.355f: can you explain why the spectra differ at low frequencies?
p.22, ll.434ff: I think the ideas and comparison to other lidar configurations are a bit superficial and not exactly based on results from this study. I recommend to skip them and focus more on the direct findings of the new variance method.
p.22, l.457: Intermittent turbulence is especially observed in stable boundary layers with strong shear. A violation of homogeneity assumptions cannot be deirectly associated with neutral and unstable conditions alone. Neutral conditions can be perfectly homogeneous over flat terrain, stable conditions can be non-homogeneous with only slightly complex terrain. I think you should be a bit more precise what you mean here.
p.23, l.477: I think this dependency on wind direction should be shown explicitly.
p.24, l.497: You should at least mention that a two-lidar setup doubles the cost at this point.Citation: https://doi.org/10.5194/wes-2025-179-RC2 -
AC2: 'Reply on RC2', Maxime Thiébaut, 29 Jan 2026
The comment was uploaded in the form of a supplement: https://wes.copernicus.org/preprints/wes-2025-179/wes-2025-179-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Maxime Thiébaut, 29 Jan 2026
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 469 | 133 | 31 | 633 | 31 | 29 |
- HTML: 469
- PDF: 133
- XML: 31
- Total: 633
- BibTeX: 31
- EndNote: 29
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1