Articles | Volume 18, issue 3
https://doi.org/10.5194/nhess-18-759-2018
https://doi.org/10.5194/nhess-18-759-2018
07 Mar 2018
 | 07 Mar 2018

Dynamic magnification factors for tree blow-down by powder snow avalanche air blasts

Perry Bartelt, Peter Bebi, Thomas Feistl, Othmar Buser, and Andrin Caviezel
Abstract

We study how short duration powder avalanche blasts can break and overturn tall trees. Tree blow-down is often used to back-calculate avalanche pressure and therefore constrain avalanche flow velocity and motion. We find that tall trees are susceptible to avalanche air blasts because the duration of the air blast is near to the period of vibration of tall trees, both in bending and root-plate overturning. Dynamic magnification factors for bending and overturning failures should therefore be considered when back-calculating avalanche impact pressures.

Dates
1 Introduction

In this paper we develop a simple method to determine the dynamic response of trees to impulsive loads. This is an important problem in natural hazards engineering where historical evidence of forest destruction or tree breakage is often used to evaluate the potential avalanche hazard. Any indication of forest damage is particularly valuable to avalanche engineers because it helps define the destructive reach of an extreme and infrequent event. Fallen tree stems delineate the spatial extent of an avalanche and create a natural vector field indicating the primary flow direction of the movement (Fig. 1). The age of the destroyed trees can be additionally used to link the historical observations to the avalanche return period (Reardon et al.2008; Schläppy et al.2014; Gadek et al.2017). In many cases observations of forest destruction are the only data the engineer has to quantify avalanche danger.

The problem with using evidence of tree destruction for avalanche mitigation planning is that a simple relationship between avalanche impact pressure and tree failure is difficult to establish. Tree-breaking depends on both the avalanche loading and tree strength. Trees fall if the bending stress exerted by the avalanche exceeds the bending strength of the tree stem (Johnson1987; Mattheck and Breloer1994; Peltola et al.1997, 1999) or if the applied torque overcomes the strength of the root-soil plate, leading to uprooting and overturning (Coutts1983; Mattheck and Breloer1994; Jonsson et al.2006). Both mechanisms depend on the local flow height of the avalanche. Recent observations by Feistl et al. (2015b) suggest that the magnitude of the avalanche impact pressure is strongly related to the avalanche flow regime. Although long recognised that dense flowing avalanches can easily break, overturn and uproot trees (Bartelt and Stöckli2001; Feistl et al.2015a), tree destruction by powder avalanche air blasts has received less attention. A mechanical understanding of how trees are blown-down by powder avalanche blasts would allow engineers to quantify powder avalanche pressures from case studies and historical records.

Here we develop a mechanical model to predict the natural frequency of trees subject to full-height air-blasts of powder snow or ice avalanches. We assume two deformation modes: stem bending and root-plate overturning, see Figs. 2 and 3. The ratio of the natural tree frequency to the frequency of the avalanche air-blast defines the dynamic magnification factor D (Clough and Penzien1975). This value is used to magnify the non-impulsive loadings D> 1 to account for the increase in stress under an impulsive load. The eigenfrequency of the tree is a function of the tree height, stiffness and mass distribution between the stem and branches. It therefore depends on forest age and tree species. We show that dynamic magnification factors for fully grown trees are large indicating that mature forests are especially vulnerable to powder snow avalanches. As we shall see, an error of up to 25 % can be made when back-calculating avalanche velocities. For example, an avalanche travelling at 35 m s−1 exerts the same pressure as an avalanche travelling at 50 m s−1 if the impulsive nature of the loading is considered. These are significant differences in hazard mitigation studies.

https://www.nat-hazards-earth-syst-sci.net/18/759/2018/nhess-18-759-2018-f01

Figure 1Tree breakage caused by the air blast of a powder avalanche, Zernez, Switzerland, 1999. The trees failed through a combination of bending and root-plate overturning. Photograph: Peter Bebi, SLF.

Download

2 Tree response to impulsive loading

Measurements on real avalanches reveal that the air-blast is intermittent and of short duration, lasting only a few seconds (Grigoryan et al.1982; Sukhanov and Kholobaeva1982; Sukhanov1982). When a powder avalanche hits a forest the ice-dust cloud is typically moving at velocities in excess of 50 m s−1 (similar to extreme wind gusts). The height of the cloud is equal, if not larger, than the height of the tree, i.e. H> 20 m. The pressure blast thus acts over the entire width and height of the tree, producing large bending moments in the stem and straining the root base plate. The impulsive character of the powder avalanche air-blast, however, magnifies the static stress state (Clough and Penzien1975). The fallen tree stems often point in the direction of the flow, indicating that the trees had little time to sway and react to blast and that the inertial effects are of considerable importance.

To calculate the dynamic magnification factor D we first make three simplifying assumptions. Firstly, the air blast can be expressed as a sine wave impulse with duration time t0. Moreover,

(1) F ( t ) = F 0 sin ω t ,

where ω is the circular frequency of the loading ω=πt0. The magnitude of the force F0 is as follows:

(2) F 0 = p 0 A = 1 2 c d ρ U max 2 A ,

where p0 is the amplitude of the avalanche pressure given by the density of the powder cloud ρ, the form drag coefficient of the tree cd and the maximum velocity of the blast Umax (Bozhinskiy and Losev1998; Feistl et al.2015a). The tree area over which the blast acts is denoted A, typically given by the tree height H and effective tree width W. Thus, if the cloud density and velocity are known as well as the tree geometry, the magnitude of the applied blast force F0 is given.

After the loading time t0, the tree vibrates freely with natural frequency ω. The natural frequency is found using the Rayleigh quotient method (Clough and Penzien1975), which assumes the deflected form is known (but not the magnitude of deformation). The assumption of a deflected shape reduces the tree to a single degree of freedom system. The frequency is found by equating the maximum strain energy Vmax to the maximum kinetic energy Tmax developed during the tree response. By calculating the strain and kinetic energy produced by the avalanche blast, we find the generalised stiffness K and generalised mass M of the tree:

(3) ω 2 = K M .

The natural frequency for two different deformation modes, stem bending ωsb and root-overturning ωro will be determined in the next sections.

In both cases the total tree height is H. Tree mass is divided into two parts: the stem mass ms (a mass per unit length of the tree kg m−1) and the total mass of the branches Mb (kg). The branch mass, including the mass of needles, is lumped at the tree centre-of-mass. The mass Mb can include the mass of snow held by the branches and thus, like the tree elasticity, have some seasonal variation. As we assume a constant stem diameter d the stem mass per unit length is,

(4) m s = ρ t A t

with

(5) A t = π 4 d 2 .

The density of the stem wood is ρt. For both the bending and overturning cases, the concentrated load F0 acts at the tree centre-of-mass, which is located a distance a from the ground (see Figs. 2 and 3).

Finally, the third assumption, the maximum response of the tree, will be reached before the damping forces can absorb the energy of the air blast. Only the undamped response to a short duration blast is considered.

https://www.nat-hazards-earth-syst-sci.net/18/759/2018/nhess-18-759-2018-f02

Figure 2A tree of height H breaks in bending. The avalanche exerts a loading p(t) of known (but short) duration. The load acts in the centre-of-mass of the tree located a distance a from the ground. The mass of the linear distributed mass of the tree stem is mt and the lumped mass of the branches is Mb. Tree deformation is given by the non-linear distribution x(z).

Download

2.1 Eigenfrequency: tree bending mode

For the case of tree bending, the deformation x(z) at height z is given by (see Fig. 2):

x1(z)=X0ψ1(z)=Fa2(3H-a)3EI3az2-z32a2(3H-a)(6)forza

and

x2(z)=X0ψ2(z)=Fa2(3H-a)3EI3za2-a32a2(3H-a)(7)forz>a,

where E is the modulus of elasticity of the tree stem and I is the moment of inertia. The functions ψ1(z) and ψ2(z) represent interpolation functions for the deformation field. These equations for lateral tree deformation are found by assuming the tree is a statically determinate cantilever-type structure fixed at the base to the ground (see Fig. 2 and Clough and Penzien1975). The largest bending moment in the tree is found at the tree base, z= 0. The quantity X0 is the static deformation under the blast load F,

(8) X 0 = F a 2 ( 3 H - a ) 3 E I .

The moment of inertia is taken for circular stem sections,

(9) I = π d 4 64 .

The maximum potential strain energy in bending is as follows (Clough and Penzien1975)

(10) V max = 1 2 X 0 2 0 a E I ( z ) x 1 2 ( z ) d z = 1 2 3 E I a 2 ( 3 H - a ) X 0 2 .

In the bending case, the tree is firmly rooted in the ground and strain energy is stored in the tree stem between the ground and the point of load application z=a. The tree stem above z>a is stress free, swaying back and forth as a rigid body. The maximum kinetic energy Tmax is composed of two parts containing the stem energy Tmaxs and the branch energy Tmaxb of the tree, Tmax=Tmaxs+Tmaxb (Clough and Penzien1975):

Tmaxs=ωsb220amsx12(z)dz+ωsb22aHmsx22(z)dz(11)=1280ms105H3-105aH2+35Ha2-2a3(3H-a)2X02,

and

(12) T max b = M b ω sb 2 2 x 1 2 ( z = a ) = M b ω sb 2 2 X 0 2 a 2 ( 3 H - a ) 2 .

The eigenfrequency ωsb2 is found by equating Tmax=Vmax:

(13) ω sb 2 = 420 E I ( 3 H - a ) a 2 m s 105 H 3 - 105 a H 2 + 35 H a 2 - 2 a 3 + 140 a 2 M b m t .
https://www.nat-hazards-earth-syst-sci.net/18/759/2018/nhess-18-759-2018-f03

Figure 3A tree of height H breaks by overturning at the root-plate. The avalanche exerts a loading p(t) of known (but short) duration. The load acts in the centre-of-mass of the tree located a distance a from the ground. The mass of the linear distributed mass of the tree stem is mt and the lumped mass of the branches is Mb. Tree deformation is given by the linear distribution v(z).

Download

2.2 Eigenfrequency: tree overturning mode

For the tree overturning case,

(14) x ( z ) = X 0 ψ ( z ) = F a H k z H ,

where k is the overturning stiffness of the root-plate. This equation is found by assuming the lateral tree deformation is governed by a torsional spring, representing the stiffness of the root-plate (see Fig. 3 and Chajes1974). The maximum potential strain energy (overturning) is then

(15) V max = 1 2 F X 0 = 1 2 k a H X 0 2 .

Similar to the bending case, the maximum kinetic energy is found by considering the stem and branch energies separately:

(16) T max s = ω ro 2 2 a H m s x 1 2 ( z ) d z = 1 6 m s a 3 H 2 X 0 2

and

(17) T max b = M b ω ro 2 2 x 2 ( z = a ) = M b ω ro 2 2 X 0 2 a 2 H 2 .

The eigenfrequency ωro2 is found by equating Tmax=Vmax:

(18) ω ro 2 = 3 m s a + 3 M b H k a 3 .
3 Dynamic magnification of avalanche blast

The equation of motion for an undamped system subjected to a harmonic loading is as follows:

(19) M x ¨ ( t ) + K x ( t ) = F ( t ) = F 0 sin ω t

which has the general solution for 0 tt0,

(20) x ( t ) = F 0 K 1 1 - β 2 ( sin t - β sin ω t )

and for t>t0:

(21) x ( t ) = x ˙ t 0 ω sin ω t - t 0 - x t 0 sin ω t - t 0 ,

where β=ωω is the ratio between the frequency of the avalanche blast and eigenfrequency of the tree. The magnitude of the dynamic response therefore depends on the ratio of the load duration to the period of vibration of the tree. For the case when β< 1 the maximum deformation occurs when the impulsive load is active. It can be shown (see Clough and Penzien1975) that the time to this peak response tmax is:

(22) ω t max = 2 π β β + 1 ,

which can be substituted into the general solution to find the dynamic magnification factor for a long duration impulse:

(23) D = 1 1 - β 2 sin ω t max - β sin ω t max β .

It can likewise be shown that the maximum response for the free vibration case occurs when β> 1, t>t0. For this case, the dynamic magnification factor for a short duration impulse is:

(24) D = 2 β 1 - β 2 cos π 2 β .

For the resonance case β= 1

(25) D = π 2 .

Table 1Numerical values for the mass distribution of spruce for different tree heights. Table is constructed from data contained in Indermühle (1978), Kalberer (2006) and Kramer (1988). The stated values represent average values for spruce trees in alpine environments. Values are approximate and will change depending on their location in forest, slope exposition, etc. Branch mass includes needle mass, which is given in parenthesis. Intercepted snow mass is not included in the calculations.

Download Print Version | Download XLSX

4 Application

To demonstrate how the dynamic magnification factor D can be found, we consider the following problem: a powder snow avalanche enters a spruce forest with considerable speed (> 50 m s−1) and exerts a short duration air-blast with frequency ω. The duration of the blast is on the order of a few seconds. The height of the trees is between 25 and 30 m, which is also the height of the powder cloud. The cloud has decoupled from the avalanche core which has stopped before reaching the forest. Moreover, the only loading on the trees is the air-blast.

Table 2Natural frequencies in bending and overturning for spruce trees of different heights. E= 10 GPa. A reduced stem diameter d= 0.5dDBH produces a good agreement to measured frequencies. Mass distribution taken from Table 1.

Download Print Version | Download XLSX

Using the measured mass values tabulated in Table 1, we set the total branch and needle mass of a single tree to be Mb= 540 kg. The stem mass per length is approximately 60 kg m−1 (wood density 480 kg m−3). The total force of the avalanche impact acts at the tree's centre-of-mass which is located a= 16.5 m above ground. This allows us to define the natural frequency in bending of the tree by Eq. (13), ωsb= 1.48 rad s−1 (0.24 Hz), see Table 2. This value is in very good agreement with the measurements (see Jonsson et al.2007). The modulus of elasticity was set to E= 10 GPa based on experimental measurements (Haines et al.1996). For the calculations, a tree diameter somewhat smaller than the diameter at breast height (DBH) is selected. In this case d= 0.2 m, which is 1∕2 of the DBH diameter (this provides the best match to the experimental frequencies).

Consider first a duration sine impulse lasting 2.50 s (ω=π6). In this case β= 0.699; that is, the maximum deformation occurs during the time the load is acting. For this case, application of Eq. (23), we find D= 1.76, a rather large magnification factor. For a shorter duration impulse lasting 1.66 s, β= 1.27 and from Eq. (24), we find D= 1.36. The primary conclusion to draw from this analysis is that the natural frequency in bending of tall trees is close to the frequency of the applied avalanche air-blast. Measurements of air-blast duration times reported by Russian researchers are within this range, lasting only a few seconds (see Grigoryan et al.1982; Sukhanov and Kholobaeva1982; Sukhanov1982).

Measurements of root plate stiffness are rare; however, values for 10–14 m high spruce reported in Neild and Wood (1998) vary between k= 80 kN m (H= 10 m) and k= 1200 kN m (H= 14 m). These values suggest a large variation in k depending on growth conditions. The application of these k stiffness values for spruce trees predicts natural frequencies for root-plate overturning in ωo> 2 Hz (Eq. 18), see Table 2. The calculated β factors for overturning are typically β< 1. This result suggests that large dynamic magnification factors can only be generated by very short duration impulses (less than t< 0.5 s). Tall trees (H> 20 m) with low root plate stiffness (k 100 kN m) are vulnerable to powder avalanche air-blasts.

5 Conclusions

We draw several conclusions from our analysis. Firstly, the natural frequency of tall trees – in bending and overturning – is close to the loading frequency of powder avalanches, ωω. Thus, tall trees are susceptible to powder avalanche blow-down. When using tree blow-down to estimate avalanche impact pressures (and therefore speed and density of the powder cloud) a dynamic magnification factor should be applied in the analysis. Moreover, powder avalanches can knock down trees with lower velocity than is presently assumed. This result is also valid for other types of tall structures, including power pylons, or buildings with long over-hanging roofs.

Secondly, both tree bending and root-plate overturning are possible tree failure modes when hit by a powder avalanche. Interestingly, the natural frequencies of tree bending and root-plate overturning are similar, when the root-plate stiffness is low (k< 100 kN m) and the tree is tall (H> 20 m). Although there is considerable data available to constrain the value of the modulus of elasticity of wood E, there is less information available to constrain the root-plate stiffness. In the future, field investigations that document forest destruction should clearly separate bending and overturning failures. This would help understand the variability of tree anchorage on mountain slopes. The field examinations should also quantify the stem diameter d at more than one location as this is necessary to accurately determine the bending eigenfrequency.

Finally, the fact that tall trees can be broken in bending and overturning indicates the nature of the avalanche air blast. It appears to be a high velocity, short duration pulse of flowing material (ice-dust), similar to a high-density gust of wind. It is not a compression wave travelling at the speed of sound.

Data availability

No data sets were used in this article.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

This work was performed within the framework of the joint Austrian-Swiss project bDFA, a study of avalanche motion beyond the dense flow avalanche regime. We thank the Austrian Academy of Science (ÖAW) for their financial support as well as the Austrian research partners (Austrian Research Centre for Forests, Torrent and Avalanche Control and the University of Innsbruck).

Edited by: Oded Katz
Reviewed by: two anonymous referees

References

Bartelt, P. and Stöckli, V.: The influence of tree and branch fracture, overturning and debris entrainment on snow avalanche flow, Ann. Glaciol., 32, 209–216, 2001. a

Bozhinskiy, A. N. and Losev, K. S.: The fundamentals of avalanche science, Mitt. Eidgenöss. Inst. Schnee- Lawinenforsch., Davos, p. 280, 1998. a

Chajes, A.: Principles of Structural Stability Theory, Prentice Hall Inc, Englewood Cliffs, p. 336, 1974. a

Clough, R. W. and Penzien, J.: Dynamics of Structures, McGraw-Hill Inc, New York, p. 634, 1975. a, b, c, d, e, f, g

Coutts, M.: Root architecture and tree stability, Plant Soil, 71, 171–188, 1983. a

Feistl, T., Bebi, P., Christen, M., Margreth, S., Diefenbach, L., and Bartelt, P.: Forest damage and snow avalanche flow regime, Nat. Hazards Earth Syst. Sci., 15, 1275-1288, https://doi.org/10.5194/nhess-15-1275-2015, 2015a. a, b

Feistl, T., Bebi, P., Teich, M., Bühler, Y., Christen, M., Thuro, K., and Bartelt, P.: Observations and modeling of the braking effect of forests on small and medium avalanches, J. Glaciol., 60, 124–138, https://doi.org/10.3189/2014JoG13J055, 2015b. a

Gadek, B., Kaczka, R. J., Raczkowska, Z., Rojan, E., Casteller, A., and Bebi, P.: Snow avalanche activity in Zleb Zandarmerii in a time of climate change (Tatra Mts., Poland), Catena, 158, 201–212, 2017. a

Grigoryan, S., Urubayev, N., and Nekrasov, I.: Experimental investigation of an avalanche air blast, Data Glaciol. Stud., 44, 87–93, 1982. a, b

Haines, D. W., Leban, J. M., and Herbe, C.: Determination of Young's modulus for spruce, fir and isotropic materials by the resonance flexure method with comparisons to static flexure and other dynamic methods, Wood Sci. Technol., 30, 253–263, 1996. a

Indermühle, M. P.: Struktur, Alters- und Zuwachsuntersuchungenin einem Fichten-Plenterwald der subalpinen Stufe, in: Beiheft Nr. 60 zur Schweiz, Z. Forstwesen, Dissertaion ETHZ-Zürich, Zürich, p. 98, 1978.  a

Johnson, E. A.: The relative importance of snow avalanche disturbance and thinning on canopy plant populations, Ecology, 68, 43–53, 1987. a

Jonsson, M. J., Foetzki, A., Kalberer, M., Lundström, T., Ammann, W., and Stöckli, V.: Root-soil rotation stiffness of norway spruce (Picea abies (L.) Karst) growing on subalpine forested slopes, Plant Soil, 285, 267–277, 2006. a

Jonsson, M. J., Foetzki, A., Kalberer, M., Lundström, T., Ammann, W., and Stöckli, V.: Natural frequencies and damping ratios of Norway spruce (Picea abies (L.) Karst) growing on subalpine forested slopes, Trees, 21, 541–548, https://doi.org/10.1007/s00468-007-0147-x, 2007. a

Kalberer, M.: Quantifizierung und Optimierung der Schutzwaldleistung gegenüber Steinschlag, Dissertaion Albert-Ludwigs-Universität, Freiburg, 2006. a

Kramer, H.: Waldwachstumslehre, Parey, Hamburg, Berlin, p. 374, 1988. a

Mattheck, C. and Breloer, H.: Handbuch der Schadenskunde von Bäumen: Der Baumbruch in Mechanik und Rechtsprechung, Rombach, Freiburg im Breisgau, 1994. a, b

Neild, S. A. and Wood, C. J.: Estimating stem and root-anchorage flexibility in trees, Tree Physiol., 19, 141–151, 1998. a

Peltola, H., Nykänen, M. L., and Kellomäki, S.: Model computations on the critical combination of snow loading and windspeed for snow damage of scots pine, norway spruce and birch at stand edge, Forest Ecol. Manage., 95, 229–241, 1997. a

Peltola, H., Kellomäki, S., Väisänen, H., and Ikonen, V.: A mechanisticmodel for assessing the risk of wind and snow damage to single trees and stands of scots pine, norway spruce, and birch, Can. J. Forest Res., 29, 647–661, 1999. a

Reardon, B. A., Pederson, G. T., Caruso, C. J., and Fagre, D. B.: Spatial Reconstructions and Comparisons of Historic Snow Avalanche Frequency and Extent Using Tree Rings in Glacier National Park, Montana, U.S.A., Arct. Antarct. Alp. Res., 40, 148–160, 2008. a

Schläppy, R., Eckert, N., Jomelli, C., Stoffel, M., Grancher, D., Brunstein, D., Naaim, M., and Deschatres, M.: Validation of extreme snow avalanches and related return periods derived from a statistical-dynamical model using tree-ring techniques, Cold Reg. Sci. Technol., 99, 12–26, 2004. a

Sukhanov, G.: The mechanism of avalanche air blast formation as derived from field measurements, Data Glaciol. Stud., 44, 94–98, 1982. a, b

Sukhanov, G. and Kholobaeva, P.: Vriability of avalanche air blast in time and space, Data Glaciol. Stud., 44, 98–105, 1982. a, b

Download
Short summary
We study how short duration powder avalanche blasts break and overturn tall trees. Tree blow-down is often used to back-calculate avalanche pressure and therefore constrain avalanche flow velocity and motion. We find that tall trees are susceptible to avalanche air blasts because the duration of the air blast is near to the period of vibration of tall trees. Dynamic magnification factors should therefore be considered when back-calculating powder avalanche impact pressures.
Altmetrics
Final-revised paper
Preprint