21 - Slope Stability and Landslides
21 - Slope Stability and Landslides
Landslides go by many different names from rockfalls to mudlows, but the point
is that they all slide and they all involve shearing resistance. In terms of lives and
property lost, landslides are among the most tragic and devastating events that
involve soils. Landslides occur naturally as streams cut valleys deeper and wider,
and as lakes and oceans erode their shores. Landslides are triggered by shaking
from an earthquake, or by willful indiscretions of man.
The most common naturally occurring slopes are stream valley walls and eroding
shorelines, but slopes also occur as exposed sides of faults and uplifted mountain
ranges. Depositional slopes include sand dunes, deltas, and volcanic cinder cones.
The most common depositional slopes are man-made—road embankments, earth
dams, piles of sand, grain, or industrial or mining waste. Other man-made slopes
include cuts made by trenching and by excavations for foundations, basements, or
roadcuts. All slopes have one thing in common: if conditions are right, they have
it in them to host a landslide.
This chapter will discuss the soil mechanics of naturally occurring slopes and
landslides, followed by design considerations for cut slopes and embankments to
prevent slope failures.
The best time to stop a landslide is before it starts, because shearing and
remolding of soil in the slip zone causes it to lose a substantial part of its shearing
strength. Therefore, after sliding starts it may speed up, depending on geometric
factors and how much strength has been lost. Sliding then continues until the
shearing resistance improves or the ground reaches a configuration that com-
pensates for the loss of soil strength. Sensitive soils move farther and much more
rapidly that those that are less sensitive, sometimes becoming devastating mud
flows that move so fast that they are virtually inescapable.
Normally for a landslide to stop it must move to a flatter surface. This often takes
place gradually, as the sliding soil moves along a concave surface that is nearly
vertical at the top and nearly horizontal at the bottom. Passive resistance of soil
building up at the bottom stops the slide. If that soil is removed, either naturally
or as a result of man’s activities, the slide continues.
One of the worst things that can be done to encourage a landslide to continue also
seems to be the most intuitive, which is to put the soil back where it belongs.
Curiously it sometimes takes two or even three attempts before owners are willing
to admit defeat. Putting the soil back not only removes passive resistance, it puts
weight where it will be most effective at renewing sliding.
Figure 21.1
Causes of this landslide were excavation at the toe of the slope and doing too much laundry at
the top (waste water went into a septic drainfield). The slide enlarged uphill in steps as sliding
removed lateral support from soil exposed in the slide scarp. This landslide was stabilized with
drilled lime.
Active landslides are more readily analyzed because the slip surface can be
located by inspection and by drilling. The analysis can determine the causes, or
if more than one causal factor is involved, prorate responsibility to the several
contributors. Analysis also is required to compare the effectiveness of different
repair methods. Homeowner insurance policies almost inevitably have a waiver
in fine print that excludes damages from ground movements, placing them in
the category with nuclear wars and volcanic eruptions.
The forces involved in an angle of repose are shown in Fig. 21.2, where W is the
weight of an element. As shown in the figure, the weight may be resolved into N,
Figure 21.2
A simple case of
borderline slope
stability: sand at
the angle of
repose.
a force normal to the slip surface, and Sa, a force acting parallel with the surface
and tending to cause sliding. The acting force is
Sa ¼ W sin i ð21:1Þ
N ¼ W cos i ð21:2Þ
From the Coulomb equation, if cohesion is zero, the resisting force is only
a function of N:
Sr ¼ N tan
ð21:3Þ
¼ W cos i tan
tan
FS ¼ ð21:5Þ
tan i
If the factor of safety is 1.0, then i ¼ (eq. 21.4).
Example 21.1
A highway slope is to be cut through a stable sand dune with a factor of safety of 1.2
against slipping. The sand is predominately quartz having an angle of internal friction of
258. What is the design slope angle assuming no influence from water?
Answer: tan i ¼ tan 258/1.2 ¼ 0.389; i ¼ 218.
Equation (21.1) can be written for a unit volume of soil oriented in the direction
of seepage (Fig. 21.3), so the soil unit weight can be substituted for W. Because
the soil is submerged, the submerged unit weight is used to calculate the
acting force. Then
Soil: Sa ¼ ð w Þsin i ð21:6Þ
Figure 21.3
Seepage force S
acting within a
unit volume of soil
having weight W.
where Sa is parallel with the slope. With seepage parallel to the slope the hydraulic
gradient is
iW ¼ hW ¼ sin i ð21:7Þ
Seepage force S ¼ w sin i. Adding this to the soil acting force gives
Sa þ S ¼ sin i ð21:8Þ
The submerged unit weight also is used for the calculation of normal force:
N ¼ ð w Þcos i ð21:9Þ
The factor of safety is
Sr N tan ð w Þcos i tan
FS ¼ ¼ ¼
Sa sin i sin i
ð w Þ tan
FS ¼ ð21:10Þ
tan i
If the factor of safety is 1.0,
ð w Þ
tan i ¼ tan ð21:11Þ
It will be noted that the acting force, in the denominator, depends on the total unit
weight of the soil, whereas the resisting force in the numerator depends on
( w), the submerged unit weight.
Example 21.2
Assume a worst moisture condition for the highway cut in the previous example and
determine the corresponding slope angle. Is this value appropriate for design?
Answer: 1.2 ¼ 0.5 tan 258/tan i; solving for i gives 118. This would be conservative because
fully saturated conditions are not likely to occur in a sand dune.
Example 21.3
Estimate the creep angle for a saturated cohesive soil with ¼ 208.
Answer: tan i ¼ 0.5 tan 208; i ¼ 108, or approximately 1 vertical to 5.5 horizontal.
Figure 21.4
Sliding of unsaturated soil on an inclined weaker layer, where W is the weight of the sliding soil.
This was the case for a tragic failure at Vaiont Dam in northern Italy, where the slip surface was an
old fault plane.
Thus, the steeper the angle, the higher the acting force. Friction depends on the
weight component acting normal to the slip surface:
N ¼ W cos ð21:14Þ
In this case the steeper the angle, the lower the normal force and the lower the
frictional resistance to sliding. The Coulomb equation for maximum resisting
force, which is the same as eq. (18.6) but on a total force basis, is
Sr ¼ cL þ N tan ð21:15Þ
where c is the soil cohesion, L is the length of the slip surface, and is the friction
angle. Hence
Sr cL þ W cos tan
FS ¼ ¼ ð21:17Þ
Sr W sin
If Sr is fully mobilized and equals Sa, the factor of safety is 1.0 and the soil mass
will slide. In this form eq. (21.17) does not include pore pressure, which normally
plays a critical role in landslides, so the equation is applicable only for
unsaturated conditions.
Example 21.4
The size of the triangular mass of soil in Fig. 21.4 is scaled from the drawing: L ¼ 35 ft
(10.7 m) and ¼ 248. The height of the triangle normal to L is h ¼ 7 ft (2.13 m). Strength
parameters for the inclined weak layer are determined to be approximately c ¼ 50 lb/ft2
(2.4 kPa) and ¼ 208. The unit weight is ¼ 100 lb/ft3 (15.7 kN/m3). What is the factor of
safety against sliding?
Question: Which is the more important component of the resisting force, friction
or cohesion?
The notation in Fig. 21.5 is the same as in Fig. 21.4, and eqs. (21.13) to (21.17)
still apply, except that is not known. It therefore is a simple matter to apply
eq. (21.17) to a series of trial values of , plot the factor of safety versus the
trial value, and pick off the most critical value. An analytical solution is more
precise:
Substituting the value of W from eq. (21.18) into eq. (21.19) and solving for cd
gives
An additional substitution is
Figure 21.5
Culmann’s
method to define
the most critical
slip angle for
planar failure in
homogeneous
soil.
Setting the derivative of cd relative to equal to zero establishes the most critical
value of , which is
c ¼ ½ ð i þ d Þ ð21:21Þ
That is, the most critical failure angle for a planar failure is the average of the slope
angle and the angle of internal friction. The steeper the slope, the higher the failure
angle. This applies only to a planar failure surface.
Example 21.5
According to eq. (21.21), what is the most critical failure angle for a vertical slope? How
does this compare with the angle for Rankine active failure?
Answer: c ¼ ð90 þ d Þ ¼ 45 þ d =2, which is identical to the Rankine active case that
assumes parallel planar failures.
Figure 21.6
Three versions of
critical height for a
vertical cut. When
the cliff lets go,
better toss the
umbrella and get
out of there.
One method for predicting the maximum height of a vertical cut can be obtained
by substituting eq. (21.21) into eq. (21.18) and setting the factor of safety equal to
1.0, but it must be emphasized that this is not an appropriate model, as discussed
later in this section. This substitution gives
4c sin i cos
Hcul ¼
½1 cos ði Þ
where Hcul is the Culmann critical height. Setting i ¼ 908 gives
4c cos
Hcul ¼
1 sin
ð21:22mathematically correct but unsafe!Þ
4c
Hcul ¼ tan 45
2
Equation (21.22) is unsafe for a vertical cut because it does not address the most
likely failure mechanism.
Terzaghi’s Correction
As shown in Fig. 21.6, Terzaghi (1943) indicated that the Culmann Hcul is unsafe
because vertical tension cracks running parallel to the faces of an open vertical cut
reduce the length and the amount of restraint along the slip plane. He therefore
reduced the critical height by one-third. While that is a substantial improvement,
it is empirical and inadequate. An analytical solution should be more accurate
if the model is successfully drawn.
Let us therefore redefine the critical height, Hc, in terms of the unconfined
compressive strength, qu. At depth Hc in the face of a vertical cut,
qu ¼ Hc
Hc ¼ qu = ð21:23Þ
It will be noted that this is the critical height at failure and does not include a
factor of safety, yet is only one-half of that obtained with eq. (21.22), emphasizing
the importance of using an appropriate model for an analytical treatment.
Equation (21.24) still can be on the unsafe side if a trench cuts below a groundwater
table because seepage forces will be directed toward the the open cut. This situation
also is unsatisfactory from a construction standpoint. Restrictive regulations now
require that workers in a trench deeper than a specific height must be protected by a
steel ‘‘trench box,’’ or the sides of the trench must be cut back to a stable angle.
Example 21.6
A trench is to be excavated in soil that is reported to have an unconfined compressive
strength of 500 lb/ft2 (23.9 kN/m2) and a unit weight of 120 lb/ft3 (18.8 kN/m3) What is Hc
with a factor of safety of 1.2?
Answer: According to the qu criterion Hc ¼ 500/120 ¼ 4.2 ft, or 23.9/18.8 ¼ 1.27 m. The
allowable height is Ha ¼ 4.2 1.2 ¼ 3.5 ft (1.1 m).
21.3 LANDSLIDES
An active landslide should be checked for sounds of running water and for any
other evidence of broken water or sewer pipes indicated by welling-up of the
groundwater table and local wet spots. Nearby residents should be asked if they
have had any problems with sewer backups. Leakage is stopped by shutting the
water off.
Landslides often are spoon-shaped with the long dimension of the spoon oriented
across the slope. A three-dimensional analysis is possible but adds complexity to
an already complicated problem. If the landslide is wide, two or more traverses
may be investigated.
Borings, cone probings, and/or in-situ shear tests, discussed in Chapter 26, are
made at selected locations along a traverse. Difficult terrain requires the use of
a track vehicle. Shallow landslides may most conveniently be investigated by
hand augering. Undisturbed soil samples also may be taken for laboratory tests,
but it may be difficult or impossible to obtain duplicate samples, in which case
tests can be performed in situ.
Ground elevations and soil boring information are plotted to scale to develop
a cross-section, as shown in Fig. 21.7. There should be no vertical exaggeration
because slip angles may be measured from the graph. The cross-section should
show soil strata, the slip surface, and the position of the groundwater table.
Average unit weight of the sliding mass of soil is needed, and cohesions and
friction angles of the slip zone soil are required for an analysis.
Landslides generally are analyzed on an effective stress basis. Pore water pressures
are determined from groundwater levels in the borings with the aid of flow nets.
Pore pressures also can be monitored in the field with piezometers, although
active sliding can destroy a piezometer. Anomalies in groundwater elevations
can be critically important clues to broken sewers or water pipes, or to aquifers
and artesian effects, all of which can be causal factors.
The factor of safety is 1.0 when a landslide starts, and immediately decreases
as soil in the slip zone is remolded and loses strength. By the time a field inves-
tigation is performed the landslide probably will temporarily have stopped so the
factor of safety still is 1.0, or it may be slightly higher if the pause has allowed
water to drain out.
The analysis will provide valuable benchmark values that will enable predictions
of the effectiveness of different repair methods. The design factor of safety will
be influenced by drainage conditions. The worst-case scenario will involve full
saturation of the sliding mass.
Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com)
Copyright © 2007 The McGraw-Hill Companies. All rights reserved.
Any use is subject to the Terms of Use as given at the website.
Slope Stability and Landslides
Figure 21.7
Landslide cross-section divided into arbitrary vertical slices for analysis. The factors of safety for
individual slices shown at the bottom define active and passive parts of the landslide—the part doing
the pushing and the part being pushed. This is an important consideration for stabilization.
The analysis did not include vertical shearing stress, which adds to the factor of safety where the
slip surface curves.
The exposed part of a landslide slip surface at the top of a landslide is the
landslide ‘‘scarp.’’ The field investigation should include a walk along and above
the scarp to examine for tension cracks running parallel to the scarp, as these
normally will delineate the next part of the hill to go. The situation is aggravated
by surface water running down into the cracks, which therefore should be sealed
by cutting in a soil wedge with a spade.
crack prior to failure. Another evidence for progressive failure is that back-
calculation of the resisting force favors the residual, remolded shear strength
instead of the peak shear strength (Duncan and Wright, 2005). Thus, as one part
of the slide moves sufficiently to reduce the soil strength to a remolded state,
additional stress is transferred to the remaining areas until the slip zone becomes
continuous and the entire mass slides. The strength reduction also is confirmed
by Borehole Shear tests that inevitably show lower strength in the shear
zone compared with the same soil outside of the zone. Conventional analysis
methods nevertheless assume simultaneous shear all along a slip surface, and
progressive failure is taken into account by using the residual shear strength in
calculations. If residual strength data are not available, cohesion may be assumed
to be zero.
Example 21.7
What average thickness of slip zone can be expected in a landslide that has moved
30 ft (10 m)?
Answer: 0.3 ft ¼ 4 in. (0.1 m ¼ 100 mm), not enough for replicate samples. Adequate for a
Borehole Shear test.
8. Consider moving structures to a different site before it is too late and they
become involved in the landslide.
More permanent repair measures then can be applied, such as removing load
from the top, adding load at the toe, installing drains or structural measures,
or chemical stabilization in responsive soils.
Because the soil mass is assumed to rotate about a common center, solution
is based on moment equilibrium—that is, the sum of the acting forces along the
failure arc times the radius to the center equals the sum of the resisting forces
along the same arc times the same radius. Since the radius is constant, in polar
coordinates this is the equivalent of a force equilibrium of tangential forces.
Moment equilibrium still remains the basis for modern computer solutions,
including those for which the slip surface is not a circle and the radius is not
constant. However, there also are limitations to a moment equilibrium. For
example, it cannot be used when part of the slip surface is flat because the radius
is infinite. Then a force equilibrium is used, as in the preceding discussions.
The first step is to draw a cross-section of a slope (or landslide) to scale, and
include the actual slip surface or, if the slip surface is not known, a trial surface.
The mass of soil above the slip surface is then divided into vertical slices as in
Fig. 21.7, and forces on each slice are calculated using the same principles
and formulas as previously presented for a triangular mass of soil slipping along
a planar surface. Analysis involves summing resisting and acting forces on all
slices and dividing to obtain a factor of safety.
The relationships in eqs. (21.13) to (21.17) may be used to define acting and
resisting forces along the base of each slice. A factor of safety then may be
calculated that does not include side forces and can define active and passive
zones, as in Fig. 21.7. Horizontal side forces are required for equilibrium, and it is
these forces and resulting friction between slices that are ignored in the simplified
method of slices. The result is to underpredict the overall factor of safety,
particularly for deep slip surfaces and long slices.
Whereas in static equilibrium the pore water pressure is simply the depth below a
groundwater table times the unit weight of water, in a dynamic situation involving
seepage, the resistance to seepage causes a reduction in pore pressure. In Fig. 21.8
this is indicated by sloping equipotential lines. Therefore the pore water pressure
at point 8 on the slip surface does not equal the depth below the groundwater
table, but a depth Hw below the intersection of the equipotential line with the
groundwater table.
Figure 21.8
Flow net in a
landslide. Upward
seepage force in
the toe area further
destabilizes that
area, which may
be too weak to
support the weight
of a person.
Figure 21.9
Dimensions and
forces in the
simplified method
of slices.
Soil Data:
1. Slice number.
2. ¼ soil unit weight.
3. 0 ¼ soil friction angle measured on an effective stress basis.
4. c0 ¼ soil cohesion measured on an effective stress basis.
Calculations:
9. Weight W ¼ xz.
10. Acting force Sa ¼ W sin .
Summations: P
17. Sum of acting forces is Sa.
P
18. Sum of resisting forces is Sr.
P P
19. Global factor of safety ¼ FSt ¼ Sr/ Sa for the slope.
Notes:
Each of the steps can be set up as a column in a spreadsheet, and appro-
priate numbers and formulas copied vertically to enable a solution for each
slice.
In Step 2 the soil unit weight is the total unit weight including water.
In Steps 3 and 4 the soil cohesion and internal friction are on an effective
stress basis, which means that they either are measured with the soil in a fully
drained condition, or stresses have been corrected by subtracting pore water
pressures measured in a test specimen at the moment of failure.
In Step 8 the pore water pressure is measured from the equipotential line as
shown in Fig. 21.9, which is not the same as the depth below the phreatic
line because of seepage restraints of the water.
In Step 14 pore water pressure acts across the bottom of each slice, so it
must be multiplied by length L to give the force acting in opposition to N.
Steps 1–17 are calculated for each slice, and Steps 18–20 are summations that
are performed only once per analysis.
Example 21.8
Sum forces and calculate a factor of safety for the following slice:
1. Slice number ¼ 5.
2. ¼ soil unit weight ¼ 125 lb/ft3 or 19.6 kN/m3.
3. 0 ¼ 208.
4. c0 ¼ 500 lb/ft2 or 23.9 kN/m2.
5. X ¼ 10 ft or 3.05 m.
6. Z ¼ 16 ft or 4.88 m.
7. ¼ 328.
8. hw ¼ 8 ft or 2.44 m.
Calculations:
9. W ¼ XZ ¼ 10(16)(125) ¼ 20,000 lb/ft length, or 3.05(4.88)(19.6) ¼ 292 kN/m.
10. Sa ¼ W sin ¼ 20,000 sin 328 ¼ 10,600 lb/ft, or 292 sin 328 ¼ 155 kN/m.
11. N ¼ W cos ¼ 20,000 cos 328 ¼ 16,960 lb/ft, or 292 cos 328 ¼ 248 kN/m.
12. L ¼ X sec ¼ 10 sec 328 ¼ 11.8 ft, or 3.05 sec 328 ¼ 3.60 m.
13. N0 ¼ N hw w L ¼ 16,960 8ð62:4Þð11:8Þ ¼ 16,960 5891 ¼ 11,070 lb=ft, or
¼ 248 2:44ð9:81Þð3:60Þ ¼ 162 kN=m:
14. Sr ¼ c0 LþN0 tan 0 ¼ ð500Þð11:8Þ þ 11, 070 tan 20 ¼ 5900 þ 4029
¼ 9930 lb=ft, or ¼ ð23:9Þð3:60Þ þ 162 tan 20 ¼ 86 þ 59 ¼ 145 kN=m:
15. FSi ¼ 9930/10,600 ¼ 0.94, or 145/155 ¼ 0.94.
Figure 21.10
Diagram
to illustrate
equivalence of
seepage force and
effective stress
analysis.
Sa ¼ sin
N ¼ cos
The normal force on an effective stress basis is
From the equipotential line u ¼ w hw, where hw ¼ ab cos and ab ¼ cos . Then
U ¼ w cos2 = cos ¼ w cos
N0 ¼ cos w cos ¼ ð w Þ cos
Sr ¼ N0 tan 0 ¼ ð w Þ cos tan 0
which is the equation that was developed on the basis of seepage forces, with
the notation on 0 indicating an effective stress basis. This illustration does not
include soil cohesion as cohesion is not affected by seepage forces or buoyancy.
Example 21.9
Use the buoyant unit weight to estimate a fully saturated factor of safety for the slice in the
preceding example, using 125 lb/ft3 (19.6 kN/m3) as the saturated unit weight.
Answer:
Sa ¼ W sin ¼ 10,600 lb=ft or 155 kN=m ðunchangedÞ
W 0 ¼ XZ 0 ¼ ð10Þð16Þð125 62:4Þ ¼ 10,020 lb=ft ðbuoyant weightÞ or
¼ ð3:05Þð4:89Þð19:6 9:81Þ ¼ 146 kN=ft
N0 ¼ 10,020 cos 32 ¼ 8490 lb=ft or 146 cos 32 ¼ 124 kN=m ðbased on buoyant
weightÞ
cL ¼ 5900 lb=ft or 86 kN=m ðunchangedÞ
Sr ¼ 5900 þ 8490 tan 20 ¼ 8990 lb=ft or 89 þ 124 tan 20 ¼ 134 kN=m
FSi ¼ 8990=10,600 ¼ 0:85 or 134=155 ¼ 0:86
compared with the previous answer of 0.94 for a partially saturated condition.
Example 21.10
Calculate a fully saturated factor of safety for the slice in the preceding example using
effective stresses.
Answer:
Sa ¼ 10,600 lb=ft or 155 kN=m ðunchangedÞ
hw ¼ 16cos2 32 ¼ 11:5 ft or 4:89 cos 232 ¼ 3:52 m
u¼ 11:5ð62:4Þ ¼ 718 lb=ft2 or 3:52ð9:81Þ ¼ 34:5 kN=m2
N0 ¼N uL¼ 16,960 ð718Þð11:8Þ ¼ 8187 lb=ft or
¼ 248 ð34:5Þð3:60Þ ¼ 124 kN=m
Sr ¼ c0 LþN0 tan 0 ¼ ð500Þð11:8Þ þ 8187 tan 20 ¼ 5900 þ 3089 ¼ 8989 lb=ft or
¼ ð23:9Þð3:60Þ þ 124 tan 20 ¼ 86:0 þ 45:1 ¼ 131 kN=m
FSi ¼ 8989=10,600 ¼ 0:85 or 131=155 ¼ 0:85
In this case the solution involving submerged unit weight is simpler because of possible
ambiguity in evaluating hw.
Table 21.1
Taylor stability numbers for various failure surfaces for soil having a friction angle ¼ 258. The higher
the number, the lower the indicated stability. Answers that are in close agreement are bracketed—but
they still only apply to simultaneous failure along the entire slip surface
Slope i, degrees Surface Method s , degrees
Because the obliquity and the factor of safety affect both the resisting force, the
solution is by trial and error. Bishop proposed a simplified method that rapidly
converges. The method was developed specifically for circular slip surfaces.
Figure 21.11
The simplified
Bishop method
assumes a
circular slip
surface and
adjusts the factor
of safety
for a partially
developed
resisting force Sr.
Sr sin
N cos ¼ W ð21:29Þ
FSi
Thus as the factor of safety increases, N decreases. By assuming that the factor
of safety is constant along the entire slip plane, Bishop obtained the following
expression:
X c0 l cos þ ðW ul cos Þtan 0 X
FS ¼ W sin ð21:30Þ
cos þ ðsin tan 0 Þ=FS
This is the same as the Ordinary Method of Slices. Bishop’s method does not
specify the radius or position of the center of the most critical circle, which are
determined by repeating the calculations for different trial surfaces and selecting
the one having the lowest factor of safety. This is quickly done with a computer
that can be programmed to find the 10 most critical circles. An example solution is
shown in Fig. 21.12, and it will be seen that a diversity of circular shear paths gives
almost the same factor of safety. The actual path therefore can be expected to be
sensitive to small variations in soil strength but can deviate from the most critical
circle without greatly influencing the factor of safety.
Figure 21.12
Trial circular
surfaces from a
STABL computer
run using the
simplified Bishop
method of slices.
All factors of safety
are in a range
1.5 percent, the
most critical being
shown with a
heavy line.
Figure 21.13
Spencer’s (1967)
charts for
estimating the
factor of safety of
a slope in
homogeneous soil
with selected pore
pressure ratios.
Incorporation of pore water pressure into the analysis was simplified by defining a
pore pressure ratio, defined as the ratio of pore water pressure at the base of each
slice to the weight of the slice:
u
ru ¼ ð21:31Þ
h
where u is the pore water pressure at the slip surface, and h represents the
overburden pressure.
Spencer’s method was derived for a circular slip surface but also can be extended
to noncircular surfaces. Duncan and Wright (2005) indicate that this is the
simplest complete solution for estimating a factor of safety.
Example 21.11
Estimate the factor of safety for a 35 ft (10.7 m) high slope at 1 vertical to 3 horizontal, in
soil having a friction angle of 258, cohesion of 300 lb/ft2 (14.4 kN/m3), and unit weight of
120 lb/ft3 (18.8 kN/m3). Assume a worst saturation condition.
Answer: The slope angle is i ¼ tan1 0.333 ¼ 18.48. For full saturation, the lower graph in
Fig. 21.11 gives
c 300 14:4
¼ ¼ 0:028, or ¼
FH Fð120Þð35Þ Fð18:8Þð10:7Þ
where F is the factor of safety. Solving for F gives 2.5.
Answer: If c ¼ 0 the failure is parallel to the soil surface and Spencer’s charts
cannot be used. The factor of safety for a continuous slope in cohesionless soil is
discussed at the beginning of this chapter. From eq. (21.5), FS ¼ tan /tan i ¼ 1.4.
and the books are allowed to rotate, there is no book-to-book friction, but if
the curvature changes, individual books must slide up or down relative to one
another. This vertical resistance to slipping adds to the factor of safety.
Figure 21.14
Janbu’s empirical
correction to
factors of safety to
account for
mobilization of
side forces.
Example 21.12
An analysis of the landslide in Fig. 21.7 by the ordinary method of slices gave a factor of
safety of 1.0 but the slide temporarily was stable. What is the Janbu correction?
Answer: d/L is approximately 0.1, so the factor of safety should be multiplied by 1.05.
Slope stability analyses require averaged cohesions and angles of internal friction
that can be obtained only with large numbers of tests that are specific for these
values. The problem is both simplified and made more difficult in the case of an
active landslide—simplified because remolding in the shear zone tends to move
shearing strengths toward an average so that fewer tests are required, and more
difficult because the shear zone usually is too thin to obtain multiple identical
samples for laboratory testing. In that case a stage test, either a triaxial stage test
or the in-situ Borehole Shear test, may be the only practical solution.
Where adequate soil data are lacking, as often is the case, analyses sometimes are
performed using assumed values for cohesion and friction angle, or values obtained
on the basis of empirical correlations such as shown in Fig. 19.2. The degree of
sophistication of a computer program then has little to do with the reliability of the
results, as any analysis based on a guess remains a guess. Chapter 27 presents
somewhat more reliable estimates of friction angle based on penetration test data.
thixotropic recovery if the landslide stops. Still, it is on the safe side to assume that
it is zero.
In a stable slope, additional unknowns include location and shape of the most
critical slip surface, location and anticipated changes in the groundwater table,
and indeterminacy of the analysis. Because of these additional unknowns the
design factor of safety usually is of the order of 1.2 or more.
As the number of computer runs increases, the computer can produce some
unrealistic quirks in the slip surface. These probably relate to the assumptions made
A unique property of the log spiral is that the angle between the radius and the
spiral is constant. If that angle is assigned to the angle of internal friction,
any radius for the spiral defines the direction of the resultant force at failure.
The corresponding equation is
Figure 21.15
Ideal slip surfaces
and (a–c) some
developments
involving cycloids.
weight of a soil mass, such as in slopes or the active side of retaining walls.
From eq. (21.32) and Fig. 21.15 will be seen that a circle is a special case of
a log spiral with ¼ 0.
21.7.3 Cycloids
A cycloid is generated by a point on the rim of a wheel that is moving
along a horizontal surface (Fig. 21.15(a)). This does not appear particularly
relevant, but the cycloid has some other features that have interesting parallels
for slope stability.
As pointed out by Fox (1981), the cycloid was shown by Bernoulli to be the
solution to the ‘‘brachistochrone’’ problem that had intrigued Galileo: what is the
shape surface that will provide the fastest and most efficient path for a ball rolling
down a slope? The problem is illustrated in Fig. 21.15(b). Bernoulli used an
analogy to solve the problem by noting a similarity to the path of a beam of light
penetrating at an angle through a medium with a gradually increasing index of
refraction (Fig. 21.15(c)). The solution was obtained by integrating Snell’s Law of
refraction. There also is a striking analogy in soil mechanics to a shear path that
encounters increasing resistance with depth because of the increasing overburden
pressure and friction. Although analogy can lead to new innovations and insights
(Koestler, 1964), it is important to remember that, regardless of the depth of
philosophical conjecture, it does not constitute a proof. However, the analogy
that can be reinforced by reasoning can come closer.
Conceivably the energy-efficient cycloidal path for a rolling ball might also
describe the most energy-efficient path for the slip surface in a landslide. This
would be consistent with a concept of progressive failure, as each new segment of
a slip surface searches for the most energy-efficient path. Fox (1981), who pointed
out many of the features in this discussion, also pointed out that catenary curves,
which according to arching theory are trajectories of major principal stress,
for practical purposes are orthogonal to cycloids.
The x and y parameters of a cycloid with a starting point at the lower end
and verticality at the upper end are
x ¼ að þ cos Þ and y ¼ að1 cos Þ ð21:33Þ
where is in radians.
21.8.2 Drainage
Drainage is one of the oldest methods for stopping a landslide, and is often
accomplished with corrugated plastic pipe that has been perforated with narrow
slots to let water in. However, the slots also can let water out, and a tile line that
crosses an active slip zone will be stretched, bent, compressed, or sheared off,
Figure 21.16
Some repair
methods for
landslides. The
one at the bottom
is too frequently
used and can only
make matters
worse.
so it may drain directly into the slip zone. Trenching to install the tile is a risky
adventure because it can reactivate the slide. Trenching can be done only during
dry periods when the water table is low, and should be done so that only short
sections are open at any one time.
Directional drilling such as shown in Fig. 21.17 is expensive, but can extend much
deeper than is possible with trenching, and avoids the dangers associated with
open trenches. Directional drilling is accomplished in soils by pushing a beveled
mandrel using water jetting. The mandrel carries electronic direction sensors so
that directions can be changed by changing the orientation of the bevel. A small-
diameter perforated plastic pipe is attached to the mandrel at the exit point, and
then pulled back and left in the boring.
Gravity drainage requires an exit point, and the groundwater table also can
be lowered by the use of well points, vacuum well points, or electro-osmosis.
The latter processes are relatively expensive but may be the only effective ways
to drain low-permeable clay soils. Vertical wells can be effective but require
intermittent pumping.
A key trench filled with granular material at the toe ideally will cross the slip zone
and act as a drain. Another advantage is that part of the soil at the toe is replaced
with soil having a higher friction angle, a factor that is readily included in stability
calculations. However, a key trench is not necessarily a perfect solution because
the landslide still can go over the top.
Figure 21.17
Directional drilling
can install
perforated plastic
drain lines without
trenching.
Example 21.13
The factor of safety of an active landslide is assumed to 1.0 and with full saturation the
friction angle in the slip zone is back-calculated to be 128. What percent replacement will be
needed to increase the factor of safety to 1.2, (a) using stone columns having a friction
angle of 358, (b) using Rammed Aggregate Piers having a friction angle of 488?
Answer:
ðaÞ Stone columns: m tan 35 þ ð1 mÞ tan 12 þ ¼ 1:2 tan 12
0:70 m þ 0:21 0:21 m ¼ 0:26
m ¼ 0:10 or 10%
ðbÞ Rammed Aggregate Piers: m tan 48 þð1 mÞ tan 12 ¼ 0:26
1:1m þ 0:21 0:21m ¼ 0:26
m ¼ 0:06 or 6%
handling the chemical to prevent personal injuries. If the soil is too soft to hold
an open boring, reverse augering is used to push lime down into the soil, and
immediately creates a sheath of stabilized soil. In extreme cases lime can be
introduced in holes created with a concrete vibrator, which has been used to save
houses that were audibly creaking from pressures from landslides.
The drilled lime process works by immediately drawing water from the soil to
hydrate the lime, and sliding immediately stops. However, this reduction in water
content of the soil is only temporary. The lime column approximately doubles in
diameter as the lime hydrates, drying and creating radial tension cracks in the
surrounding soil. The cracks then are filled by a paste of hydrated lime, greatly
enhancing its distribution into the soil (Handy and Williams, 1967). For eco-
nomical reasons borings may only cover the upper part of the active zone so long
as no other changes are made that would reduce the factor of safety. The boring
diameters and spacing are designed to introduce from 1 to 3 percent lime,
depending on the clay content of the soil. Excess lime is used in a pozzolanic
reaction with soil clay minerals that creates hydrated silicates and aluminates
that are the same products as in hydrating Portland cement.
A simple field test to determine applicability of the drilled lime method for a
particular soil is to place small quantities of soft, saturated soil from the slip zone
into two sealable plastic bags. A few percent hydrated lime is added to soil in one
of the bags, and both are sealed and the samples remolded by hand. Within a few
minutes a reactive soil will become crumbly and appear to dry out even though no
moisture has been lost, as the plastic limit increases and exceeds the soil natural
moisture content. If that does not occur, either the moisture content is too high or
the soil does not contain the necessary active clay minerals.
21.9.1 Relevance
Earth dams vary from some of the smallest to some of the largest man-made
structures in the world. Design of large dams takes on additional importance
because of safety issues. Many topics such as compaction, settlement, slope failure,
dispersive soils, piping, etc., are covered elsewhere in this text. Because of the
severity of failures and the continuing construction of new earth dams, particularly
in developing countries, the topic merits some additional consideration.
Figure 21.18
Causes of failures of earth dams. Slope failures include those occurring both in
the dam and around the reservoir; internal failures include piping and improper
compaction. (From data of Milligan, 1999.)
the consequences are not nearly so severe, and there typically is a lower level of
care in design and construction. Small dams also are much easier to repair and
restore to operating condition.
Most failures occurred before the dams were 10 years old, perhaps in part because
soil aging results in a gradual gain in shear strength. A statistical study indicates
that the odds of failure of a dam are about 1 in 7000 or 0.014 percent for any
particular year, but the odds are influenced by seismicity as well as the size of the
dam. In the typical 100-year useful life of a dam the probability is 1.4 percent,
but the odds are improved after the dam has survived its first decade.
The pie chart in Fig. 21.18 indicates that slope failures are not the major cause
of distress in earth dams, being eclipsed by foundation failures and overtopping.
Nevertheless the basic tenets are followed regardless of size: using properly
selected soil, having an adequate site preparation, and compaction of the
embankment at a moisture content that is at or above the optimum. Some
problems specific to earth dams and levees are discussed below.
Upstream faces of earth dams are covered with coarse stone, broken concrete,
stabilized soil, or grout-filled plastic mattresses to prevent wave erosion. The dam
obviously must be high enough above the reservoir level to prevent overtopping
by waves, whose height depends on wind velocity and the distance the wind blows
across water, or ‘‘fetch.’’ The influence of fetch is illustrated by a guide to
freeboard for small dams (Fig. 21.19).
Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com)
Copyright © 2007 The McGraw-Hill Companies. All rights reserved.
Any use is subject to the Terms of Use as given at the website.
Slope Stability and Landslides
Figure 21.19
Suggested
freeboards for
earth dams.
(Adapted from
U.S. Bureau of
Reclamation,
1960.)
Figure 21.20
Rate of movement
of the ground
surface along a
traverse running
downslope on the
Vaiont landslide.
(Modified after
Leonards, 1987,
p. 608.)
Figure 21.21
Some
representative
earth dam designs
to direct and
control seepage.
Riprap protects the
faces from erosion;
berm adds to
stability of the
downstream slope.
(Adapted from
Creager et al.,
1948.)
The resulting erosion and piping can have much more serious consequences than
simply a leak.
The tip of a hydraulic fracture becomes a focal point for tension that can be
calculated from a stress concentration factor from elastic theory and mechanics of
materials. The smaller the radius at the tip of the fracture, the higher the tensile
stress. Hydraulic fracturing is prevented by equalizing seepage pressure with
pore water pressure, which means slow first-filling of a reservoir. Piezometers
monitor pore water pressures within the dam and also preferably within the
abutments to determine if they are as predicted from the flow net—a low pressure
being conducive to hydraulic fracturing, and a high pressure indicating that
fracturing already may have occurred. Piezometers are monitored throughout the
life of a large dam to show any deviations from expected behavior.
This potential for leakage through soils and rocks underneath a dam is evaluated
with pressure tests performed in borings, and these are an important part of site
exploration. Underdam leakage is most effectively prevented by use of a cutoff wall
that can be concrete, clay, or driven sheet pile. Grouting is frequently used but is
difficult to direct to the desired area, and in severe cases of underdam leakage the
grout is washed away as fast as it is injected, particularly after a reservoir is filled
with water. In some cases the volume of grout pumped has exceeded the volume of
the dam itself, resulting in a huge cost overrun. Leaky reservoirs also impact an
entire neighborhood by raising the groundwater table and flooding nearby fields
and basements.
In recognition that a river often flows on a bed of sand that grades downward into
gravel, an upstream clay blanket may be used to seal off the river channel prior to
reservoir filling. If the hydraulic conductivity of the potential seepage layer is 10
times that of the average soil underneath the dam, the clay blanket can extend a
distance upstream equal to 10 times the width of the dam to achieve the same
resistance to flow. Some common methods for controlling through-dam and
under-dam seepage are illustrated in Fig. 21.21.
Vaiont Dam occupies a deep valley in the Dolomite Alps in northern Italy, and is
one of the highest concrete arch dams in the world, standing tall at 860 ft (262 m).
The landslide area was about 2 km wide by 1.5 km long (1.2 0.9 miles).
The speed of sliding was back-calculated to be about 110 km/hr (70 mph), and
a combination of wind and water from the slide blew out windows in a town
100 m above the reservoir on the opposite hillside. The failure in 1963 remains
one of the worst dam tragedies to date.
The landslide at least in part developed along a prehistoric slip zone, and as
the lake was being filled a smaller landslide developed in the same area. During
the first filling of the reservoir livestock on the hillside were observed to be very
nervous, probably as a result of rock noises that can be a prelude to shear
failure but are inaudible to the human ear. Ground surveys were conducted
and showed that the hillside was moving at a very slow rate, as shown at the left
in Fig. 21.20. Attempts were made to reduce the rate of movement by reducing
the lake level, but the lake level ominously kept rising as the hill squeezed in on the
reservoir.
The engineers then were faced with a serious dilemma: should they increase the
rate of drainage, which would create a rapid drawdown condition, or should they
try and maintain the reservoir at a constant elevation, and/or should they require
an evacuation? There were no easy answers, particularly as attempts to lower or
maintain the reservoir level were frustrated by movement of the landslide.
The data in Fig. 21.20 present a classic example of different stages of soil creep,
although this may not have been fully appreciated at the time. During primary
creep, which is the mechanism of secondary consolidation, the rate is decreasing
and may stop. During secondary creep the rate is constant, and during tertiary
Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com)
Copyright © 2007 The McGraw-Hill Companies. All rights reserved.
Any use is subject to the Terms of Use as given at the website.
Slope Stability and Landslides
creep the rate is increasing and is preliminary to rapid shear failure. From the
graph it appears that tertiary creep initiated about 2 months prior to failure and
was assured by a month prior to failure, which is when warnings would have been
appropriate.
The sudden and rapid rate of sliding was not expected, and it later was conjec-
tured that heat generated from friction might have contributed by decreasing the
viscosity. Steam pressure also was suggested but is unlikely under the huge
confining pressure. A more likely explanation for the speed of the sliding comes
from an examination of measurements at different points on the sliding mass.
These indicated that the upper part of the slide had moved 4.03 m while the lower
point moved only 2.87 m. This indicates that the rock mass between those two
points had compressed 1.16 m, confirming the thesis that a landslide has an active,
pushing part and a passive part that is being pushed. The separation distance
between the points was about 700 m, so the compressive strain was about
0.17 percent as stored strain energy. Being elastic, that energy would be subject
to instant release and conversion to acceleration.
Several approaches might have helped avert this tragedy, including a better
appreciation of the geological factors contributing to the slide. Modern
instrumentation, for example the use of inclinometers and acoustic micro-
phones, could have provided a more stringent warning. The cows had the
right idea.
Problems
21.1. Describe two situations when the factor of safety of a landslide is exactly 1.0,
and explain how this observation can be used to estimate soil sensitivity.
21.2. Explain the mobilization of shearing stresses within a sliding mass of soil,
and two situations where they theoretically will not be mobilized.
21.3. Analysis of an active landslide yields a factor of safety of 2.8. Is that pos-
sible? In your opinion is this more likely to be due to erroneous data or to the
method of analysis? What procedures should be used to check the results?
21.4. The friction angle of a soil is 258. Calculate the angle of repose on the
Moon where gravity is about one-fifth that on Earth.
21.5. A deep railroad cut intersects the plane of contact between a loess soil and
an impermeable clay paleosol. The plane of contact makes an angle of 128
with horizontal. Soil information is given in Fig. 21.22. Assume that
cohesion is zero and calculate the factor of safety with no contribution from
seepage.
21.6. Repeat Problem 21.5 assuming a worst-case scenario for seepage.
21.7. In the preceding problem even though there is a plane along which sliding
may occur, there is no assurance that a more critical curved slip surface
Figure 21.22
Conditions for
Problems 21.5 and
21.6.
may not exist in the soil above that plane. Use Spencer’s method with three
different degrees of saturation to estimate factors of safety. Which do you
feel is most appropriate for design of a roadcut along a county road? Of a
cut bank behind a house?
21.8. The stability of a slope is being analyzed by the method of slices. On a
particular curved surface through the soil mass, the shearing component of
the weight of each slice and the normal component of the reaction at the base
of each slice are presented in Table 21.2. The length of the curved surface is
29.9 m (98 ft). If the friction angle of the soil is 128 and the cohesion is
2.39 kPa (50 lb/ft2), what is the factor of safety along this particular curved
surface?
21.9. A 21 m (70 ft) roadcut has side slopes of 2 to 1. The base of the cut is soil
in limestone. The soil has the following properties: ¼ 19.6 kN/m3
(125 lb/ft2); c 0 ¼ 26.8 kPa (560 lb/ft2); 0 ¼ 208. Draw a trial circular failure
with the center at the crest of the cut and intersecting the cut at the toe. Use
the simplified method of slices to determine the minimum factor of safety
if seepage is not a factor, and apply the Janbu correction.
21.10. Recalculate the factor of safety in Problem 21.8 using Spencer’s method.
Figure 21.23
Landslide problem
for the method of
slices data in
Tables 21.3 and
21.4.
Table 21.4
Spreadsheet calculations of resisting forces in Fig. 21.23
Slice L c cL N u uL N0 N 0 tan ’ Sri FSi
No. (m) (kPa) (MN/m) (MN/m) (kPa) (MN/m) (MN/m) (deg.) (MN/m) (MN/m) (Sr/Sa)
Note : Negative FSi for individual slices results when Sa and Sr act in the same direction.
21.19. The slip surface of an active landslide is planar at the top and thickens to
more than a meter downslope. Explain.
21.20. Look to sources to review geological factors that may have influenced the
failure at Vaiont Dam.
21.21. Prepare a table with columns for quicksand, dispersive clay, and hydraulic
fracturing and entries describing (1) flow rate required, (2) flow direction,
(3) factors governing hydraulic pressure, (4) prevention, (5) high, medium,
or low potential danger to a dam, and (6) clues that the phenomenon is
occurring.
21.22. What should be looked for when inspecting an earth dam? When should
inspection be most rigorous, during high or low level of the reservoir? How
often should an earth dam be inspected (a) the first 10 years, (b) after 10
years?
Logani, K. L. (1973). ‘‘Dilatancy Model for the Failure of Rocks.’’ Unpublished. Ph.D.
dissertation, Iowa State University, Ames, Iowa.
Milligan, V. (1999). ‘‘An Historical Perspective of the Development of the Embankment
Dam.’’ Canadian Dam Association. Conference, Sudbury, Ontario.
Morgan, A. E. (1971). dams (sic) and other disasters. Porter Sargent Publisher, Boston.
Scholtz, C. H. (1967). ‘‘Wear and Gouge Formation in Brittle Faulting.’’ Geology 15,
493–495.
Seed, H. B. (1968). ‘‘Landslides During Earthquakes Due to Soil Liquefaction.’’
Fourth Terzaghi Lecture. ASCE J.Soil Mech. and Foundation Eng. Div. 94(SM5),
1053–1122.
Skempton, A. W. (1964). ‘‘Long Term Stability of Clay Slopes.’’ Geotechnique 14(2),
75–102.
Spencer, E. (1967). ‘‘A Method of Analysis of the Stability of Embankments Assuming
Parallel Interslice Forces.’’ Ge´otechnique 17(1), 11–26.
Taylor, D. W. (1948). Fundamentals of Soil Mechanics. John Wiley & Sons, New York.
Terzaghi, K. (1943). Theoretical Soil Mechanics. John Wiley & Sons, New York.
U.S. Department of Interior Bureau of Reclamation (1960). Small Dams. U.S. Govt.
Printing Office, Washington, D.C.