Physical Modeling of Snow Gliding: A Case Study in the NW Italian Alps

Snow gliding, a slow movement downhill of snow cover, is complex to forecast and model and yet is extremely important, because it drives snowpack dynamics in the pre-avalanching phase. Despite recent interest in this process and the development of some studies therein, this phenomenon is poorly understood and represents a major point of uncertainty for avalanche forecasting. This study presents a data-driven, physically based, time-dependent 1D model, Poli-Glide, able to predict the slow movement of snowpacks along a flow line at the daily scale. The objective of the work was to create a useful snow gliding model, requiring few, relatively easily available input data, by (i) modeling snowpack evolution from measured precipitation and air temperature, (ii) evaluating the rate and extent of movement of the snowpack in the gliding phase, and (iii) assessing fracture (i.e., avalanching) timing. Such a model could be then used to provide hazard assessment in areas subject to gliding, thereby, and subsequent avalanching. To do so, some simplifying assumptions were introduced, namely that (i) negligible traction stress occurs within soil, (ii) water percolation into snow occurs at a fixed rate, and (iii) the micro topography of soil is schematized according to a sinusoidal function in the absence of soil erosion. The proposed model was then applied to the “Torrent des Marais-Mont de La Saxe” site in Aosta Valley, monitored during the winters of 2010 and 2011, featuring different weather conditions. The results showed an acceptable capacity of the model to reproduce snowpack deformation patterns and the final snowpack’s displacement. Correlation analysis based upon observed glide rates further confirmed dependence against the chosen variables, thus witnessing the goodness of the model. The results could be a valuable starting point for future research aimed at including more complex parameterizations of the different processes that affect gliding.

Methodological and conceptual challenges in rare and severe event forecast-verification

There are distinctive methodological and conceptual challenges in rare and severe event (RSE) forecast-verification, that is, in the assessment of the quality of forecasts involving natural hazards such as avalanches or tornadoes. While some of these challenges have been discussed since the inception of the discipline in the 1880s, there is no consensus about how to assess RSE forecasts. This article offers a comprehensive and critical overview of the many different measures used to capture the quality of an RSE forecast and argues that there is only one proper skill score for RSE forecast-verification. We do so by first focusing on the relationship between accuracy and skill and show why skill is more important than accuracy in the case of RSE forecast-verification. Subsequently, we motivate three adequacy constraints for a proper measure of skill in RSE forecasting. We argue that the Peirce Skill Score is the only score that meets all three adequacy constraints. We then show how our theoretical investigation has important practical implications for avalanche forecasting by discussing a recent study in avalanche forecast-verification using the nearest neighbour method. Lastly, we raise what we call the “scope challenge" that affects all forms of RSE forecasting and highlight how and why the proper skill measure is important not only for local binary RSE forecasts but also for the assessment of different diagnostic tests widely used in avalanche risk management and related operations. Finally, our discussion is also of relevance to the thriving research project of designing methods to assess the quality of regional multi-categorical avalanche forecasts.

Modeling spatially distributed snow instability at a regional scale using Alpine3D

Assessing the avalanche danger level requires snow stratigraphy and instability data. As such data are usually sparse, we investigated whether distributed snow cover modeling can be used to provide information on spatial instability patterns relevant for regional avalanche forecasting. Using Alpine3D, we performed spatially distributed simulations to evaluate snow instability for the winter season 2016–17 in the region of Davos, Switzerland. Meteorological data from automatic weather stations were interpolated to 100 m horizontal resolution and precipitation was scaled with snow depth measurements from airborne laser scanning. Modeled snow instability metrics assessed for two different weak layers suggested that the weak layer closer to the snow surface was more variable. Initially, it was less stable than the weak layer closer to the ground, yet it stabilized faster as the winter progressed. In spring, the simulated snowpack on south-facing slopes stabilized faster than on north-facing slopes, in line with the regional avalanche forecast. In the winter months January to March 2017, simulated instability metrics did not suggest that the snowpack on south-facing slopes was more stable, as reported in the regional avalanche forecast. Although a validation with field data is lacking, these model results still show the potential and challenges of distributed modeling for supporting operational avalanche forecasting.

Climate drivers of large magnitude snow avalanche years in the U.S. northern Rocky Mountains

Large magnitude snow avalanches pose a hazard to humans and infrastructure worldwide. Analyzing the spatiotemporal behavior of avalanches and the contributory climate factors is important for understanding historical variability in climate-avalanche relationships as well as improving avalanche forecasting. We used established dendrochronological methods to develop a long-term (1867–2019) regional avalanche chronology for the Rocky Mountains of northwest Montana using tree-rings from 647 trees exhibiting 2134 avalanche-related growth disturbances. We then used principal component analysis and a generalized linear autoregressive moving average model to examine avalanche-climate relationships. Historically, large magnitude regional avalanche years were characterized by stormy winters with positive snowpack anomalies, with avalanche years over recent decades increasingly influenced by warmer temperatures and a shallow snowpack. The amount of snowpack across the region, represented by the first principal component, is shown to be directly related to avalanche probability. Coincident with warming and regional snowpack reductions, a decline of ~ 14% (~ 2% per decade) in overall large magnitude avalanche probability is apparent through the period 1950–2017. As continued climate warming drives further regional snowpack reductions in the study region our results suggest a decreased probability of regional large magnitude avalanche frequency associated with winters characterized by large snowpacks and a potential increase in large magnitude events driven by warming temperatures and spring precipitation.

A data efficient machine learning model for autonomous operational avalanche forecasting

Snow avalanches pose serious hazard to people and property in snow bound mountains. Snow mass sliding downslope can gain sufficient momentum to destroy buildings, uproot trees and kill people. Forecasting and in turn avoiding exposure to avalanches is a much practiced measure to mitigate hazard world over. However, sufficient snow stability data for accurate forecasting is generally difficult to collect. Hence forecasters infer snow stability largely through intuitive reasoning based upon their knowledge of local weather, terrain and sparsely available snowpack observations. Machine learning models may add more objectivity to this intuitive inference process. In this paper we propose a data efficient machine learning classifier using the technique of Random Forest. The model can be trained with significantly lesser training data compared to other avalanche forecasting models and it generates useful outputs to minimise and quantify uncertainty. Besides, the model generates intricate reasoning descriptions which are difficult to observe manually. Furthermore, the model data requirement can be met through automatic systems. The proposed model advances the field by being inexpensive and convenient for operational use due to its data efficiency and ability to describe its decisions besides the potential of lending autonomy to the process.

Tomography-based investigation of concurrent snow creep and isothermal metamorphism

<p>Once on the ground, the microstructure of snow, i.e. the three-dimensional arrangement of ice and pores, quickly evolves with metamorphism and deforms under the overburden of the overlaying snow. Understanding these concurrent processes is important to predict the evolution of the physical and mechanical properties of snow which are crucial for many applications, such as avalanche forecasting. To this end, we monitored oedometric creep tests of snow under isothermal conditions at -8.6°C for about one week with X-ray tomography. We investigated the evolution of recent snow under a constant load of around 4 kPa, where both ice matrix creep and metamorphism are active. Our time-series comprises one of the most highly-resolved images of snow microstructure evolution, with a temporal resolution of 3 h and spatial resolution of 8.5 microns and thousands of images. Interestingly, we observed distinct effects of the overburden and of the vapor transport on the microstructure evolution. In particular, the quantification of the ice bond network through the Euler characteristic and the min-cut surface shows that metamorphism progressively increases the bond size almost independently of the applied overburden, while the application of an overburden yields a rapid increase of the bond coordination number. These distinct impacts exhibit the difficulty to accurately reproduce the time evolution of recent snow by snow cover models, whose snow microstructure representation with density and snow type remains too coarse. </p>

Data-driven automatic predictions of avalanche danger in Switzerland

<p>Avalanche forecasting implies predicting current and future snow instability in time and space. In Switzerland, avalanche bulletins are issued daily during the winter season to warn the public about the avalanche hazard, described by region with one of five danger levels. Assessing avalanche danger is by large a data-driven, yet experience-based decision-making process. It involves analysing a multitude of data diverse in scale – time and space, and concluding by expert judgment on the avalanche scenario. Numerous statistical models were developed in the past, but rarely applied due to limited usefulness in operational forecasting. Modern machine learning techniques open up new possibilities for developing support tools for operational avalanche forecasting. With this aim, we developed a data-driven approach based on the supervised Random Forest (RF) classifier to automatically predict the danger level for dry-snow avalanche conditions in the Swiss Alps. A large database of more than 20 years of meteorological data and modelled snow stratigraphy data obtained with the numerical snow cover model SNOWPACK were used to train the RF algorithm. We optimized the model and selected the best set of input features that combine meteorological variables and features extracted from the simulated profiles, resampled at the same daily resolution as the forecasts. Our target variable was the regional danger level forecast in the public bulletin. We evaluated the predictive performance of the RF model with an independent test set with data of two winter seasons (2018-2019 and 2019-2020). The test set accuracy was 72 %, which is slightly lower than the accuracy estimate of the public forecasts (about 76 %). Given this uncertainty in our target variable, we trained an optimized RF model on a subset containing so-called verified avalanche danger levels. The test set accuracy then increased to 80 %. During the winter season 2020-2021, both RF models were tested in operational setting and automatically predicted a ‘nowcast’ and a ‘forecast’ in real-time.  In parallel, we also tested a deep recurrent neural network model, which used a 7-days time series with 3-hours time resolution as input and also predicted the avalanche danger level. We present a comparison of the performance of the three models. This is one of the first times that a data-driven approach is tested in real-time as a feasible tool for operational avalanche forecasting.</p>

A tracking algorithm to identify slab and weak layer combinations for assessing snow instability and avalanche problem type

<p>Snow cover models have mostly been developed to support avalanche forecasting. Recently developed snow instability metrics can help interpreting modeled snow cover data. However, presently snow cover models cannot forecast the relevant avalanche problem types – an essential element to describe avalanche danger. We present an approach to detect, track and assess weak layers in snow cover model output data to eventually assess the related avalanche problem type. We demonstrate the applicability of this approach with both, SNOWPACK and CROCUS snow cover model output for one winter season at Weissfluhjoch. We introduced a classification scheme for four commonly used avalanche problem types including new snow, wind slabs, persistent weak layers and wet snow, so different avalanche situations during a winter season can be classified based on weak layer type and meteorological conditions. According to the modeled avalanche problem types and snow instability metrics both models produced weaknesses in the modeled stratigraphy during similar periods. For instance, in late December 2014 the models picked up a non-persistent as well as a persistent weak layer that were both observed in the field and caused widespread instability in the area. Times when avalanches released naturally were recorded with two seismic avalanche detection systems, and coincided reasonably well with periods of low modeled stability. Moreover, the presented approach provides the avalanche problem types that relate to the observed natural instability which makes the interpretation of modeled snow instability metrics easier. As the presented approach is process-based, it is applicable to any model in any snow avalanche climate. It could be used to anticipate changes in avalanche problem type due to changing climate. Moreover, the presented approach is suited to support the interpretation of snow stratigraphy data for operational forecasting.</p>

Snow profile alignment and similarity assessment for aggregating, clustering, and evaluating snowpack model output for avalanche forecasting

Snowpack models simulate the evolution of the snow stratigraphy based on meteorological inputs and have the potential to support avalanche risk management operations with complementary information relevant for their avalanche hazard assessment, especially in data-sparse regions or at times of unfavorable weather and hazard conditions. However, the adoption of snowpack models in operational avalanche forecasting has been limited, predominantly due to missing data processing algorithms and uncertainty around model validity. Thus, to enhance the usefulness of snowpack models for the avalanche industry, numerical methods are required that evaluate and summarize snowpack model output in accessible and relevant ways. We present algorithms that compare and assess generic snowpack data from both human observations and models, which consist of multidimensional sequences describing the snow characteristics of grain type, hardness, and age. Our approach exploits Dynamic Time Warping, a well-established method in the data sciences, to match layers between snow profiles and thereby align them. The similarity of the aligned profiles is then evaluated by our independent similarity measure based on characteristics relevant for avalanche hazard assessment. Since our methods provide the necessary quantitative link to data clustering and aggregating methods, we demonstrate how snowpack model output can be grouped and summarized according to similar hazard conditions. By emulating aspects of the human avalanche hazard assessment process, our methods aim to promote the operational application of snowpack models so that avalanche forecasters can begin to build an understanding of how to interpret and trust operational snowpack simulations.

Brief Communication: Initializing RAMMS with High Resolution LiDAR Data for Avalanche Simulations

The Rapid Mass Movements Simulator (RAMMS) is an avalanche dynamics software tool for research and forecasting. Since the model’s conception, the sensitivity of inputs on simulation results has been well-documented. Here, we introduce a new method for initializing RAMMS that can be easily operationalized for avalanche forecasting using high resolution LiDAR data. As a demonstration, hypothetical avalanche simulations were performed while incrementally incorporating semi-automated LiDAR-derived values for snow depth, interface topography, and vegetative cover from field-collected LiDAR data. Results show considerable variation in the calculated runout extent, flow volume, pressure, and velocity of the simulated avalanches when incorporating these LiDAR-derived values.