Causes of Earth Dam Failure: Seepage, Overtopping and Slope Stability (2026 Guide)

1. Introduction to Embankment Dam Safety and Failure Statistics

Embankment dams comprise more than $75\%$ of all major water-retention structures globally. Constructed using naturally occurring soil, gravel, and rockfill materials, these gravity structures rely on self-weight and material shear strength to resist hydrostatic reservoir forces. However, unlike rigid concrete gravity dams, earthen embankments are inherently vulnerable to internal water action and external erosion.

Historical forensic records cataloged by the International Commission on Large Dams (ICOLD) reveal that understanding the causes of earth dam failure is central to preventing catastrophe. A breach releases uncontrolled reservoir volumes, generating destructive flash floods that endanger downstream populations, infrastructure, and ecological corridors.

Each mechanism develops through distinct physical transitions. Hydraulic failures progress rapidly over hours, while seepage piping can operate undetected internally for months before culminating in rapid breach formation.

3. Hydraulic Failure Modes: Overtopping, Spillway Inadequacy, and Wave Erosion

3.1 Overtopping Breach Mechanisms and Hydrodynamic Shear

Overtopping accounts for approximately $30\% – 35\%$ of all recorded earthen dam breaches. When incoming flood inflows exceed spillway discharge capacity, reservoir water rises above the impervious crest.

Water cascading down the unprotected downstream embankment slope exerts tractive shear stress $\tau_b$:

$$\tau_b = \gamma_w \cdot y \cdot S_0$$

Where:

• $\gamma_w$ = Unit weight of water ($9.81 \text{ kN/m}^3$)

• $y$ = Flow depth over the slope ($\text{ m}$)

• $S_0$ = Downstream slope gradient ($\sin\beta \approx \tan\beta$).

When bed shear $\tau_b$ exceeds the critical shear stress $\tau_{cr}$ of the downstream embankment fill, soil particles detach. A headcut knickpoint forms on the downstream face and migrates upstream toward the reservoir crest:

Reservoir Water Level > Dam Crest
~~~~~~~~~~~~~~~~~~\
                   \========= [CREST]
                    \       \
                     \       \---- Hydrodynamic Flow Cascade
                      \           \
                       \           \==== Headcut Knickpoint Formation
                        \                \---> Knickpoint migrates UPSTREAM
                         \________________\

Once the headcut breaches the reservoir crest edge, uncontrolled discharge erodes a deep trapezoidal breach channel. The outflow hydrograph peaks rapidly, causing catastrophic inundation downstream.

3.2 Wave Action, Freeboard Depletion, and Surface Scour

High wind velocities across long reservoir fetches generate continuous surface waves. Without adequate riprap stone armoring, breaking waves batter the upstream face, washing away fine cohesive particles and undercutting the slope.

additionally, prolonged wind setup combined with wave run-up reduces design net freeboard:

$$F_{net} = H_{crest} – (H_{max\_pool} + S_{setup} + R_{runup})$$

If $F_{net} \le 0$, wave spray and surging wash over the crest, softening compacted crest gravel and initiating localized progressive rilling.

4. Seepage Failure Modes: Internal Erosion and Piping Mechanics

4.1 Backward Erosion Piping and Exit Hydraulic Gradients

Internal erosion and piping represent the most insidious causes of earth dam failure, responsible for over one-third of all historical catastrophes, including the famous Teton Dam failure in 1976.

Seepage water percolating through the porous dam core exits at the downstream toe or into unprotected foundation joints. As seepage discharge concentrates, the seepage force per unit volume $j$ acts in the direction of flow:

$$j = \gamma_w \cdot i$$

Where $i = \frac{\Delta h}{\Delta L}$ represents the hydraulic gradient.

At the downstream exit face, upward seepage force directly opposes the submerged effective weight of the soil skeleton $\gamma’$. When the exit gradient reaches the critical hydraulic gradient $i_{cr}$, effective stress drops to zero:

$$\sigma’ = \sigma – u = 0$$

The critical hydraulic gradient derives directly from soil phase relations:

$$i_{cr} = \frac{\gamma’}{\gamma_w} = \frac{G_s – 1}{1 + e}$$

Where:

• $G_s$ = Specific gravity of soil solids (typically $2.65 – 2.72$)

• $e$ = Void ratio of the soil matrix.

When actual exit gradient $i_{exit} \ge i_{cr}$, soil particles boil and wash out. A continuous open pipe erodes backward along the seepage flowline toward the reservoir:

               [ Upstream Reservoir ]
                     |~~~~~~~~|
                     |        |====== [IMPERVIOUS CORE]
                     |        |     \
                     |        |      \--- Phreatic Seepage Line
                     |        |       \
                     |        |        \====== [Backward Piping Pipe] <=== Exits at Toe
                     +--------+-----------------\
                                                 \---> Sand Boils Form Here!

4.2 Internal Suffusion, Dispersive Clays, and Contact Erosion

• Suffusion: In gap-graded soils, fine particles migrate through the pore channels formed by the coarse gravel skeleton under steady seepage, increasing void ratio and local permeability without overall volume change.

• Suffosion: Suffusion accompanied by volumetric collapse of the coarse soil framework.

• Dispersive Clay Erosion: Highly sodium-saturated clays ($\text{ ESP} > 6\%$) deflocculate in low-salinity water, dissolving colloidal clay particles into suspension along cracks even under very low hydraulic gradients ($i < 0.05$).

• Conduit Contact Erosion: Inadequate compaction around rigid concrete outlet conduits creates annular voids. Seepage water exploits this smooth interface, eroding a continuous passage along the conduit exterior.

4.3 Terzaghi and Sherard Filter Design Criteria

Engineers prevent internal erosion by installing graded granular filter and drainage zones adjacent to the impervious core. Modern filter criteria follow Terzaghi and Sherard rules:

1. Retention (Piping Prevention) Criterion: Pore openings in the filter must be small enough to block coarse particles of base soil:
   $$\frac{D_{15\text{ (filter)}}}{d_{85\text{ (base)}}} \le 4 \text{ to } 5$$
2. Permeability (Drainage) Criterion: Filter must be porous enough to dissipate pore pressures without building excess hydraulic head:
   $$\frac{D_{15\text{ (filter)}}}{d_{15\text{ (base)}}} \ge 4 \text{ to } 5$$
3. Segregation Avoidance: Filter uniformity coefficient $C_u = \frac{D_{60}}{D_{10}} \le 6$ to prevent particle segregation during placement.

5. Structural and Slope Stability Failure Modes

5.1 Upstream Slope Failure Under Rapid Drawdown

When a reservoir level drops rapidly (e.g., during emergency discharge or heavy irrigation demand), water drains from the upstream reservoir faster than pore water can escape from low-permeability upstream clay shell fill.

       Rapid Drop in Reservoir Head (hw2 << hw1)
                 \
                  \     [Upstream Shell]
                   \       /~~~~~ [Unbalanced High Pore Pressure u]
                    \     /     \
                     \   /       \====== Critical Slip Surface (Downhill)
                      \_/_________\

The stabilizing hydrostatic water pressure on the upstream face vanishes instantly. However, internal pore pressure $u$ remains high. Effective normal stress $\sigma’ = \sigma – u$ plummets. Shear strength along potential circular slip surfaces drops:

$$\tau_f = c’ + (\sigma_n – u) \tan\phi’$$

Driving gravitational shear stresses exceed available shear resistance, triggering massive rotational slumps of the upstream embankment shell.

5.2 Downstream Steady-State Seepage Slope Instability

Under prolonged normal operating pool levels, steady-state seepage saturates the lower portion of the downstream shell. The phreatic surface intersects the downstream slope if internal chimney/blanket drains are absent or clogged.

High pore water pressures reduce effective stresses along deep circular and non-circular slip surfaces. Limit equilibrium formulations, such as Bishop’s Simplified Method, evaluate the Factor of Safety ($FS$):

$$FS = \frac{\sum \left[ \frac{c’ b_i + (W_i – u_i b_i) \tan\phi’}{m_\alpha} \right]}{\sum W_i \sin\alpha_i}$$

Where:

• $m_\alpha = \cos\alpha_i \left(1 + \frac{\tan\alpha_i \tan\phi’}{FS}\right)$

• $W_i$ = Slice total weight

• $u_i$ = Pore pressure at base of slice

• $\alpha_i$ = Slice base inclination angle.

Design codes mandate $FS \ge 1.50$ for steady-state downstream seepage conditions.

5.3 Foundation Shear Failure and Liquefaction Under Dynamic Shaking

Embankment dams resting on soft alluvial clays or loose, saturated silty sand deposits face foundation shear failure:
1. Spreading Along Weak Seams: High lateral earth pressures from the heavy embankment push weak foundation clay layers horizontally, initiating base extrusion.
2. Dynamic Liquefaction: Earthquake shaking induces cyclic shear strains in saturated, loose granular materials ($N_{SPT} < 10$). Excess pore water pressure rises until $\Delta u = \sigma_v$, transforming solid soil into a viscous slurry. The dam crest subsides abruptly, causing catastrophic overtopping.

6. Step-by-Step Solved Geotechnical Calculation Examples

6.1 Example 1: Critical Hydraulic Exit Gradient and Factor of Safety Against Piping

  Dam Cross-Section & Flow Net:
  Upstream Head H1 = 28.0 m ------------> Downstream Head H2 = 2.0 m
  Total Head Loss ΔH = 26.0 m
  Number of Equipotential Drops Nd = 14
  Length of Final Flow Element ΔL = 1.85 m

Problem Statement:

An earth dam retains a reservoir with total differential head $\Delta H = 26.0 \text{ m}$. A flownet analysis yields:

• Total number of equipotential drops: $N_d = 14$

• Length of the exit flow element at the downstream toe: $\Delta L = 1.85 \text{ m}$

• Soil solids specific gravity: $G_s = 2.68$

• Soil in-situ void ratio: $e = 0.65$

• Evaluate the critical hydraulic gradient, actual exit gradient, and factor of safety against backward piping.

Step 1: Compute Critical Hydraulic Gradient ($i_{cr}$)

Using the Terzaghi phase relationship:

$$i_{cr} = \frac{G_s – 1}{1 + e} = \frac{2.68 – 1.00}{1 + 0.65} = \frac{1.68}{1.65} \approx 1.0182$$

Step 2: Calculate Head Loss per Equipotential Drop ($\Delta h$)

$$\Delta h = \frac{\Delta H}{N_d} = \frac{26.00 \text{ m}}{14} = 1.8571 \text{ m}$$

Step 3: Compute Actual Exit Hydraulic Gradient ($i_{exit}$)

The exit gradient across the final element of length $\Delta L = 1.85 \text{ m}$ is:

$$i_{exit} = \frac{\Delta h}{\Delta L} = \frac{1.8571 \text{ m}}{1.85 \text{ m}} = 1.0038$$

Step 4: Calculate Factor of Safety Against Piping ($FS_{\pi\ping}$)

$$FS_{\pi\ping} = \frac{i_{cr}}{i_{exit}} = \frac{1.0182}{1.0038} = 1.014$$

Engineering Evaluation:

The calculated $FS_{\pi\ping} = 1.014$ is far below the minimum code-mandated threshold of $FS \ge 3.0 – 4.0$. Sand boiling and rapid backward piping erosion are imminent at the downstream toe. A weighted toe filter berm or inverted granular filter must be installed immediately to extend seepage path length and suppress exit gradients.

6.2 Example 2: Terzaghi Filter Granular Compatibility Verification

Problem Statement:

Verify whether a candidate crushed sand filter protects an impervious silty clay core:

• Base Core Soil Particle Distribution: $d_{15} = 0.004 \text{ mm}$, $d_{50} = 0.035 \text{ mm}$, $d_{85} = 0.080 \text{ mm}$

• Candidate Granular Filter Distribution: $D_{15} = 0.320 \text{ mm}$, $D_{50} = 1.10 \text{ mm}$, $D_{85} = 3.50 \text{ mm}$, $D_{60} = 1.40 \text{ mm}$, $D_{10} = 0.25 \text{ mm}$.

Step 1: Check Retention (Piping Prevention) Criterion

The allowable upper bound:
$$\text{Retention Ratio} = \frac{D_{15\text{ (filter)}}}{d_{85\text{ (base)}}} = \frac{0.320 \text{ mm}}{0.080 \text{ mm}} = 4.00$$

Because $4.00 \le 5.00$, the candidate filter satisfies the retention criterion and prevents base particle migration.

Step 2: Check Permeability (Drainage) Criterion

$$\text{Permeability Ratio} = \frac{D_{15\text{ (filter)}}}{d_{15\text{ (base)}}} = \frac{0.320 \text{ mm}}{0.004 \text{ mm}} = 80.00$$

Because $80.00 \ge 5.00$, the filter possesses adequate hydraulic conductivity to drain seepage without pore pressure buildup.

Step 3: Check Coefficient of Uniformity ($C_u$) to Avoid Segregation

$$C_u = \frac{D_{60}}{D_{10}} = \frac{1.40 \text{ mm}}{0.25 \text{ mm}} = 5.60$$

Because $C_u = 5.60 \le 6.00$, the granular filter will not segregate during placement. The material is fully compliant.

7. Dam Safety Standards, ICOLD Guidelines, and USBR Regulations

International engineering bodies govern embankment dam safety through rigorous design codes:

• USBR Design Standards No. 13 (Embankment Dams): Mandates specific drainage blanket details, minimum filter thicknesses ($\ge 1.0 \text{ m}$ for machine-placed layers), and rapid drawdown stability requirements.

• ICOLD Bulletin 164 (Internal Erosion of Existing Dams): Establishes the modern toolbox for assessing backward erosion, contact erosion, and suffusion risks in aging infrastructure.

• FEMA 642 / Federal Guidelines for Dam Safety: Establishes emergency action planning, probable maximum flood (PMF) routing criteria, and risk-informed decision making.

• Eurocode 7 (EN 1997-1) Geotechnical Design: Prescribes limit state verification (GEO and UPL states) for seepage buoyancy and hydraulic heave.

8. Synthesis and Engineering Takeaways on Dam Failure Prevention

Mitigating the causes of earth dam failure requires comprehensive defense-in-depth engineering. Hydraulic overtopping demands adequately sized spillways capable of passing the Probable Maximum Flood without crest submergence.

Internal seepage demands multi-stage granular filters designed under Terzaghi and Sherard retention rules to suppress backward piping. Structural slope stability requires rigorous limit equilibrium evaluations under rapid drawdown and seismic conditions. Robust geotechnical instrumentation—piezometers, seepage measuring flumes, and crest settlement monuments—ensures dam safety throughout operational service life.

References & Standards Cited

1. Fell, R., MacGregor, P., Stapledon, D., Bell, G., & Foster, M. (2015). Geotechnical Engineering of Dams (2nd Edition). London: CRC Press / Balkema.
2. International Commission on Large Dams (ICOLD). (2017). Internal Erosion of Existing Dams, Levees and Dikes, and Their Foundations. ICOLD Bulletin 164, Paris.
3. Sherard, J. L., & Dunnigan, L. P. (1989). Critical Filters for Impervious Soils. Journal of Geotechnical Engineering, ASCE, 115(7), 927-947.
4. Terzaghi, K., Peck, R. B., & Mesri, G. (1996). Soil Mechanics in Engineering Practice (3rd Edition). New York: John Wiley & Sons.
5. U.S. Army Corps of Engineers (USACE). (2004). General Design and Construction Considerations for Earth and Rock-Fill Dams. Engineer Manual EM 1110-2-2300, Washington, D.C.
6. U.S. Bureau of Reclamation (USBR). (2014). Design Standards No. 13: Embankment Dams (Chapter 5: Protective Filters). Denver, CO: USBR.

Frequently Asked Questions (FAQ)

What are the most common causes of earth dam failure?

The three primary causes are hydraulic failures (overtopping and spillway inadequacy, $approx 35%$), seepage failures (backward erosion piping and internal suffusion, $approx 38%$), and structural/slope stability failures (drawdown slips and earthquake liquefaction, $approx 27%$).

How does backward erosion piping cause a dam to breach?

Backward erosion piping initiates when high exit hydraulic gradients at the downstream toe dislodge soil particles. This forms sand boils and erodes an open conduit progressively backward into the core, eventually creating an uncontrolled release path.

What is the critical hydraulic gradient for soil piping?

The critical hydraulic gradient is calculated as $i_{cr} = (G_s - 1) / (1 + e)$, where $G_s$ is specific gravity and $e$ is void ratio. For typical soils with $G_s approx 2.65$ and $e approx 0.65$, $i_{cr} approx 1.0$.

Why does rapid drawdown destabilize the upstream slope of an earth dam?

During rapid drawdown, stabilizing reservoir water weight is removed quickly, while high pore water pressures remain trapped inside the low-permeability upstream shell. This reduces effective stress and shear strength, triggering deep rotational slips.

How do granular filters prevent internal erosion in embankment dams?

Granular filters follow Terzaghi retention criteria ($D_{15text{ filter}} / d_{85text{ base}} le 5$) to physically block moving core soil particles while maintaining high permeability ($D_{15text{ filter}} / d_{15text{ base}} ge 5$) to safely drain seepage without building excess pore pressure.