Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment

Three-dimensional, coupled variably saturated flow and biogeochemical reactive transport modeling of a 2008 in situ uranium bioremediation field experiment is used to better understand the interplay of transport and biogeochemical reactions controlling uranium behavior under pulsed acetate amendment, seasonal water table variation, spatially variable physical (hydraulic conductivity, porosity) and geochemical (reactive surface area) material properties. While the simulation of the 2008 Big Rusty acetate biostimulation field experiment in Rifle, Colorado was generally consistent with behaviors identified in previous field experiments at the Rifle IFRC site, the additional process and property detail provided several new insights. A principal conclusion from this work is that uranium bioreduction is most effective when acetate, in excess of the sulfate-reducing bacteria demand, is available to the metal-reducing bacteria. The inclusion of an initially small population of slow growing sulfate-reducing bacteria identified in proteomic analyses led to an additional source of Fe(II) from the dissolution of Fe(III) minerals promoted by biogenic sulfide. The falling water table during the experiment significantly reduced the saturated thickness of the aquifer and resulted in reactants and products, as well as unmitigated uranium, in the newly unsaturated vadose zone. High permeability sandy gravel structures resulted in locally high flow rates in the vicinity of injection wells that increased acetate dilution. In downgradient locations, these structures created preferential flow paths for acetate delivery that enhanced local zones of TEAP reactivity and subsidiary reactions. Conversely, smaller transport rates associated with the lower permeability lithofacies (e.g., fine) and vadose zone were shown to limit acetate access and reaction. Once accessed by acetate, however, these same zones limited subsequent acetate dilution and provided longer residence times that resulted in higher concentrations of TEAP reaction products when terminal electron donors and acceptors were not limiting. Finally, facies-based porosity and reactive surface area variations were shown to affect aqueous uranium concentration distributions with localized effects of the fine lithofacies having the largest impact on U(VI) surface complexation. The ability to model the comprehensive biogeochemical reaction network, and spatially and temporally variable processes, properties, and conditions controlling uranium behavior during engineered bioremediation in the naturally complex Rifle IFRC subsurface system required a subsurface simulator that could use the large memory and computational performance of a massively parallel computer. In this case, the eSTOMP simulator, operating on 128 processor cores for 12h, was used to simulate the 110-day field experiment and 50 days of post-biostimulation behavior.     

Large-scale modeling of reactive solute transport in fracture zones of granitic bedrocks

Final disposal of high-level radioactive waste in deep repositories located in fractured granite formations is being considered by several countries. The assessment of the safety of such repositories requires using numerical models of groundwater flow, solute transport and chemical processes. These models are being developed from data and knowledge gained from in situ experiments such as the Redox Zone Experiment carried out at the underground laboratory of Äspö in Sweden. This experiment aimed at evaluating the effects of the construction of the access tunnel on the hydrogeological and hydrochemical conditions of a fracture zone intersected by the tunnel. Most chemical species showed dilution trends except for bicarbonate and sulphate which unexpectedly increased with time. Molinero and Samper [Molinero, J. and Samper, J. Groundwater flow and solute transport in fracture zones: an improved model for a large-scale field experiment at Äspö (Sweden). J. Hydraul. Res., 42, Extra Issue, 157–172] presented a two-dimensional water flow and solute transport finite element model which reproduced measured drawdowns and dilution curves of conservative species. Here we extend their model by using a reactive transport which accounts for aqueous complexation, acid–base, redox processes, dissolution–precipitation of calcite, quartz, hematite and pyrite, and cation exchange between Na⁺ and Ca²⁺. The model provides field-scale estimates of cation exchange capacity of the fracture zone and redox potential of groundwater recharge. It serves also to identify the mineral phases controlling the solubility of iron. In addition, the model is useful to test the relevance of several geochemical processes. Model results rule out calcite dissolution as the process causing the increase in bicarbonate concentration and reject the following possible sources of sulphate: (1) pyrite dissolution, (2) leaching of alkaline sulphate-rich waters from a nearby rock landfill and (3) dissolution of iron monosulphides contained in Baltic seafloor sediments. Based on these results, microbially mediated processes are postulated as the most likely hypothesis to explain the measured increase of dissolved bicarbonates and sulphates after tunnel construction.     

Modeling in-situ uranium(VI) bioreduction by sulfate-reducing bacteria

We present a travel-time based reactive transport model to simulate an in-situ bioremediation experiment for demonstrating enhanced bioreduction of uranium(VI). The model considers aquatic equilibrium chemistry of uranium and other groundwater constituents, uranium sorption and precipitation, and the microbial reduction of nitrate, sulfate and U(VI). Kinetic sorption/desorption of U(VI) is characterized by mass transfer between stagnant micro-pores and mobile flow zones. The model describes the succession of terminal electron accepting processes and the growth and decay of sulfate-reducing bacteria, concurrent with the enzymatic reduction of aqueous U(VI) species. The effective U(VI) reduction rate and sorption site distributions are determined by fitting the model simulation to an in-situ experiment at Oak Ridge, TN. Results show that (1) the presence of nitrate inhibits U(VI) reduction at the site; (2) the fitted effective rate of in-situ U(VI) reduction is much smaller than the values reported for laboratory experiments; (3) U(VI) sorption/desorption, which affects U(VI) bioavailability at the site, is strongly controlled by kinetics; (4) both pH and bicarbonate concentration significantly influence the sorption/desorption of U(VI), which therefore cannot be characterized by empirical isotherms; and (5) calcium-uranyl-carbonate complexes significantly influence the model performance of U(VI) reduction.     

The addition of organic carbon and nitrate affects reactive transport of heavy metals in sandy aquifers

Organic carbon introduction in the soil to initiate remedial measures, nitrate infiltration due to agricultural practices or sulphate intrusion owing to industrial usage can influence the redox conditions and pH, thus affecting the mobility of heavy metals in soil and groundwater. This study reports the fate of Zn and Cd in sandy aquifers under a variety of plausible in-situ redox conditions that were induced by introduction of carbon and various electron acceptors in column experiments. Up to 100% Zn and Cd removal (from the liquid phase) was observed in all the four columns, however the mechanisms were different. Metal removal in column K1 (containing sulphate), was attributed to biological sulphate reduction and subsequent metal precipitation (as sulphides). In the presence of both nitrate and sulphate (K2), the former dominated the process, precipitating the heavy metals as hydroxides and/or carbonates. In the presence of sulphate, nitrate and supplemental iron (Fe(OH)(3)) (K3), metal removal was also due to precipitation as hydroxides and/or carbonates. In abiotic column, K4, (with supplemental iron (Fe(OH)(3)), but no nitrate), cation exchange with soil led to metal removal. The results obtained were modeled using the reactive transport model PHREEQC-2 to elucidate governing processes and to evaluate scenarios of organic carbon, sulphate and nitrate inputs.     

Intercomparison of reactive transport models applied to UO2 oxidative dissolution and uranium migration

Oxidative dissolution of uranium dioxide (UO(2)) and the subsequent migration of uranium in a subsurface environment and an underground waste disposal have been simulated with reactive transport models. In these systems, hydrogeological and chemical processes are closely entangled and their interdependency has been analyzed in detail, notably with respect to redox reactions, kinetics of mineralogical evolution and hydrodynamic migration of species of interest. Different codes, where among CASTEM, CHEMTRAP and HYTEC, have been used as an intercomparison and verification exercise. Although the agreement between codes is satisfactory, it is shown that the discretization method of the transport equation (i.e. finite elements (FE) versus mixed-hybrid FE and finite differences) and the sequential coupling scheme may lead to systematic discrepancies.     

Intercomparison of reactive transport models applied to UO2 oxidative dissolution and uranium migration

Oxidative dissolution of uranium dioxide (UO₂) and the subsequent migration of uranium in a subsurface environment and an underground waste disposal have been simulated with reactive transport models. In these systems, hydrogeological and chemical processes are closely entangled and their interdependency has been analyzed in detail, notably with respect to redox reactions, kinetics of mineralogical evolution and hydrodynamic migration of species of interest.Different codes, where among CASTEM, CHEMTRAP and HYTEC, have been used as an intercomparison and verification exercise. Although the agreement between codes is satisfactory, it is shown that the discretization method of the transport equation (i.e. finite elements (FE) versus mixed-hybrid FE and finite differences) and the sequential coupling scheme may lead to systematic discrepancies.     

Transport and anaerobic biodegradation of propylene glycol in gravel-rich soil materials

Continued input of airplane de-icing/anti-icing fluids (ADAF) to runway adjacent soils may result in the depletion of soil-borne terminal electron acceptors. We studied the transport and transformation of propylene glycol (PG), the major constituent of many ADAF, in topsoil and subsoil samples using saturated column experiments at 4 degrees C and 20 degrees C. The export of soil-borne DOC was generally high, non-exhaustive and rate limited. Retardation of added PG was negligible. Rapid PG degradation was observed only in topsoil materials high in organic matter at 20 degrees C. At 4 degrees C, no significant degradation was observed. Thus, under unfavorable, i.e., wet and cold conditions typical for winter de-icing operations, PG and its metabolites will be relocated to deeper soil horizons or even to the groundwater. In subsoil materials, PG degradation was very slow and incomplete. We found that subsoil degradation depended on the import of active microorganisms originating from the organic-rich topsoil material. The degradation efficiency is strongly influenced by the flow velocity, i.e., the residence time of PG in the soil column. Poorly crystalline iron(III) and manganese(IV) (hydr)oxides are used during microbial respiration acting as terminal electron acceptors. This results in the formation and effective relocation of reduced and mobile Fe and Mn species. Long-term application of ADAF to runway adjacent soil as well as the lasting consumption of Fe and Mn will tend to decrease the soil redox potential. Without proper counteractive measures, this will eventually favor the development of methanogenic conditions.     

Mechanisms of electron acceptor utilization: implications for simulating anaerobic biodegradation

Simulation of biodegradation reactions within a reactive transport framework requires information on mechanisms of terminal electron acceptor processes (TEAPs). In initial modeling efforts, TEAPs were approximated as occurring sequentially, with the highest energy-yielding electron acceptors (e.g. oxygen) consumed before those that yield less energy (e.g., sulfate). Within this framework in a steady state plume, sequential electron acceptor utilization would theoretically produce methane at an organic-rich source and Fe(II) further downgradient, resulting in a limited zone of Fe(II) and methane overlap. However, contaminant plumes often display much more extensive zones of overlapping Fe(II) and methane. The extensive overlap could be caused by several abiotic and biotic processes including vertical mixing of byproducts in long-screened monitoring wells, adsorption of Fe(II) onto aquifer solids, or microscale heterogeneity in Fe(III) concentrations. Alternatively, the overlap could be due to simultaneous utilization of terminal electron acceptors. Because biodegradation rates are controlled by TEAPs, evaluating the mechanisms of electron acceptor utilization is critical for improving prediction of contaminant mass losses due to biodegradation. Using BioRedox-MT3DMS, a three-dimensional, multi-species reactive transport code, we simulated the current configurations of a BTEX plume and TEAP zones at a petroleum-contaminated field site in Wisconsin. Simulation results suggest that BTEX mass loss due to biodegradation is greatest under oxygen-reducing conditions, with smaller but similar contributions to mass loss from biodegradation under Fe(III)-reducing, sulfate-reducing, and methanogenic conditions. Results of sensitivity calculations document that BTEX losses due to biodegradation are most sensitive to the age of the plume, while the shape of the BTEX plume is most sensitive to effective porosity and rate constants for biodegradation under Fe(III)-reducing and methanogenic conditions. Using this transport model, we had limited success in simulating overlap of redox products using reasonable ranges of parameters within a strictly sequential electron acceptor utilization framework. Simulation results indicate that overlap of redox products cannot be accurately simulated using the constructed model, suggesting either that Fe(III) reduction and methanogenesis are occurring simultaneously in the source area, or that heterogeneities in Fe(III) concentration and/or mineral type cause the observed overlap. Additional field, experimental, and modeling studies will be needed to address these questions.     

Concurrent nitrate and Fe(III) reduction during anaerobic biodegradation of phenols in a sandstone aquifer

The biodegradation of phenols (5, 60, 600 mg l⁻¹) under anaerobic conditions (nitrate enriched and unamended) was studied in laboratory microcosms with sandstone material and groundwater from within an anaerobic ammonium plume in an aquifer. The aqueous phase was sampled and analyzed for phenols and selected redox sensitive parameters on a regular basis. An experiment with sandstone material from specific depth intervals from a vertical profile across the ammonium plume was also conducted. The miniature microcosms used in this experiment were sacrificed for sampling for phenols and selected redox sensitive parameters at the end of the experiment. The sandstone material was characterized with respect to oxidation and reduction potential and Fe(II) and Fe(III) speciation prior to use for all microcosms and at the end of the experiments for selected microcosms.The redox conditions in the anaerobic microcosms were mixed nitrate and Fe(III) reducing. Nitrate and Fe(III) were apparently the dominant electron acceptors at high and low nitrate concentrations, respectively. When biomass growth is taken into account, nitrate and Fe(III) reduction constituted sufficient electron acceptor capacity for the mineralization of the phenols observed to be degraded even at an initial phenols concentration of 60 mg l⁻¹ (high) in an unamended microcosm, whereas nitrate reduction alone is unlikely to have provided sufficient electron acceptor capacity for the observed degradation of the phenols in the unamended microcosm.For microcosm systems, with solid aquifer materials, dissolution of organic substances from the solid material may occur. A quantitative determination of the speciation (mineral types and quantity) of electron acceptors associated with the solids, at levels relevant for degradation of specific organic compounds in aquifers, cannot always be obtained. Hence, complete mass balances of electron acceptor consumption for specific organic compounds degradation are difficult to confine. For aquifer materials with low initial Fe(II) content, Fe(II) determinations on solids and in aqueous phase samples may provide valuable information on Fe(III) reduction. However, in microcosms with natural sediments and where electron acceptors are associated with the sediments, complete mass-balances for substrates and electron acceptors are not likely to be obtained.     

Modelling of iron cycling and its impact on the electron balance at a petroleum hydrocarbon contaminated site in Hnevice, Czech Republic

Over a period of several decades multiple leaks of large volumes from storage facilities located near Hnevice (Czech Republic) have caused the underlying Quaternary aquifer to be severely contaminated with nonaqueous phase liquid (NAPL) petroleum hydrocarbons. Beginning in the late 1980's the NAPL plume started to shrink as a consequence of NAPL dissolution exceeding replenishment and due to active remediation. The subsurface was classified geochemically into four different zones, (i) a contaminant-free zone never occupied by NAPL or dissolved contaminants, (ii) a re-oxidation zone formerly occupied by NAPL, (iii) a zone currently occupied by NAPL, and (iv) a lower fringe zone between the overlying NAPL and the deeper underlying contaminant-free zone. The study investigated the spatial and temporal variability of the redox zonation at the Hnevice site and quantified the influence of iron-cycling on the overall electron balance. As a first step inverse geochemical modelling was carried out to identify possible reaction models and mass transfer processes. In a subsequent step, two-dimensional (forward) multi-component reactive transport modelling was performed to evaluate and quantify the major processes that control the geochemical evolution at the site. The study explains the observed enrichment of the lower fringe zone with ferrihydrite as a result of the re-oxidation of ferrous iron. It suggests that once the NAPL zone started to shrink the dissolution of previously formed siderite and FeS by oxygen and nitrate consumed a significant part of the oxidation capacity for a considerable time period and therefore limited the penetration of electron acceptors into the NAPL contaminated zone.     

Kinetic modeling of microbially-driven redox chemistry of radionuclides in subsurface environments: coupling transport, microbial metabolism and geochemistry

Microbial reactions play an important role in regulating pore water chemistry as well as secondary mineral distribution in many subsurface systems and, therefore, may directly impact radionuclide migration in those systems. This paper presents a general modeling approach to couple microbial metabolism, redox chemistry, and radionuclide transport in a subsurface environment. To account for the likely achievement of quasi-steady state biomass accumulations in subsurface environments, a modification to the traditional microbial growth kinetic equation is proposed. The conditions for using biogeochemical models with or without an explicit representation of biomass growth are clarified. Based on the general approach proposed in this paper, the couplings of uranium reactions with biogeochemical processes are incorporated into computer code BIORXNTRN Version 2.0. The code is then used to simulate a subsurface contaminant migration scenario, in which a water flow containing both uranium and a complexing organic ligand is recharged into an oxic carbonate aquifer. The model simulation shows that Mn and Fe oxyhydroxides may vary significantly along a flow path. The simulation also shows that uranium(VI) can be reduced and therefore immobilized in the anoxic zone created by microbial degradation.     

电动力学技术强化原位生物修复研究进展

介绍了利用电动力学技术强化土壤及地下水原位生物修复的原理和最新进展。电动力学强化的基本原理是利用电渗析、电迁移和电泳等电动力学效应加速污染环境中有机污染物和微生物运动,注入营养物、电子受体或活性微生物,或者利用电极反应和电流热效应为地下生物降解创造有利条件。研究表明．电动力学技术能有效地强化原位生物修复,而且该技术不破环生态环境．安装和操作简单,成本低廉．有广泛的应用前景。     

Uranium removal from groundwater via in situ biostimulation: Field-scale modeling of transport and biological processes

During 2002 and 2003, bioremediation experiments in the unconfined aquifer of the Old Rifle UMTRA field site in western Colorado provided evidence for the immobilization of hexavalent uranium in groundwater by iron-reducing Geobacter sp. stimulated by acetate amendment. As the bioavailable Fe(III) terminal electron acceptor was depleted in the zone just downgradient of the acetate injection gallery, sulfate-reducing organisms came to dominate the microbial community. In the present study, we use multicomponent reactive transport modeling to analyze data from the 2002 field experiment to identify the dominant transport and biological processes controlling uranium mobility during biostimulation, and determine field-scale parameters for these modeled processes. The coupled process simulation approach was able to establish a quantitative characterization of the principal flow, transport, and reaction processes based on the 2002 field experiment, that could be applied without modification to describe the 2003 field experiment. Insights gained from this analysis include field-scale estimates of the bioavailable Fe(III) mineral threshold for the onset of sulfate reduction, and rates for the Fe(III), U(VI), and sulfate terminal electron accepting processes.     

Bacterial interactions with uranium: An environmental perspective

The presence of actinides in radioactive wastes is of major concern because of their potential for migration from the waste repositories and long-term contamination of the environment. Studies have been and are being made on inorganic processes affecting the migration of radionuclides from these repositories to the environment but it is becoming increasingly evident that microbial processes are of importance as well. Bacteria interact with uranium through different mechanisms including, biosorption at the cell surface, intracellular accumulation, precipitation, and redox transformations (oxidation/reduction). The present study is intended to give a brief overview of the key processes responsible for the interaction of actinides e.g. uranium with bacterial strains isolated from different extreme environments relevant to radioactive repositories. Fundamental understanding of the interaction of these bacteria with U will be useful for developing appropriate radioactive waste treatments, remediation and long-term management strategies as well as for predicting the microbial impacts on the performance of the radioactive waste repositories.     

Release of trace metals, sulfate and complexed cyanide from soils contaminated with gas-purifier wastes: a microcosm study

Deposited gas-purifier wastes are commonly contaminated with trace metals, sulfate and cyanide (CN) compounds. We investigated their release from three soils contaminated with gas-purifier wastes into solution in microcosm experiments under varying redox conditions (E(H) 170-620 mV). The soils differed in pH (2.2; 4.9; 7.4) and featured low amounts of trace metals, but large amounts of total S and total CN. The pH governed trace metal release in the case of the acidic soil and CN release in the case of the slightly alkaline soil. The redox potential controlled trace metal and CN release in the case of the moderately acidic soil. Sources of dissolved SO(4)(2-) were dissolution of gypsum, desorption from Fe oxides and probably oxidation of elemental S. The geochemical behaviors of trace metals (soluble under acidic and reducing conditions) and CN (soluble under alkaline and oxidizing conditions) were diametrically opposed.     

Evaluating factors that influence microbial phenanthrene biodegradation rates by regression with categorical variables

To advance the accuracy of bioremediation measurements, it is useful before specific experiments to attribute or estimate the influence of both experimental as well as field conditions on the expected magnitudes of microbial degradation rate coefficients. This paper analyzes the numerical contribution, or influence, of categories of conditions, such as bacterial adaptive state, electron acceptor type, mixing, generalized sorption conditions, and biodegradation temperature, on published phenanthrene biodegradation rates as an example of our regression approach. A fundamental microbial degradation rate equation is transformed to an additive model, then using multiple linear regression on published data, coefficients (of categorical variables) and a linear model are presented that estimate first-order biodegradation rate coefficients to within a factor of 3. Numerical estimates of how much bacterial adaptive state and presence of a sorption phase, the two most statistically significant factors, alter the phenanthrene biodegradation rate are presented. The influence of some measurement or field conditions, for example, the influence of oxygen reduction versus optimal nitrate reduction, cannot be distinguished statistically given the available data and range. The regression model is tested using conditions from newly published papers to estimate a priori the expected rate, which compares very favorably to measurements reported in the papers. Due to limited published data and range for extreme cases, the current coefficients do not apply to degradation of very aged phenanthrene nor very low concentrations of electron acceptors. As estimating tools, however, the coefficients themselves and the regression approach have very beneficial roles in design of experiments for both laboratory and field settings. Our method can be applied to other PAHs as sufficient data become available.     

Conceptual and numerical model of uranium(VI) reductive immobilization in fractured subsurface sediments

A conceptual model and numerical simulations of bacterial U(VI) reduction in fractured subsurface sediments were developed to assess the potential feasibility of biomineralization at the fracture/matrix interface as a mechanism for immobilization of uranium in structured subsurface media. The model envisions flow of anaerobic groundwater, with or without acetate as an electron donor for stimulation of U(VI) reduction by dissimilatory metal-reducing bacteria (DMRB), within mobile macropores along a one-dimensional flow path. As the groundwater moves along the flow path, U(VI) trapped in the immobile mesopore and micropore domains (the sediment matrix) becomes desorbed and transferred to the mobile macropores (fractures) via a first-order exchange mechanism. By allowing bacterial U(VI) reduction to occur in the mesopore domain (assumed to account for 12% of total sediment pore volume) according to experimentally-determined kinetic parameters and an assumed DMRB abundance of 10(7) cells per cm3 bulk sediment (equivalent to 4 mg of cells per dm3 bulk sediment), the concentration of U(VI) in the macropore domain was reduced ca. 10-fold compared to that predicted in the absence of mesopore DMRB activity after a 6-month simulation period. The results suggest that input of soluble electron donors over a period of years could lead to a major redistribution of uranium in fractured subsurface sediments, converting potentially mobile sorbed U(VI) to an insoluble reduced phase (i.e. uraninite) in the mesopore domain that is expected to be permanently immobile under sustained anaerobic conditions.     

Biogeochemistry and isotope geochemistry of a landfill leachate plume

The biogeochemical processes were identified which improved the leachate composition in the flow direction of a landfill leachate plume (Banisveld, The Netherlands). Groundwater observation wells were placed at specific locations after delineating the leachate plume using geophysical tests to map subsurface conductivity. Redox processes were determined using the distribution of solid and soluble redox species, hydrogen concentrations, concentration of dissolved gases (N(2), Ar, and CH(4)), and stable isotopes (delta15N-NO(3), delta34S-SO(4), delta13C-CH(4), delta2H-CH(4), and delta13C of dissolved organic and inorganic carbon (DOC and DIC, respectively)). The combined application of these techniques improved the redox interpretation considerably. Dissolved organic carbon (DOC) decreased downstream in association with increasing delta13C-DOC values confirming the occurrence of degradation. Degradation of DOC was coupled to iron reduction inside the plume, while denitrification could be an important redox process at the top fringe of the plume. Stable carbon and hydrogen isotope signatures of methane indicated that methane was formed inside the landfill and not in the plume. Total gas pressure exceeded hydrostatic pressure in the plume, and methane seems subject to degassing. Quantitative proof for DOC degradation under iron-reducing conditions could only be obtained if the geochemical processes cation exchange and precipitation of carbonate minerals (siderite and calcite) were considered and incorporated in an inverse geochemical model of the plume. Simulation of delta13C-DIC confirmed that precipitation of carbonate minerals happened.     

Multi-component reactive transport modeling of natural attenuation of an acid groundwater plume at a uranium mill tailings site

Natural attenuation of an acidic plume in the aquifer underneath a uranium mill tailings pond in Wyoming, USA was simulated using the multi-component reactive transport code PHREEQC. A one-dimensional model was constructed for the site and the model included advective-dispersive transport, aqueous speciation of 11 components, and precipitation-dissolution of six minerals. Transport simulation was performed for a reclamation scenario in which the source of acidic seepage will be terminated after 5 years and the plume will then be flushed by uncontaminated upgradient groundwater. Simulations show that successive pH buffer reactions with calcite, Al(OH)3(a), and Fe(OH)3(a) create distinct geochemical zones and most reactions occur at the boundaries of geochemical zones. The complex interplay of physical transport processes and chemical reactions produce multiple concentration waves. For SO4(2-) transport, the concentration waves are related to advection-dispersion, and gypsum precipitation and dissolution. Wave speeds from numerical simulations compare well to an analytical solution for wave propagation.