Simulating airborne 'snow walls' of Antarctica using CRYOWRF v1.0

When a well-developed, high velocity katabatic flow draining down the ice sheet of Antarctica reaches the coast, it experiences an abrupt and rapid transition due to change in slope resulting in formation of a hydraulic jump. A remarkable manifestation of the hydraulic jump, given the "right" surface conditions, is the large-scale entrainment and convergence of blowing snow particles within the hydraulic jump. This can result in formation of 100-1000 m high, highly localized "walls" of snow in the air in an otherwise cloud-free sky.Recent work by Vignon et al. (2020) has described in detail, the mechanisms resulting in the formation of hydraulic jumps and excitation of gravity waves during a particularly notable event at the Dumont d"Urville (DDU) station in August 2017. They used a combination of satellite images, mesoscale simulations with WRF and station measurements (including Micro Rain Radars) in their study, notably relying on the snow wall for diagnosing and quantifying the hydraulic jump in satellite images. On the other hand, relatively less importance was given towards the surface snow processes including the transport of snow particles in the wall.In this presentation, we present results from simulations done using the recently developed CRYOWRF v1.0 to recreate the August 2017 episode at DDU and explicitly simulate the formation and the dynamics of the snow wall itself. CRYOWRF enhances the standard WRF model with the state-of-the-art surface snow modelling scheme SNOWPACK as well as a completely new blowing snow scheme. SNOWPACK essentially acts as a land surface model for the WRF atmospheric model, thus making a quantum leap over the existing snow cover models in WRF. Since SNOWPACK is a grain-scale snow model, it allows for the proper formulation of boundary conditions for simulating blowing snow dynamics.Results show the formation of the snow wall due to large scale entrainment over a wide area of the ice sheet, the mass balance of the snow wall within the hydraulic jump and finally, the destruction of the snow wall and the ultimate fate of all the entrained snow. We also show results for the influence of the snow wall on the local surface radiation at DDU. Overall, we test the capabilities of CRYOWRF to simulate such a complex phenomenon and highlight possible applications now feasible due the tight coupling of an advanced snow cover model and a multi-scale, non-hydrostatic atmospheric flow solver.

Simulating the influence of snow surface processes on soil moisture dynamics and streamflow generation in an alpine catchment

Mathias Thierry Pierre Bavay, Francesco Comola, Michael Lehning, Nander Wever

The assessment of flood risks in alpine, snow-covered catchments requires an understanding of the linkage between the snow cover, soil and discharge in the stream network. Here, we apply the comprehensive, distributed model Alpine3D to investigate the role of soil moisture in the predisposition of the Dischma catchment in Switzerland to high flows from rainfall and snowmelt. The recently updated soil module of the physics-based multilayer snow cover model SNOWPACK, which solves the surface energy and mass balance in Alpine3D, is verified against soil moisture measurements at seven sites and various depths inside and in close proximity to the Dischma catchment. Measurements and simulations in such terrain are difficult and consequently, soil moisture was simulated with varying degrees of success. Differences between simulated and measured soil moisture mainly arise from an overestimation of soil freezing and an absence of a groundwater description in the Alpine3D model. Both were found to have an influence in the soil moisture measurements. Using the Alpine3D simulation as the surface scheme for a spatially explicit hydrologic response model using a travel time distribution approach for inter-flow and baseflow, streamflow simulations were performed for the discharge from the catchment. The streamflow simulations provided a closer agreement with observed streamflow when driving the hydrologic response model with soil water fluxes at 30 cm depth in the Alpine3D model. Performance decreased when using the 2 cm soil water flux, thereby mostly ignoring soil processes. This illustrates that the role of soil moisture is important to take into account when understanding the relationship between both snowpack runoff and rainfall and catchment discharge in high alpine terrain. However, using the soil water flux at 60 cm depth to drive the hydrologic response model also decreased its performance, indicating that an optimal soil depth to include in surface simulations exists and that the runoff dynamics are controlled by only a shallow soil layer. Runoff coefficients (i.e. ratio of rainfall over discharge) based on measurements for high rainfall and snowmelt events were found to be dependent on the simulated initial soil moisture state at the onset of an event, further illustrating the important role of soil moisture for the hydrological processes in the catchment. The runoff coefficients using simulated discharge were found to reproduce this dependency, which shows that the Alpine3D model framework can be successfully applied to assess the predisposition of the catchment to flood risks from both snowmelt and rainfall events.

Copernicus GmbH2017

Evaluating a prediction system for snow management

Michael Lehning

The evaluation of snowpack models capable of accounting for snow management in ski resorts is a major step towards acceptance of such models in supporting the daily decision-making process of snow production managers. In the framework of the EU Horizon 2020 (H2020) project PROSNOW, a service to enable real-time optimization of grooming and snow-making in ski resorts was developed. We applied snow management strategies integrated in the snowpack simulations of AMUNDSEN, Crocus, and SNOWPACK-Alpine3D for nine PROSNOW ski resorts located in the European Alps. We assessed the performance of the snow simulations for five winter seasons (2015-2020) using both ground-based data (GNSS-measured snow depth) and spaceborne snow maps (Copernicus Sentinel-2). Particular attention has been devoted to characterizing the spatial performance of the simulated piste snow management at a resolution of 10 m. The simulated results showed a high overall accuracy of more than 80% for snow-covered areas compared to the Sentinel-2 data. Moreover, the correlation to the ground observation data was high. Potential sources for local differences in the snow depth between the simulations and the measurements are mainly the impact of snow redistribution by skiers; compensation of uneven terrain when grooming; or spontaneous local adaptions of the snow management, which were not reflected in the simulations. Subdividing each individual ski resort into differently sized ski resort reference units (SRUs) based on topography showed a slight decrease in mean deviation. Although this work shows plausible and robust results on the ski slope scale by all three snowpack models, the accuracy of the results is mainly dependent on the detailed representation of the real-world snow management practices in the models. As snow management assessment and prediction systems get integrated into the workflow of resort managers, the formulation of snow management can be refined in the future.

COPERNICUS GESELLSCHAFT MBH2021

A review in snow saltation dynamics and its implications for the surface mass balance

Daniela Brito Melo, Francesco Comola, Franziska Gerber, Michael Lehning, Varun Sharma, Armin Sigmund

Snow covered regions are frequently subjected to strong winds. This leads to the erosion of the snow surface and the occurrence of drifting and blowing snow events. To correctly predict and model these phenomena is of the utmost importance to assess snowpack stability and avalanche formation, as well as airborne snow sublimation and the resultant surface mass balance. Nonetheless, snow transport is frequently neglected or misrepresented in regional and mesoscale models. One of the main challenges is the accurate representation of snow transport close to the ground, where snow particles are transported by a process called saltation. This shallow layer comprises most of the horizontal mass flux and sets the lower boundary condition to model snow suspension clouds. A detailed study of snow saltation dynamics has been conducted using a Large-Eddy-Simulation flow solver coupled with a Lagrangian model for particle trajectories. The effect of particle size distribution and interparticle cohesion on particle speed, mass flux and surface friction velocity has systematically been investigated. The results show that snow cohesion and grain size heterogeneity can significantly increase saltation mass flux, specially at high friction velocities. Moreover, the agreement between simulation results and the saltation models typically used in large scale atmospheric models is highly dependent on the assumed bed characteristics. These findings can support the development of comprehensive saltation models that specifically take into account snow bed properties. These new saltation models can be implemented in mesoscale atmospheric models and have the potential to significantly improve surface mass balance predictions if grain-scale snow surface properties are available. This is possible in the newly developed CRYOWRF model which expands WRF’s surface modeling suite by including the advanced complexity, grain-scale snow model, SNOWPACK, along with a blowing snow scheme.

2021