Avalanche Front Research Articles

Phenomenology of Avalanche Recordings From Distributed Acoustic Sensing

AbstractAvalanches and other hazardous mass movements pose a danger to the population and critical infrastructure in alpine areas. Hence, understanding and continuously monitoring mass movements are crucial to mitigate their risk. We propose to use Distributed Acoustic Sensing (DAS) to measure strain rate along a fiber‐optic cable to characterize ground deformation induced by avalanches. We recorded 12 snow avalanches of various dimensions at the Vallée de la Sionne test site in Switzerland, utilizing existing fiber‐optic infrastructure and a DAS interrogation unit during the winter 2020/2021. By training a Bayesian Gaussian Mixture Model, we automatically characterize and classify avalanche‐induced ground deformations using physical properties extracted from the frequency‐wavenumber and frequency‐velocity domain of the DAS recordings. The resulting model can estimate the probability of avalanches in the DAS data and is able to differentiate between the avalanche‐generated seismic near‐field, the seismo‐acoustic far‐field, and the mass movement propagating on top of the fiber. By analyzing the mass‐movement propagation signals, we are able to identify group velocity packages within an avalanche that propagate faster than the phase velocity of the avalanche front, indicating complex internal structures. Importantly, we show that the seismo‐acoustic far‐field can be detected before the avalanche reaches the fiber‐optic array, highlighting DAS as a potential research and early warning tool for hazardous mass movements.

Detrainment and braking of snow avalanches interacting with forests

Abstract. Mountain forests provide natural protection against avalanches. They can both prevent avalanche formation in release zones and reduce avalanche mobility in runout areas. Although the braking effect of forests has been previously explored through global statistical analyses on documented avalanches, little is known about the mechanism of snow detrainment in forests for small and medium avalanches. In this study, we investigate the detrainment and braking of snow avalanches in forested terrain, by performing three-dimensional simulations using the material point method (MPM) and a large-strain elastoplastic snow constitutive model based on critical state soil mechanics. First, the snow internal friction is evaluated using existing field measurements based on the detrainment mass, showing the feasibility of the numerical framework and offering a reference case for further exploration of different snow types. Then, we systematically investigate the influence of snow properties and forest parameters on avalanche characteristics. Our results suggest that for both the cold and warm snow parameterized in our simulations, the detrainment mass decreases with the square of the avalanche front velocity before it reaches a plateau value. Furthermore, the detrainment mass significantly depends on snow properties. It can be as much as 10 times larger for warm snow compared to cold snow. By examining the effect of forest configurations, it is found that forest density and tree diameter have cubic and square relations with the detrainment mass, respectively. The outcomes of this study may contribute to the development of improved formulations of avalanche–forest interaction models in popular operational simulation tools and thus improve hazard assessment for alpine geophysical mass flows in forested terrain.

Numerical modelling of turbulent geophysical flows using a hyperbolic shear shallow water model: Application to powder snow avalanches

In this work we apply a mathematical model developed by Teshukov (2007) to simulate turbulent powder snow avalanches. The two-parameter model describes the production of turbulent energy from shearing. This energy is associated with the formation of small and large vortices which provide avalanches with their distinctive billow and cleft-like structures. The model accurately predicts the concentration of translational kinetic energy at the avalanche front and likewise the formation of an almost stationary turbulent wake. The calculation of turbulent energy can be exploited to improve air-entrainment and turbulent drag models and therefore to improve engineering calculations of powder cloud height, speed and density, an important problem in snow avalanche mitigation. In present work we focus on the one-dimensional case. The governing equations are discretized with a finite volume second order Godunov-type scheme using HLLC Riemann solver. A good agreement between numerical solution of the new model and the photogrammetric measurements (height, length and frequency of billows, depths of clefts) is observed both at the front and tail of the avalanche for two different data sets.

The mechanical origin of snow avalanche dynamics and flow regime transitions

Abstract. Snow avalanches cause fatalities and economic damage. Key to their mitigation is the understanding of snow avalanche dynamics. This study investigates the dynamic behavior of snow avalanches, using the material point method (MPM) and an elastoplastic constitutive law for porous cohesive materials. By virtue of the hybrid Eulerian–Lagrangian nature of the MPM, we can handle processes involving large deformations, collisions and fractures. Meanwhile, the elastoplastic model enables us to capture the mixed-mode failure of snow, including tensile, shear and compressive failure. Using the proposed numerical approach, distinct behaviors of snow avalanches, from fluid-like to solid-like, are examined with varied snow mechanical properties. In particular, four flow regimes reported from real observations are identified, namely, cold dense, warm shear, warm plug and sliding slab regimes. Moreover, notable surges and roll waves are observed peculiarly for flows in transition from cold dense to warm shear regimes. Each of the flow regimes shows unique flow characteristics in terms of the evolution of the avalanche front, the free-surface shape, and the vertical velocity profile. We further explore the influence of slope geometry on the behavior of snow avalanches, including the effect of slope angle and path length on the maximum flow velocity, the runout angle and the deposit height. Unified trends are obtained between the normalized maximum flow velocity and the scaled runout angle as well as the scaled deposit height, reflecting analogous rules with different geometry conditions of the slope. It is found that the maximum flow velocity is mainly controlled by the friction between the bed and the flow, the geometry of the slope, and the snow properties. We reveal the crucial effect of both flow and deposition behaviors on the runout angle. Furthermore, our MPM modeling is calibrated and tested with simulations of real snow avalanches. The evolution of the avalanche front position and velocity from the MPM modeling shows reasonable agreement with the measurement data from the literature. The MPM approach serves as a novel and promising tool to offer systematic and quantitative analysis for mitigation of gravitational hazards like snow avalanches.

Estimation of Avalanche Development and Frontal Velocities Based on the Spectrogram of the Seismic Signals Generated at the Vallée de la Sionne Test Site

The changes in the seismic signals generated by avalanches recorded at three sites along a path at the Vallée de la Sionne (VdlS) experimental site are presented. We discuss and correlate the differences in the duration, signal amplitudes, and frequency content of the sections (Signal ONset (ON), Signal Body (SBO), and Signal TAil and Signal ENd STA-SEN) of the spectrograms with the evolution of the powder, transitional and wet snow avalanches along a path. The development of the avalanche front was quantified using the exponential function in time F (t) = K’ exp (β t) fitted to the shape of the signal ONset (SON section of the spectrogram. The speed of the avalanche front is contained in β. To this end, a new method was developed. The three seismic components were converted into one seismic component (FS), when expressing the vector in polar coordinates. We linked the theoretical function of the shape of the FS-SON section of the spectrogram to the numerical coefficients of its shape after considering the spectrogram as an image. This allowed us to obtain the coefficients K’ and β. For this purpose, the Hough Transform (HT) was applied to the image. The values of the resulting coefficients K’ and β are included in different ranges in accordance with the three types of avalanche. Curves created with these coefficients enable us to estimate the development of the different avalanche types along the path. Our results show the feasibility of classifying the type of avalanche through these coefficients. Average speeds of the avalanches approaching the recording sites were estimated. The speed values of wet and transitional avalanches are consistent with those derived from GEODAR (GEOphysical Doppler radAR) measurements, when available. The absence of agreement in the speed values obtained from seismic signals and GEODAR measurements for powder snow avalanches indicates, for this type of avalanche, a different source of the measured signal. Hence, the use of the two measuring systems proves to be complementary.

A Review of Russian Snow Avalanche Models—From Analytical Solutions to Novel 3D Models

The article is a review of mathematical models of snow avalanches that have been proposed since the middle of the 20th century and are still in use. The main attention is paid to the work of researchers from the Soviet Union and Russia, since many of their works were published only in Russian and are not widely available. Mathematical models of various levels of complexity for avalanches of various types—from dense to powder-snow avalanches—are discussed. Analytical solutions including formulas for the avalanche front speed are described. The results of simulations of the movement of avalanches are given that were used to create avalanche hazard maps. The last part of the article is devoted to constructing models of a new type, in which avalanches are considered as laminar or turbulent flows of non-Newtonian fluids, using the full (not depth-averaged) equations of continuum mechanics. The results of a numerical study of the effect of non-Newtonian rheology and mass entrainment on the avalanche dynamics are presented.

Hill Slope Falls and Long-Runout Rockslides under Large-Scale Underground Blasting

The collected data on long-runout rockslides induced be large-scale underground blasts are presented. The critical conditions of the rockslides are determined as a function of the dynamic impact, rockfall volume, as well as the slope height and gradient. It is shown that the main condition for the mobile rockslides is their volume exceeding 106-108 m3. A collateral effect for rockslides in case of smaller volume of 106-107 m3 is the presence of an acceleration phase at the slope toe and the absence of a counter slope. The relations of the avalanche front range, rockslide volume and fall height are given. It is mentioned that the avalanche range is associated with their composition, in particular, rock-and-ice or rock-and-snow avalanches feature increased mobility starting from their volumes of 105 m3. A brief review of theoretical and modeling studies intended to explain the reduced friction in large avalanches is provided.

On the velocity of turbidity currents over moderate slopes

In the present article we consider the problem of underwater avalanches propagating over moderate slopes. The main goal of our work is to investigate the avalanche front velocity selection mechanism when it propagates downwards. In particular, we show that the front velocity does not depend univocally on the mass of sediments. This phenomenon is investigated and explained in our study. Moreover, we derive from the first principles a depth-averaged model. Then, we assume that sediments are uniformly distributed along the slope. In this case, they can be entrained into the flow head and a self-sustained regime can be established. One of the main findings of our study is that the avalanche front velocity is not unique due to a hysteresis phenomenon. We attempt to explain this phenomenon using dynamical systems considerations.

The Intermittency Regions of Powder Snow Avalanches

Powder snow avalanches are typically composed of several regions characterized by different flow regimes. These include a turbulent suspension cloud of fine particles, a dense basal flow, and an intermittency frontal region, which is characterized by large fluctuations in impact pressure, air pressure, velocity, and density, but whose origin remains unknown. In order to describe the physical processes governing the intermittency region, we present data from four large powder snow avalanches measured at the Vallée de la Sionne test site in Switzerland, which show that the intermittency is caused by mesoscale coherent structures. These structures have a length of 3–14 m and a height of 10 m or more. The structures can have velocities as much as 60% larger than the avalanche front speed and are characterized by an air/particle mixture whose average density can be as high as 20 kg/m3. This average density increases the drag on large granules by a factor of up to 20 compared to pure air, so that each structure can maintain denser snow clusters and single snow granules in suspension for several seconds. The intermittency region has importance for the dynamics of an avalanche, as it provides an efficient mechanism for moving snow from the dense layer to the powder cloud, but also for risk assessment, as it can cause large forces at large heights above the basal dense layer.

Breaking size-segregation waves and mobility feedback in dense granular avalanches

Through experiments and discrete particle method (DPM) simulations we present evidence for the existence of a recirculating structure, that exists near the front of dense granular avalanches, and is known as a breaking size-segregation (BSS) wave. This is achieved through the study of three-dimensional bidisperse granular flows in a moving-bed channel. Particle-size segregation gives rise to the formation of a large-particle-rich front and a small-particle-rich tail with a BSS wave positioned between the tail and front. We experimentally resolve the structure of the BSS wave using refractive-index matched scanning and find that it is qualitatively similar to the structure observed in DPM simulations. Our analysis demonstrates a relation between the concentration of small particles in the flow and the amount of basal slip, in which the structure of the BSS wave plays a key role. This leads to a feedback between the mean bulk flow velocity and the process of particle-size segregation. Ultimately, these findings shed new light on the recirculation of large and small grains near avalanche fronts and the effects of this behaviour on the mobility of the bulk flow.

Intermittency between avalanche regimes on grain piles.

We experimentally investigate discrete avalanches of grains, driven by a low inflow rate, on an erodible pile in a channel. We observe intermittency between one regime, in which avalanches are quasiperiodic and system spanning, and another, in which they pass at irregular intervals and have a power-law size distribution. Observations are robust to changes of inflow rate and grain type and require no tuning of external parameters. We demonstrate that the state of the pile's surface determines whether avalanche fronts propagate to the end of the channel or stop partway down, and we introduce a toy model for the latter case that reproduces the observed power-law size distribution. We suggest direct applications to avalanches of pharmaceutical and geophysical grains, and the possibility of reconciling the "self-organized criticality" predicted by several authors with the hysteretic behavior described by others.

An analysis of the entrainment effect of dry debris avalanches on loose bed materials.

Substrate entrainment can greatly influence the mass movement process of a debris avalanche because it can enlarge the landslide volume and change the motion characteristics of the sliding masses. To study the interaction between debris avalanches and erodible substrate, physical modeling experiments varying in the mass of granular flow and substrate thickness were performed. The experimental results show that both the entrained materials and the maximum erosion depth are increased with increasing mass of the debris avalanche and decreasing substrate thickness. During the experiment, several tests were recorded using a high-speed digital camera with a frequency of 500 frames per second, so that the process of entrainment could be clearly observed. Combined with the experiment result and results of previous studies from predecessors, the entrainment mechanism during debris avalanches are analyzed and discussed. The entrainment effect of the sliding masses on the loose bed materials include basal abrasion and impact erosion of the avalanche front, the latter of which can contribute to the former by failing or yielding the erodible bed.

On the front shape of an inertial granular flow down a rough incline

Granular material flowing on complex topographies are ubiquitous in industrial and geophysical situations. In this paper, we study the small-scale experiment of a granular layer flowing on a rough incline. The shape of the granular front is solved analytically by using depth-averaged mass and momentum equations with a fractional expression for the frictional rheology μ(I), which is a generalization of Gray and Ancey [“Segregation, recirculation and deposition of coarse particles near two-dimensional avalanche fronts,” J. Fluid Mech. 629, 387 (2009)]. Unlike previous studies where a “plug flow dynamics” is assumed, a free shape factor α describing the vertical velocity profile is taken into account. The effect of inertia and shear rate on the front profile is evidenced through the introduction of the Froude number and the shape factor α. The analytical predictions are compared to experimental results published by Pouliquen [“On the shape of granular fronts down rough inclined planes,” Phys. Fluids 11, 1956 (1999)] and with our new experimental data obtained at higher Froude numbers. A good agreement between theory and experiments is found for α = 5/4, corresponding to a Bagnold-like velocity profile. However, we observe a systematic deviation near the head of the front where the height vanishes: the theory predicts a continuous precursor layer, while a grain-free region is observed experimentally. This suggests that the vertical velocity profile is not uniform inside the front, but the shape factor α tends to 1 near the head of the front. This raises questions about the vertical velocity profile in granular flows and about the expression of the rheological function μ(I) and its calibration from experimental data.

Stick–slip motion and high speed ejecta in granular avalanches detected through a multi-sensors flume

An oscillatory motion of dense, dry granular flows is experimentally observed by means of a 5.5m-long densely instrumented flume. An array of 35 sensing devices allowed recording kinematic and dynamic data of the entire flow with high spatial and temporal resolution. Dynamic data, in addition to the kinematics, provide valuable information to unveil stick–slip motion that is proposed as the triggering mechanism for the appearance of periodic changes in compaction of the granular material called density waves. Periodic reverse graded strata after deposition of the flow accounts for such density waves. Besides, internal interactions between different parts of the avalanche, lead to the ejection of material and a high peak of the avalanche front velocity — or an “abnormal acceleration” — at the break slope and ejection of rocks from the flow was systematically observed, evidencing internal interactions among different parts of the avalanche body while descending the inclined plane. The velocity peak registered at the break in slope is interpreted as complementary evidence of these internal collisions and a non perfect inelastic-collision simple model is assumed in order to explain the counter-intuitive acceleration observed during the ejection process. High speed ejection at the break in slope, stick–slip, and expulsion from the flow are robust, ubiquitous phenomena which do not strongly depend on the granulometric distribution used.

The full-scale avalanche test-site at Lautaret Pass (French Alps)

The full-scale avalanche test site at Lautaret Pass in the southern French Alps has been used by IRSTEA-Cemagref Research Institute since 1972. Over recent years, two avalanche paths have been used routinely to release avalanches and study avalanche dynamics and interactions between avalanches and obstacles. Avalanche flows are generally dense and dry, sometimes with a powder cloud on top. Main avalanche path no. 2 is dedicated to studies on avalanche dynamics. Within the flow of the avalanche, flow height and vertical profiles of pressure and velocity are measured along a 3.5m tripod. The snow volume released in the release zone is quantified by differential analysis of laser scanning measurements performed before and after triggering. High-speed positioning of the avalanche front along the track is carried out by terrestrial oblique photogrammetry. Above the dense layer, the upper layer of the avalanche is characterized by particle and air flux measurements. Avalanche path no. 1 is smaller in size and particularly well-suited to experiments on structures exposed to small to medium-size avalanches (<1000m3). A macroscopic sensor structure consisting of a one square-meter plate supported by a 3.5m high steel cantilever beam is fixed in the ground, facing the avalanche. Impact pressures are reconstructed from the beam deformations and avalanche velocity is measured by optical sensors. For these experimental devices dedicated to improving our understanding of avalanche physics, a national and international partnership has been developed over the years, including INSA de Lyon, CNRS and Université Joseph Fourier (France), Aalto University (Finland), Nagoya University (Japan), Boku University (Austria) and IGEMA (Bolivia).

Quantum ignition of deflagration in theFe8molecular magnet

We report spatially resolved, time-dependent, magnetization reversal measurements of an Fe$_{8}$ single molecular magnet using a microscopic Hall bar array. We found that under some conditions the molecules reverse their spin direction at a resonance field in the form of an avalanche. The avalanche front velocity is of the order of $1$m/sec and is sensitive to field gradients and sweep rates. We also measured the propagation velocity of a heat pulse and found that it is much slower than the avalanche velocity. We therefore conclude that in Fe$_{8}$, the avalanche front propagates without thermal assistance.

Two-dimensional radar imaging of flowing avalanches

Radar has emerged as an important tool in avalanche research. However, existing radar sensors suffer from coarse range resolution capabilities. This limits the usefulness of the data they collect in validating models of avalanche dynamics. This paper details the development of a frequency modulated continuous wave, phased array radar, and its associated signal processing, for non-invasive measurements of entire avalanche events. The radar outperforms existing avalanche radar sensors in terms of range resolution, and it provides cross-range resolution using a phased array receiver. The radar has been operating at the Vallée de la Sionne avalanche test site in Switzerland since the 2010 winter season. It has successfully gathered measurements of entire natural avalanche events. In this paper we show two-dimensional radar images of a naturally occurring avalanche, the first of their kind, which reveal movements of layers or particles of the flowing avalanche in unparalleled detail. Furthermore, the potential of the measured data is shown with tracking of avalanche fronts in two spatial dimensions. This marks an important step towards providing a library of high-quality avalanche measurements to improve our knowledge of avalanche dynamics.

Comparison of avalanche front velocity measurements and implications for avalanche models

Besides the runout distance of an avalanche information on avalanche intensity along the path is often required for hazard zoning or planning of mitigation measures. The avalanche (front) velocity is a common intensity measure as it can be linked to expected impact pressures. Furthermore, the velocity of an avalanche determines if the avalanche stays in its usual track or if the avalanche unexpectedly deviates and thus endangers areas believed to be safe. Therefore, a reliable prediction of the expected velocities is most important. However, many of the prevailing avalanche models either tend to underestimate velocities or they overestimate the runout distances.In this paper, several avalanche front velocity measurements from the Ryggfonn test site are presented. The measurements are derived from photo and/or video analyses of mainly dry-mixed avalanches. Some of these observations are partly supplemented by Doppler radar measurements.Additionally, the measurements from the Ryggfonn test site are compared with published velocity measurements from other locations to obtain a wider spectrum of avalanche conditions. By analyzing these velocities, constraints for possible rheological models of avalanche flows are obtained.The measurements of the presented avalanches, of which the most were relatively large to their path and can be classified as dry-mixed avalanches (i.e. they were partially fluidized and accompanied by a powder cloud), suggest a lower (apparent) velocity dependency of the bottom friction term than generally believed. The measurements show that the maximum front velocities of these avalanches scale with the square root of the total fall height of the avalanche front, that is, the maximum velocity Um∝Hsc.

Comparison of avalanche front velocity measurements: Supplementary energy considerations

A series of front velocity measurements were re-analyzed to derive estimates on the mean retarding acceleration of dry-mixed avalanches (i.e. dry-snow avalanches, which were partially fluidized and accompanied by a powder cloud) using energy considerations. The obtained estimates correspond well with those derived from around 320 solely runout observations. Based on these results it seems questionable whether classical avalanche models are adequate as basis for more physically-based ones.

Formation of levees and en-echelon shear planes during snow avalanche run-out

AbstractSnow avalanches often form levees and en-echelon shear planes in the run-out zone. We describe the formation of these depositional structures using a simple model that accounts for the role of granular fluctuations in avalanche motion. A mathematical feature of this model is the existence of a bifurcation saddle point, describing how granular fluctuations control the avalanche velocity in the runout zone. The saddle point discriminates between a flowing and stopping regime and defines the physical boundary between the flow and non-flow regions of the avalanche, i.e. the location of shear planes in the avalanche deposits. The formation of a shear plane depends on the interplay between terrain slope and avalanche mass flux, which varies from avalanche head to tail. Levees can form immediately at the avalanche front or, for steep slopes and low mass fluxes, at the avalanche tail. At ravine and gully shoulders the mass flux is restricted, thus initiating levee formation. We find that the levee lines are parallel to the flow direction when the mass flux is constant; en-echelon shear lines occur when the mass flux is decreasing. We test the model using several case studies where we have accurate laser scans of avalanche deposits. Our results suggest that avalanche flow parameters can be determined from simple levee measurements or, conversely, formation of levees and flow fingers can be predicted once the parameters governing the granular fluctuations are known.