ALBERT

Description: The significant ambiguities inherent in the determination of a particular vertical rain intensity profile from a given time profile of radar echo powers measured by a downward-looking (spaceborne or airborne) radar at a single attenuating frequency are well documented. Indeed, one already knows that by appropriately varying the parameters of the reflectivity-rain rate (Z-R) and/or attenuation-rain rate (k- R) relationships one can produce several substantially different rain-rate profiles that would produce the same radar power profile. Imposing the additional constraint that the path-averaged rain rate be a given fixed number does reduce the ambiguities but falls far short of eliminating them. While formulas to generate all mutually ambiguous rain-rate profiles from a given profile of received radar reflectivities have already been derived, there remains to be produced a quantitative measure to assess how likely each of these profiles is, what the appropriate "average" profile should be, and what the "variance" of these multiple solutions is. To do this, one needs to spell out the stochastic constraints that can allow us to make sense of the words "average" and "variance" in a mathematically rigorous way. Such a quantitative approach would be particularly well suited for such systems as the planned precipitation radar of the Tropical Rainfall Measuring Mission (TRMM). Indeed, one would then be able to use the radar reflectivities measured by the TRMM radar to estimate the rain-rate profile that would most likely have produced the measurements, as well as the uncertainty in the estimated rain rates as a function of range. Such an optimal approach is described in this paper.

Description: This paper addresses the problem of finding a parametric form for the raindrop size distribution (DSD) that(1) is an appropriate model for tropical rainfall, and (2) involves statistically independent parameters. Such a parameterization is derived in this paper. One of the resulting three "canonical" parameters turns out to vary relatively little, thus making the parameterization particularly useful for remote sensing applications. In fact, a new set of r drop-size-distribution-based Z-R and k-R relations is obtained. Only slightly more complex than power laws, they are very good approximations to the exact radar relations one would obtain using Mie scattering. The coefficients of the new relations are directly related to the shape parameters of the particular DSD that one starts with. Perhaps most important, since the coefficients are independent of the rain rate itself, the relations are ideally suited for rain retrieval algorithms.

Description: This paper describes a computationally efficient nearly optimal Bayesian algorithm to estimate rain (and drop size distribution) profiles, given a radar reflectivity profile at a single attenuating wavelength. In addition to estimating the averages of all the mutually ambiguous combinations of rain parameters that can produce the data observed, the approach also calculates the n-ns uncertainty in its estimates (this uncertainty thus quantifies "the amount of ambiguity" in the "solution"). The paper also describes a more general approach that can make estimates based on a radar reflectivity profile together with an approximate measurement of the path-integrated attenuation, or a radar reflectivity profile and a set of passive microwave brightness temperatures. This more general "combined" algorithm is currently being adapted for the Tropical Rainfall Measuring Mission.

Description: The significant ambiguities inherent in the determination of a particular vertical rain intensity profile from a given time profile of radar echo powers measured by a downward-looking (spaceborne or airborne) radar at a single attenuating frequency are well-documented. Indeed, one already knows that by appropriately varying the parameters of the reflectivity-rain-rate (Z - R) and/or attenuation-rain-rate (k - R) relationships, one can produce several substantially different hypothetical rain rate profiles which would have the same radar power profile. Imposing the additional constraint that the path-averaged rain-rate be a given fixed number does reduce the ambiguities but falls far short of eliminating them. While we now know how to generate as many mutually ambiguous rain-rate profiles from a given profile of received radar reflectivities as we like, there remains to produce a quantitative measure to assess how likely each of these profiles is, what the appropriate 'average' profile should be, and what the 'variance' of these multiple solutions is. Of course, in order to do this, one needs to spell out the stochastic constraints that can allow us to make sense of the words 'average' and 'variance' in a mathematically rigorous way. Such a quantitative approach would be particularly well-suited for such systems as the proposed Precipitation Radar of the Tropical Rainfall Measuring Mission (TRMM). Indeed, one would then be able to use the radar reflectivities measured by the TRMM radar from one particular look in order to estimate the most likely rain-rate profile that would have produced the measurements, as well as the uncertainty in the estimated rain-rates as a function of range. Such an optimal approach is described in this paper.

Description: In this paper, an analytical treatment of the atmospheric remote sensing problem of determining the raindrop size distribution (DSD) with a spaceborne multifrequency microwave nadir-looking radar system is presented. It is typically assumed that with two radar measurements at different frequencies one ought to be able to calculate two state variables of the DSD: a bulk quantity, such as the rain rate, and a distribution shape parameter. To determine if this nonlinear problem can indeed be solved, the DSD is modeled as a Gamma distribution and quadratic approximations to the corresponding radar-rain relations are used to examine the invertibility of the resulting system of equations in the case of two as well as three radar frequencies. From the investigation, it is found that for regions of DSD state space multiple solutions exist for two or even three different frequency radar measurements. This should not be surprising given the nonlinear coupled nature of the problem.

Description: The precipitation radar planned for the Tropical Rainfall Measuring Mission (TRMM) will be the first of its kind to measure vertical rainfall distributions from space. The TRMM radar will scan +/- 20 degrees across the nadir track. The range-gated backscattering powers over the entire scan swath will be measured, classified (rain versus no-rain), averaged, and processed to derive the rainfall rates. With this observation scheme, there are two major reasons why it is important to know the rain-perturbed backscattering coefficient of the surface background (tilde over sigma_0)...

Description: The first problem addressed concerns passive-microwave rain retrievals. Most current approaches start by building off-line a cloud-model-derived database. Given data, the retrieval algorithms search the database for the microwave temperatures "closest" to the observed data, then after some fine-tuning (performed in different ways by different implementations) the rain is estimated to be that which corresponds to the selected (and fine-tuned) set of database temperatures. These approaches have three drawbacks: they cannot properly take into account the ambiguities which arise from the fact that several rain scenarios can produce the same observed temperatures; they are quite inefficient since they require manipulating a large database along with often complex "fine-tuning" procedures; and they cannot refine their estimates if additional data is available. This past year we have derived closed formulae relating observed microwave brightness temperatures, T(sub b), and the underlying rain rates, R: average T(sub b) =f (rain) and average rain = g (T(sub b)), along with the corresponding covariance matrices. These results are sufficient to describe the conditional probabilities p(R/T(sub b)) and p(T(sub b)/R) to second order. Progress has also been made towards deriving a robust description of the rain drop size distribution (DSD). The widespread approach consisting in parameterizing the DSD as a gamma-distribution in terms of the drop diameter D suffers from the facts that, in reality, the DSD is not a smooth function of D and that the largely arbitrary Gamma model imposes unintended behavior, which has implications on any quantities derived from the DSD model. We have therefore developed a non-parametric yet practical description of the DSD, which is particularly well-suited for use in remote-sensing applications. The diagram on the left shows a comparison between an actual DSD sample and the truncated non-parametric representation. One figure shows the relation between radar reflectivity and rain rate derived using this representation. Validation of the Tropical Rainfall Measuring Mission (TRMM) radar-radiometer combined R and DSD algorithm is underway. This algorithm was designed to make optimal use of the instantaneous reflectivity profiles measured by the TRMM radar and the microwave brightness temperatures measured by the TRMM passive radiometer. So far, it appears to be the most reliable TRMM rain algorithm.