This paper describes the results of the research fulfilled in TU-Delft by joint Ukrainian and Dutch team. It analyzes multi-instrument rain observations, using the instrument set, which includes W-band cloud radar, laser optical disdrometers, weather station, and microwave radiometer. New friendly interface software is developed, presented, and used as a tool for comparison and fusion of diverse sensors datasets. The results obtained demonstrate the synergy of multi-instrument measurements and corresponds to the overarching trends of big data analysis. The intricacies of combining data from various sources to enhance calibration and improve the accuracy of atmospheric studies is discussed. In particular, analysis of 94 GHz cloud radar calibration based on disdrometer measurements with application of additional multi-instrument measurements is performed.

@en

This paper is devoted to discussing peculiarities of W-band cloud radar calibration. After a brief overview of meteorological radar calibration methods for quantitative information retrieval, we focus on problems and their possible solutions with respect to mm-wave radar calibration. The experimental part of the research is based on multi-instrument measurements performed during several years in the Cabauw experimental meteorological site in the Netherlands. The accumulated data are used for comparison of 94 GHz radar rain measurements with non-radar droplet size distribution measurements, provided by laser disdrometers. Calculations are done taking into account data of other in situ meteorological measurements. A specialized MATLAB software tool for processing such complex data and radar calibration is developed and demonstrated.

@en

A better understanding of aerosol–cloud interaction processes is important to quantify the role of clouds and aerosols on the climate system. There have been significant efforts to explain the ways aerosols modulate cloud properties. However, from the observational point of view, it is indeed challenging to observe and/or verify some of these processes because no single instrument or platform has been proven to be sufficient. Discrimination between aerosol and cloud is vital for the quantification of aerosol–cloud interaction. With this motivation, a set of observational field campaigns named balloon-borne aerosol–cloud interaction studies (BACIS) is proposed and conducted using balloon-borne in situ measurements in addition to the ground-based (lidar; mesosphere, stratosphere and troposphere (MST) radar; lower atmospheric wind profiler; microwave radiometer; ceilometer) and space-borne (CALIPSO) remote sensing instruments from Gadanki (13.45^◦ N, 79.2^◦ E), India. So far, 15 campaigns have been conducted as a part of BACIS campaigns from 2017 to 2020. This paper presents the concept of the observational approach, lists the major objectives of the campaigns, describes the instruments deployed, and discusses results from selected campaigns. Balloon-borne measurements of aerosol and cloud backscatter ratio and cloud particle count are qualitatively assessed using the range-corrected data from simultaneous observations of ground-based and space-borne lidars. Aerosol and cloud vertical profiles obtained in multi-instrumental observations are found to reasonably agree. Apart from this, balloon-borne profiling is found to provide information on clouds missed by ground-based and/or space-borne lidar. A combination of the Compact Optical Backscatter AerosoL Detector (COBALD) and Cloud Particle Sensor (CPS) sonde is employed for the first time in this study to discriminate cloud and aerosol in an in situ profile. A threshold value of the COBALD colour index (CI) for ice clouds is found to be between 18 and 20, and CI values for coarse-mode aerosol particles range between 11 and 15. Using the data from balloon measurements, the relationship between cloud and aerosol is quantified for the liquid clouds. A statistically significant slope (aerosol–cloud interaction index) of 0.77 found between aerosol backscatter and cloud particle count reveals the role of aerosol in the cloud activation process. In a nutshell, the results presented here demonstrate the observational approach to quantifying aerosol–cloud interactions.

@en

Accurate lidar-based measurements of cloud optical extinction, even though perhaps limited to the cloud base region, are useful. Arguably, more advanced lidar techniques (e.g. Raman) should be applied for this purpose. However, simpler polarisation and backscatter lidars offer a number of practical advantages (e.g. better resolution and more continuous and numerous time series). In this paper, we present a backscatter lidar signal inversion method for the retrieval of the cloud optical extinction in the cloud base region. Though a numerically stable method for inverting lidar signals using a far-end boundary value solution has been demonstrated earlier and may be considered as being well established (i.e. the Klett inversion), the application to high-extinction clouds remains problematic. This is due to the inhomogeneous nature of real clouds, the finite range resolution of many practical lidar systems, and multiple scattering effects. We use an inversion scheme, where a backscatter lidar signal is inverted based on the estimated value of cloud extinction at the far end of the cloud, and apply a correction for multiple scattering within the cloud and a range resolution correction. By applying our technique to the inversion of synthetic lidar data, we show that, for a retrieval of up to 90g m from the cloud base, it is possible to obtain the cloud optical extinction within the cloud with an error better than 5g %. In relative terms, the accuracy of the method is smaller at the cloud base but improves with the range within the cloud until 45g m and deteriorates slightly until reaching 90g m from the cloud base.

@en

In the coming decades, the European energy system is expected to become increasingly reliant on non-dispatchable generation such as wind and solar power. Under such a renewable energy scenario, a better characterization of the extreme weather condition ‘Dunkelflaute’, which can lead to a sustained reduction of wind and solar power, is important. In this paper, we report findings from the very first climatological study of Dunkelflaute events occurring in eleven countries surrounding the North and Baltic Sea areas. By utilizing multi-year meteorological and power production datasets, we have quantified various statistics pertaining to these events and also identified their underlying meteorological drivers. It was found that almost all periods tagged as Dunkelflaute events (with a length of more than 24 h) are in November, December, and January for these countries. On average, there are 50–100 h of such events happening in each of these three months per year. The limited wind and solar power production during Dunkelflaute events is shown to be mainly driven by large-scale high-pressure systems and extensive low-cloud coverage. Even though the possibility of simultaneous Dunkelflaute events in neighboring countries can be as high as 30–40%, such events hardly occur simultaneously in all the eleven countries. Through an interconnected EU-11 power system, the mean frequency of Dunkelflaute drops from 3–9% for the individual countries to approximately 3.5% for the combined region, highlighting the importance of aggregating production over a wide area to better manage the integration of renewable energy generation.@en

In the near future, wind and solar generation are projected to play an increasingly important role in Europe's energy sector. With such fast-growing renewable energy development, the presence of simultaneous calm wind and overcast conditions could cause significant shortfalls in production with potentially serious consequences for system operators. Such events are sometimes dubbed “Dunkelflaute” events and have occurred several times in recent history. The capabilities of contemporary mesoscale models to reliably simulate and/or forecast a Dunkelflaute event are not known in the literature. In this paper, a Dunkelflaute event near the coast of Belgium is simulated utilizing the Weather Research and Forecasting (WRF) model. Comprehensive validation using measured power production data and diverse sets of meteorological data (e.g., floating lidars, radiosondes, and weather stations) indicates the potential of WRF to reproduce and forecast the boundary layer evolution during the event. Extensive sensitivity experiments with respect to grid-size, wind farm parameterization, and forcing datasets provide further insights on the reliability of the WRF model in capturing the Dunkelflaute event.

@en

Radio frequency interference (RFI) has become a growing concern for weather radar, distorting radar variable estimation. By simultaneously or alternately transmitting the horizontal and vertical polarized waves, polarimetric weather radar can be referred to as SHV radar or AHV radar. The SHV radar can mimic the AHV radar by discarding either H- or V-channel measurements, which leads to an alternating scheme. In this research, the real RFI measurements from an operational C-band SHV radar are used to characterize the RFI temporal, spectral, and polarimetric features. Then, the RFI is simulated to quantify the performance of the object-orientated spectral polarimetric (OBSPol) filter in RFI mitigation. The OBSPol filter has been previously proposed by the authors to mitigate the narrowband clutter (both stationary and moving) and noise. This work extends the application of the filter to remove the RFI for SHV radar. Specifically, by taking advantage of the low copolar correlation of the RFI signal measured in AHV radar, the RFI mitigation method is designed, and its effectiveness is proven by qualitative and quantitative analyses. In particular, in the case of RFI overlapped to weather echoes in the time domain, the RFI can be mitigated, also when the duty cycle of the RFI is high. However, this work does not provide a full evaluation of the RFI mitigation performance on all radar data outputs but a proof of concept to show the effectiveness of the proposed filter for RFI mitigation.

@en

In the coming decades, both wind and solar power production will be playing increasingly important roles in Europe's energy economy. It is absolutely essential that power grids are resilient against any unusual weather phenomena. One such meteorological phenomenon, "Dunkelflaute", is causing serious concern to the renewable energy industry, which is primarily characterized by calm winds and overcast conditions. For example, a Dunkelflaute event happened in the Netherlands on 30th April 2018 leading to a significant shortfall in renewable energy generation requiring emergency intervention by the system operator. By analyzing this case, this paper investigates the performance of a state-of-the-art mesoscale model, Weather Research and Forecasting (WRF), in forecasting a Dunkelflaute event. Multiple WRF simulations are driven using real-time Global Forecast System (GFS) operational data over a range of prediction horizons. For comparison, a benchmark run is carried out using ERA5 reanalysis data as boundary conditions. Through validation using a variety of measured data covering onshore and offshore areas, wind speed is shown to be more predictable than cloud-cover in this particular case study.

@en

The adequacy of the gamma model to describe the variability of raindrop size distributions (DSD) is studied using observations from an optical disdrometer. Model adequacy is checked using a combination of Kolmogorov–Smirnov goodness-of-fit test and Kullback–Leibler divergence and the sensitivity of the results to the sampling resolution is inves-tigated. A new adaptive DSD sampling technique capable of determining the highest possible temporal sampling resolution at which the gamma model provides an adequate representation of sampled DSDs is proposed. The results show that most DSDs at 30 s are not strictly distributed according to a gamma model, while at the same time they are not far away from it either. According to the adaptive DSD sampling algorithm, the gamma model proves to be an adequate choice for the majority (85.81%) of the DSD spectra at resolutions up to 300 s. At the same time, it also reveals a considerable number of DSD spectra (5.55%) that do not follow a gamma distribution at any resolution (up to 1800 s). These are attributed to transitional periods during which the DSD is not stationary and exhibits a bimodal shape that cannot be modeled by a gamma distribution. The proposed resampling procedure is capable of automatically identifying and flagging these periods, providing new valuable quality control mechanisms for DSD retrievals in disdrometers and weather radars.

@en

In this article, five velocity-based energy dissipation rate (EDR) retrieval techniques are assessed. The EDR retrieval techniques are applied to Doppler measurements from Transportable Atmospheric Radar (TARA)—a precipitation profiling radar—operating in the vertically fixed-pointing mode. A generalized formula for the Kolmogorov constant is derived, which gives potential for the application of the EDR retrieval techniques to any radar line of sight (LOS). Two case studies are discussed that contain rain events of about 2 and 18 h, respectively. The EDR values retrieved from the radar are compared to in situ EDR values from collocated sonic anemometers. For the two case studies, a correlation coefficient of 0.79 was found for the wind speed variance (WSV) EDR retrieval technique, which uses 3D wind vectors as input and has a total sampling time of 10 min. From this comparison it is concluded that the radar is able to measure EDR with a reasonable accuracy. Almost no correlation was found for the vertical wind velocity variance (VWVV) EDR retrieval technique, as it was not possible to sufficiently separate the turbulence dynamics contribution to the radar Doppler mean velocities from the velocity contribution of falling raindrops. An important cause of the discrepancies between radar and in situ EDR values is thus due to insufficient accurate estimation of vertical air velocities.@en

Weather radar is well recognized as an effective sensor for obtaining the microphysical and dynamical properties of precipitation at high spatial and temporal resolution. Radar calibration is one of the most important prerequisites for achieving accurate observations. In this article, a portable, cost-effective and repeatable radar calibration technique, namely, unmanned aerial vehicle (UAV)-aided radar calibration, is proposed. A UAV serves as the stable aerial platform carrying a metal sphere, flying over the radar illumination areas to complete the calibration process. The flying routine of the UAV can be pre-programmed, and thus, the antenna pattern regarding different elevation and azimuth angles can be retrieved. To obtain the position of the sphere, the real-time single-frequency precise point positioning-type global navigation satellite system solution is developed. In addition, the radar constant is calculated in the range-Doppler domain, and only the data where the metal sphere separates from clutter and other objects are selected. The S-band polarimetric Doppler transportable atmospheric radar (TARA) is used in the calibration campaign. The experiments demonstrate the following results: 1) antenna pointing calibration can be completed and 2) antenna pattern can be retrieved and weather radar constant can be accurately calculated.@en

A Sentinel-5p/TROPOMI validation campaign was held in the Netherlands based at the Cabauw Experimental Site for Atmospheric Research during September 2019. The TROpomi vaLIdation eXperiment (TROLIX) consisted of active and passive remote sensing platforms in conjunction with several balloon-borne, airborne and surface chemical measurements. The goal of this geophysical validation campaign was to make intensive observations to establish the quality of TROPOMI L2 main data products (UVAI, Aerosol Layer Height, NO2, O3, HCHO, Clouds) under realistic non-idealized conditions with varying cloud cover and a wide range of atmospheric conditions. Since TROPOMI is a hyperspectral imager with a very high spatial resolution of 3.5x7 km2, understanding local effects such as inhomogeneous sources of pollution, sub-pixel clouds and variations in ground albedo is important to interpret TROPOMI results. Therefore, the campaign included sub-pixel resolution local networks of sensors, involving Pandora and MAXDOAS instruments, around Cabauw (51.97° N, 4.93° E) and within the city of Rotterdam. Cabauw is considered rural while Rotterdam is densely populated and industrialized. These focal areas were connected through airborne as well as ground based mobile observations. Cabauw, using its comprehensive in-situ and remote sensing observation program in and around the 213 m meteorological tower, was the main site of the campaign with focus on vertical profiling using lidar instruments for aerosols, clouds, water vapor, tropospheric and stratospheric ozone, as well as balloon-borne sensors for NO2 and ozone. The data set collected can be directly compared to the TROPOMI L2 data products, while measurements of parameters related to a-priori data and auxiliary parameters that influence the quality of the L2 products such as aerosol and cloud profiles and in-situ aerosol and atmospheric chemistry were also collected. This paper gives an overview of the campaign, and an overview of the participating main and ancillary instrumentation. Furthermore, an overview of the meteorological and atmospheric conditions observed during the campaign is given from the respective perspectives of the participating instruments, including satellite observations and the support by atmospheric modeling (CAMS). @en

The growth of ice crystals in presence of supercooled liquid droplets represents the most important process for precipitation formation in the mid-latitudes. However, such mixed-phase interaction processes remain relatively unknown, as capturing the complexity in cloud dynamics and microphysical variabilities turns to be a real observational challenge. Ground-based radar systems equipped with fully polarimetric and Doppler capabilities in high temporal and spatial resolutions such as the S-band transportable atmospheric radar (TARA) are best suited to observe mixed-phase growth processes. In this paper, measurements are taken with the TARA radar during the ACCEPT campaign (analysis of the composition of clouds with extended polarization techniques). Besides the common radar observables, the 3-D wind field is also retrieved due to TARA unique three beam configuration. The novelty of this paper is to combine all these observations with a particle evolution detection algorithm based on a new fall streak retrieval technique in order to study ice particle growth within complex precipitating mixed-phased cloud systems. In the presented cases, three different growth processes of ice crystals, plate-like crystals, and needles are detected and related to the presence of supercooled liquid water. Moreover, TARA observed signatures are assessed with co-located measurements obtained from a cloud radar and radiosondes. This paper shows that it is possible to observe ice particle growth processes within complex systems taking advantage of adequate technology and state of the art retrieval algorithms. A significant improvement is made towards a conclusive interpretation of ice particle growth processes and their contribution to rain production using fall streak rearranged radar data.

@en

Long-term measurements of PM_2.5 mass concentrations and aerosol particle size distributions from 2008 to 2015, as well as hygroscopicity measurements conducted over one year (2008–2009) at Cabauw, The Netherlands, are compiled here in order to provide a comprehensive dataset for understanding the trends and annual variabilities of the atmospheric aerosol in the region. PM_2.5 concentrations have a mean value of 14.4 μg m^-3 with standard deviation 2.1 μg m^-3, and exhibit an overall decreasing trend of −0.74 μg m^-3 year^-1. The highest values are observed in winter and spring and are associated with a shallower boundary layer and lower precipitation, respectively, compared to the rest of the seasons. Number concentrations of particles smaller than 500 nm have a mean of 9.2 × 10³particles cm^-3 and standard deviation 4.9 × 10³particles cm^-3, exhibiting an increasing trend between 2008 and 2011 and a decreasing trend from 2013 to 2015. The particle number concentrations exhibit highest values in spring and summer (despite the increased precipitation) due to the high occurrence of nucleation-mode particles, which most likely are formed elsewhere and are transported to the observation station. Particle hygroscopicity measurements show that, independently of the air mass origin, the particles are mostly externally mixed with the more hydrophobic mode having a mean hygroscopic parameter κ of 0.1 while for the more hydrophilic mode κ is 0.35. The hygroscopicity of the smaller particles investigated in this work (i.e., particles having diameters of 35 nm) appears to increase during the course of the nucleation events, reflecting a change in the chemical composition of the particles.

@en

The cloud droplet number concentration (N_d) is of central interest to improve the understanding of cloud physics and for quantifying the effective radiative forcing by aerosol-cloud interactions. Current standard satellite retrievals do not operationally provide N_d, but it can be inferred from retrievals of cloud optical depth (τ_c) cloud droplet effective radius (r_e) and cloud top temperature. This review summarizes issues with this approach and quantifies uncertainties. A total relative uncertainty of 78% is inferred for pixel-level retrievals for relatively homogeneous, optically thick and unobscured stratiform clouds with favorable viewing geometry. The uncertainty is even greater if these conditions are not met. For averages over 1° ×1° regions the uncertainty is reduced to 54% assuming random errors for instrument uncertainties. In contrast, the few evaluation studies against reference in situ observations suggest much better accuracy with little variability in the bias. More such studies are required for a better error characterization. N_d uncertainty is dominated by errors in r_e, and therefore, improvements in r_e retrievals would greatly improve the quality of the N_d retrievals. Recommendations are made for how this might be achieved. Some existing N_d data sets are compared and discussed, and best practices for the use of N_d data from current passive instruments (e.g., filtering criteria) are recommended. Emerging alternative N_d estimates are also considered. First, new ideas to use additional information from existing and upcoming spaceborne instruments are discussed, and second, approaches using high-quality ground-based observations are examined.

@en

In situ measurements using unmanned aerial vehicles (UAVs) and remote sensing observations can independently provide dense vertically resolved measurements of atmospheric aerosols, information which is strongly required in climate models. In both cases, inverting the recorded signals to useful information requires assumptions and constraints, and this can make the comparison of the results difficult. Here we compare, for the first time, vertical profiles of the aerosol mass concentration derived from light detection and ranging (lidar) observations and in situ measurements using an optical particle counter on board a UAV during moderate and weak Saharan dust episodes. Agreement between the two measurement methods was within experimental uncertainty for the coarse mode (i.e. particles having radii > 0.5ĝ€μm), where the properties of dust particles can be assumed with good accuracy. This result proves that the two techniques can be used interchangeably for determining the vertical profiles of aerosol concentrations, bringing them a step closer towards their systematic exploitation in climate models.

@en

One of the most beneficial polarimetric variables may be the specific differential phase K_DP because of its independence from power attenuation and radar miscalibration. However, conventional K_DP estimation requires a substantial amount of range smoothing as a result of the noisy characteristic of the measured differential phase ψ_DP. In addition, the backscatter differential phase δ_hv component of ψ_DP, significant at C- and X-band frequency, may lead to inaccurate K_DP estimates. In this work, an adaptive approach is proposed to obtain accurate K_DP estimates in rain from noisy ψ_DP, whose δ_hv is of significance, at range resolution scales. This approach uses existing relations between polarimetric variables in rain to filter δ_hv from ψ_DP while maintaining its spatial variability. In addition, the standard deviation of the proposed K_DP estimator is mathematically formulated for quality control. The adaptive approach is assessed using four storm events, associated with light and heavy rain, observed by a polarimetric X-band weather radar in the Netherlands. It is shown that this approach is able to retain the spatial variability of the storms at scales of the range resolution. Moreover, the performance of the proposed approach is compared with two different methods. The results of this comparison show that the proposed approach outperforms the other two methods in terms of the correlation between K_DP and reflectivity, and K_DP standard deviation reduction.

@en

In radar polarimetry, the differential phase CDP consists of the propagation differential phase FDP and the backscatter differential phase dhv. While FDP is commonly used for attenuation correction (i.e., estimation of the specific attenuation A and specific differential phase KDP), recent studies have demonstrated that dhv can provide information concerning the dominant size of raindrops. However, the estimation of FDP and dhv is not straightforward given their coupled nature and the noisy behavior of CDP, especially over short paths. In this work, the impacts of estimating FDP on the estimation of A over short paths, using the extended version of the ZPHI method, are examined. Special attention is given to the optimization of the parameter a that connects KDP and A. In addition, an improved technique is proposed to compute dhv from CDP and FDP in rain. For these purposes, diverse storm events observed by a polarimetric X-band radar in the Netherlands are used. Statistical analysis based on the minimum errors associated with the optimization of a and the consistency between KDP and A showed that more accurate and stable a and A are obtained if FDP is estimated at range resolution, which is not possible by conventional range filtering techniques. Accurate dhv estimates were able to depict the spatial variability of dominant raindrop size in the observed storms. By following the presented study, the ZPHI method and its variations can be employed without the need for considering long paths, leading to localized and accurate estimation of A and dhv.@en

A novel algorithm is put forward to separate radar target and moving clutter based on a combination of the low-rank matrix optimization (LRMO) and the decision tree. Similar to the moving object detection in the field of automated video analysis, the proposed separation method, which is carried out in the range-Doppler domain, makes use of different motion variations of radar target and clutter in the spectrogram sequence. The technique is very general, but the focus of this paper is on narrowband moving clutter suppression in a weather radar. The first step in implementing this method is the generation of a range-Doppler spectrogram sequence. Then, the LRMO is applied to the obtained sequence to divide target and moving clutter into foreground and background. From the foreground sequence which is obtained by solving the LRMO, foreground frequency and spectral width are combined in a decision tree to obtain a filtering mask to mitigate the narrowband moving clutter and noise. Data collected by a polarimetric Doppler weather radar known as the IRCTR Drizzle Radar are used to validate the performance of the proposed algorithm. Moreover, its effectiveness in removing narrowband moving clutter is quantitatively assessed through comparisons with another clutter mitigation method.

@en