🎓 Lesson 9
D5
Calibrating Infiltration Parameters from Field Data
Calibrating infiltration parameters means using real field measurements—like how fast water soaks into soil—to fine-tune the numbers in hydrologic models so they accurately predict groundwater recharge and surface runoff.
🎯 Learning Objectives
- ✓ Calculate saturated hydraulic conductivity (Ksat) from double-ring infiltrometer data using the steady-state approximation
- ✓ Analyze discrepancies between observed and simulated infiltration curves to identify dominant parameter sensitivity
- ✓ Apply the Green-Ampt model to estimate time-to-ponding and cumulative infiltration for a given soil profile and storm event
- ✓ Explain how antecedent moisture conditions and soil layering affect parameter identifiability during calibration
📖 Why This Matters
In mining, inaccurate infiltration estimates lead to underdesigned surface water diversions, overestimated seepage into open pits, or failed tailings dam cover systems—costing millions in remediation and risking environmental non-compliance. Real-world examples include the 2014 Mount Polley tailings breach, where oversimplified infiltration assumptions contributed to underestimated pore pressure buildup. Calibrating infiltration parameters bridges the gap between textbook theory and site-specific geology—making it foundational for responsible hydrologic design.
📘 Core Principles
Infiltration is governed by soil matrix properties (porosity, texture, structure), initial moisture content, and driving forces (matric suction gradient, ponding depth). Empirical models like Horton (time-decay exponential) and Green-Ampt (sharp wetting front with capillary drive) simplify physics but require site-specific parameterization. Calibration isn’t curve-fitting—it’s constraint-based inverse modeling: Ksat must align with grain-size-derived pedotransfer functions; initial deficit must reflect recent antecedent rainfall; and effective porosity must be physically bounded by measured bulk density and particle density. Heterogeneity (e.g., fractured bedrock beneath colluvium) demands layered or dual-domain calibration strategies.
📐 Green-Ampt Cumulative Infiltration
The Green-Ampt model predicts cumulative infiltration (F) over time (t) assuming a sharp wetting front and constant saturated hydraulic conductivity. It’s widely used in mine drainage design due to its physical basis and robustness for coarse- to medium-textured soils.
Green-Ampt Cumulative Infiltration (Implicit Form)
F = K_sat · t + ψ_f · Δθ · ln(1 + F / (ψ_f · Δθ))Estimates total water infiltrated (F) into homogeneous soil over time (t), accounting for capillary drive and soil storage capacity.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| F | Cumulative infiltration | m | Total depth of water infiltrated up to time t |
| K_sat | Saturated hydraulic conductivity | m/s | Maximum steady-state flow rate through saturated soil |
| t | Time | s | Duration of infiltration |
| ψ_f | Average matric suction head at wetting front | m | Negative pressure (capillary head) driving infiltration |
| Δθ | Effective porosity (θ_s − θ_i) | m³/m³ | Volumetric water content difference between saturated and initial states |
Typical Ranges:
Sand: 1.0 × 10⁻⁵ to 1.0 × 10⁻³ m/s
Silt loam: 1.0 × 10⁻⁷ to 1.0 × 10⁻⁵ m/s
Clay: 1.0 × 10⁻⁹ to 1.0 × 10⁻⁷ m/s
💡 Worked Example
Problem: A double-ring infiltrometer test on glacial till at a proposed leach pad site yields steady-state infiltration rate = 0.8 mm/hr after 90 minutes. Measured saturated moisture content θs = 0.38 cm³/cm³, initial moisture content θi = 0.12 cm³/cm³, and average matric suction head ψf = −25 cm. Estimate cumulative infiltration at t = 120 min.
1.
Step 1: Compute effective porosity (Δθ) = θs − θi = 0.38 − 0.12 = 0.26 cm³/cm³
2.
Step 2: Convert ψf to consistent units: −25 cm = −0.25 m; convert Ksat = 0.8 mm/hr = 2.22 × 10⁻⁷ m/s
3.
Step 3: Apply Green-Ampt implicit solution: F = Ksat·t + ψf·Δθ·ln(1 + F/(ψf·Δθ)). Solve iteratively (or use approximate explicit form F ≈ Ksat·t + ψf·Δθ·ln(t) for t > 1 hr). Using Newton-Raphson with t = 7200 s → F ≈ 0.011 m (11 mm).
Answer:
The estimated cumulative infiltration at 120 minutes is 11 mm, which falls within the typical range of 5–25 mm for compacted glacial till under 2-hour ponding.
🏗️ Real-World Application
At the Cadia East copper-gold mine (NSW, Australia), engineers calibrated infiltration parameters for a 30-ha waste rock dump cover system using 12 double-ring infiltrometer tests across three soil layers (clay cap, sandy loam transition, gravelly base). By fitting Green-Ampt parameters to 48 hours of post-rainfall infiltration monitoring, they reduced predicted long-term seepage volume uncertainty from ±65% to ±12%, enabling redesign of the collection trench capacity and avoiding $2.3M in over-engineering.
🔧 Interactive Calculator
🔧 Open Drainage & Hydrologic Design Calculator📋 Case Connection
📋 Industrial Park Stormwater Master Plan in Indianapolis, IN
Achieving 90% annual phosphorus load reduction while accommodating truck traffic and hazardous material spill containmen...
📋 Mountainous Ski Resort Stormwater Retrofit in Aspen, CO
Preventing erosion and sediment delivery to Roaring Fork River (Class I trout stream) under high-intensity snowmelt puls...