1. Introduction
Dual-polarization radar provides a unique set of variables that can be used to infer winter cloud and precipitation characteristics. A better understanding of such characteristics is essential for accurate snow forecasting (Reeves et al. 2014) and hydrometeor classification (Thompson et al. 2014). The recent advent of ubiquitous dual-polarization capabilities in operational NEXRAD instruments provides ample opportunity for developing appropriate methods to exploit these polarimetric capabilities. As a result, operational radars can help to improve snowfall forecasting capabilities and further enhance understanding of the microphysical particle properties and processes that cause dual-polarization signatures or “fingerprints” (Kumjian et al. 2016).
To complement the recent polarimetric advances in radar meteorology, newer studies have developed mathematical theories and methodologies that exploit multiple radar variables to retrieve parameters describing particle size distributions (PSDs; Ryzhkov and Zrnić 2019), particle characteristics such as shape and orientation (e.g., Melnikov and Straka 2013; Matrosov et al. 2017; Matrosov 2020; Matrosov et al. 2020), ice water content (IWC; Ryzhkov et al. 1998), and snow precipitation rates (S; Bukovčić et al. 2018, 2020). These microphysical retrievals are either purely empirical, thereby linking radar quantities directly to IWC or S through the use of (bivariate) power laws (e.g., Hogan et al. 2006; Bukovčić et al. 2018; Nguyen et al. 2019; Bukovčić et al. 2020), or estimated from PSD parameters that, in turn, are used to calculate IWC or S (e.g., Ryzhkov and Zrnić 2019). Microphysical relations that rely on PSD characterization must make numerous approximations and assumptions about the shapes, densities, orientations, and degree of riming for all particles (see Ryzhkov and Zrnić 2019). Such assumptions are very crude since ice particles are well known to exhibit a wide variety of microphysical characteristics even within the same ice particle class (Locatelli and Hobbs 1974; Mitchell 1996) and are expected to vary spatiotemporally depending upon environmental conditions (Bailey and Hallett 2009) and growth histories (Chen and Lamb 1994; Harrington et al. 2013; Sulia et al. 2013; Jensen and Harrington 2015). Despite this variability, current radar retrieval techniques often employ the same set of assumptions and approximations including a fixed (oblate) spheroid particle shape, a fixed density–size (ρs–D) relation, and a fixed canting angle standard deviation σ. Such assumptions necessitate their own sets of errors and uncertainties (Westbrook 2014; Schrom and Kumjian 2018), yet these assumptions have not yet been rigorously evaluated through careful analyses with in situ cloud probe observations.
Radar retrievals using reflectivity Z, differential reflectivity ZDR, specific differential phase Kdp, and copolar correlation coefficient ρhv should theoretically be most accurate in regions where all four variables are most prominent.1 One such scenario can regularly occur within the dendritic growth layer, approximately from −12° to −18°C, where ZDR and Kdp are elevated (0.5 < ZDR < 5.0 dB and 0.2 < Kdp < 0.8° km−1 at S band) and ρhv is diminished (0.96 < ρhv < 0.99). For Rayleigh scattering, Kdp depends upon the inverse of radar wavelength. As such, dendritic growth layer Kdp enhancements are roughly 3 times as high for X band as for S band. Enhanced ZDR and Kdp implies the presence of anisotropic particles such as dendrites or hexagonal plates where the magnitude of ZDR and Kdp depends upon the particle aspect ratios and densities. In addition, unlike ZDR, which is not sensitive to number concentration, Kdp is dependent upon number concentration and can therefore indicate enhancements in the number of particles as well. Diminished ρhv suggests a diversity of different particle types, shapes, densities, and/or orientations. Many studies have investigated these dendritic growth layer signatures and have provided various theories that explain the underlying microphysical causes (Ryzhkov and Zrnić 1998; Kennedy and Rutledge 2011; Andrić et al. 2013; Bechini et al. 2013; Grazioli et al. 2015; Moisseev et al. 2015; Schrom et al. 2015; Schrom and Kumjian 2016, 2018; Kumjian and Lombardo 2017; Griffin et al. 2018) yet there is no overall consensus for the validity or reconciliation of these theories. Some studies such as Kennedy and Rutledge (2011), Schrom et al. (2015), and Schrom and Kumjian (2018) have hypothesized that the Kdp and ZDR dendritic growth layer enhancements result from large dendrites that grow in situ, whereas other studies have hypothesized the role of “seeder–feeder”-type processes and oblate aggregates as a potential explanation (Moisseev et al. 2015). Furthermore, other studies such as Griffin et al. (2018) have provided evidence that hexagonal plates in the dendritic growth layer can cause extremely large values of ZDR, and they attributed ZDR and Kdp enhancements to separate types of particles. Still other studies such as Andrić et al. (2013) and Sinclair et al. (2016) attribute Kdp enhancements to secondary ice production of small, anisotropic crystals. Homogeneous oblate spheroids are often used in these studies to represent various ice particle habits; however, recent discrete dipole approximation modeling results (Westbrook 2014; Schrom and Kumjian 2018) illustrate that this approximation can severely underestimate or overestimate both ZDR and Kdp for branched planar crystals depending upon the assumed effective density.
Untangling the various causes of these polarimetric signatures requires careful consideration of collocated radar and in situ observations along with appropriate, physical approximations that govern each retrieval. However, few studies have directly compared dual-polarization radar retrieval estimates with in situ observations (e.g., Ryzhkov et al. 1998; Nguyen et al. 2019; Murphy et al. 2020). Such comparison efforts are essential to ensure retrievals are consistent with observations and to investigate any biases and uncertainties in such calculations. Indeed, each retrieval relies on physical assumptions about ice particle shapes, densities, orientations, and degree of riming that not only influence the PSD parameters but the resulting quantities (e.g., IWC or S) that often follow from such assumptions through weighted integration of the retrieved PSD. Any errors in radar quantities from, for instance, calibration biases or spatiotemporal matching will necessarily propagate into any retrieved quantity; it is our understanding that few studies have actually quantified how radar variable uncertainties would affect retrievals and how the magnitude of such uncertainties compare to the retrieved quantities themselves.
Toward this end, this paper evaluates current PSD and aspect ratio retrieval methodologies through direct comparison of Multi-Radar Multi-Sensor (MRMS; Zhang et al. 2016) observations collocated with in situ aircraft observations obtained during the 2020 Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Storms (IMPACTS; McMurdie et al. 2022) campaign. While we only evaluate two flight legs of this winter storm case (constituting approximately 1 h of flight data), the MRMS in situ collocation approach and subsequent investigation into fascinating polarimetric signatures can be thought of as a first attempt at understanding such radar features through the combined use of multiple observational platforms. This study is the first to incorporate polarimetric MRMS data and MRMS uncertainties in the use of radar retrieval algorithms and the first to attempt collocation of MRMS observations with measurements obtained from in situ aircraft data.
2. The MRMS system
In this study, the MRMS mosaic of NEXRAD WSR-88D instruments is used to collocate ice microphysical radar retrievals with in situ aircraft measurements. The MRMS is an attractive system for these tasks because the system incorporates data from multiple radars. Thus, the use of multiple radars allows for a much larger spatial range than provided by a single radar and can help mitigate potential calibration biases introduced by individual radars. The MRMS is also an operational product. Therefore, it would be possible for future MRMS upgrades to include real-time CONUS-wide microphysical retrieval estimates.
The MRMS consists of 143 WSR-88Ds as well as 30 radars from Environment and Climate Change Canada seamlessly interpolated onto a common 3D CONUS grid, forming a mosaic. Each individual WSR-88D typically has 0.5° (azimuthal) by 250-m (range) resolution with ranges of 420 km at each site. However, the MRMS grid itself has a horizontal resolution of 0.01° latitude × 0.01° longitude that corresponds to about a 1.1-km resolution and 33 irregularly spaced vertical levels with 250-m resolution from 0 to 3 km MSL, 500-m resolution from 3 to 9 km MSL, and 1000-m resolution from 9 to 20 km MSL (Smith et al. 2016). This grid in general has a temporal resolution of 2 min for Z and 10 min for each polarimetric variable. However, to maintain consistent radar retrievals (which use both Z and each polarimetric variable), we use the 10-min MRMS data. The MRMS radar merging methodology comes from the work of Lakshmanan et al. (2006) who uses exponential weights in time and in distance between the individual radar data bin and the MRMS grid cell. Radar data are averaged only between radars, whereas MRMS grid points near each radar are not smoothed. The Kdp is calculated by using a local linear fitting to differential phase along the radial direction (Ryzhkov and Zrnić 1996).
MRMS uncertainties for spatiotemporal collocation with aircraft can result from a number of factors. For example, individual radars can be miscalibrated, which can affect Z and ZDR, and Kdp can be noisy. MRMS uncertainties can also result from the multiradar merging algorithm (Lakshmanan et al. 2006), which preferentially weights data according to time and distance. Therefore, while there are many sources of potential error from the MRMS, we develop a methodology for incorporating some of these errors in section 4a when we perform our radar retrievals.
3. Case description
In situ and radar observations used in this work are from the 7 February 2020 IMPACTS case. Measurements were made across the final frontal system in a series of shortwaves that moved across the Ohio Valley and mid-Atlantic region in the preceding days. The fast-moving shortwave associated with the surface low pressure system was embedded in an extremely deep and negatively tilted trough centered over the Mississippi Valley with an attendant jet streak at 300 hPa off the mid-Atlantic coast that exceeded 225 kt (116 m s−1) by 0000 UTC 8 February 2020.
At the surface, a low pressure system formed in the Tennessee Valley at 1200 UTC 6 February 2020 and moved quickly northeast, undergoing rapid cyclogenesis on 7 February 2020 and deepening to 966 hPa by 0000 UTC 8 February 2020. Concurrently, an Arctic front moved southward on the backside of the system, resulting in frontogenetic forcing within the deformation zone over the Hudson Valley. A broad swath of heavy snow fell across the interior Northeast and New England, with up to 30 cm of snow falling across parts of central New York and snowfall rates exceeding 3–5 cm h−1.
The NASA P3 aircraft flew from 1400 to 1900 UTC, making multiple passes cutting across the cold front situated in upstate New York (Fig. 1). Another aircraft, the NASA ER2, made coordinated passes with the P3 for some of these flight legs as well. The ER2 aircraft flew above the P3 aircraft at approximately 20 km MSL and utilized multiple radar and lidar instruments. In particular, the ER2 aircraft housed a nadir-pointing X-band radar (EXRAD) as well as Ku-band and Ka-band radars (HIWRAP) that could be used to collocate radar reflectivity with the P3 aircraft location.
(a) 5.0-km Z at 1510 UTC 7 Feb 2020 from the MRMS system at 5.0 km MSL. The orange, dashed line corresponds to a convergence line estimate adapted from the 1500 UTC 7 Feb 2020 Weather Prediction Center surface analysis. The P3 flight track is shown in black, with the flight direction of each of the relevant flight legs denoted with arrows. (b) As in (a), but with the GOES-East visible satellite. (c) As in (a) and (b), but with just the terrain overlay. Also shown are time series of (d) altitude and (e) temperature (black; °C) with relative humidity with respect to ice (red; %) measured by the P3 probes. The periods corresponding to each leg are additionally shown as white triangles in (a)–(c).
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
The west side of the frontal boundary was composed of snow at the surface while the east side of the frontal boundary was composed of rain at the surface. Herein, we focus on performing radar and in situ retrievals along the first two P3 legs, which correspond to time periods when the P3 was to the west of the frontal boundary (Fig. 1). Leg 1 corresponds to roughly 1508–1524 UTC when the P3 flew west to east, whereas leg 2 corresponds to roughly 1530–1551 UTC when the P3 flew east to west. During leg 1, the P3 aircraft flew at approximately 5 km MSL (T ≈ −16°C) throughout the flight leg, whereas in leg 2 the P3 aircraft flew at approximately 3.5 km MSL but with temperatures decreasing from approximately −6°C at 1530 UTC to −13°C at 1551 UTC (Figs. 1d,e).
Figure 2 shows a Cartesian representation of the 5-km MRMS horizontal cross section for Z, ZDR, and Kdp at 1510 and 1520 UTC during leg 1 along with the P3 locations for these times. At both time periods, there are two separate Kdp and ZDR enhancement regions. At each period, the P3 aircraft is in the vicinity of Kdp values over 0.2° km−1 and ZDR values over 1.0 dB. In comparison with the stratiform nature of Z, these enhancement regions are very small [O(10 km)]. Despite these features occurring on such a small scale, the P3 aircraft flew through both regions nearly 100 km apart. This implies that it is possible for us to utilize in situ measurements taken aboard the aircraft in conjunction with the MRMS to directly investigate the microphysical properties of these regions (and in between) and to connect these dual-polarization radar features with physical, causal explanations. It is important to realize that the small scale of these features shown in Fig. 2 and the even smaller scale sampling volumes swept out by the in situ measurement probes provide a host of uncertainties in such an investigation. As a result, the next section lays out a general methodology for collocating the MRMS with the P3 in situ measurements that encompasses these spatiotemporal uncertainties.
MRMS horizontal cross sections of (top) Z, (middle) ZDR, and (bottom) Kdp at 5 km for (left) 1510 and (right) 1520 UTC during leg 1. Purple dots represent the P3 aircraft locations at each time period. The x and y coordinates are centered on the P3 aircraft location at 1510 UTC. MRMS radar artifacts near the center of a nearby radar can be seen in Z in the top-right corner.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
4. Methodology
a. In situ aircraft and MRMS radar collocation
MRMS horizontal Cartesian cross sections at 5.0 km MSL with P3 locations at MRMS time stamps (a) 1510 and (b) 1520 UTC (purple dots) and a time period in between MRMS time stamps (1513 UTC; black dots). Red lines indicate the temporal contribution to the weighting function, where shorter distances indicate a higher probability for the MRMS time period to be sampled. Also shown are 5.0 km MSL MRMS cross sections (zoomed in) at (left) 1510 UTC and (right) 1520 UTC for (c),(d) reflectivity; (e),(f) ZDR; and (g),(h) Kdp. MRMS points within the gray circles represent those of the nearest 1.0 km that are selected during the bootstrap sampling.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
We then use each sampled set of MRMS radar variables to retrieve three key bulk ice microphysical quantities: mean volume diameter Dm, number concentration Nt, and IWC. The Dm for this work is defined as the ratio of the fourth and third moments of the equivolume size distribution.
b. In situ probes and data processing methodologies
The P3 aircraft housed several different cloud probe instruments including a two-dimensional stereo (2DS) that has both horizontal (2DS-H) and vertical (2DS-V) oriented photodiode arrays, two high-volume precipitation spectrometers (HVPS-A and HVPS-B), a cloud probe imager (CPI), a fast-cloud droplet probe (FCDP), and a Rosemount icing detector (RICE), as well as several other probes not used in this study (Delene 2022). Aircraft probe data were processed using the Airborne Data Processing and Analysis package (Delene 2011). The FCDP and RICE instruments were used to indicate regimes of supercooled liquid drops and the CPI was used along with 3V-CPIview software (SPEC, Inc.) to estimate habit fractions for each leg. While the liquid water content (LWC) estimated from the FCDP was artificially high (up to 200%) due to shattering effects, we use the spikes in FCDP LWC as a way to infer the presence of supercooled liquid droplets. The University of Illinois/Oklahoma OAP Processing Software (UIOOPS; McFarquhar et al. 2018) was used to derive PSDs from the 2DS-V and HVPS-A probes. All 1-Hz PSD samples from the 2DS were ignored when dead time exceeded 70% (i.e., when the 2DS sampled less than 30% of the 1-Hz period). Stricter dead time thresholds would limit the number of time periods with higher concentrations. The remaining samples from the 2DS and HVPS were then averaged in 5-s time intervals for particle sizes in the range 0.1–30 mm where sizes smaller than 1.4 mm used 2DS data, and sizes larger than 1.4 mm used HVPS data. Particle sizes were defined using the circle of minimum enclosure approach Dmax of Wu and McFarquhar (2016).
Gamma and exponential distribution fits were performed on the P3 in situ observations for each 5-s binned PSD using the form
We estimated ice particle aspect ratios by fitting ellipses to each particle image. Ellipse fits and corresponding ellipse orientation angles were computed using the Computational Geometry Algorithms Library (CGAL) C++ library functions and the methodology of Welzl (1991) and Gärtner and Schönherr (1998) to find the ellipse of the smallest area enclosing a finite set of points. Ellipse orientation angles are defined as the angle between the maximum dimension and the time dimension (i.e., flight direction). Aspect ratio statistics were only computed for particle images that were fully within the probe sampling volume. Ellipse fits were computed from the horizontally and vertically oriented HVPS and 2DS datasets. While radar retrieval calculations used the actual intrinsic (spheroid) aspect ratios, these in situ ellipse fits represent 2D projections of each ice particle. Therefore, the projected ellipse fit aspect ratio Φ can be in the range φ ≤ Φ ≤ 1.0 depending upon particle orientations as imaged by the in situ cloud probes.
IWC was estimated from the observed PSDs following two different methods. The first method adopted a modified form of the static m–D relationship from Brown and Francis (1995, hereinafter BF95) as described in Hogan et al. (2012). The second method employed scattering simulations of varying riming degrees and m–D relations from Leinonen and Szyrmer (2015, hereinafter LS15) to determine the m–D relation that best minimized the reflectivity difference between ER2 HIWRAP radar measurements at the P3 location and that simulated from the PSDs for each 5-s interval. Unless otherwise noted, the majority of results use the BF95 IWC calculations since LS15 calculations are only available when the ER2 and P3 were less than 4 km apart horizontally and when data were within 3 min of each other.
c. Radar retrieval methodology
Radar retrieval equations used in this work come primarily from two sources: Ryzhkov and Zrnić (2019, hereinafter RZ19) and Bukovčić et al. (2020, hereinafter B20), with the exception of the nonpolarimetric IWC retrieval from Hogan et al. (2006, hereinafter H06) that uses Z (dBZ) and temperature in degrees Celsius. RZ19 and B20 assume that particles are described in terms of equivolume diameters D and are represented as oblate spheroids with fixed aspect ratios:
Table of retrieval equations used in this paper: Dm is in units of millimeters, Nt is in units of number per meter cubed, and IWC is in units of grams per meter cubed; Fs is the combined shape/orientation factor that depends upon
Figure 4 shows contours of Dm, Nt, and IWC relative uncertainties for various ranges of each independent variable in Eq. (3). Figure 4a shows that RZ19 retrievals of Dm depend equally on relative uncertainties in both ZDR and Kdp with relative errors strongly dependent upon the absolute Z uncertainty and weakly dependent on the magnitude of ZDR. Overall Dm uncertainties for σZ ≤ 5.0 dB are less than unity as long as relative uncertainties of ZDR and Kdp are also less than unity. The Nt uncertainties (Fig. 4b) mirror the contours of Dm uncertainties but are 2–4 times as large. RZ19 IWC retrievals (Fig. 4c) also depend equally on ZDR and Kdp relative errors but do not depend upon Z uncertainties. If ZDR is small (i.e., 0.1 dB), the relative errors of IWC depend equally upon the relative errors in ZDR and Kdp. Greater values of ZDR stretch IWC errors along the
Relative error contours for Dm (increments of 0.05), Nt (increments of 0.2), and IWC (increments of 0.1) retrievals using Eq. (3) and Table 1 for different values of
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
Although the apparent maximum dimensions and aspect ratios from in situ observations depend upon particle orientations (McFarquhar et al. 1999; Jiang et al. 2017; Dunnavan and Jiang 2019) that can be influenced by the cloud probe itself (King 1985, 1986), the use of maximum dimension here acts as a simple way to compare radar retrievals with in situ observations. It is important to note that, neglecting any optical errors, truly spheroidal particles will have either their semimajor (oblate) or semiminor (prolate) axis lengths represented in their 2D ellipse projections due to axis symmetry. Therefore, we would expect that only prolate-like, or triaxial particles would exhibit errors in Dmax, whereas oblate-like or planar crystals would have apparent maximum dimensions that differ strictly only to morphological effects (e.g., asymmetric branch structure).
To estimate
Notice from Table 1 that the general retrieval equations for B20 depend directly upon μ, α and Fs (which itself depends on
5. Results
a. In situ measurements and MRMS collocation
Time–height profiles of Z, ZDR, Kdp, and CDRp are shown in Fig. 5 for legs 1 and 2 using the bootstrapping procedure described in section 4c. Figure 5 shows a classic dendritic growth layer enhancement region in ZDR, Kdp, and CDRp from 1508 to 1512 UTC (with a maximum at approximately 1510:30 UTC). Values of Kdp in particular are quite high and approach values of 0.5° km−1 along the P3 flight path. A weaker Kdp, but higher ZDR enhancement region can be seen later from approximately 1520 UTC, which lies just above the frontal boundary on the easternmost portion of leg 1 (Fig. 1). During leg 2, the P3 aircraft flew at a lower altitude (higher temperature). As such, enhanced aggregation at this lower level contributed to larger Z and diminished Kdp and ZDR relative to leg 1.
A GOES-16 multichannel imaging algorithm (MIA; Hu et al. 2021) and CPI particle images suggests that the first Kdp and ZDR peak at 1510 UTC was primarily dominated by irregular dry ice particles, whereas the second peak at 1520 UTC corresponds to a large number of small, supercooled liquid droplets mixed together with ice particles. In Fig. 6, “irregular” particles (i.e., polycrystalline) constitute approximately 60% of all particles throughout each leg except during regions close to the frontal boundary where “spheroidal” particles, indicative of supercooled liquid droplets, occasionally dominate the observed particle type in conjunction with peaks in FCDP LWC measurements. While the 3V-CPIview software used to perform these habit classifications is not perfect because of misclassification of planar crystals and dendrites (such particles would normally be classified in the “plates” category but were instead misclassified as “rosettes”), casual observation of CPI images indicate that irregular polycrystals and spatial planar crystals do dominate the majority of ice throughout both legs. In addition, the RICE frequency behavior during this leg shows sharp frequency decreases (indicating the presence of liquid particles) at various times. Despite the P3 flying below cloud top, both MIA and RICE observations (Fig. 6) show a sharp transition from a dry ice environment to a mixed-phase environment around 1520 UTC. The correspondence between MIA and P3 observations appears to even capture the sharp spikes of LWC that are estimated from the FCDP data at 1521 UTC in Fig. 6. In particular, FCDP size distributions for 1521 UTC illustrate a pronounced peak at 25 μm, which is indicative of liquid droplets (Cober et al. 2001; McFarquhar et al. 2007). The appearance of this peak corresponds to a FCDP number concentration of 19 cm−3 (not shown). At 1510 UTC, on the other hand, the FCDP measured number concentration values that were fewer by a factor of approximately 6 (3 cm−3) and there was no pronounced peak. Therefore, the Kdp and ZDR peaks shown at 1510 UTC and 1520 UTC in Fig. 5 suggest two different microphysical regimes where 1510 UTC represents a population of dry ice, and 1520 UTC represents a population of supercooled liquid droplets and ice particles.
(top four rows) MRMS data collocated along the P3 aircraft flight path latitude and longitude (thick, black line) for (left) leg 1 and (right) leg 2. Also shown are (bottom) the number of unique MRMS grid points used for collocation with P3 aircraft.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
(top) CPI-derived habit fractions for (left) leg 1 and (right) leg 2. Plates were not identified for either leg and thus were not included in the plot. (bottom) Time series phase fractions of ice (blue), supercooled liquid droplets (red), and liquid water (yellow) as estimated by MIA for (left) leg 1 and (right) leg 2. Also shown are P3 estimates of LWC (gray) from the FCDP instrument and the RICE (black) frequency.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
Figure 7 shows the time series of bulk parameters Dm,max and Nt from the P3 observations and the exponential and gamma PSD fits. These plots not only show the observed values of Dm,max and Nt that will be used to evaluate MRMS retrieval performance, but they also show how well the exponential PSD and gamma PSD forms can capture the measured PSDs. Exponential PSD fits (μ = 0) tend to overestimate Dm,max by at least 0.3 mm and underestimate Nt by as much as 20 L−1. Therefore, errors in the exponential fitting alone can account for at least 50% of the observed Nt values. Using a gamma distribution (with variable μ) dramatically improves the PSD fit where Nt can be captured to within a few numbers per liter for a majority of each flight leg and Dm,max can be captured well within 0.1 mm. The μ is generally positive for a majority of leg 1 and for the last part of leg 2. Negative μ values are apparent during periods of supercooled liquid droplets in leg 1 (e.g., 1521 UTC) and leg 2 before 1538 UTC.
PSD measurements of Dm,max and Nt computed directly from measured PSDs from the 2DS and HVPS (black lines), gamma distribution fits (blue lines), and exponential distribution fits to the 2DS/HVPS PSDs (orange lines) for (top) μ, (top middle) Dm,max, (middle) Dm,max bias, (bottom middle) Nt, and (bottom) Nt bias for (left) leg 1 and (right) leg 2; Dm,max and Nt biases represent P3 probe observations subtracted from fits to the P3 probe observations. Fit parameters Dm,max and Nt were calculated over the 2DS and HVPS size range (0.10 ≤ Dmax ≤ 30.0 mm) representing the P3-observed PSD.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
The ZDR and Kdp during leg 1 are positively correlated with Nt for the dry ice period (Fig. 8) where the Kdp and ZDR peak at 1510 UTC corresponds to a Nt peak, whereas the later Kdp and ZDR peak at 1520 UTC does not show this correlation. From 1508 to 1518 UTC, correlation coefficients between Nt and the MRMS radar variables are 0.77 and 0.81 for Kdp and ZDR, respectively. While leg-1 reflectivities are generally higher during the first peak at 1510 UTC than during the second peak at 1520 UTC, the correlation between Nt and Z (0.48) is not as strong as with Kdp or ZDR from 1508 to 1518 UTC. The correlation coefficients for the entirety of leg 1 are overall lower for Kdp, ZDR, and Z (0.48, 0.16, and 0.45, respectively), whereas leg-2 correlations from 1530 to 1551 UTC are lower as well (0.43, 0.24, and −0.11, respectively). Size distributions for the 1510 UTC period (not shown) illustrate that the maximum dimensions of the particles corresponding to the Nt peak are between 0.1 and 1.0 mm, with a mean maximum dimension of approximately 0.5 mm.
(left) Leg-1 and (right) leg-2 time series of (top) ZDR, (middle) Kdp, and (bottom) Z along with the P3 probe Nt calculations from Fig. 7. Shaded regions represent the 5th–95th-percentile range of MRMS bootstrap samples.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
Ellipse fit joint distributions of major dimension and aspect ratio (Fig. 9) suggest that particles are apparently very spherical at all sizes for both the 1510 UTC period and the entirety of leg 1. While dendrites or other anisotropic particles could yield high ellipse fit aspect ratios if they are preferentially viewed along their faces, the observed small sizes instead suggest that the numerous particles responsible for the first Kdp and ZDR enhancement region between 1510 and 1512 UTC are either hexagonal plates or irregular particles consistent with the hypothesis of Moisseev et al. (2015), Griffin et al. (2018) and others or secondary ice production (e.g., Andrić et al. 2013). CPI images and the CPI classification algorithm shown in Fig. 6 suggests that there is no substantial change in habit type throughout the first Kdp and ZDR enhancement region. This is somewhat surprising considering that these enhancement regions could result strictly from changes in particle properties alone. Unlike in Griffin et al. (2018), the ZDR and Kdp during leg 1 are generally correlated, which suggests perhaps that both signatures were the result of the same or similar particle types. However, ρhv (not shown) indicates a brief, sharp decrease from 1510:10 to 1510:40 UTC and reaches a low value of almost 0.96. Thus, it is possible that some or all of the 1510 UTC Kdp and ZDR peak results from particles that are, for instance, denser than the observed irregular particles from Fig. 6. For example, this decrease in ρhv could result from secondary ice production such as from the collisional breakup of ice particles. The second Kdp and ZDR peak occurred at lower Nt and Z values (Fig. 8), which suggests that lower cloud tops with less numerous but larger anisotropic particles like dendrites were not masked as much by the presence of larger, more spherical particles (Griffin et al. 2018). This leads to ZDR being larger at 1520 UTC than during the 1510 UTC peak but with lower Kdp due to lower Nt at 1520 UTC. Particle images indicate local peaks in Dm,max shown in Fig. 7 are the result of large aggregates, irregular particles, and/or dendrites.
Observed leg-1 joint probability distributions of ellipse fit aspect ratios and maximum dimension for ice particles viewed from the horizontal plane: (a) 2DS probe for entire leg, (b) HVPS probe for entire leg, (c) 2DS probe for 1510–1512 UTC, and (d) HVPS probe for 1510–1512 UTC. Only particles that are fully within the probe sampling volumes are included.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
b. MRMS radar retrieval results
Figure 10 shows the leg-1 and leg-2 distributions of relative error for Dm,max, Nt, and IWC. Simplified retrievals (Fig. 10a) consistently underestimate Dm,max by approximately 50% for both legs. The Nt is consistently overestimated where RZ19 errors are less than B20. Both Nt and IWC relative error distributions are heavily skewed to larger, more positive relative errors, especially for leg 2 where retrievals are heavily biased during the 1530–1534 UTC period of supercooled liquid droplets (Fig. 6) where in situ measured IWC values are very low. The median relative errors of RZ19 Nt and IWC for leg 1 are small (+17% and −2%, respectively) but are much larger, and more positive, for leg 2 (+110% and +75%, respectively). B20 median relative errors are consistently larger overall for both Nt (+117%) and IWC (+43%)
Box-and-whisker plots of Dm,max, Nt, and IWC relative errors for legs 1 (1508–1524 UTC) and 2 (1530–1551 UTC) using the simplified retrieval equations from RZ19 and B20 (see Table 1). Circles represent the distribution mean, and horizontal black lines represent distribution medians. Boxes encompass the 25th–75th interquartile range, and whiskers represent the smallest and largest data values within 1.5 times the interquartile range.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
Retrieved aspect ratios from the Matrosov (2020) CDRp retrieval methodology (Fig. 10b) for each leg are very low (
Figure 11 compares the IWC MRMS retrievals from the two simplified polarimetric methods shown in Table 1 and the nonpolarimetric method (H06) with the BF95 and LS15 in situ measurements. H06 retrievals using the ER2 EXRAD reflectivity are also included for comparison as well. For leg 1, RZ19 retrievals are by far the most consistent with observations, whereas B20 overestimates IWC by approximately 50% and the H06 retrievals underestimate IWC by approximately 50%. Unlike the nonpolarimetric H06 retrievals, both the RZ19 and B20 simplified retrieval methodologies capture the correct trend where IWC fluctuates at higher values from 1508 to 1515 UTC and then decreases from 1515 to 1518 UTC. B20 is overall biased more to larger IWC values than RZ19 for both legs and in situ IWC estimates and thus produces larger overall errors. The H06 methodologies generally underestimate IWC relative to observations by approximately 50% for the majority of the leg-1 dry ice half but more consistently captures the observed low IWC values during the supercooled liquid droplet period around 1520 UTC. Using the ER2 EXRAD reflectivity does improve the IWC retrievals but only by about 10%. During leg 2, both RZ19 and B20 on average overestimates IWC by well over 50%, whereas both H06 retrievals on average captures IWC to within 50% of both sets of in situ measurements.
(top) Comparison of bootstrap IWC retrievals following RZ19 (green; simplified RZ19), B20 (red; simplified B20), and H06 using the MRMS bootstraps (blue) and EXRAD reflectivity (purple) with two PSD estimates obtained from the P3 aircraft for legs (left) 1 and (right) 2. (bottom) IWC relative error boxplots for each leg using both BF95 and LS15 IWC estimates. Boxplots have the same statistical meanings as in Fig. 10.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
B20 retrievals with optimized parameters
While the previous section evaluated the simplified approaches with assumed parameters for α and σ, it is natural to wonder how well retrievals could capture the measured Dm,max, Nt, and IWC if each assumed parameter varies. To investigate such optimal parameter sets, we use MCMC runs along with the general B20 equations from Table 1 such that all three parameters:
B20 parameters and corresponding retrievals for (from top to bottom): α,
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
Unlike α and
Correlations between each variable are shown in Fig. 13 where each set of joint posterior distributions corresponds to the collated optimal parameters for the μ-dependent B20 retrievals for each 5 s time period. As shown in this figure, each pair of parameters is roughly related to one another in the form of power laws. In particular, the correlations between
Combined posterior parameter distributions for generalized B20 retrievals using MCMC methodology for (top) leg 1 and (bottom) leg 2.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
6. Discussion
a. What caused each dendritic growth layer Kdp and ZDR enhancement?
The collocated MRMS variables with in situ measurements shown in section 5a provide a complex picture for inferring which physical mechanism caused each Kdp and ZDR enhancement in leg 1. Here, we attempt to explain which physical mechanisms could plausibly explain each set of observations for each enhancement region and how each mechanism relates to previous theories and hypotheses.
For the first Kdp and ZDR peak at 1510 UTC, Fig. 8 clearly suggests that each enhancement is in some way related to a local enhancement in Nt. However, it is not entirely clear which physical mechanism corresponds to such an Nt enhancement. A local increase in the number of dendrites or planar crystals, as suggested by studies such as Kennedy and Rutledge (2011), Schrom et al. (2015), and Schrom (2018), provides perhaps the most straightforward explanation. However, the 3V-CPIview habit classifications (Fig. 6) do not suggest a prevalence of such particles at any portion of the flight, and the majority of particles tend to rather be irregular. Secondary ice production processes such as rime splintering represent another possibility. However, the lack of supercooled liquid droplets and the colder temperatures during leg 1 do not suggest any rime splintering. Another secondary ice production mechanism, collisional breakup of ice particles, represents a better possible explanation. A laboratory experiment by Lo (1983) suggests that fragment sizes of ice–ice collisions are approximately exponentially distributed where a majority of fragments (for collisions of particles 3 mm or larger) have sizes less than 0.6 mm. These fragment sizes are somewhat consistent with the 0.5-mm particles that primarily account for the Nt enhancement at 1510 UTC. Another laboratory study by Vardiman (1978) illustrates that light-to-moderately rimed spatial planar crystals, such as those observed in the current study, were the most efficient at producing fragments. However, Vardiman (1978) also found that such breakup required heavily rimed particles or graupel to produce such fragments. Thus, it is unclear whether secondary ice from collisional breakup could even efficiently take place for the first Kdp and ZDR peak given the lack of observed rimed particles and supercooled liquid droplets. Some in situ studies such as Lo and Passarelli (1982) have, however, emphasized that such rimed particles are not required for breakup to occur. Theoretical studies such as Phillips et al. (2017) have highlighted the importance of ice–ice collisions yet some researchers such as Korolev and Leisner (2020) have questioned the reliability of such theory in the context of known artificial shattering issues induced by cloud probes and differences in laboratory setups for the few experiments that do exist. Considering the apparent correlation between the increase in Nt with increases in Kdp and ZDR for the first peak, we hypothesize that either ice–ice collisional breakup or a seeder–feeder process such as that proposed by Moisseev et al. (2015) most likely explains the polarimetric enhancements at 1510 UTC.
For the second Kdp and ZDR peak at about 1520 UTC, in situ probe data illustrates that supercooled liquid droplets were present (Fig. 6). Riming therefore provides a possible explanation for the second ZDR and Kdp peak. Riming would naturally increase the density of ice particles, which, in turn, would enhance the effect anisotropic particles have on Kdp and ZDR measurements. Particles at 1520 UTC were also generally larger than at 1510 UTC as indicated by Dm,max ≈ 3.0 mm values as compared with Dm,max ≈ 2.0 mm at 1510 UTC (Fig. 7). As indicated in the reflectivity time–height profile for leg 1 (Fig. 5), cloud heights were generally lower at 1520 UTC than at 1510 UTC. As suggested by Griffin et al. (2018) and briefly described in section 5a, the lower and warmer cloud top at 1520 UTC likely led to fewer large aggregates that usually tend to mute ZDR in the dendritic growth layer. Large anisotropic crystals (as observed in some of the 2DS images) within these shallower cloud systems could therefore contribute more to ZDR and Kdp than for colder clouds with higher cloud tops. This also explains why ZDR was larger whereas Kdp was smaller for the second peak than for the first. However, unlike in Griffin et al. (2018), ZDR and Kdp are clearly correlated, rather than anticorrelated. Also, unlike in Griffin et al. (2018), it does not appear that hexagonal plates were the primary cause of the second enhancement region. Thus, our best explanation for the cause of the second ZDR and Kdp enhancement region at 1520 UTC is based on the combination of 1) riming affecting particle density, 2) the presence of a few large, anisotropic particles, and 3) a lower cloud top that limits the number of large aggregates and, in turn, reduces the ZDR masking effect common for colder clouds.
b. What are the true intrinsic ice particle aspect ratios?
While MRMS retrievals seemingly require low aspect ratios (
To help reconcile this paradox, we perform idealized simulations of rotating and projecting oblate spheroids (mimicking low-aspect-ratio, planar crystals) and triaxial ellipsoids (mimicking various “irregular” particles) under various aspect ratio and orientation assumptions. This allows us to investigate which combinations of ice particle shapes and orientations could explain the in situ projection observations as shown in Fig. 9. The results of such an analysis are shown in Fig. 14 for oblate spheroids with aspect ratios of 0.2 (green, purple, and cyan lines), which are symbolic of
Leg 1 probability density comparison of observed ellipse projections from the P3 HVPS instrument (shaded histograms) and Monte Carlo simulated estimates of various ellipsoids (colored lines) assuming various orientation distributions for test conditions given in Table 2: (a) Horizontal view ellipse aspect ratio 2D projection distributions. (b) Vertical view ellipse aspect ratio 2D projection distributions. (c) Horizontal view ellipse orientation 2D projection distributions. (d) Vertical view ellipse orientation 2D projection distributions. Also shown is the (e) Monte Carlo simulated joint distribution (shaded) of triaxial ellipsoid aspect ratios used for projection analysis with the analytic representation (contours). The oblate spheroid aspect ratio in the projection analysis is represented with a green circle at φba = 1.0 and φcb = 0.2.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
Ellipsoid test parameters used in Figs. 14 and 15, where a ≥ b ≥ c are semiaxis lengths with aspect ratios: φba ≡ b/a, φca ≡ c/a, and φcb ≡ c/b; B(a, b) is the beta function; K(x) is the complete elliptic integral of the first kind; 0 ≤ ϕ(ψ) ≤ 2π; 0 ≤ θ ≤ π/2;
For φ = 0.2 oblate spheroids to yield the observed bell-shaped 2D projected aspect ratio distributions, the spheroids need to be canted such that their faces are preferentially pointed toward the imager as suggested by some of the 2DS particle images. As shown in Fig. 14, a Gaussian distribution for θ can approximate the aspect ratio distribution for one view (i.e., horizontal view) but not for the other view (purple lines). Normal distributions for all three Euler orientation angles (cyan lines) can produce similar bell-shaped aspect ratio distributions to that of the observations but cannot capture the concave 2D ellipse orientation distribution. As stated before, retrievals are thought to potentially introduce intrinsic aspect ratio uncertainties up to 0.2 (Matrosov 2020). Despite this uncertainty, the inability for oriented oblate spheroids to produce the correct distributions of projected ellipse aspect ratios and orientations suggests that the observed ice particle shapes were more representative of irregular shapes characteristic of the triaxial ellipsoid tests in Fig. 14. This is also consistent with CPI habit classifications (Fig. 6), which indicate that irregular particles, rather than hexagonal plates or dendrites, dominate ice habits for all periods.
This projection behavior for oblate spheroids is confirmed in Fig. 15, which shows a side-by-side comparison of example 2DS images of dendrites and planar particles with the 2DS-derived projected aspect ratio and ellipse orientation distributions for each orthogonal 2DS probe view. Unlike the HVPS, the 2DS-H and 2DS-V projections have partially overlapping sampling volumes and thus the 2DS-H and 2DS-V distributions more-or-less represent the same particles being viewed along flight (2DS-H) and along the particle fall speed direction (2DS-V). Moreover, an additional Monte Carlo oblate spheroid projection test (test 7 in Table 2) is performed to compare both the resulting distributions and the example projection images. Here, oblate spheroid aspect ratios are randomly chosen between 0.1 and 0.3 to mimic potential observed dendrites and planar crystals consistent with retrievals using CDRp. Similar to Fig. 14, low-aspect-ratio oblate spheroids can explain the observed 2D ellipse orientation distributions if they exhibit some beta-like canting angle distribution with their major dimension randomly oriented in the horizontal plane. As in Fig. 14, the resulting projected aspect ratio distribution peaks at lower values than observed. Example 2DS images of dendrites and planar crystals in Fig. 15 are qualitatively consistent with the Monte Carlo simulated oblate spheroids: both viewing directions show examples where the basal facet of the planar crystals are tilted away from the imager as well as facing toward the imager. This further illustrates that a mixture of a low number concentration of large dendrites and a higher concentration of higher-aspect-ratio, irregular, and aggregate particles with triaxial shapes can explain the 2DS and HVPS observations.
(left) 2DS images of dendrites and planar particles taken during leg 1. Clear examples of dendrites and planar particles are circled in red. (right) Comparison of 2DS distributions of 2D ellipse fit aspect ratios and ellipse orientation angles (shaded histograms) and analogous distributions of beta-oriented (see Table 2) Monte Carlo simulated 3D oblate spheroids with random aspect ratios between 0.1 and 0.3 (lines). Below each distribution plot are 5 random samples of oblate spheroid views assuming the (a) 2DS-H and (b) 2DS-V orthogonal viewing angles. The 2DS-H and 2DS-V share the same sampling volume.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
c. MRMS retrieval uncertainties and biases
One reason for a possible underestimation of radar-based retrieved particle aspect ratios are potential biases in the estimates of CDRp, which are based on averaged values of Zdr and ρhv [Eq. (5)] as reported by the MRMS system. Previous studies (Matrosov et al. 2017; Matrosov 2020) with radar measurements at closer ranges [thus generally with high signal-to-noise ratios (SNR)] indicated that minimal CDRp values corresponding to quasi-spherical particles with larger aspect ratios are around −27 dB. High minimal CDRp values are generally expected from particles with lower aspect ratios (Matrosov et al. 2020). CDRp estimates are quite sensitive to biases in Zdr and ρhv data. Underestimation of ρhv values at low SNRs by about 0.006 could explain an approximate −4 dB in minimal CDRp values. Additionally, radar-based ice particle aspect ratio retrievals are weighted by larger particles because of a general Dn dependence of radar variables as shown in Eq. (A9). As a result of this dependence, particle mixtures (e.g., more numerous smaller quasi-spherical particles and larger dendritic particles) will yield biased retrieved aspect ratios relative to a simple mean aspect ratio value.
The MRMS spatiotemporal resolution could also explain some of the observed discrepancies between the in situ observed Dm,max peaks as shown in Fig. 7 and the retrieved Dm,max from the results section (section 5). Because the Dm,max peaks are rather sharp, it is possible for the MRMS collocated time series to not properly resolve these regions of enhanced Z responsible for these peaks. Figure 16 shows evidence for this hypothesis by comparing the collocated MRMS reflectivity and the airborne ER2 EXRAD X-band reflectivity with the P3-observed Dm,max. Figure 16 shows that various Dm,max peaks (e.g., 1515, 1519, 1537, and 1544 UTC) correlate with the ER2 EXRAD X-band radar reflectivity peaks, whereas the MRMS S-band radar bootstrap observations are not able to resolve these peaks. Furthermore, the MRMS Z is generally biased toward lower values than the EXRAD observations. This suggests that some of the biases in the retrievals (i.e., low Dm,max and high Nt) could result from a low bias in Z from the MRMS. The 2-min MRMS Z (blue) are very similar to the 10-min MRMS Z except from 1514 to 1521 UTC and from 1534 to 1538 UTC where the 2-min data has slightly lower reflectivity values. This suggests that the majority of the discrepancy between the EXRAD and MRMS comes from spatial uncertainties rather than temporal uncertainties.
Time series of reflectivity for the ER2 EXRAD X-band radar (black lines) and the 2-min (blue) and 10-min (gray) MRMS S-band radar along with estimates of Dm,max from the P3 aircraft. For the EXRAD reflectivity, time periods are not included if the horizontal difference in distance between the ER2 and P3 is greater than 4 km.
Citation: Journal of Applied Meteorology and Climatology 61, 11; 10.1175/JAMC-D-21-0220.1
7. Summary and conclusions
This study shows that current PSD and aspect ratio dual-polarization radar retrieval methodologies for ice particles in a winter storm can produce similar values and behaviors to those derived explicitly from in situ observations in the mean. However, there are several important limitations, uncertainties, and caveats that heavily influence these retrieved PSD and aspect ratio variables. Using simplified retrieval equations that assume typical values of
Some discrepancy between the retrieved particle aspect ratios and their in situ estimates necessarily resulted from discrepancies between the radar retrievals that represent intrinsic or 3D mean aspect ratio values and in situ ellipse fit data that were derived from 2D projections. Thus, some planar particle projections will have artificially large 2D ellipse fit aspect ratios because of their quasi-random orientations within the 2DS and HVPS sample volumes. However, numerical simulations (Figs. 14 and 15) showed that low-aspect-ratio planar particles cannot simultaneously match the observed distribution of 2D ellipse orientation angles and aspect ratios for both horizontally oriented and vertically oriented in situ cloud probes. To match the observed distribution of 2D ellipse angles for both probe orientations, planar particles must be quasi-randomly oriented as might be expected from turbulence generated within and around in situ cloud probes. However, these orientations necessarily lead to a bimodal distribution of 2D ellipse aspect ratios rather than the observed unimodal distribution. More likely, MRMS uncertainties and retrieval assumptions are responsible for the observed discrepancies. In particular, a comparison between the collocated MRMS reflectivity and that of the ER2 EXRAD radar and measurements of Dm,max taken from the P3 probes (Fig. 16) illustrates that the MRMS cannot resolve local Z peaks responsible for enhanced Dm,max values due to the MRMS horizontal resolution. Another source of discrepancy between retrievals of aspect ratios and their 2D projections is that radar-based retrievals are weighted by larger dendritic-type particles due to reflectivity dependence as approximately D4, while in situ estimates represent simple mean values, which could be dominated by smaller quasi-spherical particles.
The ability to collocate polarimetric radar variables directly with in situ aircraft data has also allowed us to investigate the various underlying causes of the often-discussed dendritic growth layer Kdp and ZDR enhancement regions associated with the onset of aggregation (Kennedy and Rutledge 2011; Moisseev et al. 2015). In situ aircraft and satellite observations (Fig. 6) illustrate that both dry ice particles and mixtures of ice particles and supercooled liquid droplets can exist in these enhancement regions. The first Kdp and ZDR peak occurred in a region with little to no observed supercooled liquid droplets and no noticeable change in ice particle habit (Fig. 6). Low temperatures (T ≈ −16°C) and lack of supercooled liquid droplets suggests that rime splintering was not active in this region. However, this first peak did show elevated values of Nt, which suggests that increases in Nt from irregular crystals drove enhancements in both Kdp and ZDR. It is possible that secondary ice production such as ice–ice collisional breakup could plausibly explain this local increase in Nt. However, various studies present conflicting views on the need for rimed particles for such breakup to occur as well as the reliability of previous experimentation for understanding the ice–ice collisional breakup process (Lo and Passarelli 1982; Vardiman 1978; Phillips et al. 2017; Korolev and Leisner 2020). On the other hand, in situ probe data illustrate that supercooled liquid droplets were present during the second Kdp and ZDR peak. Higher Dm,max and lower Nt values at 1520 UTC than at 1510 UTC suggest that large, anisotropic particles at 1520 UTC had a larger impact on Kdp and ZDR than the smaller, irregular particles. Overall, despite similar polarimetric radar variable values and tendencies, in situ observations in these regions suggest that IWC, Nt, and Dm,max in these regions can be vastly different. It is clear that some of these differences are attributable to differences in particle characteristics related to particle shape (
Herein, radar variable subscripts represent either linear units (lowercase) or logarithmic units (uppercase).
Acknowledgments.
The authors are very thankful for the NSSL MRMS team and, in particular Andrew Osborne, who graciously provided us with the MRMS datasets used for this study. We are also tremendously grateful to the IMPACTS team that collected a comprehensive in situ and remote sensing dataset enabling this study. In particular, we thank Gerry Heymsfield for his input into this study with regard to the ER2 EXRAD data. The authors also thank Peter Brechner (University of Oklahoma School of Meteorology) for his willingness to provide the incomplete gamma fitting routine used in this work as well as Paul Lawson and Sarah Woods (SPEC, Inc.) for providing us with the 3V-CPIview software we used for CPI habit classification. The authors thank the three reviewers of this paper who helped to greatly improve this study with their feedback and suggestions. The IMPACTS project was funded by the NASA Earth Venture Suborbital-3 (EVS-3) program and managed by the Earth System Science Pathfinder (ESSP) Program Office. Funding for this study was provided by NOAA/Office of Oceanic and Atmospheric Research under NOAA–University of Oklahoma Cooperative Agreement NA21OAR4320204, U.S. Department of Commerce, by the U.S. Department of Energy’s Atmospheric Science Program Atmospheric System Research, an Office of Science, Office of Biological and Environmental Research program under DOE ASR Grant DE-SC0018967, by NASA Grant 80NSSC19K0399, and by NSF Grant 1841246.
Data availability statement.
Archival MRMS data used in this paper were requested directly from the NSSL MRMS team. The IMPACTS EXRAD, CPI, RICE, FCDP, 2DS, and HVPS probes used within this analysis are available online from the NASA EOSDIS Global Hydrology Resource Center Distributed Activate Archive Center (https://doi.org/10.5067/IMPACTS/DATA101). UND particle probe data can be found online (https://doi.org/10.5067/IMPACTS/MULTIPLE/DATA101; Delene 2022). Aircraft data were processed using the Airborne Data Processing and Analysis package (Delene 2011), the University of Illinois/Oklahoma Optical Array Probe (OAP) Processing Software (https://doi.org/10.5281/zenodo.1285969), and 3V-CPIview software (SPEC, Inc.; http://www.specinc.com/3v-cpi-combo).
APPENDIX
Radar Retrieval Methodology
Notice in Eq. (A1) that ρs becomes unphysically larger than the bulk density of ice for small sizes. As a result, the CDRp retrievals described in section 4c explicitly sets the upper limit of ρs to the bulk density of ice ρi for these small sizes. For the RZ19 and B20 retrieval methodologies, this inconsistency is ignored such that the retrievals can be expressed as simple formulas. Therefore, the following equations are specifically valid for sufficiently large Dm (Dm > 1.0 mm).
REFERENCES
Andrić, J., M. R. Kumjian, D. S. Zrnić, J. M. Straka, and V. M. Melnikov, 2013: Polarimetric signatures above the melting layer in winter storms: An observational and modeling study. J. Appl. Meteor. Climatol., 52, 682–700, https://doi.org/10.1175/JAMC-D-12-028.1.
Auer, A. H., and D. L. Veal, 1970: The dimension of ice crystals in natural clouds. J. Atmos. Sci., 27, 919–926, https://doi.org/10.1175/1520-0469(1970)027<0919:TDOICI>2.0.CO;2.
Bailey, M. P., and J. Hallett, 2009: A comprehensive habit diagram for atmospheric ice crystals: Confirmation from the laboratory, AIRS II, and other field studies. J. Atmos. Sci., 66, 2888–2899, https://doi.org/10.1175/2009JAS2883.1.
Bechini, R., L. Baldini, and V. Chandrasekar, 2013: Polarimetric radar observations in the ice region of precipitating clouds at C-band and X-band radar frequencies. J. Appl. Meteor. Climatol., 52, 1147–1169, https://doi.org/10.1175/JAMC-D-12-055.1.
Brandes, E. A., K. Ikeda, G. Zhang, M. Schönhuber, and R. M. Rasmussen, 2007: A statistical and physical description of hydrometeor distributions in Colorado snowstorms using a video disdrometer. J. Appl. Meteor. Climatol., 46, 634–650, https://doi.org/10.1175/JAM2489.1.
Brown, P. R. A., and P. N. Francis, 1995: Improved measurements of the ice water content in cirrus using a total-water probe. J. Atmos. Oceanic Technol., 12, 410–414, https://doi.org/10.1175/1520-0426(1995)012<0410:IMOTIW>2.0.CO;2.
Bukovčić, P., A. Ryzhkov, D. Zrnić, and G. Zhang, 2018: Polarimetric radar relations for quantification of snow based on drisdrometer data. J. Appl. Meteor. Climatol., 57, 103–120, https://doi.org/10.1175/JAMC-D-17-0090.1.
Bukovčić, P., A. Ryzhkov, and D. Zrnić, 2020: Polarimetric relations for snow estimation – radar verification. J. Appl. Meteor. Climatol., 59, 991–1009, https://doi.org/10.1175/JAMC-D-19-0140.1.
Chase, R. J., S. W. Nesbitt, and G. M. McFarquhar, 2021: A dual-frequency radar retrieval of two parameters of the snowfall particle size distribution using a neural network. J. Appl. Meteor. Climatol., 60, 341–359, https://doi.org/10.1175/JAMC-D-20-0177.1.
Chen, J.-P., and D. Lamb, 1994: The theoretical basis for the parameterization of ice crystal habit: Growth by vapor deposition. J. Atmos. Sci., 51, 1206–1222, https://doi.org/10.1175/1520-0469(1994)051<1206:TTBFTP>2.0.CO;2.
Cober, S. G., G. A. Isaac, A. V. Korolev, and J. W. Strapp, 2001: Assessing cloud-phase conditions. J. Appl. Meteor., 40, 1967–1983, https://doi.org/10.1175/1520-0450(2001)040<1967:ACPC>2.0.CO;2.
Delene, D. J., 2011: Airborne data processing and analysis software package. Earth Sci. Inform., 4, 29–44, https://doi.org/10.1007/s12145-010-0061-4.
Delene, D. J., 2022: UND Cloud Microphysics IMPACTS 2020, version 1. NASA Global Hydrometeorology Resource Center DAAC, accessed 21 February 2021, https://doi.org/10.5067/IMPACTS/MULTIPLE/DATA101.
Dunnavan, E. L., and Z. Jiang, 2019: A general method for estimating bulk 2d projections of ice particle shape: Theory and applications. J. Atmos. Sci., 76, 305–332, https://doi.org/10.1175/JAS-D-18-0177.1.
Dunnavan, E. L., Z. Jiang, J. Y. Harrington, J. Verlinde, K. Fitch, and T. J. Garrett, 2019: The shape and density evolution of snow aggregates. J. Atmos. Sci., 76, 3919–3940, https://doi.org/10.1175/JAS-D-19-0066.1.
Finlon, J. A., G. M. McFarquhar, S. W. Nesbitt, R. M. Rauber, H. Morrison, W. Wu, and P. Zhang, 2019: A novel approach for characterizing the variability in mass-dimension relationships: Results from MC3E. Atmos. Chem. Phys., 19, 3621–3643, https://doi.org/10.5194/acp-19-3621-2019.
Fitch, K. E., C. Hang, A. Talaei, and T. J. Garrett, 2021: Arctic observations and numerical simulations of surface wind effects on Multi-Angle Snowflake Camera measurements. Atmos. Meas. Tech., 14, 1127–1142, https://doi.org/10.5194/amt-14-1127-2021.
Garrett, T. J., S. E. Yuter, C. Fallgatter, K. Shukurko, S. R. Rhodes, and J. L. Endries, 2015: Orientations and aspect ratios of falling snow. Geophys. Res. Lett., 42, 4617–4622, https://doi.org/10.1002/2015GL064040.
Gärtner, B., and S. Schönherr, 1998: Exact primitives for smallest enclosing ellipses. Inf. Process. Lett., 68, 33–38, https://doi.org/10.1016/S0020-0190(98)00132-X.
Goodman, J., and J. Weare, 2010: Ensemble samplers with affine invariance. Comm. App. Math. Comp. Sci., 5, 65–80, https://doi.org/10.2140/camcos.2010.5.65.
Grazioli, J., G. Lloyd, L. Panziera, C. R. Hoyle, P. J. Connolly, J. Henneberger, and A. Berne, 2015: Polarimetric radar and in situ observations of riming and snowfall microphysics during CLACE 2014. Atmos. Chem. Phys., 15, 13 787–13 802, https://doi.org/10.5194/acp-15-13787-2015.
Griffin, E. M., T. J. Schuur, and A. V. Ryzhkov, 2018: A polarimetric analysis of ice microphysical processes in snow, using quasi-vertical profiles. J. Appl. Meteor. Climatol., 57, 31–50, https://doi.org/10.1175/JAMC-D-17-0033.1.
Grinsted, A., 2022: Ensemble MCMC sampler. Github, accessed 18 April 2022, https://github.com/grinsted/gwmcmc.
Harrington, J. Y., K. Sulia, and H. Morrison, 2013: A method for adaptive habit prediction in bulk microphysical models. Part II: Parcel model corroboration. J. Atmos. Sci., 70, 365–376, https://doi.org/10.1175/JAS-D-12-0152.1.
Heymsfield, A. J., A. Bansemer, C. Schmitt, C. Twohy, and M. R. Poellot, 2004: Effective ice particle densities derived from aircraft data. J. Atmos. Sci., 61, 982–1003, https://doi.org/10.1175/1520-0469(2004)061<0982:EIPDDF>2.0.CO;2.
Hogan, R. J., M. P. Mittermaier, and A. J. Illingworth, 2006: The retrieval of ice water content from radar reflectivity factor and temperature and its use in evaluating a mesoscale model. J. Appl. Meteor. Climatol., 45, 301–317, https://doi.org/10.1175/JAM2340.1.
Hogan, R. J., L. Tian, P. R. A. Brown, C. D. Westbrook, A. J. Heymsfield, and J. D. Eastment, 2012: Radar scattering from ice aggregates using the horizontally aligned oblate spheroid approximation. J. Appl. Meteor. Climatol., 51, 655–671, https://doi.org/10.1175/JAMC-D-11-074.1.
Hu, J., D. Rosenfeld, Y. Zhu, X. Lu, and J. Carlin, 2021: Multi-channel Imager Algorithm (MIA): A novel cloud-top phase classification algorithm. Atmos. Res., 261, 105767, https://doi.org/10.1016/j.atmosres.2021.105767.
Jensen, A. A., and J. Y. Harrington, 2015: Modeling ice crystal aspect ratio evolution during riming: A single-particle growth model. J. Atmos. Sci., 72, 2569–2590, https://doi.org/10.1175/JAS-D-14-0297.1.
Jiang, Z., M. Oue, J. Verlinde, E. E. Clothiaux, K. Aydin, G. Botta, and Y. Lu, 2017: What can we conclude about the real aspect ratios of ice particle aggregates from two-dimensional images? J. Appl. Meteor. Climatol., 56, 725–734, https://doi.org/10.1175/JAMC-D-16-0248.1.
Kennedy, P. C., and S. A. Rutledge, 2011: S-band dual-polarization radar observations of winter storms. J. Appl. Meteor. Climatol., 50, 844–858, https://doi.org/10.1175/2010JAMC2558.1.
King, W. D., 1985: Air flow and particle trajectories around aircraft fuselages. III: Extensions to particles of arbitrary shape. J. Atmos. Oceanic Technol., 2, 539–547, https://doi.org/10.1175/1520-0426(1985)002<0539:AFAPTA>2.0.CO;2.
King, W. D., 1986: Air flow and particle trajectories around aircraft fuselages. IV: Orientation of ice crystals. J. Atmos. Oceanic Technol., 3, 433–439, https://doi.org/10.1175/1520-0426(1986)003<0433:AFAPTA>2.0.CO;2.
Korolev, A., and G. Isaac, 2003: Roundness and aspect ratio of particles in ice clouds. J. Atmos. Sci., 60, 1795–1808, https://doi.org/10.1175/1520-0469(2003)060<1795:RAAROP>2.0.CO;2.
Korolev, A., and T. Leisner, 2020: Review of experimental studies of secondary ice production. Atmos. Chem. Phys., 20, 11 767–11 797, https://doi.org/10.5194/acp-20-11767-2020.
Kumjian, M. R., and K. A. Lombardo, 2017: Insights into the evolving microphysical and kinematic structure of northeastern U.S. winter storms from dual-polarization Doppler radar. Mon. Wea. Rev., 145, 1033–1061, https://doi.org/10.1175/MWR-D-15-0451.1.
Kumjian, M. R., S. Mishra, S. E. Giangrande, T. Toto, A. V. Ryzhkov, and A. Bansemer, 2016: Polarimetric radar and aircraft observations of saggy bright bands during MC3E. J. Geophys. Res. Atmos., 121, 3584–3607, https://doi.org/10.1002/2015JD024446.
Lakshmanan, V., T. Smith, K. Hondl, G. J. Stumpf, and A. Witt, 2006: A real-time, three dimensional, rapidly updating, heterogeneous radar merger technique for reflectivity, velocity, and derived products. Wea. Forecasting, 21, 802–823, https://doi.org/10.1175/WAF942.1.
Lamb, D., and W. D. Scott, 1972: Linear growth rates of ice crystals grown from the vapor phase. J. Cryst. Growth, 12, 21–31, https://doi.org/10.1016/0022-0248(72)90333-8.
Leinonen, J., and W. Szyrmer, 2015: Radar signature of snowflake riming: A modeling study. Earth Space Sci., 2, 346–358, https://doi.org/10.1002/2015EA000102.
Lo, K. K., 1983: Growth processes of snow. Ph.D. thesis, Massachusetts Institute of Technology, 193 pp.
Lo, K. K., and R. E. Passarelli, 1982: The growth of snow in winter storms: An airborne observational study. J. Atmos. Sci., 39, 697–706, https://doi.org/10.1175/1520-0469(1982)039<0697:TGOSIW>2.0.CO;2.
Locatelli, J. D., and P. V. Hobbs, 1974: Fall speeds and masses of solid precipitation particles. J. Geophys. Res., 79, 2185–2197, https://doi.org/10.1029/JC079i015p02185.
Magono, C., and T. Nakamura, 1965: Aerodynamic studies of falling snowflakes. J. Meteor. Soc. Japan, 43, 139–147, https://doi.org/10.2151/jmsj1965.43.3_139.
Matrosov, S. Y., 2020: Ice hydrometeor shape estimations using polarimetric operational and research radar measurements. Atmosphere, 11, 97, https://doi.org/10.3390/atmos11010097.
Matrosov, S. Y., C. G. Schmitt, M. Maahn, and G. de Boer, 2017: Atmospheric ice particle shape estimates from polarimetric radar measurements and in situ observations. J. Atmos. Oceanic Technol., 34, 2569–2587, https://doi.org/10.1175/JTECH-D-17-0111.1.
Matrosov, S. Y., A. V. Ryzhkov, M. Maahn, and