The Goddard profiling algorithm has evolved from a pseudoparametric algorithm used in the current TRMM operational product (GPROF 2010) to a fully parametric approach used operationally in the GPM era (GPROF 2014). The fully parametric approach uses a Bayesian inversion for all surface types. The algorithm thus abandons rainfall screening procedures and instead uses the full brightness temperature vector to obtain the most likely precipitation state. This paper offers a complete description of the GPROF 2010 and GPROF 2014 algorithms and assesses the sensitivity of the algorithm to assumptions related to channel uncertainty as well as ancillary data. Uncertainties in precipitation are generally less than 1%–2% for realistic assumptions in channel uncertainties. Consistency among different radiometers is extremely good over oceans. Consistency over land is also good if the diurnal cycle is accounted for by sampling GMI product only at the time of day that different sensors operate. While accounting for only a modest amount of the total precipitation, snow-covered surfaces exhibit differences of up to 25% between sensors traceable to the availability of high-frequency (166 and 183 GHz) channels. In general, comparisons against early versions of GPM’s Ku-band radar precipitation estimates are fairly consistent but absolute differences will be more carefully evaluated once GPROF 2014 is upgraded to use the full GPM-combined radar–radiometer product for its a priori database. The combined algorithm represents a physically constructed database that is consistent with both the GPM radars and the GMI observations, and thus it is the ideal basis for a Bayesian approach that can be extended to an arbitrary passive microwave sensor.

Rainfall retrieval algorithms often assume a gamma-shaped raindrop size distribution (DSD) with three mathematical parameters N w , D m , and m. If only two independent measurements are available, as with the dualfrequency precipitation radar on the Global Precipitation Measurement (GPM) mission core satellite, then retrieval algorithms are underconstrained and require assumptions about DSD parameters. To reduce the number of free parameters, algorithms can assume that m is either a constant or a function of D m . Previous studies have suggested m-L constraints [where L 5 (4 1 m)/D m ], but controversies exist over whether m-L constraints result from physical processes or mathematical artifacts due to high correlations between gamma DSD parameters. This study avoids mathematical artifacts by developing joint probability distribution functions (joint PDFs) of statistically independent DSD attributes derived from the raindrop mass spectrum. These joint PDFs are then mapped into gamma-shaped DSD parameter joint PDFs that can be used in probabilistic rainfall retrieval algorithms as proposed for the GPM satellite program. Surface disdrometer data show a high correlation coefficient between the mass spectrum mean diameter D m and mass spectrum standard deviation s m . To remove correlations between DSD attributes, a normalized mass spectrum standard deviation s 0 m is constructed to be statistically independent of D m , with s 0 m representing the most likely value and std(s 0 m ) representing its dispersion. Joint PDFs of D m and m are created from D m and s 0 m . A simple algorithm shows that rain-rate estimates had smaller biases when assuming the DSD breadth of s 0 m than when assuming a constant m.

In this paper, the operational Global Precipitation Measurement (GPM) mission combined radar–radiometer algorithm is thoroughly described. The operational combined algorithm is designed to reduce uncertainties in GPM Core Observatory precipitation estimates by effectively integrating complementary information from the GPM Dual-Frequency Precipitation Radar (DPR) and the GPM Microwave Imager (GMI) into an optimal, physically consistent precipitation product. Although similar in many respects to previously developed combined algorithms, the GPM combined algorithm has several unique features that are specifically designed to meet the GPM objectives of deriving, based on GPM Core Observatory information, accurate and physically consistent precipitation estimates from multiple spaceborne instruments, and ancillary environmental data from reanalyses. The algorithm features an optimal estimation framework based on a statistical formulation of the Gauss–Newton method, a parameterization for the nonuniform distribution of precipitation within the radar fields of view, a methodology to detect and account for multiple scattering in Ka-band DPR observations, and a statistical deconvolution technique that allows for an efficient sequential incorporation of radiometer information into DPR precipitation retrievals.

Retrievals of falling snow from space-based observations represent key inputs for understanding and linking Earth’s atmospheric, hydrological, and energy cycles. This work quantifies and investigates causes of differences among the first stable falling snow retrieval products from the Global Precipitation Measurement (GPM) Core Observatory satellite and CloudSat’s Cloud Profiling Radar (CPR) falling snow product. An important part of this analysis details the challenges associated with comparing the various GPM and CloudSat snow estimates arising from different snow–rain classification methods, orbits, resolutions, sampling, instrument specifications, and algorithm assumptions. After equalizing snow–rain classification methodologies and limiting latitudinal extent, CPR observes nearly 10 (3) times the occurrence (accumulation) of falling snow as GPM’s Dual-Frequency Precipitation Radar (DPR). The occurrence disparity is substantially reduced if CloudSat pixels are averaged to simulate DPR radar pixels and CPR observations are truncated below the 8-dBZ reflectivity threshold. However, even though the truncated CPR- and DPR-based data have similar falling snow occurrences, average snowfall rate from the truncated CPR record remains significantly higher (43%) than the DPR, indicating that retrieval assumptions (microphysics and snow scattering properties) are quite different. Diagnostic reflectivity (Z)–snow rate (S) relationships were therefore developed at Ku and W band using the same snow scattering properties and particle size distributions in a final effort to minimize algorithm differences. CPR–DPR snowfall amount differences were reduced to ~16% after adopting this diagnostic Z–S approach.

No abstract

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.

The estimation of the drop size distribution parameter [median volume diameter (D 0 )] and rain rate (R) from the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) as well as from combined PR-TRMM Microwave Imager (TMI) algorithms are considered in this study for two TRMM satellite overpasses near the Kwajalein Atoll. An operational dual-polarized S-band radar (KPOL) located in Kwajalein is central as the only TRMM ground validation site for measurement of precipitation over the open ocean. The accuracy of the TRMM PR in retrieving D 0 and R is better for precipitation over the ocean based on a more stable surface reference technique for estimating the path-integrated attenuation. Also, combined PR-TMI methods are more accurate over the open ocean because of better knowledge of the surface microwave emissivity. Using Z h (horizontal polarized radar reflectivity) and Z dr (differential reflectivity) data for the two TRMM overpass events over Kwajalein, D 0 and R from KPOL are retrieved. Herein, the main objective is to see if the D 0 retrieved from either PR or the combined PR-TMI algorithms are in agreement with KPOL-derived values. Also, the variation of D 0 versus R is compared for convective rain pixels from KPOL, PR, and PR-TMI. It is shown that the PR-TMI optimal estimation scheme does indeed adjust the D 0 in the ''correct'' direction, on average, from the a priori state if the KPOL data are considered to be the ground truth. This correct adjustment may be considered as evidence of the value added by the TMI brightness temperatures in the combined PR-TMI variational scheme, at least for the two overpass events considered herein.

Precipitation, including rain and snow, is a critical part of the Earth's energy and hydrology cycles. Precipitation impacts latent heating profiles locally while global circulation patterns distribute precipitation and energy from the equator to the poles. For the hydrological cycle, falling snow is a primary contributor in northern latitudes during the winter seasons. Falling snow is the source of snow pack accumulations that provide fresh water resources for many communities in the world. Furthermore, falling snow impacts society by causing transportation disruptions during severe snow events. In order to collect information on the complete global precipitation cycle, both liquid and frozen precipitation must be collected.The challenges of estimating falling snow from space still exist though progress is being made. These challenges include weak falling snow signatures with respect to background (surface, water vapor) signatures for passive sensors over land surfaces, unknowns about the spherical and non-spherical shapes of the snowflakes, their particle size distributions (PSDs) and how the assumptions about the unknowns impact observed brightness temperatures or radar reflectivities, differences in near surface snowfall and total column snow amounts, and limited ground truth to validate against. While these challenges remain, knowledge of their impact on expected retrieval results is an important key for understanding falling snow retrieval estimations.Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types. For example, can a lake effect snow system with low (-2.5 km) cloud tops having an ice water content (IWC) at the surface of 0.25 g m· 3 and dendrite snowflakes be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results.Here, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases since simulations provide all the information to determine the measurements from space and the ground truth. Sensitivity analyses were performed to better ascertain the relationships between multifrequency microwave and millimeter-wave sensor observations and the falling snow/underlying field of view. In addition, thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types were studied. Results will be presented for active radar at Ku, Ka, and W-band and for passive radiometer channels from 10 to 183 GHz.https://ntrs.nasa.gov/search