skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Retrieval of Drop Size Distribution Parameters Using a Dual-Polarimetric Radar
The raindrop size distribution (DSD) is vital for applications such as quantitative precipitation estimation, understanding microphysical processes, and validation/improvement of two-moment bulk microphysical schemes. We trace the history of the DSD representation and its linkage to polarimetric radar observables from functional forms (exponential, gamma, and generalized gamma models) and its normalization (un-normalized, single/double-moment scaling normalized). The four-parameter generalized gamma model is a good candidate for the optimal representation of the DSD variability. A radar-based disdrometer was found to describe the five archetypical shapes (from Montreal, Canada) consisting of drizzle, the larger precipitation drops and the ‘S’-shaped curvature that occurs frequently in between the drizzle and the larger-sized precipitation. Similar ‘S’-shaped DSDs were reproduced by combining the disdrometric measurements of small-sized drops from an optical array probe and large-sized drops from 2DVD. A unified theory based on the double-moment scaling normalization is described. The theory assumes the multiple power law among moments and DSDs are scaling normalized by the two characteristic parameters which are expressed as a combination of any two moments. The normalized DSDs are remarkably stable. Thus, the mean underlying shape is fitted to the generalized gamma model from which the ‘optimized’ two shape parameters are obtained. The other moments of the distribution are obtained as the product of power laws of the reference moments M3 and M6 along with the two shape parameters. These reference moments can be from dual-polarimetric measurements: M6 from the attenuation-corrected reflectivity and M3 from attenuation-corrected differential reflectivity and the specific differential propagation phase. Thus, all the moments of the distribution can be calculated, and the microphysical evolution of the DSD can be inferred. This is one of the major findings of this article.  more » « less
Award ID(s):
1901585
PAR ID:
10406094
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Remote Sensing
Volume:
15
Issue:
4
ISSN:
2072-4292
Page Range / eLocation ID:
1063
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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. 
    more » « less
  2. Lupo, Anthony (Ed.)
    We examine several different features of DSDs based on data and observations from two mid-latitude coastal locations: (a) the Delmarva peninsula, USA, and (b) Incheon, South Korea. In each case, the full DSD spectra were obtained from two collocated disdrometers. Two events from location (a) and one event from location (b) are presented. For (a), observations and retrievals from NASA’s S-band polarimetric radar are included in the analyses as well as retrieved DSD parameters from the dual-wavelength precipitation radar onboard the Global Precipitation Measurement satellite. For (b), the disdrometer-based DSD data are compared with measurements from another sensor. Our main aim is to examine the underlying shape of the DSDs and their representation by the generalized gamma model. 
    more » « less
  3. Abstract Polarimetric coastal radar data are used to compare the rainfall characteristics of Hurricanes Harvey (2017) and Florence (2018). Intense rainfall was an infrequent yet important contributor to the total rainfall in Harvey, but its relative contribution varied spatially. The total rainfall over land maximized near the coast over Beaumont, TX, due to intense convection resulting from prolonged onshore flow downshear from the circulation center. Overall, polarimetric radar observations in Harvey show a dominance of high concentrations of small‐to‐medium drops, consistent with prior tropical cyclone studies. The microphysical characteristics were spatially and temporally inhomogeneous however, with larger drops more frequent on 27 August and higher number concentrations more frequent on 28 and 30 August. The polarimetric variables and raindrop characteristics observed during Florence share broad similarities to Harvey, but had reduced variability, fewer observations of stronger reflectivity and differential reflectivity, and a lower frequency of high number concentrations and medium‐sized drops. The radar data indicate Florence had reduced coverage of stronger convection compared to Harvey. We hypothesize that differences in storm motion, intensity decay rates, and vertical wind shear produce the distinct precipitation structures and microphysical differences seen in Harvey and Florence. 
    more » « less
  4. null (Ed.)
    Abstract The time preceding supercell tornadogenesis and tornadogenesis “failure” has been studied extensively to identify differing attributes related to tornado production or lack thereof. Studies from the Verification of the Origins of Rotation in Tornadoes Experiment (VORTEX) found that air in the rear-flank downdraft (RFD) regions of non- and weakly tornadic supercells had different near-surface thermodynamic characteristics than that in strongly tornadic supercells. Subsequently, it was proposed that microphysical processes are likely to have an impact on the resulting thermodynamics of the near-surface RFD region. One way to view proxies to microphysical features, namely drop size distributions (DSDs), is through use of polarimetric radar data. Studies from the second VORTEX used data from dual-polarization radars to provide evidence of different DSDs in the hook echoes of tornadic and non-tornadic supercells. However, radar-based studies during these projects were limited to a small number of cases preventing result generalizations. This study compiles 68 tornadic and 62 non-tornadic supercells using Weather Surveillance Radar–1988 Doppler (WSR-88D) data to analyze changes in polarimetric radar variables leading up to, and at, tornadogenesis and tornadogenesis failure. Case types generally did not show notable hook echo differences in variables between sets, but did show spatial hook echo quadrant DSD differences. Consistent with past studies, differential radar reflectivity factor (Z DR ) generally decreased leading up to tornadogenesis and tornadogenesis failure; in both sets, estimated total number concentration increased during the same times. Relationships between DSDs and the near-storm environment, and implications of results for nowcasting tornadogenesis, also are discussed. 
    more » « less
  5. null (Ed.)
    The Remote sensing of Electrification, Lightning, And Meso-scale/micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) and the Cloud, Aerosol, and Complex Terrain Interactions Experiment Proposal (CACTI) field campaigns provided an unprecedented thirteen-disdrometer dataset in Central Argentina during the Intensive (IOP, 15 November to 15 December 2018) and Extended (EOP, 15 October 2018 to 30 April 2019) Observational Periods. The drop size distribution (DSD) parameters and their variability were analyzed across the region of interest, which was divided into three subregions characterized by the differing proximity to the Sierras de Córdoba (SDC), in order to assess the impact of complex terrain on the DSD parameters. A rigorous quality control of the data was first performed. The frequency distributions of DSD-derived parameters were analyzed, including the normalized intercept parameter (logNw), the mean volume diameter (D0), the mean mass diameter (Dm), the shape parameter (μ), the liquid water content (LWC), and the rain rate (R). The region closest to the SDC presented higher values of logNw, lower D0, and higher μ, while the opposite occurred in the farthest region, i.e., the concentration of small drops decreased while the concentration of bigger drops increased with the distance to the east of the SDC. Furthermore, the region closest to the SDC showed a bimodal distribution of D0: the lower values of D0 were associated with higher values of logNw and were found more frequently during the afternoon, while the higher D0 were associated with lower logNw and occurred more frequently during the night. The data were analyzed in comparison to the statistical analysis of Dolan et al. 2018 and sorted according to the classification proposed in the cited study. The logNw-D0 and LWC-D0 two-dimensional distributions allowed further discussion around the applicability of other mid-latitude and global precipitation classification schemes (startiform/convection) in the region of interest. Finally, three precipitation case studies were analyzed with supporting polarimetric radar data in order to relate the DSD characteristics to the precipitation type and the microphysical processes involved in each case. 
    more » « less