The global precipitation measurement mission (GPM) has been in operation for seven years and continues to provide a vast quantity of global precipitation data at finer temporospatial resolutions with improved accuracy and coverage. GPM’s signature algorithm, the integrated multisatellite retrievals for GPM (IMERG) is a next-generation of precipitation product expected for wide variety of research and operational applications. This study evaluates the latest version (V06B) of IMERG and its predecessor, the tropical rainfall measuring mission (TRMM) multisatellite precipitation (TMPA) 3B42 (V7) using ground-based and gauge-corrected multiradar multisensor system (MRMS) precipitation products over the conterminous United States (CONUS). The spatial distributions of all products are analyzed. The error characteristics are further examined for 3B42 and IMERG in winter and summer by an error decomposition approach, which partitions total bias into hit bias, biases due to missed precipitation and false precipitation. The volumetric and categorical statistical metrics are used to quantitatively evaluate the performance of the two satellite-based products. All products show a similar precipitation climatology with some regional differences. The two satellite-based products perform better in the eastern CONUS than in the mountainous Western CONUS. The evaluation demonstrates the clear improvement in IMERG precipitation product in comparison with its predecessor 3B42, especially in reducing missed precipitation in winter and summer, and hit bias in winter, resulting in better performance in capturing lighter and heavier precipitation.

Remote-sensing observations are needed to estimate the regional and global impacts of snow. However, to retrieve accurate estimates of snow mass and rate, these observations require augmentation through additional information and assumptions about hydrometeor properties. The Precipitation Imaging Package (PIP) provides information about precipitation characteristics and can be utilized to improve estimates of snowfall rate and accumulation. Here, the goal is to demonstrate the quality and utility of two higher-order PIP-derived products: liquid water equivalent snow rate and an approximation of volume-weighted density called equivalent density. Accuracy of the PIP snow rate and equivalent density is obtained through intercomparison with established retrieval methods and through evaluation with colocated ground-based observations. The results confirm the ability of the PIP-derived products to quantify properties of snow rate and equivalent density, and demonstrate that the PIP produces physically realistic snow characteristics. When compared to the National Weather Service (NWS) snow field measurements of six-hourly accumulation, the PIP-derived accumulations were biased only +2.48% higher. Additionally, this work illustrates fundamentally different microphysical and bulk features of low and high snow-to-liquid ratio events, through assessment of observed particle size distributions, retrieved mass coefficients, and bulk properties. Importantly, this research establishes the role that PIP observations and higher-order products can serve for constraining microphysical assumptions in ground-based and spaceborne remotely sensed snowfall retrievals.

A unique capability of the Global Precipitation Measurement (GPM) mission is its ability to better estimate the raindrop size distribution (DSD) on a global scale. To validate the GPM DSD retrievals, a network of more than 100 ground-based polarimetric radars from across the globe are utilized within the broader context of the GPM Validation Network (VN) processing architecture. The GPM VN ensures quality controlled dual-polarimetric radar moments for use in providing reference estimates of the DSD. The VN DSD estimates are carefully geometrically matched with the GPM core satellite measurements for evaluation of the GPM algorithms. We use the GPM VN to compare the DSD retrievals from the GPM’s Dual-frequency Precipitation Radar (DPR) and combined DPR–GPM Microwave Imager (GMI) Level-2 algorithms. Results suggested that the Version 06A GPM core satellite algorithms provide estimates of the mass-weighted mean diameter (D_m) that are biased 0.2 mm too large when considered across all precipitation types. In convective precipitation, the algorithms tend to overestimate D_m by 0.5–0.6 mm, leading the DPR algorithm to underestimate the normalized DSD intercept parameter (N_w) by a factor of two, and introduce a significant bias to the DPR retrievals of rainfall rate for DSDs with large D_m. The GPM Combined algorithm performs better than the DPR algorithm in convection but provides a severely limited range of N_w estimates, highlighting the need to broaden its a priori database in convective precipitation.

Spaceborne microwave sensors provide critical rain information used in several global multisatellite rain products, which in turn are used for a variety of important studies, including landslide forecasting, flash flood warning, data assimilation, climate studies, and validation of model forecasts of precipitation. This study employs 4 yr (2003–06) of satellite data to assess the relative performance and skill of the Special Sensor Microwave Imager [SSM/I (F13, F14, and F15], Advanced Microwave Sounding Unit [AMSU-B (N15, N16, and N17)], Advanced Microwave Scanning Radiometer for Earth Observing System [AMSR-E (Aqua)], and the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) in estimating surface rainfall based on direct instantaneous comparisons with ground-based rain estimates from the TRMM Ground Validation (GV) sites at Kwajalein, Republic of the Marshall Islands (KWAJ), and Melbourne, Florida (MELB). The relative performance of each of these satellite estimates is examined via comparisons with space- and time-coincident GV radar-based rain-rate estimates. Because underlying surface terrain is known to affect the relative performance of the satellite algorithms, the data for MELB were further stratified into ocean, land, and coast categories using a 0.25° terrain mask. Of all the satellite estimates compared in this study, TMI and AMSR-E exhibited considerably higher correlations and skills in estimating–observing surface precipitation. While SSM/I and AMSU-B exhibited lower correlations and skills for each of the different terrain categories, the SSM/I absolute biases trended slightly lower than AMSR-E over ocean, where the observations from both emission and scattering channels were used in the retrievals. AMSU-B exhibited the least skill relative to GV in all of the relevant statistical categories, and an anomalous spike was observed in the probability distribution functions near 1.0 mm h−1. This statistical artifact appears to be related to attempts by algorithm developers to include some lighter rain rates, not easily detectable by its scatter-only frequencies. AMSU-B, however, agreed well with GV when the matching data were analyzed on monthly scales. These results signal to developers of global rainfall products, such as the TRMM Multisatellite Precipitation Analysis (TMPA) and the Climate Data Center’s Morphing (CMORPH) technique, that care must be taken when incorporating data from these input satellite estimates to provide the highest-quality estimates in their products.

Surface precipitation phase is a fundamental meteorological property with immense importance. Accurate classification of phase from satellite remotely sensed observations is difficult. This study demonstrates the ability of the Precipitation Imaging Package (PIP), a ground-based, in situ precipitation imager, to distinguish precipitation phase. The PIP precipitation phase identification capabilities are compared to observer records from the National Weather Service (NWS) office in Marquette, Michigan, as well as co-located observations from profiling and scanning radars, disdrometer data, and surface meteorological measurements. Examined are 13 events with at least one precipitation phase transition. The PIP-determined onsets and endings of the respective precipitation phase periods agree to within 15 min of NWS observer records for the vast majority of the events. Additionally, the PIP and NWS liquid water equivalent accumulations for 12 of the 13 events were within 10%. Co-located observations from scanning and profiling radars, as well as reanalysis-derived synoptic and thermodynamic conditions, support the accuracy of the precipitation phases identified by the PIP. PIP observations for the phase transition events are compared to output from a parameterization based on wet bulb and near-surface lapse rates to produce a probability of solid precipitation. The PIP phase identification and the parameterization output are consistent. This work highlights the ability of the PIP to properly characterize hydrometeor phase and provide dependable precipitation accumulations under complicated mixed-phase and rain and snow (or vice versa) transition events.

Hurricane rainbands are very efficient rain producers, but details on drop size distributions are still lacking. This study focuses on the rainbands of hurricane Dorian as they traversed the densely instrumented NASA precipitation-research facility at Wallops Island, VA, over a period of 8 h. Drop size distribution (DSD) was measured using a high-resolution meteorological particle spectrometer (MPS) and 2D video disdrometer, both located inside a double-fence wind shield. The shape of the DSD was examined using double-moment normalization, and compared with similar shapes from semiarid and subtropical sites. Dorian rainbands had a superexponential shape at small normalized diameter values similar to those of the other sites. NASA’s S-band polarimetric radar performed range height-indicator (RHI) scans over the disdrometer site, showing some remarkable signatures in the melting layer (bright-band reflectivity peaks of 55 dBZ, a dip in the copolar correlation to 0.85 indicative of 12–15 mm wet snow, and a staggering reflectivity gradient above the 0 °C level of −10 dB/km, indicative of heavy aggregation). In the rain layer at heights < 2.5 km, polarimetric signatures indicated drop break-up as the dominant process, but drops as large as 5 mm were detected during the intense bright-band period.