ALBERT

Description: The paper examines how well an Eddington approximation can reproduce brightness temperatures obtained from a more complete, N-stream discrete ordinate solution in the microwave regime. Radiation propagation through a plane parallel medium is considered. Although model discrepancies are complicated functions of the cloud constituents, the differences between an eight-stream discrete ordinate solution and an analytical Eddington solution were found to be generally small, ranging from 0 to 6 K when only one uniform layer of hydrometeors was considered. When realistic multilayered cloud hydrometeor profiles were used, the differences between these two models never exceeded 3 K over the entire range of microwave frequencies considered (6.6-183 GHz). The models agreed to within 0.2 K in the absence of scattering constituents.

Description: This study presents the first quantitative retrievals of vertical profiles of precipitation derived from multispectral passive microwave radiometry. Measurements of microwave brightness temperature (Tb) obtained by a NASA high-altitude research aircraft are related to profiles of rainfall rate through a multichannel piecewise-linear statistical regression procedure. Statistics for Tb are obtained from a set of cloud radiative models representing a wide variety of convective, stratiform, and anvil structures. The retrieval scheme itself determines which cloud model best fits the observed meteorological conditions. Retrieved rainfall rate profiles are converted to equivalent radar reflectivity for comparison with observed reflectivities from a ground-based research radar. Results for two case studies, a stratiform rain situation and an intense convective thunderstorm, show that the radiometrically derived profiles capture the major features of the observed vertical structure of hydrometer density.

Description: A multichannel statistical approach is used to retrieve rainfall rates from the brightness temperature T(B) observed by passive microwave radiometers flown on a high-altitude NASA aircraft. T(B) statistics are based upon data generated by a cloud radiative model. This model simulates variabilities in the underlying geophysical parameters of interest, and computes their associated T(B) in each of the available channels. By further imposing the requirement that the observed T(B) agree with the T(B) values corresponding to the retrieved parameters through the cloud radiative transfer model, the results can be made to agree quite well with coincident radar-derived rainfall rates. Some information regarding the cloud vertical structure is also obtained by such an added requirement. The applicability of this technique to satellite retrievals is also investigated. Data which might be observed by satellite-borne radiometers, including the effects of nonuniformly filled footprints, are simulated by the cloud radiative model for this purpose.

Description: Algorithms that are currently being developed for the Tropical Rainfall Measuring Mission for the retrieval of the surface rainfall and of the rain-system vertical structure from passive microwave measurements are described. The utility of the procedure is illustrated using data from the Special Sensor Microwave/Imager. The results, which are still of preliminary nature, show the advantages of incorporating the rain-system vertical structure into the precipitation retrievals.

Description: This paper describes a multichannel physical approach for retrieving rainfall and vertical structure information from satellite-based passive microwave observations. The algorithm makes use of statistical inversion techniques based upon theoretically calculated relations between rainfall rates and brightness temperatures. Potential errors introduced into the theoretical calculations by the unknown vertical distribution of hydrometeors are overcome by explicity accounting for diverse hydrometeor profiles. This is accomplished by allowing for a number of different vertical distributions in the theoretical brightness temperature calculations and requiring consistency between the observed and calculated brightness temperatures. This paper will focus primarily on the theoretical aspects of the retrieval algorithm, which includes a procedure used to account for inhomogeneities of the rainfall within the satellite field of view as well as a detailed description of the algorithm as it is applied over both ocean and land surfaces. The residual error between observed and calculated brightness temperatures is found to be an important quantity in assessing the uniqueness of the solution. It is further found that the residual error is a meaningful quantity that can be used to derive expected accuracies from this retrieval technique. Examples comparing the retrieved results as well as the detailed analysis of the algorithm performance under various circumstances are the subject of a companion paper.

Description: Because passive microwave instruments are confined to polar-orbiting satellites, rainfall estimates must interpolate across long time periods, during which no measurements are available. In this paper the authors discuss a technique that allows one to partially overcome the sampling limitations by using frequent infrared observations from geosynchronous platforms. To accomplish this, the technique compares all coincident microwave and infrared observations. From each coincident pair, the infrared temperature threshold is selected that corresponds to an area equal to the raining area observed in the microwave image. The mean conditional rainfall rate as determined from the microwave image is then assigned to pixels in the infrared image that are colder than the selected threshold. The calibration is also applied to a fixed threshold of 235 K for comparison with established infrared techniques. Once a calibration is determined, it is applied to all infrared images. Monthly accumulations for both methods are then obtained by summing rainfall from all available infrared images. Two examples are used to evaluate the performance of the technique. The first consists of a one-month period (February 1988) over Darwin, Australia, where good validation data are available from radar and rain gauges. For this case it was found that the technique approximately doubled the rain inferred by the microwave method alone and produced exceptional agreement with the validation data. The second example involved comparisons with atoll rain gauges in the western Pacific for June 1989. Results here are overshadowed by the fact that the hourly infrared estimates from established techniques, by themselves, produced very good correlations with the rain gauges. The calibration technique was not able to improve upon these results.

Description: The second intercomparison project of the Global Precipitation Climatology Project (GPCP) examined the estimation of midlatitude, cool-season precipitation. As part of that effort, the authors report here on the results of two microwave techniques, the Goddard scattering algorithm and the physical retrieval algorithm of Kummerow. Results from the estimation of instantaneous rain rate for five overpasses of the Special Sensor Microwave/Imager (SSM/I) are presented in a case study mode to illustrate both the strong and weak points of each technique. These five cases represent a sampling of the various types of precipitating systems observed. Results for the complete set of 20 swaths chosen by the United Kingdom Meteorological Office (UKMO) are then categorized by scatterplots and statistics of instantaneous radar versus microwave-estimated rain rate, rain no-rain contingency tables, and scatterplots of areal coverage of rainfall. Neither algorithm produced a good statistical correlation with the radar data, yet in general, both did well at determining rainy areas. Two reasons are suggested for the low correlation coefficients between both algorithms and the radar data. Time differences between the SSM/I overpass and the radar observations can occasionally account for some of the differences. The primary reason for the low correlations, however, appears to be the predominance of very light rain in the area of interest during the winter. Both algorithms are in good spatial agreement with the radar when the radar data are restricted to rates above 1 mm/h. When all radar rain rates are included, the radar areal coverage increases by as much as a factor of 10 in some cases. Because the Kummerow algorithm does not handle such low rain rates over land very well, and because the Goddard scattering algorithm uses 1 mm/h as the minimum reliably detectable rain rate, regimes that contain large areas of very light rain present inherent difficulties for these retrieval methods. Therefore, the proliferation of low rain rates observed during the experiment is the main contributor to low correlation coefficients and high root-mean-square differences. Mididentification of cold surface (e.g., snow cover) as precipitation was also a problem in several instances.

Description: A procedure for the retrieval of cloud vertical structure from passive microwave radiometry is demonstrated by using passive microwave radiometry observations made during the Tropical Rainfall Measuring Mission. The procedure uses a set of cloud radiative models, with each model consisting of five vertical layers, specifying a distinct cloud vertical structure in terms of the near-surface parameters. The retrieval procedure is separated into two tasks (1) retrieving a set of geophysical parameters for each cloud radiative model and (2) finding which of the cloud radiative models and its associated retrieved parameters best fit the observed geophysical conditions. It is shown that this retrieval technique can detect differences and similarities between precipitating systems.

Description: Recent advances in cloud microphysical models have led to realistic three-dimensional distributions of cloud constituents. Radiative transfer schemes can make use of this detailed knowledge in order to study the effects of horizontal as well as vertical inhomogeneities within clouds. This study looks specifically at the differences between three-dimensional radiative transfer results and those obtained by plane parallel, independent pixel approximations in the microwave spectrum. A three-dimensional discrete ordinates method as well as a backward Monte Carlo method are used to calculate realistic radiances emerging from the cloud. Analyses between these models and independent pixel approximations reveal that plane parallel approximations introduce two distinct types of errors. The first error is physical in nature and is related to the fact that plane parallel approximations do not allow energy to leak out of dense areas into surrouding areas. In general, it was found that these errors are quite small for emission-dominated frequencies (37 GHz and lower) and that physical errors are highly pronounced only at scattering frequencies (85 GHz) where large deviations and biases up to 8 K averaged over the entire cloud were found. The second error is more geometric in nature and is related to the fact that plane parallel approximations cannot accommodate physical boundaries in the horizontal dimension for off-nadir viewing angles. The geometric errors were comparable in magnitude for all frequencies. Their magnitude, however, depends on a number of factors including the scheme used to deal with the edge, the nature of the surface, and the viewing angle.

Description: A multi channel physical approach for retrieving rainfall and its vertical structure from Special Sensor Microwave/Imager (SSM/I) observations is examined. While a companion paper was devoted exclusively to the description of the algorithm, its strengths, and its limitations, the main focus of this paper is to report on the results, applicability, and expected accuraciesfrom this algorithm. Some examples are given that compare retrieved results with ground-based radar data from different geographical regions to illustrate the performance and utility of the algorithm under distinct rainfall conditions. More quantitative validation is accomplished using two months of radar data from Darwin, Australia, and the radar network over Japan. Instantaneous comparisons at Darwin indicate that root-mean-square errors for 1.25 deg areas over water are 0.09 mm/h compared to the mean rainfall value of 0.224 mm/h while the correlation exceeds 0.9. Similar results are obtained over the Japanese validation site with rms errors of 0.615 mm/h compared to the mean of 0.0880 mm/h and a correlation of 0.9. Results are less encouraging over land with root-mean-square errors somewhat larger than the mean rain rates and correlations of only 0.71 and 0.62 for Darwin and Japan, respectively. These validation studies are further used in combination with the theoretical treatment of expected accuracies developed in the companion paper to define error estimates on a broader scale than individual radar sites from which the errors may be analyzed. Comparisons with simpler techniques that are based on either emission or scattering measurements are used to illustrate the fact that the current algorithm, while better correlated with the emission methods over water, cannot be reduced to either of these simpler methods.