1. Introduction and Setting of The Problem
Each winter’s contest between stresses and strengths in an inclined snow-pack leads to many avalanches in moutainous areas with a serious death toll and extensive damage to property. Therefore, the fracture mechanism of avalanches and the mechanical stability of the snow cover are very important for avalanche forecasting and other applications, e.g. estimating the time for closing roads or ski runs, for choosing the correct route for skiers, or for artificial triggering of avalanches.
The fracture of a slab avalanche (Fig. 1) is a multistage, progressive, and very fast-running process, details of which are still clouded with many secrets. Among the unknowns are the locally and temporally varying tensile, compressional, and shear stresses, and also the resisting strengths, as they cannot be measured directly because of the specific properties of the snow and the severe danger involved. Because of the topographic relief and the specific build-up of the snow cover in the course of a winter, its stability across a potential avalanche area is variable both in space and time. Therefore, idealized models, limiting considerations, and separate investigations of each stage of this complex process are needed. The present simplified model deals only with the fracture process itself in a geotechnical manner. The preceding avalanche “maturing” will be touched on briefly, where necessary.
2. Stages of Fracture
In response to the variety of forces along its boundaries (Figs 1 and 2), a snow-slab avalanche fractures in several stages.
2.1. Shear fracture
The shear fracture along the shear interface, which is more or less parallel to the surface, is the most important as regards the area, and is therefore considered to be the primary fracture by most experts. The development leading up to the failure of the entire shear surface is characterized by the large deformability of snow under shear stresses, and the resulting dissipative processes. Therefore, the strength depends heavily (i.e. inversely) on the strain-rate. Following Reference GublerGubler (1977, 1988), the initial fracture along a limited weak zone of the shear surface requires that a minimal strain-rate be exceeded such that the ductile initial strength may be overcome. Further propagation of the fracture then occurs as a brittle fracture quasi-elastically at high deformation rates. In the course of this progression, therefore, two quite different strengths play a role, one being the higher initial strength (at low strain-rate), and the other the lower brittle strength (at a very fast deformation speed in the course of the crack’s expansion). The question to what extent the actual shear strain is important – with granular media, as we know, a certain shearing displacement is required to activate the shear strength – appears to be unresolved.
Weakening and failure of the shear surface may be caused by:
Increasing stresses and consequent higher strain-rates, caused by loading (snow precipitation, skiers), or redistribution of stresses.
Local strain softening of the snow caused by heterogeneous deformations.
Decreasing grain connections caused by TG metamorphism or melting processes.
Temperature rising with strength decreasing.
Collapse of a depth-hoar layer (especially in areas with a continental climate).
2.2. Tensile fracture
The tensile fracture at the crown of the snow-slab avalanche – the secondary fracture in most cases – is a cleavage fracture, which is determined by the rather brittle deformation characteristics of snow under tension. This fact gives many reasons for believing that here could be primary fracture (e.g.Reference Haefeli and Kingery Haefeli, 1963, 1967; Reference SommerfeldSommerfeld, 1969). As opposed to the shearing zone, which is often only a few millimetres thick, here the properties of the snow cover in its entire depth, consisting of layers with varying tensile strengths and deformability, come into play. Stress concentrations are built up in stiffer layers, also causing the fracture process to be progressive in this case. This effect is superimposed by the fact that the crown’s height is frequently not uniform.
2.3. Flank fracture
During flank fracture, which is the third phase, shear-fracture mechanisms and tensile-fracture mechanisms overlap. As a consequence, we find different aspect angles between the lateral flanks and the down-slope direction (Fig. 2).
2.4. Compressive fracture
Finally, following considerable compression of the snow at the stauchwall, compressive fracture will occur as shearing across the whole snow cover. Here again, what matters are the strength properties of the various snow layers involved.
3. Stress Redistributions and Propagation of The Failure
3.1. Primary fracture of a thin shear layer
Before the initial fracture (see section 2.1) occurs, there is an equilibrium state between the driving stresses and the shear, tensile, and compressive strengths, of which different fractions are activated. When the shear strength ßs drops by Δβs or ceases altogether in the initial fracture zone (of slope-parallel length Δx), the fraction which is no longer transmitted down into the deeper snow layers (for each metre width of the slab)
must be supported by larger circumferential forces in the snow layers above the shear surface. If these are tensile forces, the tensile stress σx , on average, is increased by
where D is thickness of snow-pack (normal to the surface).
This leads to increased stresses near the initial fracture zone, which may either be transmitted down into deeper snow layers when the shear strength in the surrounding zone is high enough (in which case the shear fracture is stopped), or they may also cause a shear fracture here. Because of this propagating shear fracture, the tensile stresses accumulate to the point where the tensile strength of the snow-pack is reached and the tensile fracture is initiated. It will be shown later that the developing slab avalanche is greater the stronger the snow layers are above the shear surface.
3.2.Collapse of a thick layer (e.g. depth hoar)
Depth hoar at the bottom has a texture with coarse pores, large snow crystals, and few points of contact for transmitting forces. At higher stress, collapse of the texture may occur due to a local compressive failure, immediately preceding the avalanche fracture itself (seeReference Bucher Bucher, 1948; or Reference Bradley, Brown and WilliamsBradley and others, 1977). In the case of depth-hoar collapse, which causes compression of the depth hoar along with a decrease of the shear strength, we may imagine a variety of complex failure mechanisms:
(a) The collapse of texture leads to a loss of cohesion (mostly small to start with) in the basal layer, which cannot be compensated for by increased friction due to denser packing of the snow grains. Shear strength decreases and, consequently, the shear fracture widens across the basal layer much as in the case of thin shear layers, until circumferential fractures set in.
(b) The collapse of texture is accompanied by compression. If there is a stiff upper layer, the compression leads to a state in which the upper layer locally becomes detached from the basal layer, being subjected to bending forces. The area of bending widens along with the lateral propagation of the collapse. A certain part of the upper layer may remain hollow without breaking, but its size depends on its bending tensile strength, its thickness, the load, and the support conditions. Before the simultaneous failure of the shear support leads to fracture propagation, as in case (a), the bending strength is exceeded and cracks across the upper layer interrupt the lateral and longitudinal force transmission. Thus, the circumferential support, being already severely stressed by the failure of the shear support, collapses entirely. Fracture of the slab avalanche becomes unavoidable.
Preliminary investigations of a circular slab supported along the edges have shown that the bending tensile stresses at the snow surface along the edges due to the slab’s own weight and the effect of clamping are probably larger than those in the middle of the lower side. Therefore, the bending tensile cracks begin at the snow surface. No values have been found in the literature for the bending strength of snow. An investigation of the fracture mechanics of these processes is being planned for the future.
4. Stability Against Shear Fracture
As a rule, especially for practical purposes, the slope stability is investigated solely with respect to the possible shear fracture in the neutral zone, disregarding the boundary conditions along the entire snow-slab avalanche, as well as the fracture’s progressive character. This means that, for the potential shear surface, a stability factor is formed (really a safety factor in the technical sense). s is the ratio of the shear strength ßs to driving shear stress parallel to the slope τzx (Reference RochRoch, 1966a;Reference FÖhn Föhn, 1981, 1987; Reference Conway and AbrahamsonConway and Abrahamson, 1984):
Normal and shear stresses are obtained from the dead load in the special case after Rankine (activated friction angle equals slope angle):
Here, we use the notation: ρ is mean density of the snow-pack, D is thickness of the snow-pack (normal to the surface), ψ is slope inclination, and g = 10 m/s2.
If there is an additional load (component parallel to the slope Δτzx ), caused by snowfall, snowdrift, skiers, snowmobiles, or avalanche blasting, the stability changes as follows:
Usually ßs is determined from in situ shear-frame tests, the values therefore being rather inaccurate. Also, ßs may change, e.g. due to loads which increase σz :
In Equation (equ.6), a possible change in the strength Δβs has been included. The Mohr diagram (Fig. 3) is a graphical
representation of Equations (5) and (6) for varying loads. The shear-strength line (Coulomb-Mohr) is taken to be invariable.
Reference FÖhnFöhn (1981, 1987) investigated mainly the amount of additional tangential stress (load case Β in Figure 3) and the changes in stability, particularly in the case of avalanche blasting, but at the same time disregarding the changes of strength due to the increased normal stresses or to dynamic effects in the snow-pack.
As proof of the fact that the method based on the stability factor is unsatisfactory, we take the results that Roch (1966a), as well as Conway and Abrahamson (1984), with their improved experimental technique obtained shear stabilities up to 4 (with a very large scatter) even when the avalanche had fractured. This can partly be explained by the fact that strengths are particularly variable across short distances. A similarly large variance is, no doubt, true for the stresses. Föhn (1987) stated: “ …, the method reveals still too many deficiencies.” Apparently, locally restricted considerations of strength and stress are insufficient for describing the overall stability of a slope. Therefore, in the following, we will consider the total potential avalanche.
5. Supporting Forces Acting On A Snow Slab
5.1. The “standard” snow-slab avalanche
Although we are aware of the complex and heterogeneous conditions, we are trying to estimate the balance of forces acting on a snow slab by adding up individual contributions. By this means, we are neglecting the sequence of the above-described fracture phases and their exact course, especially progressive processes and stress redistributions. We will try to model important factors such as the suspension of the slab at the crown, its support at the stauchwall, and boundary influences from the flanks in relation to the forces along the shear surface by using a simplified “standard avalanche”. The following forces act along the down-slope direction (Fig. 2): Τ is the driving component of the total weight G, Fz is the tension force at the crown, FD is the compression force at the stauchwall, FF is flank force, Fs is shear force along the shear surface.
According to Reference PerlaPerla (1977, 1980) and Föhn (1981), the standard avalanche is characterized by the following values: Β = 50 m (width of avalanche), L = 50 m (length of avalanche), D = 0.7 m (thickness of avalanche), ψ = 38° (slope inclination), and ρ = 220 kg/m3 (density of snow (valid also below the sliding surface)).
The following active forces may be derived:
Following Reference MellorMellor (1975, figs 17 and 18), we obtain limiting values of the tensile, compressive, and shear strengths for snow density of 220 kg/m3. These values are shown in Table I below. This table also contains the
computation of minimal and maximal supporting forces of the standard snow slab.
The equation for the uniaxial tensile strength βz is accepted, being justified by the mechanism of cleavage fracture. The simplified computation for the stauchwall with the uniaxial compressive strength ßD has been checked by investigating a shear fracture across the snow-pack, similar to passive Earth pressure. Following Reference ZieglerZiegler (1963) and applying the Tresca criterion for fracture
we compute the slope-parallel supporting force at the stauchwall as
with limiting values 193 < FD < 774 kN. Compared to Table I with 175 < FD < 700 kN (calculated with ßD ), there is a difference of only roughly 10%. The complete solution for the inclined half-space including the internal friction is too cumbersome in the present context and changes the force FD only slightly (Reference KupperKupper, 1967). The Tresca criterion was also applied to the cohesion ß F of the flanks.
If, in place of the Tresca criterion (as in Table I: ßs = c = constant), we use the Coulomb-Mohr equation for obtaining the shear strength of the shear surface, we obtain for the chosen standard avalanche the limits for ßs as 0.4 and 2.6 kPa, respectively, following Föhn’s (1981) analyses based on papers by Reference RochRoch (1966a, b). The limits of strengths shown in Table I may, therefore, be considered as the right order of magnitude.
The computed minimal and maximal forces from Table I may now be related to different avalanche scenarios, such that the restraining shear forces Fs may combine in various ways with the supporting circumferential forces Fu
We distinguish between an internally strong snow slab (e.g. formed by wind compaction) and a snow slab of low strength consisting of loose snow, and further between a solid substrate well connected with the snow slab, and a weak shear surface (e.g. hoar frost being snowed in). The various possible combinations of suporting forces and the respective integral stabilities:
are put together in Table II.
Taking the combination with minimal shear force in the second-to-last line of this table, a case which may occur with heavy snowfall, we see from s = 0.72 that fracturing of the snow-slab avalanche is unavoidable unless
* FU = Fz + FF .
additional supporting forces can be activated, for example vaulting with additional lateral support or support at protrusions in the relief.
The last line of Table II refers to a combination with minimal shear force, whereby the pressure force at the stauchwall has been disregarded because much deformation is needed to activate it. Here, despite the large tensile and flank forces, integral stability barely exceeds 1. This reinforces the overwhelming importance of the shear strength.
It has frequently been assumed that the influence of the boundaries is negligible due to their small contribution to the total area; for the standard avalanche, the circumferential surface is only 5.6% of the shear surface, but for many real avalanches it is even less. The above conclusion, however, may be quite erroneous; taking the ratio FU /Fs as a measure of the boundary contribution to the support of a snow-slab avalanche, this ratio (from Table II) is indeed only 4–15% at high stability (which tends to prevent avalanche fracture), but increases rapidly to 36–154% in the critical cases of low stability which are of particular interest.
5.2. The parameter study
Generalizing the computations (Tables I and II), a parameter study has been undertaken using the following constants and variables:
Constants
Density ρ = 220 kg/m3 (see note (a))
Circumferential strengths ßz , ßD , ßF … minima and maxima for each (from Table I)
Shear strength βs = 0.5 kN/m2 … minimum (see note (b))
Slope inclination ψ = 38° (see note (c))
Variables
Width Β = length L (10 m ⩽ Β ⩽ 110 m)
Thickness D (0.4 m ⩽ D ⩽ 1.5 m)
Note (a). In a comparative investigation for higher density and disproportionately higher strengths (see Mellor, 1975), stabilities 5 worked out to be many times higher; these values are therefore not relevant. For lower densities and
consequently lower strengths, applicable particularly to new-snow avalanches, there are hardly any physical parameters in the literature.
Note (b). As shown already in Table II, no avalanche releases may be expected with high shear strengths.
Note (c). Varying the slope inclination angle by ±6°, and retaining all other parameters and Η = D/cos ψ = constant, resulted in a change in stability s of only roughly ±10% for the “standard avalanche”.
Overall stability for Β = L is obtained from the forces according to Table I and Equation (equ.13):
Varying the width or the shear surface area (B = (As )½), and keeping the slab’s thickness D = 0.7 m constant, gives the results shown in Figure 4. It can be seen from the approximately hyperbolic curve shape for s as well as for FU /Fs that, with increasing size of the slab, the circumferential forces’ supporting role weakens considerably; their contribution, however, remains above 70% of the shearing force in this case with a weak shear surface. Even so, the circumferential forces can prevent fracture only up to a certain critical width , because beyond this point-force transmission becomes impossible due to the preceding circumferential fracture. BK increases considerably with increasing circumferential strengths (see Figs 4 and 6). Stronger, wind-blown slab avalanches with proportionately larger circumferential forces (90–300% of Fs ) may, therefore, fracture as a whole large area, even when the shear surface is weak.
The influence of varying depths on s retaining the width constant, is shown in Figure 5. Stability decreases with increasing depth, this decrease being more pronounced the larger and weaker the slab. The distance between the two lines for a strong and a weak slab, respectively, is a measure of how the circumferential forces matter more for smaller slabs. Upon reaching the critical depth DK the slab avalanche fractures; this point is reached sooner both with weak and with large slabs.
Setting (s = 1) ., we can compute DK for constant Β and BK for constant D from Equation (equ.14) as follows:
They are shown in one diagram (Fig. 6).
The hyperbolic shape of the curves means that, with the assumed low shear strength and at small depths, a small increase in depth leads to a very noticeable drop in the maximum possible slab size or, alternatively, in the region with (s = 1) . in an inclined snow-pack. This region is smaller the lower the circumferential strength (line b, Fig. 6). This favors the spontaneous release of new-snow avalanches. Conversely, hardening of the snow enables larger slabs to remain stable (line a, Fig. 6). When, however, the avalanche is released, e.g. by a skier’s forces, it has a larger area.
On the other hand, if a slab avalanche’s size is constrained by the topographic relief, a certain depth (the critical depth DK ) cannot be exceeded.
6. Conclusions
The first result of the parameter study carried out is confirmation of the overwhelming importance of the shear strength of surfaces parallel to the slope for maintaining the snow-pack’s stability. Particularly low strengths should be the preferred objective of intensive future studies, even though they are the most difficult to control in experiments.
The second result, however, is that the boundary conditions play an important role in the process of stress redistribution, particularly when the resisting shear forces are weak. Even with low circumferential strengths, circumferential forces make a considerable contribution to overall stability of a slab; this fraction, however, drops rapidly with increasing avalanche size. Boundary conditions must be increasingly included in theoretical and practical investigations, although, naturally, this will not facilitate the solution of problems connected with avalanche fracture.
Taking note of the great variety of snow types, and of each winter’s different time sequence of the weather, both extensive observations and measurements of snow-slab avalanches (geometry, snow data, most of all density, stratification, strengths), as well as further mechanical models of the force interplay, the stresses, deformations, and failure mechanisms will be needed in order to understand better the important and complex problem of snow-slab avalanche release, particularly regarding the processes initiating fracture.