Relations for retrieving precipitation and attenuation information from radar measurements play a key role in radar meteorology. The uncertainty in such relations highly affects the precipitation and attenuation estimates. Weather radar algorithms are often derived by applying regression methods to precipitation measurements and radar observables simulated from datasets of drop size distributions (DSD) using microphysical and electromagnetic assumptions. DSD datasets can be derived from theoretical considerations or obtained from experimental measurements collected throughout the years by disdrometers. Although the relations obtained from experimental disdrometer datasets can be generally considered more representative of a specific climatology, the measuring errors, which depend on the specific type of disdrometer used, introduce an element of uncertainty to the final retrieval algorithms. Eventually, data quality checks and filtering procedures applied to disdrometer measurements play an important role. In this study, we pursue two main goals: (i) evaluate two different techniques for establishing weather radar algorithms from measured DSD, and (ii) investigate to what extent dual-polarization radar algorithms derived from experimental DSD datasets are influenced by the different error structures introduced by the various disdrometer types (namely 2D video disdrometer, first and second generation of OTT Parsivel disdrometer, and Thies Clima disdrometer) used to collect the data. Furthermore, weather radar algorithms optimized for Italian climatology are presented and discussed.

A multi-sensor analysis of convective precipitation events that occurred in central Italy in autumn 2012 during the HyMeX (Hydrological cycle in the Mediterranean experiment) Special Observation Period (SOP) 1.1 is presented. Various microphysical properties of liquid and solid hydrometeors are examined to assess their relationship with lightning activity. The instrumentation used consisted of a C-band dual-polarization weather radar, a 2-D video disdrometer, and the LINET lightning network. Results of T-matrix simulation for graupel were used to (i) tune a fuzzy logic hydrometeor classification algorithm based on Liu and Chandrasekar (2000) for the detection of graupel from C-band dual-polarization radar measurements and (ii) to retrieve graupel ice water content. Graupel mass from radar measurements was related to lightning activity. Three significant case studies were analyzed and linear relations between the total mass of graupel and number of LINET strokes were found with different slopes depending on the nature of the convective event (such as updraft strength and freezing level height) and the radar observational geometry. A high coefficient of determination (R² = 0.856) and a slope in agreement with satellite measurements and model results for one of the case studies (15 October 2012) were found. Results confirm that one of the key features in the electrical charging of convective clouds is the ice content, although it is not the only one. Parameters of the gamma raindrop size distribution measured by a 2-D video disdrometer revealed the transition from a convective to a stratiform regime. The raindrop size spectra measured by a 2-D video disdrometer were used to partition rain into stratiform and convective classes. These results are further analyzed in relation to radar measurements and to the number of strokes. Lightning activity was not always recorded when the precipitation regime was classified as convective rain. High statistical scores were found for relationships relating lightning activity to graupel aloft.

Snow plays a crucial role in the hydrological cycle and energy budget of the Earth, and remote sensing instruments with the necessary spatial coverage, resolution, and temporal sampling are essential for snowfall monitoring. Among such instruments, ground-radars have scanning capability and a resolution that make it possible to obtain a 3D structure of precipitating systems or vertical profiles when used in profiling mode. Radars from space have a lower spatial resolution, but they provide a global view. However, radar-based quantitative estimates of solid precipitation are still a challenge due to the variability of the microphysical, geometrical, and electrical features of snow particles. Estimations of snowfall rate are usually accomplished using empirical, long-term relationships between the equivalent radar reflectivity factor (Ze) and the liquid-equivalent snowfall rate (SR). Nevertheless, very few relationships take advantage of the direct estimation of the microphysical characteristics of snowflakes. In this work, we used a K-band vertically pointing radar collocated with a laser disdrometer to develop Ze-SR relationships as a function of snow classification. The two instruments were located at the Italian Antarctic Station Mario Zucchelli. The K-band radar probes the low-level atmospheric layers, recording power spectra at 32 vertical range gates. It was set at a high vertical resolution (35 m), with the first trusted range gate at a height of only 100 m. The disdrometer was able to provide information on the particle size distribution just below the trusted radar gate. Snow particles were classified into six categories (aggregate, dendrite aggregate, plate aggregate, pristine, dendrite pristine, plate pristine). The method was applied to the snowfall events of the Antarctic summer seasons of 2018–2019 and 2019–2020, with a total of 23,566 min of precipitation, 15.3% of which was recognized as showing aggregate features, 33.3% dendrite aggregate, 7.3% plates aggregate, 12.5% pristine, 24% dendrite pristine, and 7.6% plate pristine. Applying the appropriate Ze-SR relationship in each snow category, we calculated a total of 87 mm water equivalent, differing from the total found by applying a unique Ze-SR. Our estimates were also benchmarked against a colocated Alter-shielded weighing gauge, resulting in a difference of 3% in the analyzed periods.

This work proposes a multi-parameter method for the detection of cloud-to-ground stroke rate (SR_CG) associated to convective cells, based on the measurements of a low-cost single-polarization X-band weather radar. To train and test our procedure, we built up a multi-year dataset, collecting 1575 radar reflectivity volumes that were acquired in the pilot study area of Naples metropolitan environment matched with the LIghtning NETwork (LINET) strokes and meteorological in-situ data. Three radar-based variables are extracted simultaneously for each rain cell and properly merged together, using “ad hoc” classification methods, to produce an estimation of the expected lightning activity for each rain cell. These variables, proxies of mixed-phase particles and ice amount into a convective cell, are combined into a single label to cluster the SR_CG into two categories: SR_CG = 0 (no production of strokes) or SR_CG > 0 (stroke production), respectively. Overall, the main results are comparable with those that were obtained from more advanced radar systems, showing a Critical Success Index of 0.53, an Equitable Threat Score of 0.34, a Frequency Bias Index of 1.00, a Heidke Skill Score of 0.42, a Hanssen-Kuiper Skill Score of 0.42, and an area under the curve of probability of detection as a function of false alarm rate (usually referred as ROC curve) equal to 0.78. The developed technique, although with some limitations, outperforms those based on the use of single stroke proxy parameters.