**Stability and Run-out Analysis of Earthquake-induced Landslides**

Yingbin Zhang

[32] Wieland, M., 2007, The seismic performance of concrete face rockfill dams under strong ground Shaking*.* International Water Power and Dam Construction, April,

202 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

[33] Wieland, M., 2013, Seismic design of major components. International Water Power

[34] Wieland, M., 2014, Seismic hazard and seismic design and safety aspects of large dam projects. Second European Conference on Earthquake Engineering and Seismol‐

18-20.

and Dam Construction, February, 16-19.

ogy, Istanbul, Aug. 25-29, 2014.

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/59439

## **1. Introduction**

A large number of landslides can be caused by a strong earthquake and they have been the source of significant damage and loss of people and property. Therefore, it is very important to predict the stability of slope and the movement behaviors of a potential landslide under an earthquake loading, i.e., stability and run-out analysis (Figure 1).

Earthquake-induced landslides have been the source of significant damage and loss of people and property. One of the most serious event is the 1970 Peru earthquake. This event caused a huge rock avalanche that killed almost 54,000 people and buried two cities [143]. Another example is, in the 1920 Haiyuan earthquake, a large number of landslides caused widespread damage to infrastructure and buildings and killed at least 100,000 people, almost half of the total earthquake deaths [82].

Therefore, it is very important to predict the earthquake-induced landslides and to take countermeasures for potential landslides.

Main topics of earthquake-induced landslides are the following:


management of earthquake-related landslides, and (iii) hazard map and early warning system for earthquake-related landslides

**Figure 1.** Stability and run-out analysis of earthquake-induced landslides

This chapter focuses on the prediction of potential earthquake-induced landslides. The prediction of potential landslide can be carried out using detailed geotechnical investigations and stability calculations. (i) Failure mechanism and stability analysis of seismic slopes, i.e. seismic slope stability analysis and (ii) movement mechanism and behaviours of earthquakeinduced landslides, i.e. landside run-out analysis are outlined firstly, and then the merits and demerits of each method are clarified in this chapter.

## **2. Seismic slope stability analysis**

So far, methods developed to analyze the stability of earthquake slopes can be divided into three types: (1) pseudo-static methods, (2) dynamic sliding block methods, and (3) stress-strain methods. These three types of methods can be applied in different cases due to each of them has merit and demerit [73].

#### **2.1. Pseudo-static methods**

[166] first presented the pseudo-static method, which is a simple method for evaluating of seismic stability of a slope. This type of method can be used to man-made or natural slopes based on either analytical method or numerical method. The earthquake force, acting on the an element or whole of the slope, is writed by a horizontal force and/or a vertical volum force equal to the gravitation force multiple a coefficient *k*, called the pseudo-static coefficient as shown in Figure 2 and Equation (1).

**Figure 2.** Forces acting on a slope in pseudo-static slope stability analysis

management of earthquake-related landslides, and (iii) hazard map and early warning

204 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

This chapter focuses on the prediction of potential earthquake-induced landslides. The prediction of potential landslide can be carried out using detailed geotechnical investigations and stability calculations. (i) Failure mechanism and stability analysis of seismic slopes, i.e. seismic slope stability analysis and (ii) movement mechanism and behaviours of earthquakeinduced landslides, i.e. landside run-out analysis are outlined firstly, and then the merits and

So far, methods developed to analyze the stability of earthquake slopes can be divided into three types: (1) pseudo-static methods, (2) dynamic sliding block methods, and (3) stress-strain methods. These three types of methods can be applied in different cases due to each of them

[166] first presented the pseudo-static method, which is a simple method for evaluating of seismic stability of a slope. This type of method can be used to man-made or natural slopes

system for earthquake-related landslides

**Figure 1.** Stability and run-out analysis of earthquake-induced landslides

demerits of each method are clarified in this chapter.

**2. Seismic slope stability analysis**

has merit and demerit [73].

**2.1. Pseudo-static methods**

Thus, *k* times the gravitational acceleration *g*, i.e. *a=kg* forms the assumed seismic acceleration *a*. The assumed pseudo-static forces acting on a potential sliding mass of weight *W* will be

$$\begin{aligned} f\_h &= \frac{a\_h W}{\mathcal{g}} = k\_h W \\ f\_v &= \frac{a\_v W}{\mathcal{g}} = k\_v W \end{aligned} \tag{1}$$

where *ah* and *av* are horizontal and vertical pseudo-static accelerations, respectively, *kh* and *kv* are horizontal and vertical pseudo-static coefficients, respectively. The factor of safety(FOS) is represented as the ratio of the resisting force to the driving force, Equation (2).

$$\text{FOS} = \frac{\tau\_r}{\tau\_d} \tag{2}$$

From Equation (1), the pseudo-static force is determined by the seismic coefficient. The key problem for the pseudo-static procedure is how to select an appropriate seismic coefficient under an acceptable FOS. There have been studies for determining the most appropriate pseudo-static coefficient by a matter of experience and judgment.

[166] classical paper made the original suggestion to use of *kh*=0.1 for severe earthquakes, *kh*=0.2 for violent and/or destructive earthquakes, and of *kh*=0.5 for catastrophic earthquakes.

[103] presented a minimum pseudo-static FOS of 1.5 based on a slope material strength reduction factor (SRF) of 0.8 and the following acceleration values associated with two different earthquake magnitudes *M*. The same values of seismic coefficients for magnitude 6.5 and 8.25 earthquakes are recommended by [154], but with an acceptable FOS of 1.15.

$$\begin{aligned} a &= 0.1 \text{g for} \, M = 6.5 \text{ implying} \, k = 0.1\\ a &= 0.15 \, \text{g for} \, M = 8.25 \text{ implying} \, k = 0.15 \end{aligned} \tag{3}$$

[137] also presented the pseudo-static coefficient related to earthquake magnitude. In detail, for an 8.25, 7.5, 7.0 and 6.5 magnitude earthquakes, if the seismic coefficients equal to 1/2, 1/3, 1/4 and 1/5 of the PGA, respectively, the computed FOSs are larger than 1.0, the accumulated displacements of slope are likely to be acceptably small.

In the report published by the International Commission of Large Dams (ICOLD), [154] shows a list of the minimum FOS value and horizontal seismic coefficients for 14 large dams world‐ wide, in which the minimum FOSs range from 1.0 to 1.5 and the horizontal earthquake coefficients range from 0.1 to 0.15. The Corps of Engineers Manual recommended a earthquake coefficient of 0.1 or 0.15 for areas where major and great earthquake threats are estimated, respectively, and a FOS of no larger than 1.0 for all magnitude earthquakes.

Some references related the earthquake coefficient value to the peak ground acceleration (PGA) [10, 67, 108]. [108] related a pseudo-static coefficient of 1/3 to 1/2 of the PGA at the top of a double-side slope (a dam in the source reference), whereas [67] related a pseudo-static coefficient of 1/2 of the PGA of bedrock (PGArock) with a FOS of no larger than 1.0 and a SRF of 20%. And, [10] recommended the pseudo-static coefficient of 0.6 or 0.75 times of the PGA of bedrock (0.6 or 0.75PGArock). It should be noted that the value given by [10] is conservative because the original study is designed for solid-waste landfills, where the allowable deforma‐ tion are relatively small. [89] pointed that although engineering judgment is required for all cases, the criteria of [67] should be appropriate for most slopes.

[91] suggested one-half of PGA to use in an area of low seismicity (peak acceleration <0.15g) for the stability of earth embankments. This can be obtained from the peak horizontal motion (mean) from Modified Mercalli Intensity (MMI), magnitude-distance attenuation and the probability of a 50-year, 90% nonexceedance. However, in an area of moderate to strong seismicity (0.15*g*≤PGA≤0.40*g*), PGA is obtained from the peak horizontal motion, from MMI, magnitude-distance attenuation and probability of 250-year, 90% nonexceedance.

[76] suggested a minimum FOS of 1.0, also based on a slope material SRF of 0.8 and the following values of pseudo-static coefficient: *a* equals to 0.17PGA or 0.5PGA for the dynamic response analysis is to be performed for the slope or earthquake structure or not.

[163] developed an expression for the earthquake coefficient in terms of characters of ground motion and magnitude of earthquake based on the data of [10].

It is almost common that only the horizontal acceleration is considered in evaluating the stability and deformation of a slope because the horizontal acceleration is the principal destabilizing force that acts on earth structures as well as the principal source of damage observed in earthquakes [4].

From Figure 2, the horizontal force clearly increases the driving force and decreases the FOS. The vertical pseudo-static force generally has less influence on the FOS than the horizontal pseudo-static force does because the vertical pseudo-static reduces both the driving and resisting forces. Hence, the effects of vertical seismic loading are frequently omitted in pseudostatic analysis [89].

Several investigators performed some analyses and have shown that the inclination of seismic loading have a significant influence on the seismic stability of slope by coupling the vertical and horizontal components of seismic force [20, 100].

In summary, pseudo-static method can be simply and directly used to identify the FOS and the critical seismic coefficient *kc*. In addition, performance of slope is closely related to permanent displacement, but the results of pseudo-static method are difficult to interpret the performance of slope after a seismic event because this method provides no information about permanent displacement. Because the pseudo-static analysis method provides only a rough assessment of seismic slope, it should be only used for the preliminary procedures. More accurate methods can be used to the followed process [73, 163, 170].

#### **2.2. Dynamic sliding block methods**

[166] classical paper made the original suggestion to use of *kh*=0.1 for severe earthquakes, *kh*=0.2 for violent and/or destructive earthquakes, and of *kh*=0.5 for catastrophic earthquakes.

[103] presented a minimum pseudo-static FOS of 1.5 based on a slope material strength reduction factor (SRF) of 0.8 and the following acceleration values associated with two different earthquake magnitudes *M*. The same values of seismic coefficients for magnitude 6.5 and 8.25

[137] also presented the pseudo-static coefficient related to earthquake magnitude. In detail, for an 8.25, 7.5, 7.0 and 6.5 magnitude earthquakes, if the seismic coefficients equal to 1/2, 1/3, 1/4 and 1/5 of the PGA, respectively, the computed FOSs are larger than 1.0, the accumulated

In the report published by the International Commission of Large Dams (ICOLD), [154] shows a list of the minimum FOS value and horizontal seismic coefficients for 14 large dams world‐ wide, in which the minimum FOSs range from 1.0 to 1.5 and the horizontal earthquake coefficients range from 0.1 to 0.15. The Corps of Engineers Manual recommended a earthquake coefficient of 0.1 or 0.15 for areas where major and great earthquake threats are estimated,

Some references related the earthquake coefficient value to the peak ground acceleration (PGA) [10, 67, 108]. [108] related a pseudo-static coefficient of 1/3 to 1/2 of the PGA at the top of a double-side slope (a dam in the source reference), whereas [67] related a pseudo-static coefficient of 1/2 of the PGA of bedrock (PGArock) with a FOS of no larger than 1.0 and a SRF of 20%. And, [10] recommended the pseudo-static coefficient of 0.6 or 0.75 times of the PGA of bedrock (0.6 or 0.75PGArock). It should be noted that the value given by [10] is conservative because the original study is designed for solid-waste landfills, where the allowable deforma‐ tion are relatively small. [89] pointed that although engineering judgment is required for all

[91] suggested one-half of PGA to use in an area of low seismicity (peak acceleration <0.15g) for the stability of earth embankments. This can be obtained from the peak horizontal motion (mean) from Modified Mercalli Intensity (MMI), magnitude-distance attenuation and the probability of a 50-year, 90% nonexceedance. However, in an area of moderate to strong seismicity (0.15*g*≤PGA≤0.40*g*), PGA is obtained from the peak horizontal motion, from MMI,

[76] suggested a minimum FOS of 1.0, also based on a slope material SRF of 0.8 and the following values of pseudo-static coefficient: *a* equals to 0.17PGA or 0.5PGA for the dynamic

magnitude-distance attenuation and probability of 250-year, 90% nonexceedance.

response analysis is to be performed for the slope or earthquake structure or not.

*a gM <sup>k</sup>* (3)

earthquakes are recommended by [154], but with an acceptable FOS of 1.15.

206 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

== = == = 0.1 for 6.5 implying 0.1 0.15 for 8.25 implying 0.15

*a gM k*

respectively, and a FOS of no larger than 1.0 for all magnitude earthquakes.

cases, the criteria of [67] should be appropriate for most slopes.

displacements of slope are likely to be acceptably small.

Displacement-based dynamic slidng block method is another alternative approach to evaluate the seismic slope stability, as permanent displacement is a useful index of slope performance, especially for those man-made slopes constructed for special purposes such as dams, em‐ bankments et al. This method has been widely used in earthquake geotechnical engineering.

In 1965, [119] proposed the dynamic sliding block method for estmating the permanent displacement of embankment affected by a seismic loading. In this method, sliding would be induced once the seismic loading exceed the critical seismic force of a potential failure surface as shown in Figure 3. The sliding would be accumulated until the end of seismic loading. We can evaluate the accumulated permanent displacement to assess the seismic stability of a slope.

Newmark's method showes that the yield acceleration of a potential block is a function of the FOS and slope angle, as:

$$a\_c = (\text{FOS} - \text{l}) \text{g} \sin \alpha \tag{4}$$

where *ac* is in terms of the gravity acceleration *g*; FOS is the static factor of safety; and *α* is the slope angle.

Since then, the method has been numerous extensions and applications. The section 2.2.1 and section 2.2.2 will give reviews for these two aspects, respectively. In addition, a regional scale application of the dynamic sliding block method is reviewed in section 2.2.3.

**Figure 3.** Illustration of the original Newmark's method

#### *2.2.1. Extensions*

More attention has been focused over the last decades on developing methods to more accurately analyze the seismic stability of a slope for dams, embankments or other important structures by modeling the dynamic response of a slope more rigorously.

After the first dynamic sliding rigid block method, [155 and 97] published more sophisticated methods to account for the un-rigid block. Similar studies also given by [103]. As the classifi‐ cation given by [73], methods for estimating the permanent displacement of a sliding system induced by earthquake loading can been grouped into: (1) rigid-block model [119], (2) decoupled model [10, 104], and (3) coupled model [11, 97, 139].

**Figure 4.** Applications of dynamic sliding block method in geotechnical engineering

#### *2.2.2. Applications*

where *ac* is in terms of the gravity acceleration *g*; FOS is the static factor of safety; and *α* is the

Since then, the method has been numerous extensions and applications. The section 2.2.1 and section 2.2.2 will give reviews for these two aspects, respectively. In addition, a regional scale

More attention has been focused over the last decades on developing methods to more accurately analyze the seismic stability of a slope for dams, embankments or other important

After the first dynamic sliding rigid block method, [155 and 97] published more sophisticated methods to account for the un-rigid block. Similar studies also given by [103]. As the classifi‐ cation given by [73], methods for estimating the permanent displacement of a sliding system induced by earthquake loading can been grouped into: (1) rigid-block model [119], (2)

structures by modeling the dynamic response of a slope more rigorously.

decoupled model [10, 104], and (3) coupled model [11, 97, 139].

application of the dynamic sliding block method is reviewed in section 2.2.3.

208 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

**Figure 3.** Illustration of the original Newmark's method

*2.2.1. Extensions*

slope angle.

Since the rigid-block method was published in 1965 by Newmark, it has seen numerous applications, four of which are shown in Figure 4. The applications in recent years include (1) the seismic deformation analysis of earth dams and embankments [1, 22, 48, 49, 89, 90, 97, 103, 138, 144, 145, 150, 155, 179, 180]; (2) the displacements associated with landslides [34, 53, 70, 171]; (3) the seismic deformation of landfills with geosynthetic liners [10, 181]; (4) the seismic settlement of surface foundations [141]; and (5) the potential sliding of concrete gravity dams [32, 47, 95]. The extension of the analogue by [140] to gravity retaining walls has met worldwide acceptance, and has found its way into seismic codes of practice. Several other generalised applications have also appeared (e.g. [2, 3, 45, 99, 139, 162, 169]).

#### *2.2.3. Regional scale analysis*

Except a single slope analysis, where the landslides are likely to occur and what kind of seismic conditions will cause it failure are two important topics in seismic hazard assessment, i.e.regional scale analysis [59].

For a regional scale analysis, slope stability analysis menthds will be not suitable [143, 168]. With the development of Geographic Information Systems (GIS) tools in recent years, regional scale analyses by the dynamic sliding block method have been proposed, in which ground shaking characteristic parameters, geotechnical material and topographic data are considered (e.g. [34, 71, 75, 106, 114, 151, 155]).

The Newmark analysis (which combines slope stability calculations with seismic groundmotion records) is widely used to evaluate the potential for landslides that could be triggered by earthquake shaking [70, 71, 72, 74, 113].

#### **2.3. Stress-strain methods**

With the developments of the simulation approach and computer technology in recent years, the stress-strain method is becoming increasingly used in seismic slope stability analysis. These methods can be grouped into continuous methods, e.g. finite element method (FEM) [21], finite difference method (FDM) [116], boundary element method (BEM) [12], and discontinuous methods, e.g. rigid block spring method (RBSM) [77; 80], discontinuous deformation analysis (DDA) [159, 160] and discrete element method (DEM) [31].

#### *2.3.1. Continuous methods*

[21] developed and named FEM of engineering analysis, in which the studied system is meshed into small many elements. This method can be applied to estimate the slope stability including dynamic stability analysis.

Some applications of the continuous methods have been proposed, e.g. [89, 94, 153] and [156]. Recently, nonlinear in-elastic soil models have been developed and implemented in twodimension (2-D) and three-dimension (3-D) models (e.g., [42, 50, 135, 164]). In addition, [93] and [183] studied the seismic slope stability by using FDM.

#### *2.3.2. Discontinuous methods*

For the analysis of a potential failure mass consisting of multiple blocks as shown in Figure 5, the discontinuous methods are more applicable [120]. Some applications of RBSM and DEM can be found in some literature (e.g. [8, 52, 77, 79, 80, 85121, 127128, 129, 130, 131, 136, 182]).

**Figure 5.** A jointed rock slope (modified from Bhasin and Kaynia, 2004)

DDA is also a discontinuous method developed for the modeling of the behaviors of multiple block systems. Since the novel formulation and the numerical code of DDA were presented, DDA draws more and more attention and many extensions and modifications to the original method have been proposed to overcome some limitations [19, 37, 38, 81, 87, 98] and make it more suitable, practical and efficient to seismic slope stability.

The DDA can be used both to static rock slope engineering (e.g. [17, 81, 102, 123, 176, 187]) and the seismic rock slope stability analysis [56, 57, 54].

In summary, stress-strain method represents a powerful alternative approach for seismic slope stability analysis which is accurate, versatile and requires fewer a priori assumptions, espe‐ cially, regarding the shape of failure surface.

## **3. Landslide run-out analysis**

It is important to estimate the movement behaviour of a potential landslide. For example, the movement distance is an important parameter in risk assessment and measure design. There are many run-out analysis methods, which can fall into four categories: (1) experimental methods, (2) empirical methods, (3) analytical methods, and (4) numerical simulation methods. The states of the art of these methods are reviewed in the following four subsections 3.1 - 3.4.

### **3.1. Experiment methods**

scale analyses by the dynamic sliding block method have been proposed, in which ground shaking characteristic parameters, geotechnical material and topographic data are considered

210 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

The Newmark analysis (which combines slope stability calculations with seismic groundmotion records) is widely used to evaluate the potential for landslides that could be triggered

With the developments of the simulation approach and computer technology in recent years, the stress-strain method is becoming increasingly used in seismic slope stability analysis. These methods can be grouped into continuous methods, e.g. finite element method (FEM) [21], finite difference method (FDM) [116], boundary element method (BEM) [12], and discontinuous methods, e.g. rigid block spring method (RBSM) [77; 80], discontinuous deformation analysis

[21] developed and named FEM of engineering analysis, in which the studied system is meshed into small many elements. This method can be applied to estimate the slope stability including

Some applications of the continuous methods have been proposed, e.g. [89, 94, 153] and [156]. Recently, nonlinear in-elastic soil models have been developed and implemented in twodimension (2-D) and three-dimension (3-D) models (e.g., [42, 50, 135, 164]). In addition, [93]

For the analysis of a potential failure mass consisting of multiple blocks as shown in Figure 5, the discontinuous methods are more applicable [120]. Some applications of RBSM and DEM can be found in some literature (e.g. [8, 52, 77, 79, 80, 85121, 127128, 129, 130, 131, 136, 182]).

(e.g. [34, 71, 75, 106, 114, 151, 155]).

**2.3. Stress-strain methods**

*2.3.1. Continuous methods*

dynamic stability analysis.

*2.3.2. Discontinuous methods*

by earthquake shaking [70, 71, 72, 74, 113].

(DDA) [159, 160] and discrete element method (DEM) [31].

and [183] studied the seismic slope stability by using FDM.

**Figure 5.** A jointed rock slope (modified from Bhasin and Kaynia, 2004)

Physical modelling typically involves using scale models to capture the motion of landslides. Physical experiments are usually preferred to models because models require more assump‐ tions than direct measurements. But for landslides, direct experiment is difficult, dangerous, expensive, and of limited utility. Based on laboratory experiments and filed investigation data, there are many different available models developed for calculating run-out zones.

Some full-scale direct experiments with artificial landslides have been completed [118, 122124, 125, 126] and others). However, since landslides are frequently heterogeneous and single event cannot be repeated carefully through adjusting only one factor, direct experiment is difficult, dangerous, expensive and of limited utility. And observing conditions are complicated by the danger of being in close proximity to a landslide and the difficulty of measuring a material with properties that change when observed in-situ or when isolated for measurement. But laboratory experiments are still the first qualitative and quantitative observations on the obtained results became fundamental for a better understanding of movement runout behaviour.

#### **3.2. Empirical methods**

Several empirical methods for assessing landslide travel distance and velocity have been developed based on historical data and on the analysis of the relationship between parameters characterizing both the landslide, e.g. the volume of the landslide mass, and the path, e.g. local morphology, and the distance travelled by the failure mass [65]. Regression model-based methods and geomorphology-based methods are two kinds of common methods.

#### *3.2.1. Regression model-based methods*

The regression model-based methods are developed on an apparent inverse relationship between landslide volume and angle of reach (also called as fahrböschund by [58]). Several linear regression equations have been proposed [25, 96, 153]. Introduced by [58], the angle of reach is the inclination of the line connecting the crest of the source with the toe of the deposit, as measured along the approximate streamline of motion. The angle of reach is considered an index of the efficiency of energy dissipation, and so is inversely related to mobility. Similar correlations between volume and other simple mobility indices have been proposed [33, 60, 142]. Given estimated source location, volume and path direction, these methods provide estimates of the distal limit of motion [111].

Improved empirical model notable performing regressions on subsets with varying scopes were presented by [13, 25, 69] and others.

Regression model-based models play a valuable role in landslide run-out analysis due to the regression model-based methods are simple. But the regression model-based methods are difficult to apply in practice with a high degree of certainty. For example, the correlation coefficients for some of regression models are 0.7-0.8, while a value of larger than 0.95 generally indicates a strong correlation. And it is difficult in this method to take account of influences of the ground condition, the micro-topography, the degree of saturation of the landslide mass and et al. For this point, geomorphology-based method is another alternative approach to predict the run-out of landslide.

#### *3.2.2. Geomorphology-based methods*

Field work and photo interpretation are the main sources of the geomorphological analysis for determining the travel distance of landslides [65]. The outer margin of the landslide deposits give an appraisal of the maximum distances that landslides have been able to reach during the present landscape (Figure 6). Several authors have provided these studies (e.g. [23, 24, 26, 88]).

The geomorphological approach does not give any clue of the emplacement mechanism. Furthermore, the slope geometry and the circumstances responsible for past landslides might have changed. Therefore, results obtained in a given place cannot be easily exported to other localities.

In summary, empirical methods, both regression model-based methods and geomorphologybased methods, typically predict travel distances, while the deformation characteristics or the slide velocities of the landslide are not predicted. These models may be applied to establish initial hazard characteristics for preliminary run-out analysis, which may be later refined by other models.

**Figure 6.** Boundary of the potential rockfall run-out area in Santa Coloma (Principality of Andorra), defined by the line that links the farthest fallen blocks observed in the field [24]. Arrows indicate historical rockfall paths and solid circles are large fallen boulders [65]

#### **3.3. Analytical methods**

morphology, and the distance travelled by the failure mass [65]. Regression model-based

The regression model-based methods are developed on an apparent inverse relationship between landslide volume and angle of reach (also called as fahrböschund by [58]). Several linear regression equations have been proposed [25, 96, 153]. Introduced by [58], the angle of reach is the inclination of the line connecting the crest of the source with the toe of the deposit, as measured along the approximate streamline of motion. The angle of reach is considered an index of the efficiency of energy dissipation, and so is inversely related to mobility. Similar correlations between volume and other simple mobility indices have been proposed [33, 60, 142]. Given estimated source location, volume and path direction, these methods provide

Improved empirical model notable performing regressions on subsets with varying scopes

Regression model-based models play a valuable role in landslide run-out analysis due to the regression model-based methods are simple. But the regression model-based methods are difficult to apply in practice with a high degree of certainty. For example, the correlation coefficients for some of regression models are 0.7-0.8, while a value of larger than 0.95 generally indicates a strong correlation. And it is difficult in this method to take account of influences of the ground condition, the micro-topography, the degree of saturation of the landslide mass and et al. For this point, geomorphology-based method is another alternative approach to

Field work and photo interpretation are the main sources of the geomorphological analysis for determining the travel distance of landslides [65]. The outer margin of the landslide deposits give an appraisal of the maximum distances that landslides have been able to reach during the present landscape (Figure 6). Several authors have provided these studies (e.g. [23,

The geomorphological approach does not give any clue of the emplacement mechanism. Furthermore, the slope geometry and the circumstances responsible for past landslides might have changed. Therefore, results obtained in a given place cannot be easily exported to other

In summary, empirical methods, both regression model-based methods and geomorphologybased methods, typically predict travel distances, while the deformation characteristics or the slide velocities of the landslide are not predicted. These models may be applied to establish initial hazard characteristics for preliminary run-out analysis, which may be later refined by

methods and geomorphology-based methods are two kinds of common methods.

212 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

*3.2.1. Regression model-based methods*

estimates of the distal limit of motion [111].

were presented by [13, 25, 69] and others.

predict the run-out of landslide.

*3.2.2. Geomorphology-based methods*

24, 26, 88]).

localities.

other models.

In contrast to empirical methods, analytical methods are based on mechanics and involve the solution of motion equations [111]. The simplest analytical model is the classical sliding block model as shown in Figure 7, which is based on work-energy theory [6, 9, 43, Müller-Bernet in 58, 63, 83, 84, 132, 147]. Internal deformation and it's associated energy dissipation are neglected and the landslides is treated as a lumped mass. At any position along the path, the sum of the energies including the potential energy, kinetic energy and net energy loss equals the initial potential energy. This energy balance can be visualized using the concept of energy grade lines, as shown in Figure 9. The concept of energy grade lines is useful for visualizing the energy balance. *v* is the velocity of the block, *g* is the vertical acceleration due to gravity and *v*<sup>2</sup> */*2*g* is known as the velocity head, which is the kinetic energy of the block normalized by the product of its mass and *g*. The same normalization of net energy loss is known as head loss. Note that the positions of the energy lines are referenced to the centre of mass of the block and that the true energy line and mean energy line do not necessarily coincide. Given the initial position of the center of mass and a suitable relationship to approximate the energy losses, the position and velocity of the block can be determined at any given time.

Three-dimensional analysis for investigating runout of a slope were also proposed [36, 40, 51, 92 and 109]. These models require a high resolution Digital Elevation Model (DEM).

Generally speaking, the use of analytical methods is somewhat motivated by the limitations of purely empirical methods, as the unique geometry and materials involved in each case can be accounted for explicitly and a statistically-significant database of previous events is not

**Figure 7.** The classical sliding block model, based on work-energy theory [111]

necessarily required. The simplicity of a lumped mass allows analytical solutions, fast and effectively [66]. However, because the landslide is reduced to a single point, lumped mass models cannot provide the exact maximum runout distance, but only the displacement concerning the centre of mass [44, 62].

#### **3.4. Numerical simulation methods**

The single-block model should be only applied to the motion of the center of mass of a rigid body, but more complex continuum deformable mass or multi-block system is often appeared in practice. Some numerical simulation methods have been developed to account explicitly for deformation during motion.

#### *3.4.1. Continuous methods*

When considering that the dimensions of a typical particle is much smaller than the depth and length of the debris, the debris mass is treated as continuum. According to depth averaged Saint Venant approach, the material is assumed to be incompressible and the mass and momentum equations are written in a depth-averaged form. Many numerical methods now exist to investigate the run-out process of landslide (e.g. [18, 27, 28, 29, 30, 35, 62, 111, 112, 133, 149, 161, 165]). These methods are usually based on continuum mechanics and assume that the avalanche thickness is very much smaller than its extent parallel to the bed, i.e. thin layer depth-averaged models. The primary differences are their representation of basal resistance force and the constitutive relations describing the mechanical behaviour of the considered material. These models can accurately take account of detailed topography effects, shown to be significant, with a reasonable computational time, making it possible to perform sensitivity studies of the parameters used in the model. They can provide effective properties that make it possible to roughly reproduce not only the deposit shape but also the dynamic as shown in [46) and 117] for examples. However, conventional continuum approaching models, which neglects the contact between rocks, makes it impossible to trace the position of individual rock during a landslide.

#### *3.4.2. Discontinuous methods*

When the landslide mass consists of large fragments and boulders, the run-out mass is modelled as an assembly of blocks moving down a surface. Some authors take circular shape models in their run-out analysis to evaluate maximum runout and final deposit position of past or potential events (e.g., [134]). Although polygonal shapes have the disadvantages due to the complexity of the contact patterns and penalty in computational time, methods using non-circular shapes will be required for more real-world problems. It is more appropriate when problems are limited in finite blocks. Discontinuous numerical simulation methods are powerful tools in simulation of failure and run-out process of rock avalanche controlled by weakness surface. DEM [31] and DDA [159, 157] are two of the most commonly used methods.

Both DEM and DDA employ the equations of dynamic motion which are solved at finite points in time, in a series of time steps, but there are some subtle but significant differences in their formulations of the solution schemes and contact mechanics. In the solution schemes, equa‐ tions of motion in DDA are derived using the principle of minimization of the total potential energy of the system, while the equations of motion as implemented in DEM are derived directly from the force balance equations, which still resultant unbalanced force after a time step and damping is necessarily used to dissipate energy. In the contact mechanics, the DDA used a penalty method in which the contact is assumed to be rigid. No overlapping or interpenetration of the blocks is allowed as the same as real physical cases, whereas soft contact approach is used in DEM. The soft contact approach requires laboratory or field measured joint stiffness, which may be difficult to obtain in many cases. Many comparisons of basic models (sliding, colliding and rolling models) between the DEM and DDA were carried out and show that the results from DDA are more close to the analytical values than that from DEM [188]. Compared to DEM, DDA has a simpler and more straightforward physical meaning [172].

Applications of DEM can be found in some literatures, such as [85, 128, 129, 131, 136, 182].

DDA can be used for estimating the affected area of an earthquake-induced landslide.

[55] first validated the applicability of DDA for the dynamic behaviour of block sliding on an slope. Based on the same inputs model of seismic loadings, [7576886105, 158, 167, 178] studied the dynamic response or/and stability analysis of tunnel, slope, dam, foundation or ancient masonry structure by using DDA. Alternatively, the seismic loadings also can be applied to the base block [146, 147], which is different from the original DDA. Later, [173, 174, 175, 177] applied DDA to simulate the kinematic behavior of sliding rock blocks in the Tsaoling landslide and the Chiu-fen-erh-shan landslide induced by the 1999 Chi-Chi earthquake. Recently, [184] applied newest DDA program to simulate the largest landslides induced by the 2008 Wen‐ chuan earthquake.

## **4. Comparisons of various methods**

necessarily required. The simplicity of a lumped mass allows analytical solutions, fast and effectively [66]. However, because the landslide is reduced to a single point, lumped mass models cannot provide the exact maximum runout distance, but only the displacement

214 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

The single-block model should be only applied to the motion of the center of mass of a rigid body, but more complex continuum deformable mass or multi-block system is often appeared in practice. Some numerical simulation methods have been developed to account explicitly for

When considering that the dimensions of a typical particle is much smaller than the depth and length of the debris, the debris mass is treated as continuum. According to depth averaged Saint Venant approach, the material is assumed to be incompressible and the mass and momentum equations are written in a depth-averaged form. Many numerical methods now exist to investigate the run-out process of landslide (e.g. [18, 27, 28, 29, 30, 35, 62, 111, 112, 133, 149, 161, 165]). These methods are usually based on continuum mechanics and assume that the avalanche thickness is very much smaller than its extent parallel to the bed, i.e. thin layer depth-averaged models. The primary differences are their representation of basal resistance force and the constitutive relations describing the mechanical behaviour of the considered material. These models can accurately take account of detailed topography effects, shown to be significant, with a reasonable computational time, making it possible to perform sensitivity studies of the parameters used in the model. They can provide effective properties that make it possible to roughly reproduce not only the deposit shape but also the dynamic as shown in [46) and 117] for examples. However, conventional continuum approaching models, which neglects the contact between rocks, makes it impossible to trace the position of individual rock

When the landslide mass consists of large fragments and boulders, the run-out mass is modelled as an assembly of blocks moving down a surface. Some authors take circular shape

concerning the centre of mass [44, 62].

**Figure 7.** The classical sliding block model, based on work-energy theory [111]

**3.4. Numerical simulation methods**

deformation during motion.

*3.4.1. Continuous methods*

during a landslide.

*3.4.2. Discontinuous methods*

The studies in the field of the earthquake-induced landslides are generally reviewed. Two parts of contents, i) seismic stability analysis and ii) run-out analysis are reviewed and compared. Some conclusions can be drawn:


## **5. A case study: The Daguangbao landslide [185, 186]**

#### **5.1. Background**

The Daguangbao landslide is located in the hanging wall only 6.5 km away from the Yingxiu-Beichuan fault. It is a typical bedding landslide.

Figure 8 gives a pre- and post-earthquake 3-D topographies, from which cross-section of the Daguangbao landslide can be obtained (Figure 9). The extent of the damage caused by the Daguangbao landslide is reflected in the following statistics [61]:


**Figure 8.** Pre- and post-earthquake 3-D topographies of the Daguangbao landslide

**Figure 9.** Cross-section of the Daguangbao landslide before and after the Wenchuan earthquake.

#### **5.2. Material properties and ground motion**

#### *5.2.1. Material properties*

**1.** Three categories methods can be used to analyse the seismic stability of a slope. Each of these types of methods has strengths and weaknesses and each can be appropriately applied in different situations. In detail, pseudo-static methods can simply and directly determine the FOS and the critical coefficient *kc* of a slope, while the widely used New‐ mark's methods and its extensions can determine the co-seismic deformation of a slope. And the Newmark's methods can be used to estimate where earthquake-induced landslides are likely to occur and what kind of shaking conditions will trigger them based on the GIS technology. More sophisticated analysis for real dynamic process of a seismic slope should be carried out by stress-strain methods, including both continuous methods

216 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

**2.** Four kinds of methods can be used to analyse the run-out of a landslide. In detail, experiment method can provide the qualitative and quantitative observations on the obtained results although this method is difficult, dangerous, expensive and of limited utility. Empirical method can be directly used for assessing landslide travel distance and velocity based on historical data and on the analysis of the relationship between param‐ eters characterizing both landslide and the path. Analytical method can be more directly used without the need of statistically-significant database of previous events. Numerical simulation method can be used to provide more information for the landslide composed

The Daguangbao landslide is located in the hanging wall only 6.5 km away from the Yingxiu-

Figure 8 gives a pre- and post-earthquake 3-D topographies, from which cross-section of the Daguangbao landslide can be obtained (Figure 9). The extent of the damage caused by the

;

;

by the complex continuum deformable mass or multi-blocks.

**5. A case study: The Daguangbao landslide [185, 186]**

Daguangbao landslide is reflected in the following statistics [61]:

**3.** Estimated volume of collapsed rock mass is 750~840 million m3

Beichuan fault. It is a typical bedding landslide.

**1.** The affected area covered 7.3~10 km2

**4.** The failure zone is more than 1 km;

**5.** The failure mass moved about 4.5 km;

**6.** Formed a 600m high landslide dam.

**2.** The accumulation body width is 2.2 km;

and discontinuous methods.

**5.1. Background**

The Daguangbao landslide is so huge that the size effect must be considered. To account for this discrepancy, experience equations based on Hoek-Brown failure criterion, which size effect can be considered, is used to back calculate the material strength. Table 1 lists the material properties of the Daguangbao landslide.


**Table 1.** Material properties of the Daguangbao landslide in FLAC3D and DDA

#### *5.2.2. Ground motion*

The horizontal earthquake wave is the projection combination in the main sliding direction (N60ºE) use the MZQP acceleration records in E-W and N-S directions as Equation (5). The inputted vertical earthquake wave is the MZQP acceleration records in U-D direction.

$$a\_H = a\_{E-W} \cdot \sin 60^\circ + a\_{N-S} \cdot \cos 60^\circ \tag{5}$$

Figure 10 shows the input combined acceleration records. Velocity and displacement time histories can be obtained by first and second integration from acceleration record. The duration of earthquake wave is 60s.

**Figure 10.** Input post-corrected horizontal and vertical ground records projected to N60ºE direction

### **5.3. Numerical simulations-Run-out analysis**

Seismic DDA can successfully simulate the movement of earthquake induced landslide. Two main features determine the Daguangbao landslide is a unique case, one is near-fault location (≈6.5 km) and the other one is huge scale (≈800×10<sup>6</sup> m3 ). The near-fault location determines the Daguangbao landslide must be shocked by the extreme ground motion from the strong Wenchuan earthquake. And the Daguangbao landslide located on the meizoseismal area where the vertical seismic component is very large. In addition, the landslide is so huge that the size effects must be considered. The friction coefficient measured in the laboratory is no longer suitable for stability and run-out analysis.

To these two features, the Daguangbao landslide is simulated by the newest seismic DDA code in which multi-direction seismic forces can be applied in the base block directly, and experience equations based on Hoek-Brown failure criterion is applied to back-calculate the material strength by trying to consider the size effect.


#### **Table 2.** Control parameters for DDA

*5.2.2. Ground motion*

Tensile strength of discontinuities (*σ* <sup>t</sup>

**Table 1.** Material properties of the Daguangbao landslide in FLAC3D and DDA

of earthquake wave is 60s.

Acceleration[g]

Displacement [m]

Velocity [m/sec]

The horizontal earthquake wave is the projection combination in the main sliding direction (N60ºE) use the MZQP acceleration records in E-W and N-S directions as Equation (5). The inputted vertical earthquake wave is the MZQP acceleration records in U-D direction.

) : kPa 12 32 556

Density (*ρ*): g/cm3 2.5 2.6 2.6 260000 Unit weight of rock (*γ*) : kN/m3 25 26 26 0 Elastic modulus (*E*) : Gpa 1.86 2.63 14.76 Poisson's ratio (*ν*) 0.2 0.2 0.1 Friction angle of discontinuities (*φ*) : º 10.8 12.18 23.53 Cohesion of discontinuities (*c*) : Mpa 1.276 1.576 4.052

218 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

Figure 10 shows the input combined acceleration records. Velocity and displacement time histories can be obtained by first and second integration from acceleration record. The duration

Time [sec]

**Figure 10.** Input post-corrected horizontal and vertical ground records projected to N60ºE direction

Time [sec]

Time [sec]

sin 60 cos60 *H EW N S aa a* = × °+ × ° - - (5)

**Material 1 Material 2**

(a)

**Materia 3**

**FLAC3D DDA**

(b)

(c)

#### *5.3.1. Geometry of sliding blocks*

The main sliding direction of the Daguangbao landslide, N60ºE, is selected as analysis profile. The DDA model is depicted in Figure 11. In this simulation, based on the shape of failure surface and the character of slope topography, the whole slope is divided into three parts: base block, upper sliding mass, and lower sliding mass. Then two sliding masses are divided into the smaller discrete deformable blocks based on pre-existing discontinuities.

**Figure 11.** DDA model of the Daguangbao landslide

#### *5.3.2. Results*

Figure 15 shows the post-failure behavior of the Daguangbao landslide simulated by the seismic DDA code. Simulated results show that the sliding blocks climb over the Pingliangzi. After overlapping the final step of DDA calculation with the topographic cross-section at the Daguangbao landslide, the deposit pattern of the simulated Daguangbao landslide under horizontal-and-vertical situation coincides well with local topography.

**Figure 12.** Simulation results of the Daguangbao landslide

### **6. Conclusions**

Five cases are performed using finite difference program FLAC3D, under the real seismic waves near the study site. The results show that the seismic conditions cause a significant reduction in factor of safety than static situation. It also found that the vertical seismic has a significant influence on tension failure of block, although it has an insignificant influence on change of the factor of safety. Another important conclusion is the effect of vertical seismic force on relative displacement of potential sliding mass is significant. In addition, large area of tension failure caused by the combined seismic forces at back edge of the slope applies the evidence of effect of vertical seismic force on failure mechanism of slope stability.

A comparison of simulation results from three situations, static, only-horizontal and horizon‐ tal-and-vertical, is carried out. Seismic force has a significant influence on the arrival distance, and shape of post-failure. Arrival distance from horizontal-and-vertical situation is larger than that from only-horizontal situation. In addition, the deposit pattern of the simulated Dag‐ uangbao landslide under horizontal-and-vertical situation coincides well with local topogra‐ phy. The vertical seismic force should be considered for landslide assessment and management, especially in the situation that the studied site located on the meizoseismal area during the earthquake.

## **Acknowledgements**

*5.3.2. Results*

Figure 15 shows the post-failure behavior of the Daguangbao landslide simulated by the seismic DDA code. Simulated results show that the sliding blocks climb over the Pingliangzi. After overlapping the final step of DDA calculation with the topographic cross-section at the Daguangbao landslide, the deposit pattern of the simulated Daguangbao landslide under

220 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

Five cases are performed using finite difference program FLAC3D, under the real seismic waves near the study site. The results show that the seismic conditions cause a significant reduction in factor of safety than static situation. It also found that the vertical seismic has a significant influence on tension failure of block, although it has an insignificant influence on change of the factor of safety. Another important conclusion is the effect of vertical seismic force on relative displacement of potential sliding mass is significant. In addition, large area of tension

horizontal-and-vertical situation coincides well with local topography.

**Figure 12.** Simulation results of the Daguangbao landslide

**6. Conclusions**

The authors gratefully acknowledge the financial support provided by the National Natural Science Foundation of China (No. 51408511), the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry and the opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection (Chengdu Univer‐ sity of Technology) (No. SKLGP2014K015).

## **Author details**

#### Yingbin Zhang\*

Address all correspondence to: yingbinz516@126.com

Department of Geotechnical Engineering, School of Civil Engineering, Southwest Jiaotong University, Chengdu, China

## **References**


[20] Chopra, A.K. 1966. The importance of the vertical component of earthquake motions. Bulletin of the Seismological Society of America, 56(5): 1163-1175.

[4] Anderson, D., and Kavazanjian Jr, E. 1995. Performance of landfills under seismic loading. In Proc., 3rd Int. Conf. on Recent Advances in Geotechnical Earthquake En‐

[6] Azzoni, A., La Barbera, G., and Zaninetti, A. 1995. Analysis and prediction of rock‐ falls using a mathematical model. In International Journal of Rock Mechanics and

[7] Bakun‐Mazor, D., Hatzor, Y., and Glaser, S. 2012. Dynamic sliding of tetrahedral wedge: The role of interface friction. International Journal for Numerical and Analyt‐

[8] Bhasin R, Kaynia AM. 2004. Static and dynamic simulation of a 700-m high rock

[9] Bozzolo, D., and Pamini, R. 1986. Simulation of rock falls down a valley side. Acta

[10] Bray, J.D., and Rathje, E.M. 1998. Earthquake-induced displacements of solid-waste landfills. Journal of Geotechnical and Geoenvironmental Engineering, 124(3):

[11] Bray, J.D., and Travasarou, T. 2007. Simplified procedure for estimating earthquakeinduced deviatoric slope displacements. Journal of Geotechnical and Geoenviron‐

[12] Brebbia, C.A., and Wrobel, L. 1980. The boundary element method. Computer meth‐

[13] Cannon, S.H. 1993. An empirical model for the volume-change behavior of debris

[14] CDCDMG:California Department of Conversation, Division of Mines and Geology, 1997. Guidelines for evaluating and mitigating seismic hazards in California. CDMG

[15] Chang, K.-T., Lin, M.-L., Dong, J.-J., and Chien, C.-H. 2011. The Hungtsaiping land‐

[16] Chen, L., and Zhao, W. 1979. Longling earthquake, 1976. Earthquake press, Beijing.

[17] Chen, G., and Ohnishi, Y. 1999. Slope stability analysis using Discontinuous Defor‐

[18] Chen, H., and Lee, C. 2000. Numerical simulation of debris flows. Canadian Geotech‐

[19] Cheng, Y.M. 1998. Advancements and improvement in discontinuous deformation

ods in fluids.(A 81-28303 11-34) London, Pentech Press, Ltd., 1980: 26-48.

gineering and Soil Dynamics. Univ. of Missouri Rolla, MO, Vol.3, pp. 277-306.

[5] Army, U.C.o.E. 1960. Stability of Earth and Rockfill Dams, EM 1110-2-1902.

222 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

Mining Sciences & Geomechanics Abstracts. Elsevier, Vol.32, pp. 709-724.

slope in western Norway. Engineering Geology, 71(3-4): 213-226.

ical Methods in Geomechanics, 36(3): 327-343.

Mechanica, 63(1-4): 113-130.

mental Engineering, 133(4): 381-392.

flows. In Hydraulic Engineering '93, San Francisco.

slides: from ancient to recent. Landslides, 9(2): 205-214.

analysis. Computers and Geotechnics, 22(2): 153-163.

mation Analysis method. Rock Mecganics for Industry: 535-541.

242-253.

Special Publication.

nical Journal, 37(1): 146-160.


[48] Franklin, A.G., and Chang, F.K. 1977. Permanent displacements of earth embank‐ ments by Newmark sliding block analysis.

[34] Del Gaudio, V., Pierri, P., and Wasowski, J. 2003. An approach to time-probabilistic evaluation of seismically induced landslide hazard. Bulletin of the Seismological So‐

224 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

[35] Denlinger, R.P., and Iverson, R.M. 2001. Flow of variably fluidized granular masses across three‐dimensional terrain: 2. Numerical predictions and experimental tests.

[36] Descoeudres, F., and Zimmermann, T. 1987. Three-dimensional dynamic calculation

[37] Doolin, D.M. 2005. Unified displacement boundary constraint formulation for dis‐ continuous deformation analysis (DDA). International Journal for Numerical and

[38] Doolin, D.M., and Sitar, N. 2004. Time Integration in Discontinuous Deformation

[39] Dorren, L.K. 2003. A review of rockfall mechanics and modelling approaches. Prog‐

[40] Dorren, L., and Heuvelink, G.B. 2004. Effect of support size on the accuracy of a dis‐ tributed rockfall model. International Journal of Geographical Information Science,

[41] Dreyfus, D., Rathje, E.M., and Jibson, R.W. 2013. The Influence of Different Simpli‐ fied Sliding-Block Models and Input Parameters on Regional Predictions of Seismic

[42] Elgamal, A.-W.M., Scott, R.F., Succarieh, M.F., and Yan, L. 1990. La Villita dam re‐ sponse during five earthquakes including permanent deformation. Journal of Geo‐

[43] Evans, S., and Hungr, O. 1993. The assessment of rockfall hazard at the base of talus

[44] Evans, S., Hungr, O., and Enegren, E. 1994. The Avalanche Lake rock avalanche, Mackenzie mountains, northwest territories, Canada: description, dating, and dy‐

[45] Fardis, M.N. 2009. Seismic design, assessment and retrofitting of concrete buildings:

[46] Favreau, P., Mangeney, A., Lucas, A., Crosta, G., and Bouchut, F. 2010. Numerical

[47] Fenves, G., and Chopra, A.K. 1986. Simplified analysis for earthquake resistant de‐ sign of concrete gravity dams. Earthquake Engineering Research Center, University

modeling of landquakes. Geophysical Research Letters, 37(15): L15305.

Landslides Triggered by the Northridge Earthquake. Engineering Geology.

Journal of Geophysical Research: Solid Earth (1978–2012), 106(B1): 553-566.

ciety of America, 93(2): 557-569.

of rockfalls. In 6th ISRM Congress.

ress in Physical Geography, 27(1): 69-87.

technical Engineering, 116(10): 1443-1462.

based on EN-Eurocode 8. Springer.

slopes. Canadian Geotechnical Journal, 30(4): 620-636.

namics. Canadian Geotechnical Journal, 31(5): 749-768.

18(6): 595-609.

of California.

Analytical Methods in Geomechanics, 29(12): 1199-1207.

Analysis. Journal of Engineering Mechanics, 130(3): 249-258.


[77] Kawai, T. 1977. A new discrete analysis of nonlinear solid mechanics problems in‐ volving stability, plasticity and crack. In the Symposium on Applications of Comput‐ er Methods in Engineering, Los Angeles, USA, pp. 1029-1038.

[63] Hungr, O., and Evans, S. 1988. Engineering evaluation of fragmental rockfall haz‐ ards. In Proceedings of the Fifth International Symposium on Landslides, Lausanne,

226 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

[64] Hungr, O., and Evans, S. 2004. Entrainment of debris in rock avalanches: An analysis of a long run-out mechanism. Geological Society of America Bulletin, 116(9-10):

[65] Hungr, O., Corominas, J., and Eberhardt, E. 2004. Estimating landslide motion mech‐

[66] Hürlimann, M., Rickenmann, D., Medina, V., and Bateman, A. 2008. Evaluation of approaches to calculate debris-flow parameters for hazard assessment. Engineering

[67] Hynes-Griffin, M.E., and Franklin, A.G. 1984. Rationalizing the seismic coefficient

[68] Ishikawa, T., Sekine, E., and Ohnishi, Y. 2002. Shaking table tests of coarse granular materials with discontinuous analysis. Proc. of ICADD-5, BALKEMA: 181-187. [69] Jakob, M., and Hungr, O. 2005. Debris-flow hazards and related phenomena. Spring‐

[70] Jibson, R.W. 1993. Predicting earthquake-induced landslide displacements using

[71] Jibson, R. 2000. A method for producing digital probabilistic seismic landslide haz‐

[72] Jibson, R.W. 2007. Regression models for estimating coseismic landslide displace‐

[73] Jibson, R.W. 2011. Methods for assessing the stability of slopes during earthquakes—

[74] Jibson, R.W., and Jibson, M.W. 2003. Java programs for using Newmark's method and simplified decoupled analysis to model slope performance during earthquakes.

[75] Jibson, R.W., Harp, E.L., and Michael, J.A. 1998. A Method for Producing Digital Probabilistic Seismic Landslide Hazard Maps: An Example from the Los Angeles,

[76] Kavazanjian, E., and Consultants, G. 1997. Design Guidance: Geotechnical Earth‐ quake Engineering for Highways. Design Principles. Federal Highway Administra‐

Newmark's sliding block analysis. Transportation Research Record: 9-9.

AA Balkema, Rotterdam, Netherlands, pp. 685-690.

method. Defense Technical Information Center.

anism, travel distance and velocity.

ard maps. Engineering Geology.

California, Area.

tion.

ment. Engineering Geology, 91(2-4): 209-218.

A retrospective. Engineering Geology, 122(1-2): 43-50.

US Department of the Interior, US Geological Survey.

Geology, 102(3): 152-163.

1240-1252.

er.


[105] Makris, N., and Roussos, Y. 2000. Rocking response of rigid blocks under near-source ground motions. Geotechnique, 50(3): 243-262.

[91] Krinitzsky, E.L. 1993. Fundamentals of earthquake-resistant construction. Wiley.

[92] Lan, H., Derek Martin, C., and Lim, C. 2007. RockFall analyst: A GIS extension for three-dimensional and spatially distributed rockfall hazard modeling. Computers &

[93] Latha, G.M., and Garaga, A. 2010. Seismic Stability Analysis of a Himalayan Rock

[94] Lee, K.L. 1974. Seismic permanent deformation in earth dams, University of Califor‐

[95] Leger, P., and Katsouli, M. 1989. Seismic stability of concrete gravity dams. Earth‐

[96] Li, T. 1983. A mathematical model for predicting the extent of a major rockfall. Zeits‐

[97] Lin, J.S., and Whitman, R.V. 1983. Decoupling approximation to the evaluation of earthquake‐induced plastic slip in earth dams. Earthquake Engineering & Structural

[98] Lin, C.T., Amadei, B., Jung, J., and Dwyer, J. 1996. Extensions of discontinuous defor‐ mation analysis for jointed rock masses. international Journal of Rock Mechanics and

[99] Ling, H.I. 2001. Recent applications of sliding block theory to geotechnical design.

[100] Ling, H., and Leshchinsky, D. 1998. Effects of vertical acceleration on seismic design

[101] Lomnitz, C. 1970. Casualties and behavior of populations during earthquakes. Bulle‐

[102] Luan, M., Li, Y., and Yang, Q. 2000. Discontinuous deformation computational me‐ chanics model and its application in stability analysis of rock slope. Chinese Journal

[103] Makdisi, F.I., and Seed, H.B. 1977. Simplified procedure for estimating dam and em‐ bankment earthquake-induced deformations. In ASAE Publication No. 4-77. Pro‐ ceedings of the National Symposium on Soil Erosion and Sediment by Water,

[104] Makdisi, F.I., and Seed, H.B. 1978. Simplified procedure for estimating dam and em‐ bankment earthquake-induced failures. Journal of the Geotechnical Division, ASCE,

of geosynthetic-reinforced soil structures. Geotechnique, 48(3): 347-373.

Slope. Rock Mechanics and Rock Engineering, 43(6): 831-843.

228 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

quake Engineering & Structural Dynamics, 18(6): 889-902.

Mining Sciences & Geomechanics Abstracts, 33(7): 671-694.

Soil Dynamics and Earthquake Engineering, 21: 189-197.

tin of the Seismological Society of America, 60(4): 1309-1313.

of Rock Mechanics and Engineering, 3: 006.

Chicago, Illinois, December 12-13, 1977.

104: 849–861.

chrift Fur Geomorphologie, 24: 473-482.

com.

Geosciences, 33(2): 262-279.

nia, Los Angeles, CA.

Dynamics, 11(5): 667-678.


[134] Poisel, R., Preh, A., and Hungr, O. 2008. Run Out of Landslides–Continuum Mechan‐ ics versus Discontinuum Mechanics Models. Geomechanics and Tunnelling, 1(5): 358-366.

[119] Newmark, N.M. 1965. Effects of earthquakes on dams and embankments. Géotechni‐

230 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

[120] Ning, Y., and Zhao, Z. 2012. A detailed investigation of block dynamic sliding by the discontinuous deformation analysis. International Journal for Numerical and Analyt‐

[121] Niwa, K., Kawai, T., Ikeda, M., and Takeda, T. 1984. Application of a new discrete method to fracture analysis of brittle materials. In the 3rd International Conference

[122] Ochiai, H., Okada, Y., Furuya, G., Okura, Y., Matsui, T., Sammori, T., Terajima, T., and Sassa, K. 2004. A fluidized landslide on a natural slope by artificial rainfall.

[123] Ohnishi, Y., Chen, G., and Miki, S. 1995. Recent development of DDA in rock me‐

[124] Okura, Y., Kitahara, H., and Sammori, T. 2000. Fluidization in dry landslides. Engi‐

[125] Okura, Y., Kitahara, H., Sammori, T., and Kawanami, A. 2000. The effects of rockfall

[126] Okura, Y., Kitahara, H., Ochiai, H., Sammori, T., and Kawanami, A. 2002. Landslide fluidization process by flume experiments. Engineering Geology, 66(1): 65-78.

[127] Pal, S., Kaynia, A.M., Bhasin, R.K., and Paul, D.K. 2011. Earthquake Stability Analysis

[128] Papaloizou, L., and Komodromos, P. 2009. Planar investigation of the seismic re‐ sponse of ancient columns and colonnades with epistyles using a custom-made soft‐

[129] Papantonopoulos, C., Psycharis, I.N., Papastamatiou, D.Y., Lemos, J.V., and Mouza‐ kis, H.P. 2002. Numerical prediction of the earthquake response of classical columns using the distinct element method. Earthquake Engineering & Structural Dynamics,

[130] Pekau, O.A, and Yuzhu, C. 2004. Failure analysis of fractured dams during earth‐

[131] Pekau, O.A, and Yuzhu, C. 2004. Seismic collapse behaviour of Damaged dams. In 13 WCEE: 13 th World Conference on Earthquake Engineering Conference Proceedings.

[132] Pfeiffer, T.J., and Bowen, T. 1989. Computer simulation of rockfalls. Bulletin of the

[133] Pirulli, M. 2005. Numerical modelling of landslide runout. A continuum mechanics

quakes by DEM. Engineering Structures, 26(10): 1483-1502.

Association of Engineering Geologists, 26(1): 135-146.

approach, Politecnico di Torino.

of Rock Slopes: a Case Study. Rock Mechanics and Rock Engineering.

ware. Soil Dynamics and Earthquake Engineering, 29(11-12): 1437-1454.

volume on runout distance. Engineering Geology, 58(2): 109-124.

on Numerical Methods in Fracture Mechanics, Swansea, U.K., pp. 13-27.

que, 15: 139-159.

ical Methods in Geomechanics: 1-21.

Landslides, 1(3): 211-219.

31(9): 1699-1717.

chanics. Proc. ICADD, 1: 26-47.

neering Geology, 56(3): 347-360.


[148] Sassa, K. 1988. Motion of Landslides and Debris Flows: Prediction of Hazard Area: Report for Grant-in-aid for Scientific Research by Japanese Ministry on Education, Science and Culture (project No. 61480062). Disaster Prevention Research Institute.

[149] Savage, S., and Hutter, K. 1989. The motion of a finite mass of granular material

[150] Sawada, T., Chen, W.F., and Nomachi, S.G. 1993. Assessment of seismic displace‐ ments of slopes. Soil Dynamics and Earthquake Engineering, 12: 357-362.

[151] Saygili, G., and Rathje, E.M. 2009. Probabilistically based seismic landslide hazard maps: An application in Southern California. Engineering Geology, 109(3): 183-194.

[152] Scheidegger, A.E. 1973. On the prediction of the reach and velocity of catastrophic

[153] Seed, H.B. 1973. Analysis of the Slides in the San Fernando Dams During the Earth‐ quake of Feb. 9, 1971: Report to State of California Department of Water Resources, Los Angeles Department of Water and Power, National Science Foundation. College

[154] Seed, H.B. 1979. Considerations in the earthquake-resistant design of earth and rock‐

[155] Seed, H.B., and Martin, G.R. 1966. The seismic coefficient in earth dam design. Jour‐

[156] Serff, N. 1976. Earthquake induced deformations of earth dams. College of Engineer‐

[157] Shi, G.-H. 1988. Discontinuous Deformation Analysis A New Numerical Model for the Statics and Dynamics of Block Systems, University of California, Berkeley.

[158] Shi, G. 2002. Single and multiple block limit equilibrium of key block method and discontinuous deformation analysis. In Proceedings of the 5th International Confer‐ ence on Analysis of Discontinuous Deformation. Rotterdam: AA Balkema, pp. 3-43.

[159] Shi, G.-H., and Goodman, R.E. 1985. Two dimensional discontinuous deformation analysis. International Journal for Numerical and Analytical Methods in Geomechan‐

[160] Shi, G.-H., and Goodman, R.E. 1989. Generalization of two-dimensional discontinu‐ ous deformation analysis for forward modelling. International Journal for Numerical

[161] Sousa, J., and Voight, B. 1991. Continuum simulation of flow failures. Geotechnique,

[162] Stamatopoulos, C. 1996. Sliding system predicting large permanent co-seismic move‐ ments of slopes. Earthquake Engineering & Structural Dynamics, 25(10): 1075-1093.

and Analyrucal Methods in Geomechanics, 13: 359-380.

nal of Soil Mechanics & Foundations Div, 92(Proc. Paper 4824).

down a rough incline. Journal of Fluid Mechanics, 199(1): 177-215.

232 Earthquake Engineering - From Engineering Seismology to Optimal Seismic Design of Engineering Structures

landslides. Rock Mechanics, 5(4): 231-236.

of Engineering, University of California.

fill dams. Geotechnique, 29(3): 13-41.

ing, University of California.

ics, 9: 541-556.

41(4): 515-538.


**Chapter 9**
