Ann. Geophys., 27, 555-567, 2009
http://www.ann-geophys.net/27/555/2009/ doi:10.5194/angeo-27-555-2009 © Author(s) 2009. This work is distributed under the Creative Commons Attribution 3.0 License. |

04 Feb 2009

Cooperative Institute for Research in Environmental Sciences (CIRES), Univ. Colorado, Boulder, CO 80309-0216, USA

Earth System Research Laboratory (ESRL), NOAA, Boulder, CO 80305, USA

Abstract. Before radar estimates of the raindrop size distribution (DSD) can be assimilated into numerical weather prediction models, the DSD estimate must also include an uncertainty estimate. Ensemble statistics are based on using the same observations as inputs into several different models with the spread in the outputs providing an uncertainty estimate. In this study, Doppler velocity spectra from collocated vertically pointing profiling radars operating at 50 and 920 MHz were the input data for 42 different DSD retrieval models. The DSD retrieval models were perturbations of seven different DSD models (including exponential and gamma functions), two different inverse modeling methodologies (convolution or deconvolution), and three different cost functions (two spectral and one moment cost functions).

Two rain events near Darwin, Australia, were analyzed in this study producing 26 725 independent ensembles of mass-weighted mean raindrop diameter*D*_{m} and rain rate *R*. The mean and the standard deviation (indicated by the symbols <*x*> and σ*x*) of *D*_{m} and *R* were estimated for each ensemble. For small ranges of <*D*_{m}> or <*R*>, histograms of σ*D*_{m} and σ*R* were found to be asymmetric, which prevented Gaussian statistics from being used to describe the uncertainties. Therefore, 10, 50, and 90 percentiles of σ*D*_{m} and σ*R* were used to describe the uncertainties for small intervals of <*D*_{m}> or <*R*>. The smallest *D*_{m} uncertainty occurred for <*D*_{m}> between 0.8 and 1.8 mm with the 90th and 50th percentiles being less than 0.15 and 0.11 mm, which correspond to relative errors of less than 20% and 15%, respectively. The uncertainty increased for smaller and larger <*D*_{m}> values. The uncertainty of *R* increased with <*R*>. While the 90th percentile uncertainty approached 0.6 mm h^{−1} for a 2 mm h^{−1} rain rate (30% relative error), the median uncertainty was less than 0.15 mm h^{−1} at the same rain rate (less than 8% relative error). This study addresses retrieval error and does not attempt to quantify absolute or representativeness errors.

**Citation:**
Williams, C. R. and Gage, K. S.: Raindrop size distribution variability estimated using ensemble statistics, Ann. Geophys., 27, 555-567, doi:10.5194/angeo-27-555-2009, 2009.

Earth System Research Laboratory (ESRL), NOAA, Boulder, CO 80305, USA

Abstract. Before radar estimates of the raindrop size distribution (DSD) can be assimilated into numerical weather prediction models, the DSD estimate must also include an uncertainty estimate. Ensemble statistics are based on using the same observations as inputs into several different models with the spread in the outputs providing an uncertainty estimate. In this study, Doppler velocity spectra from collocated vertically pointing profiling radars operating at 50 and 920 MHz were the input data for 42 different DSD retrieval models. The DSD retrieval models were perturbations of seven different DSD models (including exponential and gamma functions), two different inverse modeling methodologies (convolution or deconvolution), and three different cost functions (two spectral and one moment cost functions).

Two rain events near Darwin, Australia, were analyzed in this study producing 26 725 independent ensembles of mass-weighted mean raindrop diameter