skip to main content

Explore Research Products in the NSF-PAR

Title: How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates
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.  more » « less
Award ID(s):
1841246
NSF-PAR ID:
10209864
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
Volume:
78
Issue:
1
ISSN:
0022-4928
Page Range / eLocation ID:
51 to 73
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ( , Journal of Hydrometeorology)
    null (Ed.)
    Abstract This study investigates the shapes and fall speeds of freezing and frozen raindrops through field observations using an instrument called the high-speed optical disdrometer (HOD) that we developed recently. Our field observations showed that while the shapes of all of the observed freezing raindrops and a portion of the frozen raindrops (39% of the frozen raindrops that are larger than 1.0 mm in volume equivalent diameter D ) resemble the shapes of warm raindrops, majority of frozen raindrops (61% of the frozen raindrops with D > 1.0 mm) exhibited a distinct feature such as a spicule, bulge, cavity, or aggregation. Field observations of axis ratios (i.e., ratio of the vertical to horizontal chord) and fall speeds were compared with the predictions of available models. Separate empirical axis ratio parameterizations were developed for the freezing and frozen raindrops using the HOD field observations and extensions to an available shape model were also incorporated. For the fall speeds of freezing and frozen raindrops, field observations demonstrated a good agreement with the predictions of the available parameterizations. Frozen raindrops showed a larger scatter of fall speeds around the mean fall speed of a given drop size than those of the freezing raindrops due to the shape variety among the frozen raindrops with the aforementioned distinct features. The drag coefficients for the observed hydrometeors were compared with the predictions of the available drag coefficient models. Separate “drag coefficient–Reynolds number” relationships for freezing and frozen raindrops were developed. 
    more » « less
  2. ; ; ; ; ; ; ( , Atmospheric Chemistry and Physics)
    Abstract. Cirrus clouds that form in the tropical tropopause layer(TTL) can play a key role in vertical transport through the uppertroposphere and lower stratosphere, which can significantly impact theradiative energy budget and stratospheric chemistry. However, the lack ofrealistic representation of natural ice cloud habits in microphysicalparameterizations can lead to uncertainties in cloud-related processes andcloud–climate feedbacks. The main goal of this study is to investigate therole of different cloud regimes and the associated ice habits in regulatingthe properties of the TTL. We compare aircraft measurements from theStratoClim field campaign to a set of numerical experiments at the scale of large-eddy simulations (LESs) for the same case study that employ differentmicrophysics schemes. Aircraft measurements over the southern slopes of theHimalayas captured high ice water content (HIWC) up to 2400 ppmv and iceparticle aggregates exceeding 700 µm in size with unusually longresidence times. The observed ice particles were mainly of liquid origin,with a small amount formed in situ. The corresponding profile of ice water content (IWC) fromthe ERA5 reanalysis corroborates the presence of HIWC detrained from deep-convective plumes in the TTL but underestimates HIWC by an order ofmagnitude. In the TTL, only the scheme that predicts ice habits canreproduce the observed HIWC, ice number concentration, and bimodal iceparticle size distribution. The lower range of particle sizes is mostlyrepresented by planar and columnar habits, while the upper range isdominated by aggregates. Large aggregates with sizes between 600 and 800 µm have fall speeds of less than 20 cm s−1, which explains thelong residence time of the aggregates in the TTL. Planar ice particles ofliquid origin contribute substantially to HIWC. The columnar and aggregatehabits are in the in situ range with lower IWC and number concentrations. Forall habits, the ice number concentration increases with decreasingtemperature. For the planar ice habit, relative humidity is inverselycorrelated with fall speed. This correlation is less evident for the othertwo ice habits. In the lower range of supersaturation with respect to ice,the columnar habit has the highest fall speed. The difference in ice numberconcentration across habits can be up to 4 orders of magnitude, withaggregates occurring in much smaller numbers. We demonstrate and quantifythe linear relationship between the differential sedimentation of pristineice crystals and the size of the aggregates that form when pristine crystalscollide. The slope of this relationship depends on which pristine ice habitsediments faster. Each simulated ice habit is associated with distinctradiative and latent heating rates. This study suggests that a modelconfiguration nested down to LES scales with a microphysicalparameterization that predicts ice shape evolution is crucial to provide anaccurate representation of the microphysical properties of TTL cirrus andthus the associated (de)hydration process. 
    more » « less
  3. ; ; ; ; ; ( , Entropy)
    null (Ed.)
    We develop a simple Quantile Spacing (QS) method for accurate probabilistic estimation of one-dimensional entropy from equiprobable random samples, and compare it with the popular Bin-Counting (BC) and Kernel Density (KD) methods. In contrast to BC, which uses equal-width bins with varying probability mass, the QS method uses estimates of the quantiles that divide the support of the data generating probability density function (pdf) into equal-probability-mass intervals. And, whereas BC and KD each require optimal tuning of a hyper-parameter whose value varies with sample size and shape of the pdf, QS only requires specification of the number of quantiles to be used. Results indicate, for the class of distributions tested, that the optimal number of quantiles is a fixed fraction of the sample size (empirically determined to be ~0.25–0.35), and that this value is relatively insensitive to distributional form or sample size. This provides a clear advantage over BC and KD since hyper-parameter tuning is not required. Further, unlike KD, there is no need to select an appropriate kernel-type, and so QS is applicable to pdfs of arbitrary shape, including those with discontinuous slope and/or magnitude. Bootstrapping is used to approximate the sampling variability distribution of the resulting entropy estimate, and is shown to accurately reflect the true uncertainty. For the four distributional forms studied (Gaussian, Log-Normal, Exponential and Bimodal Gaussian Mixture), expected estimation bias is less than 1% and uncertainty is low even for samples of as few as 100 data points; in contrast, for KD the small sample bias can be as large as −10% and for BC as large as −50%. We speculate that estimating quantile locations, rather than bin-probabilities, results in more efficient use of the information in the data to approximate the underlying shape of an unknown data generating pdf. 
    more » « less
  4. ; ( , Journal of the Atmospheric Sciences)
    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
  5. ; ( , Journal of Geophysical Research: Atmospheres)
    Abstract

    Detailed ground‐based observations of snow are scarce in remote regions, such as the Arctic. Here, Multi‐Angle Snowflake Camera measurements of over 55,000 solid hydrometeors—obtained during a two‐year period from August 2016 to August 2018 at Oliktok Point, Alaska—are analyzed and compared to similar measurements from an earlier experiment at Alta, Utah. In general, distributions of hydrometeor fall speed, fall orientation, aspect ratio, flatness, and complexity (i.e., riming degree) were observed to be very similar between the two locations, except that Arctic hydrometeors tended to be smaller. In total, the slope parameter defining a negative exponential of the size distribution was approximately 50% steeper in the Arctic as at Alta. Sixty‐six percent of particles were observed to be rimed or moderately rimed with some suggestion that riming is favored by weak boundary layer stability. On average, the fall speed of rimed particles was not notably different from aggregates. However, graupel density and fall speed increase as cloud temperatures approach the melting point.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email