🎓 Lesson 9
D5
Saturation Indices and Mineral Stability Fields: What They Tell You (and Don’t Tell You)
Saturation indices tell you whether a mineral is likely to dissolve or form in water, based on how 'full' the water is with its dissolved ingredients.
🎯 Learning Objectives
- ✓ Calculate saturation indices for key minerals (e.g., gypsum, ferrihydrite, calcite) using aqueous speciation output from PHREEQC or similar batch models
- ✓ Interpret mineral stability fields on Pourbaix or activity–activity diagrams to predict dominant solid phases under varying pH and redox conditions
- ✓ Explain limitations of saturation indices—including kinetic inhibition, metastable phases, and absence of kinetic or microbial controls—in predicting actual mineral formation or dissolution in mine waste environments
- ✓ Apply saturation index trends across a geochemical time-series dataset to diagnose evolving acid rock drainage (ARD) or neutral mine drainage (NMD) behavior
📖 Why This Matters
In mine waste management, predicting whether iron hydroxides will precipitate—or whether sulfates will crystallize and heave containment liners—is not academic: it determines whether a tailings facility remains stable for decades or fails catastrophically. Saturation indices are your first-line diagnostic tool—they flag *thermodynamic potential* for mineral reactions before field evidence appears. But misreading them as guarantees of reaction leads to costly oversights: e.g., assuming jarosite precipitation means ARD is controlled, when in reality it may be kinetically stalled or bypassed by faster pathways. This lesson bridges the gap between model output and engineering judgment.
📘 Core Principles
Saturation indices derive from equilibrium thermodynamics: they compare the *actual* ion activity product (Q) in solution to the *theoretical* equilibrium constant (K) for a mineral’s dissolution reaction. However, natural systems are rarely at equilibrium—especially in dynamic mine waste environments where flow rates, microbial activity, surface area, and aging alter reaction pathways. Mineral stability fields—graphical representations of SI = 0 contours across axes like pH–Eh or log[Fe²⁺]–pH—show *where* minerals *could* be stable *if* equilibrium were achieved and all relevant phases were considered. Critically, stability fields assume idealized systems: no solid-solution effects, no amorphous or nanocrystalline phases, and no competing nucleation barriers. Real-world interpretation requires overlaying field observations (e.g., SEM-EDS mineral ID, porewater trends) to test whether predicted phases actually occur—and if not, why.
📐 Saturation Index Calculation
The saturation index quantifies departure from equilibrium for any mineral reaction. It is computed directly from aqueous speciation modeling software (e.g., PHREEQC) but can be verified manually using activity-based calculations. A positive SI does *not* guarantee precipitation—it only signals thermodynamic favorability; nucleation kinetics, inhibitor ions (e.g., Si, PO₄³⁻), and lack of seed surfaces often prevent it.
Saturation Index (SI)
SI = log₁₀(Q / K)Quantifies thermodynamic tendency for a mineral to precipitate (SI > 0), dissolve (SI < 0), or remain at equilibrium (SI = 0).
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| Q | Ion activity product | dimensionless | Product of aqueous species activities, each raised to stoichiometric coefficient in dissolution reaction. |
| K | Solubility product constant | dimensionless | Equilibrium constant for mineral dissolution at specified temperature and pressure; tabulated in thermodynamic databases. |
Typical Ranges:
Acid rock drainage systems: -5.0 to +3.0
Neutral mine drainage systems: -4.0 to +2.5
💡 Worked Example
Problem: Given PHREEQC output for a mine seepage sample (pH = 2.8, T = 25°C): [Fe³⁺]ₐcₜ = 1.2×10⁻⁴ mol/kg, [SO₄²⁻]ₐcₜ = 3.7×10⁻² mol/kg, [H⁺]ₐcₜ = 1.6×10⁻³ mol/kg. Calculate SI for jarosite (K = 10⁻²²·⁵⁶, using the simplified dissolution reaction: KFe₃(SO₄)₂(OH)₆ ⇌ K⁺ + Fe³⁺ + 2SO₄²⁻ + 6H⁺ + 3H₂O). Assume unit activity for solids and water; use γ ≈ 1 for dilute estimate.
1.
Step 1: Write Q = (a_K⁺)(a_Fe³⁺)(a_SO₄²⁻)²(a_H⁺)⁶. With γ ≈ 1 and negligible K⁺ (assume [K⁺] = 10⁻⁶ M, below detection), approximate Q ≈ (1.2×10⁻⁴)(3.7×10⁻²)²(1.6×10⁻³)⁶
2.
Step 2: Compute: (3.7×10⁻²)² = 1.37×10⁻³; (1.6×10⁻³)⁶ = (1.6)⁶ × 10⁻¹⁸ ≈ 16.8 × 10⁻¹⁸ = 1.68×10⁻¹⁷; then Q ≈ (1.2×10⁻⁴)(1.37×10⁻³)(1.68×10⁻¹⁷) ≈ 2.76×10⁻²⁴
3.
Step 3: SI = log₁₀(Q/K) = log₁₀(2.76×10⁻²⁴ / 10⁻²²·⁵⁶) = log₁₀(2.76×10⁻²⁴) − (−22.56) = (−23.56) + 22.56 = −1.00
Answer:
SI = −1.00, indicating jarosite is undersaturated—consistent with absence of jarosite crusts in this low-pH, high-sulfate seepage. This supports ARD persistence and informs selection of alkaline amendments rather than relying on passive jarosite attenuation.
🏗️ Real-World Application
At the Mount Polley tailings storage facility (British Columbia), pre-closure PHREEQC batch modeling predicted SI > +1.5 for schwertmannite across porewater samples from the upper oxidized zone. Field validation via XRD and synchrotron μ-XANES confirmed schwertmannite dominance—validating model predictions and supporting long-term ARD forecasts. In contrast, SI for goethite was +0.8, yet goethite abundance was low: kinetic inhibition due to rapid sulfate adsorption on Fe-oxyhydroxide surfaces suppressed transformation. This case underscores that SI > 0 is necessary—but not sufficient—for observed mineral occurrence; engineering decisions (e.g., cover design, water routing) must integrate both SI trends *and* mineralogical monitoring.
🔧 Interactive Calculator
🔧 Open Mine Waste Characterization & Geochemical Modeling Calculator📋 Case Connection
📋 Gold Tailings Geochemical Stabilization at Granny Smith Mine (WA)
Arsenic-rich tailings (up to 120 mg/kg As) exhibiting elevated As leaching under oxidizing conditions
📋 Limestone Mine Neutral Drainage Management at Mount Read Complex (Tasmania)
Historic waste dumps containing carbonate-hosted Pb-Zn mineralization generating neutral metal leachate (Zn >15 mg/L, Cd...
📋 Iron Ore Mine Waste Rock Long-Term Stability at Brockman 4 (Pilbara)
Massive hematite-goethite waste rock (low sulfide but high Mn/Al) showing delayed acidity and Al leaching post-construct...
📋 Coal Mine Spoil Geochemical Capping at Hunter Valley Reclamation Project
Spoil with pyritic shale interbeds generating ARD despite initial alkaline overburden; inconsistent capping led to local...