scholarly journals Sensitivity of modeled snow stability data to meteorological input uncertainty

2020 ◽  
Vol 20 (11) ◽  
pp. 2873-2888 ◽  
Author(s):  
Bettina Richter ◽  
Alec van Herwijnen ◽  
Mathias W. Rotach ◽  
Jürg Schweizer

Abstract. To perform spatial snow cover simulations for numerical avalanche forecasting, interpolation and downscaling of meteorological data are required, which introduce uncertainties. The repercussions of these uncertainties on modeled snow stability remain mostly unknown. We therefore assessed the contribution of meteorological input uncertainty to modeled snow stability by performing a global sensitivity analysis. We used the numerical snow cover model SNOWPACK to simulate two snow instability metrics, i.e., the skier stability index and the critical crack length, for a field site equipped with an automatic weather station providing the necessary input for the model. Simulations were performed for a winter season, which was marked by a prolonged dry period at the beginning of the season. During this period, the snow surface layers transformed into layers of faceted and depth hoar crystals, which were subsequently buried by snow. The early-season snow surface was likely the weak layer of many avalanches later in the season. Three different scenarios were investigated to better assess the influence of meteorological forcing on snow stability during (a) the weak layer formation period, (b) the slab formation period, and (c) the weak layer and slab formation period. For each scenario, 14 000 simulations were performed, by introducing quasi-random uncertainties to the meteorological input. Uncertainty ranges for meteorological forcing covered typical differences observed within a distance of 2 km or an elevation change of 200 m. Results showed that a weak layer formed in 99.7 % of the simulations, indicating that the weak layer formation was very robust due to the prolonged dry period. For scenario a, modeled grain size of the weak layer was mainly sensitive to precipitation, while the shear strength of the weak layer was sensitive to most input variables, especially air temperature. Once the weak layer existed (case b), precipitation was the most prominent driver for snow stability. The sensitivity analysis highlighted that for all scenarios, the two stability metrics were mostly sensitive to precipitation. Precipitation determined the load of the slab, which in turn influenced weak layer properties. For cases b and c, the two stability metrics showed contradicting behaviors. With increasing precipitation, i.e., deep snowpacks, the skier stability index decreased (became less stable). In contrast, the critical crack length increased with increasing precipitation (became more stable). With regard to spatial simulations of snow stability, the high sensitivity to precipitation suggests that accurate precipitation patterns are necessary to obtain realistic snow stability patterns.

2020 ◽  
Author(s):  
Bettina Richter ◽  
Alec van Herwijnen ◽  
Mathias W. Rotach ◽  
Jürg Schweizer

Abstract. To perform spatial snow cover simulations for numerical avalanche forecasting, interpolation and downscaling of meteorological data are required, which introduce uncertainties. The repercussions of these uncertainties on modeled snow stability remain mostly unknown. We therefore assessed the contribution of meteorological input uncertainty on modeled snow stability by performing a global sensitivity analysis. We used the numerical snow cover model SNOWPACK to simulate two snow instability metrics, i.e. the skier stability index and the critical crack length, for a field site equipped with an automatic weather station providing the necessary input for the model. Uncertainty ranges for meteorological forcing covered typical differences observed within a distance of 2 km and an elevation change of 200 m. Three different scenarios were investigated to better assess the influence of meteorological forcing on snow stability during a) the weak layer formation period, b) the slab formation period, and c) the weak layer and slab formation period. For each scenario, 14 000 simulations were performed, by introducing quasi-random uncertainties to the meteorological input. Results showed that a weak layer formed in 99.7 % of the simulations, indicating that the weak layer formation was very robust due to the prolonged dry period. For scenario a), modeled grain size of the weak layer was mainly sensitive to precipitation, while the shear strength of the weak layer was sensitive to most input variables, especially air temperature. Once the weak layer existed (case b), precipitation was the most prominent driver for snow stability. The sensitivity analysis highlighted that for all scenarios, the two stability metrics were mostly sensitive precipitation. Precipitation determined the load of the slab, which in turn influenced weak layer properties. For case b) and c), the two stability metrics showed contradicting behaviors. With increasing precipitation, i.e. deep snowpacks, the skier stability index decreased (less stable). In contrast, the critical crack length increased with increasing precipitation. With regard to spatial simulations of snow stability, the high sensitivity on precipitation suggests that accurate precipitation patterns are necessary to obtain realistic snow stability patterns. With regard to spatial simulations of snow stability, the high sensitivity on precipitation suggests that accurate precipitation patterns are necessary to obtain realistic snow stability patterns.


2017 ◽  
Vol 11 (1) ◽  
pp. 217-228 ◽  
Author(s):  
Johan Gaume ◽  
Alec van Herwijnen ◽  
Guillaume Chambon ◽  
Nander Wever ◽  
Jürg Schweizer

Abstract. The failure of a weak snow layer buried below cohesive slab layers is a necessary, but insufficient, condition for the release of a dry-snow slab avalanche. The size of the crack in the weak layer must also exceed a critical length to propagate across a slope. In contrast to pioneering shear-based approaches, recent developments account for weak layer collapse and allow for better explaining typical observations of remote triggering from low-angle terrain. However, these new models predict a critical length for crack propagation that is almost independent of slope angle, a rather surprising and counterintuitive result. Based on discrete element simulations we propose a new analytical expression for the critical crack length. This new model reconciles past approaches by considering for the first time the complex interplay between slab elasticity and the mechanical behavior of the weak layer including its structural collapse. The crack begins to propagate when the stress induced by slab loading and deformation at the crack tip exceeds the limit given by the failure envelope of the weak layer. The model can reproduce crack propagation on low-angle terrain and the decrease in critical length with increasing slope angle as modeled in numerical experiments. The good agreement of our new model with extensive field data and the ease of implementation in the snow cover model SNOWPACK opens a promising prospect for improving avalanche forecasting.


2020 ◽  
Author(s):  
Stephanie Mayer ◽  
Alec van Herwijnen ◽  
Mathias Bavay ◽  
Bettina Richter ◽  
Jürg Schweizer

<p>Numerical snow cover models enable simulating present or future snow stratigraphy based on meteorological input data from automatic weather stations, numerical weather prediction or climate models. To assess avalanche danger for short-term forecasts or with respect to long-term trends induced by a warming climate, the modeled vertical layering of the snowpack has to be interpreted in terms of mechanical instability. In recent years, improvements in our understanding of dry-snow slab avalanche formation have led to the introduction of new metrics describing the fracture processes leading to avalanche release. Even though these instability metrics have been implemented into the detailed snow cover model SNOWPACK, validated threshold values that discriminate rather stable from rather unstable snow conditions are not readily available. To overcome this issue, we compared a comprehensive dataset of almost 600 manual snow profiles with simulations. The manual profiles were observed in the region of Davos over 17 different winters and include stability tests such as the Rutschblock test as well as observations of signs of instability. To simulate snow stratigraphy at the locations of the manual profiles, we obtained meteorological input data by interpolating measurements from a network of automatic weather stations. By matching simulated snow layers with the layers from traditional snow profiles, we established a method to detect potential weak layers in the simulated profiles and determine the degree of instability. To this end, thresholds for failure initiation (skier stability index) and crack propagation criteria (critical crack length) were calibrated using the observed stability test results and signs of instability incorporated in the manual observations. The resulting instability criteria are an important step towards exploiting numerical snow cover models for snow instability assessment.</p>


2019 ◽  
Vol 13 (12) ◽  
pp. 3353-3366 ◽  
Author(s):  
Bettina Richter ◽  
Jürg Schweizer ◽  
Mathias W. Rotach ◽  
Alec van Herwijnen

Abstract. Observed snow stratigraphy and snow stability are of key importance for avalanche forecasting. Such observations are rare and snow cover models can improve the spatial and temporal resolution. To evaluate snow stability, failure initiation and crack propagation have to be considered. Recently, a new stability criterion relating to crack propagation, namely the critical crack length, was implemented into the snow cover model SNOWPACK. The critical crack length can also be measured in the field with a propagation saw test, which allows for an unambiguous comparison. To validate and improve the parameterization for the critical crack length, we used data from 3 years of field experiments performed close to two automatic weather stations above Davos, Switzerland. We monitored seven distinct weak layers and performed in total 157 propagation saw tests on a weekly basis. Comparing modeled to measured critical crack length showed some discrepancies stemming from model assumption. Hence, we replaced two variables of the original parameterization, namely the weak layer shear modulus and thickness, with a fit factor depending on weak layer density and grain size. With these adjustments, the normalized root-mean-square error between modeled and observed critical crack lengths decreased from 1.80 to 0.28. As the improved parameterization accounts for grain size, values of critical crack lengths for snow layers consisting of small grains, which in general are not weak layers, become larger. In turn, critical weak layers appear more prominently in the vertical profile of critical crack length simulated with SNOWPACK. Hence, minimal values in modeled critical crack length better match observed weak layers. The improved parameterization of critical crack length may be useful for both weak layer detection in simulated snow stratigraphy and also providing more realistic snow stability information – and hence may improve avalanche forecasting.


2021 ◽  
Author(s):  
Bertil Trottet ◽  
Ron Simenhois ◽  
Gregoire Bobillier ◽  
Alec van Herwijnen ◽  
Chenfanfu Jiang ◽  
...  

<p>Snow slab avalanche release can be separated in four distinct phases : (i) failure initiation in a weak snow layer buried below a cohesive snow slab, (ii) the onset and, (iii) dynamic phase of crack propagation within the weak layer across the slope and (iv) the slab release. The highly porous character of the weak layer implies volumetric collapse during failure which leads to the closure of crack faces followed by the onset of frictional contact. To better understand the mechanisms of dynamic crack propagation, we performed numerical simulations, snow fracture experiments, and analyzed the release of full scale avalanches. Simulations of crack propagation are based on the Material Point Method (MPM) and finite strain elastoplasticity. Experiments consist of the so-called Propagation Saw Test (PST). Concerning full scale measurements, an algorithm is applied to detect changes in image pixel intensity induced by slab displacements. We report the existence of a transition from sub-Rayleigh anticrack to supershear crack propagation following the Burridge-Andrews mechanism. In detail, after reaching the critical crack length, self-propagation starts in a sub-Rayleigh regime and is driven by slab bending induced by weak layer collapse. If the slope angle is larger than a critical value, and if a so-called super critical crack length is reached, supershear crack propagation occurs. The corresponding critical angle may be lower than the weak layer friction angle due to the loss of frictional resistance during volumetric collapse. The sub-Rayleigh regime is driven by mixed mode anticrack propagation while the supershear regime corresponds to a pure mode II propagation with intersonic crack speeds (v: crack speed, c<sub>s</sub>: shear wave speed, c<sub>p</sub>: longitudinal wave speed, E: slab Young's modulus and ρ: slab density). This intersonic regime of crack propagation thus leads to pure tensile slab fractures initiating from the bottom of the slab as opposed to top initiations induced by slab bending in the sub-Rayleigh regime. Key ingredients for the existence of this transition are discussed such as the role played by friction angle, collapse height and slab secondary fractures. </p>


2016 ◽  
Vol 16 (3) ◽  
pp. 775-788 ◽  
Author(s):  
Fabiano Monti ◽  
Johan Gaume ◽  
Alec van Herwijnen ◽  
Jürg Schweizer

Abstract. The process of dry-snow slab avalanche formation can be divided into two phases: failure initiation and crack propagation. Several approaches tried to quantify slab avalanche release probability in terms of failure initiation based on shear stress and strength. Though it is known that both the properties of the weak layer and the slab play a major role in avalanche release, most previous approaches only considered slab properties in terms of slab depth, average density and skier penetration. For example, for the skier stability index, the additional stress (e.g. due to a skier) at the depth of the weak layer is calculated by assuming that the snow cover can be considered a semi-infinite, elastic, half-space. We suggest a new approach based on a simplification of the multi-layered elasticity theory in order to easily compute the additional stress due to a skier at the depth of the weak layer, taking into account the layering of the snow slab and the substratum. We first tested the proposed approach on simplified snow profiles, then on manually observed snow profiles including a stability test and, finally, on simulated snow profiles. Our simple approach reproduced the additional stress obtained by finite element simulations for the simplified profiles well – except that the sequence of layering in the slab cannot be replicated. Once implemented into the classical skier stability index and applied to manually observed snow profiles classified into different stability classes, the classification accuracy improved with the new approach. Finally, we implemented the refined skier stability index into the 1–D snow cover model SNOWPACK. The two study cases presented in this paper showed promising results even though further verification is still needed. In the future, we intend to implement the proposed approach for describing skier-induced stress within a multi-layered snowpack into more complex models which take into account not only failure initiation but also crack propagation.


2010 ◽  
Vol 24 (15n16) ◽  
pp. 3161-3166 ◽  
Author(s):  
LIYONG NI ◽  
CHAO LIU ◽  
CHUNGEN ZHOU

The durability and reliability of thermal barrier coatings(TBCs) have become a major concern of hot-section components due to lack of a reliable life prediction model. In this paper, it is found that the failure location of TBCs is at the TBC/TGO interface by a sequence of crack propagation and coalescence process. The critical crack length of failure samples is 8.8mm. The crack propagation rate is 3-10µm/cycle at the beginning and increases largely to 40µm/cycle near coating failure. A life prediction model based a simple fracture mechanics approach is proposed.


2019 ◽  
Vol 29 (5) ◽  
pp. 810-830 ◽  
Author(s):  
Dao-Hang Li ◽  
De-Guang Shang ◽  
Jin Cui ◽  
Luo-Jin Li ◽  
Ling-Wan Wang ◽  
...  

A fatigue–oxidation–creep damage model that can take into account the effect of multiaxial cyclic feature on the damage mechanism is proposed under axial-torsional thermo-mechanical fatigue loading. In the proposed model, the effects of non-proportional additional hardening on fatigue, oxidation, and creep damages are considered, and the variation of oxidation damage under different high temperature loading conditions is also described. Moreover, the intergranular creep damage needs to be equivalent to the transgranular damage before accumulating with the fatigue and oxidation damages. The fatigue, oxidation, and creep damages can be expressed as the fractions of fatigue life, critical crack length, and creep rupture time, respectively, which allows the linear accumulation of different types of damages on the basis of life fraction rule. In addition, the proposed model is validated by various fatigue experimental results, including uniaxial thermo-mechanical fatigue, axial-torsional thermo-mechanical fatigue, and isothermal axial-torsional fatigue under proportional and non-proportional loadings. The results showed that the errors are within a factor of 2.


1998 ◽  
Vol 26 ◽  
pp. 307-312 ◽  
Author(s):  
H. Conway

Measurements and observations by others indicate that a potential slab avalanche consists of a relatively cohesive slab of snow overlying a thin weak layer that coniains flaws where locally the shear stress from the overburden is not fully supported. Under favorable conditions, snow will shear strain-soften, which provides the basis for applying a slip-weakening model to examine the size of flaw needed to initiate sub-critical crack propagation along the weak layer. Using typical values for snow properties, the model predicts sub-critical crack growth can initiate from a relatively small flaw well before the shear stress from the overburden approaches the peak shear strength at tin-bed. The occurrence of small flaws or imperfections in the basal layer would explain field measurements which usually indicate that avalanching occurs before the applied shear stress exceeds the shear strength at the basal layer.Widespread slab-avalanche activity often increases significantly soon after the onset of rain on new snow. Measurements of temperature and mechanical properties show that only the upper 0.15 m or less of the slab has been altered at the time of avalanching; alterations at the sliding layer have not yet been detected. Results from the slip-weakening model indicate that the rain-induced alterations would reduce the size of flaw needed to initiate sub-critical crack growth by 10–20%. The observations and model results show clearly the importance of the slab properties; it is evident that both the slab and the weak layer act together to control slope stability. A further implication is that the stability of freshly deposited snow is often close to critical, because a relatively small surface perturbation is often sufficient to cause avalanching. This is not surprising, because it is well known from field observations that new snow on slopes should be treated with caution.


Sign in / Sign up

Export Citation Format

Share Document