Articles | Volume 15, issue 21
https://doi.org/10.5194/amt-15-6507-2022
https://doi.org/10.5194/amt-15-6507-2022
Research article
 | 
11 Nov 2022
Research article |  | 11 Nov 2022

Complementarity of wind measurements from co-located X-band weather radar and Doppler lidar

Jenna Ritvanen, Ewan O'Connor, Dmitri Moisseev, Raisa Lehtinen, Jani Tyynelä, and Ludovic Thobois
Abstract

Accurate wind profile measurements are important for applications ranging from aviation to numerical weather prediction. The spatial pattern of winds can be obtained with ground-based remote sensing instruments, such as weather radars and Doppler lidars. As the return signal in weather radars is mostly due to hydrometeors or insects, and in Doppler lidars due to aerosols, the instruments provide wind measurements in different weather conditions. However, the effect of various weather conditions on the measurement capabilities of these instruments has not been previously extensively quantified. Here we present results from a 7-month measurement campaign that took place in Vantaa, Finland, where a co-located Vaisala WRS400 X-band weather radar and WindCube 400S Doppler lidar were employed continuously to perform wind measurements. Both instruments measured plan position indicator (PPI) scans at 2.0 elevation from the horizontal. Direct comparison of radial Doppler velocities from both instruments showed good agreement with R2=0.96. We then examined the effect of horizontal visibility, cloud base height, and precipitation intensity on the measurement availability of each instrument. The Doppler lidar displayed good availability in clear air situations and the X-band radar in precipitation. Both instruments exhibited high availability in clear air conditions in summer when insects were present. The complementary performance in the measurement availability of the two instruments means that their combination substantially increases the spatial coverage of wind observations across a wide range of weather conditions.

Dates
1 Introduction

Accurate wind profile measurements are required by numerous applications, ranging from assimilation in numerical weather prediction (NWP) models to monitoring weather conditions in weather-sensitive industries such as aviation. The World Meteorological Organization (WMO) Statement of Guidance for Global NWP states that “wind profiles at all levels outside the main populated areas are a top priority among variables that are not adequately measured by current or planned systems” (Andersson2017).

Remote sensing instruments, such as Doppler lidars and Doppler weather radars, are appealing technologies for measuring wind, since they can provide wide spatial coverage of near-surface wind measurements. Doppler lidars measure the radial wind velocity based on the Doppler effect on the signal backscattered from aerosols (Werner2005). Given a suitable concentration of atmospheric aerosols, Doppler lidars can measure the radial wind velocity out to ranges of 10 to 20 km, depending on the instrument. However, strong attenuation of the lidar beam by precipitation, fog, and cloud will often limit the measurement availability to much closer ranges (Guo et al.2015; Liu et al.2019).

Operating at longer wavelengths than Doppler lidar, Doppler weather radars measure the radial velocity from scattering by hydrometeors (Rauber and Nesbitt2018). However, their measurement availability is often limited in clear air conditions. Biological scatterers, such as insects in the boundary layer during the warm season, can provide radial velocity measurements in clear air conditions, as can Bragg scattering, caused by fluctuations in air refractive index, which can be observed at the S and C band (Achtemeier1991; Wilson et al.1994; Contreras and Frasier2008; Franck et al.2021). Operating at higher frequencies, such as X band, offers several advantages over S and C band (McLaughlin et al.2009), including ease of co-location with Doppler lidar in critical environments such as airports.

The measurement limitations suggest that, for all-weather capability, neither instrument is sufficient on its own as there will be significant gaps in data availability for particular situations (Oude Nijhuis et al.2018). However, with different scattering media responsible for Doppler lidar and weather radar observations, the combined data availability from both instruments may close the gaps in the data availability significantly. Therefore, we quantify the performance of both instruments in various weather conditions in order to demonstrate the complementarity of radar and lidar wind observations. The study uses observations collected during a 7-month measurement campaign in Vantaa, Finland, where a co-located Vaisala WindCube 400S Doppler lidar and Vaisala WRS400 X-band weather radar were deployed to perform simultaneous wind measurements. These observations were supplemented by surface station measurements of horizontal visibility, cloud base height, and precipitation intensity, in order to characterise the weather conditions.

The rest of this article is structured as follows. Section 2 describes the measurement campaign and instruments. The analysis methods are described in Sect. 3, and the results comparing the Doppler lidar and X-band weather radar Doppler velocity measurements and their availability are presented in Sect. 4. Finally, Sect. 5 summarises the main findings.

2 Data

2.1 Measurement campaign and site

A Vaisala WindCube 400S Doppler lidar (Dolfi-Bouteyre et al.2008; Thobois et al.2019) and a Vaisala WRS400 X-band weather radar were deployed a few metres apart (see Fig. 1) on the roof of a building at Vaisala Oyj headquarters in Vantaa, Finland (601656.6′′ N 245233.9′′ E, 55 m above mean sea level). Both instruments were configured to measure plan position indicator (PPI) scans with an elevation angle of 2.0 from the horizontal. The campaign period was from 1 May to 30 November 2021; however, the WindCube 400S lidar was out of operation due to a broken component from 22 June to 17 August. The lidar and X-band radar measurements are described in detail in Sect. 2.2 and 2.3, respectively.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f01

Figure 1Vaisala WRS400 X-band radar and WindCube 400S Doppler lidar installed on the roof of a building at Vaisala Oyj headquarters in Vantaa, Finland. Photo credit: Raisa Lehtinen.

Download

The measurement location was chosen for convenience reasons; however, it is not optimal. For both instruments, nearby trees blocked specific azimuths, and additionally, for the lidar, a sector to the east was blocked by the X-band radar. The blocked rays are clearly seen in Fig. 2, which presents the fraction of available 2.0 PPI measurements during the entire measurement campaign for each instrument. For analysis purposes, any ray with a measurement availability of less than 5 % was considered to be blocked.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f02

Figure 2Fraction of available measurements during the entire measurement campaign for (a) lidar and (b) X-band radar. The location of the surface weather station is marked with a black cross. Background map tiles by Stamen Design, under CC BY 3.0. Map data by OpenStreetMap, under Open Database License (ODbL).

Since the instruments were located at the same location and altitude, and the analysis presented here is limited to be within the lidar measurement range (i.e. within 14.3 km), we considered the lidar beam and the centre of the X-band radar beam to follow the same path. Additionally, the lidar and X-band radar measurements were temporally synchronised so that the scan starting times were as close as possible. Therefore, this measurement campaign provided a unique opportunity to compare Doppler velocity measurements from both instruments.

2.2 Doppler wind lidar

The Vaisala WindCube 400S Doppler wind lidar is a coherent pulsed Doppler lidar operating at a wavelength of 1.54 µm (Dolfi-Bouteyre et al.2008). Coherent Doppler lidars typically use high pulse repetition rates and accumulate many pulses in order to achieve the desired sensitivity. The high pulse repetition rate limits the maximum range of the system; however, this is not a major issue, since signals close to the maximum range are normally far below the detection threshold. Radial velocities are obtained via the heterodyne technique (Frehlich and Kavaya1991), and the maximum unambiguous velocity depends on the sampling rate. The lidar beam itself is very narrow, and when scanning, the angular resolution of a ray is a function of the pulse accumulation time.

The instrument-provided data were ingested and processed using the open-source package wradlib (Heistermann et al.2013), and the radial velocity data then filtered using the quality flags provided by the instrument and a carrier-to-noise ratio (CNR) threshold of −30 dB. No clutter filtering was applied to the Doppler lidar data.

Other relevant lidar specifications are summarised in Table 1.

Table 1Vaisala WindCube 400S Doppler wind lidar specifications.

Download Print Version | Download XLSX

2.3 Weather radar

The Vaisala WRS400 X-band weather radar uses solid-state transmitters for horizontal and vertical polarisation channels. To achieve the required sensitivity, a 44 µs pulse with frequency modulation and corresponding pulse compression is used (Mudukutore et al.1998; Bharadwaj and Chandrasekar2012; O'Hora and Bech2007). Because the use of such a long pulse leads to a blind region in the vicinity of the radar, a shorter pulse of 1 µs is used in addition (Bharadwaj and Chandrasekar2012). The transition from the short pulse to the long pulse measurements occurs at 13.5 km.

The Doppler observations are performed using a dual pulse repetition frequency (dual-PRF) pulsing scheme with pulse repetition frequencies (PRFs) of 2100 and 1400 Hz. The corresponding maximum unambiguous velocity for this pulsing scheme is 32.66 m s−1. A Doppler clutter filter (Siggia and Passarelli2004) was applied and a post-processing algorithm was used to correct spurious Doppler velocities caused by the dual-PRF velocity unfolding (Holleman and Beekhuis2003; implementation from Altube2020, https://github.com/meteocat/vcor_dual_prf, last access: 2 November 2022). The radar data were ingested and processed using the open-source package Py-ART (Helmus and Collis2016). Other relevant radar specifications are summarised in Table 2.

Table 2Vaisala WRS400 X-band weather radar specifications.

Download Print Version | Download XLSX

2.4 Weather station measurements

The surface station measurements used in the analysis were measured at the Helsinki Airport weather station (WMO code 02974; 601936.1′′ N 245724.3′′ E; 51 m above mean sea level) operated by the Finnish Meteorological Institute (FMI). The station is located 6.7 km northeast from the instrument site (see Fig. 2).

Horizontal visibility, taken as a 1 min average of the measured values, was provided by a Vaisala FS11P present weather sensor. The cloud base altitude was measured by the Vaisala CL31 ceilometer, which has a temporal resolution of 30 s, but the cloud base reported by the instrument is an average from several profiles, with the amount of averaging increasing with increasing altitude. A Vaisala FS11P present weather sensor measured the precipitation intensity, which is taken as the 10 min average reported every 10 min. The measurements are summarised in Table 3.

Table 3Surface weather station measurements used in the analysis.

Download Print Version | Download XLSX

When comparing to the measurement times of the lidar and radar measurements, the nearest (in time) reported surface station measurement value was always used. That is, for horizontal visibility and cloud base height, the measurements were taken from the same minute, while for the precipitation intensity the measurement was taken from the nearest complete 10 min period.

3 Methods

Although the Doppler lidar and X-band radar are co-located and observe along the same path, the measurements themselves have different radial and azimuthal resolutions. Therefore, to compare the Doppler velocity measurements, the measurements were interpolated onto a shared Cartesian grid with 250 m × 250 m pixel size. The grid was limited to 14.5 km from the site (i.e. the maximum range of the lidar). The interpolation was performed with the nearest-neighbour interpolation implemented in the wradlib Python package (Heistermann et al.2013). After interpolation, a median filter was applied in the radar measurements to remove any incorrectly unfolded dual-PRF velocities that were not corrected by the dual-PRF correction algorithm (see Sect. 2.3).

The pixel-wise agreement of the measurements from the closest-in-time lidar and X-band radar scans was then studied using an ordinary least squares regression model. Given the gridded X-band radar measurement ViX and lidar measurement ViL, the regression model is defined as follows:

(1) V i X = β 0 + β 1 V i L + ϵ i ,

where ϵi represents the measurement error, assumed to be independent and identically normally distributed, and the coefficients β0, β1 are estimated by minimising the residual,

(2) S ( β 0 , β 1 ) = i = 1 N V i X - β 0 - β 1 V i L 2 ,

where N is the number of paired measurements over the whole measurement campaign.

The goodness of fit of the linear regression model can be estimated using the coefficient of determination R2, as follows:

(3) R 2 = 1 - SS res SS tot ,

where the residual sum of squares SSres=S(β0,β1), and the total sum of squares is

(4) SS tot = i = 1 N V i X - 1 N i = 1 N V i X 2 .

The R2 term indicates the proportion of variance in the data predicted by the model, with values from 0 to 1. The root mean squared difference (RMSD) and mean error (ME) were also calculated, with the RMSD being defined as

(5) RMSD = 1 N i = 1 N V i X - V i L 2 ,

and ME as

(6) ME = 1 N i = 1 N V i X - V i L .

The regression model and the statistics were calculated with the statsmodels Python package (Seabold and Perktold2010).

Additionally, for each scan, the fraction of measurement bins with valid measurements compared to the total number of bins in the scan was calculated. The fraction was calculated taking into account the blocked beams separately for each instrument, so that a value of 1 indicates that every non-blocked bin had a valid measurement value. These fractions were then used to compare the measurement performance in different weather conditions.

Finally, the performance of the instruments as a function of the measurement range was studied by calculating the fraction f(r) of valid measurements at range r, as follows:

(7) f ( r ) = 1 N scans N rays i = 1 N scans j = 1 N rays I ( r , j , i ) ,

where Nscans is the number of scans, and Nrays is the number of non-blocked rays for the instrument. I(r,j,i) is an indicator function for whether the measurement Vr,j,i in ray j at range r in scan i is valid, with I(r,j,i) defined as follows:

(8) I ( r , j , i ) = 1 , V r , j , i valid 0 , V r , j , i not valid .
4 Results

4.1 Comparison of measured radial velocities

Figure 3 presents a scatterplot comparing gridded Doppler lidar and X-band radar Doppler velocity measurements. The results of the statistical comparison are given in Table 4. A value of R2=0.96 suggests that the agreement between the two observations is rather good. We should point out, however, that some artefacts can be seen in Fig. 3. For example, there are some points where the lidar is measuring a zero velocity, but the X-band radar is not. This is due to differences in the processing and clutter filtering. Additionally, the scattered points located approximately symmetrically on both sides of the main point cloud could be caused by errors in the radar dual-PRF unfolding algorithm when the measurement has a low signal-to-noise ratio (SNR) that were not removed by our postprocessing.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f03

Figure 3Scatterplot of Doppler lidar versus X-band radar radial Doppler velocity measurements. The orange line indicates the linear fit to the data, and the black dashed line indicates the one-to-one agreement.

Download

Table 4Statistics on the agreement of Doppler lidar and X-band radar radial Doppler velocity measurements.

Download Print Version | Download XLSX

The statistics in Table 4 show that the mean error between the measurements is only −0.042 m s−1. Additionally, there is no significant correlation between the linear regression residual and lidar CNR or X-band SNR (not shown). Therefore, we can see that the measurements by the two instruments agree well, and possible differences in the measurements are better explained by the different measurement volumes and scatterers, and the time difference between the measurements, than instrument error or bias between the instruments.

Note that Fig. 3 and the statistics in Table 4 have been calculated from measurements where both Doppler lidar and X-band radar measurements were available, which in our data correspond to approximately 5 % of all measurements. The differences in measurement availability for the two instruments are discussed in the following sections.

4.2 Measurement performance

The data availability, computed as the fraction of measurements available at each range bin over the entire measurement campaign period, provides a description of the overall performance of each instrument. This is shown in Fig. 4 as the fraction of measurements available with range over the entire measurement campaign period. The availability as the function of range depends on both the specific instrument and local conditions, e.g. aerosol content for a Doppler lidar; therefore, when examining the measurement performance in different conditions, we need to consider the effect relative to Fig. 4.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f04

Figure 4Fraction of available measurements, i.e. data availability, as a function of measurement range for Doppler lidar (dashed line) and X-band radar (solid line) measurements.

Download

The Doppler lidar displays very high availability at short ranges, with availability decreasing almost linearly with range. At 14.3 km in range, the measurement availability is about 12 %.

For the X-band radar, the overall availability is much lower than for the Doppler lidar. Again, availability is better at short ranges and decreases with range. The change in pulse length at 13.5 km is clearly visible, with the 14.6 dB greater sensitivity of the long pulse responsible for the increase (see Table 2). Note that, here, there has been no blending of short and long pulse.

4.3 Measurement performance as function of horizontal visibility

Next we study the effect of horizontal visibility on the measurement performance. Figure 5 presents the distribution of the fraction of valid measurements for both instruments as a function of horizontal visibility, separated into summer (May to September) and autumn (October to November) months. Only measurements with a corresponding valid horizontal visibility measurement were included, which was approximately 90 % of all measurements. The Doppler lidar measurements exhibit similar distributions, even though the majority of the cases in the summer months have visibility above 35 km (Fig. 5a). Figure 5a and c show a linear increase in Doppler lidar measurement availability with horizontal visibility up until a maximum at about 40 km, after which the availability then decreases, presumably due to the decreasing atmospheric aerosol number concentrations responsible for increasing visibility.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f05

Figure 5Fraction of valid measurements as a function of horizontal visibility from May to September and October to November 2021 for (a, c) Doppler lidar and (b, d) X-band radar.

Download

A very different response is seen for the X-band radar. In the summer months (Fig. 5b), insects are responsible for a large portion of the distribution towards the lower right quadrant (fraction values from 0 to 0.9 and horizontal visibility above 30 km), which then mostly disappear during the autumn (Fig. 5d). The portion of high fractions with low visibility in the top left corner are attributed to precipitation and are present during both seasons.

These effects can also be seen in Fig. 6, which shows the measurement availability as a function of range for different horizontal visibility intervals. For low horizontal visibilities, the Doppler lidar measurement availability decreases rapidly with increasing range (Fig. 6a). As visibility increases, the Doppler lidar measurement availability increases up until visibility reaches about 40 to 50 km, after which the availability decreases. The relatively larger decreases seen beyond a distance in range of 4 km for some curves is likely due to the shape of the telescope focus function that the Doppler lidars are configured with (Pentikäinen et al.2020).

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f06

Figure 6Fraction of available measurements as a function of measurement range for (a) Doppler lidar and (b) X-band radar calculated separately for different horizontal visibility classes. Each line corresponds to a horizontal visibility class with a bin width of 5000 m.

Download

For the X-band radar (Fig. 6b), the highest measurement availability at all ranges is achieved at low horizontal visibility due to precipitation. Note that, since the low horizontal visibility cases will also include situations with fog where the X-band radar will not measure well, the average measurement availability of the cases here is at most 0.7. The availability decreases as the horizontal visibility increases, and the lowest availability is obtained at highest horizontal visibility values. The jump in availability at the pulse change is clearly visible, but the horizontal visibility affects the magnitude of the jump, with the largest jump occurring at low visibility.

Two example cases demonstrating the discussed effects are shown in Figs. 7 and 8. Both cases are clear air cases with no precipitation. During the first case in Fig. 7 from 28 August 2021, the X-band radar had a strong insect echo and thus was able to measure throughout almost the entire measurement area (excluding beam blockage). The horizontal visibility during the measurements was approximately 45 km, resulting in the Doppler lidar being able to measure out to close to the maximum range.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f07

Figure 7Data example from 28 August 2021 at 12:28 UTC, showing PPI scans at 2 of X-band (a) radar reflectivity and (b) Doppler velocity and Doppler lidar (c) carrier-to-noise ratio (CNR) and (d) Doppler velocity. The Doppler velocity colour map is from Crameri (2021).

The second case (Fig. 8) on 17 June 2021, on the other hand, occurred on a day with few insects in the boundary layer, so the X-band radar had no measurements beyond the first few kilometres. The horizontal visibility during the case was approximately 63 km, and correspondingly, the Doppler lidar was not able to measure consistently beyond 10 km in range.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f08

Figure 8Data example from 17 June 2021 at 16:28 UTC, showing PPI scans at 2 of X-band (a) radar reflectivity and (b) Doppler velocity and Doppler lidar (c) carrier-to-noise ratio (CNR) and (d) Doppler velocity. The Doppler velocity colour map is from Crameri (2021).

Since the radar measurement performance in clear air is determined by whether there are insects in the air or not, it is important to note that insects are not always an unbiased target for measuring Doppler velocity (Martin and Shapiro2007; Rennie2012; Hannesen et al.2014). To investigate whether insects should be excluded due to bias in our location, we selected the data from Fig. 3, where the X-band radar differential reflectivity ZDR≥5 dB and co-polar correlation coefficient ρHV≤0.9, as measurements with those values are expected to be from insect scatterers (Zrnic and Ryzhkov1999; Stepanian et al.2016; Jatau et al.2021). These measurements correspond to approximately 10 % of the original data. The resulting scatterplot is shown in Fig. 9.

For this subset of the data, the coefficient of determination R2=0.95 is slightly decreased and bias (ME =0.078 m s−1) slightly increased, but the RMSD =1.13 m s−1 is decreased. Visually, the artefacts seen in the original scatterplot (Fig. 3) are not present in the new scatterplot (Fig. 9). This leads us to conclude that, for our location, insects do not cause significant errors in the radar Doppler velocity measurements. This may be by virtue of smaller insects that act as passive tracers dominating in our location compared to larger insects. However, this is not the case everywhere, so the issue should be investigated separately for each location.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f09

Figure 9Scatterplot of Doppler lidar versus X-band radar radial Doppler velocity measurements, where X-band radar ZDR≥5 dB and ρHV≤0.9. The orange line indicates the linear fit to the data, and the black dashed line indicates the one-to-one agreement.

Download

4.4 Measurement performance as function of cloud base height

Similar analysis as for the horizontal visibility was performed with respect to cloud base height. The cloud base height analysis was limited to the lowest 3 km, as this was considered most relevant for scanning at low elevation angles. Only measurements with a corresponding valid cloud base height measurement were included, which was approximately 72 % of all measurements.

Figure 10 shows how the fraction of valid measurements varies with respect to cloud base height. The distributions show the impact of the severe in-cloud attenuation experienced by Doppler lidar. Assuming a reasonably homogeneous cloud layer, the valid Doppler lidar measurements are limited to the range at which the lidar beam impinges the cloud layer. A cloud layer at 250 m would be intercepted at about 7 km in range by a lidar beam scanning at 2 in elevation, implying a maximum availability of 0.5, since the maximum range is about double this; this indeed corresponds to the value we observe for a cloud base height of 250 m. This is the case for both summer and autumn. When the cloud layers are at higher altitudes, and not intercepted by the lidar beam, summer months still show reasonable availability, which is not dependent on cloud base height. It is not certain whether the different response seen in cloud base height above 1 km in autumn months is due to a lack of clouds at these heights in our dataset.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f10

Figure 10Fraction of valid measurements as a function of cloud base height from May to September and October to November 2021 for (a, c) lidar and (b, d) X-band radar.

Download

The X-band radar exhibits a more bimodal distribution with respect to cloud base height (Fig. 10b, d). For all cloud base heights shown, the availability of radar measurements is either mostly low, which we attribute to non-precipitating clouds, or mostly high, which we attribute to precipitating clouds. Cloudy cases with intermittent precipitation are presumably responsible for the low cloud situations with valid measurement fractions between 0.1 and 0.9, and these appear to be more common in winter months. The increase in scatter observed during the summer months (Fig. 10b) for cloud bases at higher altitudes is likely due to the presence of insects in the boundary layer below clouds, especially likely below non-precipitating clouds.

Figure 11 shows the measurement availability as a function of range and cloud base height. For Doppler lidar below 500 m in cloud base height, the availability is very dependent on cloud base height, as shown previously. Once the cloud base height is above 500 m, the lidar beam will no longer be impacted by the cloud itself, and there is a less rapid increase in availability with cloud base height. This suggests that there is still some reduction in availability due to the attenuation by precipitation falling from these clouds, the amount of which reduces as the cloud base height increases further.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f11

Figure 11Same as Fig. 6 but as a function of cloud base height separated into different cloud base height classes. Each line corresponds to a cloud base height class with a bin width of 100 m below 500 m and a bin width of 500 m above. The line colour indicates the upper limit of the bin. Note the logarithmic scale for the colour bar.

Download

The performance of the X-band radar again shows different behaviour to the Doppler lidar. For the X-band radar, a low cloud base height results in high availability, with the availability reducing as the cloud base height increases, especially for cloud base heights above 1 km. As for Doppler lidar, this can be attributed to the same reason because the amount of precipitation falling into the radar beam reduces as the cloud base height increases, except that now this reduces rather than increases the data availability. The relative reduction in availability with range is more rapid for high cloud base heights than for low cloud base heights.

4.5 Measurement performance as function of precipitation intensity

Figure 12 shows the distributions of the fraction of valid measurements as a function of precipitation intensity for the Doppler lidar and X-band radar. Only measurements with a corresponding valid precipitation intensity measurement were included, which was approximately 56 % of all measurements. For Doppler lidar (Fig. 12a), the presence of precipitation reduces measurement availability due to attenuation of the lidar beam. However, precipitation is also a strong scatterer at lidar wavelengths, and attenuation is not as rapid as for liquid cloud; hence, the Doppler lidar is still capable of measuring up to a kilometre or more into precipitation – and further if the precipitation is patchy. For the X-band radar (Fig. 12b), as can be expected, the presence of any precipitation will improve the availability.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f12

Figure 12Fraction of valid measurements as a function of precipitation intensity (10 min average; optical) for (a) Doppler lidar and (b) X-band radar Doppler velocity measurements.

Download

Figure 13 shows the effect of precipitation intensity on the measurement availability as a function of range. For Doppler lidar (Fig. 13a), increasing precipitation intensity decreases the availability, but this effect is only evident at ranges beyond 3 km.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f13

Figure 13Same as Fig. 6 but as a function of precipitation intensity separated into different precipitation classes. Each line corresponds to a precipitation intensity class with a bin width of 0.25 mm h−1.

Download

The X-band radar has high availability in precipitation (Fig. 13b), with slight reductions in availability for precipitation intensities less than 0.5 mm h−1. Note that, even for high precipitation intensities, the increase in availability for the longer pulse length is still evident.

Figure 14 shows an example of the Doppler lidar and X-band weather radar measurements from a precipitation case from 17 May 2021. During this time, a thunderstorm was passing over the measurement area. The X-band radar provided excellent measurement coverage inside the precipitation in the eastern half of the area but hardly any measurements in the clear air area in the west. On the contrary, the Doppler wind lidar obtained measurements in the clear air area, but the signal was attenuated quickly inside the heavy rainfall.

https://amt.copernicus.org/articles/15/6507/2022/amt-15-6507-2022-f14

Figure 14Data example from 17 May 2021 at 11:13 UTC, showing PPI scans at 2 of X-band (a) radar reflectivity and (b) Doppler velocity and Doppler lidar (c) carrier-to-noise ratio (CNR) and (d) Doppler velocity. The Doppler velocity colour map is from Crameri (2021).

5 Discussion and conclusions

The aim of this study was to compare Doppler velocities measured by Doppler lidar and X-band weather radar to see, first, how well the co-located measurements agree and, second, how the availability of the measurements varies in different weather conditions. When Doppler velocity measurements are available from both Doppler lidar and X-band weather radar, the agreement is good with R2=0.96. Some departures in agreement were apparent for specific Doppler velocities, indicating that, in some situations, it might be necessary to apply clutter filtering to Doppler lidar measurements at low elevation angles.

The two instruments operate at very different wavelengths and have very different responses to the scattering media present in the atmosphere. Therefore, as expected, the two instruments exhibited different data availabilities. If their different capabilities in specific weather situations can complement one another, then their combination would increase the spatial coverage of wind observations across a wide range of weather conditions.

To better understand the response of each instrument to different meteorological conditions, we investigated their measurement availability in terms of horizontal visibility, cloud base height, and precipitation intensity. In general, Doppler lidar and X-band radar displayed contrasting measurement availabilities in the meteorological conditions examined here. Precipitation, often accompanied by low horizontal visibility and low cloud base height, lead to high radar measurement availability, whereas the Doppler lidar measurement availability was reduced due to the signal being attenuated. Our study showed that the attenuation of the Doppler lidar signal by precipitation was significant but not severe, with reasonably good measurement availability out to at least 2 km, which is in agreement with Oude Nijhuis et al. (2018). Attenuation by cloud, however, is very rapid, and the Doppler lidar does not measure much more than 300 m beyond cloud base.

In contrast, in clear air situations, with high horizontal visibility and high cloud bases, the Doppler lidar is able to measure out to longer ranges. The only hindrance for lidar measurements in clear air is when the atmospheric aerosol concentration falls so low that there is no longer sufficient signal, as in cases with extremely high horizontal visibility. For the radar, there may not be enough scatterers for the radar to measure consistently in clear air. One exception to this is insects, which can increase the radar measurement availability substantially in daytime during summer. This leads to a seasonality in the availability of radar clear air measurements, as seen in Fig. 5, where in summer the availability can be quite high, but in other seasons no measurements are obtained. Note that whether insects are acceptable targets in the radar measurements is location dependent, as using migrating insects as targets can cause a bias in the Doppler velocity measurements compared to the lidar. However, for our location, insects were not found to cause a significant bias in the measurements.

Our measurements show that the conditions where both Doppler lidar and the X-band radar are able to obtain Doppler velocity measurements are the summer clear air situations in which insects and aerosol are present. During our measurement campaign, the optimal horizontal visibility in these cases was approximately 40…60 km and the cloud base height higher than 1 km. However, optimal conditions depend on the measurement location as, for example, horizontal visibility can be decreased by poor air quality, which will also attenuate the Doppler lidar signal more.

It is also of interest is to determine whether there are conditions under which neither instrument is able to measure Doppler velocity and provide winds. For example, conditions with low cloud base height but with no precipitation, including fog, can be responsible for poor data availability from both instruments, especially for measurements at long distances. However, our analysis did not reveal a significant proportion of situations where measurements from both instruments were not available. As the data availability for radar and lidar exhibits contrasting behaviour, it is likely that at least one of the instruments will provide Doppler velocity observations.

Code and data availability

The code used for the analysis is available at https://doi.org/10.5281/zenodo.7234135 (Ritvanen2022). Numerical versions for the data availability curves in Figs. 4, 6, 11, and 13, and the observations from the FMI automatic weather station used in this study can be found at https://doi.org/10.23728/fmi-b2share.8353436be6e444dea496c0c01e6a3658 (Ritvanen et al.2022).

Author contributions

JR, EOC, DM, RL, JT, and LT conceptualised the study. JR did the analysis, led the investigation, developed the methodology and software, and visualized and validated the results. JR and RL curated the data. JR, RL, and JT acquired the funding for the project, and RL and JT managed the project. JT acquired the data storage, and RL and LT acquired the measurement instruments. EOC, DM, RL, JT, and LT supervised the study. JR, EOC, and DM wrote the draft and reviewed and edited the paper.

Competing interests

Raisa Lehtinen and Ludovic Thobois are employees of Vaisala Oyj. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

The authors thank Pekka Rossi, for his input in early stages of the study. We also thank Sebastian Kauczok, for providing feedback that helped improve the article, along with the two reviewers.

Financial support

This research has been supported by the European Space Agency (grant no. 4000132768/20/UK/ND).

Open-access funding was provided by the Helsinki University Library.

Review statement

This paper was edited by Ulla Wandinger and reviewed by Melanie Eickmeier and one anonymous referee.

References

Achtemeier, G. L.: The Use of Insects as Tracers for “Clear-Air” Boundary-Layer Studies by Doppler Radar, J. Atmos. Ocean. Tech., 8, 746–765, https://doi.org/10.1175/1520-0426(1991)008<0746:TUOIAT>2.0.CO;2, 1991. a

Altube, P.: meteocat/vcor_dual_prf, Github [code], https://github.com/meteocat/vcor_dual_prf (last access: 2 November 2022), 2020. a

Andersson, E.: How to Evolve Global Observing Systems, ECMWF newsletter, 153, 37–40, https://doi.org/10.21957/9fxea2, 2017. a

Bharadwaj, N. and Chandrasekar, V.: Wideband Waveform Design Principles for Solid-State Weather Radars, J. Atmos. Ocean. Tech., 29, 14–31, https://doi.org/10.1175/JTECH-D-11-00030.1, 2012. a, b

Contreras, R. F. and Frasier, S. J.: High-Resolution Observations of Insects in the Atmospheric Boundary Layer, J. Atmos. Ocean. Tech., 25, 2176–2187, https://doi.org/10.1175/2008JTECHA1059.1, 2008. a

Crameri, F.: Scientific Colour Maps (7.0.1), Zenodo, https://doi.org/10.5281/zenodo.5501399, 2021. a, b, c

Dolfi-Bouteyre, A., Augère, B., Besson, C., Canat, G., Fleury, D., Gaudo, T., Goular, D., Lombard, L., Planchat, C., Valla, M., Cariou, J.-P., Petilon, O., and Lawson-Daku, J.: 1.5 µm All Fiber Pulsed Lidar For Wake Vortex Monitoring, in: Conference on Lasers and Electro-Optics/Quantum Electronics and Laser Science Conference and Photonic Applications Systems Technologies (2008), San Jose, California, United States, 4–9 May 2008, Paper CMQ3, p. CMQ3, Optica Publishing Group, https://opg.optica.org/abstract.cfm?uri=CLEO-2008-CMQ3 (last access: 2 November 2022), 2008. a, b

Franck, A., Moisseev, D., Vakkari, V., Leskinen, M., Lampilahti, J., Kerminen, V.-M., and O'Connor, E.: Evaluation of convective boundary layer height estimates using radars operating at different frequency bands, Atmos. Meas. Tech., 14, 7341–7353, https://doi.org/10.5194/amt-14-7341-2021, 2021. a

Frehlich, R. G. and Kavaya, M. J.: Coherent laser radar performance for general atmospheric refractive turbulence, Appl. Optics, 30, 5325–5352, https://doi.org/10.1364/AO.30.005325, 1991. a

Guo, J., Zhang, H., and Zhang, X.-j.: Propagating Characteristics of Pulsed Laser in Rain, Int. J. Antenn. Propag., 2015, 1–7, https://doi.org/10.1155/2015/292905, 2015. a

Hannesen, R., Kauczok, S., and Weipert, A.: Quality of Clear-Air Radar Radial Velocity Data: Do Insects Matter?, in: Proc. Eighth European Conf. on Radar in Meteorology and Hydrology, Garmisch-Partenkirchen, Germany, 1–5 September 2014, ERAD, 2014. a

Heistermann, M., Jacobi, S., and Pfaff, T.: Technical Note: An open source library for processing weather radar data (wradlib), Hydrol. Earth Syst. Sci., 17, 863–871, https://doi.org/10.5194/hess-17-863-2013, 2013. a, b

Helmus, J. and Collis, S.: The Python ARM Radar Toolkit (Py-ART), a Library for Working with Weather Radar Data in the Python Programming Language, Journal of Open Research Software, 4, e25, https://doi.org/10.5334/jors.119, 2016. a

Holleman, I. and Beekhuis, H.: Analysis and Correction of Dual PRF Velocity Data, J. Atmos. Ocean. Tech., 20, 443–453, https://doi.org/10.1175/1520-0426(2003)20<443:AACODP>2.0.CO;2, 2003. a

Jatau, P., Melnikov, V., and Yu, T.-Y.: A Machine Learning Approach for Classifying Bird and Insect Radar Echoes with S-Band Polarimetric Weather Radar, J. Atmos. Ocean. Tech., 38, 1797–1812, https://doi.org/10.1175/JTECH-D-20-0180.1, 2021. a

Liu, Z., Barlow, J. F., Chan, P.-W., Fung, J. C. H., Li, Y., Ren, C., Mak, H. W. L., and Ng, E.: A Review of Progress and Applications of Pulsed Doppler Wind LiDARs, Remote. Sens., 11, 2522, https://doi.org/10.3390/rs11212522, 2019. a

Martin, W. J. and Shapiro, A.: Discrimination of Bird and Insect Radar Echoes in Clear Air Using High-Resolution Radars, J. Atmos. Ocean. Techn., 24, 1215–1230, https://doi.org/10.1175/JTECH2038.1, 2007. a

McLaughlin, D., Pepyne, D., Chandrasekar, V., Philips, B., Kurose, J., Zink, M., Droegemeier, K., Cruz-Pol, S., Junyent, F., Brotzge, J., Westbrook, D., Bharadwaj, N., Wang, Y., Lyons, E., Hondl, K., Liu, Y., Knapp, E., Xue, M., Hopf, A., Kloesel, K., DeFonzo, A., Kollias, P., Brewster, K., Contreras, R., Dolan, B., Djaferis, T., Insanic, E., Frasier, S., and Carr, F.: Short-wavelength technology and the potential for distributed networks of small radar systems, B. Am. Meteorol. Soc., 90, 1797–1818, https://doi.org/10.1175/2009BAMS2507.1, 2009. a

Mudukutore, A. S., Chandrasekar, V., and Keeler, R. J.: Pulse compression for weather radars, IEEE T. Geosci. Remote, 36, 125–142, 1998. a

Oude Nijhuis, A. C. P., Thobois, L. P., Barbaresco, F., Haan, S. D., Dolfi-Bouteyre, A., Kovalev, D., Krasnov, O. A., Vanhoenacker-Janvier, D., Wilson, R., and Yarovoy, A. G.: Wind hazard and turbulence monitoring at airports with lidar, radar and mode-S downlinks: the UFO project, B. Am. Meteorol. Soc., 99, 2275–2294, https://doi.org/10.1175/BAMS-D-15-00295.1, 2018. a, b

O'Hora, F. and Bech, J.: Improving weather radar observations using pulse-compression techniques, Meteorol. Appl., 14, 389–401, 2007.  a

Pentikäinen, P., O'Connor, E. J., Manninen, A. J., and Ortiz-Amezcua, P.: Methodology for deriving the telescope focus function and its uncertainty for a heterodyne pulsed Doppler lidar, Atmos. Meas. Tech., 13, 2849–2863, https://doi.org/10.5194/amt-13-2849-2020, 2020. a

Rauber, R. M. and Nesbitt, S. W.: Radar Meteorology: A First Course, Advancing Weather and Climate Science Series, John Wiley & Sons, Hoboken, NJ, USA, 1st edn., ISBN 9781118432655, 2018. a

Rennie, S. J.: Doppler Weather Radar in Australia, Technical Report No. 0.55, CAWCR, ISBN 9780643109483, https://www.cawcr.gov.au/technical-reports/CTR_055.pdf (last access: 2 November 2022), 2012. a

Ritvanen, J.: ritvje/lidar-xband-article-scripts: Article version, Zenodo [code], https://doi.org/10.5281/zenodo.7234135, 2022. a

Ritvanen, J., O'Connor, E., Moisseev, D., Lehtinen, R., Tyynelä, J., and Thobois, L.: Data for Ritvanen et al. (2022) “Complementarity of Wind Measurements from Co-located X-band Weather Radar and Doppler Lidar”, Finnish Meteorological Institute [data set], https://doi.org/10.23728/fmi-b2share.8353436be6e444dea496c0c01e6a3658, 2022. a

Seabold, S. and Perktold, J.: Statsmodels: Econometric and Statistical Modeling with Python, Proceedings of the 9th Python in Science Conference (SciPy 2010), Austin, Texas, 28 June–3 July 2010, 92–96, https://doi.org/10.25080/Majora-92bf1922-011, 2010. a

Siggia, A. D. and Passarelli, R. E.: Gaussian model adaptive processing (GMAP) for improved ground clutter cancellation and moment calculation, in: Proc. Third European Conf. on Radar Meteorology, Visby, Sweden, 6–10 September 2004, ERAD, 67–73, 2004. a

Stepanian, P. M., Horton, K. G., Melnikov, V. M., Zrnić, D. S., and Gauthreaux Jr., S. A.: Dual-Polarization Radar Products for Biological Applications, Ecosphere, 7, e01539, https://doi.org/10.1002/ecs2.1539, 2016. a

Thobois, L., Cariou, J. P., and Gultepe, I.: Review of Lidar-Based Applications for Aviation Weather, Pure Appl. Geophys., 176, 1959–1976, https://doi.org/10.1007/s00024-018-2058-8, 2019. a

Werner, C.: Doppler Wind Lidar, in: Lidar: Range-Resolved Optical Remote Sensing of the Atmosphere, edited by: Weitkamp, C., Springer Series in Optical Sciences, Springer, New York, NY, USA, 325–354, https://doi.org/10.1007/0-387-25101-4_12, 2005. a

Wilson, J. W., Weckwerth, T. M., Vivekanandan, J., Wakimoto, R. M., and Russell, R. W.: Boundary Layer Clear-Air Radar Echoes: Origin of Echoes and Accuracy of Derived Winds, J. Atmos. Ocean. Tech., 11, 1184–1206, https://doi.org/10.1175/1520-0426(1994)011<1184:BLCARE>2.0.CO;2, 1994. a

Zrnic, D. S. and Ryzhkov, A. V.: Polarimetry for Weather Surveillance Radars, B. Am. Meteorol. Soc., 80, 389–406, https://doi.org/10.1175/1520-0477(1999)080<0389:PFWSR>2.0.CO;2, 1999. a

Download
Short summary
Doppler lidars and weather radars provide accurate wind measurements, with Doppler lidar usually performing better in dry weather conditions and weather radar performing better when there is precipitation. Operating both instruments together should therefore improve the overall performance. We investigate how well a co-located Doppler lidar and X-band radar perform with respect to various weather conditions, including changes in horizontal visibility, cloud altitude, and precipitation.