This numerical study investigates how different types of snow avalanches behave, how key factors affect their dynamics and flow regime transitions, and what are the underpinning rules. According to the unified trends obtained from the simulations, we are able to quantify the complex interplay between bed friction, slope geometry and snow mechanical properties (cohesion and friction) on the maximum velocity, runout distance and deposit height of the avalanches.

Based on DEM simulations we developed a new model for the onset of crack propagation in snow slab avalanche release. The model reconciles past approaches by considering the complex interplay between slab elasticity and the mechanical behavior of the weak layer including its structural collapse. The model agrees with extensive field data and can reproduce crack propagation on low-angle terrain and the decrease in critical crack length with increasing slope angle observed in numerical experiments.

We propose 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. The method was tested on simplified snow profiles, then on manually observed snow profiles including a stability test and, finally, on simulated snow profiles, thereby showing the promise of our approach.

We proposed a new approach to characterize the dynamic phase of crack propagation in weak snowpack layers as well as fracture arrest propensity by means of numerical "propagation saw test" simulations based on the discrete element method. Crack propagation speed and distance before fracture arrest were derived from the simulations for different snowpack configurations and mechanical properties. Numerical and experimental results were compared and the mechanical processes at play were discussed.

Slab tensile failure propensity is examined using a mechanical--statistical model of the slab–-weak layer (WL) system based on the finite element method. This model accounts for WL heterogeneity, stress redistribution by elasticity of the slab and the slab possible tensile failure. For realistic values of the parameters, the tensile failure propensity is mainly driven by slab properties. Hard and thick snow slabs are more prone to wide–scale crack propagation and thus lead to larger avalanches.