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

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 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

Accurate characterization of winter precipitation using multi-angle snowflake camera, visual hull, advanced scattering methods and polarimetric radar

Includes bibliographical references (pages 28-31).This article proposes and presents a novel approach to the characterization of winter precipitation and modeling of radar observables through a synergistic use of advanced optical disdrometers for microphysical and geometrical measurements of ice and snow particles (in particular, a multi-angle snowflake camera-MASC), image processing methodology, advanced method-of-moments scattering computations, and state-of-the-art polarimetric radars. The article also describes the newly built and established MASCRAD (MASC + Radar) in-situ measurement site, under the umbrella of CSU-CHILL Radar, as well as the MASCRAD project and 2014/2015 winter campaign. We apply a visual hull method to reconstruct 3D shapes of ice particles based on high-resolution MASC images, and perform "particle-by-particle" scattering computations to obtain polarimetric radar observables. The article also presents and discusses selected illustrative observation data, results, and analyses for three cases with widely-differing meteorological settings that involve contrasting hydrometeor forms. Illustrative results of scattering calculations based on MASC images captured during these events, in comparison with radar data, as well as selected comparative studies of snow habits from MASC, 2D video-disdrometer, and CHILL radar data, are presented, along with the analysis of microphysical characteristics of particles. In the longer term, this work has potential to significantly improve the radar-based quantitative winter-precipitation estimation.Published with support from the Colorado State University Libraries Open Access Research and Scholarship Fund

Spatial Correlation of Rain Drop Size Distribution from Polarimetric Radar and 2D-Video Disdrometers

Spatial correlations of two of the main rain drop-size distribution (DSD) parameters - namely the median-volume diameter (Do) and the normalized intercept parameter (Nw) - as well as rainfall rate (R) are determined from polarimetric radar measurements, with added information from 2D video disdrometer (2DVD) data. Two cases have been considered, (i) a widespread, long-duration rain event in Huntsville, Alabama, and (ii) an event with localized intense rain-cells within a convection line which occurred during the MC3E campaign. For the first case, data from a C-band polarimetric radar (ARMOR) were utilized, with two 2DVDs acting as ground-truth , both being located at the same site 15 km from the radar. The radar was operated in a special near-dwelling mode over the 2DVDs. In the second case, data from an S-band polarimetric radar (NPOL) data were utilized, with at least five 2DVDs located between 20 and 30 km from the radar. In both rain event cases, comparisons of Do, log10(Nw) and R were made between radar derived estimates and 2DVD-based measurements, and were found to be in good agreement, and in both cases, the radar data were subsequently used to determine the spatial correlations For the first case, the spatial decorrelation distance was found to be smallest for R (4.5 km), and largest fo Do (8.2 km). For log10(Nw) it was 7.2 km (Fig. 1). For the second case, the corresponding decorrelation distances were somewhat smaller but had a directional dependence. In Fig. 2, we show an example of Do comparisons between NPOL based estimates and 1-minute DSD based estimates from one of the five 2DVDs

Polarization Weather Radar Development from 1970–1995: Personal Reflections

The modern era of polarimetric radar begins with radiowave propagation research starting in the early 1970s with applications to measurement and modeling of wave attenuation in rain and depolarization due to ice particles along satellite–earth links. While there is a rich history of radar in meteorology after World War II, the impetus provided by radiowave propagation requirements led to high-quality antennas and feeds. Our journey starts by describing the key institutions and personnel responsible for development of weather radar polarimetry. The early period was dominated by circularly polarized radars for propagation research and at S band (frequency near 3 GHz) for hail detection. By the mid to late 70s, a paradigm shift occurred which led to the dominance of linear polarizations with applications to slant path attenuation prediction as well as estimation of rain rates and inferences of precipitation physics. The period from the early 1980s to 1995 can be considered as the “golden” period of rapid research that brought in meteorologists, cloud physicists, and hydrologists. This article describes the evolution of this technology from the vantage point of the authors. Their personal reflections and “behind the scenes” descriptions offer a glimpse into the inner workings at several key institutions which cannot be found elsewhere

Measurements of Rainfall Rate, Drop Size Distribution, and Variability at Middle and Higher Latitudes: Application to the Combined DPR-GMI Algorithm

The Global Precipitation Measurement mission is a major U.S.–Japan joint mission to understand the physics of the Earth’s global precipitation as a key component of its weather, climate, and hydrological systems. The core satellite carries a dual-precipitation radar and an advanced microwave imager which provide measurements to retrieve the drop size distribution (DSD) and rain rates using a Combined Radar-Radiometer Algorithm (CORRA). Our objective is to validate key assumptions and parameterizations in CORRA and enable improved estimation of precipitation products, especially in the middle-to-higher latitudes in both hemispheres. The DSD parameters and statistical relationships between DSD parameters and radar measurements are a central part of the rainfall retrieval algorithm, which is complicated by regimes where DSD measurements are abysmally sparse (over the open ocean). In view of this, we have assembled optical disdrometer datasets gathered by research vessels, ground stations, and aircrafts to simulate radar observables and validate the scattering lookup tables used in CORRA. The joint use of all DSD datasets spans a large range of drop concentrations and characteristic drop diameters. The scaling normalization of DSDs defines an intercept parameter NW, which normalizes the concentrations, and a scaling diameter Dm, which compresses or stretches the diameter coordinate axis. A major finding of this study is that a single relationship between NW and Dm, on average, unifies all datasets included, from stratocumulus to heavier rainfall regimes. A comparison with the NW–Dm relation used as a constraint in versions 6 and 7 of CORRA highlights the scope for improvement of rainfall retrievals for small drops (Dm < 1 mm) and large drops (Dm > 2 mm). The normalized specific attenuation–reflectivity relationships used in the combined algorithm are also found to match well the equivalent relationships derived using DSDs from the three datasets, suggesting that the currently assumed lookup tables are not a major source of uncertainty in the combined algorithm rainfall estimates

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

研究船「みらい」で観測された海洋上の雨滴粒径分布のモデル化の試みと、その降雨特性による差異

http://www.godac.jamstec.go.jp/darwin/cruise/mirai/mr15-04/

Testing the Drop-Size Distribution-Based Separation of Stratiform and Convective Rain Using Radar and Disdrometer Data from a Mid-Latitude Coastal Region

Stratiform and convective rain are associated with different microphysical processes and generally produce drop-size distributions (DSDs) with different characteristics. Previous studies using data from (a) a tropical coastal location, (b) a mid-latitude continental location with semi-arid climate, and (c) a sub-tropical continental location, found that the two rain types could be separated in the NW–Dm space, where Dm is the mass-weighted mean diameter and NW is the normalized intercept parameter. In this paper, we investigate the same separation technique using data and observations from a mid-latitude coastal region. Three-minute DSDs from disdrometer measurements are used for the NW- versus Dm-based classification and are compared with simultaneous observations from an S-band polarimetric radar 38 km away from the disdrometer site. Specifically, RHI (range-height indicator) scans over the disdrometer were used for confirmation. Results show that there was no need to modify the separation criteria from previous studies. Three-minute DSDs from the same location were used as input to scattering calculations to derive retrieval equations for NW and Dm for the S-band radar using an improved technique and applied to the RHI scans to identify convective and stratiform rain regions. Two events are shown as illustrative examples

Dual-wavelength radar technique development for snow rate estimation : a case study from GCPEx

Quantitative precipitation estimation (QPE) of snowfall has generally been expressed in power-law form between equivalent radar reflectivity factor (Z(e)) and liquid equivalent snow rate (SR). It is known that there is large variability in the prefactor of the power law due to changes in particle size distribution (PSD), density, and fall velocity, whereas the variability of the exponent is considerably smaller. The dual-wavelength radar reflectivity ratio (DWR) technique can improve SR accuracy by estimating one of the PSD parameters (characteristic diameter), thus reducing the variability due to the prefactor. The two frequencies commonly used in dual-wavelength techniques are Ku- and Kabands. The basic idea of DWR is that the snow particle size-to-wavelength ratio is falls in the Rayleigh region at Ku-band but in the Mie region at Ka-band. We propose a method for snow rate estimation by using NASA D3R radar DWR and Ka-band reflectivity observations collected during a long-duration synoptic snow event on 30-31 January 2012 during the GCPEx (GPM Cold-season Precipitation Experiment). Since the particle mass can be estimated using 2-D video disdrometer (2DVD) fall speed data and hydrodynamic theory, we simulate the DWR and compare it directly with D3R radar measurements. We also use the 2DVD-based mass to compute the 2DVD-based SR. Using three different mass estimation methods, we arrive at three respective sets of Z-SR and SR(Z(h), DWR) relationships. We then use these relationships with D3R measurements to compute radar-based SR. Finally, we validate our method by comparing the D3R radar-retrieved SR with accumulated SR directly measured by a well-shielded Pluvio gauge for the entire synoptic event.Peer reviewe

Measurements and Modeling of the Full Rain Drop Size Distribution

The raindrop size distribution (DSD) is fundamental for quantitative precipitation estimation (QPE) and in numerical modeling of microphysical processes. Conventional disdrometers cannot capture the small drop end, in particular the drizzle mode which controls collisional processes as well as evaporation. To overcome this limitation, the DSD measurements were made using (i) a high-resolution (50 microns) meteorological particle spectrometer to capture the small drop end, and (ii) a 2D video disdrometer for larger drops. Measurements were made in two climatically different regions, namely Greeley, Colorado, and Huntsville, Alabama. To model the DSDs, a formulation based on (a) double-moment normalization and (b) the generalized gamma (GG) model to describe the generic shape with two shape parameters was used. A total of 4550 three-minute DSDs were used to assess the size-resolved fidelity of this model by direct comparison with the measurements demonstrating the suitability of the GG distribution. The shape stability of the normalized DSD was demonstrated across different rain types and intensities. Finally, for a tropical storm case, the co-variabilities of the two main DSD parameters (normalized intercept and mass-weighted mean diameter) were compared with those derived from the dual-frequency precipitation radar onboard the global precipitation mission satellite