To characterize atmospheric turbulence, the Doppler moments are estimated by weather radars. However, moment accuracy is highly sensitive to radar transmission parameters such as pulse repetition time (T_s) and number of pulses (N_p), which affect Doppler ambiguity and estimation variance. Traditional fixed-parameter radars face trade-offs between aliasing and measurement precision. This paper proposes an adaptive radar framework that dynamically adjusts T_s and N_p on a per-scan basis to improve Doppler moment estimation at a single resolution cell level. Inspired by the Fully Adaptive Radar (FAR) concept, the method also includes a novel multi-lag Doppler width estimation scheme. Results demonstrate enhanced estimation accuracy, enabling better responsiveness to localized and non-stationary weather conditions.

A novel disaggregation algorithm for commercial microwave links (CMLs), named CLEAR (CML Segments with Equal Amounts of Rain), is proposed. CLEAR utilizes a multiplicative random cascade generator to control the splitting of link segments, with the generator's standard deviation dependent on the rain rate and segment length. Spatial consistency during the splitting process is maintained using rain rate information from neighboring CMLs. CLEAR is evaluated on a network of 77 CMLs in Prague. The performance is assessed first using simulated rainfall fields and second through a case study with real attenuation data from the network to demonstrate its applicability in real-world scenarios. Results from the virtual rainfall fields indicate good overall performance, including the generation of realistic spatial patterns and effective estimation of maximal and minimal rain rates along CML paths. The stochastic nature of CLEAR allows it to represent uncertainty as an ensemble of rain rate distributions along CML paths. However, the generated ensembles significantly underestimate overall variability along the paths. Additionally, the case study on real data highlights challenges associated with uncertainties in CML quantitative precipitation estimates, which are common across all methods. In conclusion, CLEAR contributes to generating more representative rainfall distributions along CMLs, which is critical for spatial reconstruction of rainfall fields from path-integrated CML data. It also has the potential to reduce errors in CML quantitative precipitation estimates caused by assuming uniform rain rates along CML paths.

Precipitation nowcasting – the short-term prediction of rainfall using recent radar observations – is critical for weather-sensitive sectors such as transportation, agriculture, and disaster mitigation. While recent deep learning models have shown promise in improving nowcasting skill, most approaches rely solely on 2D radar reflectivity fields, discarding valuable vertical information available in the full 3D radar volume. In this work, we explore the use of echo top height (ETH), a 2D projection indicating the maximum altitude of radar reflectivity above a given threshold, as an auxiliary input variable for deep learning-based nowcasting. We examine the relationship between ETH and radar reflectivity, confirming its relevance for predicting rainfall intensity. We implement a single-pass 3D U-Net that processes both the radar reflectivity and ETH as separate input channels. While our models are able to leverage ETH to improve skill at low rain-rate thresholds, results are inconsistent at higher intensities and the models with ETH systematically underestimate precipitation intensity. Three case studies are used to illustrate how ETH can help in some cases, but also confuse the models and increase the error variance. Nonetheless, the study serves as a foundation for critically assessing the potential contribution of additional variables to nowcasting performance.

Radar rainfall nowcasting has mostly been applied to relatively large (often rural) domains (e.g., river basins), although rainfall nowcasting in small urban areas is expected to be more challenging. Here, we selected 80 events with high rainfall intensities (at least one 1-km ² grid cell experiences precipitation >15 mm h ²¹ for 1-h events or 30 mm day ²¹ for 24-h events) in five urban areas (Maastricht, Eindhoven, The Hague, Amsterdam, and Groningen) in the Netherlands. We evaluated the performance of 9060 probabilistic nowcasts with 20 ensemble members by applying the short-term ensemble prediction system (STEPS) from Pysteps to every 10-min issue time for the selected events. We found that nowcast errors increased with decreasing (urban) areas especially when below 100 km ². In addition, at 30-min lead time, the underestimation of nowcasts was 38% larger and the discrimination ability was 11% lower for 1-h events than for 24-h events. A set of gridded correction factors for the Netherlands, CARROTS (Climatology-based Adjustments for Radar Rainfall in an Operational Setting) could adjust the bias in real-time QPE and nowcasts by 70%. Yet, nowcasts were still found to underestimate rainfall more than 50% above 40-min lead time relative to the reference, which indicates that this error originates from the nowcasting model itself. Also, CARROTS did not adjust the rainfall spatial distribution in urban areas much. In summary, radar-based nowcasting for urban areas (between 67 and 213 km ²) in the Netherlands exhibits a short skillful lead time of about 20 min, which can only be used for last-minute warning and preparation.

Quantifying the magnitude and frequency of extreme precipitation events is key in translating climate observations to planning and engineering design. Past efforts have mostly focused on the estimation of daily extremes using gauge observations. Recent development of high-resolution global precipitation products, now allow estimation of global extremes. This research aims to quantitatively characterize the spatiotemporal behavior of precipitation extremes, by calculating extreme precipitation return levels for multiple durations on the global domain using the Multi-Source Weighted-Ensemble Precipitation (MSWEP) dataset. Both classical and novel extreme value distributions are used to provide insight into the spatial patterns of precipitation extremes. Our results show that the traditional Generalized Extreme Value (GEV) distribution and Peak-Over-Threshold (POT) methods, which only use the largest events to estimate precipitation extremes, are not spatially coherent. The recently developed Metastatistical Extreme Value (MEV) distribution, that includes all precipitation events, leads to smoother spatial patterns of local extremes. For durations of 5 and 10 days, however, there are less events per year to fit the distribution (37 and 22 on average, respectively), leading to larger inter-annual variability and possible overestimation of the extremes. While the GEV and POT methods predict a consistent shift from heavy to thin tails with increasing duration, the MEV method predicts a relatively constant heaviness of the tail for any precipitation duration, opening up an important research question on what is the ‘correct’ tail behavior of extreme precipitation for different durations. The generated extreme precipitation return levels and corresponding parameters are provided as the Global Precipitation EXtremes (GPEX) dataset. These data can be useful for studying the underlying physical processes causing the spatiotemporal variations of the heaviness of extreme precipitation distributions.

Raindrop size distributions (DSDs) play a crucial role in quantitative rainfall estimation using weather radar. Thanks to dual polarization capabilities, crucial information about the DSD in a given volume of air can be retrieved. One popular retrieval method assumes that the DSD can be modeled by a constrained gamma distribution in which the shape (μ) and rate (Λ) parameters are linked together by a deterministic relationship. In the literature, μ-Λ relationships are often taken for granted and applied without much critical discussion. In this study, we take another look at this important issue by conducting a detailed analysis of μ-Λ relations in stratiform rain and quantifying the accuracy of the associated DSD retrievals. Crucial aspects of our research include the sensitivity of μ-Λ relations to the temporal aggregation scale, drop concentration, inter-event variability, and adequacy of the gamma distribution model. Our results show that μ-Λ relationships in stratiform rain are surprisingly robust to the choice of the sampling resolution, sample size, and adequacy of the gamma model. Overall, the retrieved DSDs are in a rather decent agreement with ground observations (correlation coefficient of 0.57 and 0.74 for μ and Dm). The main sources of errors and uncertainty during the retrievals are calibration offsets in reflectivity (Zhh) and differential reflectivity (Zdr). Measurement noise and differences in scale between radars and disdrometers also play a minor role. The raindrop concentration (NT) remains the most difficult parameter to retrieve, which can be off by several orders of magnitude. After careful data filtering and removal of problematic Zhh/Zdr pairs, the correlation coefficient for the retrieved NT values remained low, only slightly increasing from 0.12 into 0.24.

Commercial microwave links (CMLs) in telecommunication networks can provide relevant information for remote sensing of precipitation and other environmental variables, such as path-averaged drop size distribution, evaporation, or humidity. The CoMMon field experiment (COmmercial Microwave links for urban rainfall MONitoring) mainly focused on the rainfall observations by monitoring a 38ĝ€¯GHz dual-polarized CML of 1.85ĝ€¯km path length at a high temporal resolution (4ĝ€¯s), as well as a co-located array of five disdrometers and three rain gauges over 1 year. The dataset is complemented with observations from five nearby weather stations. Raw and pre-processed data, which can be explored with a custom static HTML viewer, are available at 10.5281/zenodo.4923125 . The data quality is generally satisfactory for further analysis, and potentially problematic measurements are flagged to help the analyst identify relevant periods for specific study purposes. Finally, we encourage potential applications and discuss open issues regarding future remote sensing with CMLs.

Conventionally, Micro Rain Radars (MRRs) have been used as a tool to calibrate reflectivity from weather radars, estimate the relation between rainfall rate and reflectivity, and study microphysical processes in precipitation. However, limited attention has been given to the reliability of the retrieved drop size distributions (DSDs) from MRRs. This study sheds more light on this aspect by examining the sensitivity of retrieved DSDs to the assumptions made to map Doppler spectra into size distributions, and investigates the capability of an MRR to assess polarimetric observations from operational weather radars. For that, an MRR was installed near the Cabauw observatory in the Netherlands, between the International Research Center for Telecommunications and Radar (IRCTR) Drizzle Radar (IDRA) X-band radar and the Herwijnen operational C-band radar. The measurements of the MRR from November 2018 to February 2019 were used to retrieve DSDs and simulate horizontal reflectivity Z_e, differential reflectivity Z_DR, and specific differential phase K_DP in rain. Attention is given to the impact of aliased spectra and right-hand-side truncation on the simulation of polarimetric variables. From a quantitative assessment, the correlations of Z_e and Z_DR between the MRR and Herwijnen radar were 0.93 and 0.70, respectively, while those between the MRR and IDRA were 0.91 and 0.69. However, Z_e and Z_DR from the Herwijnen radar showed slight biases of 1.07 and 0.25 dB. For IDRA, the corresponding biases were 2.67 and-0.93 dB. Our results show that MRR measurements are advantageous to inspect the calibration of scanning radars and validate polarimetric estimates in rain, provided that the DSDs are correctly retrieved and controlled for quality assurance.

Spatial downscaling of rainfall fields is a challenging mathematical problem for which many different types of methods have been proposed. One popular solution consists of redistributing rainfall amounts over smaller and smaller scales by means of a discrete multiplicative random cascade (DMRCs). This works well for slowly varying homogeneous rainfall fields but often fails in the presence of intermittency (i.e., large amounts of zero rainfall values). The most common workaround in this case is to use two separate cascade models, namely one for the occurrence and another for the intensity. In this paper, a new and simpler approach based on the notion of equal-volume areas (EVAs) is proposed. Unlike classical cascades where rainfall amounts are redistributed over grid cells of equal size, the EVA cascade splits grid cells into areas of different sizes, with each of them containing exactly half of the original amount of water. The relative areas of the subgrid cells are determined by drawing random values from a logit-normal cascade generator model with scale and intensity-dependent standard deviation (SD). The process ends when the amount of water in each subgrid cell is smaller than a fixed-bucket capacity, at which point the output of the cascade can be resampled over a regular Cartesian mesh. The present paper describes the implementation of the EVA cascade model and gives some first results for 100 selected events in the Netherlands. Performance is assessed by comparing the outputs of the EVA model to bilinear interpolation and to a classical DMRC model based on fixed grid cell sizes. Results show that, on average, the EVA cascade outperforms the classical method, producing fields with more realistic distributions, small-scale extremes and spatial structures. Improvements are mostly credited to the higher robustness of the EVA model in the presence of intermittency and to the lower variance of its generator. However, both approaches have their advantages and weaknesses. For example, while the classical cascade tends to overestimate small-scale variability and extremes, the EVA model tends to produce fields that are slightly too smooth and block shaped compared to the observations. The complementary nature of the two approaches, and the fact that they produce errors of opposite signs, opens up new possibilities for quality control and bias corrections of downscaled fields.

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.

Weather radar has become an invaluable tool for monitoring rainfall and studying its link to hydrological response. However, when it comes to accurately measuring small-scale rainfall extremes responsible for urban flooding, many challenges remain. The most important of them is that radar tends to underestimate rainfall compared to gauges. The hope is that by measuring at higher resolutions and making use of dual-polarization radar, these mismatches can be reduced. Each country has developed its own strategy for addressing this issue. However, since there is no common benchmark, improvements are hard to quantify objectively. This study sheds new light on current performances by conducting a multinational assessment of radar's ability to capture heavy rain events at scales of 5 min up to 2 h. The work is performed within the context of the joint experiment framework of project MUFFIN (Multiscale Urban Flood Forecasting), which aims at better understanding the link between rainfall and urban pluvial flooding across scales. In total, six different radar products in Denmark, the Netherlands, Finland and Sweden were considered. The top 50 events in a 10- year database of radar data were used to quantify the overall agreement between radar and gauges as well as the bias affecting the peaks. Results show that the overall agreement in heavy rain is fair (correlation coefficient 0.7-0.9), with apparent multiplicative biases on the order of 1.2-1.8 (17 %- 44 % underestimation). However, after taking into account the different sampling volumes of radar and gauges, actual biases could be as low as 10 %. Differences in sampling volumes between radar and gauges play an important role in explaining the bias but are hard to quantify precisely due to the many post-processing steps applied to radar. Despite being adjusted for bias by gauges, five out of six radar products still exhibited a clear conditional bias, with intensities of about 1 %-2% per mmh-1. As a result, peak rainfall intensities were severely underestimated (factor 1.8-3.0 or 44 %- 67 %). The most likely reason for this is the use of a fixed Z-R relationship when estimating rainfall rates (R) from reflectivity (Z), which fails to account for natural variations in raindrop size distribution with intensity. Based on our findings, the easiest way to mitigate the bias in times of heavy rain is to perform frequent (e.g., hourly) bias adjustments with the help of rain gauges, as demonstrated by the Dutch C-band product. An even more promising strategy that does not require any gauge adjustments is to estimate rainfall rates using a combination of reflectivity (Z) and differential phase shift (Kdp), as done in the Finnish OSAPOL product. Both approaches lead to approximately similar performances, with an average bias (at 10 min resolution) of about 30% and a peak intensity bias of about 45 %.