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.

A comprehensive and accurate global water vapor dataset is critical to the adequate understanding of water vapor's role in the earth's climate system. To begin to satisfy this need, the authors have produced a blended dataset made up of global, 5-yr (1988-92), l°x 1° spatial resolution, atmospheric water vapor (WV) and liquid water path products. These new products consist of both the daily total column-integrated composites and a multilayered WV product at three layers (1000-700, 700-500, 500-300 mb). The analyses combine WV retrievals from the Television and Infrared Operational Satellite (TIROS) Operational Vertical Sounder (TOVS), the Special Sensor Microwave/Imager, and radiosonde observations. The global, vertical-layered water vapor dataset was developed by slicing the blended total column water vapor using layer information from TOVS and radiosonde. Also produced was a companion, over oceans only, liquid water path dataset. Satellite observations of liquid water path are growing in importance since many of the global climate models are now either incorporating or contain liquid water as an explicit variable. The complete dataset (all three products) has been named NVAP, an acronym for National Aeronautics and Space Administration Water Vapor Project. This paper provides examples of the new dataset as well as scientific analysis of the observed annual cycle and the interannual variability of water vapor at global, hemispheric, and regional scales. A distinct global annual cycle is shown to be dominated by the Northern Hemisphere observations. Planetary-scale variations are found to relate well to recent independent estimates of tropospheric temperature variations. Maps of regional interannual variability in the 5-yr period show the effect of the 1992 ENSO and other features.

The combination of active and passive microwave sensors on board the Tropical Rainfall Measuring Mission (TRMM) satellite have been used to construct observationally constrained databases of precipitation profiles for use in passive microwave rainfall retrieval algorithms over oceans. The method uses a very conservative approach that begins with the operational TRMM precipitation radar algorithm and adjusts its solution only as necessary to simultaneously match the radiometer observations. Where the TRMM precipitation radar (PR) indicates no rain, an optimal estimation procedure using TRMM Microwave Imager (TMI) radiances is used to retrieve nonraining parameters. The optimal estimation methodology ensures that the geophysical parameters are fully consistent with the observed radiances. Within raining fields of view, cloud-resolving model outputs are matched to the liquid and frozen hydrometeor profiles retrieved by the TRMM PR. The profiles constructed in this manner are subsequently used to compute brightness temperatures that are immediately compared to coincident observations from TMI. Adjustments are made to the rainwater and ice concentrations derived by PR in order to achieve agreement at 19 and 85 GHz, vertically polarized brightness temperatures at monthly time scales. The database is generated only in the central 11 pixels of the PR radar scan, and the rain adjustment is performed independently for distinct sea surface temperature (SST) and total precipitable water (TPW) values. Overall, the procedure increases PR rainfall by 4.2%, but the adjustment is not uniform across all SST and TPW regimes. Rainfall differences range from a minimum of −57% for SST of 293 K and TPW of 13 mm to a maximum of +53% for SST of 293 K and TPW of 45 mm. These biases are generally reproduced by a TMI retrieval algorithm that uses the observationally generated database. The algorithm increases rainfall by 5.0% over the PR solution with a minimum of −99% for SST of 293 K and TPW of 14 mm to a maximum of +11.8% for an SST of 294 K and TPW of 50 mm. Some differences are expected because of the algorithm mechanics.

The estimation of precipitation across the globe from satellite sensors provides a key resource in the observation and understanding of our climate system. Estimates from all pertinent satellite observations are critical in providing the necessary temporal sampling. However, consistency in these estimates from instruments with different frequencies and resolutions is critical. This paper details the physically based retrieval scheme to estimate precipitation from cross-track (XT) passive microwave (PM) sensors on board the constellation satellites of the Global Precipitation Measurement (GPM) mission. Here the Goddard profiling algorithm (GPROF), a physically based Bayesian scheme developed for conically scanning (CS) sensors, is adapted for use with XT PM sensors. The present XT GPROF scheme utilizes a model-generated database to overcome issues encountered with an observational database as used by the CS scheme. The model database ensures greater consistency across meteorological regimes and surface types by providing a more comprehensive set of precipitation profiles. The database is corrected for bias against the CS database to ensure consistency in the final product. Statistical comparisons over western Europe and the United States show that the XT GPROF estimates are comparable with those from the CS scheme. Indeed, the XT estimates have higher correlations against surface radar data, while maintaining similar root-mean-square errors. Latitudinal profiles of precipitation show the XT estimates are generally comparable with the CS estimates, although in the southern midlatitudes the peak precipitation is shifted equatorward while over the Arctic large differences are seen between the XT and the CS retrievals.

An intercomparison of radiation codes used in retrieving upper-tropospheric humidity (UTH) from observations in the v2 (6.3 /xm) water vapor absorption band was performed. This intercomparison is one part of a coordinated effort within the Global Energy and Water Cycle Experiment Water Vapor Project to assess our ability to monitor the distribution and variations of upper-tropospheric moisture from spaceborne sensors. A total of 23 different codes, ranging from detailed line-by-line (LBL) models, to coarser-resolution narrowband (NB) models, to highly parameterized single-band (SB) models participated in the study. Forward calculations were performed using a carefully selected set of temperature and moisture profiles chosen to be representative of a wide range of atmospheric conditions. The LBL model calculations exhibited the greatest consistency with each other, typically agreeing to within 0.5 K in terms of the equivalent blackbody brightness temperature (T b ). The majority of NB and SB models agreed to within ±1 K of the LBL models, although a few older models exhibited systematic T b biases in excess of 2 K. A discussion of the discrepancies between various models, their association with differences in model physics (e.g., continuum absorption), and their implications for UTH retrieval and radiance assimilation is presented.