Numerical investigation of the mixed-mode failure of snow

The failure of a weak snow layer underlying a cohesive slab is the primary step in the release process of a dry snow slab avalanche. The complex and heterogeneous microstructure of snow limits our understanding of failure initiation inside the weak layer, especially under mixed-mode shear-compression loading. Further complication arises from the dependence of snow strength on the loading rate induced by the balance between bond breaking and bond formation (sintering) during the failure process. Here, we use the discrete element method to investigate the influence of mixed-mode loading and fast sintering on the failure of a weak layer generated using cohesive ballistic deposition. Both fast and slow loading simulations resulted in a mixed-mode failure envelope in good agreement with laboratory experiments. We show that the number of broken bonds at failure and the weak layer strength significantly decreases with increasing loading angle, regardless of the loading rate. While the influence of loading rate appears negligible in shear-dominant loading (for loading angles above $$ ), simulations suggest a significant increase in the weak layer strength at low loading angles and low loading rates, characteristic of natural avalanches, due to the presence of an active sintering mechanism.