skip to main content

Title: Retrieval of lower-order moments of the drop size distribution using CSU-CHILL X-band polarimetric radar: a case study
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 more » 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. « less
Authors:
; ; ; ;
Award ID(s):
1901585
Publication Date:
NSF-PAR ID:
10232025
Journal Name:
Atmospheric Measurement Techniques
Volume:
13
Issue:
9
Page Range or eLocation-ID:
4727 to 4750
ISSN:
1867-8548
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 quadrantmore »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.« less
  2. Commonly used disdrometers tend not to accurately measure concentrations of very small drops in the raindrop size distribution (DSD), either through truncation of the DSD at the small-drop end or because of large uncertainties on these measurements. Recent studies have shown that, as a result of these inaccuracies, many if not most ground-based disdrometers do not capture the “drizzle mode” of precipitation, which consists of large concentrations of small drops and is often separated from the main part of the DSD by a shoulder region. We present a technique for reconstructing the drizzle mode of the DSD from “incomplete” measurements in which the drizzle mode is not present. Two statistical moments of the DSD that are well measured by standard disdrometers are identified and used with a double-moment normalized DSD function that describes the DSD shape. A model representing the double-moment normalized DSD is trained using measurements of DSD spectra that contain the drizzle mode obtained using collocated Meteorological Particle Spectrometer and 2D video disdrometer instruments. The best-fitting model is shown to depend on temporal resolution. The result is a method to estimate, from truncated or uncertain measurements of the DSD, a more complete DSD that includes the drizzle mode.more »The technique reduces bias on low-order moments of the DSD that influence important bulk variables such as the total drop concentration and mass-weighted mean drop diameter. The reconstruction is flexible and often produces better rain-rate estimations than a previous DSD correction routine, particularly for light rain.

    « less
  3. Abstract

    Radar retrievals of drop size distribution (DSD) parameters are developed and evaluated over the mountainous Olympic Peninsula of Washington State. The observations used to develop retrievals were collected during the 2015/16 Olympic Mountain Experiment (OLYMPEX) and included the NASA S-band dual-polarimetric (NPOL) radar and a collection of second-generation Particle Size and Velocity (PARSIVEL2) disdrometers over the windward slopes of the barrier. Nonlinear and random forest regressions are applied to the PARSIVEL2 data to develop retrievals for median volume diameter, liquid water content, and rain rate. Improvement in DSD retrieval accuracy, defined by the mean error of the retrieval relative to PARSIVEL2 observations, was achieved when using the random forest model when compared with nonlinear regression. Evaluation of disdrometer observations and the retrievals from NPOL indicate that the radar retrievals can accurately reproduce observed DSDs in this region, including the common wintertime regime of small but numerous raindrops that is important there. NPOL retrievals during the OLYMPEX period are further evaluated using two-dimensional video disdrometers (2DVD) and vertically pointing Micro Rain Radars. Results indicate that radar retrievals using random forests may be skillful in capturing DSD characteristics in the lowest portions of the atmosphere.

  4. Two case studies of marine stratocumulus (one nocturnal and drizzling, the other daytime and nonprecipitating) are simulated by the UCLA large-eddy simulation model with bin microphysics for comparison with aircraft in situ observations. A high-bin-resolution variant of the microphysics is implemented for closer comparison with cloud drop size distribution (DSD) observations and a turbulent collision–coalescence kernel to evaluate the role of turbulence on drizzle formation. Simulations agree well with observational constraints, reproducing observed thermodynamic profiles (i.e., liquid water potential temperature and total moisture mixing ratio) as well as liquid water path. Cloud drop number concentration and liquid water content profiles also agree well insofar as the thermodynamic profiles match observations, but there are significant differences in DSD shape among simulations that cause discrepancies in higher-order moments such as sedimentation flux, especially as a function of bin resolution. Counterintuitively, high-bin-resolution simulations produce broader DSDs than standard resolution for both cases. Examination of several metrics of DSD width and percentile drop sizes shows that various discrepancies of model output with respect to the observations can be attributed to specific microphysical processes: condensation spuriously creates DSDs that are too wide as measured by standard deviation, which leads to collisional production of too manymore »large drops. The turbulent kernel has the greatest impact on the low-bin-resolution simulation of the drizzling case, which exhibits greater surface precipitation accumulation and broader DSDs than the control (quiescent kernel) simulations. Turbulence effects on precipitation formation cannot be definitively evaluated using bin microphysics until the artificial condensation broadening issue has been addressed.

    « less
  5. 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, themore »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.« less