Search

Data from an S-band polarimetric radar located at a mid-latitude, coastal location are used to compare two different methods for identifying stratiform and convective rain regions. The first method entails the retrievals of two (main) parameters of the rain drop size distributions using the radar reflectivity and the differential reflectivity. The second technique is a well-known texture-based method which utilizes the radar reflectivity and its spatial variability. A widespread event with embedded line convection was used as a test case. The two methods were compared using 500 m by 500 m pixel resolution gridded data constructed from the radar volume scans. Only 12% of the pixels showed disagreement between the two methods.

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.

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.

On 9 September 2019, rain-bands of category-1 Hurricane Dorian passed over a ground instrumentation site in Delmarva peninsula, USA. Drop shapes derived from 2D Video Disdrometer measurements at this site were used to compute the S-band radar cross sections (RCS) for horizontal and vertical polarizations for each drop with equi-volume diameter > 2 mm. These are combined with RCS for the smaller drops assuming equilibrium shapes. Radar reflectivity (Z H ) and differential reflectivity (Z DR ) are calculated for each of the 3 minutes throughout the event which lasted for more than 8 hours. These are compared with simultaneous observations from an S-band polarimetric radar 38 km away. The comparisons highlight the impact of large amplitude drop oscillations on Z DR .

Stratiform and convective rain are associated with different microphysical processes and generally produce drop size distributions (DSDs) with different characteristics. Identification of these two rain types is also important for estimating rainfall rates from ground-based polarimetric radars as well as spaceborne radars. Previous studies have investigated DSD characteristics using disdrometer data along with radar observations and/or vertically pointing Doppler radar observations. One such study, using data from Darwin, Australia (a tropical coastal 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. Since then, the separation method has been tested using data and observations from Greeley, Colorado, a mid-latitude continental location with semi-arid climate, and Huntsville, Alabama, a sub-tropical continental location. In this paper, we investigate the same separation technique using data and observations from a mid-latitude coastal region, situated in the Delmarva peninsula in Virginia at NASA Wallops Flight Facility (WFF). Two different types of disdrometers were used to construct the full DSD spectra and the NW versus Dm based classification is compared with simultaneous observations from a S-band polarimetric radar 38 km away from the disdrometer site. Three-minute DSDs were used for the classification and RHI (range-height indicator) radar scans over the disdrometer were used for validation. Events which occurred on 14 and 16 October 2019 were chosen. Results show, surprisingly, that there was no need to modify the separation criteria from those used in Darwin, AU. Also considered were the outer rain bands of Hurricane Dorian (as Category-1) which occurred on 06 September 2019. However, in this case, it was largely stratiform rain during the major part of the storm affecting WFF. Scattering (T-matrix) calculations using the 3-minute DSD spectra were used to derive retrieval equations for NW and Dm for the S-band radar data. These were applied to the radar scans to identify convective and stratiform rain regions as well as mixed or transition regions. Chosen RHI scans from two events (14 October 2019 and 06 September 2019) will be used as illustrative examples. Vertical profiles of reflectivity, differential reflectivity and copolar correlation coefficient over the disdrometer site will be extracted to establish whether or not the melting layer can be clearly distinguished. Specific times will be chosen from the two events and compared against the disdrometer data based classification as well for those corresponding times.

Hurricane rainbands are very efficient rain producers, but details on drop size distributions are still lacking. This study focuses on the rainbands of hurricane Dorian as they traversed the densely instrumented NASA precipitation-research facility at Wallops Island, VA, over a period of 8 h. Drop size distribution (DSD) was measured using a high-resolution meteorological particle spectrometer (MPS) and 2D video disdrometer, both located inside a double-fence wind shield. The shape of the DSD was examined using double-moment normalization, and compared with similar shapes from semiarid and subtropical sites. Dorian rainbands had a superexponential shape at small normalized diameter values similar to those of the other sites. NASA&rsquo, s S-band polarimetric radar performed range height-indicator (RHI) scans over the disdrometer site, showing some remarkable signatures in the melting layer (bright-band reflectivity peaks of 55 dBZ, a dip in the copolar correlation to 0.85 indicative of 12&ndash, 15 mm wet snow, and a staggering reflectivity gradient above the 0 &deg, C level of &ndash, 10 dB/km, indicative of heavy aggregation). In the rain layer at heights &lt, 2.5 km, polarimetric signatures indicated drop break-up as the dominant process, but drops as large as 5 mm were detected during the intense bright-band period.

Tropical storm Nate, which was a powerful hurricane prior to landfall along the US Gulf coast, traversed north and weakened considerably to a tropical depression as it moved near an instrumented site in Hunstville, AL. The outer rain bands lasted 18 h (03:00 to 21:00 UTC on 08 October 2017) and a 2D-video disdrometer (2DVD) captured the event which was shallow at times and indicative of pure warm rain processes. The 2DVD measurements are used for 3D reconstruction of drop shapes (including the rotationally asymmetric drops) and the drop-by-drop scattering matrix has been computed using Computer Simulation Technology integral equation solver for drop sizes &gt, 2.5 mm. From the scattering matrix elements, the polarimetric radar observables are simulated by integrating over 1 min consecutive segments of the event. These simulated values are compared with dual-polarized C-band radar data located at 15 km range from the 2DVD site to evaluate the contribution of the asymmetric drop shapes, specifically to differential reflectivity. The drop fall velocities and drop horizontal velocities in terms of magnitude and direction, all being derived from each drop image from two orthogonal cameras of the 2DVD, are also considered.

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

Two-dimensional video disdrometer (2DVD) data from a line convection rain event are analyzed using the method of moments surface integral equation (MoM-SIE) via drop-by-drop polarimetric scattering calculations at C band that are compared with radar measurements. Drop geometry of asymmetric drop shapes is reconstructed from 2DVD measurements, and the MoM-SIE model is created by meshing the surface of the drop. The differential reflectivity Zdr calculations for an example asymmetric drop are validated against an industry standard code solution at C band, and the azimuthal dependence of results is documented. Using the MoM-SIE analysis on 2DVD drop-by-drop data (also referred to as simply MoM-SIE), the radar variables [Zh, Zdr, Kdp, ρhv] are computed as a function of time (with 1-min resolution) and compared to C-band radar measurements. The importance of shape variability of asymmetric drops is demonstrated by comparing with the traditional (or “bulk”) method, which uses 1-min averaged drop size distributions and equilibrium oblate shapes. This was especially pronounced for ρhv, where the MoM-SIE method showed lowered values (dip) during the passage of the line convection consistent with radar measurements, unlike the bulk method. The MoM-SIE calculations of [Zh, Zdr, Kdp] agree very well with the radar measurements, whereas linear depolarization ratio (LDR) calculations from the drop-by-drop method are found to be larger than the values from the bulk method, which is consistent with the dip in simulated and radar-measured ρhv. Our calculations show the importance of the variance of shapes resulting from asymmetric drops in the calculation of ρhv and LDR.

Dual-polarization radar rainfall estimation relationships have been extensively tested in continental and subtropical coastal rain regimes, with little testing over tropical oceans where the majority of rain on Earth occurs. A 1.5-yr Indo-Pacific warm pool disdrometer dataset was used to quantify the impacts of tropical oceanic drop-size distribution (DSD) variability on dual-polarization radar variables and their resulting utility for rainfall estimation. Variables that were analyzed include differential reflectivity Zdr; specific differential phase Kdp; reflectivity Zh; and specific attenuation Ah. When compared with continental or coastal convection, tropical oceanic Zdr and Kdp values were more often of low magnitude (−1) and Zdr was lower for a given Kdp or Zh, consistent with observations of tropical oceanic DSDs being dominated by numerous, small, less-oblate drops. New X-, C-, and S-band R estimators were derived: R(Kdp), R(Ah), R(Kdp, ζdr), R(z, ζdr), and R(Ah, ζdr), which use linear versions of Zdr and Zh, namely ζdr and z. Except for R(Kdp), convective/stratiform partitioning was unnecessary for these estimators. All dual-polarization estimators outperformed updated R(z) estimators derived from the same dataset. The best-performing estimator was R(Kdp, ζdr), followed by R(Ah, ζdr) and R(z, ζdr). The R error was further reduced in an updated blended algorithm choosing between R(z), R(z, ζdr), R(Kdp), and R(Kdp, ζdr) depending on Zdr > 0.25 dB and Kdp > 0.3° km−1 thresholds. Because of these thresholds and the lack of hail, R(Kdp) was never used. At all wavelengths, R(z) was still needed 43% of the time during light rain (R < 5 mm h−1, Zdr < 0.25 dB), composing 7% of the total rain volume. As wavelength decreased, R(Kdp, ζdr) was used more often, R(z, ζdr) was used less often, and the blended algorithm became increasingly more accurate than R(z).

Analysis of drop size distributions (DSD) measured by collocated Meteorological Particle Spectrometer (MPS) and a third-generation, low-profile, 2D-video disdrometer (2DVD) are presented. Two events from two different regions (Greeley, Colorado, and Huntsville, Alabama) are analyzed. While the MPS, with its 50-μm resolution, enabled measurements of small drops, typically for drop diameters below about 1.1 mm, the 2DVD provided accurate measurements for drop diameters above 0.7 mm. Drop concentrations in the 0.7–1.1-mm overlap region were found to be in excellent agreement between the two instruments. Examination of the combined spectra clearly reveals a drizzle mode and a precipitation mode. The combined spectra were analyzed in terms of the DSD parameters, namely, the normalized intercept parameter NW, the mass-weighted mean diameter Dm, and the standard deviation of mass spectrum σM. The inclusion of small drops significantly affected the NW and the ratio σM/Dm toward higher values relative to using the 2DVD-based spectra alone. For each of the two events, polarimetric radar data were used to characterize the variation of radar-measured reflectivity Zh and differential reflectivity Zdr with Dm from the combined spectra. In the Greeley event, this variation at S band was well captured for small values of Dm (Zdr tended to 0 dB but Zh showed a noticeable decrease with decreasing Dm. For the Huntsville event, an overpass of the Global Precipitation Measurement mission Core Observatory satellite enabled comparison of satellite-based dual-frequency radar retrievals of Dm with ground-based DSD measurements. Small differences were found between the satellite-based radar retrievals and disdrometers.

Data analyses for the mobile Iowa X-band polarimetric (XPOL) radar from a long-duration rain event that occurred during the NASA Iowa Flood Studies (IFloodS) field campaign are presented. A network of six 2D video disdrometers (2DVDs) is used to derive four rain-rate estimators for the XPOL-5 radar. The rain accumulation validations with a collocated network of twin and triple tipping-bucket rain gauges have highlighted the need for combined algorithms because no single estimator was found to be sufficient for all cases considered. A combined version of weighted and composite algorithms is introduced, including a new R(Ah, Zdr) rainfall estimator for X band, where Ah is the specific attenuation for horizontal polarization and Zdr is the differential reflectivity. Based on measurement and algorithm errors, the weights are derived to be as piecewise constant functions over reflectivity values. The weights are later turned into continuous functions using smoothing splines. A methodology to derive the weights in near–real time is proposed for the composite-weighted algorithm. Comparisons of 2-h accumulations and 8-h event totals obtained from the XPOL-5 with 12 rain gauges have shown 10%–40% improvement in normalized bias over individual rainfall estimators. The analyses have enabled the development of rain-rate estimators for the Iowa XPOL.

Comprehensive analysis of an unusual graupel-shower event recorded by an S-band polarimetric radar and two optical-imaging surface instruments is presented. The primary radar characteristic was negative differential reflectivity Zdr values along a vertical column. During the afternoon hours of 16 February 2015, a sequence of three showers that were composed primarily of small (8–15-mm diameter) graupel affected the ground instrumentation site that was established for the Multi-Angle Snowflake Camera and Radar (MASCRAD) experiment in the high plains of Colorado. While these showers passed the instrumentation site, the CSU–CHILL radar conducted high-time-resolution (~2.5-min cycle time) range–height indicator (RHI) scans from a range of 13 km. The RHI data show that the negative Zdr values extended vertically through much of the reflectivity cores, implying that the reflectivity-weighted mean axis ratios of the graupel particles in this event remained somewhat prolate throughout their lifetime. To be specific, the cores of the convective showers only extended to heights of ~3.5 km AGL and had fractionally negative (from ~−0.3 to −0.7 dB) Zdr levels in those cores. Particle-image data obtained by the MASC system and by a collocated 2D video disdrometer measured the diameters, shapes, and fall speeds of the graupel particles as they reached the surface. The graupel particles were found to be primarily of the lump type with a slightly prolate mean shape (especially for the larger-diameter particles). Microwave backscatter calculations confirm that the graupel-particle shape and orientation characteristics are consistent with the observed slightly, but consistently, negative Zdr values.

Two-dimensional video disdrometer (2DVD) data were analyzed from two equatorial Indian (Gan) and west Pacific Ocean (Manus) islands where precipitation is primarily organized by the intertropical convergence zone and the Madden–Julian oscillation (MJO). The 18 (3.5) months of 2DVD data from Manus (Gan) Island show that 1) the two sites have similar drop size distribution (DSD) spectra of liquid water content, median diameter, rain rate R, radar reflectivity z, normalized gamma number concentration Nw, and other integral rain parameters; 2) there is a robust Nw-based separation between convective (C) and stratiform (S) DSDs at both sites that produces consistent separation in other parameter spaces. The 2DVD data indicate an equatorial, maritime average C/S rainfall accumulation fraction (frequency) of 81/19 (41/59) at these locations. It is hypothesized that convective fraction and frequency estimates are slightly higher than previous radar-based studies, because the ubiquitous weak, shallow convection (

Measurements using the 2-D video disdrometer (2DVD) taken during a heavy rainfall event in Huntsville, Alabama, are analysed. The 2DVD images were processed to derive the rain microstructure parameters for each individual drop, which in turn were used as input to the T-matrix method to compute the forward and back scatter amplitudes of each drop at C-band. The polarimetric radar variables were then calculated from the individual drop contribution over a finite time period, e.g., 1 min. The calculated co-polar reflectivity, differential reflectivity, specific differential propagation phase and the co-polar correlation coefficient were compared with measurements from a C-band polarimetric radar located 15 km away. An attenuation-correction method based on the specific differential propagation phase was applied to the co-polar and differential reflectivity data from the C-band radar, after ensuring accurate radar calibration. Time series comparisons of the parameters derived from the 2DVD and C-band radar data show very good agreement for all four quantities, the agreement being sometimes better than the computations using the 1-min drop size distribution and bulk assumptions on rain microstructure (such as mean shapes and model-based assumptions for drop orientation). The agreement is particularly improved in the case of co-polar correlation coefficient since this parameter is very sensitive to variation of shapes as well as orientation angles. The calculations mark the first attempt at utilizing experimentally derived "drop- by-drop" rain microstructure information to compute the radar polarimetric parameters and to demonstrate the value of utilizing the 2-D video disdrometer for studying rain microstructure under various precipitation conditions. Histograms of drop orientation angles as well as the most probable drop shapes and the corresponding variations were also derived and compared with prior results from the 80 m fall "artificial rain" experiment.

Polarimetric radar data obtained at high spatial and temporal resolutions offer a distinct advantage in estimating the spatial correlation function of drop size distribution (DSD) parameters and rain rate compared with a fixed gauge–disdrometer network. On two days during the 2011 Midlatitude Continental Convective Clouds Experiment (MC3E) campaign in Oklahoma, NASA’s S-band polarimetric radar (NPOL) performed repeated PPI scans every 40 s over six 2D video disdrometer (2DVD) sites, located 20–30 km from the radar. The two cases were 1) a rapidly evolving multicell rain event (with large drops) and 2) a long-duration stratiform rain event. From the time series at each polar pixel, the Pearson correlation coefficient is computed as a function of distance along each radial in the PPI scan. Azimuthal dependence is found, especially for the highly convective event. A pseudo-1D spatial correlation is computed that is fitted to a modified-exponential function with two parameters (decorrelation distance R0 and shape F). The first event showed significantly higher spatial variability in rain rate (shorter decorrelation distance R0 = 3.4 km) compared with the second event with R0 = 10.2 km. Further, for the second event, the spatial correlation of the DSD parameters and rain rate from radar showed good agreement with 2DVD-based spatial correlations over distances ranging from 1.5 to 7 km. The NPOL also performed repeated RHI scans every 40 s along one azimuth centered over the 2DVD network. Vertical correlations of the DSD parameters as well as the rainwater content were determined below the melting level, with the first event showing more variability compared with the second event.

We present observations and analyses for several cases with contrasting hydrometeor forms, recorded by several ground instruments at the MASCRAD (MASC + Radar) snow field site, radiosondes, and polarimetric radars. The instrumentation includes the multi-angle snowflake camera (MASC) and 2D-video disdrometer. The study focuses on microphysical characteristics analysis and is aimed at improving the radar-based winter-precipitation estimation.

Accurate measurements of rain drop size distributions (DSD), with particular emphasis on small and tiny drops, are presented. Measurements were conducted in two very different climate regions, namely Northern Colorado and Northern Alabama. Both datasets reveal a combination of (i) a drizzle mode for drop diameters less than 0.7 mm and (ii) a precipitation mode for larger diameters. Scattering calculations using the DSDs are performed at S and X bands and compared with radar observations for the first location. Our accurate DSDs will improve radar-based rain rate estimates as well as propagation predictions.

We present scattering models of snow and ice hydrometeors and computation of full polarimetrie radar variables for winter precipitation using a higher order method of moments (MoM) in the surface integral equation (SIE) formulation. The studies of winter precipitation are based primarily on measurements by a multi-angle snowflake camera (MASC), reconstruction of 3D hydrometeor shapes by means of the visual hull method, MoM-SIE scattering computations, and measured polarimetric observables by CSU-CHILL Radar.

The effect of ship motion on shipborne polarimetric radar measurements is considered at C band. Calculations are carried out by (i) varying the “effective” mean canting angle and (ii) separately examining the elevation dependence. Scattering from a single oblate hydrometeor is considered at first. Equations are derived (i) to convert the measured differential reflectivity for nonzero mean canting angles to those for zero mean canting angle and (ii) to do the corresponding corrections for nonzero elevation angles. Scattering calculations are also performed using the T-matrix method with measured drop size distributions as input. Dependence on mean volume diameter is examined as well as variations of the four main polarimetric parameters. The results show that as long as the ship movement is limited to a roll of less than about 10°–15°, the effects are tolerable. Furthermore, the results from the scattering simulations have been used to provide equations for correction factors that can be applied to compensate for the “apparent” nonzero canting angles and nonzero elevation angles, so that drop size distribution parameters and rainfall rates can be estimated without any bias.

A new 10-category, polarimetric-based hydrometeor identification algorithm (HID) for C band is developed from theoretical scattering simulations including wet snow, hail, and big drops/melting hail. The HID is applied to data from seven wet seasons in Darwin, Australia, using the polarimetric C-band (C-POL) radar, to investigate microphysical differences between monsoon and break periods. Scattering simulations reveal significant Mie effects with large hail (diameter > 1.5 cm), with reduced reflectivity and enhanced differential reflectivity Zdr and specific differential phase Kdp relative to those associated with S band. Wet snow is found to be associated with greatly depreciated correlation coefficient ρhv and moderate values of Zdr. It is noted that large oblate liquid drops can produce the same electromagnetic signatures at C band as melting hail falling quasi stably, resulting in some ambiguity in the HID retrievals. Application of the new HID to seven seasons of C-POL data reveals that hail and big drops/melting hail occur much more frequently during break periods than during monsoon periods. Break periods have a high frequency of vertically aligned ice above 12 km, suggesting the presence of strong electric fields. Reflectivity and mean drop diameter D0 statistics demonstrate that convective areas in both monsoon and break periods may have robust coalescence or melting precipitation ice processes, leading to enhanced reflectivity and broader distributions of D0. Conversely, for stratiform regions in both regimes, mean reflectivity decreases below the melting level, indicative of evaporative processes. Break periods also have larger ice water path fractions, indicating substantial mixed-phase precipitation generation as compared with monsoonal periods. In monsoon periods, a larger percentage of precipitation is produced through warm-rain processes.

Analysis of data from three snow events during the MASCRAD winter campaign in Greeley, Colorado, USA, are presented. The campaign entailed multi-angle snowflake camera (MASC), 2D-video disdrometer, precipitation occurrence sensor system, Pluvio gauge, weather station, all installed at the well-established MASCRAD (MASC + Radar) snow field site, as well as state-of-the-art S- and X-band polarimetric radars. We present here microphysical characteristics analysis for three heavy snowfall cases, on (a) 21–22 Feb 2015, (b) 15 Dec 2015, and (c) 01–02 Feb 2016.

We present our continued studies of winter precipitation within the MASCRAD (MASC + Radar) project, using multi-angle snowflake camera (MASC), 2D-video disdrometer, computational electromagnetic scattering methods, and state-of-the-art polarimetric radar. We also introduce some recent advancements to the observation and analysis process, and discuss new illustrative results.

We present our ongoing studies of rain precipitation synergistically using a 2D-video disdrometer, particle spectrometer, precipitation occurrence sensor system, Pluvio precipitation gauge, state-of-the-art polarimetric CSU-CHILL radar, and a higher order electromagnetic scattering method. We present and discuss measurements and analyses for several rain events in 2015 at MASCRAD Instrumentation Site in Colorado.

Propagation effects relating to rainfall remote sensing at X-band are examined using data from an event which occurred during a six-month rain measurement campaign in Greeley, Colorado. The campaign involved three disdrometers, a rain gauge, and polarimetric radar measurements at S and X bands. Attenuation correction schemes using well-established procedures are applied to the radar data. Use of the specific attenuation derived from the correction procedures for estimating rainfall rates is examined.

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

The algorithms used to estimate rainfall from polarimetric radar variables show significant variance in error characteristics over the range of naturally occurring rain rates. As a consequence, to improve rainfall estimation accuracy using polarimetric radar, it is necessary to optimally combine a number of different algorithms. In this study, a new composite method is proposed that weights the algorithms by the inverse of their theoretical error. A number of approaches are discussed and are investigated using simulated radar data calculated from disdrometer measurements. The resultant algorithms show modest improvement over composite methods based on decision-tree logic—in particular, at rain rates above 20 mm h−1.

Rain drop size distributions retrieved from polarimetric radar measurements over regularly occurring thunderstorms over the islands north of Darwin, Australia, are used to test if aerosol contributions to the probability distributions of the drop size distribution parameters (median volume diameter and normalized intercept parameter) are detectable. The observations reported herein are such that differences in cloud properties arising from thermodynamic differences are minimized but even so may be a factor. However, there is a clear signature that high aerosol concentrations are correlated with smaller number concentrations and larger drops. This may be associated with enhanced ice multiplication processes for low aerosol concentration storms or other processes such as invigoration of the updrafts.

In this paper raindrop shapes from laboratory, field and model investigations are examined in order to distill a consistent picture of raindrop axis ratios as a function of size. The first realistic measurements of raindrop shape were made in the UCLA vertical wind tunnel that was built under the supervision of Hans R. Pruppacher at UCLA in 1968. Subsequent research has investigated raindrop shape in more detail by laboratory studies, field measurements and modeling. Measurements of rain at night in backscattered light have revealed that raindrops larger than 1 mm diameter oscillate in conditions where winds and collisions are too weak to excite oscillations. One apparent source of these oscillations is a resonance with eddy shedding. Subsequent laboratory studies have demonstrated that such oscillations increase the average axis ratio of small raindrops. Details of past theoretical research show that raindrop axis ratios decrease with size in the manner of a sessile drop because of its increasing weight, but also because of an increasing differential in the fore-aft aerodynamic pressure. Modern analytical and numerical models predicted similar axis ratios for water drops falling at terminal velocity and these results are generally consistent with laboratory and field measurements. Although the axis ratios first measured in the UCLA wind tunnel and the more recent laboratory measurements have led to a broader understanding of raindrop axis ratios, new information about raindrop shape is available from recent field measurements using improved techniques. Raindrop axis ratios, orientation and canting angles using 2-D video disdrometer (2DVD) are also reported. The results include measurements from an artificial rain experiment conducted under calm wind conditions where the drops (over 115,000) were allowed to fall a distance of 80 m. Axis ratios and drop shapes show good agreement with the Mainz wind-tunnel data and images, as well as canting angle distributions symmetric about zero degrees with a standard deviations of 7–8° for all drops > 2 mm. Measurements of 2DVD raindrop shape has been collected from Okinawa, Sumatra, Toronto, Huntsville and Brisbane. An unusually wide axis ratio distribution was observed in the D = 1.6–1.8 mm region, which may be evidence for drop oscillations caused by eddy-shedding. Some subtle differences can also be seen for larger drops (e.g. D = 4 mm), the vertical chord being slightly larger than the horizontal chord leading to slightly larger axis ratios compared with the 80 m fall data. Additionally, an example using individual drop information from 2DVD measurements is considered, which enhances agreement with C-band polarimetric radar measurements. In this particular case, a considerable shift in the axis ratio distribution for 3.5 mm drops is observed during the initial strong convective phase of the event. Further ongoing work will combine the 2DVD measurements with polarimetric radar observations in a synergistic manner in order to identify the meteorological or microphysical conditions under which significant changes in drop shapes may occur in natural rain. Improved polarimetric radar estimates of canting angles are reported from an S-band system with high cross-polar performance antenna. The histogram of the standard deviation derived from the radar data obtained in light stratiform rain event with embedded convection shows a canting angle mode of about 7°, with a significant positive skewness. An increase in stability with drop diameter is also inferred from the radar measurements, which is in agreement with direct estimates from the 2DVD. Ongoing experiments in the Mainz wind tunnel are evaluating how the raindrop shapes are affected by collisions with drizzle drops within rain shafts.

This note builds on prior technique development related to the classification of rain types utilizing C-band polarimetric (CPOL) radar measurements. While the prior work was preliminary and limited in scope, the authors elaborate here on the basis of the drop size distribution (DSD)-based indexing technique for rain-type classification (convective/stratiform/mixed), and place it on firmer footing by testing the methodology against texture- and disdrometer-based methods as applied to Darwin datasets. A microphysical-based methodology is attractive as it links more directly to the underlying rainfall physical processes. Statistics of the DSD parameters, namely, histograms of log10(Nw) and D0, for convective and stratiform rain types across the premonsoon buildup and monsoon regimes were derived and further separated for over land and over ocean regions. The maximum value for mean D0 (1.64 mm) and the largest histogram standard deviation (0.32 mm) occurred for convective rain over land during the buildup regime. The largest differences in D0 and NW histograms were found to be for convective rain between the buildup and monsoon regimes (independent of land or ocean areas). Stratiform rain histograms were found to be very similar during the buildup regime with little land–ocean differences. However, somewhat larger land–ocean differences were found for the monsoon stratiform rain. The main histogram characteristics of the “mixed” or “uncertain” rain type were closer to the convective rain type than to stratiform, across both regimes and land–ocean areas. Additionally, the Nw versus D0 cluster of points (mean ±1σ) for convective rain agrees very well with the previously published range of values for maritime convective (equilibrium-like) DSDs.

Comparisons are made between the reflectivity Z, median volume diameter D0, and rain rate R from a dual-frequency profiler and the C-band polarimetric radar (C-POL), which are both located near Darwin, Australia. Examples from the premonsoon “buildup” regime and the monsoon (oceanic) regime are used to illustrate the excellent agreement between the dual-profiler retrievals and the polarimetric radar-based retrievals. This work builds on similar works that were limited in scope to shallow tropical showers and predominantly stratiform rain events. The dual-frequency profiler retrievals of D0 and R herein are based on ensemble statistics, whereas the polarimetric radar retrievals are based on algorithms derived by using one season of disdrometer data from Darwin along with scattering simulations. The latest drop shape versus D relation is used as well as the canting angle distribution results obtained from the 80-m fall bridge experiment in the scattering simulations. The scatterplot of D0 from dual-frequency profiler versus Zdr measurements from C-POL is shown to be consistent not only with the theoretical simulations and prior data but also within prior predicted error bars for both stratiform rain as well as convective rain. Based on dual-frequency profiler–retrieved gamma drop size distribution parameters, a new smoothly varying “separator” indexing scheme has been developed that classifies between stratiform and convective rain types, including a continuous “transition” region between the two. This indexing technique has been applied on a number of low-elevation-angle plan position indicator (PPI) sweeps with the C-POL from the two regime examples, to construct “unconditioned” histograms of D0 in stratiform and convective rain (to within the sensitivity of the radar). With reference to the latter, it is demonstrated that the distribution of D0 is different in the buildup example than in the monsoon example, because of the differences in both the microphysical and kinematic features between the two regimes. In particular, (i) the mean D0 is significantly larger in the buildup example than in the monsoon example, irrespective of rain type; (ii) the histogram width (or standard deviation) is much larger for the buildup example than the monsoon example, irrespective of rain type; and (iii) the histogram skewness is negative for the monsoon regime example because of a lack of larger D0 values, whereas the buildup histogram is positively skewed irrespective of rain type.

A decrease in copolar correlation coefficient (ρco) at C band has been observed for several rain events with broad drop size distributions (DSDs). Observational evidence comes from simultaneous measurements with a C-band dual-polarization radar and a 2D video disdrometer. The possibility of utilizing the ρco decrease for DSD retrievals is discussed. A preliminary method using the copolar reflectivity, differential reflectivity, and ρco is given for estimating the DSD parameters. Validation is carried out by deriving the differential propagation phase (Φdp) from the estimated DSD parameters and comparing against measurements. The method presented here shows potential but needs to be further assessed in different rain climatologies.

Drop shapes derived from a previously conducted artificial rain experiment using a two-dimensional video disdrometer (2DVD) are presented. The experiment involved drops falling over a distance of 80 m to achieve their terminal velocities as well as steady-state oscillations. The previous study analyzed the measured axis ratios (i.e., ratio of maximum vertical to maximum horizontal chord) as a function of equivolumetric spherical drop diameter (Deq) for over 115 000 drops ranging from 1.5 to 9 mm. In this paper, the actual contoured shapes of the drops are reported, taking into account the finite quantization limits of the instrument. The shapes were derived from the fast line-scanning cameras of the 2DVD. The drops were categorized into Deq intervals of 0.25-mm width and the smoothed contours for each drop category were superimposed on each other to obtain their most probable shapes and their variations due to drop oscillations. The most probable shapes show deviation from oblate spheroids for Deq > 4 mm, the larger drops having a more flattened base, in good agreement with the equilibrium (nonoblate) shape model of Beard and Chuang. Deviations were noted from the Beard and Chuang model shapes for diameters larger than 6 mm. However, the 2DVD measurements of the most probable contour shapes are the first to validate the Beard and Chuang model shapes for large drops, and further to demonstrate the differences from the equivalent oblate shapes. The purpose of this paper is to document the differences in radar polarization parameters and the range of error incurred when using the equivalent oblate shapes versus the most probable contoured shapes measured with the 2DVD especially for drop size distributions (DSDs) with large median volume diameters (>2 mm). The measured contours for Deq > 1.5 mm were fitted to a modified conical equation, and scattering calculations were performed to derive the complex scattering amplitudes for forward and backscatter for H and V polarizations primarily at 5.34 GHz (C band) but also at 3 GHz (S band) and 9 GHz (X band). Calculations were also made to derive the relevant dual-polarization radar parameters for measured as well as model-based drop size distributions. When comparing calculations using the contoured shapes against the equivalent oblate spheroid shapes, good agreement was obtained for cases with median volume diameter (D0) less than around 2 mm. Small systematic differences in the differential reflectivity (Zdr) values of up to 0.3 dB were seen for the larger D0 values when using the oblate shapes, which can be primarily attributed to the shape differences in the resonance region, which occurs in the 5.5–7-mm-diameter range at C band. Lesser systematic differences were present in the resonance region at X band (3–4 mm). At S band, the impact of shape differences in the polarimetric parameters were relatively minor for D0 up to 2.5 mm. Unusual DSDs with very large D0 values (>3 mm) (e.g., as can occur along the leading edge of severe convective storms or aloft due localized “big drop” zones) can accentuate the Zdr difference between the contoured shape and the oblate spheroid equivalent, especially at C band. For attenuation-correction schemes based on differential propagation phase, it appears that the equivalent oblate shape approximation is sufficient using a fit to the axis ratios from the 80-m fall experiment given in this paper. For high accuracy in developing algorithms for predicting D0 from Zdr, it is recommended that the fit to the most probable contoured shapes as given in this paper be used especially at C band.

We present our ongoing studies of winter precipitation using multi-angle snowflake camera (MASC), 2D-video disdrometer, computational electromagnetic scattering methods, and state-of-the-art polarimetric radars. The newly built and established MASCRAD (MASC + Radar) Snow Field Site is one of the currently best instrumented and most sophisticated field sites for winter precipitation measurements and analysis in the nation. We present and discuss MASCRAD measurements for the snow event on Nov 15, 2014 in La Salle, Colorado.

This paper presents an inter-comparison of rainfall parameters (median volume diameter and rain rate) using C-band polarimetric radar, a 2D-video disdrometer and a 400 MHz profiler for the Baiu front event of 8-9 June 2005 in Okinawa, Japan. These instruments are part of the Okinawa Sub-Tropical Environment Remote Sensing Center, operated by the National Institute of Information and Communications Technology (NICT). The 2D-video disdrometer is used to derive the mean axis ratio of raindrops versus drop diameter, as well as the drop size distribution for the Baiu event. The data are then used to simulate various relations between polarimetric scattering parameters such as: specific attenuation (Ah), and specific differential attenuation (Adp), versus specific differential phase (Kdp) which are required to correct the measured reflectivity at horizontal polarization (Zh), and the differential reflectivity (Zdr) for rain attenuation. The 2D-video disdrometer data are also used to arrive at retrieval formulas for median volume diameter (D0) from radar Zdr and rain rate from radar Kdp.The intense Baiu event of 8-9 June 2005 was composed of heavy convective rain cells embedded in large areas of stratiform rain. The inter-comparison of D0 and rain rate (R) between instruments was conducted for 12 hours (03:00-07:00, 11:00-19:00 UTC on 8th June 2004). The C-band radar retrievals were found to be in excellent agreement with the 2D-video disdrometer for the entire period. The 400 MHz profiler retrievals of D0 and R were in good agreement with 2D-video disdrometer during the more steady rain periods, with more scatter observed during the heavier convective rain periods. These inter-comparisons demonstrate the accuracy of C-band polarimetric radar to retrieve important rainfall parameters, as well as the accurate correction for rain attenuation using differential propagation phase.

Results from an experiment to measure the drop shapes using a 2D video disdrometer (2DVD) are reported. Under calm conditions, drops were generated from a hose located on a bridge 80 m above ground, this height being sufficient to allow drop oscillations to reach a steady state. The disdrometer data had to be carefully processed so as to eliminate the drops mismatched by the instrument and to remove the system spreading function. The total number of drops analyzed was around 115 000. Their axis ratio distributions were obtained for diameters ranging from 1.5 to 9 mm. The mean axis ratio decreases with increasing drop diameter, in agreement with the upper bound of the Beard and Chuang equilibrium shape model. The inferred mode of oscillation appears to be dominated by the oblate–prolate axisymmetric mode for the diameter range of 1.5 to 9 mm. The mean axis ratio agrees well with two empirically fitted formulas reported in earlier studies. In addition, a linear fit was applied to the data for radar applications relating to rain retrievals from dual-polarization measurements. The 2DVD data taken in moderate stratiform rain were also analyzed in a similar way and the results agree with the artificially generated drop experiment, at least up to 4 mm. No data for larger diameters were available for stratiform precipitation. Finally, the fall velocity was examined in terms of drop diameter. The results closely follow an empirical formula fitted to the Gunn and Kinzer data as well as the Beard and Pruppacher data including a slight decrease in the terminal velocity with a diameter beyond 7 mm.

The calculation of fade margins required to meet the user-specified availability criteria needs to take into consideration a number of meteorological factors specific to the earth-station location. One such factor is the annual average of the freezing level height. Information on this height is available in ITU-R Recommendation P. 839-3, which provides contours that are generated on a 1.5° by 1.5° latitude by longitude resolution grid. This paper compares these heights with the bright-band heights (BBH) obtained from the precipitation radar on-board the TRMM (Tropical Rainfall Measuring Mission) satellite. Four years of TRMM radar data have been analysed, and their averages are compared within the latitude range covered by the satellite orbit (35°N to 35°S). Comparisons show that the radar bright-band heights typically occur 300 m below the 0°C isotherm heights. However, results also indicate that this difference may be latitude dependent. Also examined are the year-to-year variability and the seasonal variation. In the former case, only 1998 showed BBH values which are somewhat higher in tropical regions, this year having had one of the strongest El Nino events recorded. In the latter case, results show significant seasonal variation, becoming more pronounced at higher latitudes. Copyright © 2003 John Wiley & Sons, Ltd.

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.

Meteorological radars are unique instruments for the observation and measurement of precipitation, of which they can provide extensive information depending on the equipment characteristics. In particular, all radars provide reliable measurements over large areas (100 km in range and even more) and, at the same time, with high spatial detail (from some tenths of meters to some hundreds of meters).

[1] This paper presents the results and analyses of Doppler radar measurements taken as part of a long-term campaign conducted in Singapore. Eight months of data taken in the zenith-pointing mode have been analyzed, in order to examine the validity of the assumptions underlying the retrieval algorithm proposed for the precipitation radar on board the TRMM satellite. Stratiform and convective events were analyzed separately. The former exhibits a clear melting layer in the 4–5 km height region, defined by an enhancement in the radar reflectivity, an increase in the mean fall velocity from 1–2 m/s in the snow region just above the melting layer to over 8 m/s in the rain region just below it. The corresponding width increases from less than 1 m/s to greater than 1 m/s. The variation in all three parameters can be accurately represented by the Noncoalescence-Nonbreakup melting layer model used for the satellite radar retrieval algorithms. Also examined was the monthly variation of the reflectivity-height profiles. In the rain region below the melting layer, a consistent decrease in the radar reflectivity with decreasing height is seen, indicating that the drop size distribution (DSD) cannot be assumed constant over all heights. However, if only the parameter N0 of the DSD is assumed constant, then good agreement is obtained between the mean velocity estimated from the reflectivity data and the measured velocity throughout the height profile. The agreement is also obtained in the spectrum width profile in the rain region. Convective precipitation on the other hand, does not exhibit a melting layer, and its averaged dBZ and velocity profiles cannot be characterized even under the constant N0 assumption. Its averaged reflectivity-height profiles show more variability than stratiform precipitation. For such cases, factors such as up and down drafts and nonstandard atmosphere need to be taken into account. The implication for the so-called N0 adjustment method is that while it can be applied to stratiform precipitation, the convective events need additional information, such as dual-frequency radar measurements, in order to improve the accuracy of the rainfall estimates.

A range profiling method proposed for nadir-looking rain-radar has been investigated using aircraft measurements of typhoon at 10 GHz and 35 GHz. In order to take into account the effects of change in the hydrometeor phase along the radar beam, it was necessary to modify the original method. Instead of using an analytic expression for the retrieval of radar reflectivity, a numerical solution is described which enables one to include such effects as the melting layer contribution and temperature change in the rain. Results show a marked improvement in the estimated rainfall rates even when a fairly simple bright band model is used. The improvement occurs mainly in the retrieved Ka-band results and indicates the importance of the melting layer, particularly for high frequencies. >

Most prediction methods relating to propagation in rain and other hydrometeors take into consideration the height of the melting layer in order to separate the effects of rain from ice. Here, we present the annual variability of this height, from five years' of data from the precipitation radar on board the Tropical Rainfall Measurement Mission (TRMM) satellite. The variability is derived for latitudes ranging from −30° to +30°, and compared with the corresponding freezing heights derived from the ECMWF reanalysis data. The two sets of heights show a strong variation with latitude, with the tropics showing the narrowest distribution. Their differences also show a latitude dependence, and indicate that the melting layer thickness is more variable for the subtropics. The implications are examined with reference to hydrometeor scatter interference. The recommended global standard prediction procedure for this mode of interference is evaluated using the variable distribution of the melting layer heights and their effects are examined in terms of the overall interference signal levels at 11 GHz and 30 GHz. When compared with the results using the fixed distribution specified in the current recommendation, small but significant differences at the 0.01 time percentages are observed, particularly for the 30 GHz case. The evaluation of an enhanced version of the procedure shows that significant differences can be observed even at 11 GHz, thus demonstrating the need to use a latitude-dependent rain height variability for the full evaluation of rain scatter interference levels.

Meteorological radars that are capable of making both polarisation diversity and phase measurements require an accurate method for absolute power calibration to achieve their full potential. A technique is described which uses a combination of differential reflectivity and differential phase measurements in rain to achieve this. The method is demonstrated using data obtained with the 3 GHz Chilbolton radar, which was recently upgraded to measure differential phase and Doppler velocity. The sensitivity of this technique to the assumed drop size distribution is estimated, and shown to have a minor effect on the accuracy of the calibration.