ALBERT

Description: The spatial variability of parameters of the raindrop size distribution and its derivatives is investigated through a field study where collocated Particle Size and Velocity (Parsivel2) and two-dimensional video disdrometers were operated at six sites at Wallops Flight Facility, Virginia, from December 2013 to March 2014. The three-parameter exponential function was employed to determine the spatial variability across the study domain where the maximum separation distance was 2.3 km. The nugget parameter of the exponential function was set to 0.99 and the correlation distance d0 and shape parameter s0 were retrieved by minimizing the root-mean-square error, after fitting it to the correlations of physical parameters. Fits were very good for almost all 15 physical parameters. The retrieved d0 and s0 were about 4.5 km and 1.1, respectively, for rain rate (RR) when all 12 disdrometers were reporting rainfall with a rain-rate threshold of 0.1 mm h−1 for 1-min averages. The d0 decreased noticeably when one or more disdrometers were required to report rain. The d0 was considerably different for a number of parameters (e.g., mass-weighted diameter) but was about the same for the other parameters (e.g., RR) when rainfall threshold was reset to 12 and 18 dBZ for Ka- and Ku-band reflectivity, respectively, following the expected Global Precipitation Measurement mission’s spaceborne radar minimum detectable signals. The reduction of the database through elimination of a site did not alter d0 as long as the fit was adequate. The correlations of 5-min rain accumulations were lower when disdrometer observations were simulated for a rain gauge at different bucket sizes.

Description: A network of seven two-dimensional video disdrometers (2DVD), which were operated during the Midlatitude Continental Convective Clouds Experiment (MC3E) in northern Oklahoma, are employed to investigate the spatial variability of raindrop size distribution (DSD) within the footprint of the dual-frequency precipitation radar (DPR) on board the National Aeronautics and Space Administration’s Global Precipitation Measurement (GPM) mission core satellite. One-minute 2DVD DSD observations were interpolated uniformly to 13 points distributed within a nearly circular DPR footprint through an inverse distance weighting method. The presence of deep continental showers was a unique feature of the dataset resulting in a higher mean rain rate R with respect to previous studies. As a measure of spatial variability for the interpolated data, a three-parameter exponential function was applied to paired correlations of three parameters of normalized gamma DSD, R, reflectivity, and attenuation at Ka- and Ku-band frequencies of DPR (Z_Ka, Z_Ku, k_Ka, and k_Ku, respectively). The symmetry of the interpolated sites allowed quantifying the directional differences in correlations at the same distance. The correlation distances d0 of R, k_Ka, and k_Ku were approximately 10 km and were not sensitive to the choice of four rain thresholds used in this study. The d0 of Z_Ku, on the other hand, ranged from 29 to 20 km between different rain thresholds. The coefficient of variation (CV) remained less than 0.5 for most of the samples for a given physical parameter, but a CV of greater than 1.0 was also observed in noticeable samples, especially for the shape parameter and Z_Ku.

Description: One of the main goals of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM) mission is to retrieve parameters of the raindrop size distribution (DSD) globally. As a standard product of the Dual-Frequency Precipitation Radar (DPR) on board the GPM Core Observatory satellite, the mass-weighted mean diameter Dm and the normalized intercept parameter Nw are estimated in three dimensions at the resolution of the radar. These are two parameters of the three-parameter gamma model DSD adopted by the GPM algorithms. This study investigates the accuracy of the Dm retrieval through a comparative study of C-band ground radars (GRs) and GPM products over Italy. The reliability of the ground reference is tested by using two different approaches to estimate Dm. The results show good agreement between the ground-based and spaceborne-derived Dm, with an absolute bias being generally lower than 0.5 mm over land in stratiform precipitation for the DPR algorithm and the combined DPR–GMI algorithm. For the DPR–GMI algorithm, the good agreement extends to convective precipitation as well. Estimates of Dm from the DPR high-sensitivity (HS) Ka-band data show slightly worse results. A sensitivity study indicates that the accuracy of the Dm estimation is independent of the height above surface (not shown) and the distance from the ground radar. On the other hand, a nonuniform precipitation pattern (interpreted both as high variability and as a patchy spatial distribution) within the DPR footprint is usually associated with a significant error in the DPR-derived estimate of Dm.

Description: In this paper, precipitation estimates derived from the Italian ground radar network (IT GR) are used in conjunction with Spinning Enhanced Visible and InfraRed Imager (SEVIRI) measurements to develop an operational oriented algorithm (RAdar INfrared Blending algorithm for Operational Weather monitoring (RAINBOW)) able to provide precipitation pattern and intensity. The algorithm evaluates surface precipitation over five geographical boxes (in which the study area is divided). It is composed of two main modules that exploit a second-degree polynomial relationship between the SEVIRI brightness temperature at 10.8 µm TB10.8 and the precipitation rate estimates from IT GR. These relationships are applied to each acquisition of SEVIRI in order to provide a surface precipitation map. The results, based on a number of case studies, show good performance of RAINBOW when it is compared with ground reference (precipitation rate map from interpolated rain gauge measurements), with high Probability of Detection (POD) and low False Alarm Ratio (FAR) values, especially for light to moderate precipitation range. At the same time, the mean error (ME) values are about 0 mmh−1, while root mean square error (RMSE) is about 2 mmh−1, highlighting a limited variability of the RAINBOW estimations. The precipitation retrievals from RAINBOW have been also compared with the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Support to Operational Hydrology and Water Management (H SAF) official microwave (MW)/infrared (IR) combined product (P-IN-SEVIRI). RAINBOW shows better performances than P-IN-SEVIRI, in terms of both detection and estimates of precipitation fields when they are compared to the ground reference. RAINBOW has been designed as an operational product, to provide complementary information to that of the national radar network where the IT GR coverage is absent, or the quality (expressed in terms of Quality Index (QI)) of the RAINBOW estimates is low. The aim of RAINBOW is to complement the radar and rain gauge network supporting the operational precipitation monitoring.

Description: Precipitation modifies atmospheric column thermodynamics through the process of evaporation and serves as a proxy for latent heat modulation. For this reason, a correct precipitation parameterization (especially for low-intensity precipitation) within global scale models is crucial. In addition to improving our modeling of the hydrological cycle, this will reduce the associated uncertainty of global climate models in correctly forecasting future scenarios, and will enable the application of mitigation strategies. In this manuscript we present a proof of concept algorithm to automatically detect precipitation from lidar measurements obtained from the National Aeronautics and Space Administration Micropulse lidar network (MPLNET). The algorithm, once tested and validated against other remote sensing instruments, will be operationally implemented into the network to deliver a near real time (latency

Description: Coalescence and breakup of drops are recognized as the main mechanisms determining raindrop size distributions on the ground. Full knowledge of these processes is hindered by the challenging difficulties both in the laboratory and tunnel experiments and during observations in the open air. In real rain breakup is mainly due to collision between drops of different sizes (collisional breakup) and occurs when the collisional kinetic energy (CKE) is not absorbed by the colliding drops. In this work, the authors observe and measure the dependence on altitude of the occurrence of collisional breakup in real rainfall events, and then estimate the corresponding limit terminal velocities of drops and their size when breakup significantly takes place. Data from Pludix, an X-band microwave disdrometer, were collected at three locations at different elevations: collisional breakup position in the power spectrum of Pludix increases toward higher frequencies with increasing altitude. Terminal velocities and sizes of the drops at breakup were determined consequently, with drop sizes resulting in 4.55 ± 0.35, 4.02 ± 0.32, and 3.16 ± 0.3 mm for altitudes of 15, 950, and 3300 m MSL, respectively. The authors computed the CKE of the colliding drops at the breakup, finding an upper limiting value of about 1.22 × 10−5 J for all three altitudes. This shows that most dominant collisional breakup signature occurs at similar CKE values for all three locations, corresponding to different drop diameters at different altitudes because of the effect of air density on the drop terminal velocity.

Description: Numerous laboratory and numerical studies have been dedicated to understanding collisional breakup as one of the most important processes in rain formation. The present study aims to identify when, in natural rain, collisional breakup is dominant and thus able to modify the shape of the raindrop size distribution (DSD), up to the equilibrium DSD. To this end, an automated objective algorithm has been developed and applied to a total of more than 6000 two-minute-averaged DSDs. Since breakup is mostly observed in heavy precipitation, the method was applied to the DSDs where rain rate was above 5 mm h−1. The selected breakup DSDs had good agreement with those predicted to be the equilibrium DSD by different theoretical models. The equilibrium DSD was found in a variable fraction of the total samples (0%–7%), confirming that the onset of equilibrium is a rare event in natural rain. The occurrence of a DSD in which breakup is dominant and modifies the DSD but the equilibrium DSD is not reached is higher (15%–47%). The gamma distribution, which has been widely used in the parameterization of observed size spectra, had a poor fitting in breakup-induced DSD, especially in the 1.0–2.6-mm-diameter interval. This can impact applications for which the parameterization of DSD is needed, such as in the retrieval of a DSD integral parameter (such as rain rate) from active remote sensor data.

Description: The National Aeronautics and Space Administration Global Precipitation Measurement (GPM) mission ground validation program uses dual-polarization radar moments to estimate raindrop size distribution (DSD) parameters, the mass-weighted mean drop diameter Dmass, and normalized intercept parameter NW, to validate the GPM Core Observatory–derived DSD parameters. The disdrometer-based Dmass and NW are derived through empirical relationships between Dmass and differential reflectivity ZDR, and between NW, reflectivity ZH, and Dmass. This study employs large datasets collected from two-dimensional video disdrometers (2DVD) during six different field studies to derive the requisite empirical relationships. The uncertainty of the derived Dmass(ZDR) relationship is evaluated through comparisons of 2DVD-calculated and ZDR-estimated Dmass, where ZDR is calculated directly from 2DVD observations. Similarly, the uncertainty of the NW(ZH, Dmass) relationship is evaluated through 2DVD-calculated and Dmass and ZH-estimated NW, where Dmass and ZH are directly calculated from 2DVD observations. This study also presents the sensitivity of Dmass(ZDR) relationships to climate regime and to disdrometer type after developing three additional Dmass(ZDR) relationships from second-generation Particle Size Velocity (PARSIVEL2) disdrometer (P2) observations collected in the Pacific Northwest, in Iowa, and at Kwajalein Atoll in the tropical Pacific Ocean. The application of P2-derived Dmass(ZDR) relationship based on precipitation in the northwestern United States to P2 observations collected over the tropical ocean resulted in the highest error among comparisons of the three datasets.

Description: This study shows how satellite-based passive and active microwave (MW) sensors can be used in conjunction with high-resolution Numerical Weather Prediction (NWP) simulations to provide insights of the precipitation structure of the tropical-like cyclone (TLC) Numa, which occurred on 15–19 November 2017. The goal of the paper is to characterize and monitor the precipitation at the different stages of its evolution from development to TLC phase, throughout the storm transition over the Mediterranean Sea. Observations by the NASA/JAXA Global Precipitation Measurement Core Observatory (GPM-CO) and by the GPM constellation of MW radiometers are used, in conjunction with the Regional Atmospheric Modeling System (RAMS) simulations. The GPM-CO measurements are used to analyze the passive MW radiometric response to the microphysical structure of the storm, while the comparison between successive MW radiometer overpasses shows the evolution of Numa precipitation structure from its early development stage on the Ionian Sea into its TLC phase, as it persists over southern coast of Italy (Apulia region) for several hours. Measurements evidence stronger convective activity at the development phase compared to the TLC phase, when strengthening or weakening phases in the eye development, and the occurrence of warm rain processes in the areas surrounding the eye, are identified. The weak scattering and polarization signal at and above 89 GHz, the lack of scattering signal at 37 GHz, and the absence of electrical activity in correspondence of the rainbands during the TLC phase, indicate weak convection and the presence of supercooled cloud droplets at high levels. RAMS high-resolution simulations support what inferred from the observations, evidencing Numa TLC characteristics (closed circulation around a warm core, low vertical wind shear, intense surface winds, heavy precipitation), persisting for more than 24 h. Moreover, the implementation of DPR 3D reflectivity field in the RAMS data assimilation system shows a small (but non negligible) impact on the precipitation forecast over the sea up to a few hours after the DPR overpass.