skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. 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. Abstract Hail research and forecasting models necessarily involve explicit or implicit—and uncertain—physical assumptions regarding hailstones’ shape, tumbling behavior, fall speed, and thermal energy transfer. Whereas most models assume spherical hailstones, we relax this assumption by using hailstone shape data from field observations to establish empirical size–shape relationships with reasonable degrees of randomness considering hailstones’ natural shape variability, capturing the observed distribution of triaxial ellipsoidal shapes. We also incorporate explicit, random tumbling of individual hailstones during their growth to simulate their free-falling behavior and the resultant changes in cross-sectional area (which affects growth by hydrometeor collection). These physical attributes are incorporated in calculating hailstones’ fall speeds, using either empirical relationships or analytical relationships based on each hailstone’s Best and Reynolds numbers. Options for drag coefficient modification are added to emulate hailstones’ rough surfaces (lobes), which then modifies their thermal energy and vapor exchange with the environment. We investigate how applying these physical assumptions about nonspherical hail to the Penn State hail growth trajectory model, coupled with Cloud Model 1 supercell simulations, impacts hail production and examine the reasons behind the resulting variability in hail statistics. The choice of hailstone size–mass relation and fall speed scheme have the strongest influence on hail sizes. Using nonspherical, tumbling hailstones reduces the number of large hailstones produced. Applying shape-specific thermal energy transfer coefficients subtly increases sizes; the effects of lobes vary depending on the fall speed scheme used. These physical assumptions, although adding complexity to modeling, can be parameterized efficiently and potentially used in microphysics schemes. Significance StatementIn numerical modeling of hailstorms, we usually consider hailstones to be spherical to simplify calculations, but in nature, hailstones generally are not spheres and can be rather lumpy and have spikes. The purpose of this study is to examine how the model result would change when nonspherical hailstone shape is implemented. We examine the relationship between hailstone shape and physical processes during hail growth in effort to explain why these changes occur and offer insights on how nonspherical hailstone shape may be parameterized in bulk microphysics schemes. 
    more » « less
  3. 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
  4. null (Ed.)
    The effect of turbulence on snow precipitation is not incorporated into present weather forecasting models. Here we show evidence that turbulence is in fact a key influence on both fall speed and spatial distribution of settling snow. We consider three snowfall events under vastly different levels of atmospheric turbulence. We characterize the size and morphology of the snow particles, and we simultaneously image their velocity, acceleration and relative concentration over vertical planes approximately $$30\ \textrm {m}^2$$ in area. We find that turbulence-driven settling enhancement explains otherwise contradictory trends between the particle size and velocity. The estimates of the Stokes number and the correlation between vertical velocity and local concentration are consistent with the view that the enhanced settling is rooted in the preferential sweeping mechanism. When the snow vertical velocity is large compared to the characteristic turbulence velocity, the crossing trajectories effect results in strong accelerations. When the conditions of preferential sweeping are met, the concentration field is highly non-uniform and clustering appears over a wide range of scales. These clusters, identified for the first time in a naturally occurring flow, display the signature features seen in canonical settings: power-law size distribution, fractal-like shape, vertical elongation and large fall speed that increases with the cluster size. These findings demonstrate that the fundamental phenomenology of particle-laden turbulence can be leveraged towards a better predictive understanding of snow precipitation and ground snow accumulation. They also demonstrate how environmental flows can be used to investigate dispersed multiphase flows at Reynolds numbers not accessible in laboratory experiments or numerical simulations. 
    more » « less
  5. null (Ed.)
    Abstract. Ground-based measurements of frozen precipitation are heavily influenced by interactions of surface winds with gauge-shield geometry. The Multi-Angle Snowflake Camera (MASC), which photographs hydrometeors in free-fall from three different angles while simultaneously measuring their fall speed, has been used in the field at multiple midlatitude and polar locations both with and without wind shielding. Here, we present an analysis of Arctic field observations – with and without a Belfort double Alter shield – and compare the results to computational fluid dynamics (CFD) simulations of the airflow and corresponding particle trajectories around the unshielded MASC. MASC-measured fall speeds compare well with Ka-band Atmospheric Radiation Measurement (ARM) Zenith Radar (KAZR) mean Doppler velocities only when winds are light (≤5ms-1) and the MASC is shielded. MASC-measured fall speeds that do not match KAZR-measured velocities tend to fall below a threshold value that increases approximately linearly with wind speed but is generally <0.5ms-1. For those events with wind speeds ≤1.5ms-1, hydrometeors fall with an orientation angle mode of 12∘ from the horizontal plane, and large, low-density aggregates are as much as 5 times more likely to be observed. Simulations in the absence of a wind shield show a separation of flow at the upstream side of the instrument, with an upward velocity component just above the aperture, which decreases the mean particle fall speed by 55 % (74 %) for a wind speed of 5 m s−1 (10 m s−1). We conclude that accurate MASC observations of the microphysical, orientation, and fall speed characteristics of snow particles require shielding by a double wind fence and restriction of analysis to events where winds are light (≤5ms-1). Hydrometeors do not generally fall in still air, so adjustments to these properties' distributions within natural turbulence remain to be determined. 
    more » « less