Search

The influence of drifting and blowing snow on surface mass and energy exchange is difficult to quantify due to limitations in both measurements and models, but is still potentially very important over large areas with seasonal or perennial snow cover. We present a unique set of measurements that make possible the calculation of turbulent moisture, heat, and momentum fluxes during conditions of drifting and blowing snow. From the data, Monin–Obukhov estimation of bulk fluxes is compared to eddy-covariance-derived fluxes. In addition, large-eddy simulations with sublimating particles are used to more completely understand the vertical profiles of the fluxes. For a storm period at the Syowa S17 station in East Antarctica, the bulk parametrization severely underestimates near-surface heat and moisture fluxes. The large-eddy simulations agree with the eddy-covariance fluxes when the measurements are minimally disturbed by the snow particles. We conclude that overall exchange over snow surfaces is much more intense than current models suggest, which has implications for the total mass balance of the Antarctic ice sheet and the cryosphere. Supplementary Information The online version supplementary material available at 10.1007/s10546-021-00653-x.

While turbulent bursts are considered critical for blowing-snow transport and initiation, the interaction of the airflow with the snow surface is not fully understood. To better characterize the coupling of turbulent structures and blowing-snow transport, observations collected in natural environments at the necessary high-resolution time scales are needed. To address this, high-frequency measurements of turbulence, blowing-snow density and particle velocity were made in the Canadian Rockies. During blowing-snow storms, modified variable-interval time averaging enabled identification of periods of near-surface blowing-snow coupling with shear-stress-producing motions in the lowest 2 m of the atmospheric surface layer. The identification of those turbulent motions responsible for blowing snow yields a better understanding of the event-driven mechanics of initiation and sustained transport. The type of coherent structures generating the Reynolds stress are just as important as the magnitude of the Reynolds stress in initiating and sustaining near-surface blowing snow. Our results suggest that blowing-snow models driven by merely the time-averaged shear stress lack physical realism in the near-surface region. The next phase of the development of blowing-snow models should incorporate parametrizations of coherent turbulent structures.

A model is proposed to determine the electric field strength in blowing snow. To test this model, the electric field strength was measured over an 80-day period during the Canadian Arctic Shelf Exchange Study (CASES) in 2004. The electric field strength at 0.5 m correlates well with the difference between 10-m wind speed and a threshold wind speed, although there is a large amount of variation between the electric fields generated during different blowing snow events. Although the model predicts that the electric field should be proportional to particle number density, the correlation is weak. The correlation of wind speed and electric field strength suggests that particles become charged primarily due to friction-induced temperature difference as they impact upon the surface. The strength of the electric field is likely influenced by a large number of other factors that are difficult to measure. However, the model predicts electric field strengths in excess of 25 kV m−1 near the surface, which would have a significant effect on particle motion.

This paper presents a new triple-moment blowing snow model PIEKTUK-T by including predictive equations for three moments of the gamma size distribution. Specifically, predictive equations for the total number concentration, total mass mixing ratio, and total radar reflectivity for blowing snow are included. Tests in the context of idealized experiments and observed case studies demonstrate that the triple-moment model performs better than the double-moment model PIEKTUK-D in predicting the evolution of the number concentration, mixing ratio, shape parameter, and visibility in blowing snow, provided that the fall velocities for the total number concentration, mass mixing ratio, and radar reflectivity are weighted by the same order of the respective moments in both models. The power law relationship between the radar reflectivity factor and particle extinction coefficient found in PIEKTUK-T is consistent with one observed in snow storms. Coupling of the triple-moment blowing snow model to an atmospheric model would allow realistic studies of the effect of blowing snow on weather and climate.

A linear theory is used to model windflow over a tundra landscape, taken as an example of complex topography. The autocorrelation between gradients at different points in this landscape has an exponential dependence on separation, characterized by a correlation length. Areal variances of simulated windspeeds are proportional to the variance of the topographic gradient. Simulations with synthetic topographies having the same correlation structure show that the constants of proportionality in these relationships depend on the correlation length. Statistics of snow sublimation and transport calculated using simulated windspeed distributions in a distributed blowing-snow model are found to be in good agreement for real and synthetic topographies.

We describe in this paper the development of a double-moment modelof blowing snow and its application to the Canadian Arctic. Wefirst outline the formulation of the numerical model, whichsolves a prognostic equation for both the blowing snow mixingratio and total particle numbers, both moments of particles thatare gamma-distributed. Under idealized simulations, the modelyields realistic evolutions of the blowing snow particledistributions, transport and sublimation rates as well as the thermodynamic fields at low computational costs. A parametrizationof the blowing snow sublimation rate is subsequently derived. The model and parametrization are then applied to a Canadian Arctictundra site prone to frequent blowing snow events. Over a period of210 days during the winter of 1996/1997, the near-surfacerelative humidity consistently approaches saturationwith respect to ice. These conditions limit snowpack erosion byblowing snow sublimation to ≈3 mm snow water equivalent (swe)with surface sublimation removing an additional 7 mm swe.We find that our results are highly sensitiveto the proper assimilation of the humidity measurements and the evolving thermodynamic fields in the atmospheric boundary layer during blowingsnow. These factors may explain the lower values of blowing snow sublimationreported in this paper than previously published for the region.

A mesoscale atmospheric numerical model is coupled with a physically based snow-pack model and with a snow-drift model. The snow model is verified for the French Alps by comparing its simulations to observations performed at the Col de Forte in the Chartreuse Massif. The snow erosion threshold depends on snow-pack properties such as density, dendricity, sphericity and particle size. The atmospheric turbulence scheme is modified in order to take into account stabilization effects due to airborne blown snow particles. In particular, vertically integrated stability functions for the stable boundary layer are completed by including the threshold friction velocity for snow erosion. The snow-drift model is calibrated by simulating the conditions observed during the Byrd snow project, held in West Antarctica in 1962. Finally, sensitivity experiments to the snow-surface properties show the importance of their accurate representation when modelling the contribution of deflation to the Antarctic surface mass balance.

Four one-dimensional, time-dependent blowing snow models areintercompared. These include three spectral models, PIEKTUK-T,WINDBLAST, SNOWSTORM, and the bulk version of PIEKTUK-T,PIEKTUK-B. Although the four models are based on common physicalconcepts, they have been developed by different research groups. Thestructure of the models, numerical methods, meteorological field treatmentand the parameterization schemes may be different. Under an agreed standardcondition, the four models generally give similar results for the thermodynamic effects of blowing snow sublimation on the atmospheric boundary layer, including an increase of relative humidity and a decrease of the ambient temperature due to blowing snow sublimation. Relative humidity predicted by SNOWSTORM is lower than the predictions of the other models, which leads to a larger sublimation rate in SNOWSTORM. All four models demonstrate that sublimation rates in a column of blowing snow have a single maximum in time, illustrating self-limitation of the sublimation process of blowing snow. However, estimation of the eddy diffusioncoefficient for momentum (Km), and thereby the diffusion coefficients for moisture (Kw) and for heat (Kh), has a significant influence on the process. Sensitivitytests with PIEKTUK-T show that the sublimation rate can be approximately constant with time after an initial phase, if Km is a linear function with height. In order to match the model results with blowing snow observations, some parameters in the standard run, such as settling velocity of blowing snow particles in these models, may need to be changed to more practical values.

We present in this paper a simple and computationally efficient numerical model that depicts a column of sublimating, blowing snow. This bulk model predicts the mixing ratio of suspended snow by solving an equation that considers the diffusion, settling and sublimation of blowing snow in a time-dependent mode. The bulk model results compare very well with those of a previous spectral version of the model, while increasing its computational efficiency by a factor of about one hundred. This will allow the use of the model to estimate the effects of blowing snow upon the atmospheric boundary layer and to the mass balance of such regions as the Mackenzie River Basin of Canada.

Separation of electrostatic charge during the transport of particles by wind adds a force to the gravitational and fluid forces that determine trajectories of particles moving by saltation. Evaluating this electrostatic force requires the electric field strength very near the saltation surface, and charge-to-mass ratios for the moving particles. Field mill readings 4 cm above the surface in a moderate blizzard showed electric field strength as high as +30 kV m-1. Another experiment gave charge-to-mass ratios of individual saltation particles in low-level drifting that ranged between +72 μC kg-1 and -208 μC kg-1. From these measurements, we estimated electrostatic forces as large as the gravitational force on some saltating particles. Including forces of this magnitude in the equations of motion significantly alters predicted saltation trajectories from those for uncharged particles. Predictions appear reasonable that for some saltating particles, the electrostatic force prevents new surface impacts. These results should help improve models of energy transfer in the planetary boundary layer during blizzards and sandstorms.

A seasonal snowcover blankets much of Canada during wintertime. In such an environment, the frequency of blowing snow events is relatively high and can have important meteorological and hydrological impacts. Apart from the transport of snow, the thermodynamic impact of sublimating blowing snow in air near the surface can be investigated. Using a time or fetch-dependent blowing snow model named 'PIEKTUK' that incorporates prognostic equations for a spectrum of sublimating snow particles, plus temperature and humidity distributions, it is found that the sublimation of blowing snow can lead to temperature decreases of the order of 0.5 °C and significant water vapour increases in the near-surface air. Typical predicted snow removal rates due to sublimation of blowing snow are several millimetres snow water equivalent per day over open Arctic tundra conditions. The model forecast sublimation rates are most sensitive to humidity, as well as wind speed, temperature and particle distributions, with a maximum value in sublimation typically found approximately 1 km downstream from blowing snow initiation. This suggests that the sublimation process is self-limiting despite ongoing transport of snow by wind, yielding significantly lower values of blowing snow sublimation rates (nearly two-thirds less) compared to situations where the thermodynamic feedbacks are neglected. The PIEKTUK model may provide the necessary thermodynamic inputs or blowing snow parameterizations for mesoscale models, allowing the assessment of the contribution of blowing snow fluxes, in more complex situations, to the moisture budgets of high-latitude regions.

A density current model was extended for use in katabatic flow over the steep slopes of Antarctica through the inclusion of the Coriolis effect and weight flux terms corresponding to blowing snow and cooling caused by sublimation. The model was calibrated and tested against data obtained during two flights in Adelie Land, Antarctica, along a trajectory starting about 170 km inland and extending to Dumont d'Urville. The predicted trend in water vapor flux agrees with measurements of this flux, lending support to empirical formulae for both snow flux and sublimation rate. Model predictions of velocity were in good agreement with measured quantities when reasonable estimates of radiation divergence and surface heat exchange were provided as input to the model. The potential temperature gradient above the katabatic layer was found to play a major role in flow stability for high velocity and deep katabatic flows. Velocity predictions were in better agreement with the data when a locally determined value was used for the coefficient in the empirical snow flux expression.

To estimate the sublimation rate of snow during relocation by wind, sizes and concentration of ice crystal fragments were measured at 6 levels in the lowest 1 m, during ten 10-min runs, in a nocturnal blizzard. Power-law functions of height described the decrease in mean particle diameter and concentration. The vertical gradient of water vapor, measured with a thermocouple psychrometer, was approximately linear from 0.2 to 1.0m above the surface. Evaporation of blowing snow over 3 km of transport distance was estimated to be 39% of transport rate, under conditions of the experiment.