Antimony and arsenic migration in a heterogeneous subsurface at an abandoned antimony smelter under rainfall

While antimony (Sb) and arsenic (As) co-contamination in subsurface soil systems due to the legacy of Sb smelting wastes has been documented, the role of inherent heterogeneity on pollutant migration is largely overlooked. Herein this study investigated Sb and As migration in a slag impacted, vertically stratified subsurface at an abandoned Sb smelter. A 2-dimensional flume was assembled as a lab-scale analogue of the site and subject to rainfall and stop-rain events. Reactive transport modeling was then performed by matching the experimental observations to verify the key factors and processes controlling pollutant migration. Results showed that rainfall caused Sb and As release from the shallow slag layer and promoted their downward movement. Nevertheless, the less permeable deeper layers limited physical flow and transport, which led to Sb and As accumulation at the interface. The re-adsorption of Sb and As onto iron oxides in the deeper, more acidic layers further retarded their migration. Because of the large difference between Sb and As concentrations, Sb re-adsorption was much less effective, which led to higher mobility. Our findings overall highlight the necessity of understanding the degree and impacts of physicochemical heterogeneity for risk exposure assessment and remediation of abandoned Sb smelting sites.

Bioremediation of chlorinated ethenes contaminated groundwater and the reactive transport modeling - A review

Improper disposal of chlorinated ethenes (CEs), a class of widely used solvents in chemical manufacturing and cleaning industries, often leads to severe groundwater contamination. In situ bioremediation of CE-contaminated groundwater has received continuous attention in recent years. The reactive transport simulation is a valuable tool for planning and designing in situ bioremediation systems. This paper presents a detailed and comprehensive review on the main biotransformation pathways of CEs in aquifers, the mathematical modeling of bioremediation processes, and the available software developed for the simulation of reactive transport of CEs over past three decades. The aim of this research is to provide guidance on the selection of appropriate models and software suitable for systems of varying scales, and to discern prevailing research trends while identifying areas worthy of further study. This paper provides a detailed summary of the equations, parameters, and applications of existing biotransformation models from literature studies, highlighting the operation, benefits, and limitations of software available for CEs reactive transport simulations. Lastly, the support of reactive transport simulation programs for the design of full-scale in situ bioremediation systems was elucidated. Further research is needed for incorporating the effects of key subsurface environmental factors on biodegradation processes into models and balancing model complexity with computer data processing power to better support the development and application of reactive transport modeling software.

Temperature-dependent dynamics of electrokinetic conservative and reactive transport in porous media: A model-based analysis

Electrokinetic techniques employ direct current electric fields to enhance the transport of amendments in low permeability porous media and have been demonstrated effective for in situ remediation of both organic contaminants and heavy metals. The application of electric potential gradients give rise to coupled chemical, hydraulic and electric fluxes, which are at the basis of the main transport mechanisms: electromigration and electroosmosis. Previous research has highlighted the significant impacts of charge interactions and fluid composition, including temperature-dependent properties such as electrolyte conductivity and density, on these transport phenomena. However, current models of electrokinetic applications often assume isothermal conditions and overlook the production of heat resulting from Joule heating. This study provides a detailed model-based investigation, systematically exploring the effects of temperature on electrokinetic conservative and reactive transport in porous media. By incorporating temperature-dependent material properties and progressively investigating the impact of temperature on each transport mechanism, we analyze the effects of temperature variations in both 1D and 2D systems. The study reveals how temperature dynamically influences the physical, chemical and electrostatic processes controlling electrokinetic transport. A temperature increase results in a higher speed of amendments delivery by both electromigration and electroosmosis and increases the kinetics of degradation reactions. The simulations also reveal a feedback mechanism in which higher aqueous conductivity results in increased Joule heating, leading to a faster temperature rise and, subsequently, to higher electrolyte conductivity. Finally, we estimate the electric energy requirements of the system at varying temperatures and show how these changes impact the rate of contaminant removal.

Release and mobility of hexavalent chromium in contaminated soil with chemical factory waste: Experiments, Cr isotope analysis and reactive transport modeling

Our study focused on the leaching processes in soil contaminated with hexavalent chromium (Cr(VI)), traced to industrial waste from a disused site and resulting in groundwater contamination. Mineral and geochemical characterization of the soil by means of XRD, SEM-EDS, total digestion and sequential extractions revealed that the main Cr content was from solid waste located in the upper meter of the soil profile. Flow-through and column experiments were carried out to investigate the processes responsible for Cr(VI) release. Cr(VI) mobility along the soil profile was also assessed. Moreover, Cr isotope signatures were used to evaluate a potential Cr(VI) reduction process, which preferably could immobilize toxic Cr(VI) complexes. One-dimensional (1D) numerical simulations reproduced the Cr(VI) release from the flow-through experiment containing the Cr(VI) rich-solid waste and also the Cr(VI) mobility along the column experiment. These results enabled us to interpret quantitatively the processes resulting in Cr(VI) contamination and mobility along the soil profile. Cr(VI) was released from dissolving Cr(VI)-rich phases (e.g., sodium chromate, Cr(VI)-hydrocalumite and Cr(VI)-ettringite) of the solid waste layer. Cr(VI) reduction and Cr(VI) adsorption did not take place along the column. Such accurate characterization of these processes is necessary for the mitigation of Cr(VI) mobility in contaminated soils.

Electrokinetic delivery of permanganate in clay inclusions for targeted contaminant degradation

The use of electrokinetics (EK) has great potential to deliver reactants in impervious porous media, thus overcoming some of the challenges in the remediation of contaminants trapped in low-permeability zones. In this work we experimentally investigate electrokinetic transport in heterogeneous porous media consisting of a sandy matrix with a target clay inclusion. We demonstrate the efficient EK-delivery of permanganate in the target clay zone (transport velocity 0.3-0.5 m day^-1) and its reactivity with Methylene Blue, a positively charged contaminant trapped within the inclusion. The delivery method was optimized using a KH₂PO₄/K₂HPO₄ buffer to attenuate the effect of electrolysis reactions in the electrode chambers, thus mitigating the propagation of pH fronts and preventing the phenomenon of permanganate stalling. The experiments showed that the buffer electrical conductivity greatly impacts the potential gradient in the heterogeneous porous medium with implications on the observed rates of electrokinetic transport (variation up to 40%). The reactive experiments provided direct evidence of the permanganate penetration within the clay and of its capability to degrade the target immobilized contaminant. The experimental results were analyzed using a process-based model, elucidating the governing transport mechanisms and highlighting the effect of different mass transfer processes on conservative and reactive electrokinetic transport.

Leaching and transport of technetium from reducing cementitious waste forms in field lysimeters

Field lysimeters tests examined leaching of technetium-99 (⁹⁹Tc) from two types of cementitious waste forms and found that the presence of blast furnace slag reduced the overall leaching of ⁹⁹Tc from the waste form. The two cementitious waste forms were a slag-grout 45%/45%/10% mixture of fly ash, blast furnace slag, and cement, respectively, referred to as slag-grout or a 55%/45% mixture of cement and fly ash, respectively, referred to here simply as cement. Duplicate sources of each composition were buried in four lysimeters for approximately 10 months to evaluate leaching characteristics under natural meteorological conditions in South Carolina, USA. Effluent samples were collected four times during the experiment, and the distribution of ⁹⁹Tc in the sediment was determined by destructively segmenting the lysimeters at the end of the experiment. The transport of Tc within the lysimeter was simulated by assuming advection, dispersion, and sorption in partially saturated porous media, and by using a shrinking-core type approximation for the release of Tc from the source. The shrinking-core model predicted that the oxidation front created by the oxygenated infiltrating groundwater moved into the cementitious source at a rate of 14 μm/day. As this front moved through the source, Tc(IV) was oxidized to the highly mobile Tc(VII) (as TcO₄^-) species, which then was transported through the sediment primarily via advection due to a small partitioning coefficient (C_solid/C_aq; K_d = 0.14 mL/g). The simulations predicted a cycle of accumulation of Tc in sediment at the source between rainfall events, followed by downward advection due to infiltration during rainfall events. The anomalous upward movement of Tc peak was predicted to be due to upward flux caused by evaporation after the experiment was terminated by capping the lysimeter. These experiments demonstrate that Tc leaching from cementitious waste forms under simulated vadose zone oxidizing conditions can be reasonably approximated by the shrinking core model, and the migration of Tc through the sediment is profoundly influenced by the presence of slag in the grout formulation and hydraulic conditions due to the low sorption affinity of TcO₄^-.

Geochemical audit of a historical tailings storage facility in Japan: Acid mine drainage formation, zinc migration and mitigation strategies

Historical tailings storage facilities (TSFs) are either abandoned or sparsely rehabilitated promoting acid mine drainage (AMD) formation and heavy metal release. To sustainably manage these sites, a geochemical audit coupled with numerical simulation to predict AMD flow paths and heavy metal migration are valuable. In this study, a 40-year-old TSF in Hokkaido, Japan was investigated. Tailings in this historical TSF contain pyrite (FeS₂) while its copper (Cu) and zinc (Zn) contents were 1400-6440 mg/kg and 2800-22,300 mg/kg, respectively. Copper and Zn were also easily released in leaching tests because they are partitioned with the exchangeable phase (29% of Zn; 15% of Cu) and oxidizable fraction (25% of Zn; 33% of Cu). Kinetic modeling results attributed AMD formation to the interactions of pyrite and soluble phases in the tailings with oxygenated groundwater, which is supported by the sequential extraction and leaching results. Calibrations of groundwater/AMD flow and solute transport in the 2D reactive transport model were successfully done using hydraulic heads measured on-site and leaching results, respectively. The model forecasted the quality of AMD to deteriorate with time and AMD formation to continue for 1000 years. It also predicted ~24% AMD flux reduction, including lower Zn release with time when recharge reduction interventions are implemented on-site.

On the interplay between electromigration and electroosmosis during electrokinetic transport in heterogeneous porous media

Electrokinetic techniques represent a valuable approach to enhance solute transport, reactant delivery and contaminant degradation in complex environmental matrices, such as contaminated soil and groundwater, and have a great potential for the remediation of many organic and inorganic pollutants. This study investigates the complex interplay between the key electrokinetic transport mechanisms, electromigration and electroosmosis, in physically heterogeneous porous media and its impact on tracer distribution, reactant mixing and degradation efficiency. We perform experiments in a multidimensional setup, considering different types of heterogeneities, injected tracers and reactants, as well as background electrolyte pore water with different chemical composition and pH. We show that EK transport is significantly affected by the physical heterogeneities, due to the interaction between electrokinetic and hydraulic processes, and by the pore water chemistry that plays a key role on the magnitude and spatial distribution of electroosmotic fluxes. The latter affect the overall transport of charged and non-charged species, including the migration velocity of injected plumes, their spatial patterns, spreading and mixing with the background groundwater, and the extent of degradation and the spatio-temporal evolution of reactive zones in the heterogeneous porous media. Process-based numerical modeling allowed us to interpret the experimental observations and to disentangle the coupled effects of physical, chemical and electrostatic processes in the multidimensional, heterogeneous setups. Besides elucidating the mechanisms controlling electrokinetic transport, the results of this study have also important implications for practical field implementation of EK approaches in intrinsically heterogeneous subsurface systems.

Modeling of soil gas radon as an in situ partitioning tracer for quantifying LNAPL contamination

In the last decades radon (Rn) has been widely proposed as a naturally occurring tracer for non-aqueous phase liquids (NAPL) in the soil. This work examines the feasibility of using soil gas data collected at some distance from the source zone for the application of the Rn deficit technique for the identification and quantification of NAPL contamination. To this end, we used a steady-state 1-D analytical solution that is based on a 3-layer model that allows to simulate the transport and distribution of Rn in the source zone, capillary fringe and overlying unsaturated soil. The analytical solution was first validated against a more detailed numerical model available in the literature. Then, a series of simulations were carried out to evaluate the vertical concentration profiles of Rn in soil gas above the source zone and in background location not impacted by NAPL. Simulation results showed that the parameters that most influence the migration and distribution of Rn in the subsurface are the distance of the soil gas probe from the source zone and, to a lower extent, the type of contamination (e.g. diesel or gasoline) and the soil type. On the basis of these results, we developed some easy-to-use nomographs to estimate the residual NAPL phase based on the observed radon deficit in soil gas and on the probe to source distance and soil and NAPL characteristics. According to the obtained results, the radon deficit technique results a feasible method for a qualitative identification of residual NAPL when radon in soil gas is measured at distances lower than 2 m from the contaminated zone. However, for an accurate quantitative estimation of the NAPL phase content, soil gas probes should be preferably located at distances lower than 1 m from the source zone.

On the role of rock matrix to heat transfer in a fracture-rock matrix system

In this study, a fully coupled analytical model is developed for thermal energy transfer in a single fracture-rock matrix system where the coupling implies that the governing equations of thermal transfer in the fracture and rock matrix are supplemented with the continuity conditions of temperature and thermal flux at the interfaces of the fracture-rock matrix. The proposed model accounts for thermal convection, longitudinal thermal conduction and thermal dispersion in the fracture, and transverse thermal conduction in the rock matrix. The fully coupled two-dimensional model is established to investigate the thermal energy transfer processes, assess the spatiotemporal temperature distribution in the fracture and rock matrix system and investigate the role of the rock matrix. The solutions are verified with the existing studies and proven to be accurate and robust. The present study demonstrates that: 1) thermal dispersion in the fracture plays an important role in the temperature distribution in the fracture and rock matrix domains, and longitudinal thermal conduction in the fracture has minor effects on the temperature distribution in the system; 2) transverse thermal conduction in the rock matrix is a critical parameter that determines the spatiotemporal temperature distribution in both the fracture and the rock matrix domains. Ignoring thermal conduction in the rock matrix will lead to a significant overestimation of temperature in the short and long terms; 3) the sensitivity analysis implies that thermal energy transfer in the system is sensitive to the fluid velocity in the fracture, thermal dispersivity in the fracture and thermal conductivity in the rock matrix, and less sensitive to thermal conductivity in the fracture.

Evaluating dual-domain models for upscaling multicomponent reactive transport in mine waste rock

Reactive transport models have proven abilities to simulate the quantity and quality of drainage from mine waste rock. Tracer experiments indicate the presence of fast and slow flow regimes in many heterogeneous waste-rock piles. Although multidomain models have been developed specifically for systems with such distinctive hydrodynamics, there have been limited applications of multidomain reactive transport models to simulate composite drainage chemistries from waste-rock piles to date. This work evaluated the ability of dual-domain multicomponent reactive transport models (DDMRTMs) to reproduce breakthrough curves of conservative (chloride) and reactive (molybdenum) solutes observed at a well-characterized experimental waste-rock pile at the Antamina Mine, Peru. We found that the DDMRTM simulations quantitatively matched eight-year-long records of conservative transport through the waste-rock pile when parameterized mainly with field-measured properties obtained from the site and limited calibration. The DDMRTM model also provided a reasonable match to field observations of the reactive solute. The limited calibrated parameters are physically realistic, corroborating the ability of these multidomain models to reproduce the complex reactive-transport processes governing polluted rock drainage from large-scale waste-rock piles.

Impact of solute charge and diffusion coefficient on electromigration and mixing in porous media

The application of electrokinetic techniques in porous media has great potential to enhance mass transfer rates and, thus, to mobilize contaminants and effectively deliver reactants and amendments. However, the transport mechanisms induced by the application of an external electric field are complex and entail the coupling of physical, chemical and electrostatic processes. In this study we focus on electromigration and we provide experimental evidence of the impact of compound-specific properties, such as the aqueous diffusivity and the valence of charged species, on the macroscopic electrokinetic transport. We performed a series of multidimensional experiments considering the displacement of three different tracer plumes (i.e., permanganate, allura red and new coccine) in different background electrolyte solutions. The outcomes of the experiments clearly show that both the compound-specific diffusivity and the charge of the injected and resident ions impact the transport of the selected color tracer plumes, whose evolution was monitored with image analysis. The investigated experimental scenarios led to distinct plume behavior characterized by different mass distribution, average displacement velocities, longitudinal and lateral plume spreading, shape of the invading and receding fronts, as well as dilution of the injected solutes. A numerical simulator, based on the Nernst-Planck-Poisson equations and on aqueous speciation reactions in the pore water, allowed us to quantitatively interpret the experimental results, to capture the observed patterns of plume evolution, and to illuminate the coupling between the governing physico-chemical mechanisms and the controlling role of small scale compound-specific and electrostatic properties. Finally, the model was also extended to a typical configuration of in situ electrokinetic remediation of contaminated groundwater to show the impact of such mechanisms at larger scale.

Impact of compost reactive layer on hydraulic transport and C & N cycles: Biogeochemical modeling of infiltration column experiments

Managed Aquifer Recharge (MAR) is a key strategy to increase freshwater resources in many regions facing water scarcity. MAR issues are related to both quantity and quality of the infiltrating water. In most countries, very high quality of the infiltrating water is required, to limit the impact on the aquifer geochemistry. In this paper, the possibility of injecting water of lower quality in the aquifer and letting the biogeochemical reactions take place in order to enhance its quality is explored. Here, we present the fate of nutrients (C, N) in the biogeochemical system of a reactive barrier formed by mixture of different proportions of sand and compost, supplied with treated wastewater to mimic MAR. An integrated conceptual model involving the nutrient cycles and biomass dynamics (auto- and heterotrophic) was developed, and then tested with a number of solute transport experiments in columns with different compost fraction in the column filling. The model incorporated both saturation and inhibition processes (regarding the nutrients and their byproducts) to provide a comprehensive picture of the nutrient dynamics within the column. The model developed (three if considering the 3 column setups) allowed to discriminate the processes that govern the fate of nutrients in relation with the compost enhancing long-term nutrient degradation, yet hindering hydraulic parameters that affect infiltration rates.

Arsenic release and transport during oxidative dissolution of spatially-distributed sulfide minerals

The oxidative dissolution of sulfide minerals, naturally present in the subsurface, is one of the major pathways of arsenic mobilization. This study investigates the release and fate of arsenic from arsenopyrite and löllingite oxidation under dynamic redox conditions. We performed multidimensional flow-through experiments focusing on the impact of chemical heterogeneity on arsenic mobilization and reactive transport. In the experimental setups the As-bearing sulfide minerals were embedded, with different concentrations and spatial distributions, into a sandy matrix under anoxic conditions. Oxic water flushed in the flow-through setups triggered the oxidative dissolution of the reactive minerals, the release of arsenic, as well as changes in pore water chemistry, surface-solution interactions and mineral precipitation. We developed a reactive transport model to quantitatively interpret the experimental results. The simulation outcomes showed that 40% of the arsenic released was reincorporated into a freshly precipitated iron-arsenate phase that created a coating on the mineral surface limiting the dissolution reactions. The faster dissolution rate of löllingite compared to arsenopyrite was responsible for sustaining the continuous release of As-contaminated plumes. The model also allowed shedding light on the spatial distribution, on the temporal dynamics, and on the interactions between arsenic sources (As-bearing minerals) and sinks (freshly formed secondary phases) in flow-through systems.

Modeling of microbial kinetics and mass transfer in bioreactors simulating the natural attenuation of arsenic and iron in acid mine drainage

Natural attenuation in acid mine drainage (AMD) due to biological iron and arsenic oxidation offers a promising strategy to treat As-rich AMD in passive bioreactors. A reactive transport model is developed in order to identify the main controlling factors. It simulates batch and flow-through experiments that reproduce natural attenuation in a high-As AMD. The 2-D model couples second-order microbial kinetics (Fe- and As- oxidation) and geochemical reactions to hydrodynamic transport. Oxidation only occurrs in the biofilm with an oxygen transfer from the air through the water column. The model correctly simulates the Fe(II)-Fe(III) and As(III)-As(V) concentrations in the outlet waters and the precipitates, over hydraulic retention times from 30 min to 800 min. It confirms that the natural attenuation at 20 °C is driven by the fast Fe(II) oxidation and slow As(III) oxidation that favors arsenite trapping by schwertmannite over amorphous ferric arsenate (AFA) formation. The localization of iron oxidation in the biofilm limits the attenuation of arsenic and iron as the water column height increases. The change in the composition of the bacterial iron-oxidizer community of the biofilm at the lowest pH boundary seems to control the Fe(II) oxidation kinetic rate besides the bacterial concentration.

Seasonal dynamics modifies fate of oxygen, nitrate, and organic micropollutants during bank filtration - temperature-dependent reactive transport modeling of field data

Bank filtration is considered to improve water quality through microbially mediated degradation of pollutants and is suitable for waterworks to increase their production. In particular, aquifer temperatures and oxygen supply have a great impact on many microbial processes. To investigate the temporal and spatial behavior of selected organic micropollutants during bank filtration in dependence of relevant biogeochemical conditions, we have set up a 2D reactive transport model using MODFLOW and PHT3D under the user interface ORTI3D. The considered 160-m-long transect ranges from the surface water to a groundwater extraction well of the adjacent waterworks. For this purpose, water levels, temperatures, and chemical parameters were regularly measured in the surface water and groundwater observation wells over one and a half years. To simulate the effect of seasonal temperature variations on microbial mediated degradation, we applied an empirical temperature factor, which yields a strong reduction of the degradation rate at groundwater temperatures below 11 °C. Except for acesulfame, the considered organic micropollutants are substantially degraded along their subsurface flow paths with maximum degradation rates in the range of 10^-6 mol L^-1 s^-1. Preferential biodegradation of phenazone, diclofenac, and valsartan was found under oxic conditions, whereas carbamazepine and sulfamethoxazole were degraded under anoxic conditions. This study highlights the influence of seasonal variations in oxygen supply and temperature on the fate of organic micropollutants in surface water infiltrating into an aquifer.

Identifying remedial solutions through optimal bioremediation design under real-world field conditions

Over more than a century of intense industrial production and associated accidental release, petroleum products (e.g., gasoline, diesel, fuel oil) have contaminated a significant portion of the world's groundwater resources. Groundwater remediation is generally a complex task, especially where aquifers and the associated contaminant distribution are highly heterogeneous. The ability to predict the efficiency of such remediation is of crucial importance, as the costs are strongly linked to the treatment design and duration. In this study, a coupled simulation-optimization (S/O) framework, consisting of a process-based reactive transport simulation model linked with particle swarm optimization (PSO) was developed. It was subsequently applied for the design of a real-world in situ bio-treatment of a BTEX contaminated aquifer in France. In the application, the optimization framework was used to simultaneously determine optimal well locations and their optimal injection rates, both constituting key elements of the enhanced biodegradation design problem. The optimization of the treatment efficiency was examined in terms of three different regulatory objectives, (1) minimization of the residual NAPL mass of the key contaminant, i.e., benzene, in the source zone, (2) reduction of the maximum concentration of benzene in groundwater, and (3) minimization of the time required to reduce the benzene concentration in groundwater to below a threshold value. Our analysis of potential, optimal remediation strategies showed that: (i) the complexity of the biodegradation behavior at real sites may favor very different remediation options as a result of varying remediation targets, (ii) the long term behavior of the contaminants after the end of the active treatment period, which is often neglected, showed to have a significant influence on remediation design that requires increased attention, (iii) PSO has shown to be a very efficient algorithm in the context of the present study. The insights that can be gained from such a framework will provide decision support to select the most suitable remediation strategy while facing different regulatory objectives.

Impact of diffuse layer processes on contaminant forward and back diffusion in heterogeneous sandy-clayey domains

Low-permeability aquitards can significantly affect the transport, distribution, and persistence of contaminant plumes in subsurface systems. Although such low-permeability materials are often charged, the key role of charge-induced electrostatic processes during contaminant transport has not been extensively studied. This work presents a detailed investigation exploring the coupled effects of heterogeneous distribution of physical, chemical and electrostatic properties on reactive contaminant transport in field-scale groundwater systems including spatially distributed clay zones. We performed an extensive series of numerical experiments in three distinct heterogeneous sandy-clayey domains with different levels of complexity. The flow and reactive transport simulations were performed by explicitly resolving the complex velocity fields, the small-scale electrostatic processes, the compound-specific diffusive/dispersive fluxes and the chemical processes utilizing a multi-continua based reactive transport code (MMIT-Clay). In each particular domain, numerical experiments were performed focusing on both the forward and back diffusion through the sandy-clayey interfaces. The results illuminate the control of microscopic electrostatic mechanisms on macroscopic mass transfer. Coulombic interactions in the clay's diffuse layer can significantly accelerate or retard a particular species depending on its charge. Furthermore, the chemical heterogeneity plays a major role in mass storage and release during reactive transport. Neglecting such processes can lead to substantial over- or underestimation of the overall transport behavior, which underlines the need for integrated physical, chemical and electrostatic approaches to accurately describe mass transfer processes in systems including low-permeability inclusions.

Experiment and multicomponent model based analysis on the effect of flow rate and nitrate concentration on denitrification in low-permeability media

To better understand the combined effects of flow rate and NO₃^- concentration on denitrification rate and NO₃^- removal efficiency in the low-permeability media, a set of column experiments with different flow rates and injected NO₃^- concentrations were conducted. Denitrification processes under these different conditions were simulated using the PHREEQC code that couples the biogeochemical reactions and hydrological transport processes. In these reactive transport models, Monod kinetics were applied to describe the denitrification process. It was found that, among the experiments conducted in this study, the low flow rate (0.023 m/d) resulted in the low denitrification rate but high NO₃^- removal efficiency. Meanwhile, NO₃^- removal efficiency was the highest (85%) under moderate NO₃^- concentration of 1.29 mmol/L, although denitrification rate increased in response to the increase of NO₃^- concentration. The model results also indicated that NO₃^- removal efficiency of 97% can be achieved with relatively low flow rate and high influent NO₃^- concentration. The results in this study provide insights into NO₃^- remediation, and the temporal and spatial flow rate, as well as NO₃^- concentration distribution, should be pre-evaluated for the effective removal strategies.

Process-based modeling of electrokinetic-enhanced bioremediation of chlorinated ethenes

This study presents a process-based modeling analysis of electrokinetic-enhanced bioremediation (EK-Bio) to illuminate the complex interactions between physical, electrostatic and biogeochemical processes occurring during the application of this remediation technique. The features of the proposed model include: (i) multidimensional electrokinetic transport in saturated porous media by electromigration and electroosmosis, (ii) charge interactions, (iii) degradation kinetics, (iv) microbial populations dynamics of indigenous and specialized exogenous degraders, (v) mass transfer limitations, and (vi) geochemical reactions. A scenario modeling investigation is presented, which was inspired by an EK-Bio pilot application conducted in a clayey aquitard at the Skuldelev site (Denmark) contaminated by chlorinated ethenes. Lactate and specialized degraders are delivered under conservative and reactive transport conditions. In the considered setup, transport of lactate using electrokinetics results in more than fourfold increase in the distribution efficiency with respect to a diffusion-only scenario. Moreover, EK transport by electromigration and electroosmosis yields fluxes at least two orders of magnitude larger than diffusive fluxes. Quantitative metrics are also defined and used to assess the amendment distribution and the enhanced contaminant biodegradation in the different conservative and reactive transport scenarios.

Spatiotemporal evolution of iron and sulfate concentrations during riverbank filtration: Field observations and reactive transport modeling

Riverbank filtration is a commonly-used technology that improves water quality by passing river water through aquifers. In this study, a riverbank filtration site in Busan, South Korea, was investigated to understand the spatiotemporal evolution of high iron and sulfate concentrations observed in the riverbank-filtered water. Discrepancies between the nonreactive transport results and field measurements suggest that iron-sulfate-related geochemical reactions play a major role in the spatiotemporal evolution of the hydrochemical properties. Pyrite oxidation was hypothesized to be the main process driving the release of iron and sulfate. To test this hypothesis, a reactive transport model was developed, that implemented pyrite oxidation as a kinetic process and subsequent ferrous iron oxidation and ferric iron precipitation as equilibrium processes. The model accurately captured the temporal evolution of sulfate; however, iron concentrations were underestimated. Sensitivity tests revealed that adjusting reaction constants significantly improved the prediction of iron concentrations. The results of this study suggest that pyrite oxidation can affect the hydrochemistry of riverbank-filtered water and highlight the potential limitations of using theoretical reaction constants in field modeling applications.

Model-based interpretation of methane oxidation and respiration processes in landfill biocovers: 3-D simulation of laboratory and pilot experiments

Landfill biocovers are an efficient strategy for the mitigation of greenhouse gas emissions from landfills. A complex interplay between key physical and reactive processes occurs in biocovers and affects the transport of gas components. Therefore, numerical models can greatly help the understanding of these systems, their design and optimal operation. In this study, we developed a 3-D multicomponent modeling approach to quantitatively interpret experimental datasets measured in the laboratory and in pilot-scale landfill biocovers. The proposed model is able to reproduce the observed spatial and temporal dynamics of CH₄, O₂ and CO₂ migration in biocovers under different operating conditions and demonstrates the importance of dimensionality in understanding the propagation of gas flow and migration of gas components in such porous media. The model allowed us to capture the coupled transport behavior of gas components, to evaluate the exchange of gas fluxes at the interface between the biocover surface and free air flow, and to investigate the effects of different gas injection patterns on the distribution of gas components within biocovers. The model also helps elucidating the dynamics and competition between methane oxidation and respiration processes observed in the different experimental setups. The simulation outcomes reveal that increasing availability of methane (i.e., higher injection flow rates or higher fractions of CH₄ in the injected gas composition) results in progressive dominance of methane oxidation in the biocovers and moderates the impact of respiration.

Modeling the relationship of aeration, oxygen transfer and treatment performance in aerated horizontal flow treatment wetlands

Mechanical aeration is commonly used to improve the overall treatment efficacy of constructed wetlands. However, the quantitative relationships of air flow rate (AFR), water temperature, field oxygen transfer and treatment performance have not been analyzed in detail until today. In this study, a reactive transport model based on dual-permeability flow and biokinetic formulations of the Constructed Wetland Model No. 1 (CWM1) was developed and extented to 1) simulate oxygen transfer and treatment performance for organic carbon and nitrogen of two pilot-scale horizontal flow (HF) aerated wetlands (Test and Control) treating domestic sewage, and, 2) to investigate the dependence of oxygen transfer and treatment performance on AFR and water temperature. Both pilot-scale wetlands exhibited preferential flow patters and high treatment performance for chemical oxygen demand (COD) and NH₄-N at AFRs of 128-700 L m^-2 h^-1. A reduction of the AFR in the Test system from 128 to 72 L h^-1 m^-2 substantially inhibited NH₄-N removal. Conservative tracer transport as well as reactive transport of dissolved oxygen (DO), soluble and total chemical oxygen demand (COD_s, COD_t), NH₄-N and NO_x-N measured in pilot-scale experiments were simulated with acceptable accuracy (E₁¯=0.39±0.26). An equation to estimate the volumetric oxygen transfer coefficient was found to be: k_La,20=0.511ln(AFR). Simulated treatment performance depended on k_La,20 in a non-linear manner. A local sensitivity analysis of the calibrated parameters revealed porosity, hydraulic permeability and dispersion length of the fast flow field as well as k_La,20 as most important. An optimal AFR for a spatially and temporally continuous aeration pattern for treatment wetlands treating similar influent was estimated to 150-200 L h^-1 m^-2. This study provides insights into aeration mechanisms of aerated treatment wetlands and highlights the benefits of process modeling for in-depth system analysis.

Applying aluminum-organic matter precipitates to reduce soil permeability in-situ: A field and modeling study

Using naturally occurring processes to modify the engineering properties of the subsurface has received increasing attention from industry and research communities as they aid in the development of cost-effective, robust and sustainable engineering technologies. In line with this trend, we propose to use precipitates of aluminum (Al) and organic matter (OM) to reduce soil permeability in-situ. This process is inspired by podzolization: a soil stratification process where a layer with low permeability is developed at depth via the precipitation of metal-OM complexes. In this study, the concept of applying Al-OM precipitates for in-situ soil permeability reduction was for the first time applied in the field. The aim of the field experiment was to create a cylindrical flow barrier in a sand layer at depth. In order to design and engineer the field application, we performed a series of scenario analyses with a site-specific 3D reactive transport model. This led to an in-situ engineering approach where a flow barrier was created by separate injection of Al and OM using a combined injection/extraction strategy. During the field application, the local variation of soil conditions required significant modifications to the design. Further scenario analyses with the model were conducted to adapt the original design and to understand the consequences of these modifications. The results show that a cylindrical flow barrier was created after an injection period of 8 days. The precipitation of Al-OM is a highly localized process, where large amount of precipitates is formed in the close vicinity of the injection filter screens. Evaluation of pumping tests that were performed after the injection activities revealed that the permeability of the treated sand was reduced to 2% of its original value. This first full-scale field test demonstrates that applying Al-OM precipitates is a suitable bio-based engineering tool to reduce soil permeability in-situ.

Scale dependence of surface complexation capacity and rates in heterogeneous media

Reactive solutes flow, transport, and transform in natural subsurface that is ubiquitously heterogeneous. How much are the heterogeneity effects on solute reactive transport? At what spatial scale heterogeneity effects diminish? Answers to these questions are important to understand controls of and to predict solute reactive transport in natural environments. Although these questions have been explored for non-reactive tracers and solutes that go through relatively simple single-component sorption, heterogeneity effects remain poorly understood for solutes that go through surface complexation. These reactions typically include multi-component and are subject to influences of local geochemical conditions. This work aims to address these questions and to understand how and how much effective capacity (C_sc,m), early (k_c1) and late (k_c2) rates of Cr(VI) surface complexation are influenced by heterogeneity characteristics and length scales. A two-dimensional model with Cr(VI) surface complexation were calibrated using column data. Numerical experiments were carried out in heterogeneous media characterized by different permeability variance σ²_lnκ (0.2, 4.5), correlation length λ_L (0-6.0 cm), and domain length (0.1-1.0 m). Results show that long λ_L and high σ²_lnκ facilitate the formation of preferential flow and elongate the diffusion length, leading to high connectivity that minimizes C_sc,m and k_c1. In highly connected media (>2.0), C_sc,m and k_c1 can be >2 orders of magnitude lower than those in corresponding homogeneous media. Connectivity however is not a good indicator for k_c2 that is predominantly controlled by diffusion. The heterogeneity effects are scale dependent. As the domain length increases, C_sc,m and k_c1 values approach those of homogeneous media due to longer mixing lengths and continued pH increase during Cr(VI) surface complexion. The "critical length" at which heterogeneity effects disappear occurs at about 20λ_L and differs only slightly at different σ²_lnκ values. The heterogeneity effects however linger at much longer lengths for k_c2, the values of which are >2 orders of magnitude lower in heterogeneous media even at 20λ because of the longer characteristic times of diffusive transport. These results underscore the strong regulation of solute reactive transport by mineral spatial patterns, suggesting the importance of subsurface characteristics and spatial scales at which they are important in estimating natural attenuation, residence times, and natural system capacity to sorb chemicals.

Mineralogy controls on reactive transport of Marcellus Shale waters

Produced or flowback waters from Marcellus Shale gas extraction (MSWs) typically are highly saline and contain chemicals including trace metals, which pose significant concerns on water quality. The natural attenuation of MSW chemicals in groundwater is poorly understood due to the complex interactions between aquifer minerals and MSWs, limiting our capabilities to monitor and predict. Here we combine flow-through experiments and process-based reactive transport modeling to understand mechanisms and quantify the retention of MSW chemicals in a quartz (Qtz) column, a calcite-rich (Cal) column, and a clay-rich (Vrm, vermiculite) column. These columns were used to represent sand, carbonate, and clay-rich aquifers. Results show that the types and extent of water-rock interactions differ significantly across columns. Although it is generally known that clay-rich media retard chemicals and that quartz media minimize water-rock interactions, results here have revealed insights that differ from previous thoughts. We found that the reaction mechanisms are much more complex than merely sorption and mineral precipitation. In clay rich media, trace metals participate in both ion exchange and mineral precipitation. In fact, the majority of metals (~50-90%) is retained in the solid via mineral precipitation, which is surprising because we typically expect the dominance of sorption in clay-rich aquifers. In the Cal column, trace metals are retained not only through precipitation but also solid solution partitioning, leading to a total of 75-99% retention. Even in the Qtz column, trace metals are retained at unexpectedly high percentages (~20-70%) due to precipitation. The reactive transport model developed here quantitatively differentiates the relative importance of individual processes, and bridges a limited number of experiments to a wide range of natural conditions. This is particularly useful where relatively limited knowledge and data prevent the prediction of complex rock-contaminant interactions and natural attenuation.

Benthic nitrite exchanges in the Seine River (France): An early diagenetic modeling analysis

Nitrite is a toxic intermediate compound in the nitrogen (N) cycle. Elevated concentrations of nitrite have been observed in the Seine River, raising questions about its sources and fate. Here, we assess the role of bottom sediments as potential sources or sinks of nitrite along the river continuum. Sediment cores were collected from two depocenters, one located upstream, the other downstream, from the largest wastewater treatment plant (WWTP) servicing the conurbation of Paris. Pore water profiles of oxygen, nitrate, nitrite and ammonium were measured. Ammonium, nitrate and nitrite fluxes across the sediment-water interface (SWI) were determined in separate core incubation experiments. The data were interpreted with a one-dimensional, multi-component reactive transport model, which accounts for the production and consumption of nitrite through nitrification, denitrification, anammox and dissimilatory nitrate reduction to ammonium (DNRA). In all core incubation experiments, nitrate uptake by the sediments was observed, indicative of high rates of denitrification. In contrast, for both sampling locations, the sediments in cores collected in August 2012 acted as sinks for nitrite, but those collected in October 2013 released nitrite to the overlying water. The model results suggest that the first step of nitrification generated most pore water nitrite at the two locations. While nitrification was also the main pathway consuming nitrite in the sediments upstream of the WWTP, anammox dominated nitrite removal at the downstream site. Sensitivity analyses indicated that the magnitude and direction of the benthic nitrite fluxes most strongly depend on bottom water oxygenation and the deposition flux of labile organic matter.

Kinetic and reactive transport modeling

Predicting the fates of arsenic (As) and selenium (Se) in natural geologic media like rocks and soils necessitates the understanding of how their various oxyanionic species behave and migrate under dynamic conditions. In this study, geochemical factors and processes crucial in the leaching and transport of arsenite (As^III), arsenate (As^V), selenite (Se^IV) and selenate (Se^VI) in tunnel-excavated rocks of marine origin were investigated using microscopic/extraction techniques, column experiments, dissolution-precipitation kinetics and one-dimensional reactive transport modeling. The results showed that evaporite salts were important because aside from containing As and Se, they played crucial roles in the evolution of pH and concentrations of coexisting ions, both of which had strong effects on adsorption-desorption reactions of As and Se species with iron oxyhydroxide minerals/phases. The observed leaching trends of As^V, As^III, Se^IV and Se^VI were satisfactorily simulated by one-dimensional reactive transport models, which predict that preferential adsorptions of As^V and Se^IV were magnified by geochemical changes in the columns due to water flow. Moreover, our results showed that migrations of As^III, Se^IV and Se^VI could be predicted adequately by 1D solute transport with simple activity-K'_d approach, but surface complexation was more reliable to simulate adsorption-desorption behavior of As^V.

Modeling 3D-CSIA data: Carbon, chlorine, and hydrogen isotope fractionation during reductive dechlorination of TCE to ethene

Reactive transport modeling of multi-element, compound-specific isotope analysis (CSIA) data has great potential to quantify sequential microbial reductive dechlorination (SRD) and alternative pathways such as oxidation, in support of remediation of chlorinated solvents in groundwater. As a key step towards this goal, a model was developed that simulates simultaneous carbon, chlorine, and hydrogen isotope fractionation during SRD of trichloroethene, via cis-1,2-dichloroethene (and trans-DCE as minor pathway), and vinyl chloride to ethene, following Monod kinetics. A simple correction term for individual isotope/isotopologue rates avoided multi-element isotopologue modeling. The model was successfully validated with data from a mixed culture Dehalococcoides microcosm. Simulation of Cl-CSIA required incorporation of secondary kinetic isotope effects (SKIEs). Assuming a limited degree of intramolecular heterogeneity of δ³⁷Cl in TCE decreased the magnitudes of SKIEs required at the non-reacting Cl positions, without compromising the goodness of model fit, whereas a good fit of a model involving intramolecular CCl bond competition required an unlikely degree of intramolecular heterogeneity. Simulation of H-CSIA required SKIEs in H atoms originally present in the reacting compounds, especially for TCE, together with imprints of strongly depleted δ²H during protonation in the products. Scenario modeling illustrates the potential of H-CSIA for source apportionment.

Nitrate reduction and its effects on trichloroethylene degradation by granular iron

Laboratory column experiments and reactive transport modeling were performed to evaluate the reduction of nitrate and its effects on trichloroethylene (TCE) degradation by granular iron. In addition to determining degradation kinetics of TCE in the presence of nitrate, the columns used in this study were equipped with electrodes which allowed for in situ measurements of corrosion potentials of the iron material. Together with Raman spectroscopic measurements the mechanisms of decline in iron reactivity were examined. The experimental results showed that the presence of nitrate resulted in an increase in corrosion potential and the formation of thermodynamically stable passive films on the iron surface which impaired iron reactivity. The extent of the decline in iron reactivity was proportional to the nitrate concentration. Consequently, significant decreases in TCE and nitrate degradation rates and migration of degradation profiles for both compounds occurred. Furthermore, the TCE degradation kinetics deviated from the pseudo-first-order model. The results of reactive transport modeling, which related the amount of a passivating iron oxide, hematite (α-Fe₂O₃), to the reactivity of iron, were generally consistent with the patterns of migration of TCE and nitrate profiles observed in the column experiments. More encouragingly, the simulations successfully demonstrated the differences in performances of three columns without changing model parameters other than concentrations of nitrate in the influent. This study could be valuable in the design of iron permeable reactive barriers (PRBs) or in the development of effective maintenance procedures for PRBs treating TCE-contaminated groundwater with elevated nitrate concentrations.

Reactive transport of uranium in fractured crystalline rock: Upscaling in time and distance

Batch adsorption and breakthrough column experiments were conducted to evaluate uranium transport through altered material that fills fractures in a granite rock system at the Grimsel Test Site in Switzerland at pH 6.9 and 7.9. The role of adsorption and desorption kinetics was evaluated with reactive transport modeling by comparing one-, two-, and three-site models. Emphasis was placed on describing long desorption tails that are important for upscaling in time and distance. The effect of increasing pH in injection solutions was also evaluated. For pH 6.9, a three-site model with forward rate constants between 0.07 and 0.8 ml g(-1) h(-1), reverse rate constants between 0.001 and 0.06 h(-1), and site densities of 1.3, 0.104, and 0.026 μmol g(-1) for 'weak/fast', 'strong/slow', and 'very strong/very slow' sites provided the best fits. For pH 7.9, a three-site model with forward rate constants between 0.05 and 0.8 mL g(-1) h(-1), reverse rate constants between 0.001 and 0.6 h(-1), and site densities of 1.3, 0.039, and 0.013 μmol g(-1) for a 'weak/fast', 'strong/slow', and 'very strong/very slow' sites provided the best fits. Column retardation coefficients (Rd) were 80 for pH 6.9 and 10.3 for pH 7.9. Model parameters determined from the batch and column experiments were used in 50 year large-scale simulations for continuous and pulse injections and indicated that a three-site model is necessary at pH 6.9, although a Kd-type equilibrium partition model with one-site was adequate for large scale predictions at pH 7.9. Batch experiments were useful for predicting early breakthrough times in the columns while column experiments helped differentiate the relative importance of sorption sites and desorption rate constants on transport.

Reactive transport modeling of subsurface arsenic removal systems in rural Bangladesh

Subsurface Arsenic Removal (SAR) is a technique for in-situ removal of arsenic from groundwater. Extracted groundwater is aerated and re-injected into an anoxic aquifer, where the oxygen in the injected water reacts with ferrous iron in the aquifer to form hydrous ferric oxide (HFO). Subsequent extraction of groundwater contains temporarily lower As concentrations, because As sorbs onto the HFO. Injection, storage, and extraction together is called a cycle. A reactive transport model (RTM) was developed in PHREEQC to determine the hydrogeochemical processes responsible for As (im)mobilization during experimental SAR operation performed in Bangladesh. Oxidation of Fe(II) and As(III) were modeled using kinetic-rate expressions. Cation exchange, precipitation of HFO, and surface complexation, were modeled as equilibrium processes. A best set of surface complexation reactions and corresponding equilibrium constants was adopted from previous studies to simulate all 20 cycles of a SAR experiment. The model gives a reasonable match with observed concentrations of different elements in the extracted water (e.g., the r(2) value of As was 0.59 or higher). As concentrations in the extracted water are governed by four major processes. First, As concentration decreases in response to the elevated pH of injection water and likewise increases when native neutral pH groundwater flows in. Second, the sorption capacity for As increases due to the gradual buildup of HFO. Third, As sorption is enhanced by preferential removal of As(V). Fourth, competitive sorption of Si limits the capacity of freshly precipitated HFO for As sorption. Transferability of the developed reactive transport model was demonstrated through successful application of the model, without further calibration, to two additional SAR sites in Bangladesh. This gives confidence that the model could be useful to assess potential SAR performance at locations in Bangladesh based on local hydrogeochemical conditions.

Evolution of carbon isotope signatures during reactive transport of hydrocarbons in heterogeneous aquifers

Compound-specific isotope analysis (CSIA) of organic pollutants has become a well-established tool for assessing the occurrence and extent of biodegradation processes in contaminated aquifers. However, the precision of CSIA is influenced by the degree to which assumptions underlying CSIA data interpretation hold under realistic field-scale conditions. For the first time this study demonstrates how aquifer analogs combined with reactive transport models offer an underexplored way to develop generic process understanding, evaluate monitoring and quantification strategies in highly heterogeneous subsurface settings. Data from high-resolution aquifer analogs were used in numerical experiments to track the propagation of a representative oxidizable organic compound (toluene) within a variety of realistic heterogeneous aquifers and to investigate its detailed fate. The simulations were used to analyze (1) the effects of physical aquifer heterogeneities on spatiotemporal patterns of contaminant concentrations and isotope signatures, (2) the performance of the commonly applied Rayleigh equation and (3) the applicability of an extension of the Rayleigh equation for complex hydrogeological conditions. The results indicate that if field-derived enrichment factors are applied without corrections for dilution, the conventional Rayleigh equation is inaccurate and estimates for biodegradation are typically overestimated and unreliable in heterogeneous aquifers. Underestimations can occur due to the partial source zone depletion. In contrast, if dilution can be accurately accounted for, field-derived enrichment factors comprise a suitable alternative to laboratory-derived and redox-specific enrichment factors. The study also examines to what extent variations in monitoring/sampling strategies influence the obtained results. Especially measurements from long-screened wells (>1 m) reveal to be inappropriate for the application of the Rayleigh equation in the investigated aquifer analogs, as low resolution data sampled from the simulated scenarios only enable a qualitative assessment of biodegradation. Measurements from both long- and short-screened wells employing the Rayleigh equation streamline approach are only partly viable for in situ biodegradation measurements in heterogeneous systems.

Reactive transport modeling of ⁹⁰Sr sorption in reactive sandpacks

Strontium-90 ((90)Sr) is one of the most problematic radioactive contaminants in groundwater at nuclear sites. Although (90)Sr is retarded relative to groundwater flow, it is sufficiently mobile and long-lived to require treatment in many hydrogeological settings. A detailed study was performed on the practicality of using granular clinoptilolite as a sandpack around groundwater wells where groundwater is contaminated with (90)Sr and the water table must be lowered. The effectiveness of the reactive sandpack concept and the mechanisms controlling (90)Sr attenuation was investigated by numerical analysis of data obtained from four in situ column experiments. The experiments spanned the range of pore-water velocities that would occur during radial flow through granular clinoptilolite sandpacks. A kinetic sorption model was required to adequately reproduce the experimentally observed (90)Sr behavior. Calibrated first-order kinetic rates were correlated with pore-water velocities. After calibration, three sorption models were used to simulate (90)Sr attenuation for four hypothetical pumping scenarios. Results show that a velocity-dependent kinetic model accurately simulates the observed early breakthrough for high pore-water velocities. The results indicate (1) that reactive sandpacks have good potential for in situ remediation and construction dewatering and (2) that quantitative modeling can aid in the design and application of this novel technique.

Transport of selenium oxyanions through TiO2porous media: column experiments and multi-scale modeling

The present work deals with the modeling of selenium oxyanion (selenite/selenate) retention in TiO2 rutile porous media. A set of chemical interactions was elaborated from spectroscopic measurements and adsorption experiments in batch and column reactors, and a model of transient transport of the selenium species through laboratory column was developed. The adsorption model considered that both forms of selenium (Se) compete for the same sorption sites, hydroxyl groups, allowing taking into account a competitive adsorption. Stoichiometry and equilibrium constants of adsorption reactions were determined on the basis of spectrometric measurement and adsorption isotherm curve fitting. This approach led to a model of Sips type isotherm including a pH-dependence. It offers an excellent fitting compared to the classical Langmuir equation and provides a unique set of parameters for both oxyanions. IMPACT code and associated modeling method were then used to couple transport and chemical reactions. The obtained numerical results showed a reasonable prediction of the shape and the time location of selenium oxyanions and pH breakthrough curves.

Influence of calcite on uranium(VI) reactive transport in the groundwater-river mixing zone

Calcite is an important, relatively soluble mineral phase that can affect uranium reactive transport in subsurface sediments. This study was conducted to investigate the distribution of calcite and its influence on uranium adsorption and reactive transport in the groundwater-river mixing zone of the Hanford 300A site, Washington State. Simulations using a two-dimensional (2D) reactive transport model under field-relevant hydrological and hydrogeochemical conditions revealed the development of a calcite reaction front through the mixing zone as a result of dynamic groundwater-river interactions. The calcite concentration distribution, in turn, affected the concentrations of aqueous carbonate and calcium, and pH through dissolution, as river waters intruded and receded from the site at different velocities in response to stage changes. The composition variations in groundwater subsequently influenced uranium mobility and discharge rates into the river in a complex fashion. The results implied that calcite distribution and concentration are important variables that need to be quantified for accurate reactive transport predictions of uranium, especially in dynamic groundwater-river mixing zones.

Modeling phosphate transport and removal in a compact bed filled with a mineral-based sorbent for domestic wastewater treatment

Phosphorus filter units containing mineral-based sorbents with a high phosphate (PO4) binding capacity have been shown to be appropriate for removing PO4 in the treatment of domestic wastewater in on-site facilities. However, a better understanding of their PO4 removal mechanisms, and reactions that could lead to the formation of PO4 compounds, is required to evaluate the potential utility of candidate sorbents. Models based on data obtained from laboratory-scale experiments with columns of selected materials can be valuable for acquiring such understanding. Thus, in this study the transport and removal of PO4 in experiments with a laboratory-scale column filled with a commercial silicate-based sorbent were modeled, using the hydro-geochemical transport code PHREEQC. The resulting models, that incorporated the dissolution of calcite, kinetic constrains for the dissolution of calcium oxide (CaO) and wollastonite (CaSiO3), and the precipitation of amorphous tricalcium phosphate, Ca3(PO4)2, successfully simulated the removal of PO4 observed in the experiments.

Estimating pathway-specific contributions to biodegradation in aquifers based on dual isotope analysis: theoretical analysis and reactive transport simulations

At field sites with varying redox conditions, different redox-specific microbial degradation pathways contribute to total contaminant degradation. The identification of pathway-specific contributions to total contaminant removal is of high practical relevance, yet difficult to achieve with current methods. Current stable-isotope-fractionation-based techniques focus on the identification of dominant biodegradation pathways under constant environmental conditions. We present an approach based on dual stable isotope data to estimate the individual contributions of two redox-specific pathways. We apply this approach to carbon and hydrogen isotope data obtained from reactive transport simulations of an organic contaminant plume in a two-dimensional aquifer cross section to test the applicability of the method. To take aspects typically encountered at field sites into account, additional simulations addressed the effects of transverse mixing, diffusion-induced stable-isotope fractionation, heterogeneities in the flow field, and mixing in sampling wells on isotope-based estimates for aerobic and anaerobic pathway contributions to total contaminant biodegradation. Results confirm the general applicability of the presented estimation method which is most accurate along the plume core and less accurate towards the fringe where flow paths receive contaminant mass and associated isotope signatures from the core by transverse dispersion. The presented method complements the stable-isotope-fractionation-based analysis toolbox. At field sites with varying redox conditions, it provides a means to identify the relative importance of individual, redox-specific degradation pathways.

Identifying key controls on the behavior of an acidic-U(VI) plume in the Savannah River Site using reactive transport modeling

Acidic low-level waste radioactive waste solutions were discharged to three unlined seepage basins at the F-Area of the Department of Energy (DOE) Savannah River Site (SRS), South Carolina, USA, from 1955 through 1989. Despite many years of active remediation, the groundwater remains acidic and contaminated with significant levels of U(VI) and other radionuclides. Monitored Natural Attenuation (MNA) is a desired closure strategy for the site, based on the premise that regional flow of clean background groundwater will eventually neutralize the groundwater acidity, immobilizing U(VI) through adsorption. An in situ treatment system is currently in place to accelerate this in the downgradient portion of the plume and similar measures could be taken upgradient if necessary. Understanding the long-term pH and U(VI) adsorption behavior at the site is critical to assess feasibility of MNA along with the in-situ remediation treatments. This paper presents a reactive transport (RT) model and uncertainty quantification (UQ) analyses to explore key controls on the U(VI)-plume evolution and long-term mobility at this site. Two-dimensional numerical RT simulations are run including the saturated and unsaturated (vadose) zones, U(VI) and H(+) adsorption (surface complexation) onto sediments, dissolution and precipitation of Al and Fe minerals, and key hydrodynamic processes are considered. UQ techniques are applied using a new open-source tool that is part of the developing ASCEM reactive transport modeling and analysis framework to: (1) identify the complex physical and geochemical processes that control the U(VI) plume migration in the pH range where the plume is highly mobile, (2) evaluate those physical and geochemical parameters that are most controlling, and (3) predict the future plume evolution constrained by historical, chemical and hydrological data. The RT simulation results show a good agreement with the observed historical pH and concentrations of U(VI), nitrates and Al concentrations at multiple locations. Mineral dissolution and precipitation combined with adsorption reactions on goethite and kaolinite (the main minerals present with quartz) could buffer pH at the site for long periods of time. UQ analysis using the Morris one-at-a-time (OAT) method indicates that the model/parameter is most sensitive to the pH of the waste solution, discharge rates, and the reactive surface area available for adsorption. However, as a key finding, UQ analysis also indicates that this model (and parameters) sensitivity evolves in space and time, and its understanding could be crucial to assess the temporal efficiency of a remediation strategy in contaminated sites. Results also indicate that residual U(VI) and H(+) adsorbed in the vadose zone, as well as aquifer permeability, could have a significant impact on the acidic plume long-term mobility.

Biogeochemical processes controlling the mobility of major ions and trace metals in aquitard sediments beneath an oil sand tailing pond: laboratory studies and reactive transport modeling

Increased production and expansion of the oil sand industry in Alberta are of great benefit to the economy, but they carry major environmental challenges. The volume of fluid fine tailings requiring storage is 840×10(6) m(3) and growing, making it imperative that we better understand the fate and transport of oil sand process-affected water (OSPW) seepage from these facilities. Accordingly, the current study seeks to characterize both a) the potential for major ion and trace element release, and b) the principal biogeochemical processes involved, as tailing pond OSPW infiltrates into, and interacts with, underlying glacial till sediments prior to reaching down gradient aquifers or surface waters. Objectives were addressed through a series of aqueous and solid phase experiments, including radial diffusion cells, an isotope analysis, X-ray diffraction, and sequential extractions. The diffusion cells were also simulated in a reactive transport framework to elucidate key reaction processes. The experiments indicate that the ingress and interaction of OSPW with the glacial till sediment-pore water system will result in: a mitigation of ingressing Na (retardation), displacement and then limited precipitation of exchangeable Ca and Mg (as carbonates), sulfate reduction and subsequent precipitation of the produced sulfides, as well as biodegradation of organic carbon. High concentrations of ingressing Cl (~375 mg L(-1)) and Na (~575 mg L(-1)) (even though the latter is delayed, or retarded) are expected to migrate through the till and into the underlying sand channel. Trace element mobility was influenced by ion exchange, oxidation-reduction, and mineral phase reactions including reductive dissolution of metal oxyhydroxides - in accordance with previous observations within sandy aquifer settings. Furthermore, although several trace elements showed the potential for release (Al, B, Ba, Cd, Mn, Pb, Si, Sr), large-scale mobilization is not supported. Thus, the present results suggest that in addition to the commonly cited naphthenic acids, remediation of OSPW-impacted groundwater will need to address high concentrations of major ions contributing to salinization.