Ilka Wallis

Process-Based Reactive Transport Model To Quantify Arsenic Mobility during Aquifer Storage and Recovery of Potable Water

Aquifer storage and recovery (ASR) is an aquifer recharge technique in which water is injected in an aquifer during periods of surplus and withdrawn from the same well during periods of deficit. It is a critical component of the long-term water supply plan in various regions, including Florida, USA. Here, the viability of ASR as a safe and cost-effective water resource is currently being tested at a number of sites due to elevated arsenic concentrations detected during groundwater recovery. In this study, we developed a process-based reactive transport model of the coupled physical and geochemical mechanisms controlling the fate of arsenic during ASR. We analyzed multicycle hydrochemical data from a well-documented affected southwest Floridan site and evaluated a conceptual/numerical model in which (i) arsenic is initially released during pyrite oxidation triggered by the injection of oxygenated water (ii) then largely complexes to neo-formed hydrous ferric oxides before (iii) being remobilized during recovery as a result of both dissolution of hydrous ferric oxides and displacement from sorption sites by competing anions.

Evaluation of Conceptual and Numerical Models for Arsenic Mobilization and Attenuation during Managed Aquifer Recharge

Managed Aquifer Recharge (MAR) is promoted as an attractive technique to meet growing water demands. An impediment to MAR applications, where oxygenated water is recharged into anoxic aquifers, is the potential mobilization of trace metals (e.g., arsenic). While conceptual models for arsenic transport under such circumstances exist, they are generally not rigorously evaluated through numerical modeling, especially at field-scale. In this work, geochemical data from an injection experiment in The Netherlands, where the introduction of oxygenated water into an anoxic aquifer mobilized arsenic, was used to develop and evaluate conceptual and numerical models of arsenic release and attenuation under field-scale conditions. Initially, a groundwater flow and nonreactive transport model was developed. Subsequent reactive transport simulations focused on the description of the temporal and spatial evolution of the redox zonation. The calibrated model was then used to study and quantify the transport of arsenic. In the model that best reproduced field observations, the fate of arsenic was simulated by (i) release via codissolution of arsenopyrite, stoichiometrically linked to pyrite oxidation, (ii) kinetically controlled oxidation of dissolved As(III) to As(V), and (iii) As adsorption via surface complexation on neo-precipitated iron oxides.

Simulating MODFLOW-based reactive transport under radially symmetric flow conditions.

Radially symmetric flow and solute transport around point sources and sinks is an important specialized topic of groundwater hydraulics. Analysis of radial flow fields is routinely used to determine heads and flows in the vicinity of point sources or sinks. Increasingly, studies also consider solute transport, biogeochemical processes, and thermal changes that occur in the vicinity of point sources/sinks. Commonly, the analysis of hydraulic processes involves numerical or (semi-) analytical modeling methods. For the description of solute transport, analytical solutions are only available for the most basic transport phenomena. Solving advanced transport problems numerically is often associated with a significant computational burden. However, where axis-symmetry applies, computational cost can be decreased substantially in comparison with full three-dimensional (3D) solutions. In this study, we explore several techniques of simulating conservative and reactive transport within radial flow fields using MODFLOW as the flow simulator, based on its widespread use and ability to be coupled with multiple solute and reactive transport codes of different complexity. The selected transport simulators are MT3DMS and PHT3D. Computational efficiency and accuracy of the approaches are evaluated through comparisons with full 2D/3D model simulations, analytical solutions, and benchmark problems. We demonstrate that radial transport models are capable of accurately reproducing a wide variety of conservative and reactive transport problems provided that an adequate spatial discretization and advection scheme is selected. For the investigated test problems, the computational load was substantially reduced, with the improvement varying, depending on the complexity of the considered reaction network.

Southern South Australian groundwater microbe diversity

Groundwater is increasingly used globally for domestic, industrial and agricultural production. While many studies have focused on groundwater as a resource, the diverse ecosystems within are often ignored. Here, we assess 54 Southern South Australian groundwater microbial communities from the populated part of the state to assess their status and dynamics in isolated groundwater systems. We observed a strong site-to-site individuality in groundwater bacterial communities, likely due to the isolated nature of groundwater bodies leading to unique ecosystems. Rank abundance analysis indicates bacterial diversity is maintained even at low abundances and that the distribution fits classical ecological models for strong competition in resource-limited environments. Combined, our data indicates that despite overrepresentation of pollutant-associated bacterial orders in and around the Adelaide metropolitan area, microbial communities remain diverse and show little evidence of converging on a common pollutant-effected community.

Trace metal behavior during in-situ iron removal tests in Leuven, Belgium

Subsurface iron removal (SIR) is an in-situ technique to lower the iron content of extracted groundwater. Through cyclic injection of oxygenated water ferrous iron oxidises and precipitates as iron hydroxide in a zone surrounding the extraction well, enhancing the sorptive capacity of the aquifer. During subsequent pumping phases, groundwater traverses the oxidation zone and ferrous iron sorbs to available and newly formed exchange and sorption sites, thereby retarding the breakthrough of dissolved iron. The process is well-understood in regards to the retardation of iron. Less well understood, however, is the behavior of a number of trace metals and metalloids during SIR operations, foremost arsenic (As). In this study, we analyse major and minor ion and trace metal concentrations from a number of SIR tests in a sand aquifer near Leuven, Belgium. We use reactive transport modelling to evaluate conceptual models of trace metal release and arrest. The test data, underpinned by model results, show that metal release, namely arsenic and barium (Ba), occurs through the oxidation of trace amounts of sulphide minerals during the injection phase. Sorption through cation exchange retards Ba while complexation lowers dissolved As concentrations. Arsenic is mobilized again during the pumping phase through varying phosphate concentrations in the native groundwater, despite available sorption surfaces, while Ba remains adsorbed. Concentrations, however, remain below WHO guideline values for As and Ba (10 μg/l and 0.7 mg/l), respectively. The developed conceptual model of As fate reveals a high propensity for As mobility during SIR due to desorption reactions and delivers an explanation as to why many SIR operations fail to show substantial As removal, despite efficient iron removal. Other monitored trace elements showed no mobilisation, including Zn, Al, Cd, Cr, Cu, F, Hg, Ni, Pb, Sb and Se.

Corrigendum to “Palaeohydrogeology and Transport Parameters Derived from 4He and Cl Profiles in Aquitard Pore Waters in a Large Multilayer Aquifer System, Central Australia”

A study of chloride and 4He profiles through an aquitard that separates the Great Artesian Basin from the underlying Arckaringa Basin in central Australia is presented. The aquitard separates two aquifers with long water residence times, due to low recharge rates in the arid climate. One-dimensional solute transport models were used to determine the advective flux of groundwater across the aquitard as well as establish any major changes in past hydrological conditions recorded by variations of the pore water composition. This in situ study showed that both diffusion and slow downward advection (  mm/yr) control solute transport. Numerical simulations show that an increase in chloride concentration in the upper part of the profile is due to a reduction in recharge in the upper aquifer for at least 3000 years. Groundwater extraction since 2008 has likely increased chloride and 4He concentrations in the lower aquifer by pulling up water from deeper layers; however, there has been insufficient time for upward solute transport into the pore water profile by diffusion against downward advection. The transport model of 4He and chloride provides insight into how the two aquifers interact through the aquitard and how climate change is being recorded in the aquitard profile.