skip to main content


Title: Polarimetric Observations and Simulations of Sublimating Snow: Implications for Nowcasting
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.  more » « less
Award ID(s):
1841246
NSF-PAR ID:
10316689
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Applied Meteorology and Climatology
ISSN:
1558-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Quasi-vertical profiles (QVPs) of polarimetric radar data have emerged as a powerful tool for studying precipitation microphysics. Various studies have found enhancements in specific differential phase K dp in regions of suspected secondary ice production (SIP) due to rime splintering. Similar K dp enhancements have also been found in regions of sublimating snow, another proposed SIP process. This work explores these K dp signatures for two cases of sublimating snow using nearly collocated S- and Ka-band radars. The presence of the signature was inconsistent between the radars, prompting exploration of alternative causes. Idealized simulations are performed using a radar beam-broadening model to explore the impact of nonuniform beam filling (NBF) on the observed reflectivity Z and K dp within the sublimation layer. Rather than an intrinsic increase in ice concentration, the observed K dp enhancements can instead be explained by NBF in the presence of sharp vertical gradients of Z and K dp within the sublimation zone, which results in a K dp bias dipole. The severity of the bias is sensitive to the Z gradient and radar beamwidth and elevation angle, which explains its appearance at only one radar. In addition, differences in scanning strategies and range thresholds during QVP processing can constructively enhance these positive K dp biases by excluding the negative portion of the dipole. These results highlight the need to consider NBF effects in regions not traditionally considered (e.g., in pure snow) due to the increased K dp fidelity afforded by QVPs and the subsequent ramifications this has on the observability of sublimational SIP. Significance Statement Many different processes can cause snowflakes to break apart into numerous tiny pieces, including when they evaporate into dry air. Purported evidence of this phenomenon has been seen in data from some weather radars, but we noticed it was not seen in data from others. In this work we use case studies and models to show that this signature may actually be an artifact from the radar beam becoming too big and there being too much variability of the precipitation within it. While this breakup process may actually be occurring in reality, these results suggest we may have trouble observing it with typical weather radars. 
    more » « less
  2. null (Ed.)
    Abstract In a 2018 paper by Bukovčić et al., polarimetric bivariate power-law relations for estimating snowfall rate S and ice water content (IWC), and , were developed utilizing 2D video disdrometer snow measurements in Oklahoma. Herein, these disdrometer-based relations are generalized for the range of particle aspect ratios from 0.5 to 0.8 and the width of the canting angle distribution from 0° to 40° and are validated via analytical/theoretical derivations and simulations. In addition, a novel S ( K DP , Z dr ) polarimetric relation utilizing the ratio between specific differential phase K DP and differential reflectivity Z dr , , is derived. Both K DP and are proportionally affected by the ice particles’ aspect ratio and width of the canting angle distribution; therefore, the variables’ ratio tends to be almost invariant to the changes in these parameters. The S ( K DP , Z ) and S ( K DP , Z dr ) relations are applied to the polarimetric S-band WSR-88D data obtained from three geographical locations in Virginia, Oklahoma, and Colorado, and their performance is compared with estimations from the standard S ( Z ) relations and ground snow measurements. The polarimetric estimates of snow accumulations from the three cases exhibit smaller bias in comparison with the S ( Z ), indicating good potential for more reliable radar snow measurements. 
    more » « less
  3. null (Ed.)
    Abstract A polarimetric radar–based method for retrieving atmospheric ice particle shapes is applied to snowfall measurements by a scanning K a -band radar deployed at Oliktok Point, Alaska (70.495°N, 149.883°W). The mean aspect ratio, which is defined by the hydrometeor minor-to-major dimension ratio for a spheroidal particle model, is retrieved as a particle shape parameter. The radar variables used for aspect ratio profile retrievals include reflectivity, differential reflectivity, and the copolar correlation coefficient. The retrievals indicate that hydrometeors with mean aspect ratios below 0.2–0.3 are usually present in regions with air temperatures warmer than approximately from −17° to −15°C, corresponding to a regime that has been shown to be favorable for growth of pristine ice crystals of planar habits. Radar reflectivities corresponding to the lowest mean aspect ratios are generally between −10 and 10 dB Z . For colder temperatures, mean aspect ratios are typically in a range between 0.3 and 0.8. There is a tendency for hydrometeor aspect ratios to increase as particles transition from altitudes in the temperature range from −17° to −15°C toward the ground. This increase is believed to result from aggregation and riming processes that cause particles to become more spherical and is associated with areas demonstrating differential reflectivity decreases with increasing reflectivity. Aspect ratio retrievals at the lowest altitudes are consistent with in situ measurements obtained using a surface-based multiangle snowflake camera. Pronounced gradients in particle aspect ratio profiles are observed at altitudes at which there is a change in the dominant hydrometeor species, as inferred by spectral measurements from a vertically pointing Doppler radar. 
    more » « less
  4. 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. 
    more » « less
  5. Abstract

    This study derives simple analytical expressions for the theoretical height profiles of particle number concentrations (Nt) and mean volume diameters (Dm) during the steady-state balance of vapor growth and collision–coalescence with sedimentation. These equations are general for both rain and snow gamma size distributions with size-dependent power-law functions that dictate particle fall speeds and masses. For collision–coalescence only,Nt(Dm) decreases (increases) as an exponential function of the radar reflectivity difference between two height layers. For vapor deposition only,Dmincreases as a generalized power law of this reflectivity difference. Simultaneous vapor deposition and collision–coalescence under steady-state conditions with conservation of number, mass, and reflectivity fluxes lead to a coupled set of first-order, nonlinear ordinary differential equations forNtandDm. The solutions to these coupled equations are generalized power-law functions of heightzforDm(z) andNt(z) whereby each variable is related to one another with an exponent that is independent of collision–coalescence efficiency. Compared to observed profiles derived from descending in situ aircraft Lagrangian spiral profiles from the CRYSTAL-FACE field campaign, these analytical solutions can on average capture the height profiles ofNtandDmwithin 8% and 4% of observations, respectively. Steady-state model projections of radar retrievals aloft are shown to produce the correct rapid enhancement of surface snowfall compared to the lowest-available radar retrievals from 500 m MSL. Future studies can utilize these equations alongside radar measurements to estimateNtandDmbelow radar tilt elevations and to estimate uncertain microphysical parameters such as collision–coalescence efficiencies.

    Significance Statement

    While complex numerical models are often used to describe weather phenomenon, sometimes simple equations can instead provide equally good or comparable results. Thus, these simple equations can be used in place of more complicated models in certain situations and this replacement can allow for computationally efficient and elegant solutions. This study derives such simple equations in terms of exponential and power-law mathematical functions that describe how the average size and total number of snow or rain particles change at different atmospheric height levels due to growth from the vapor phase and aggregation (the sticking together) of these particles balanced with their fallout from clouds. We catalog these mathematical equations for different assumptions of particle characteristics and we then test these equations using spirally descending aircraft observations and ground-based measurements. Overall, we show that these mathematical equations, despite their simplicity, are capable of accurately describing the magnitude and shape of observed height and time series profiles of particle sizes and numbers. These equations can be used by researchers and forecasters along with radar measurements to improve the understanding of precipitation and the estimation of its properties.

     
    more » « less