Friday, 24 April, 2026
Case Study: Cr(VI) Plume
This case study illustrates the coupling between redox speciation and transport within a groundwater chromium plume. It addresses a practical question: an industrial site has a history of chromium release, some of which has migrated into the aquifer—how can we predict its downgradient extension and identify areas where natural attenuation through redox reduction is effective?
Context
The site involves a former surface treatment facility that ceased operations fifteen years ago. Several monitoring wells show total chromium concentrations between 200 and 1,200 µg/L in the shallow aquifer, with a Cr(VI) to total chromium ratio that varies from point to point. The project owner wishes to evaluate whether natural attenuation is sufficient to consider the site stabilized, or if active remediation is required.
The typical challenges in this type of situation boil down to three questions: where is the plume front currently located, how will it evolve over the next ten years, and what fraction of the migrated chromium has already been reduced to immobile Cr(III)?
Scene Setup
The selected footprint covers 800 meters along the flow direction and 400 meters transversely, centered on the suspected source zone. The depth ranges from the surface to 15 meters, including the shallow aquifer and the clay layer separating it from the deep aquifer. Seven observation boreholes provide sampling points, with three in the source zone and four downgradient.
| Scene Element | Configuration |
|---|---|
| Horizontal Footprint | 800 × 400 m, aligned with the flow axis |
| Depth | 15 m, referenced to ground level |
| Calculation Grid | 80 × 40 × 30 = 96,000 cells (intermediate resolution) |
| Surface Layer (0 to 2 m) | Sandy loam, porosity 0.20, K = 5e-5 m/s |
| Shallow Aquifer (2 to 12 m) | Coarse sand, porosity 0.25, K = 5e-4 m/s |
| Clay Aquitard (12 to 15 m) | Clay, porosity 0.04, K = 1e-9 m/s |
| Sampling Points | 7 boreholes, all equipped with multi-level piezometers |
Engine Parameters
Chromium is simulated in its two primary forms: mobile Cr(VI) with a low Kd in sand (0.5 L/kg) and practically immobile Cr(III) (Kd 200 L/kg). The transformation from Cr(VI) to Cr(III) is modeled as a kinetic reaction where the rate depends on local Eh and the presence of electron donors (ferrous iron, organic matter). Eh measurements across the seven boreholes range from +250 mV to -50 mV, suggesting significant redox heterogeneity.
| Parameter | Value | Source |
|---|---|---|
| Source | 50 × 30 m zone, depth 2 to 6 m, estimated historical flux | Facility operational records |
| Hydraulic Gradient | 0.004 m/m, resolved by 7-point piezometry | Borehole measurements |
| Mean Effective Velocity in Sand | Approx. 25 m/year | Darcy's Law with 0.25 porosity |
| Kd Cr(VI) | 0.5 L/kg | Literature (low-adsorption sandy aquifer) |
| Kd Cr(III) | 200 L/kg | Literature (hydroxide precipitation) |
| Cr(VI) to Cr(III) Reduction Rate | Variable 0.001 to 0.05 /day | Function of local Eh, calibrated to Cr(VI)/Total Cr ratios |
| Simulated Duration | 20 years (15 years historical + 5 years prospective) | Covers period since cessation of activity |
The critical calibration point is the Cr(VI) to Cr(III) reduction rate. It varies by several orders of magnitude depending on Eh, organic matter, and microbial biomass. Calibrating it against the measured Cr(VI)/Total Cr ratios at the boreholes transforms this "fuzzy" parameter into one constrained by field data.
Results
After 20 years of simulation, the total chromium concentration field exhibits three distinct zones. The source zone still accumulates a significant residual mass, largely in the form of precipitated Cr(III). An intermediate zone spanning a few hundred meters downgradient is dominated by Cr(VI) in transit with signs of partial reduction. A distant downgradient zone shows only traces of chromium at concentrations near analytical detection limits.
A vertical cross-section along the flow axis reveals clear vertical stratification: chromium flows preferentially through the most permeable sandy layers and accumulates at the bottom of the aquifer, just above the clay aquitard. The δ⁵³Cr isotopic signature follows this structure: the most enriched values (indicating advanced reduction) are found at the bottom of the aquifer in the intermediate zone, where the lowest Eh was measured.
| Zone | Simulated Total Cr Concentration | Cr(VI) Fraction | Simulated δ⁵³Cr |
|---|---|---|---|
| Source (0 to 100 m) | 500 to 1500 µg/L | 20 to 40 % | +0.5 to +1.2 ‰ |
| Intermediate (100 to 400 m) | 100 to 500 µg/L | 40 to 70 % | +1.5 to +3.5 ‰ |
| Downgradient Front (400 to 700 m) | 5 to 50 µg/L | 60 to 90 % | +0.8 to +2.5 ‰ |
Interpretation
The model generally reproduces the measurements at the seven boreholes: concentrations and Cr(VI)/Total Cr ratios are accurate in the source and near-downgradient zones, while predicted concentrations in the far-downgradient zone fall below detection limits, consistent with field data. The simulated δ⁵³Cr isotopic signature at the three measured points falls within the confidence interval of the measurements, validating the reduction parameter settings.
A further five-year projection shows an overall stabilization of the plume. While the source zone continues to provide Cr(VI) through slow desorption, the mobilized mass is largely offset by reduction to Cr(III) in the intermediate zone. The downgradient front is no longer advancing significantly, and residual concentrations continue to decrease slowly.
This interpretation must be qualified by several identified uncertainties. First, the reduction rate is calibrated to current conditions; if redox conditions change (e.g., dewatering or changes in upstream land use), the dynamics could shift. Second, the source inventory remains poorly quantified, and slower-than-expected desorption could extend plume activity beyond current predictions. Finally, the regular grid representation does not capture potential preferential pathways linked to local heterogeneities.
The interpretation of "sufficient natural attenuation" suggested by the simulation remains conditional on the parameter assumptions. Active remediation (such as ZVI barriers or reductant injection) can accelerate reduction by a factor of 5 to 10 and remains an option to consider if the predicted stability is deemed too slow relative to site objectives.
Extensions
This case can be extended in several useful ways to deepen the understanding of the site or to test potential actions.
- "No Reduction" variant: Disable the Cr(VI) to Cr(III) reaction to see what the plume would look like without natural attenuation. This comparison quantifies the actual impact of reduction on plume extension.
- "Reactive Barrier" variant: Add a zone with high reducing capacity perpendicular to the flow to simulate a ZVI (Zero Valent Iron) installation. The comparison guides sizing (length, thickness, position).
- "Dry Climate" variant: Decrease meteorological recharge by 30%, which slows down the flow. This comparison shows whether the plume is controlled by transport (beneficial acceleration) or by reduction kinetics (neutral or harmful deceleration).
- "Persistent Source" variant: Extend the duration of active source input to test the hypothesis of a longer-than-expected release.
These variants can be saved as distinct scenarios in the Hypothesis Comparison module and compared directly within the same scene.
Learn More
- Speciation and Propagation: Details of the coupling used in this case.
- Redox Speciation: Cr(VI) vs. Cr(III) in general geochemistry.
- Multi-Hypothesis: Compare the variants discussed in the extensions.
- Interpreting Results: General guide to simulation outputs.