Christoph Wanner

Assessing the Cr(VI) reduction efficiency of a permeable reactive barrier using Cr isotope measurements and 2D reactive transport modeling

In Thun, Switzerland, a permeable reactive barrier (PRB) for Cr(VI) reduction by gray cast iron was installed in May 2008. The PRB is composed of a double array of vertical piles containing iron shavings and gravel. The aquifer in Thun is almost saturated with dissolved oxygen and the groundwater flow velocities are ca. 10-15m/day. Two years after PRB installation Cr(VI) concentrations still permanently exceed the Swiss threshold value for contaminated sites downstream of the barrier at selected localities. Groundwater δ(53/52)Cr(SRM979) measurements were used to track Cr(VI) reduction induced by the PRB. δ(53/52)Cr(SRM979) values of two samples downstream of the PRB showed a clear fractionation towards more positive values compared to four samples from the hotspot, which is clear evidence of Cr(VI) reduction induced by the PRB. Another downstream sample did not show a shift to more positive δ(53/52)Cr(SRM979) values. Because this latter location correlates with the highest downstream Cr(VI) concentration it is proposed that a part of the Cr(VI) plume is bypassing the barrier. Using a Rayleigh fractionation model a minimum present-day overall Cr(VI) reduction efficiency of ca. 15% was estimated. A series of 2D model simulations, including the fractionation of Cr isotopes, confirm that only a PRB bypass of parts of the Cr(VI) plume can lead to the observed values. Additionally, the simulations revealed that the proposed bypass occurs due to an insufficient permeability of the individual PRB piles. It is concluded that with this type of PRB a complete and long-lasting Cr(VI) reduction is extremely difficult to achieve for Cr(VI) contaminations located in nearly oxygen and calcium carbonate saturated aquifer in a regime of high groundwater velocities. Additional remediation action would limit the environmental impact and allow to reach target concentrations.

Benchmarking the simulation of Cr isotope fractionation

A benchmark problem set consisting of four problem levels was developed for the simulation of Cr isotope fractionation in 1D and 2D domains. The benchmark is based on a recent field study where Cr(VI) reduction and accompanying Cr isotope fractionation occurs abiotically by an aqueous reaction with dissolved Fe 2+ (Wanner et al., 2012., Appl. Geochem., 27, 644–662). The problem set includes simulation of the major processes affecting the Cr isotopic composition such as the dissolution of various Cr(VI) bearing minerals, fractionation during abiotic aqueous Cr(VI) reduction, and non-fractionating precipitation of Cr(III) as sparingly soluble Cr-hydroxide.Accuracy of the presented solutions was ensured by running the problems with four well-established reactive transport modeling codes: TOUGHREACT, MIN3P, CRUNCHFLOW, and FLOTRAN. Results were also compared with an analytical Rayleigh-type fractionation model. An additional constraint on the correctness of the results was obtained by comparing output from the problem levels simulating Cr isotope fractionation with the corresponding ones only simulating bulk concentrations. For all problem levels, model to model comparisons showed excellent agreement, suggesting that for the tested geochemical processes any code is capable of accurately simulating the fate of individual Cr isotopes.

A benchmark for microbially mediated chromium reduction under denitrifying conditions in a biostimulation column experiment

Bioremediation efforts in aquifers contaminated with redox-sensitive contaminants often rely on in situ reductive immobilization. The bioremediation treatment usually involves injection of organic carbon into the subsurface (e.g., acetate) to stimulate the growth of indigenous bacteria that mediate the relevant redox processes that immobilize the target contaminant. Batch and flow-through column experimental studies are conducted to elucidate reaction networks associated with specific electron acceptor pathways and/or specific bacterial isolates. The proposed benchmark involves the simulation of microbially mediated chromium reduction under denitrifying conditions in biostimulated batch and flow-through column experiments. Simulated reactive processes include multicomponent aqueous complexation, kinetically controlled mineral precipitation and dissolution, biologically mediated reactions, and biomass growth and decay. The focus of the benchmark problem set is on the simulation of microbially mediated redox reactions with the explicit inclusion of the microbial community dynamics and the impacts on reaction rates. Rate expressions for microbially mediated redox reactions include kinetic limitations (Monod and inhibition terms) as well as thermodynamic limitations. Both catabolic (energy) and anabolic pathways (biomass growth) are considered in the microbially mediated reactions. Microbial biomass is assumed to be bound to the sediment (non-planktonic). Any reactive transport model used to reproduce results of this benchmark problem must be capable of simulating multicomponent aqueous complexation, kinetically controlled mineral precipitation and dissolution and kinetically controlled aqueous reactions. Though convenient, it is not necessary to allow for specific stoichiometric relationships for catabolic and anabolic pathways; only the overall reaction stoichiometry is used. Rate expressions for microbially mediated reaction must include a rate constant, the biomass concentration, and a number of Monod and inhibition terms. To ensure that the results presented in this paper were the correct solutions to the problems posed, the general-purpose reactive transport codes CrunchFlow, PHT3D, ToughReact, and MIN3P were used to perform the simulations. In general, results obtained with all codes show excellent agreement.

A reactive transport benchmark on modeling biogenic uraninite re-oxidation by Fe(III)-(hydr)oxides

A reactive transport benchmark on uranium (U) bioreduction and concomitant reoxidation has been developed based on the multicomponent biogeochemical reaction network presented by Spycher et al. (Geochim Cosmochim Acta 75:4426–4440, 2011). The benchmark problem consists of a model inter-comparison starting with the numerical simulations of the original batch experiments of Sani et al. (Geochim Cosmochim Acta 68:2639–2648, 2004). The batch model is then extended to 1D and 2D reactive transport models, designed to evaluate the model results for the key biogeochemical reaction processes and their coupling with physical transport. Simulations are performed with four different reactive transport simulators: PHREEQC, PHT3D, MIN3P, and TOUGHREACT. All of the simulators are able to capture the complex biogeochemical reaction kinetics and the coupling between transport and kinetic reaction network successfully in the same manner. For the dispersion-free variant of the problem, a 1D-reference solution was obtained by PHREEQC, which is not affected by numerical dispersion. PHT3D using the sequential non-iterative approach (SNIA) with an explicit TVD scheme and MIN3P using the global implicit method (GIM) with an implicit van Leer flux limiter provided the closest approximation to the PHREEQC results. Since the spatial weighting schemes for the advection term and numerical dispersion played an important role for the accuracy of the results, the simulators were further compared using different solution schemes. When all codes used the same spatial weighting scheme with finite-difference approximation, the simulation results agreed very well among all four codes. The model intercomparison for the 2D-case demonstrated a high level of sensitivity to the mixing of different waters at the dispersive front. Therefore this benchmark problem is well-suited to assess code performance for mixing-controlled reactive transport models in conjunction with complex reaction kinetics.

Numerical modeling of cold magmatic CO2 flux measurements for the exploration of hidden geothermal systems

The most accepted conceptual model to explain surface degassing of cold magmatic CO2 in volcanic-geothermal systems involves the presence of a gas reservoir. In this study, numerical simulations using the TOUGH2-ECO2N V2.0 package are performed to get quantitative insights into how cold CO2 soil flux measurements are related to reservoir and fluid properties. Although the modeling is based on flux data measured at a specific geothermal site, the Acoculco caldera (Mexico), some general insights have been gained. Both the CO2 fluxes at the surface and the depth at which CO2 exsolves are highly sensitive to the dissolved CO2 content of the deep fluid. If CO2 mainly exsolves above the reservoir within a fracture zone, the surface CO2 fluxes are not sensitive to the reservoir size but depend on the CO2 dissolved content and the rock permeability. For gas exsolution below the top of the reservoir, surface CO2 fluxes also depend on the gas saturation of the deep fluid as well as the reservoir size. The absence of thermal anomalies at the surface is mainly a consequence of the low enthalpy of CO2. The heat carried by CO2 is efficiently cooled down by heat conduction and to a certain extent by isoenthalpic volume expansion depending on the temperature gradient. Thermal anomalies occur at higher CO2 fluxes (>37,000 g m−2 d−1) when the heat flux of the rising CO2 is not balanced anymore. Finally, specific results are obtained for the Acoculco area (reservoir depth, CO2 dissolved content, and gas saturation state).