Estimation of Rainfall Drop Size Distributions from Dual-Frequency Wind Profiler Spectra Using Deconvolution and a Nonlinear Least Squares Fitting Technique

Robert Schafer Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, Colorado

Search for other papers by Robert Schafer in
Current site
Google Scholar
PubMed
Close
,
Susan Avery Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, Colorado

Search for other papers by Susan Avery in
Current site
Google Scholar
PubMed
Close
,
Peter May Bureau of Meteorology Research Centre, Melbourne, Australia

Search for other papers by Peter May in
Current site
Google Scholar
PubMed
Close
,
Deepak Rajopadhyaya Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, Colorado

Search for other papers by Deepak Rajopadhyaya in
Current site
Google Scholar
PubMed
Close
, and
Christopher Williams Cooperative Institute for Research in Environmental Sciences, University of Colorado, and NOAA/Aeronomy Laboratory, Boulder, Colorado

Search for other papers by Christopher Williams in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

This paper compares two different analysis approaches for dual-frequency (VHF and UHF) wind profiler precipitation retrievals. The first technique is based on a general deconvolution process to remove broadening effects from the Doppler spectra due to turbulence and finite radar beamwidth. The second technique is based on a nonlinear least squares fitting approach, where an analytical function that describes the drop size distribution is convolved with a clear air distribution to model the broadening of the observed spectra. These techniques are tested with simulated data, where the true drop size distribution is known and represented by a gamma or exponential distribution, and with observations. The median volume diameter D0, which is based on the third moment of the drop size number distribution, is used to test these retrieval techniques. This parameter is closely related to the rain rate, which is approximately the 3.67th moment of the drop size number distribution.

Methods for extracting D0 from deconvolved and convolved spectra performed equally well up to D0 = 2.25 mm and a spectral width σ of about 0.5 m s−1, while for D0 > 2.25 mm a nonlinear least squares fit to the spectra gave better results. At σ = 1.5 m s−1, a direct retrieval of D0 from the deconvolved spectra, the method of moments applied to the deconvolved spectra, and a nonlinear least squares fit to the convolved spectra gave the best results, with retrieval success rates of about 80% for most D0. At the largest simulated spectral width, σ = 3 m s−1, the deconvolution methods performed best up to about D0 = 2 mm, with a success rate of about 60%, while the nonlinear least squares fit gave better results for D0 > 2 mm. An application of these dual-frequency retrieval techniques to real measurements taken during rain events shows a median absolute error (MAE) of less than 0.15 mm between D0 calculated by deconvolution methods, and D0 calculated by the convolution fitting approach. Comparison of retrieved rain rate and surface rain gauge measurements showed an MAE of less than 1.25 mm h−1 between the retrieval methods and the rain gauge, and less than 0.5 mm h−1 between the dual-frequency retrieval methods. These results suggest that the gamma distribution is a good representation of the drop size distribution during this rain event.

Corresponding author address: Dr. Robert Schafer, University of Colorado, CIRES Bldg., Room 318, Boulder, CO 80309-0216. Email: robert.schafer@colorado.edu

Abstract

This paper compares two different analysis approaches for dual-frequency (VHF and UHF) wind profiler precipitation retrievals. The first technique is based on a general deconvolution process to remove broadening effects from the Doppler spectra due to turbulence and finite radar beamwidth. The second technique is based on a nonlinear least squares fitting approach, where an analytical function that describes the drop size distribution is convolved with a clear air distribution to model the broadening of the observed spectra. These techniques are tested with simulated data, where the true drop size distribution is known and represented by a gamma or exponential distribution, and with observations. The median volume diameter D0, which is based on the third moment of the drop size number distribution, is used to test these retrieval techniques. This parameter is closely related to the rain rate, which is approximately the 3.67th moment of the drop size number distribution.

Methods for extracting D0 from deconvolved and convolved spectra performed equally well up to D0 = 2.25 mm and a spectral width σ of about 0.5 m s−1, while for D0 > 2.25 mm a nonlinear least squares fit to the spectra gave better results. At σ = 1.5 m s−1, a direct retrieval of D0 from the deconvolved spectra, the method of moments applied to the deconvolved spectra, and a nonlinear least squares fit to the convolved spectra gave the best results, with retrieval success rates of about 80% for most D0. At the largest simulated spectral width, σ = 3 m s−1, the deconvolution methods performed best up to about D0 = 2 mm, with a success rate of about 60%, while the nonlinear least squares fit gave better results for D0 > 2 mm. An application of these dual-frequency retrieval techniques to real measurements taken during rain events shows a median absolute error (MAE) of less than 0.15 mm between D0 calculated by deconvolution methods, and D0 calculated by the convolution fitting approach. Comparison of retrieved rain rate and surface rain gauge measurements showed an MAE of less than 1.25 mm h−1 between the retrieval methods and the rain gauge, and less than 0.5 mm h−1 between the dual-frequency retrieval methods. These results suggest that the gamma distribution is a good representation of the drop size distribution during this rain event.

Corresponding author address: Dr. Robert Schafer, University of Colorado, CIRES Bldg., Room 318, Boulder, CO 80309-0216. Email: robert.schafer@colorado.edu

Save
  • Atlas, D., Srivastava R. S. , and Sekhon R. S. , 1973: Doppler radar characteristics of precipitation at vertical incidence. Rev. Geophys. Space Sci., 11 , 135.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Beard, K. V., 1985: Simple altitude adjustments to raindrop velocities for Doppler radar analysis. J. Atmos. Oceanic Technol., 2 , 468471.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Cifelli, R., and Rutledge S. A. , 1994: Vertical motion structure in Maritime Continent mesoscale convective systems: Results from a 50-MHz profiler. J. Atmos. Sci., 51 , 26312652.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Cifelli, R., . 1998: Vertical motion, diabatic heating, and rainfall characteristics in N. Australia convective systems. Quart. J. Roy. Meteor. Soc., 124 , 11331162.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Currier, P. E., Avery S. K. , Balsley B. B. , and Gage K. S. , 1992: Use of two wind profilers for precipitation studies. Geophys. Res. Lett., 19 , 10171020.

    • Search Google Scholar
    • Export Citation
  • Doviak, R. J., and Zrnić D. S. , 1984: Doppler Radar and Weather Observations. Academic Press, 458 pp.

  • Gossard, E. E., 1988: Measuring drop-size distributions in clouds with a clear-air sensing Doppler radar. J. Atmos. Oceanic Technol., 5 , 640649.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hildebrand, P. H., and Sekhon R. S. , 1974: Objective determination of the noise level in Doppler spectra. J. Appl. Meteor., 13 , 808811.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Maguire, W. B. I. I,, and Avery S. K. , 1994: Retrieval of raindrop size distributions using two Doppler wind profilers: Model sensitivity testing. J. Appl. Meteor., 33 , 16231635.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • May, P. T., and Rajopadhyaya D. K. , 1996: Wind profiler observations of vertical motion and precipitation microphysics of a tropical squall line. Mon. Wea. Rev., 124 , 621633.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • May, P. T., Jameson A. R. , Keenan T. D. , and Johnston P. E. , 2001: A comparison between polarimetric radar and wind profiler observations of precipitation in tropical showers. J. Appl. Meteor., 40 , 17021717.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Press, W. H., Flannery B. P. , Teukolsky S. A. , and Vetterling W. T. , 1992: Numerical Recipes: The Art of Scientific Computing. Cambridge University Press, 963 pp.

    • Search Google Scholar
    • Export Citation
  • Rajopadhyaya, D. K., May P. T. , and Vincent R. A. , 1993: A general approach to the retrieval of rain drop-size distributions from VHF wind profiler Doppler spectra: Modeling results. J. Atmos. Oceanic Technol., 10 , 710717.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rajopadhyaya, D. K., Cifelli R. , Avery S. K. , Williams C. , and Ecklund W. L. , 1998: The effect of vertical wind motions on rain rates and median volume diameter determined from combined UHF and VHF wind profiler measurements and comparisons with rain gauge measurements. J. Atmos. Oceanic Technol., 15 , 13061319.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rajopadhyaya, D. K., Avery S. A. , May P. T. , and Cifelli R. C. , 1999: Comparison of precipitation estimation using single and dual-frequency wind profilers. J. Atmos. Oceanic Technol., 16 , 165173.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rogers, R. R., 1967: Doppler radar investigations of Hawaiian rain. Tellus, 19 , 432455.

  • Rogers, R. R., Baumgardner D. , Ether S. A. , Carter D. A. , and Ecklund W. L. , 1993: Comparison of raindrop size distributions measured by radar wind profiler and by airplane. J. Appl. Meteor., 32 , 694699.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Sauvageot, H., 1992: Radar Meteorology, Artech House, 366 pp.

  • Sato, T., Hiroshi D. , Hisato I. , Iwane K. , Shoichiro F. , Mamoru Y. , Toshitaka T. , and Susumu K. , 1990: Computer processing for deriving drop-size distributions and vertical air velocities from VHF Doppler radar spectra. Radio Sci., 25 , 961973.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ulbrich, C. W., 1983: Natural variations in the analytical form of the raindrop size distribution. J. Climate Appl. Meteor., 22 , 17641775.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ulbrich, C. W., and Atlas D. , 1998: Rainfall microphysics and radar properties: Analysis methods for drop size spectra. J. Appl. Meteor., 37 , 912923.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wakasugi, K., Mizutani A. , Matsuo M. , Fukao S. , and Kato S. , 1986: A direct method for deriving drop-size distribution and vertical air velocities from VHF Doppler radar spectra. J. Atmos. Oceanic Technol., 3 , 623629.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zrnić, D. S., 1975: Simulation of weatherlike Doppler spectra and signals. J. Appl. Meteor., 14 , 619620.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 540 300 58
PDF Downloads 153 50 13