Search

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter Dmass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM) mission Core Observatory satellite-based retrieved DSD. Both Dmass and NW are obtained from radar-measured reflectivity ZH and differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of ZH and ZDR and derived Dmass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, Dmass, and log(NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOL ZH–ZDR pairs were well outside the envelope of 2DVD-estimated ZH–ZDR pairs. The exclusion of these pairs improved the comparisons noticeably.

In this article, we report the first investigation over time of the atmospheric conditions around terrestrial gamma-ray flash (TGF) occurrences, using GPS sensors in combination with geostationary satellite observations and ERA5 reanalysis data. The goal is to understand which characteristics are favorable to the development of these events and to investigate if any precursor signals can be expected. A total of 9 TGFs, occurring at a distance lower than 45 km from a GPS sensor, were analyzed and two of them are shown here as an example analysis. Moreover, the lightning activity, collected by the World Wide Lightning Location Network (WWLLN), was used in order to identify any links and correlations with TGF occurrence and precipitable water vapor (PWV) trends. The combined use of GPS and the stroke rate trends identified, for all cases, a recurring pattern in which an increase in PWV is observed on a timescale of about two hours before the TGF occurrence that can be placed within the lightning peak. The temporal relation between the PWV trend and TGF occurrence is strictly related to the position of GPS sensors in relation to TGF coordinates. The life cycle of these storms observed by geostationary sensors described TGF-producing clouds as intense with a wide range of extensions and, in all cases, the TGF is located at the edge of the convective cell. Furthermore, the satellite data provide an added value in associating the GPS water vapor trend to the convective cell generating the TGF. The investigation with ERA5 reanalysis data showed that TGFs mainly occur in convective environments with unexceptional values with respect to the monthly average value of parameters measured at the same location. Moreover, the analysis showed the strong potential of the use of GPS data for the troposphere characterization in areas with complex territorial morphologies. This study provides indications on the dynamics of con-vective systems linked to TGFs and will certainly help refine our understanding of their production, as well as highlighting a potential approach through the use of GPS data to explore the lightning activity trend and TGF occurrences. publishedVersion

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.

Since early March 2014, the NASA/JAXA Global Precipitation Measurement Core- Observatory (GPM-CO) satellite has allowed analysis of precipitation systems around the globe, thanks to the capabilities of the GPM Microwave Imager (GMI) and Dual-Frequency Precipitation Radar (DPR). In this work, we demonstrate how GPM-CO measurements obtained from 4 years of observations over the Mediterranean area can be used as an extremely effective tool to study the main climatological characteristics of the most intense Mediterranean storm structures. DPR and GMI-based Precipitation Features (PFs) parameters are used as proxies of the vertical structure and microphysical properties of these events, and their statistical distribution is analyzed to identify extremes. The analysis of annual, seasonal and geographical distribution of the identified deep convective systems highlights substantial differences in their diurnal cycle and in the distribution between land-sea and summer-winter. There is a general shift of the convective systems from the south (mostly over the sea) in the cold season, to the north (mostly over land) in the warm season. The analysis shows also that the inferred convective intensity is not always related to heavy precipitation. Known DPR and GMI-based criteria were adopted to identify overshooting top events and potential hailstorms, identify extreme deep convection signatures, like those observed for tropical and subtropical systems, and the most intense occur mostly over the sea. Although the analysis is limited to four years, the results show that the GPM-CO offers unprecedented measurements to identify and characterize extreme weather events in the Mediterranean region, with unique potentials for future long-term climatology and interannual variability analysis.

Lightning data assimilation (LDA) is a powerful tool to improve the weather forecast of convective events and has been widely applied with this purpose in the past two decades. Most of these applications refer to events hitting coastal and land areas, where people live. However, a weather forecast over the sea has many important practical applications, and this paper focuses on the impact of LDA on the precipitation forecast over the central Mediterranean Sea around Italy. The 3 h rapid update cycle (RUC) configuration of the weather research and forecasting (WRF) model) has been used to simulate the whole month of November 2019. Two sets of forecasts have been considered: CTRL, without lightning data assimilation, and LIGHT, which assimilates data from the LIghtning detection NETwork (LINET). The 3 h precipitation forecast has been compared with observations of the Integrated Multi-satellitE Retrievals for Global Precipitation Mission (GPM) (IMERG) dataset and with rain gauge observations recorded in six small Italian islands. The comparison of CTRL and LIGHT precipitation forecasts with the IMERG dataset shows a positive impact of LDA. The correlation between predicted and observed precipitation improves over wide areas of the Ionian and Adriatic Seas when LDA is applied. Specifically, the correlation coefficient for the whole domain increases from 0.59 to 0.67, and the anomaly correlation (AC) improves by 5% over land and by 8% over the sea when lightning is assimilated. The impact of LDA on the 3 h precipitation forecast over six small islands is also positive. LDA improves the forecast by both decreasing the false alarms and increasing the hits of the precipitation forecast, although with variability among the islands. The case study of 12 November 2019 (time interval 00–03 UTC) has been used to show how important the impact of LDA can be in practice. In particular, the shifting of the main precipitation pattern from land to the sea caused by LDA gives a much better representation of the precipitation field observed by the IMERG precipitation product.

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.

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.

Both numerical modeling and laboratory experiments document the possibility of a raindrop size distribution (DSD) to evolve to an equilibrium stage (EDSD), where all the principal processes occur at steady rates. The aim of this work is to observe the temporal behavior of the DSD and to directly investigate the conditions favorable to the onset of the EDSD in natural rain. We exploited a large disdrometer dataset collected in the framework of the Ground Validation activities related to the NASA Global Precipitation Measurement mission. More than 200,000 one-minute data of two-dimensional video disdrometer (2DVD) are collected over USA to represent a wide range of precipitation types. The original data are averaged over 2 min and an automatic algorithm is used on a selected subset to identify samples with EDSD. Results show that the EDSD occurs mainly in convective events and lasts for very short time intervals (2 to 4 min). It is more frequent for rain rate between 20 and 40 mm h − 1 and it mostly occurs during sharp increase of precipitation rates.

Continuous estimates of the vertical integrated precipitable water vapor content from the tropospheric delay of the signal received by the antennas of the global positioning system (GPS) are used in this paper, in conjunction with the measurements of the Meteosat Second Generation (MSG) spinning enhanced visible and infrared imager (SEVIRI) radiometer and with the lightning activity, collected here by the ground-based lightning detection network (LINET), in order to identify links and recurrent patterns useful for improving nowcasting applications. The analysis of a couple of events is shown here as an example of more general behavior. Clear signs appear before the peak of lightning activity on a timescale from 2 to 3 h. In particular, the lightning activity is generally preceded by a period in which the difference between SEVIRI brightness temperature (TB) at channel 5 and channel 6 (i.e., ∆TB) presents quite constant values around 0 K. This trend is accompanied by an increase in precipitable water vapor (PWV) values, reaching a maximum in conjunction with the major flash activity. The results shown in this paper evidence good potentials of using radiometer and GPS measurements together for predicting the abrupt intensification of lightning activity in nowcasting systems.

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.

In this paper we illustrate a new, simple and complementary ground-based methodology to retrieve the vertically resolved atmospheric precipitation intensity through a synergy between measurements from the National Aeronautics and Space Administration (NASA) Micropulse Lidar network (MPLNET), an analytical model solution and ground-based disdrometer measurements. The presented results are obtained at two mid-latitude MPLNET permanent observational sites, located respectively at NASA Goddard Space Flight Center, USA, and at the Universitat Politècnica de Catalunya, Barcelona, Spain. The methodology is suitable to be applied to existing and/or future lidar/ceilometer networks with the main objective of either providing near-real time (3h latency) rainfall intensity measurements and/or to validate satellite missions, especially for critical light precipitation (

In this paper, we illustrate a new, simple and complementary ground-based methodology to retrieve the vertically resolved atmospheric precipitation intensity through a synergy between measurements from the National Aeronautics and Space Administration (NASA) Micropulse Lidar network (MPLNET), an analytical model solution and ground-based disdrometer measurements. The presented results are obtained at two mid-latitude MPLNET permanent observational sites, located respectively at NASA Goddard Space Flight Center, USA, and at the Universitat Politècnica de Catalunya, Barcelona, Spain. The methodology is suitable to be applied to existing and/or future lidar/ceilometer networks with the main objective of either providing near real-time (3 h latency) rainfall intensity measurements and/or to validate satellite missions, especially for critical light precipitation (

The Ka–Ku Dual-Frequency Precipitation Radar (DPR) and the Microwave Imager on board the Global Precipitation Measurement (GPM) mission core satellite have been collecting data for more than 3 years, providing precipitation products over the globe, including oceans and remote areas where ground-based precipitation measurements are not available. The main objective of this work is to validate the GPM-DPR products over a key climatic region with complex orography such as the Italian territory. The performances of the DPR precipitation rate products are evaluated over an 18-month period (July 2015–December 2016) using both radar and rain gauge data. The ground reference network is composed of 22 weather radars and more than 3000 rain gauges. DPR dual-frequency products generally show better performance with respect to the single-frequency (i.e., Ka- or Ku-band only) products, especially when ground radar data are taken as reference. A sensitivity analysis with respect to season and rainfall intensity is also carried out. It was found that the normal scan (NS) product outperforms the high-sensitivity scan (HS) and matched scan (MS) during the summer season. A deeper analysis is carried out to investigate the larger discrepancies between the DPR-NS product and ground reference data. The most relevant improvement of the DPR products’ performance was found by limiting the comparison to the upscaled radar data with a higher quality index. The resulting scores in comparison with ground radars are mean error (ME) = −0.44 mm h−1, RMSE = 3.57 mm h−1, and fractional standard error (FSE) = 142%, with the POD = 65% and FAR = 1% for rainfall above 0.5 mm h−1.

This paper presents the results of a recent flight campaign conducted over the Great Lakes region and reports the first observations of the W-band normalized backscattered cross section ( $\sigma _{0}$ ) for V and H polarization and the linear depolarization ratios (LDRs) from different types of surfaces at moderate incidence angles ( $\sigma _{0}$ behaves as previously reported at small incidence angles, it features a marked decrease with increasing incidence angles between 20° and 50°. There is a strong dependence of normalized backscattered cross sections both on the wind speed and on the wind direction, with larger values found in the presence of higher wind speeds and when the radar antenna is looking upwind. This is in line with theoretical models (though models tend to overpredict the range of variability at a given incidence angle) and with observations at lower frequencies. The LDRs are steadily increasing from values certainly lower than −30 dB, at vertical incidence, to the values of about −10 dB, at the incidence angles of about 60°–70°, with a good matching between observations and theoretical predictions. On the other hand, land surface backscattering properties are not characterized by a strong angular dependence: $\sigma _{0}$ and LDR values typically range between −20 and 0 dB and between −15 and −5 dB, respectively. This paper is relevant for spaceborne concepts of W-band radars, which envisage moderate incidence angles to achieve a broad swath needed for global coverage.

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.

This paper describes a new rainfall rate retrieval algorithm, developed within the EUMETSAT H SAF program, based on the Passive microwave Neural network Precipitation Retrieval approach (PNPR v3), designed to work with the conically scanning Global Precipitation Measurement (GPM) Microwave Imager (GMI). A new rain/no-rain classification scheme, also based on the NN approach, which provides different rainfall masks for different minimum thresholds and degree of reliability, is also described. The algorithm is trained on an extremely large observational database, built from GPM global observations between 2014 and 2016, where the NASA 2B-CMB (V04) rainfall rate product is used as reference. In order to assess the performance of PNPR v3 over the globe, an independent part of the observational database is used in a verification study. The good results found over all surface types (CC > 0.90, ME < −0.22 mm h−1, RMSE < 2.75 mm h−1 and FSE% < 100% for rainfall rates lower than 1 mm h−1 and around 30–50% for moderate to high rainfall rates), demonstrate the good outcome of the input selection procedure, as well as of the training and design phase of the neural network. For further verification, two case studies over Italy are also analysed and a good consistency of PNPR v3 retrievals with simultaneous ground radar observations and with the GMI GPROF V05 estimates is found. PNPR v3 is a global rainfall retrieval algorithm, able to optimally exploit the GMI multi-channel response to different surface types and precipitation structures, that provide global rainfall retrieval in a computationally very efficient way, making the product suitable for near-real time operational applications.

Numerous fields of atmospheric and hydrological sciences require the parametric form of the raindrop size distribution (DSD) to estimate the rainfall rate from radar observables as well as in cloud resolving and weather forecasting models. This study aims to investigate the capability of the Johnson SB distribution (JSB) in estimating rain integral parameters. Specifically, rainfall rate (R), reflectivity factor (Z) and mean mass diameter (Dmass) estimated by JSB are compared with those estimated by a three-parameter Gamma distribution, widely used by radar meteorologists and atmospheric physicists to model natural DSD. A large dataset consisting of more than 155,000 one-minute DSD, from six field campaigns of Ground Validation (GV) program of NASA/JAXA Global Precipitation Measurement mission (GPM), is used to test the performance of both JSB and Gamma distribution. The available datasets cover a wide range of rain regimes because of the field campaigns were carried out in different seasons and locations. Correlation coefficient, bias, root mean square error (RMSE) and fractional standard error (FSE) between estimated and measured integral parameters are calculated to compare the performances of the two distributions. The capability of JSB in estimating the integral parameters, especially R and Z, resulted very close to that of Gamma distribution. In particular, for light precipitation, JSB is superior to Gamma distribution in estimating R with FSE of 11% with respect to values ranging between 25% and 37% about for Gamma. Comparison of the estimated and measured DSDs shows that the JSB distribution reproduces the natural DSD quite accurately.