A variational scheme for retrieving rainfall rate and hail reflectivity fraction from polarization radar

Polarization radar offers the promise of much more accurate rainfall-rate R estimates than are possible from radar reflectivity factor Z alone, not only by better characterization of the drop size distribution, but also by more reliable correction for attenuation and the identification of hail. However, practical attempts to implement retrieval algorithms have been hampered by the difficulty in coping with the inherent noise in the polarization parameters. In this paper, a variational retrieval scheme is described that overcomes these problems by employing a forward model for differential reflectivity Z(dr) and differential phase shift phi(dp) and iteratively refining the coefficient a in the relationship Z = aR(b) such that the difference between the forward model and the measurements is minimized in a least squares sense. Two methods are used to ensure that a varies smoothly in both range and azimuth. In range, a is represented by a set of cubic-spline basis functions; in azimuth, the retrieval at one ray is used as a constraint on the next. The result of this smoothing is that the retrieval is tolerant of random errors in Zdr of up to 1 dB and in phi(dp) of up to 5 degrees. Correction for attenuation is achieved simply and effectively by including its effects in the forward model. If hail is present then the forward model is unable to match the observations of Z(dr) and phi(dp) simultaneously. This enables a first pass of the retrieval to be used to identify the radar pixels that contain hail, followed by a second pass in which the fraction of the Z in those gates that is due to hail is retrieved, this time with the scheme being able to forward-model both Z(dr) and phi(dp) accurately. The scheme is tested on S-band radar data from southern England in cases of rain, spherical hail, oblate hail, and mixtures of rain and hail. It is found to be robust and stable, even in the presence of differential phase shift on backscatter.