Search for: All records

Abstract On 7 February 2020, precipitation within the comma-head region of an extratropical cyclone was sampled remotely and in situ by two research aircraft, providing a vertical cross section of microphysical observations and fine-scale radar measurements. The sampled region was stratified vertically by distinct temperature layers and horizontally into a stratiform region on the west side, and a region of elevated convection on the east side. In the stratiform region, precipitation formed near cloud top as side-plane, polycrystalline, and platelike particles. These habits occurred through cloud depth, implying that the cloud-top region was the primary source of particles. Almost no supercooled water was present. The ice water content within the stratiform region showed an overall increase with depth between the aircraft flight levels, while the total number concentration slightly decreased, consistent with growth by vapor deposition and aggregation. In the convective region, new particle habits were observed within each temperature-defined layer along with detectable amounts of supercooled water, implying that ice particle formation occurred in several layers. Total number concentration decreased from cloud top to the −8°C level, consistent with particle aggregation. At temperatures > −8°C, ice particle concentrations in some regions increased to >100 L −1 , suggesting secondary ice production occurred at lower altitudes. WSR-88D reflectivity composites during the sampling period showed a weak, loosely organized banded feature. The band, evident on earlier flight legs, was consistent with enhanced vertical motion associated with frontogenesis, and at least partial melting of ice particles near the surface. A conceptual model of precipitation growth processes within the comma head is presented. Significance Statement Snowstorms over the northeast United States have major impacts on travel, power availability, and commerce. The processes by which snow forms in winter storms over this region are complex and their snowfall totals are hard to forecast accurately because of a poor understanding of the microphysical processes within the clouds composing the storms. This paper presents a case study from the NASA IMPACTS field campaign that involved two aircraft sampling the storm simultaneously with radars, and probes that measure the microphysical properties within the storm. The paper examines how variations in stability and frontal structure influence the microphysical evolution of ice particles as they fall from cloud top to the surface within the storm. 

Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting. 

Abstract A novel polarimetric radar algorithm for melting-layer (ML) detection and determination of its height has been developed and tested for a large number of cold-season weather events. The algorithm uses radial profiles of the cross-correlation coefficient ( ρ hv or CC) at the lowest elevation angles (<5°–6°). The effects of beam broadening on the spatial distribution of CC have been taken into account via theoretical simulations of the radial profiles of CC assuming intrinsic vertical profiles of polarimetric radar variables within the ML with varying heights and depths of the ML. The model radial profiles of CC and their key parameters are stored in lookup tables and compared with the measured CC profiles. The matching of the model and measured CC radial profiles allows the algorithm to determine the “true” heights of the top and bottom of the ML, H t and H b , at distances up to 150 km from the radar. Integrating the CC information from all available antenna elevations makes it possible to produce accurate maps of H t and H b over large areas of radar coverage as opposed to the previous ML detection methods including the existing algorithm implemented on the U.S. network of the WSR-88Ds. The initial version of the algorithm has been implemented in C++ and tested for a multitude of cold-season weather events characterized by a low ML with different degrees of spatial nonuniformity including cases with sharp frontal boundaries and rain–snow transitions. The new ML detection algorithm (MLDA) exhibits robust performance, demonstrating spatial and temporal continuity, and showing general consistency of the ML designations matching those obtained from the regional model and the quasi-vertical profiles (QVP) methodology output. 

Abstract Snow sublimating in dry air is a forecasting challenge and can delay the onset of surface snowfall and affect storm-total accumulations. Despite this, it remains comparatively less studied than other microphysical processes. Herein, the characteristics of sublimating snow and the potential for nowcasting snowfall reaching the surface are explored through the use of dual-polarization radar. Twelve cases featuring prolific sublimation were analyzed using range-defined quasi-vertical profiles (RDQVPs) and compared with environmental model analyses. Overall, reflectivity Z significantly decreases, differential reflectivity Z DR slightly decreases, and copolar-correlation coefficient ρ hv remains nearly constant through the sublimation layer. Regions of enhanced specific differential phase K dp were frequently observed in the sublimation layer and are believed to be polarimetric evidence of secondary ice production via sublimation. A 1D bin model was initialized using particle size distributions retrieved from the RDQVPs using numerous novel polarimetric snowretrieval relations for a wide range of forecast lead times, with the model environment evolving in response to sublimation. It was found that the model was largely able to predict the snowfall start time up to six hours in advance, with a 6-h median bias of just -18.5 minutes. A more detailed case study of the 08 December 2013 snowstorm in the Philadelphia region was also performed, demonstrating good correspondence with observations and examples of model fields (e.g., cooling rate) hypothetically available from such a tool. The proof-of-concept results herein demonstrate the potential benefits of incorporating spatially averaged radar data in conjunction with simple 1D models into the nowcasting process. 

Abstract The intrinsic uncertainty of radar-based retrievals in snow originates from a large diversity of snow growth habits, densities, and particle size distributions, all of which can make interpreting radar measurements of snow very challenging. The application of polarimetric radar for snow measurements can mitigate some of these issues. In this study, a novel polarimetric method for quantification of the extinction coefficient and visibility in snow, based on the joint use of radar reflectivity at horizontal polarization Z and specific differential phase K DP , is introduced. A large 2D-video-disdrometer snow dataset from central Oklahoma is used to derive a polarimetric bivariate power-law relation for the extinction coefficient, . The relation is derived for particle aspect ratios ranging from 0.5 to 0.8 and the width of the canting angle distribution ranging from 0° to 40°, values typical of aggregated snow, and validated via theoretical and analytical derivations/simulations. The multiplier of the relation is sensitive to variations in particles’ densities, the width of the canting angle distribution, and particles’ aspect ratios, whereas the relation’s exponents are practically invariant to changes in the latter two parameters. This novel approach is applied to polarimetric S-band WSR-88D data and verified against previous studies and in situ measurements of the extinction coefficient for four snow events in the eastern United States. The polarimetric radar estimates of the extinction coefficient exhibit smaller biases in comparison to previous studies concerning the ground measurements. The results indicate that there is good potential for reliable radar estimates of visibility from polarimetric weather radars, a parameter inversely proportional to the extinction coefficient. 

Abstract. Radar dual-wavelength ratio (DWR) measurements from the Stony Brook RadarObservatory Ka-band scanning polarimetric radar (KASPR, 35 GHz), a W-bandprofiling radar (94 GHz), and a next-generation K-band (24 GHz) micro rainradar (MRRPro) were exploited for ice particle identification using triple-frequency approaches. The results indicated that two of the radarfrequencies (K and Ka band) are not sufficiently separated; thus, thetriple-frequency radar approaches had limited success. On the other hand, ajoint analysis of DWR, mean Doppler velocity (MDV), andpolarimetric radar variables indicated potential in identifying ice particletypes and distinguishing among different ice growth processes and even inrevealing additional microphysical details. We investigated all DWR pairs in conjunction with MDV from the KASPRprofiling measurements and differential reflectivity (ZDR) and specificdifferential phase (KDP) from the KASPR quasi-vertical profiles. TheDWR-versus-MDV diagrams coupled with the polarimetric observables exhibiteddistinct separations of particle populations attributed to different rimedegrees and particle growth processes. In fallstreaks, the 35–94 GHz DWRpair increased with the magnitude of MDV corresponding to the scatteringcalculations for aggregates with lower degrees of riming. The DWR valuesfurther increased at lower altitudes while ZDR slightly decreased,indicating further aggregation. Particle populations with higher rimedegrees had a similar increase in DWR but a 1–1.5 m s−1 largermagnitude of MDV and rapid decreases in KDP and ZDR. The analysisalso depicted the early stage of riming where ZDR increased with theMDV magnitude collocated with small increases in DWR. This approach willimprove quantitative estimations of snow amount and microphysical quantitiessuch as rime mass fraction. The study suggests that triple-frequencymeasurements are not always necessary for in-depth ice microphysical studiesand that dual-frequency polarimetric and Doppler measurements cansuccessfully be used to gain insights into ice hydrometeor microphysics. 

Abstract Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s −1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models. 

Abstract A novel way to process polarimetric radar data collected via plan position indicator (PPI) scans and display those data in a time–height format is introduced. The columnar vertical profile (CVP) methodology uses radar data collected via multiple elevation scans, limited to data within a set region in range and azimuth relative to the radar, to create vertical profiles of polarimetric radar data representative of that limited region in space. This technique is compared to others existing in the literature, and various applications are discussed. Polarimetric ice microphysical retrievals are performed on CVPs created within the stratiform rain region of two mesoscale convective systems sampled during two field campaigns, where CVPs follow the track of research aircraft. Aircraft in situ data are collocated to microphysical retrieval data, and the accuracy of these retrievals is tested against other retrieval techniques in the literature.