The Translatory Wave Model for Landslides

*Jónas Elíasson and Þorsteinn Sæmundsson*

## **Abstract**

The Saint-Venant equations are usually the basis of numerical models for landslide flows. They are nonstationary and nonlinear. The theory for translatory waves in a prismatic channel and a funneling channel can be used for landslides using the assumption of either turbulent or laminar flow in the slide. The mathematics of translatory waves traveling over dry land or superimposed on another flow are developed. This results in a new slope factor controlling the flow velocity, together with the Chezy coefficient used in previous applications of the translatory wave theory. Flow times for the slide to reach a given destination, slide depth, and velocity can be calculated using the initial magnitude of the flow in the slide. The instabilities of the wave tail are discussed. Three case studies are presented: a submarine slide that started the Tohoku tsunami in Japan, the Morsárjökull rock avalanche in SE Iceland, and the Móafellshyrna slide in central N Iceland.

**Keywords:** landslide model, translatory wave, landslides, rock avalanche, debris slides, Iceland

## **1. Introduction**

Landslides are a natural phenomenon that more often than not can be described as an environmental disaster. Landslides are responsible for countless human lives, and there is no remedy against their occurrence. The only way to counteract the danger is preparedness, knowing the danger, assessing the hazard, operating a warning system, and keeping people out of the way of the slides, a task all governments and rescue services, civil and military, deal with every day.

Landslides are very different in character, so in the civil protection work, everything matters, area, weather, soil, topography, surface structure, and what triggers the slide. Consequently, the classification system for landslides is very complicated. A very good overview of the modern trend of classification can be found in [1] where types of landslides are discussed separately. A quick overview on what this classification is all about can be found in the course material in [2], where classification systems from 1938, 1958, and 1978 are listed plus a unified classification system.

The classification has to take into account all the danger factors previously mentioned and is therefore quite complicated and in constant development, and a lot has happened since it appeared. For instance, today statistical methods are used to assess the danger [3], and sometimes they bring about scaling laws known from

the fractal theory [4]. The statistical methods have not found their way into the classification system.

**2.2 Translatory waves in a flume**

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

*q* channel flow per unit width

*x* streamwise coordinate

*g* acceleration of gravity

The equations of Saint-Venant describing flow down an inclining plane can be

The system obviously admits a solution where *v = v(x* � *ct)* and *y = y(x* � *ct)* with *c* as the constant wave celerity. The momentum equation degenerates to

*gy* <sup>¼</sup> *<sup>I</sup>*<sup>0</sup> *<sup>c</sup>*

*x* � *ct y*0

This is the form of the approaching wave. The wave passes the origin at time

As time progresses y approaches the normal depth *y0*, while the discharge *q*

*<sup>q</sup>*<sup>0</sup> <sup>¼</sup> *cy*<sup>0</sup> <sup>¼</sup> *<sup>C</sup>* ffiffiffiffiffiffi

Call the shear stress *τ* and the pressure on the bottom (normal stress) *σ* in uniform flow of a soil in a landslide down a plane inclined by an angle *Ѳ* to the

from which the normal depth *y0* and *c* can be calculated if *q0* is known. Equation

*gI*0 p *y* 3*=*2

It is easily seen that the wave will continue to run if *I > 0*, i.e., the flow will not stop until the surface has become horizontal. This means that the landslide is of liquefied, or fluidized, soil. This means that the water content of the slide is rather high. This does not have to be the case. The slide's material may very well have

where C is the nondimensional Chezy constant and ξ = x � ct. This equation can

<sup>2</sup> <sup>¼</sup> *<sup>C</sup>*<sup>2</sup>

*δx* þ *∂q*

*<sup>∂</sup><sup>x</sup>* <sup>¼</sup> <sup>0</sup> *<sup>q</sup>* <sup>¼</sup> *yv* (1)

*<sup>∂</sup><sup>x</sup> <sup>g</sup>*ð*<sup>I</sup>* � *<sup>I</sup>*0Þ ¼ <sup>0</sup> (2)

*gIy* (3)

; *y*<sup>0</sup> ¼ 0; *x* � *ct*≤0 (4)

<sup>0</sup> (5)

taken from many textbooks on fluid mechanics or hydrology; they are

*δt* þ *v ∂v ∂x* þ *g ∂v*

Equation of continuity: *<sup>δ</sup><sup>y</sup>*

Momentum eq*:*: *<sup>δ</sup><sup>v</sup>*

*dy dξ* þ *c*2 *C*2

*y*0 � �

¼ *I*<sup>0</sup>

*2.2.1 The translatory wave*

where: *y* water depth

*v* velocity

*I* friction slope *I0* channel slope

readily be integrated.

zero: *x=t=0=> y = 0.*

approaches the value

**143**

*y y*0 <sup>þ</sup> ln 1 � *<sup>y</sup>*

finite shear strength. Then it stops at a finite slope.

*2.2.2 Basic fluid mechanics of landslide flows*

(4) can also be used for landslides, as will be demonstrated.

*t* time

One of the largest uncertainties and in many instances the very reason for the use of statistics in the science of landslides is the instability problem. Instability can trigger a landslide quite suddenly, and then again, small cracks may be discovered, and they keep widening for a long time; it is quite evident that a landslide will occur at some time, but there is no way to tell when it happens.

Instability problems are well known in soil mechanics, e.g., liquefaction and wave formation in hydrodynamics. Both can be triggered by an external force like earthquake. A lot of research effort is devoted to the problem of landslide triggering [5], sometimes with stochastic input routines [6, 7] and others with models describing the aftermath, for hazard assessment [8], or modeling of the flow of the slide itself [9], sometimes with assistance of remote sensing data [10].

Landslide models can be numerical or analytical. Both are very dependent of the boundary conditions. Both need a physical basis to describe the dynamics of the slide, especially the particle velocity. The simplest, but very efficient, method is to put the velocity proportional to the gradient of the elevation of the surface on the one hand and to the shear force in excess of the shear strength of the soil on the another. An example of the mathematical formulation this assumption leads to is Eqs. 1 and 2 in [5]. The basis of this approach is Newton's law of viscosity.

A necessary condition for this approach to hold is that the inertial terms, other than acceleration, in the full flow equation system can be neglected. This means that all nonlinear terms can be left out. This will be the case for many types of landslides, due to the high viscosity in them. However, many slides will be running faster than that, and they are the most dangerous ones.

In the following we shall discuss the translatory wave landslide model. This model can be used for the front of the landslide, the velocity of it, and how far it flows. The criterion for using the model is that the slide has to be flowing in a regular manner, not just a fall of soil or rock fragment down a steep slope, so many of the types defined in the classification systems and discussions of different types mentioned above will be excluded from this approach. For the others, i.e., regular slides that do actually flow, the translatory theory can in most cases be used.

The translatory model involves analyses that are complicated in nature but lead to a simple system when fully derived. It can be very useful in studies one has to do prior to numerical model studies in defining the area, extent of boundaries, and limits to the various parameters involved. Thus, it can save a lot of trial and error work in constructing the model and elimination of numerical instabilities. And sometimes it produces results that are good enough as demonstrated in the case studies.

## **2. Theoretical considerations**

#### **2.1 Theory of translatory waves**

A translatory wave runs down a channel slope with flow that is not constant in time but varies slowly. In treatment of such flow, it is natural to look for waves of fixed shape progressing down a channel. Such waves in shallow water have been investigated earlier; a further treatment is presented in [11] using Chezy's equation for the flow resistance instead of the more complicated Manning's equation. This simplifies the mathematics of the problem. The starting point is the continuity and momentum equations for an arbitrary prismatic channel. For the sake of simplicity, only two cases are treated, that of a wide rectangular inclining flume and a rectangular inclining funnel.

## **2.2 Translatory waves in a flume**

## *2.2.1 The translatory wave*

The equations of Saint-Venant describing flow down an inclining plane can be taken from many textbooks on fluid mechanics or hydrology; they are

$$\text{Equation of continuity: } \frac{\delta y}{\delta \mathbf{x}} + \frac{\partial q}{\partial \mathbf{x}} = \mathbf{0} \neq \mathbf{y}v \tag{1}$$

$$\text{Momentum eq.: } \frac{\delta v}{\delta t} + \nu \frac{\partial v}{\partial \mathbf{x}} + \mathbf{g} \frac{\partial v}{\partial \mathbf{x}} \text{ g} (I - I\_0) = \mathbf{0} \tag{2}$$

where:

the fractal theory [4]. The statistical methods have not found their way into the

at some time, but there is no way to tell when it happens.

that, and they are the most dangerous ones.

**2. Theoretical considerations**

**2.1 Theory of translatory waves**

gular inclining funnel.

**142**

One of the largest uncertainties and in many instances the very reason for the use of statistics in the science of landslides is the instability problem. Instability can trigger a landslide quite suddenly, and then again, small cracks may be discovered, and they keep widening for a long time; it is quite evident that a landslide will occur

Instability problems are well known in soil mechanics, e.g., liquefaction and wave formation in hydrodynamics. Both can be triggered by an external force like earthquake. A lot of research effort is devoted to the problem of landslide triggering

describing the aftermath, for hazard assessment [8], or modeling of the flow of the

Landslide models can be numerical or analytical. Both are very dependent of the boundary conditions. Both need a physical basis to describe the dynamics of the slide, especially the particle velocity. The simplest, but very efficient, method is to put the velocity proportional to the gradient of the elevation of the surface on the one hand and to the shear force in excess of the shear strength of the soil on the another. An example of the mathematical formulation this assumption leads to is Eqs. 1 and 2 in [5]. The basis of this approach is Newton's law of viscosity.

A necessary condition for this approach to hold is that the inertial terms, other than acceleration, in the full flow equation system can be neglected. This means that all nonlinear terms can be left out. This will be the case for many types of landslides, due to the high viscosity in them. However, many slides will be running faster than

In the following we shall discuss the translatory wave landslide model. This model can be used for the front of the landslide, the velocity of it, and how far it flows. The criterion for using the model is that the slide has to be flowing in a regular manner, not just a fall of soil or rock fragment down a steep slope, so many of the types defined in the classification systems and discussions of different types mentioned above will be excluded from this approach. For the others, i.e., regular slides that do actually flow, the translatory theory can in most cases be used.

The translatory model involves analyses that are complicated in nature but lead to a simple system when fully derived. It can be very useful in studies one has to do prior to numerical model studies in defining the area, extent of boundaries, and limits to the various parameters involved. Thus, it can save a lot of trial and error work in constructing the model and elimination of numerical instabilities. And sometimes it

A translatory wave runs down a channel slope with flow that is not constant in time but varies slowly. In treatment of such flow, it is natural to look for waves of fixed shape progressing down a channel. Such waves in shallow water have been investigated earlier; a further treatment is presented in [11] using Chezy's equation for the flow resistance instead of the more complicated Manning's equation. This simplifies the mathematics of the problem. The starting point is the continuity and momentum equations for an arbitrary prismatic channel. For the sake of simplicity, only two cases are treated, that of a wide rectangular inclining flume and a rectan-

produces results that are good enough as demonstrated in the case studies.

[5], sometimes with stochastic input routines [6, 7] and others with models

slide itself [9], sometimes with assistance of remote sensing data [10].

classification system.

*Landslides - Investigation and Monitoring*

*y* water depth *q* channel flow per unit width *v* velocity *x* streamwise coordinate *t* time *g* acceleration of gravity *I* friction slope *I0* channel slope

The system obviously admits a solution where *v = v(x* � *ct)* and *y = y(x* � *ct)* with *c* as the constant wave celerity. The momentum equation degenerates to

$$\frac{d\mathbf{y}}{d\xi} + \frac{c^2}{\mathbf{C}^2 \mathbf{y}} = I\_0 \quad c^2 = \mathbf{C}^2 \mathbf{g} \mathbf{J} \tag{3}$$

where C is the nondimensional Chezy constant and ξ = x � ct. This equation can readily be integrated.

$$\frac{\mathcal{Y}}{\mathcal{Y}\_0} + \ln\left(1 - \frac{\mathcal{Y}}{\mathcal{Y}\_0}\right) = I\_0 \frac{\mathcal{X} - ct}{\mathcal{Y}\_0}; \quad \mathcal{Y}\_0 = \mathbf{0}; \quad \mathcal{X} - ct \le \mathbf{0} \tag{4}$$

This is the form of the approaching wave. The wave passes the origin at time zero: *x=t=0=> y = 0.*

It is easily seen that the wave will continue to run if *I > 0*, i.e., the flow will not stop until the surface has become horizontal. This means that the landslide is of liquefied, or fluidized, soil. This means that the water content of the slide is rather high. This does not have to be the case. The slide's material may very well have finite shear strength. Then it stops at a finite slope.

As time progresses y approaches the normal depth *y0*, while the discharge *q* approaches the value

$$q\_0 = c \mathbf{y}\_0 = \mathbf{C} \sqrt{\mathbf{g} I\_0} \mathbf{y}\_0^{3/2} \tag{5}$$

from which the normal depth *y0* and *c* can be calculated if *q0* is known. Equation (4) can also be used for landslides, as will be demonstrated.

## *2.2.2 Basic fluid mechanics of landslide flows*

Call the shear stress *τ* and the pressure on the bottom (normal stress) *σ* in uniform flow of a soil in a landslide down a plane inclined by an angle *Ѳ* to the horizontal, and let *y*´ be the distance from the bottom and *y* the total depth. Then we have from any textbook in hydraulics

$$
\pi = \rho \text{g} \sin \theta (y - y'); \quad \sigma = \rho \text{g} \cos \theta (y - y'); \quad \rightarrow \quad \pi/\sigma = \tan \theta \tag{6}
$$

laminar friction. Our implicit relation would have to possess that property too, but otherwise it remains unknown until sufficient experimental data for landslides is

We call *Sf = sinѲ (1 - tanφ´/tanѲ)* the slope factor. It can be found in other models; see, e.g., (5) in [8]. We take the value 0.5. and list in **Table 1** corresponding

*C* is proportional with the velocity but varies much more slowly with y. It is therefore a possibility that approximate average regional values for *C* can exist in regions with similar geological settings. *C* for other slope factors *Sf* can easily be

When *C* is known, we can turn back to (4) and modify it to fit a slide where the flow resistance is the excess shear stress. In Section 2.3.1, we defined the slope factor for uniform flow from the Chezy coefficient *C* and the bottom slope. Now *C* is defined in the same way but for local surface slope *Sf*, not bottom slope *Sf0.* With

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

*x* � *ct y*0

The mathematical properties of Eq. (11) are discussed in [11] and the function

It is obvious from (11) that *y* approaches *yo* very fast with time. According to (4) the flood reaches a depth that is 90% of *yo* in *x=0* in time *t =* 1.4 *yo/(Sfo c).* This is only a few minutes according to the values in **Table 1**; the slower the slide, the

/s per m width in the slide is found by modifying (5) to (12):

*y* 3*=*2

ffiffiffiffiffiffiffiffi *gSf*<sup>0</sup> q

1 � tan *ϕ*<sup>0</sup> ð Þ *=* tan *θ*

¼ *C*

ffiffiffiffiffiffiffiffi *gySf* q

; *y*<sup>0</sup> ¼ 0; *x* � *ct*≤0 (12)

<sup>0</sup> (13)

(11)

available.

**Example 1.**

*C* values for various slides, calculated from (10).

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

*2.2.3 Properties of the slide function*

*ξ = x – ct* as before, the local velocity will be

s

<sup>þ</sup> ln 1 � *<sup>y</sup>*

steeper is the front of the approaching slide.

*c* ¼ *C*

And (4) is modified to

plotted in **Figure 1** in [12].

The flow *q0* in m<sup>3</sup>

**Table 1.**

**145**

*C for y in m, c in m/s, and slope factor 0.5.*

*y y*0 *gy* sin *θ*

*y*0 � � *dy dξ*

¼ *Sf*<sup>0</sup>

*q*<sup>0</sup> ¼ *cy*<sup>0</sup> ¼ *C*

*y/c* **0.5 1 2 3** 0.22 0.45 0.89 1.34 0.16 0.32 0.63 0.95 0.13 0.26 0.52 0.77 0.11 0.22 0.45 0.67 0.10 0.20 0.40 0.60 0.09 0.18 0.37 0.55 0.08 0.17 0.34 0.51 0.08 0.16 0.32 0.47

found in **Table 1** by multiplying the *C* values with *a* = √(*0.5/Sf*).

If the shear resistance of the soil is purely frictional with friction angle *φ*, there is no flow for *φ > Ѳ.* For cohesive soils this is a little more complicated. Motion starts when *τ > σ tan φ + c´*, where *c´* is the cohesion. In the no flow situation, the cohesion can be a substantial part of the total shear strength of the soil; however, the cohesion is not a material constant for the soil; upon deformation and start of flow, *c´* is reduced, and after that we can approximate the shear strength of the soil with a dynamical friction angle *φ´* and a zero cohesion. Then we will have *φ´ < Ѳ* as condition for flow. This means that the slide will stop when it flows on a slope *Ѳ < φ´,* and then the shear stress ratio is *τcr/σ = tanφ ´,* which is the critical shear stress ratio. The shear stress in excess of the critical *τcr* is then

$$(\tau - \tau\_{cr})/\sigma = \tan \theta - \tan \phi'; \quad \rightarrow \quad \tau - \tau\_{cr} = \rho \text{g} \sin \theta (1 - \tan \phi'/\tan \theta)(y - y') \tag{7}$$

This can be combined with Newton's law of viscosity.

$$(\pi - \pi\_{cr})/\rho = \text{g}\,\sin\theta (y - y') (\mathbf{1} - \tan\phi' / \tan\theta) = \mathbf{C}\_f \frac{dv}{dy'} \tag{8}$$

$$\text{Here we have:} \begin{cases} \text{Laminar flow:} & \text{C}\_{\text{f}} = \text{kinematic viscosity} \\ \text{Turbulent flow:} & \text{C}\_{\text{f}} = \text{eddy viscosity} \end{cases} \tag{9}$$

Kinematic viscosity is normally assumed constant, but eddy viscosity changes with the distance from the wall in accordance with Prandtl's theory of wall turbulence. This makes the difference that the velocity profile will be parabolic in the laminar case and the maximum velocity in a fixed distance in the flow path, i.e., *vmax = 1.5 vaverage* for a fixed *x* in (4). In the turbulent case, *vmax* is closer to *vaverage* than that.

Combining this with (4) and (5), we put *c=vav*, and the flow becomes *q=cy*. The actual variation of the velocity with *y´* will for all practical purposes affect the wave form (4) only slightly when the celerity in the landslide is *c* < < √(*gy*) which is the celerity of the surface gravity wave.

Equations (8) and (9) are similar to those for ordinary channel flow except (6) is multiplied by *(1 - tan φ´/tanѲ)* to get (8), which gives the excess shear stress or the resistance to the flow. Then (9) provides the relation to the velocity through the unknown resistance coefficient *Cf*. In the laminar case, it is assumed constant, and then we have a flow proportional to the slope, which in a free landscape gives a flow in the direction of the elevation gradient. This is a popular idea for landslide model (see, e.g., [6, 7]).

*Cf* does not need to be constant, it can depend on the velocity, giving a turbulence-like flow resistance. A simple way to include this possibility is to use Chezy's coefficient defined as

$$\mathcal{L} = \nu\_{av} = \mathcal{C}\sqrt{(\tau - \tau\_{cr})/\rho} = \mathcal{C}\sqrt{(\text{gy}\sin\theta(\mathbf{1} - \tan\phi'/\tan\theta))}\tag{10}$$

*C* is not a constant as originally proposed by Chezy but a slowly varying implicit function of the velocity and is thus constant while *c* is constant. To get an idea about the variation of *C* with the resistance factor in turbulent flows, one can study Section 15.2 in [12]. There is an implicit relation proposed by ASCE in 1963; it is based on a large volume of experimental data. For low velocity flows, it gives

laminar friction. Our implicit relation would have to possess that property too, but otherwise it remains unknown until sufficient experimental data for landslides is available.

## **Example 1.**

horizontal, and let *y*´ be the distance from the bottom and *y* the total depth. Then we

*τ* ¼ *ρg* sin *θ*ð Þ *y* � *y*´ ; *σ* ¼ *ρg* cos *θ*ð Þ *y* � *y*´ ; ! *τ=σ* ¼ tan *θ* (6)

; ! *τ* � *τcr* ¼ *ρg* sin *θ* 1 � tan *ϕ*<sup>0</sup> ð Þ *=* tan *θ* ð Þ *y* � *y*´

*dv*

*dy*<sup>0</sup> (8)

(7)

(9)

(10)

If the shear resistance of the soil is purely frictional with friction angle *φ*, there is no flow for *φ > Ѳ.* For cohesive soils this is a little more complicated. Motion starts when *τ > σ tan φ + c´*, where *c´* is the cohesion. In the no flow situation, the cohesion can be a substantial part of the total shear strength of the soil; however, the cohesion is not a material constant for the soil; upon deformation and start of flow, *c´* is reduced, and after that we can approximate the shear strength of the soil with a dynamical friction angle *φ´* and a zero cohesion. Then we will have *φ´ < Ѳ* as condition for flow. This means that the slide will stop when it flows on a slope *Ѳ < φ´,* and then the shear stress ratio is *τcr/σ = tanφ ´,* which is the critical shear

stress ratio. The shear stress in excess of the critical *τcr* is then

This can be combined with Newton's law of viscosity.

�

is the celerity of the surface gravity wave.

*<sup>c</sup>* <sup>¼</sup> *vav* <sup>¼</sup> *<sup>C</sup>* ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

ð Þ *<sup>τ</sup>* � *<sup>τ</sup>cr <sup>=</sup><sup>ρ</sup>* <sup>p</sup> <sup>¼</sup> *<sup>C</sup>*

(see, e.g., [6, 7]).

**144**

Chezy's coefficient defined as

ð Þ *τ* � *τcr =ρ* ¼ *g* sin *θ y* � *y*<sup>0</sup> ð Þ 1 � tan *ϕ*<sup>0</sup> ð *=* tan *θ*Þ ¼ *Cf*

Here we have: Laminar flow: Cf <sup>¼</sup> kinematic viscosity

Kinematic viscosity is normally assumed constant, but eddy viscosity changes with the distance from the wall in accordance with Prandtl's theory of wall turbulence. This makes the difference that the velocity profile will be parabolic in the laminar case and the maximum velocity in a fixed distance in the flow path, i.e., *vmax = 1.5 vaverage* for a fixed *x* in (4). In the turbulent case, *vmax* is closer to *vaverage* than that. Combining this with (4) and (5), we put *c=vav*, and the flow becomes *q=cy*. The actual variation of the velocity with *y´* will for all practical purposes affect the wave form (4) only slightly when the celerity in the landslide is *c* < < √(*gy*) which

Equations (8) and (9) are similar to those for ordinary channel flow except (6) is multiplied by *(1 - tan φ´/tanѲ)* to get (8), which gives the excess shear stress or the resistance to the flow. Then (9) provides the relation to the velocity through the unknown resistance coefficient *Cf*. In the laminar case, it is assumed constant, and then we have a flow proportional to the slope, which in a free landscape gives a flow in the direction of the elevation gradient. This is a popular idea for landslide model

> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *gy* sin *θ* 1 � tan *ϕ*<sup>0</sup> ð Þ ð Þ *=* tan *θ*

*Cf* does not need to be constant, it can depend on the velocity, giving a turbulence-like flow resistance. A simple way to include this possibility is to use

q

*C* is not a constant as originally proposed by Chezy but a slowly varying implicit function of the velocity and is thus constant while *c* is constant. To get an idea about the variation of *C* with the resistance factor in turbulent flows, one can study Section 15.2 in [12]. There is an implicit relation proposed by ASCE in 1963; it is based on a large volume of experimental data. For low velocity flows, it gives

Turbulent flow: Cf ¼ eddy viscosity

have from any textbook in hydraulics

*Landslides - Investigation and Monitoring*

ð Þ *τ* � *τcr =σ* ¼ tan *θ* � tan *ϕ*<sup>0</sup>

We call *Sf = sinѲ (1 - tanφ´/tanѲ)* the slope factor. It can be found in other models; see, e.g., (5) in [8]. We take the value 0.5. and list in **Table 1** corresponding *C* values for various slides, calculated from (10).

*C* is proportional with the velocity but varies much more slowly with y. It is therefore a possibility that approximate average regional values for *C* can exist in regions with similar geological settings. *C* for other slope factors *Sf* can easily be found in **Table 1** by multiplying the *C* values with *a* = √(*0.5/Sf*).

## *2.2.3 Properties of the slide function*

When *C* is known, we can turn back to (4) and modify it to fit a slide where the flow resistance is the excess shear stress. In Section 2.3.1, we defined the slope factor for uniform flow from the Chezy coefficient *C* and the bottom slope. Now *C* is defined in the same way but for local surface slope *Sf*, not bottom slope *Sf0.* With *ξ = x – ct* as before, the local velocity will be

$$\mathcal{L} = \mathcal{C}\sqrt{\text{gy}\sin\theta \frac{d\mathcal{Y}}{d\xi}(\mathbf{1} - \tan\phi'/\tan\theta)} = \mathcal{C}\sqrt{\text{gy}\mathcal{S}\_f} \tag{11}$$

And (4) is modified to

$$\frac{\mathcal{Y}}{\mathcal{Y}\_0} + \ln\left(1 - \frac{\mathcal{Y}}{\mathcal{Y}\_0}\right) = \mathcal{S}\_{f0}\frac{\mathcal{X} - ct}{\mathcal{Y}\_0}; \quad \mathcal{Y}\_0 = 0; \quad \mathcal{X} - ct \le 0 \tag{12}$$

The mathematical properties of Eq. (11) are discussed in [11] and the function plotted in **Figure 1** in [12].

It is obvious from (11) that *y* approaches *yo* very fast with time. According to (4) the flood reaches a depth that is 90% of *yo* in *x=0* in time *t =* 1.4 *yo/(Sfo c).* This is only a few minutes according to the values in **Table 1**; the slower the slide, the steeper is the front of the approaching slide.

The flow *q0* in m<sup>3</sup> /s per m width in the slide is found by modifying (5) to (12):

$$q\_0 = \mathfrak{c}\mathfrak{y}\_0 = \mathbb{C}\sqrt{\mathbf{g}\mathbf{S}\_{f0}}\mathfrak{y}\_0^{3/2} \tag{13}$$


#### **Table 1.**

*C for y in m, c in m/s, and slope factor 0.5.*

#### *2.2.4 Change of slope or channel form*

If the slope changes, the velocity will adjust to a new depth and velocities in a region just before and just after the slope change as long as the flow is subcritical in the fluid mechanical sense. When the slide flows downhill, the slope will in most cases be diminished. The flow will adapt to the new depth on a short distance.

Equation (5) can be used to find the new depth by using the fact that the flow is the same, upstream and downstream of the slope change. The same approach is used if the channel changes form.

#### *2.2.5 Flow in a funnel-shaped flume*

The continuity equation in a channel shaped as a funnel flume is

$$\frac{\partial y}{\partial t} + \frac{1}{r} \frac{\partial (rq)}{\partial r} = 0; \quad q = \mathcal{y}v \tag{14}$$

that can be approximated as a series of flumes with constant width and slowly changing funnels. It is therefore quite clear that a translatory wave can flow down any channel, adapt to the changing environment, and still be a translatory wave.

The most severe limitation of the translatory model is the very heart of the theory, the progressive wave. The wave starts at some point; for a while there will be either increasing or diminishing flow in this point. While we have increasing or constant flow, we are in the head of the wave. For constant flow (11) and (12) apply, but a flow increase will come as a new translatory wave riding upon the first

Diminishing flow is more complicated. Then we are in the tail of the wave head, the wave is recessing, and it is no longer progressing with a constant speed. If it is necessary to calculate this flow, numerical simulation using the full Saint-Venant equations, modified according to Section 2.2.3, may be tried. This involves more advanced programming than for a water waves and has not been tried for landslides using nonlinear friction. In [13] a very advanced model of this kind is presented, using linear friction, but allowing different rheologies, both viscous and non-Newtonian flow, which should give the same possibilities, as nonlinear friction coupled with a slope factor can be used as a tool to imitate non-Newtonian behavior. The head of the wave is the most damaging part, and the head also determines the length of the damaged area. Therefore, the head flow is the most important, and

If it is clear from the beginning that the velocity is low enough, so all inertial and nonlinear terms can be neglected, and the gradient flow model can be applied in a way like [6, 7]. Using a notation similar to theirs, but using *Δy* instead of *ΔH* (with *Δ* meaning the gradient/divergence operator), we would have a slope factor as in

Ѳ is the bottom slope in the direction of the flow. Using the laminar model from

*Sf y*<sup>3</sup>

Δ*y=*j j Δ*y*

*∂y ∂t*

<sup>Δ</sup>*y=*j j <sup>Δ</sup>*<sup>y</sup>* (19)

¼ Δ*q* (18)

*Sf y*<sup>3</sup>

¼ Δ 2*=* 3*Cf*

There are several ways to linearize (18), either piecewise or totally, and use it to find the flow and the slide area. This is probably the most popular slide

*Sf* ¼ sin *θ* 1 � tan *ϕ*<sup>0</sup> ð *=* tan *θ*Þ ¼ sin *θ*ð Þ 1 � j j Δ*y=*ð Þ j j tan *θ* (17)

[11], but (11) will apply, with local *y* and *Sf* instead of *y0* and *Sf0*.

the flow in the tail is not as interesting in disaster prevention studies.

**3. Limitations of the theory**

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

**3.1 Limitations of the translatory model**

**3.2 The recessing translatory wave**

*3.2.2 The laminar wave*

Example 1 and (10).

**147**

(8) and (9), flow equations will be

*q* ¼ *vy* ¼ 2*=* 3*Cf*

*∂y ∂t*

*3.2.1 The recessing wave and the flow head*

Here we shall only consider the translatory wave when the wave celerity is equal to the water velocity (*v=c*) and the momentum equation degenerates to the form (3). In this case we get a direct solution. But as the flow spreads out, the equilibrium between surface slope and bed slope will come even quicker than before. We can therefore assume the kinematic wave equation to be valid, a few minutes after the front has passed. The kinematic wave routing equation (valid when changes in velocity head are small) reads [11]

$$d\frac{\partial Q}{\partial t} + c\frac{\partial Q}{\partial r} = d\frac{\partial^2 Q}{\partial r^2}; \quad d = \frac{q}{2S\_{f0}}\tag{15}$$

Using *Q = Q(r - ct)* as before cancels the left-hand terms, and we get the expression *Q = α ryv* where *α* is the opening angle of the funnel.

$$Q = \mathbf{C}\_1(r - ct) + \mathbf{C}\_2 \quad \text{or} \quad y = (a(r - ct) + b)/r \tag{16}$$

In (15) a and b are unknown constants. They are needed because Eq. 11 cannot fulfill the boundary condition *Q=0* when *r – ct=0*. Equation 11 is valid with sufficient accuracy after *y* has reached 90 % of yo or more. Many natural channels have a shape

**Figure 1.** *Chart of the R(Z).*

that can be approximated as a series of flumes with constant width and slowly changing funnels. It is therefore quite clear that a translatory wave can flow down any channel, adapt to the changing environment, and still be a translatory wave.

## **3. Limitations of the theory**

*2.2.4 Change of slope or channel form*

*Landslides - Investigation and Monitoring*

used if the channel changes form.

*2.2.5 Flow in a funnel-shaped flume*

velocity head are small) reads [11]

**Figure 1.** *Chart of the R(Z).*

**146**

If the slope changes, the velocity will adjust to a new depth and velocities in a region just before and just after the slope change as long as the flow is subcritical in the fluid mechanical sense. When the slide flows downhill, the slope will in most cases be diminished. The flow will adapt to the new depth on a short distance.

Equation (5) can be used to find the new depth by using the fact that the flow is the same, upstream and downstream of the slope change. The same approach is

Here we shall only consider the translatory wave when the wave celerity is equal to the water velocity (*v=c*) and the momentum equation degenerates to the form (3). In this case we get a direct solution. But as the flow spreads out, the equilibrium between surface slope and bed slope will come even quicker than before. We can therefore assume the kinematic wave equation to be valid, a few minutes after the front has passed. The kinematic wave routing equation (valid when changes in

*Q*

Using *Q = Q(r - ct)* as before cancels the left-hand terms, and we get the

In (15) a and b are unknown constants. They are needed because Eq. 11 cannot fulfill the boundary condition *Q=0* when *r – ct=0*. Equation 11 is valid with sufficient accuracy after *y* has reached 90 % of yo or more. Many natural channels have a shape

*<sup>∂</sup>r*<sup>2</sup> ; *<sup>d</sup>* <sup>¼</sup> *<sup>q</sup>*

*Q* ¼ *C*1ð Þþ *r* � *ct C*<sup>2</sup> *or y* ¼ ð Þ *a r*ð Þþ � *ct b =r* (16)

2*Sf*<sup>0</sup>

(15)

*<sup>∂</sup><sup>r</sup>* <sup>¼</sup> 0; *<sup>q</sup>* <sup>¼</sup> *yv* (14)

The continuity equation in a channel shaped as a funnel flume is

*∂y ∂t* þ 1 *r <sup>∂</sup>*ð Þ *rq*

*∂Q ∂t* þ *c ∂Q <sup>∂</sup><sup>r</sup>* <sup>¼</sup> *<sup>d</sup> <sup>∂</sup>*<sup>2</sup>

expression *Q = α ryv* where *α* is the opening angle of the funnel.

## **3.1 Limitations of the translatory model**

The most severe limitation of the translatory model is the very heart of the theory, the progressive wave. The wave starts at some point; for a while there will be either increasing or diminishing flow in this point. While we have increasing or constant flow, we are in the head of the wave. For constant flow (11) and (12) apply, but a flow increase will come as a new translatory wave riding upon the first [11], but (11) will apply, with local *y* and *Sf* instead of *y0* and *Sf0*.

## **3.2 The recessing translatory wave**

## *3.2.1 The recessing wave and the flow head*

Diminishing flow is more complicated. Then we are in the tail of the wave head, the wave is recessing, and it is no longer progressing with a constant speed. If it is necessary to calculate this flow, numerical simulation using the full Saint-Venant equations, modified according to Section 2.2.3, may be tried. This involves more advanced programming than for a water waves and has not been tried for landslides using nonlinear friction. In [13] a very advanced model of this kind is presented, using linear friction, but allowing different rheologies, both viscous and non-Newtonian flow, which should give the same possibilities, as nonlinear friction coupled with a slope factor can be used as a tool to imitate non-Newtonian behavior.

The head of the wave is the most damaging part, and the head also determines the length of the damaged area. Therefore, the head flow is the most important, and the flow in the tail is not as interesting in disaster prevention studies.

## *3.2.2 The laminar wave*

If it is clear from the beginning that the velocity is low enough, so all inertial and nonlinear terms can be neglected, and the gradient flow model can be applied in a way like [6, 7]. Using a notation similar to theirs, but using *Δy* instead of *ΔH* (with *Δ* meaning the gradient/divergence operator), we would have a slope factor as in Example 1 and (10).

$$S\_f = \sin\theta(1 - \tan\phi'/\tan\theta) = \sin\theta(1 - |\Delta y/(|\tan\theta|)|) \tag{17}$$

Ѳ is the bottom slope in the direction of the flow. Using the laminar model from (8) and (9), flow equations will be

$$
\overline{q} = \overline{v}\overline{\nu} = 2/(\Im \mathcal{C}\_f) \mathcal{S}\_f y^3 \Delta y / |\Delta y| \qquad \frac{\partial \overline{\nu}}{\partial t} = \Delta \overline{q} \tag{18}
$$

$$\frac{\partial \mathbf{y}}{\partial t} = \Delta \left( 2 / (\mathbf{\hat{3C}\_f}) \mathbf{S}\_f \mathbf{y}^3 \Delta \mathbf{y} / |\Delta \mathbf{y}| \right) \tag{19}$$

There are several ways to linearize (18), either piecewise or totally, and use it to find the flow and the slide area. This is probably the most popular slide model, but it demands a linear relationship between the flow resistance *(τ τcr*) in (8) and the surface gradient.

In field investigations of actual slides, it may be difficult to find the actual investigations begin after the slide stops. Then the total flow of the slide is an

consider Eq. (11) at time t = 0; the tip of the slide is in x = 0. Let x < 0 be any

*<sup>Q</sup>* <sup>∗</sup> <sup>¼</sup> *ZL* � *R Z*ð Þ*=Sf* ; *R Z*ð Þ¼ <sup>1</sup> <sup>þ</sup> ð Þ <sup>1</sup> <sup>þ</sup> *<sup>Z</sup>* <sup>2</sup> � ð Þ <sup>1</sup> � *<sup>Z</sup>* ln 1ð Þ � *<sup>Z</sup>* (21)

2

To find the Chezy coefficient *C*, there must be some observations of flow times. If the time tx that it took the slide to move from – *x* in (19) to *x* = 0, we simply have

tions, (20) is a rather robust formula to find *Sf*. As seen from **Figure 1**, Z does not vary very much except in the range 0.5 *y0 < y <* 0.9 *y0*. Outside this range we have 1 or 0.5 for *R(Z)*. If *R(Z)* is estimated to be 0,5–0,6, we are in the *Z* range

This gives the possibility to use an estimated *y0* value in (20) and find *Sf*. Example 1 shows how Sf relates to the hillslope *tanѲ* and the critical shear stress

But even without an estimate for *C*, the estimate for *y0* and (11) give us an

As stated in the introduction, the translatory theory can be used for any landslide flowing in a regular manner, which means that the bottom friction and bed slope control the flow. The velocity of vertical falls, turfs, frozen blocks of soil, ice, or snow is not governed by the friction in the flowing phase only at the onset of motion and when the sliding stops. At high velocities a lot of air can be mixed into the flow; this is especially true for avalanches of dry fine-grained snow. In such cases material properties such as density and viscosity may become indeterminable. For flowing slides, the translatory wave theory can be used for any event, debris flows, or mudflows. In [5–8] slides that this theory can be used for are discussed. But it is also possible to apply it for unusual cases as submarine slides and rock

The main application of the translatory theory is to use the information available

In such cases a danger area has to be estimated for use by the civil protection authorities. They can declare the area off limits (all visits forbidden), close roads, and order areas to be evacuated. It is especially in the case of evacuations and closing roads for all traffic that an overcautious and too large a danger area estimate

of source area material properties and magnitudes, rundown slopes, and flow depths to estimate the flow coefficients *Sf* and *C*. A good database for these coefficients will make it possible to estimate them at sites where a slide danger has been

estimate of how large an area a slide with a known total mass will cover.

position, and the relation for the dimensionless total flow past it is

*Z* ¼ *y=y*<sup>0</sup> and *L* ¼ �*x=y*<sup>0</sup> (20)

*Sf* <sup>¼</sup> *R Z*ð Þ*<sup>=</sup> ZL* � *<sup>Q</sup>* <sup>∗</sup> ð Þ*:* (22)

. When Q is found from field investiga-

important factor. In a point where we have

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

The dimensional value of Q is Q = Q\*y0

ratio *tanφ*´ that also relates to the soil parameters.

avalanches as will be demonstrated in the next chapter.

identified and the slide is expected to fall anytime.

can come under severe criticism.

**149**

*c = x/ tx*. Then (12) gives an estimate for *C*.

**5.2 Discussion of the method**

From this we deduce

0.9 *< y/y0 <* 1.

## **4. The onset of motion**

### **4.1 The initial state**

#### *4.1.1 Development of instability*

Instability occurs when the shear stress in the soil becomes larger than the critical strength, then a landslide may start. The causes of this can be geological, morphological, physical, or human (see [1] or [2]) or any combination of these causes. A misjudgment of the geological or morphological conditions of the site where a slope is being constructed, or altered, may lead to a landslide. This would be listed as a human cause but can have the simple physical reason that the soil could not withstand the stress increase imposed by the construction work.

Instability of this kind can take time to develop. The soil may grow thicker by time, a bulge can develop due to swelling where the horizontal stress is at its maximum, and a slide will start when the horizontal stress exceeds the passive pressure limit. During the following plastic deformation of the soil, the cohesive strength (*c´*) of the soil is diminished, the plastic deformation is accelerated, the instability has set in, and we have the landslide running possibly a debris flow.

The instability can also set in because *c´* is reduced by water uptake of the soil, either by infiltration from increased groundwater pressure or soaking of the soil by rain. This reduces the φ´ at once and sets in the instability when we have *φ´ < Ѳ*.

A regional study of slopes and lithology, and analysis together with other factors, can bring about important statistics concerning landslide susceptibility and probability of occurrence of landslides caused by the above described instability [3]. This may even be done globally using satellite data [10]. Predicting the timing and location of landslide initiated by developing instability of this kind is very difficult, but considerable success is reported in this field [5, 9].

#### *4.1.2 Slope instability by external causes*

External forces can cause slope failure that starts a landslide. These are earthquakes, other tectonic movements, and explosions that rupture the soil and reduce the shear strength of the soil by one quick movement. An instability of this kind depends on the strength of the external force and is even more difficult to predict.

Large earthquakes usually start landslides, both on land and in the sea (submarine slides).

There may be any combination of events that cause landslides. A landslide due to developing instability may be triggered by a small earth tremor. This will cause the landslide to happen sooner than otherwise expected.

## **5. Practical aspects of the translatory model**

#### **5.1 The estimation and use of the Chezy coefficient**

The translatory model is based on the Chezy coefficient that can be relatively stable regionally, but this must be established by extensive field investigations.

In field investigations of actual slides, it may be difficult to find the actual investigations begin after the slide stops. Then the total flow of the slide is an important factor. In a point where we have

$$Z = \mathfrak{y}/\mathfrak{y}\_0 \qquad \text{and} \qquad L = -\mathfrak{x}/\mathfrak{y}\_0 \tag{20}$$

consider Eq. (11) at time t = 0; the tip of the slide is in x = 0. Let x < 0 be any position, and the relation for the dimensionless total flow past it is

$$\mathbf{Q}^{\*} = \mathbf{Z}\mathbf{L} - \mathbf{R}(\mathbf{Z})/\mathbf{S}\_{\mathbf{f}}; \quad \mathbf{R}(\mathbf{Z}) = \left(\mathbf{1} + \left(\mathbf{1} + \mathbf{Z}\right)^{2}\right) - \left(\mathbf{1} - \mathbf{Z}\right)\ln\left(\mathbf{1} - \mathbf{Z}\right) \tag{21}$$

From this we deduce

model, but it demands a linear relationship between the flow resistance *(τ τcr*)

Instability occurs when the shear stress in the soil becomes larger than the critical strength, then a landslide may start. The causes of this can be geological, morphological, physical, or human (see [1] or [2]) or any combination of these causes. A misjudgment of the geological or morphological conditions of the site where a slope is being constructed, or altered, may lead to a landslide. This would be listed as a human cause but can have the simple physical reason that the soil could

Instability of this kind can take time to develop. The soil may grow thicker by

External forces can cause slope failure that starts a landslide. These are earthquakes, other tectonic movements, and explosions that rupture the soil and reduce the shear strength of the soil by one quick movement. An instability of this kind depends on the strength of the external force and is even more difficult to predict. Large earthquakes usually start landslides, both on land and in the sea (subma-

There may be any combination of events that cause landslides. A landslide due to developing instability may be triggered by a small earth tremor. This will cause the

The translatory model is based on the Chezy coefficient that can be relatively stable regionally, but this must be established by extensive field investigations.

time, a bulge can develop due to swelling where the horizontal stress is at its maximum, and a slide will start when the horizontal stress exceeds the passive pressure limit. During the following plastic deformation of the soil, the cohesive strength (*c´*) of the soil is diminished, the plastic deformation is accelerated, the instability has set in, and we have the landslide running possibly a debris flow. The instability can also set in because *c´* is reduced by water uptake of the soil, either by infiltration from increased groundwater pressure or soaking of the soil by rain. This reduces the φ´ at once and sets in the instability when we have *φ´ < Ѳ*. A regional study of slopes and lithology, and analysis together with other factors, can bring about important statistics concerning landslide susceptibility and probability of occurrence of landslides caused by the above described instability [3]. This may even be done globally using satellite data [10]. Predicting the timing and location of landslide initiated by developing instability of this kind is very difficult,

not withstand the stress increase imposed by the construction work.

but considerable success is reported in this field [5, 9].

landslide to happen sooner than otherwise expected.

**5. Practical aspects of the translatory model**

**5.1 The estimation and use of the Chezy coefficient**

*4.1.2 Slope instability by external causes*

rine slides).

**148**

in (8) and the surface gradient.

*Landslides - Investigation and Monitoring*

**4. The onset of motion**

*4.1.1 Development of instability*

**4.1 The initial state**

$$\mathbf{S}\_{\hat{f}} = \mathbf{R}(\mathbf{Z})/(\mathbf{Z}\mathbf{L} - \mathbf{Q}^\*). \tag{22}$$

The dimensional value of Q is Q = Q\*y0 2 . When Q is found from field investigations, (20) is a rather robust formula to find *Sf*. As seen from **Figure 1**, Z does not vary very much except in the range 0.5 *y0 < y <* 0.9 *y0*. Outside this range we have 1 or 0.5 for *R(Z)*. If *R(Z)* is estimated to be 0,5–0,6, we are in the *Z* range 0.9 *< y/y0 <* 1.

This gives the possibility to use an estimated *y0* value in (20) and find *Sf*.

Example 1 shows how Sf relates to the hillslope *tanѲ* and the critical shear stress ratio *tanφ*´ that also relates to the soil parameters.

To find the Chezy coefficient *C*, there must be some observations of flow times. If the time tx that it took the slide to move from – *x* in (19) to *x* = 0, we simply have *c = x/ tx*. Then (12) gives an estimate for *C*.

But even without an estimate for *C*, the estimate for *y0* and (11) give us an estimate of how large an area a slide with a known total mass will cover.

### **5.2 Discussion of the method**

As stated in the introduction, the translatory theory can be used for any landslide flowing in a regular manner, which means that the bottom friction and bed slope control the flow. The velocity of vertical falls, turfs, frozen blocks of soil, ice, or snow is not governed by the friction in the flowing phase only at the onset of motion and when the sliding stops. At high velocities a lot of air can be mixed into the flow; this is especially true for avalanches of dry fine-grained snow. In such cases material properties such as density and viscosity may become indeterminable.

For flowing slides, the translatory wave theory can be used for any event, debris flows, or mudflows. In [5–8] slides that this theory can be used for are discussed. But it is also possible to apply it for unusual cases as submarine slides and rock avalanches as will be demonstrated in the next chapter.

The main application of the translatory theory is to use the information available of source area material properties and magnitudes, rundown slopes, and flow depths to estimate the flow coefficients *Sf* and *C*. A good database for these coefficients will make it possible to estimate them at sites where a slide danger has been identified and the slide is expected to fall anytime.

In such cases a danger area has to be estimated for use by the civil protection authorities. They can declare the area off limits (all visits forbidden), close roads, and order areas to be evacuated. It is especially in the case of evacuations and closing roads for all traffic that an overcautious and too large a danger area estimate can come under severe criticism.

When the size of a danger area is estimated, there has to be a safety margin. In construction it is decided by the safety margin to put on material strength and design loads. Using *Sf* and *C* to find a danger area, the safety margin can be decided in the same way. The most likely values for S*<sup>f</sup>* and *C* have to be found first and the safe values to be used in the estimation decided afterwards.

This shows that in assessing the tsunami risk in the Pacific coastal regions of Japan, the landslide risk must be considered. This fact may result in that considerably larger events than the Tohoku tsunami are possible if a larger slide than this

The assessment of this possibility of larger slides can be difficult; there is considerable uncertainty in estimating the height (*y0*) of possible landslide, but this may be done using drilled exploration wells, like the one shown in Figure 2 in Ref. [14]. The flow time will be the estimated duration of the earthquake, so using the *C* found here, there may be a possibility to estimate the volume moved by a slide

**6.2 Case study 2: the rock avalanche on the Morsárjökull outlet glacier, SE**

On March 20, 2007, a large rock avalanche fell on the Morsárjökull outlet glacier, in the southern part of the Vatnajökull ice cap, in Southeast Iceland (**Figure 2**). This rock avalanche is considered to be one of the largest rock avalanches which have

.

The Morsárjökull glacier is a small outlet glacier on the southern side of the Vatnajökull ice cap in Southeastern Iceland, located between the Þorsteinshöfði Mountain (1343 m) and Skarðatindur Mountain (1385 m) to the east and the Austurhnúta peak (1084 m) in the Miðfell Mountain (1430 m) to the west (**Figure 2**). The glacier is about 4 km long and 750–900 m wide in the uppermost part and around 800–900 wide in the lower part. Its surface is around 3.3 km<sup>2</sup>

The glacier is surrounded by steep mountain slopes and is categorized as an icefall glacier with two ice falls. In the uppermost part of the glacier, below the ice falls, the glacier surface is steep, but below the icefalls the surface inclination is around 85/

*Location map of the Morsárjökull outlet glacier in SE Iceland. The fracture zone and the accumulation lobe marked with gray color. Legend: (1) 20 m contour lines on land, (2) 20 m contour lines on glacier, (3)*

*proglacial lake, (4) medial moraine, and (5) equilibrium line (ELA).*

using (20)–(21) and thereby the strength of the corresponding tsunami.

8-m-thick slide is released.

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

**Iceland**

**Figure 2.**

**151**

*6.2.1 Location and event*

occurred in Iceland during the last decades [18].

## **6. Case studies**

### **6.1 Case study 1: the Tohoku submarine slide**

#### *6.1.1 Location and event*

The most famous tsunami event of recent years hit Honshu, the big island in Japan, on Friday 11 March 2011 at 05:46 UTC. Its location was off the Pacific coast of Tohoku, Japan. The earthquake event and the devastation caused by this tsunami is very well documented; here we shall follow [14].

The cause was a Mw 9.0 (magnitude moment) undersea megathrust earthquake. The epicenter was approximately 70 km east of the Oshika Peninsula of Tohoku with the hypocenter in approximately 30 km deep water.

As most earthquakes above Mw 7.0 do, it caused a tsunami. The tsunami was much bigger than what could be expected and caused a massive devastation. Tens of thousands of people suffered loss of life or all their possessions. Houses and other structures were swept away from large areas. The town of Sendai was very badly hit and 13 other Japanese cities.

The most serious incident, however, was a meltdown in the Fukushima Daiichi Nuclear Power Plant, when the tsunami proved too high for the coastal defenses and flooded the emergency generators.

#### *6.1.2 Properties and data*

In [14] it is suggested that the bottom deformations that caused the very strong Tohoku tsunami in the Pacific Ocean can be explained by a submarine landslide. This suggests that the coseismic slip of the earthquake triggers a sliding of the surface sediments. In combination they cause the bottom deformation.

Finally, it can thus be concluded that the coseismic slip and the landslide are both responsible for the Tohoku tsunami in March 2011, not the coseismic slip alone.

From the data obtained in the exploratory drilling at the site shown in Figure 2 in Ref. [14], it was concluded that the tsunami was caused by the mass movement shown in Figure 3 in Ref. [14]. The slide triggered by the earthquake amplifies the tsunami.

This landslide is not very high compared to what has happened elsewhere. But it is quite possible that the coastal defenses of Fukushima would have saved the emergency generators and prevented the nuclear disaster had the landslide not happened.

The main data for the slide were originally deduced from [15–17]. In [14] we find the following:

Bottom slope is 9 %, *c* = 2 m/s, and *y0* = 8 m. The length of the slide is 117.5 km, the width of the slide is 110 km, but the length of the flow path is only 60 m, and the flow time is thus 30 s. The volume moved by the slide is 96.8 km<sup>3</sup> .

#### *6.1.3 Evaluation*

Using **Table 1** and the bottom slope, we find a *C* value of 0.32. Corrected for *Sf* ≈ 0.1, we get the value of the Chezy coefficient C = 0.71.

## *The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

This shows that in assessing the tsunami risk in the Pacific coastal regions of Japan, the landslide risk must be considered. This fact may result in that considerably larger events than the Tohoku tsunami are possible if a larger slide than this 8-m-thick slide is released.

The assessment of this possibility of larger slides can be difficult; there is considerable uncertainty in estimating the height (*y0*) of possible landslide, but this may be done using drilled exploration wells, like the one shown in Figure 2 in Ref. [14]. The flow time will be the estimated duration of the earthquake, so using the *C* found here, there may be a possibility to estimate the volume moved by a slide using (20)–(21) and thereby the strength of the corresponding tsunami.

## **6.2 Case study 2: the rock avalanche on the Morsárjökull outlet glacier, SE Iceland**

## *6.2.1 Location and event*

When the size of a danger area is estimated, there has to be a safety margin. In construction it is decided by the safety margin to put on material strength and design loads. Using *Sf* and *C* to find a danger area, the safety margin can be decided in the same way. The most likely values for S*<sup>f</sup>* and *C* have to be found first and the

The most famous tsunami event of recent years hit Honshu, the big island in Japan, on Friday 11 March 2011 at 05:46 UTC. Its location was off the Pacific coast of Tohoku, Japan. The earthquake event and the devastation caused by this tsunami

The cause was a Mw 9.0 (magnitude moment) undersea megathrust earthquake. The epicenter was approximately 70 km east of the Oshika Peninsula of Tohoku

As most earthquakes above Mw 7.0 do, it caused a tsunami. The tsunami was much bigger than what could be expected and caused a massive devastation. Tens of thousands of people suffered loss of life or all their possessions. Houses and other structures were swept away from large areas. The town of Sendai was very badly hit

The most serious incident, however, was a meltdown in the Fukushima Daiichi Nuclear Power Plant, when the tsunami proved too high for the coastal defenses

In [14] it is suggested that the bottom deformations that caused the very strong Tohoku tsunami in the Pacific Ocean can be explained by a submarine landslide. This suggests that the coseismic slip of the earthquake triggers a sliding of the surface sediments. In combination they cause the bottom deformation.

Finally, it can thus be concluded that the coseismic slip and the landslide are both

From the data obtained in the exploratory drilling at the site shown in Figure 2 in Ref. [14], it was concluded that the tsunami was caused by the mass movement shown in Figure 3 in Ref. [14]. The slide triggered by the earthquake amplifies the tsunami. This landslide is not very high compared to what has happened elsewhere. But it is quite possible that the coastal defenses of Fukushima would have saved the emergency generators and prevented the nuclear disaster had the landslide not happened. The main data for the slide were originally deduced from [15–17]. In [14] we

Bottom slope is 9 %, *c* = 2 m/s, and *y0* = 8 m. The length of the slide is 117.5 km, the width of the slide is 110 km, but the length of the flow path is only 60 m, and the

Using **Table 1** and the bottom slope, we find a *C* value of 0.32. Corrected for

.

flow time is thus 30 s. The volume moved by the slide is 96.8 km<sup>3</sup>

*Sf* ≈ 0.1, we get the value of the Chezy coefficient C = 0.71.

responsible for the Tohoku tsunami in March 2011, not the coseismic slip alone.

safe values to be used in the estimation decided afterwards.

**6.1 Case study 1: the Tohoku submarine slide**

*Landslides - Investigation and Monitoring*

is very well documented; here we shall follow [14].

with the hypocenter in approximately 30 km deep water.

**6. Case studies**

*6.1.1 Location and event*

and 13 other Japanese cities.

*6.1.2 Properties and data*

find the following:

*6.1.3 Evaluation*

**150**

and flooded the emergency generators.

On March 20, 2007, a large rock avalanche fell on the Morsárjökull outlet glacier, in the southern part of the Vatnajökull ice cap, in Southeast Iceland (**Figure 2**). This rock avalanche is considered to be one of the largest rock avalanches which have occurred in Iceland during the last decades [18].

The Morsárjökull glacier is a small outlet glacier on the southern side of the Vatnajökull ice cap in Southeastern Iceland, located between the Þorsteinshöfði Mountain (1343 m) and Skarðatindur Mountain (1385 m) to the east and the Austurhnúta peak (1084 m) in the Miðfell Mountain (1430 m) to the west (**Figure 2**). The glacier is about 4 km long and 750–900 m wide in the uppermost part and around 800–900 wide in the lower part. Its surface is around 3.3 km<sup>2</sup> . The glacier is surrounded by steep mountain slopes and is categorized as an icefall glacier with two ice falls. In the uppermost part of the glacier, below the ice falls, the glacier surface is steep, but below the icefalls the surface inclination is around 85/

#### **Figure 2.**

*Location map of the Morsárjökull outlet glacier in SE Iceland. The fracture zone and the accumulation lobe marked with gray color. Legend: (1) 20 m contour lines on land, (2) 20 m contour lines on glacier, (3) proglacial lake, (4) medial moraine, and (5) equilibrium line (ELA).*

1000 m. A proglacial glacier lake occurs in front of the eastern side of the glacier snout. The lake, which is today around 1.8 km long and 350–400 m wide, begun to form around 1945. This lake has increased considerably in size during the last decades. About 1.6 to 1.8 km from the glacier snout, a system of end moraines occurs crossing the valley floor, marking the maximum extent of the glacier around 1910. The glacier is drained by the Morsá Glacier River, which merged into the old Skeiðará River gorge approximately 12–13 km down glacier.

the debris lobe, approximately 800 m from the scar. On the surface of the accumulation lobe flow line structures and longitudinal ridges, up to 3–4 m high were observed. Similar phenomena have been described from other rock avalanches on glaciers, e.g., the Sherman glacier in Alaska [20] and the Little Tahoma Peak glacier

It is evident that the Morsárjökull glacier has retreated considerably during the last century, and during the last decade, the melting has been very rapid. Therefore, it is thought that undercutting of the mountain slope by glacial erosion and the retreat of the glacier are the main contributing factors leading to the rock avalanche. The glacial erosion has destabilized the slope, which is chiefly composed of palagonite and dolerite rocks, affected by geothermal alteration. Hence a subsequent fracture formation has weakened the bedrock. The exact triggering factor, however, is not known. No seismic activity or meteorological signal, which could be interpreted as triggering factors, such as heavy rainfall or intensive snowmelt, was recorded prior to the rock avalanche.

A quick calculation using (11) sh0ws that there is very little difference between the average thickness and *y0*. The upper section is around 500 m long, making the lower section 1000 m. The glacial surface just below where the slide stops inclines

It is necessary to estimate the flow time of the slide in order to find the Chezy coefficient *C*; this is found from a tremor recorded in a nearby earthquake station to be 100 seconds, but this estimate is very uncertain. However, it affects only *c* and *C*. **Table 2** shows the data and the calculation of the main parameters of the slide.

The main results show similar values for the upper and lower *C*´s. Note the large difference from the underwater slide in Case 1. Had the *C* values been known, the

*y0* Data 5 7 m *Ѳ Arctan(Sfo)* 9.0 5.1 Deg *Sfo* Data 0.16 0.09 NA *Sf sinѲ (1 tanφ´/tanѲ))* 0.10 0.03 NA Time 100 totals 25 75 sec Length 1500 total 500 1000 m *c* Length/time 20 13.3 m/s *C C = c/ (g y0 Sf)* 8.8 8.7 NA

**Upper section Lower section Source Value Value Dim**

**6.3 Case study 3: the debris slide in the Móafellshyrna Mountain in the**

On September 20, 2012, a large debris slide occurred at the northwest facing slope of the Móafellshyrna Mountain in the Fljótin area situated in the Tröllaskagi

**Tröllaskagi peninsula, central north Iceland**

*6.3.1 Location and event*

*Upper section of the Morsárjökull slide.*

**Table 2.**

**153**

in the Washington state [21].

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

*6.2.3 Evaluation*

by 5.6 %; this is put as *tanϕ*'.

time could have been computed.

It is believed that the rock avalanche fell in two separate stages; the main part fell on March 20 and the second and smaller one on April 17, 2007. In both instances surface tremors were detected on seismographs, and on April 17 workers at the Skeiðarársandur sandur plain heard a rumbling noise from the glacier.

The rock avalanche was initiated from the north face of the headwall above the uppermost part of the glacier (**Figure 3**). The scar in the rock-face is around 330 m high, reaching from about 620 m up to 950 m and about 480 m wide on average, showing that the main part of the headwall collapsed. From a digital elevation model of the headwall, based on aerial photographs from 2003 and 2007, the total volume of debris is estimated to be around 4 to 4.5 million m<sup>3</sup> or about 10–12 million tons. The rock avalanche debris fell first straight onto the glacier and then turned 90° downward. No indication of any deformation or erosion of the glacier surface was observed, which is similar to the description of the rock avalanche which fell on the Sherman glacier in Alaska in 1964 [19, 20].

## *6.2.2 Properties and data*

The accumulation lobe is on average 1,5 km long, with a maximum length of 1.6 km, reaching from 520 m a.s.l. down to about 350 m a.s.l. Its width is from 125 m to 650 m or on average 480 m. The total area which the lobe covers is about 720.000 m<sup>2</sup> which was approximately 1/5 of the glacier surface in 2007. The thickness of the debris lobe was measured more than 8 m in places, but its mean thickness was somewhere between 5.5 and 6.2 m. In 2007 the debris mass is coarse-grained and boulder-rich with little fine material. Blocks over 5 to 8 m in diameter are common on the edges of the lobe up to 1.6 km from the source. The largest block which was observed in the accumulation lobe is around 800 tons located on the western side of

#### **Figure 3.**

*The headwall above the eastern icefall in the uppermost part of the Morsárjökull outlet glacier and the accumulation lobe (Photo Þorsteinn Sæmundsson 2007).*

## *The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

the debris lobe, approximately 800 m from the scar. On the surface of the accumulation lobe flow line structures and longitudinal ridges, up to 3–4 m high were observed. Similar phenomena have been described from other rock avalanches on glaciers, e.g., the Sherman glacier in Alaska [20] and the Little Tahoma Peak glacier in the Washington state [21].

It is evident that the Morsárjökull glacier has retreated considerably during the last century, and during the last decade, the melting has been very rapid. Therefore, it is thought that undercutting of the mountain slope by glacial erosion and the retreat of the glacier are the main contributing factors leading to the rock avalanche. The glacial erosion has destabilized the slope, which is chiefly composed of palagonite and dolerite rocks, affected by geothermal alteration. Hence a subsequent fracture formation has weakened the bedrock. The exact triggering factor, however, is not known. No seismic activity or meteorological signal, which could be interpreted as triggering factors, such as heavy rainfall or intensive snowmelt, was recorded prior to the rock avalanche.

## *6.2.3 Evaluation*

1000 m. A proglacial glacier lake occurs in front of the eastern side of the glacier snout. The lake, which is today around 1.8 km long and 350–400 m wide, begun to form around 1945. This lake has increased considerably in size during the last decades. About 1.6 to 1.8 km from the glacier snout, a system of end moraines occurs crossing the valley floor, marking the maximum extent of the glacier around 1910. The glacier is drained by the Morsá Glacier River, which merged into the old

It is believed that the rock avalanche fell in two separate stages; the main part fell on March 20 and the second and smaller one on April 17, 2007. In both instances surface tremors were detected on seismographs, and on April 17 workers at the

The rock avalanche was initiated from the north face of the headwall above the uppermost part of the glacier (**Figure 3**). The scar in the rock-face is around 330 m high, reaching from about 620 m up to 950 m and about 480 m wide on average, showing that the main part of the headwall collapsed. From a digital elevation model of the headwall, based on aerial photographs from 2003 and 2007, the total volume of debris is estimated to be around 4 to 4.5 million m<sup>3</sup> or about 10–12 million tons. The rock avalanche debris fell first straight onto the glacier and then turned 90° downward. No indication of any deformation or erosion of the glacier surface was observed, which is similar to the description of the rock avalanche

The accumulation lobe is on average 1,5 km long, with a maximum length of 1.6 km, reaching from 520 m a.s.l. down to about 350 m a.s.l. Its width is from 125 m to 650 m or on average 480 m. The total area which the lobe covers is about 720.000 m<sup>2</sup> which was approximately 1/5 of the glacier surface in 2007. The thickness of the debris lobe was measured more than 8 m in places, but its mean thickness was somewhere between 5.5 and 6.2 m. In 2007 the debris mass is coarse-grained and boulder-rich with little fine material. Blocks over 5 to 8 m in diameter are common on the edges of the lobe up to 1.6 km from the source. The largest block which was observed in the accumulation lobe is around 800 tons located on the western side of

*The headwall above the eastern icefall in the uppermost part of the Morsárjökull outlet glacier and the*

*accumulation lobe (Photo Þorsteinn Sæmundsson 2007).*

Skeiðará River gorge approximately 12–13 km down glacier.

*Landslides - Investigation and Monitoring*

which fell on the Sherman glacier in Alaska in 1964 [19, 20].

*6.2.2 Properties and data*

**Figure 3.**

**152**

Skeiðarársandur sandur plain heard a rumbling noise from the glacier.

A quick calculation using (11) sh0ws that there is very little difference between the average thickness and *y0*. The upper section is around 500 m long, making the lower section 1000 m. The glacial surface just below where the slide stops inclines by 5.6 %; this is put as *tanϕ*'.

It is necessary to estimate the flow time of the slide in order to find the Chezy coefficient *C*; this is found from a tremor recorded in a nearby earthquake station to be 100 seconds, but this estimate is very uncertain. However, it affects only *c* and *C*. **Table 2** shows the data and the calculation of the main parameters of the slide.

The main results show similar values for the upper and lower *C*´s. Note the large difference from the underwater slide in Case 1. Had the *C* values been known, the time could have been computed.


#### **Table 2.**

*Upper section of the Morsárjökull slide.*

## **6.3 Case study 3: the debris slide in the Móafellshyrna Mountain in the Tröllaskagi peninsula, central north Iceland**

#### *6.3.1 Location and event*

On September 20, 2012, a large debris slide occurred at the northwest facing slope of the Móafellshyrna Mountain in the Fljótin area situated in the Tröllaskagi

## *Landslides - Investigation and Monitoring*

peninsula, an area in central north Iceland. The local residents of the farm Þrasastaðir witnessed the release of the landslide. The area was heavily eroded during the last glaciation period Iceland. **Figure 4**, modified from [22], shows the site, the location, and an aerial picture of the extent of the slide. The onset of the slide is traced to the special geology of the mountain (**Figure 4**) and thawing of permafrost in the sediments at an approximate elevation of 900 m (meters elevation above sea level) at the top of the Móafellshyrna Mountain. The thawing was proven to be the primary reason for the weakening of the sediments in the source area, but the triggering of the event was a combination of heavy rains and three earthquakes of Richter magnitude 4.2–4.5 on September 19 and 20.

## *6.3.2 Properties and data*

The slide originates (at the headscarp) at an elevation of 880 m and runs all the way down to 330 m (**Figures 4** and **5**). The total debris volume is estimated to be 312,000–480,000 m3 [22].

There are three particularly interesting problems associated with this slide. Firstly, in it are large blocks of frozen sedimentary rocks; the biggest one is 12 meters wide and 4 meters high (**Figure 6**). This large block is an integral part of the debris flow, it is carried almost halfway down the mountain, and it stops at the elevation 495–505 m and has carved an erosion channel into the mountainside, shown as black streak in the talus (rock mass at the foot of a cliff) slope in **Figure 5**, but it stops at around 500 m a.s.l.

Secondly, the massive frozen blocks inside the main debris flow will affect the flow velocity and the velocity distribution within the flow. The erosion channel

#### **Figure 4.**

*The geographic setting of the Móafellshyrna site,Tröllaskagi peninsula, central northern Iceland. (A) Hillshaded digital elevation model and contours (in green, meters above sea level) of the Móafellshyrna region. Elevation data are from EU-DEM from GMES RDA. (B) Aerial photograph of the Móafellshyrna site taken before the slide in 2012 from the website Samsyn. (C) Oblique photo of the Móafellshyrna debris slide in July 2015 taken by Costanza Morino. In panels B and C, the perimeter of the slide is marked with red dashed line taken from the trimble data for the deposits and reconstructed from photographs for the upper part; in C this line has been manually traced onto the oblique image. Modified from [22].*

shows that the part of the slide flowing in it has lower velocity than the rest of the flow, which means that the main flow is providing a part of its momentum to the blocks. This may result in a thicker slide that does not run as long out as a slide of

*Landslide debris and an ice-cemented block resting at 495–505 m a.s.l. The block is around 12 m wide and 4 m*

*The Móafellshyrna landslide (photo taken few hours after the slide on September 20, 2012). (A) The topographical profile of the Móafellshyrna landslide. (B) Distribution map of the landslide (photo by Þorsteinn*

Thirdly, the blocks run over a 100-meter-high cliff below the source area at 690. Here they pick up momentum and carve the channel before they stop, 170 m higher up than the rest of the slide which runs about 700 meters farther out (**Figure 5**). The force of the flow on the blocks is unable to transport them any further in spite

uniform material would have done.

*high (photo by Þorsteinn Sæmundsson September 29, 2012).*

**Figure 5.**

**Figure 6.**

**155**

*Sæmundsson 2012). Modified from [22, 23].*

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

#### **Figure 5.**

peninsula, an area in central north Iceland. The local residents of the farm Þrasastaðir witnessed the release of the landslide. The area was heavily eroded during the last glaciation period Iceland. **Figure 4**, modified from [22], shows the site, the location, and an aerial picture of the extent of the slide. The onset of the slide is traced to the special geology of the mountain (**Figure 4**) and thawing of permafrost in the sediments at an approximate elevation of 900 m (meters elevation above sea level) at the top of the Móafellshyrna Mountain. The thawing was proven to be the primary reason for the weakening of the sediments in the source area, but the triggering of the event was a combination of heavy rains and three

earthquakes of Richter magnitude 4.2–4.5 on September 19 and 20.

The slide originates (at the headscarp) at an elevation of 880 m and runs all the way down to 330 m (**Figures 4** and **5**). The total debris volume is estimated to be

There are three particularly interesting problems associated with this slide. Firstly, in it are large blocks of frozen sedimentary rocks; the biggest one is 12 meters wide and 4 meters high (**Figure 6**). This large block is an integral part of the debris flow, it is carried almost halfway down the mountain, and it stops at the elevation 495–505 m and has carved an erosion channel into the mountainside, shown as black streak in the talus (rock mass at the foot of a cliff) slope in **Figure 5**,

Secondly, the massive frozen blocks inside the main debris flow will affect the flow velocity and the velocity distribution within the flow. The erosion channel

*The geographic setting of the Móafellshyrna site,Tröllaskagi peninsula, central northern Iceland. (A) Hillshaded digital elevation model and contours (in green, meters above sea level) of the Móafellshyrna region. Elevation data are from EU-DEM from GMES RDA. (B) Aerial photograph of the Móafellshyrna site taken before the slide in 2012 from the website Samsyn. (C) Oblique photo of the Móafellshyrna debris slide in July 2015 taken by Costanza Morino. In panels B and C, the perimeter of the slide is marked with red dashed line taken from the trimble data for the deposits and reconstructed from photographs for the upper part; in C this*

*line has been manually traced onto the oblique image. Modified from [22].*

*6.3.2 Properties and data*

*Landslides - Investigation and Monitoring*

312,000–480,000 m3 [22].

but it stops at around 500 m a.s.l.

**Figure 4.**

**154**

*The Móafellshyrna landslide (photo taken few hours after the slide on September 20, 2012). (A) The topographical profile of the Móafellshyrna landslide. (B) Distribution map of the landslide (photo by Þorsteinn Sæmundsson 2012). Modified from [22, 23].*

## **Figure 6.**

*Landslide debris and an ice-cemented block resting at 495–505 m a.s.l. The block is around 12 m wide and 4 m high (photo by Þorsteinn Sæmundsson September 29, 2012).*

shows that the part of the slide flowing in it has lower velocity than the rest of the flow, which means that the main flow is providing a part of its momentum to the blocks. This may result in a thicker slide that does not run as long out as a slide of uniform material would have done.

Thirdly, the blocks run over a 100-meter-high cliff below the source area at 690. Here they pick up momentum and carve the channel before they stop, 170 m higher up than the rest of the slide which runs about 700 meters farther out (**Figure 5**). The force of the flow on the blocks is unable to transport them any further in spite

## *Landslides - Investigation and Monitoring*

of the momentum gain. This acts as flow resistance on the part of the slide flowing by and may be a significant factor in determining the flow properties down streams, e.g., that it stops where it does and does not flow any further.

**Figure 5A** and **B** explains the properties of this very complicated slide visually.

The following description of the slide shows where erosion and deposition occurred in the slope (**Figure 5**) [22, 23].


heard in the nearby farm. The velocity from 400 to 330 m is still high; it takes a

*y0* 2/3 of the max. value 6 6 m Width **Figure 5B** 80 130 m *Ѳ Arctan(Sfo*) 19 10 Deg *Sfo* **Figure 5A** 0,33 0,17 NA *Sf sinѲ (1 - tanφ´/tanѲ))* 0,34 0,17 NA Time Length/*c*. Total flow 28 8,5 25,5 sec Length **Figure 5A** 263 360 m *c c= C (g y0 Sf)* 31 14 m/s *C* **Table 2** 8,8 8,7 NA

**495–400 m 400–330 m Source Value Value Dim**

The translatory wave theory is an old tool from the science of hydrodynamics that can be adapted to landslides. However, it is not possible to say that it may be used for certain classes or groups of landslides in the geological classification system of landslides. The theory can be used for all slides that flow downhill with a flow governed by gravitational forces, where the slide has a wavelike front and the slide

The translatory wave equation can be modified to describe this form of a landslide by introducing a new parameter, the slope factor *Sf* into the flow equations as to control the excess shear force in the landslide and its relation to the actual shear strength of the flowing mass. The slope factor controls how great a part of the actual shear force is in excess of the shear strength and thereby the velocity of the flow. The slope factor has to be assumed piecewise constant, and the flow stops when it becomes zero. This model can give a very good picture of the flow and can save valuable time and effort, when a more detailed numerical treatment is necessary. The translatory model is used for moving slides; it cannot predict the onset of a slide. But when in motion, it can predict the length and thickness of the slide from the knowledge of the ground slope and the soil parameters, the Chezy coefficient, and the slope factor. And vice versa, for a mapped slide, it can predict the soil

There are three studied cases presented, first the famous Tohoku tsunami in Japan, which became very strong because of the earthquake that started the tsunami and a landslide that amplified the tsunami. Both Chezy's coefficient and the slope

The second case study is a landslide down the Morsárjökull outlet glacier in Iceland, starting as a rock slide, but it turns downhill into a glacier and turns into a flow slide, in clear contrast to the Tohoku submarine slide. This slide flows first down a steep upper slope and then down a milder lower slope and then finally stops on an even milder slope with a flow factor zero. The results for the upper and lower sections show different slope factors, the lower one belonging to the lower section. The lower thickness value belongs to the upper section, but the Chezy coefficients

slope of 70 meters to stop it.

*Two sections in the Móafellshyrna slide.*

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

follows as a frontal wave with the same velocity.

**7. Conclusion**

**Table 3.**

parameters.

**157**

factor are found.


## *6.3.3 Evaluation*

The slide may be translatory from the resting position of the big block at around 500 m a.s.l. until it stops. The section around it is then the starting point of the downstream flow, and its initial conditions are the time history of the velocity and flow depth in that section.

From **Figure 5** we see that the slide stops at the bottom slope of *tanφ*´= 50/400 = 0,125 or *φ´*= 7°. The width of the slide can be approximately estimated from **Figures 4** and **5**. According to Section 6.2.2, the total flow mass to run from the talus slope from 700 to approximately 500 m a.s.l. is estimated to be 192,000–300,000 m3 . This provides an estimate that the flow channels from 500 down to 330 m has an average depth of 6 m. If the flow channels are supposed to have a parabolic shape, this corresponds to 2/3 of their maximum depth. The profiles and location of these channels are shown in **Figures 5A** and **B**. By assuming the same Chezy coefficients as in **Table 2**, we can calculate **Table 3**.

The results indicate that the translatory theory can model the slide quite well. The flow velocity *c* is found and used to calculate the flow times from 505 to 400 m height. It takes the flow 25,5 s to reach the 330 m elevation but only 8,5 s to run from 495 to 400 m. However, it takes the total mass 28 seconds (total flow 28) to flow from 495 to 400 m. This is longer time than the run from 400 to 330 m, so some material will be left between 495 and 400 m when the downstream end stops, making some accumulation to take at that height interval.

This is an example of computing slide properties from known *C* values in a fast and effective manner. The velocity *c* in the steeper part of the lower part of the slide from 495 to 400 m is high, certainly high enough to produce the rumbling noise

## *The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*


### **Table 3.**

of the momentum gain. This acts as flow resistance on the part of the slide flowing by and may be a significant factor in determining the flow properties down streams,

**Figure 5A** and **B** explains the properties of this very complicated slide visually. The following description of the slide shows where erosion and deposition

• The debris slide was detached from a main scarp at an elevation of 880 m a.s.l. The source material of the landslide was a colluvial cone. Part of the cones front, slides off a steep 100 m high rockwall, at around 800 m a.s.l. It is estimated that the volume of deposit which broke the frontal part of the

• Below the rockwall at around 700 m a.s.l., a 300-m-long, 80–100-m-wide, and 8–10-m-deep channel was carved into the 200-m-high talus slope deposits, entraining around 192,000–300,000 m<sup>3</sup> of material. The average slope of the

• At 495–505 m a.s.l. at the foot of the talus slope, some of the landslide deposits

• The rest of the landslide material traversed down the mountain slope from 495 m and finally came to a rest at 330 m. The uppermost 100 m of the slope are steeper with an averaged slope angle of around 20–25° and in the lowermost 70

The slide may be translatory from the resting position of the big block at around

From **Figure 5** we see that the slide stops at the bottom slope of *tanφ*´= 50/400 =

. This provides an estimate that the flow channels from 500 down to 330 m has an average depth of 6 m. If the flow channels are supposed to have a parabolic shape, this corresponds to 2/3 of their maximum depth. The profiles and location of these channels are shown in **Figures 5A** and **B**. By assuming the same Chezy

The results indicate that the translatory theory can model the slide quite well. The flow velocity *c* is found and used to calculate the flow times from 505 to 400 m height. It takes the flow 25,5 s to reach the 330 m elevation but only 8,5 s to run from 495 to 400 m. However, it takes the total mass 28 seconds (total flow 28) to flow from 495 to 400 m. This is longer time than the run from 400 to 330 m, so some material will be left between 495 and 400 m when the downstream end stops,

This is an example of computing slide properties from known *C* values in a fast and effective manner. The velocity *c* in the steeper part of the lower part of the slide from 495 to 400 m is high, certainly high enough to produce the rumbling noise

500 m a.s.l. until it stops. The section around it is then the starting point of the downstream flow, and its initial conditions are the time history of the velocity and

0,125 or *φ´*= 7°. The width of the slide can be approximately estimated from **Figures 4** and **5**. According to Section 6.2.2, the total flow mass to run from the talus slope from 700 to approximately 500 m a.s.l. is estimated to be 192,000–300,000

came to a rest, where the average slope angle is less than 10°.

. The headscarp has an average

e.g., that it stops where it does and does not flow any further.

colluvial cone was around 120,000–180,000 m<sup>3</sup>

occurred in the slope (**Figure 5**) [22, 23].

*Landslides - Investigation and Monitoring*

slope of around 30–35°.

talus is around 30°.

*6.3.3 Evaluation*

m3

**156**

flow depth in that section.

m of the slope around 10–15° angle.

coefficients as in **Table 2**, we can calculate **Table 3**.

making some accumulation to take at that height interval.

*Two sections in the Móafellshyrna slide.*

heard in the nearby farm. The velocity from 400 to 330 m is still high; it takes a slope of 70 meters to stop it.

## **7. Conclusion**

The translatory wave theory is an old tool from the science of hydrodynamics that can be adapted to landslides. However, it is not possible to say that it may be used for certain classes or groups of landslides in the geological classification system of landslides. The theory can be used for all slides that flow downhill with a flow governed by gravitational forces, where the slide has a wavelike front and the slide follows as a frontal wave with the same velocity.

The translatory wave equation can be modified to describe this form of a landslide by introducing a new parameter, the slope factor *Sf* into the flow equations as to control the excess shear force in the landslide and its relation to the actual shear strength of the flowing mass. The slope factor controls how great a part of the actual shear force is in excess of the shear strength and thereby the velocity of the flow.

The slope factor has to be assumed piecewise constant, and the flow stops when it becomes zero. This model can give a very good picture of the flow and can save valuable time and effort, when a more detailed numerical treatment is necessary.

The translatory model is used for moving slides; it cannot predict the onset of a slide. But when in motion, it can predict the length and thickness of the slide from the knowledge of the ground slope and the soil parameters, the Chezy coefficient, and the slope factor. And vice versa, for a mapped slide, it can predict the soil parameters.

There are three studied cases presented, first the famous Tohoku tsunami in Japan, which became very strong because of the earthquake that started the tsunami and a landslide that amplified the tsunami. Both Chezy's coefficient and the slope factor are found.

The second case study is a landslide down the Morsárjökull outlet glacier in Iceland, starting as a rock slide, but it turns downhill into a glacier and turns into a flow slide, in clear contrast to the Tohoku submarine slide. This slide flows first down a steep upper slope and then down a milder lower slope and then finally stops on an even milder slope with a flow factor zero. The results for the upper and lower sections show different slope factors, the lower one belonging to the lower section. The lower thickness value belongs to the upper section, but the Chezy coefficients

are similar for the two sections. There is reason to believe that further research can uncover a pattern for the Chezy coefficient for similar sites. Then runout lengths and risk area can be approximately predicted for slides of known volume, known from developing cracks or other discernable geological evidence.

**References**

[1] Highland, LM, Bobrowsky P. The Landslide Handbook—A Guide to Understanding Landslides, United States Geological Survey, Landslide Program and National Landslide Information Center, Mail Stop 966, Box 25046, Denver Federal Center, Denver, Colorado, 80225 USA. 2008. Available from: https://pubs.usgs.gov/circ/1325/

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

> assessment using a physically based model and digital elevation data. Environmental Geology. 1998;**35**(2-3):

[9] Tsai TL, Yang JC. Modeling of rainfall-triggered shallow landslide. Environmental Geology. 2006;**50**(4):

[10] Farahmand A, AghaKouchak A. A satellite-based global landslide model. Natural Hazards and Earth System Sciences. 2013;**13**(5):1259-1267

[11] Elíasson J, Kjaran SP, Holm SL, Gudmundsson MT, Larsen G. Large hazardous floods as translatory waves. Environmental Modelling and Software.

2007 Oct 1;**22**(10):1392-1399

[12] Douglas JF, Gasiorek JM,

**41**(6):1084-1097

IntechOpen; 2019

England: Longman Scientific &

Swaffield JA. Fluid Mechanics. 3rd Edn.

Technical, Essex CM20 2JE; 1995. p. 819

[13] McDougall S, Hungr O. A model for the analysis of rapid landslide motion across three-dimensional terrain. Canadian Geotechnical Journal. 2004;

[14] Eliasson J. Earthquake-generated

[15] JAMSTEC. Causal Mechanisms of

Earthquake of 2011. Revealed through Hydraulic Analysis of Fault Drilling Samples from the Deep–Sea Scientific Drilling Vessel Chikyu; Japan Agency

landslides and tsunamis. In: Earthquakes-Impact, Community Vulnerability and Resilience. Rijeka:

Large Slip during the Tohoku

for Marine-Earth Science and Technology. 2013. Available from: http://www.jamstec.go.jp/e/about/

[16] Land and Sea Areas of Crustal Movement and Slip Distribution Model of the Tohoku-Pacific Ocean

press\_release/20131008/

81-88

525-534

[2] van Westen C. Introduction to Landslides, Part 1: Types and Causes. Available from: http://www.adpc.net/ casita/Course%20Modules/Landslide% 20hazard%20assessment/Landslides\_

[3] Komac M. A landslide susceptibility model using the analytical hierarchy process method and multivariate statistics in perialpine Slovenia. Geomorphology. 2006;**74**(1-4):17-28

Turcotte DL, Reichenbach P. Power-law correlations of landslide areas in central Italy. Earth and Planetary Science Letters. 2002;**195**(3-4):169-183

[5] Casadei M, Dietrich WE, Miller NL. Testing a model for predicting the timing and location of shallow landslide initiation in soil-mantled landscapes. Earth Surface Processes and Landforms:

[6] Crozier MJ. Prediction of rainfalltriggered landslides: A test of the antecedent water status model. Earth Surface Processes and Landforms. 1999;

[7] Hergarten S, Neugebauer HJ. Selforganized criticality in a landslide model. Geophysical Research Letters.

[8] Borga M, Dalla Fontana G, Da Ros D, Marchi L. Shallow landslide hazard

types\_and\_causes.pdf

[4] Guzzetti F, Malamud BD,

The Journal of the British Geomorphological Research Group. 2003;**28**(9):925-950

**24**(9):825-833

**159**

1998;**25**(6):801-804

The third case study is a landslide down the Móafellshyrna Mountain in North Iceland. It has previously been shown that this slide can be modeled with the RAMS model. Here it is demonstrated that it can also be modeled using the Chezy coefficients estimated in case study 2. The possibility to keep a data set containing Chezy coefficients that can be used in future estimates is a considerable help in hazard assessments.

## **Author details**

Jónas Elíasson\* and Þorsteinn Sæmundsson Faculty of Engineering and Natural Sciences, University of Iceland, Reykjavík, Iceland

\*Address all correspondence to: jonase@hi.is

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

*The Translatory Wave Model for Landslides DOI: http://dx.doi.org/10.5772/intechopen.90264*

## **References**

are similar for the two sections. There is reason to believe that further research can uncover a pattern for the Chezy coefficient for similar sites. Then runout lengths and risk area can be approximately predicted for slides of known volume, known

The third case study is a landslide down the Móafellshyrna Mountain in North Iceland. It has previously been shown that this slide can be modeled with the RAMS model. Here it is demonstrated that it can also be modeled using the Chezy coefficients estimated in case study 2. The possibility to keep a data set containing Chezy coefficients that can be used in future estimates is a considerable help in hazard

from developing cracks or other discernable geological evidence.

*Landslides - Investigation and Monitoring*

assessments.

**Author details**

Iceland

**158**

Jónas Elíasson\* and Þorsteinn Sæmundsson

\*Address all correspondence to: jonase@hi.is

provided the original work is properly cited.

Faculty of Engineering and Natural Sciences, University of Iceland, Reykjavík,

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

[1] Highland, LM, Bobrowsky P. The Landslide Handbook—A Guide to Understanding Landslides, United States Geological Survey, Landslide Program and National Landslide Information Center, Mail Stop 966, Box 25046, Denver Federal Center, Denver, Colorado, 80225 USA. 2008. Available from: https://pubs.usgs.gov/circ/1325/

[2] van Westen C. Introduction to Landslides, Part 1: Types and Causes. Available from: http://www.adpc.net/ casita/Course%20Modules/Landslide% 20hazard%20assessment/Landslides\_ types\_and\_causes.pdf

[3] Komac M. A landslide susceptibility model using the analytical hierarchy process method and multivariate statistics in perialpine Slovenia. Geomorphology. 2006;**74**(1-4):17-28

[4] Guzzetti F, Malamud BD, Turcotte DL, Reichenbach P. Power-law correlations of landslide areas in central Italy. Earth and Planetary Science Letters. 2002;**195**(3-4):169-183

[5] Casadei M, Dietrich WE, Miller NL. Testing a model for predicting the timing and location of shallow landslide initiation in soil-mantled landscapes. Earth Surface Processes and Landforms: The Journal of the British Geomorphological Research Group. 2003;**28**(9):925-950

[6] Crozier MJ. Prediction of rainfalltriggered landslides: A test of the antecedent water status model. Earth Surface Processes and Landforms. 1999; **24**(9):825-833

[7] Hergarten S, Neugebauer HJ. Selforganized criticality in a landslide model. Geophysical Research Letters. 1998;**25**(6):801-804

[8] Borga M, Dalla Fontana G, Da Ros D, Marchi L. Shallow landslide hazard

assessment using a physically based model and digital elevation data. Environmental Geology. 1998;**35**(2-3): 81-88

[9] Tsai TL, Yang JC. Modeling of rainfall-triggered shallow landslide. Environmental Geology. 2006;**50**(4): 525-534

[10] Farahmand A, AghaKouchak A. A satellite-based global landslide model. Natural Hazards and Earth System Sciences. 2013;**13**(5):1259-1267

[11] Elíasson J, Kjaran SP, Holm SL, Gudmundsson MT, Larsen G. Large hazardous floods as translatory waves. Environmental Modelling and Software. 2007 Oct 1;**22**(10):1392-1399

[12] Douglas JF, Gasiorek JM, Swaffield JA. Fluid Mechanics. 3rd Edn. England: Longman Scientific & Technical, Essex CM20 2JE; 1995. p. 819

[13] McDougall S, Hungr O. A model for the analysis of rapid landslide motion across three-dimensional terrain. Canadian Geotechnical Journal. 2004; **41**(6):1084-1097

[14] Eliasson J. Earthquake-generated landslides and tsunamis. In: Earthquakes-Impact, Community Vulnerability and Resilience. Rijeka: IntechOpen; 2019

[15] JAMSTEC. Causal Mechanisms of Large Slip during the Tohoku Earthquake of 2011. Revealed through Hydraulic Analysis of Fault Drilling Samples from the Deep–Sea Scientific Drilling Vessel Chikyu; Japan Agency for Marine-Earth Science and Technology. 2013. Available from: http://www.jamstec.go.jp/e/about/ press\_release/20131008/

[16] Land and Sea Areas of Crustal Movement and Slip Distribution Model of the Tohoku-Pacific Ocean Earthquake: Coseismic Slip Distribution Model. Report by the Geographical Survey Institute, Japan. 2011. Available from: http://www.gsi.go.jp/cais/topic 110520-index-e.html

[17] Fujii Y, Satake K, Sakai SI, Shinohara M, Kanazawa T. Tsunami Earthquakes - Impact, Community Vulnerability and Resilience source of the 2011 off the Pacific coast of Tohoku earthquake. Earth, Planets and Space. 2011;**63**(7):55

[18] Sæmundsson Þ, Sigurðsson IA, Pétursson HG, Jónsson HP, Decaulne A, Roberts MJ, et al. Bergflóðið sem féll á Morsárjökul 20. mars 2007 –hverjar hafa afleiðingar þess orðið? Náttúrufræðingurinn. 2011;**81**(3–4): 131-141 (in Icelandic)

[19] Post AS. Alaskan Glaciers: recent observations in respect to the earthquake-advance theory. Science. 1965;**148**:366-368

[20] Shreve RL. Sherman landslide, Alaska. Science. 1966;**154**(3757): 1639-1643

[21] Fahnestock RK. Little Tahoma Peak rockfalls and avalanches, Mount Rainier, Washington, USA. In: Voight B, editor. Developments in Geotechnical Engineering. Vol. 14. Part A. New York: Elsevier; 1978. pp. 181-196

[22] Sæmundsson Þ, Morino C, Helgason JK, Conway SJ, Pétursson HG, Barceló D, et al. The triggering factors of the Móafellshyrna debris slide in northern Iceland: Intense precipitation, earthquake activity and thawing of mountain permafrost. In: Barceló D, Gan J, editors. Science of the Total Environment. 2018;**621**:1163-1175

[23] Morino C, Conway SJ, Sæmundsson Þ, Helgason JK, Hiller J, Butcher FEG, et al. Molards as an indicator of permafrost degradation and landslide processes. Earth and Planetary Science Letters. 2019;**516**:136-147

**161**

**Chapter 8**

*Herman Peiffer*

**Abstract**

rehabilitation

**1. Introduction**

stress at rest K0 [1].

viscoplastic models.

depends on the actual stress conditions.

Techniques to Evaluate and

Overconsolidated Clay

Remediate the Slope Stability in

In this chapter, the origin and remediation of an important sliding in the overconsolidated Boom Clay in Kruibeke (Belgium) is discussed. Local and

be used to evaluate the stability of the slope after remediation.

**Keywords:** slope stability, monitoring, in situ testing, DMT (dilatometer),

A natural or an artificial slope in a clay deposit results in a change of the morphological properties and the stress state in the soil This evolution of the stress state can be measured by in situ tests considering the variations over time of the lateral

K0 is an important parameter in the design and monitoring in many geotechnical applications (retaining walls, building foundations, and definitely the stability of slopes). For the evaluation of the stability of slopes, it can be estimated out of the results of field tests (Pressuremeter Tests [1], Vane Tests and Piezometer Cone Tests [2], Dilatometer Test [3–5] or numerical modeling, usually considering elasto-

Although the analysis of the lateral stress, in case of slope stability, also can be made based on the results of laboratory tests [6–8], one of the main problems is that K0, due to the natural variations of horizontal and vertical stresses, at each location

environmental factors caused an unstable slope about 30 m deep, A larger sensitivity to erosion resulted finally in the instability of the slope. Because of the formation of fine cracks in the soil there was a possibility for the water to penetrate in the clay close to the surface, resulting in the presence of higher water pressures. Also, the presence of the excavator on top of the slope during exploitation had an important impact on the stress state of the soil. Both an analytical and numerical approach were used to estimate the factor of safety. Because of the change of the soil characteristics, the factor of safety decreases, which can be estimated through a numerical analysis (using the Strength Reduction Method). This chapter also discussed the applied techniques for the remediation using numerical analysis. Also, the importance of the field test is discussed. An integrated approach, using numerical analysis and field tests in combination, is capable of predicting the instability. This approach can also

## **Chapter 8**

Earthquake: Coseismic Slip Distribution Model. Report by the Geographical Survey Institute, Japan. 2011. Available from: http://www.gsi.go.jp/cais/topic

*Landslides - Investigation and Monitoring*

110520-index-e.html

2011;**63**(7):55

[17] Fujii Y, Satake K, Sakai SI, Shinohara M, Kanazawa T. Tsunami Earthquakes - Impact, Community Vulnerability and Resilience source of the 2011 off the Pacific coast of Tohoku earthquake. Earth, Planets and Space.

[18] Sæmundsson Þ, Sigurðsson IA, Pétursson HG, Jónsson HP, Decaulne A, Roberts MJ, et al. Bergflóðið sem féll á Morsárjökul 20. mars 2007 –hverjar

Náttúrufræðingurinn. 2011;**81**(3–4):

[19] Post AS. Alaskan Glaciers: recent

[20] Shreve RL. Sherman landslide, Alaska. Science. 1966;**154**(3757):

[21] Fahnestock RK. Little Tahoma Peak rockfalls and avalanches, Mount Rainier, Washington, USA. In: Voight B, editor. Developments in Geotechnical Engineering. Vol. 14. Part A. New York:

Elsevier; 1978. pp. 181-196

Letters. 2019;**516**:136-147

**160**

[22] Sæmundsson Þ, Morino C,

the Móafellshyrna debris slide in northern Iceland: Intense precipitation, earthquake activity and thawing of mountain permafrost. In: Barceló D, Gan J, editors. Science of the Total Environment. 2018;**621**:1163-1175

Helgason JK, Conway SJ, Pétursson HG, Barceló D, et al. The triggering factors of

[23] Morino C, Conway SJ, Sæmundsson Þ, Helgason JK, Hiller J, Butcher FEG, et al. Molards as an indicator of permafrost degradation and landslide processes. Earth and Planetary Science

hafa afleiðingar þess orðið?

observations in respect to the earthquake-advance theory. Science.

131-141 (in Icelandic)

1965;**148**:366-368

1639-1643
