Understanding drop size distribution (DSD) variability has important implications for remote sensing and numerical modeling applications. Twelve disdrometer datasets across three latitude bands are analyzed in this study, spanning a broad range of precipitation regimes: light rain, orographic, deep convective, organized midlatitude, and tropical oceanic. Principal component analysis (PCA) is used to reveal comprehensive modes of global DSD spatial and temporal variability. Although the locations contain different distributions of individual DSD parameters, all locations are found to have the same modes of variability. Based on PCA, six groups of points with unique DSD characteristics emerge. The physical processes that underpin these groups are revealed through supporting radar observations. Group 1 (group 2) is characterized by high (low) liquid water content (LWC), broad (narrow) distribution widths, and large (small) median drop diameters D₀. Radar analysis identifies group 1 (group 2) as convective (stratiform) rainfall. Group 3 is characterized by weak, shallow radar echoes and large concentrations of small drops, indicative of warm rain showers. Group 4 identifies heavy stratiform precipitation. The low latitudes exhibit distinct bimodal distributions of the normalized intercept parameter N _w, LWC, and D₀ and are found to have a clustering of points (group 5) with high rain rates, large N _w, and moderate D₀, a signature of robust warm rain processes. A distinct group associated with ice-based convection (group 6) emerges in the midlatitudes. Although all locations exhibit the same covariance of parameters associated with these groups, it is likely that the physical processes responsible for shaping the DSDs vary as a function of location.

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 Z_dr and specific differential phase K_dp relative to those associated with S band. Wet snow is found to be associated with greatly depreciated correlation coefficient ρ_hv and moderate values of Z_dr. 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 D₀ 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 D₀. 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.

We have developed an algorithm that retrieves the size, number concentration and density of falling snow from multifrequency radar observations. This work builds on previous studies that have indicated that three-frequency radars can provide information on snow density, potentially improving the accuracy of snow parameter estimates. The algorithm is based on a Bayesian framework, using lookup tables mapping the measurement space to the state space, which allows fast and robust retrieval. In the forward model, we calculate the radar reflectivities using recently published snow scattering databases. We demonstrate the algorithm using multifrequency airborne radar observations from the OLYMPEX–RADEX field campaign, comparing the retrieval results to hydrometeor identification using ground-based polarimetric radar and also to collocated in situ observations made using another aircraft. Using these data, we examine how the availability of multiple frequencies affects the retrieval accuracy, and we test the sensitivity of the algorithm to the prior assumptions. The results suggest that multifrequency radars are substantially better than single-frequency radars at retrieving snow microphysical properties. Meanwhile, triple-frequency radars can retrieve wider ranges of snow density than dual-frequency radars and better locate regions of high-density snow such as graupel, although these benefits are relatively modest compared to the difference in retrieval performance between dual- and single-frequency radars. We also examine the sensitivity of the retrieval results to the fixed a priori assumptions in the algorithm, showing that the multifrequency method can reliably retrieve snowflake size, while the retrieved number concentration and density are affected significantly by the assumptions.

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud‐resolving models (CRMs). POLARRIS is composed of forward (POLARRIS‐f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables (Z_h, Z_dr, K_dp, ρ_hv) but also radar‐consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single‐moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables (Z_dr, K_dp) on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. The simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

Microphysical and kinematic characteristics of two storm populations, based on their macroscale charge structures, are investigated in an effort to increase our understanding of the processes that lead to anomalous (or inverted charge) structures. Nine normal polarity cases (midlevel negative charge) with dual‐Doppler and polarimetric coverage that occurred in northern Alabama and six anomalous polarity cases (midlevel positive charge) that occurred in northeastern Colorado are included in this study. The results show that even though anomalous polarity storms formed in environments with similar instability, they had significantly larger and stronger updrafts. Moreover, the anomalous polarity storms evidently have more robust mixed‐phase microphysics, based on a variety of metrics. We infer positively charged graupel, and therefore high supercooled water contents, in the midlevels of the anomalous storms based on the relationship between colocations of graupel and inferred positive charge from Lightning Mapping Array data. Anomalous polarity storms in Colorado have much higher cloud base heights and shallower warm cloud depths in this study, leading us to hypothesize that anomalous polarity storms have lower amounts of dilution and entrainment. Using representative updraft speeds and warm cloud depths, the time required for a parcel to traverse from cloud base to the freezing level was estimated for each storm observation. We suggest that this metric is the key discriminator between the two storm populations and leads us to hypothesize that it strongly influences the amount of supercooled water and the probability of positive charge in the midlevels, leading to an anomalous charge structure.

National Aeronautics and Space Administration's Investigations of Convective Updrafts (INCUS) mission aims to document convective mass flux through changes in the radar reflectivity (ΔZ) in convective cores captured by a constellation of three Ka‐band radars sampling the same convective cells over intervals of 30, 90, and 120 s. Here, high spatiotemporal resolution observations of convective cores from surface‐based radars that use agile sampling techniques are used to evaluate aspects of the INCUS measurement approach using real observations. Analysis of several convective cells confirms that large coherent ΔZ structure with measurable signal (>5 dB) can occur in less than 30 s and are correlated with underlying convective motions. The analysis indicates that the INCUS mission radar footprint and along track sampling are adequate to capture most of the desirable ΔZ signals. This unique demonstration of reflectivity time‐lapse provides the framework for estimating convective mass flux independent from Doppler techniques with future radar observations.