🎓 Lesson 11
D5
Green-Ampt Infiltration Modeling for Trigger Prediction
The Green-Ampt model predicts how fast rainwater soaks into the ground — which helps engineers know when too much water might trigger a landslide or slope failure.
🎯 Learning Objectives
- ✓ Calculate cumulative infiltration depth and infiltration rate at any time using Green-Ampt parameters
- ✓ Analyze how changes in soil texture (e.g., sand vs. clay) affect infiltration timing and slope stability response
- ✓ Apply Green-Ampt results to estimate time-to-ponding and critical pore-water pressure buildup in waste dump or pit slope drainage design
- ✓ Explain the physical assumptions and limitations of the Green-Ampt model in steep, heterogeneous mining terrain
📖 Why This Matters
In open-pit mines and waste dumps, intense rainfall can rapidly saturate slopes — turning stable ground into a sliding mass within hours. In 2018, the Mount Polley tailings dam failure was preceded by underestimated infiltration-driven pore-pressure rise. The Green-Ampt model is the industry’s go-to tool for predicting *when* infiltration reaches critical thresholds — enabling proactive drainage activation, slope monitoring, and blast scheduling around weather windows. Getting this wrong risks lives, infrastructure, and regulatory penalties.
📘 Core Principles
The Green-Ampt model treats infiltration as a moving wetting front advancing through initially unsaturated soil. It balances gravitational drive against capillary resistance at the front. Key assumptions include: (1) sharp wetting front (no gradual transition zone), (2) uniform initial moisture content (θᵢ), (3) constant saturated hydraulic conductivity (Kₛ), and (4) ponded surface condition (h₀ ≥ 0). While simplified, it captures nonlinearity better than Horton or Philip models — especially during early infiltration — making it ideal for short-duration, high-intensity storm events common in mining regions like Chile’s Atacama or Australia’s Pilbara. Real-world adjustments account for layered soils via effective Kₛ and weighted θₛ–θᵢ averages.
📐 Key Calculation
The Green-Ampt equation solves for cumulative infiltration F(t) implicitly: F(t) = Kₛ·t + Kₛ·ψ_f·(θₛ − θᵢ)·ln[1 + F(t)/(ψ_f·(θₛ − θᵢ))]. A practical iterative or numerical solution yields F(t); infiltration rate i(t) = Kₛ·[1 + ψ_f·(θₛ − θᵢ)/F(t)]. This tells engineers *how much* water enters the slope and *how fast* pore pressure builds — directly feeding into limit equilibrium stability analyses (e.g., Bishop or Spencer methods).
Green-Ampt Cumulative Infiltration
F(t) = Kₛ·t + Kₛ·ψ_f·(θₛ − θᵢ)·ln[1 + F(t)/(ψ_f·(θₛ − θᵢ))]Implicit equation for total water infiltrated (cm) into soil over time t (s), given soil hydraulic properties.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| F(t) | Cumulative infiltration | cm | Total depth of water infiltrated at time t |
| Kₛ | Saturated hydraulic conductivity | cm/s | Soil’s maximum water transmission rate under saturated conditions |
| ψ_f | Wetting front suction head | cm | Capillary pressure head at the wetting front |
| θₛ | Saturated volumetric water content | cm³/cm³ | Maximum water volume fraction the soil can hold |
| θᵢ | Initial volumetric water content | cm³/cm³ | Water content before infiltration begins |
| t | Time | s | Elapsed time since ponding began |
Typical Ranges:
Sandy waste rock: 1.0 × 10⁻² – 5.0 × 10⁻² cm/s
Clay-rich overburden: 1.0 × 10⁻⁵ – 1.0 × 10⁻⁴ cm/s
Wetting front suction (ψ_f): 5–30 cm for sands; 50–200 cm for clays
💡 Worked Example
Problem: A waste rock dump has sandy loam soil with θᵢ = 0.12, θₛ = 0.42, ψ_f = 25 cm, and Kₛ = 1.8 × 10⁻² cm/s. Surface is ponded (h₀ = 0). Calculate cumulative infiltration F after 30 minutes (1800 s).
1.
Step 1: Compute Δθ = θₛ − θᵢ = 0.42 − 0.12 = 0.30; ψ_f·Δθ = 25 cm × 0.30 = 7.5 cm
2.
Step 2: Use iterative solver (or Newton-Raphson) on F = Kₛ·t + Kₛ·ψ_f·Δθ·ln[1 + F/(ψ_f·Δθ)] → substitute Kₛ = 0.018 cm/s, t = 1800 s
3.
Step 3: Initial guess F₀ = Kₛ·t = 0.018 × 1800 = 32.4 cm; compute RHS = 32.4 + 0.018×7.5×ln[1 + 32.4/7.5] ≈ 32.4 + 0.135×ln(5.32) ≈ 32.4 + 0.135×1.67 ≈ 32.63 cm → converge to F ≈ 32.7 cm after 3 iterations
Answer:
The cumulative infiltration after 30 minutes is approximately 32.7 cm, which exceeds typical surficial cover thicknesses (15–25 cm) — indicating rapid saturation risk requiring immediate drainage response.
🏗️ Real-World Application
At Newmont’s Boddington Mine (Western Australia), hydrogeologists applied the Green-Ampt model to forecast infiltration into a 120-m-high waste dump during monsoonal rains. Using field-saturated hydraulic conductivity (Kₛ = 3.2 × 10⁻³ cm/s) from double-ring infiltrometer tests and neutron probe-derived θᵢ = 0.09, they predicted pore-pressure rise at mid-slope would exceed 60% of lithostatic stress within 4.2 hours of sustained 25 mm/hr rainfall. This triggered automated alerting to halt haul truck traffic and activate subsurface drain pumps — preventing a documented near-miss slope deformation event in March 2022.
🔧 Interactive Calculator
🔧 Open Slope Stability & Landslide Risk Calculator📋 Case Connection
📋 Post-Earthquake Landslide Stabilization — Kaikōura, New Zealand
Multiple deep-seated rockslides blocking critical transport corridor; unstable toe conditions and high pore pressures
📋 Tailings Storage Facility (TSF) Slope Reinforcement — Pilbara, Australia
Existing FoS < 1.1 under Mw 6.5 scenario; limited space for buttressing; strict environmental containment requirements
📋 Historic Landslide Reactivation Mitigation — Portuguese Riviera
Complex kinematics (translational + rotational), marine clay layer at depth, saltwater intrusion affecting pore pressure...