scholarly journals Retrieval of lower-order moments of the drop size distribution using CSU-CHILL X-band polarimetric radar: a case study

2020 ◽  
Vol 13 (9) ◽  
pp. 4727-4750
Author(s):  
Viswanathan Bringi ◽  
Kumar Vijay Mishra ◽  
Merhala Thurai ◽  
Patrick C. Kennedy ◽  
Timothy H. Raupach

Abstract. The lower-order moments of the drop size distribution (DSD) have generally been considered difficult to retrieve accurately from polarimetric radar data because these data are related to higher-order moments. For example, the 4.6th moment is associated with a specific differential phase and the 6th moment with reflectivity and ratio of high-order moments with differential reflectivity. Thus, conventionally, the emphasis has been to estimate rain rate (3.67th moment) or parameters of the exponential or gamma distribution for the DSD. Many double-moment “bulk” microphysical schemes predict the total number concentration (the 0th moment of the DSD, or M0) and the mixing ratio (or equivalently, the 3rd moment M3). Thus, it is difficult to compare the model outputs directly with polarimetric radar observations or, given the model outputs, forward model the radar observables. This article describes the use of double-moment normalization of DSDs and the resulting stable intrinsic shape that can be fitted by the generalized gamma (G-G) distribution. The two reference moments are M3 and M6, which are shown to be retrievable using the X-band radar reflectivity, differential reflectivity, and specific attenuation (from the iterative correction of measured reflectivity Zh using the total Φdp constraint, i.e., the iterative ZPHI method). Along with the climatological shape parameters of the G-G fit to the scaled/normalized DSDs, the lower-order moments are then retrieved more accurately than possible hitherto. The importance of measuring the complete DSD from 0.1 mm onwards is emphasized using, in our case, an optical array probe with 50 µm resolution collocated with a two-dimensional video disdrometer with about 170 µm resolution. This avoids small drop truncation and hence the accurate calculation of lower-order moments. A case study of a complex multi-cell storm which traversed an instrumented site near the CSU-CHILL radar is described for which the moments were retrieved from radar and compared with directly computed moments from the complete spectrum measurements using the aforementioned two disdrometers. Our detailed validation analysis of the radar-retrieved moments showed relative bias of the moments M0 through M2 was <15 % in magnitude, with Pearson’s correlation coefficient >0.9. Both radar measurement and parameterization errors were estimated rigorously. We show that the temporal variation of the radar-retrieved mass-weighted mean diameter with M0 resulted in coherent “time tracks” that can potentially lead to studies of precipitation evolution that have not been possible so far.

2020 ◽  
Author(s):  
Viswanathan Bringi ◽  
Kumar Vijay Mishra ◽  
Merhala Thurai ◽  
Patrick C. Kennedy ◽  
Timothy H. Raupach

Abstract. The lower order moments of the drop size distribution (DSD) have generally been considered as difficult to retrieve accurately from polarimetric radar data because these are related to higher order moments. For example, the 4.5th moment is associated with specific differential phase, 6th moment with reflectivity and ratio of high order moments with differential reflectivity. Thus, conventionally, the emphasis has been to estimate rain rate (3.67th moment) or parameters of the exponential or gamma distribution. Many double-moment bulk microphysical schemes predict the total number concentration (the 0th moment or M0) and the mixing ratio (or equivalently, the 3rd moment M3). Thus, it is difficult to compare the model outputs directly with polarimetric radar observations or, given the model outputs, to forward model the radar observables. This article describes the use of double-moment normalization of DSDs and the resulting stable intrinsic shape that can be fitted to the generalized gamma (G-G) distribution. The two reference moments are M3 and M6 which are shown to be retrievable using the X-band radar reflectivity, differential reflectivity and specific attenuation (from the iterative ZPHI method). Along with the climatological shape parameters of the G-G fit to the scaled/normalized DSDs, the lower order moments are then retrieved more accurately than possible hitherto. The importance of measuring the complete DSD from 0.1 mm onwards is emphasized using, in our case, an optical array probe with 50 µm resolution collocated with a two-dimensional video disdrometer with 170 µm resolution. This avoids small drop truncation and hence the accurate calculation of lower order moments. A case study of a complex multi-cell storm which traversed an instrumented site near the CSU-CHILL radar is described for which the moments were retrieved and compared with directly computed moments from the complete spectrum measurements using the aforementioned two disdrometers. Our detailed validation analysis of the radar-retrieved moments showed relative bias of the moments M0 through M2 was  0.9. Both radar measurement and parameterization errors were estimated rigorously. We show that the temporal variation of the radar-retrieved characteristic diameter with M0 resulted in coherent time tracks that can potentially lead to studies of precipitation evolution that have not been possible so far.


2008 ◽  
Vol 25 (5) ◽  
pp. 729-741 ◽  
Author(s):  
Eugenio Gorgucci ◽  
V. Chandrasekar ◽  
Luca Baldini

Abstract The recent advances in attenuation correction methodology are based on the use of a constraint represented by the total amount of the attenuation encountered along the path shared over each range bin in the path. This technique is improved by using the inner self-consistency of radar measurements. The full self-consistency methodology provides an optimization procedure for obtaining the best estimate of specific and cumulative attenuation and specific and cumulative differential attenuation. The main goal of the study is to examine drop size distribution (DSD) retrieval from X-band radar measurements after attenuation correction. A new technique for estimating the slope of a linear axis ratio model from polarimetric radar measurements at attenuated frequencies is envisioned. A new set of improved algorithms immune to variability in the raindrop shape–size relation are presented for the estimation of the governing parameters characterizing a gamma raindrop size distribution. Simulations based on the use of profiles of gamma drop size distribution parameters obtained from S-band observations are used for quantitative analysis. Radar data collected by the NOAA/Earth System Research Laboratory (ESRL) X-band polarimetric radar are used to provide examples of the DSD parameter retrievals using attenuation-corrected radar measurements. Retrievals agree fairly well with disdrometer data. The radar data are also used to observe the prevailing shape of raindrops directly from the radar measurements. A significant result is that oblateness of drops is bounded between the two shape models of Pruppacher and Beard, and Beard and Chuang, the former representing the upper boundary and the latter the lower boundary.


2016 ◽  
Author(s):  
Timothy H. Raupach ◽  
Alexis Berne

Abstract. A new technique for estimating the raindrop size distribution (DSD) from polarimetric radar data is proposed. Two statistical moments of the DSD are estimated from polarimetric variables, and the DSD is reconstructed. The technique takes advantage of the relative invariance of the double-moment normalised DSD. The method was tested using X-band radar data and networks of disdrometers in three different climatic regions. Radar-derived estimates of the DSD compare reasonably well to observations. In the three tested domains, the proposed method performs similarly to and often better than a state-of-the-art DSD-retrieval technique. The approach is flexible because no specific double-normalised DSD model is prescribed. In addition, a method is proposed to treat noisy radar data to improve DSD-retrieval performance with radar measurements.


2012 ◽  
Vol 13 (3) ◽  
pp. 1066-1079 ◽  
Author(s):  
M. Thurai ◽  
V. N. Bringi ◽  
L. D. Carey ◽  
P. Gatlin ◽  
E. Schultz ◽  
...  

Abstract The accuracy of retrieving the two drop size distribution (DSD) parameters, median volume diameter (D0), and normalized intercept parameter (NW), as well as rain rate (R), from polarimetric C-band radar data obtained during a cool-season, long-duration precipitation event in Huntsville, Alabama, is examined. The radar was operated in a special “near-dwelling” mode over two video disdrometers (2DVD) located 15 km away. The polarimetric radar–based retrieval algorithms for the DSD parameters and rain rate were obtained from simulations using the 2DVD measurements of the DSD. A unique feature of this paper is the radar-based estimation of the spatial correlation functions of the two DSD parameters and rain rate that are used to estimate the “point-to-area” variance. A detailed error variance separation is performed, including the aforementioned point-to-area variance, along with variance components due to the retrieval algorithm error, radar measurement error, and disdrometer sampling error. The spatial decorrelation distance was found to be smallest for the R (4.5 km) and largest for D0 (8.24 km). For log10(NW), it was 7.22 km. The proportion of the variance of the difference between radar-based estimates and 2DVD measurements that could be explained by the aforementioned errors was 100%, 57%, and 73% for D0, log10(NW), and R, respectively. The overall accuracy of the radar-based retrievals for the particular precipitation event quantified in terms of the fractional standard deviation were estimated to be 6.8%, 6%, and 21% for D0, log10(NW), and R, respectively. The normalized bias was &lt;1%. These correspond to time resolution of ~3 min and spatial resolution of ~1.5 km.


Sign in / Sign up

Export Citation Format

Share Document