Grant W. Petty

The Evolution of the Goddard Profiling Algorithm (GPROF) for Rainfall Estimation from Passive Microwave Sensors

Abstract This paper describes the latest improvements applied to the Goddard profiling algorithm (GPROF), particularly as they apply to the Tropical Rainfall Measuring Mission (TRMM). Most of these improvements, however, are conceptual in nature and apply equally to other passive microwave sensors. The improvements were motivated by a notable overestimation of precipitation in the intertropical convergence zone. This problem was traced back to the algorithms poor separation between convective and stratiform precipitation coupled with a poor separation between stratiform and transition regions in the a priori cloud model database. In addition to now using an improved convective–stratiform classification scheme, the new algorithm also makes use of emission and scattering indices instead of individual brightness temperatures. Brightness temperature indices have the advantage of being monotonic functions of rainfall. This, in turn, has allowed the algorithm to better define the uncertainties needed by the sc...

Intercomparison of global precipitation products : The third Precipitation Intercomparison Project (PIP-3)

Abstract A set of global, monthly rainfall products has been intercompared to understand the quality and utility of the estimates. The products include 25 observational (satellite based), four model, and two climatological products. The results of the intercomparison indicate a very large range (factor of 2 or 3) of values when all products are considered. The range of values is reduced considerably when the set of observational products is limited to those considered quasi–standard. The model products do significantly poorer in the Tropics, but are competitive with satellite–based fields in midlatitudes over land. Over ocean, products are compared to frequency of precipitation from ship observations. The evaluation of the observational products points to merged data products (including rain gauge information) as providing the overall best results.

Physical retrievals of over-ocean rain rate from multichannel microwave imagery. Part I: Theoretical characteristics of normalized polarization and scattering indices

SummaryMicrowave rain rate retrieval algorithms have most often been formulated in terms of the raw brightness temperatures observed by one or more channels of a satellite radiometer. Taken individually, single-channel brightness temperatures generally represent a near-arbitrary combination of positive contributions due to liquid water emission and negative contributions due to scattering by ice and/or visibility of the radiometrically cold ocean surface. Unfortunately, for a given rain rate, emission by liquid water below the freezing level and scattering by ice particles above the freezing level are rather loosely coupled in both a physical and statistical sense. Furthermore, microwave brightness temperatures may vary significantly (∼30–70 K) in response to geophysical parameters other than liquid water and precipitation. Because of these complications, physical algorithms which attempt to directly invert observed brightness temperatures have typically relied on the iterative adjustment of detailed microphysical profiles or cloud models, guided by explicit forward microwave radiative transfer calculations.In support of an effort to develop a significantly simpler and more efficient inversion-type rain rate algorithm, the physical information content of two linear transformations of single-frequency, dual-polarization brightness temperatures is studied: thenormalized polarization difference P of Petty and Katsaros (1990, 1992), which is intended as a measure of footprint-averaged rain cloud transmittance for a given frequency; and ascattering index S (similar to the polarization corrected temperature of Spencer et al., 1989) which is sensitive almost exclusively to ice. A reverse Monte Carlo radiative transfer model is used to elucidate the qualitative response of these physically distinct single-frequency indices to idealized 3-dimensional rain clouds and to demonstrate their advantages over raw brightness temperatures both as stand-alone indices of precipitation activity and as primary variables in physical, multichannel rain rate retrieval schemes.As a byproduct of the present analysis, it is shown that conventional plane-parallel analyses of the well-known footprint-filling problem for emission-based algorithms may in some cases give seriously misleading results.

Results of WetNet PIP-2 Project

The second WetNet Precipitation Intercomparison Project (PIP-2) evaluates the performance of 20 satellite precipitation retrieval algorithms, implemented for application with Special Sensor Microwave/Imager (SSM/I) passive microwave (PMW) measurements and run for a set of rainfall case studies at full resolution‐instantaneous space‐timescales. The cases are drawn from over the globe during all seasons, for a period of 7 yr, over a 608N‐ 178S latitude range. Ground-based data were used for the intercomparisons, principally based on radar measurements but also including rain gauge measurements. The goals of PIP-2 are 1) to improve performance and accuracy of different SSM/I algorithms at full resolution‐instantaneous scales by seeking a better understanding of the relationship between microphysical signatures in the PMW measurements and physical laws employed in the algorithms; 2) to evaluate the pros and cons of individual algorithms and their subsystems in order to seek optimal ‘‘front-end’’ combined algorithms; and 3) to demonstrate that PMW algorithms generate acceptable instantaneous rain estimates. It is found that the bias uncertainty of many current PMW algorithms is on the order of 630%. This level is below that of the radar and rain gauge data specially collected for the study, so that it is not possible to objectively select a best algorithm based on the ground data validation approach. By decomposing the intercomparisons into effects due to rain detection (screening) and effects due to brightness temperature‐rain rate conversion, differences among the algorithms are partitioned by rain area and rain intensity. For ocean, the screening differences mainly affect the light rain rates, which do not contribute significantly to area-averaged rain rates. The major sources of differences in mean rain rates between individual algorithms stem from differences in how intense rain rates are calculated and the maximum rain rate allowed by a given algorithm. The general method of solution is not necessarily the determining factor in creating systematic rain-rate differences among groups of algorithms, as we find that the severity of the screen is the dominant factor in producing systematic group differences among land algorithms, while the input channel selection is the dominant factor in producing systematic group differences among ocean algorithms. The significance of these issues are examined through what is called ‘‘fan map’’ analysis. The paper concludes with a discussion on the role of intercomparison projects in seeking improvements to algorithms, and a suggestion on why moving beyond the ‘‘ground truth’’ validation approach by use of a calibration-quality forward model would be a step forward in seeking objective evaluation of individual algorithm performance and optimal algorithm design.

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry. Part I: Improved Method and Uncertainties

Abstract A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and nonconvective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subse...

Physical retrievals of over-ocean rain rate from multichannel microwave imagery. Part II: Algorithm implementation

SummaryA new physical inversion-based algorithm for retrieving rain rate over the ocean with the Special Sensor Microwave Imager (SSM/I) is described. In a departure from other rain rate retrieval algorithms, the satellite observables inverted in the present algorithm are not the raw brightness temperatures but rather normalized polarizations for 19.35, 37.0, and 85.5 GHz, plus an 85.5 GHz scattering index which is sensitive primarily to ice particles aloft. The normalized polarizations are interpreted as footprint-averages of theoretically derived analytic functions of the column optical depth associated primarily with liquid water. The effective vertical depth of the rain layer is specified as a function of the SSM/I estimated column water vapor.The retrieval algorithm performs an iterative search for a high resolution (12.5 km) rain field which is simultaneously consistent with the 19.35 and 37.0 GHz normalized polarizations. The first-guess rain rate field is supplied by the 85.5 GHz scattering index. At gridpoints for which the rain column optical depth exceeds the dynamic range of the attenuation-based indices, the first-guess field is left essentially unmodified; elsewhere, the required consistency with the 19 and 37 GHz indices usually results in significant modification of the scattering-based rain rate estimates.The algorithm as described here is a prototype implementation which was developed with reference only to idealized theoretical models; empirical improvements to the numerical scheme and the model coefficients will be made in the future as results from the first Precipitation [algorithm] Intercomparison Project 1 (PIP-1) and the second phase of the GPCP (Global Precipitation Climatology Project) algorithm Intercomparison Project (AIP/2) are analyzed, as well as data from individual validation efforts. Although the algorithm is physically based and uses all SSM/I channels, it is computationally much less demanding than cloud/radiative transfer model-based inversion algorithms published else-where.

The Sensitivity of Microwave Remote Sensing Observations of Precipitation to Ice Particle Size Distributions

Abstract This study investigates the effect of variable size distribution and density of precipitation ice particles on microwave brightness temperatures. For this purpose, a set of self-consistent relationships among rain rate, size parameters of an exponential drop size distribution, and the radar reflectivity–rain rate relations for frozen and liquid precipitation was derived. Further, a scaling factor was introduced that is the ratio between the average melted diameter of the frozen and liquid precipitation and allows the specification of different sizes of the frozen particles. For given radar observations, this method allows size distributions of frozen and liquid precipitation to be derived, which are then used as input for a radiative transfer model. These relationships were used to perform Mie calculations for different precipitation rates and different types of hydrometeors (snow, graupel, and hail) and to investigate the dependence of their respective optical properties on rain rate as well as ...

Microwave Backscatter and Extinction by Soft Ice Spheres and Complex Snow Aggregates

Abstract Coupled-dipole approximation (CDA) calculations of microwave extinction and radar backscatter are presented for nonhomogeneous (soft) ice spheres and for quasi-realistic aggregates of elementary ice crystal forms, including both simple needles and real dendrites. Frequencies considered include selections from the Dual-Frequency Precipitation Radar (DPR; 13.4 and 35.6 GHz) and the Global Precipitation Measurement (GPM) Microwave Imager (GMI; 18.7, 36.5, and 89.0 GHz), both slated for orbit on the GPM mission. The computational method is first validated against Mie theory using dipole structures representing solid ice spheres as well as stochastically generated “soft” ice spheres of variable ice–air ratio. Neither the traditional Bruggeman nor Maxwell Garnett dielectric mixing formula is found to correctly predict the full range of CDA results for soft spheres. However, an excellent fit is found using the generalized mixing rule of Sihvola with ν = 0.85. The suitability of the soft-sphere approxima...

The status of satellite-based rainfall estimation over land

During the past decade, substantial progress has been made in the area of rainfall estimation from space. The purpose of this paper is to provide land surface scientists with a broad overview of that progress, to highlight what is known about the performance of current rainfall estimation techniques, and to draw attention to some of the practical developments and instrument improvements anticipated during the remainder of the present decade.

Frequencies and Characteristics of Global Oceanic Precipitation from Shipboard Present-Weather Reports

Abstract Ship reports of present weather obtained from the Comprehensive Ocean-Atmosphere Data Set are analyzed for the period 1958–91 in order to elucidate regional and seasonal variations in the climatological frequency, phase. intensity, and character of oceanic precipitation. Specific findings of note include the following: 1) The frequency of thunderstorm reports, relative to all precipitation reports, is a strong function of location, with thunderstorm activity being favored within 1000–3000 km of major tropical and subtropical landmasses, while being quite rare at other locations, even within the intertropical convergence zone. 2) The latitudinal frequency of precipitation over the southern oceans increases steadily toward the Antarctic continent and shows relatively little seasonal variation. The frequency of convective activity, however, shows considerable seasonal variability, with sharp winter maxima occurring near 381 latitude in both hemispheres. 3) Drizzle is the preferred form of precipitat...