The Impact of Diffusive Water Vapor Transport on Snow Profiles in Deep and Shallow Snow Covers and on Sea Ice

Water vapor transport has been highlighted as a critical process in Arctic snowpacks, shaping the snow cover structure in terms of density, thermal conductivity, and temperature profile among others. Here, we present an attempt to describe the thermally-induced vertical diffusion of water vapor in the snow cover and its effects of the snowpack structure using the SNOWPACK model. Convection, that may also constitute a significant part of vapor transport, is not addressed. Assuming saturated conditions at the upper boundary of the snowpack and as initial condition, the vapor flux between snow layers is expressed by a 1-dimensional transient diffusion equation, which is solved with a finite difference routine. The implications on the snowpack of this vertical diffusive flux, are analyzed using metrics such as the cumulative density change due to diffusive vapor transport, the degree of over- or undersaturation, the instantaneous snow density change rate, and the percentage of snow density change. We present results for four different regions sampling the space of natural snow cover variability: Alpine, Subarctic, Arctic, and Antarctic sea ice. The largest impact of diffusive water vapor transport is observed in snow on sea ice in the Weddell Sea and the shallow Arctic snowpack. The simulations show significant density reductions upon inclusion of diffusive water vapor transport: cumulative density changes from diffusive vapor transport can reach −62 and −66 kg m−3 for the bottom layer in the shallow Arctic snowpack and snow on sea ice, respectively. For comparison, in deeper snow covers, they rarely exceed −40 kg m−3. This leads to changes in density for shallow snowpacks at the soil-snow interface in the range of −5 to −21%. Mirroring the density decease at depth is a thicker deposition layer above it with increase in density around 7.5%. Similarly, for the sea ice, the density decreased at the sea ice-snow interface by −20%. We acknowledge that vapor transport by diffusion may in some snow covers—such as in thin tundra snow—be small compared to convective transport, which will have to be addressed in future work.

Convection of water vapour in snowpacks

Mahdi Jafari, Michael Lehning, Varun Sharma

This paper studies numerically the convection of water vapour in snowpacks using an Eulerian-Eulerian two-phase approach. The convective water vapour transport in snow and its effects on snow density are often invoked to explain observed density profiles, e.g. of thin Arctic snow covers, but this process has never been numerically simulated and analysed in a systematic manner. Here, the impact of convection on the thermal and phase change regimes as a function of different snowpack depths, thermal boundary conditions and Rayleigh numbers is analysed. We find considerable impact of natural convection on the snow density distribution with a layer of significantly lower density at the bottom of the snowpack and a layer of higher density located higher in the snowpack or at the surface. We find that emergent heterogeneity in the snow porosity results in a feedback effect on the spatial organization of convection cells causing their horizontal displacement. Regions where the snowpack is most impacted by phase changes are found to be horizontally extended and vertically thin, `pancake'-like layers at the top and bottom of the snowpack. We further demonstrate that among the parameters important for natural convection, the snowpack depth has the strongest influence on the heat and mass transfer. Despite our simplifying assumptions, our study represents a significant improvement over the state of the art and a first step to accurately simulate snowpack dynamics in diverse regions of the cryosphere.

CAMBRIDGE UNIV PRESS2022

Version 1 of a sea ice module for the physics-based, detailed, multi-layer SNOWPACK model

Michael Lehning, Katherine Colby Leonard, Nander Wever

Sea ice is an important component of the global climate system. The presence of a snowpack covering sea ice can strongly modify the thermodynamic behavior of the sea ice, due to the low thermal conductivity and high albedo of snow. The snowpack can be stratified and change properties (density, water content, grain size and shape) throughout the seasons. Melting snow provides freshwater which can form melt ponds or cause flushing of salt out of the underlying sea ice, while flooding of the snow layer by saline ocean water can strongly impact both the ice mass balance and the freezing point of the snow. To capture the complex dynamics from the snowpack, we introduce modifications to the physics-based, multi-layer SNOWPACK model to simulate the snow-sea-ice system. Adaptations to the model thermodynamics and a description of water and salt transport through the snow-sea-ice system by coupling the transport equation to the Richards equation were added. These modifications allow the snow microstructure descriptions developed in the SNOWPACK model to be applied to sea ice conditions as well. Here, we drive the model with data from snow and ice mass-balance buoys installed in the Weddell Sea in Antarctica. The model is able to simulate the temporal evolution of snow density, grain size and shape, and snow wetness. The model simulations show abundant depth hoar layers and melt layers, as well as superimposed ice formation due to flooding and percolation. Gravity drainage of dense brine is underestimated as convective processes are so far neglected. Furthermore, with increasing model complexity, detailed forcing data for the simulations are required, which are difficult to acquire due to limited observations in polar regions.

COPERNICUS GESELLSCHAFT MBH2020