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. 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.
- Award ID(s):
- 1661799
- NSF-PAR ID:
- 10387882
- Date Published:
- Journal Name:
- Atmosphere
- Volume:
- 12
- Issue:
- 8
- ISSN:
- 2073-4433
- Page Range / eLocation ID:
- 979
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
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.
-
Abstract Current bulk microphysical parameterization schemes underpredict precipitation intensities and drop size distributions (DSDs) during warm rain periods, particularly upwind of coastal terrain. To help address this deficiency, this study introduces a set of modifications, called RCON, to the liquid-phase (warm rain) parameterization currently used in the Thompson–Eidhammer microphysical parameterization scheme. RCON introduces several model modifications, motivated by evaluating simulations from a bin scheme, which together result in more accurate precipitation simulations during periods of warm rain. Among the most significant changes are 1) the use of a wider cloud water DSD of lognormal shape instead of the gamma DSD used by the Thompson–Eidhammer parameterization and 2) enhancement of the cloud-to-rain autoconversion parameterization. Evaluation of RCON is performed for two warm rain events and an extended period during the Olympic Mountains Experiment (OLYMPEX) field campaign of winter 2015/16. We show that RCON modifications produce more realistic precipitation distributions and rain DSDs than the default Thompson–Eidhammer configuration. For the multimonth OLYMPEX period, we show that rain rates, rainwater mixing ratios, and raindrop number concentrations were increased relative to the Thompson–Eidhammer microphysical parameterization, while concurrently decreasing raindrop diameters in liquid-phase clouds. These changes are consistent with an increase in simulated warm rain. Finally, real-time evaluation of the scheme from August 2021 to August 2022 demonstrated improved precipitation prediction over coastal areas of the Pacific Northwest.
Significance Statement Although the accurate simulation of warm rain is critical to forecasting the hydrology of coastal areas and windward slopes, many warm rain parameterizations underpredict precipitation in these locations. This study introduces and evaluates modifications to the Thompson–Eidhammer microphysics parameterization scheme that significantly improve the accuracy of rainfall prediction in those regions.