Influence of salinity gradients on the diffusion of water and ionic species in dual porosity clay samples.

Most of the available data on diffusion in natural clayey rocks consider tracer diffusion in the absence of a salinity gradient despite the fact that such gradients are frequently found in natural and engineered subsurface environments. To assess the role of such gradients on the diffusion properties of clayey materials, through-diffusion experiments were carried out in the presence and absence of a salinity gradient using salt-diffusion and radioisotope tracer techniques. The experiments were carried out with vermiculite samples that contained equal proportions of interparticle and interlayer porosities so as to assess also the role played by the two types of porosities on the diffusion of water and ions. Data were interpreted using both a classical Fickian diffusion model and with a reactive transport code, CrunchClay that can handle multi-porosity diffusion processes in the presence of charged surfaces. By combining experimental and simulated data, we demonstrated that (i) the flux of water diffusing through vermiculite interlayer porosity was minor compared to that diffusing through the interparticle porosity, and (ii) a model considering at least three types of porous volumes (interlayer, interparticle diffuse layer, and bulk interparticle) was necessary to reproduce consistently the variations of neutral and charged species diffusion as a function of salinity gradient conditions.

Resolving experimental biases in the interpretation of diffusion experiments with a user-friendly numerical reactive transport approach.

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (De), clay porosities (ε), and adsorption distribution coefficients (KD) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (De-ε-KD) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a De-ε-KD modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface-subsurface interactions.

We computationally explore the relationship between surface-subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface-subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface-subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface-subsurface exchange rates lead to a nonlinear change in the degree of connectivity-quantified through the number of disconnected clusters of ponding water-in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

Modeling geogenic and atmospheric nitrogen through the East River Watershed, Colorado Rocky Mountains.

There is a growing understanding of the role that bedrock weathering can play as a source of nitrogen (N) to soils, groundwater and river systems. The significance is particularly apparent in mountainous environments where weathering fluxes can be large. However, our understanding of the relative contributions of rock-derived, or geogenic, N to the total N supply of mountainous watersheds remains poorly understood. In this study, we develop the High-Altitude Nitrogen Suite of Models (HAN-SoMo), a watershed-scale ensemble of process-based models to quantify the relative sources, transformations, and sinks of geogenic and atmospheric N through a mountain watershed. Our study is based in the East River Watershed (ERW) in the Upper Colorado River Basin. The East River is a near-pristine headwater watershed underlain primarily by an N-rich Mancos Shale bedrock, enabling the timing and magnitude of geogenic and atmospheric contributions to watershed scale dissolved N-exports to be quantified. Several calibration scenarios were developed to explore equifinality using >1600 N concentration measurements from streams, groundwater, and vadose zone samples collected over the course of four years across the watershed. When accounting for recycling of N through plant litter turnover, rock weathering accounts for approximately 12% of the annual dissolved N sources to the watershed in the most probable calibration scenario (0-31% in other scenarios), and 21% (0-44% in other scenarios) when considering only "new" N sources (i.e. geogenic and atmospheric). On an annual scale, instream dissolved N elimination, plant turnover (including cattle grazing) and atmospheric deposition are the most important controls on N cycling.

Dissolved Carbonate and pH Control the Dissolution of Uranyl Phosphate Minerals in Flow-Through Porous Media.

Uranyl phosphate minerals represent an important secondary source of uranium release at contaminated sites. In flow-through column experiments with background porewater (BPW) of typical freshwater aquifer composition (pH 7.0, ∼0.2 mM total carbonate (TC)), dissolution of K-ankoleite (KUO2PO4·3H2O), Na-autunite (NaUO2PO4·3H2O), and Ca-autunite (Ca(UO2)2(PO4)2·6H2O) was controlled by mineral solubility at steady-state U release. Effluent concentrations indicated exchange with BPW cations, and postreaction characterization showed alteration of the initial mineral composition, changes in structure (decreased crystallinity, increased disorder, and distortion of U-P mineral sheets) and possible neoformation of phases of similar structure. Increasing the BPW pH and TC to 8.1-8.2 and 2.2-3.7 mM, respectively, resulted in mineral undersaturation and produced ca. 2 orders-of-magnitude higher U and P release without reaching steady state. Minerals incorporated less BPW cations into their structures compared to low carbonate BPW experiments but showed structural disorder and distortion. Faster dissolution rates were attributed to the formation of binary and ternary uranyl carbonate complexes that accelerate the rate-determining step of uranyl detachment from the uranyl-phosphate layered structure. Calculated dissolution rates (log Rs between -8.95 and -10.32 mol m-2 s-1), accounting for reaction and transport in porous media, were similar to dissolution rates of other classes of uranyl minerals. In undersaturated solutions, dissolution rates for uranyl phosphate, oxyhydroxide, and silicate minerals can be predicted within 1-2 orders-of-magnitude from pH ∼5-10 on the basis of pH/carbonate concentration.

Incorporating Nanoscale Effects into a Continuum-Scale Reactive Transport Model for CO2-Deteriorated Cement.

Wellbore cement deterioration is critical for wellbore integrity and the safety of CO2 storage in geologic formations. Our previous experimental work highlighted the importance of the portlandite (CH)-depleted zone and the surface dissolution zone in the CO2-attacked cement. In this study, we simulated numerically the evolution of the CH-depleted zone and the dissolution of the cement surfaces utilizing a reduced-dimension (1D) reactive transport model. The approach shows that three nanoscale effects are important and had to be incorporated in a continuum-scale model to capture experimental observations: First, it was necessary to account for the fact that secondary CaCO3 precipitation does not fill the pore space completely, with the result that acidic brine continues to diffuse through the carbonated zone to form a CH-depleted zone. Second, secondary precipitation in brine begins via nucleation kinetics, and this could not be described with previous models using growth kinetics alone. Third, our results suggest that the CaCO3 precipitates in the confined pore space are more soluble than those formed in brine. This study provides a new platform for a reduced dimension model for CO2 attack on cement that captures the important nanoscale mechanisms influencing macroscale phenomena in subsurface environments.

Reoxidation of Chromium(III) Products Formed under Different Biogeochemical Regimes.

Hexavalent chromium, Cr(VI), is a widespread and toxic groundwater contaminant. Reductive immobilization to Cr(III) is a treatment option, but its success depends on the long-term potential for reduced chromium precipitates to remain immobilized under oxidizing conditions. In this unique long-term study, aquifer sediments subjected to reductive Cr(VI) immobilization under different biogeochemical regimes were tested for their susceptibility to reoxidation. After reductive treatment for 1 year, sediments were exposed to oxygenated conditions for another 2 years in flow-through, laboratory columns. Under oxidizing conditions, immobilized chromium reduced under predominantly denitrifying conditions was mobilized at low concentrations (≪1 µM Cr(VI); ∼ 3% of Cr(III) deposited) that declined over time. A conceptual model of a limited pool of more soluble Cr(III), and a larger pool of relatively insoluble Cr(III), is proposed. In contrast, almost no chromium was mobilized from columns reduced under predominantly fermentative conditions, and where reducing conditions persisted for several months after introduction of oxidizing conditions, presumably due to the presence of a reservoir of reduced species generated during reductive treatment. The results from this 3-year study demonstrate that biogeochemical conditions present during reductive treatment, and the potential for buildup of reducing species, will impact the long-term sustainability of the remediation effort.

Water Table Dynamics and Biogeochemical Cycling in a Shallow, Variably-Saturated Floodplain.

Three-dimensional variably saturated flow and multicomponent biogeochemical reactive transport modeling, based on published and newly generated data, is used to better understand the interplay of hydrology, geochemistry, and biology controlling the cycling of carbon, nitrogen, oxygen, iron, sulfur, and uranium in a shallow floodplain. In this system, aerobic respiration generally maintains anoxic groundwater below an oxic vadose zone until seasonal snowmelt-driven water table peaking transports dissolved oxygen (DO) and nitrate from the vadose zone into the alluvial aquifer. The response to this perturbation is localized due to distinct physico-biogeochemical environments and relatively long time scales for transport through the floodplain aquifer and vadose zone. Naturally reduced zones (NRZs) containing sediments higher in organic matter, iron sulfides, and non-crystalline U(IV) rapidly consume DO and nitrate to maintain anoxic conditions, yielding Fe(II) from FeS oxidative dissolution, nitrite from denitrification, and U(VI) from nitrite-promoted U(IV) oxidation. Redox cycling is a key factor for sustaining the observed aquifer behaviors despite continuous oxygen influx and the annual hydrologically induced oxidation event. Depth-dependent activity of fermenters, aerobes, nitrate reducers, sulfate reducers, and chemolithoautotrophs (e.g., oxidizing Fe(II), S compounds, and ammonium) is linked to the presence of DO, which has higher concentrations near the water table.

Divergent aquifer biogeochemical systems converge on similar and unexpected Cr(VI) reduction products.

In this study of reductive chromium immobilization, we found that flow-through columns constructed with homogenized aquifer sediment and continuously infused with lactate, chromate, and various native electron acceptors diverged to have very different Cr(VI)-reducing biogeochemical regimes characterized by either denitrifying or fermentative conditions (as indicated by effluent chemical data, 16S rRNA pyrotag data, and metatranscriptome data). Despite the two dramatically different biogeochemical environments that evolved in the columns, these regimes created similar Cr(III)-Fe(III) hydroxide precipitates as the predominant Cr(VI) reduction product, as characterized by micro-X-ray fluorescence and micro-X-ray absorption near-edge structure analysis. We discuss two conflicting scenarios of microbially mediated formation of Cr(III)-Fe(III) precipitates, each of which is both supported and contradicted by different lines of evidence: (1) enzymatic reduction of Cr(VI) to Cr(III) followed by coprecipitation of Cr(III) and Fe(III) and (2) both regimes generated at least small amounts of Fe(II), which abiotically reduced Cr(VI) to form a Cr-Fe precipitate. Evidence of zones with different levels of Cr(VI) reduction suggest that local heterogeneity may have confounded interpretation of processes based on bulk measurements. This study indicates that the bulk redox status and biogeochemical regime, as categorized by the dominant electron-accepting process, do not necessarily control the final product of Cr(VI) reduction.

Pore-scale controls on calcite dissolution rates from flow-through laboratory and numerical experiments.

A combination of experimental, imaging, and modeling techniques were applied to investigate the pore-scale transport and surface reaction controls on calcite dissolution under elevated pCO2 conditions. The laboratory experiment consisted of the injection of a solution at 4 bar pCO2 into a capillary tube packed with crushed calcite. A high resolution pore-scale numerical model was used to simulate the experiment based on a computational domain consisting of reactive calcite, pore space, and the capillary wall constructed from volumetric X-ray microtomography images. Simulated pore-scale effluent concentrations were higher than those measured by a factor of 1.8, with the largest component of the discrepancy related to uncertainties in the reaction rate model and its parameters. However, part of the discrepancy was apparently due to mass transport limitations to reactive surfaces, which were most pronounced near the inlet where larger diffusive boundary layers formed around grains and in slow-flowing pore spaces that exchanged mass by diffusion with fast flow paths. Although minor, the difference between pore- and continuum-scale results due to transport controls was discernible with the highly accurate methods employed and is expected to be more significant where heterogeneity is greater, as in natural subsurface materials.

Upscaling calcite growth rates from the mesoscale to the macroscale.

Quantitative prediction of mineral reaction rates in the subsurface remains a daunting task partly because a key parameter for macroscopic models, the reactive site density, is poorly constrained. Here we report atomic force microscopy (AFM) measurements on the {1014} calcite surface of monomolecular step densities, treated as equivalent to the reactive site density, as a function of aqueous calcium-to-carbonate ratio and saturation index. Data for the obtuse step orientation are combined with existing step velocity measurements to generate a model that predicts overall macroscopic calcite growth rates. The model is quantitatively consistent with several published macroscopic rates under a range of alkaline solution conditions, particularly for two of the most comprehensive data sets, without the need for additional fit parameters. The model reproduces peak growth rates, and its functional form is simple enough to be incorporated into reactive transport or other macroscopic models designed for predictions in porous media. However, it currently cannot model equilibrium or pH effects and it may overestimate rates at high aqueous calcium-to-carbonate ratios. The discrepancies in rates at high calcium-to-carbonate ratios may be due to differences in pretreatment, such as exposing the seed material to SI ≥ 1.0 to generate/develop growth hillocks, or other factors.

Biostimulation induces syntrophic interactions that impact C, S and N cycling in a sediment microbial community.

Stimulation of subsurface microorganisms to induce reductive immobilization of metals is a promising approach for bioremediation, yet the overall microbial community response is typically poorly understood. Here we used proteogenomics to test the hypothesis that excess input of acetate activates complex community functioning and syntrophic interactions among autotrophs and heterotrophs. A flow-through sediment column was incubated in a groundwater well of an acetate-amended aquifer and recovered during microbial sulfate reduction. De novo reconstruction of community sequences yielded near-complete genomes of Desulfobacter (Deltaproteobacteria), Sulfurovum- and Sulfurimonas-like Epsilonproteobacteria and Bacteroidetes. Partial genomes were obtained for Clostridiales (Firmicutes) and Desulfuromonadales-like Deltaproteobacteria. The majority of proteins identified by mass spectrometry corresponded to Desulfobacter-like species, and demonstrate the role of this organism in sulfate reduction (Dsr and APS), nitrogen fixation and acetate oxidation to CO2 during amendment. Results indicate less abundant Desulfuromonadales, and possibly Bacteroidetes, also actively contributed to CO2 production via the tricarboxylic acid (TCA) cycle. Proteomic data indicate that sulfide was partially re-oxidized by Epsilonproteobacteria through nitrate-dependent sulfide oxidation (using Nap, Nir, Nos, SQR and Sox), with CO2 fixed using the reverse TCA cycle. We infer that high acetate concentrations, aimed at stimulating anaerobic heterotrophy, led to the co-enrichment of, and carbon fixation in Epsilonproteobacteria. Results give an insight into ecosystem behavior following addition of simple organic carbon to the subsurface, and demonstrate a range of biological processes and community interactions were stimulated.

Timing the onset of sulfate reduction over multiple subsurface acetate amendments by measurement and modeling of sulfur isotope fractionation.

Stable isotope fractionations of sulfur are reported for three consecutive years of acetate-enabled uranium bioremediation at the US Department of Energy's Rifle Integrated Field Research Challenge (IFRC) site. The data show a previously undocumented decrease in the time between acetate addition and the onset of sulfate reducing conditions over subsequent amendments, from 20 days in the 2007 experiment to 4 days in the 2009 experiment. Increased sulfide concentrations were observed at the same time as δ(34)S of sulfate enrichment in the first year, but in subsequent years elevated sulfide was detected up to 15 days after increased δ(34)S of sulfate. A biogeochemical reactive transport model is developed which explicitly incorporates the stable isotopes of sulfur to simulate fractionation during the 2007 and 2008 amendments. A model based on an initially low, uniformly distributed population of sulfate reducing bacteria that grow and become spatially variable with time reproduces measured trends in solute concentration and δ(34)S, capturing the change in onset of sulfate reduction in subsequent years. Our results demonstrate a previously unrecognized hysteretic effect in the spatial distribution of biomass growth during stimulated subsurface bioremediation.

Physicochemical heterogeneity controls on uranium bioreduction rates at the field scale.

It has been demonstrated in laboratory systems that U(VI) can be reduced to immobile U(IV) by bacteria in natural environments. The ultimate efficacy of bioreduction at the field scale, however, is often challenging to quantify and depends on site characteristics. In this work, uranium bioreduction rates at the field scale are quantified, for the first time, using an integrated approach. The approach combines field data, inverse and forward hydrological and reactive transport modeling, and quantification of reduction rates at different spatial scales. The approach is used to explore the impact of local scale (tens of centimeters) parameters and processes on field scale (tens of meters) system responses to biostimulation treatments and the controls of physicochemical heterogeneity on bioreduction rates. Using the biostimulation experiments at the Department of Energy Old Rifle site, our results show that the spatial distribution of hydraulic conductivity and solid phase mineral (Fe(III)) play a critical role in determining the field-scale bioreduction rates. Due to the dependence on Fe-reducing bacteria, field-scale U(VI) bioreduction rates were found to be largely controlled by the abundance of Fe(III) minerals at the vicinity of the injection wells and by the presence of preferential flow paths connecting injection wells to down gradient Fe(III) abundant areas.

Strontium and cesium release mechanisms during unsaturated flow through waste-weathered Hanford sediments.

Leaching behavior of Sr and Cs in the vadose zone of Hanford site (Washington) was studied with laboratory-weathered sediments mimicking realistic conditions beneath the leaking radioactive waste storage tanks. Unsaturated column leaching experiments were conducted using background Hanford pore water focused on first 200 pore volumes. The weathered sediments were prepared by 6 months reaction with a synthetic Hanford tank waste leachate containing Sr and Cs (10(-5) and 10(-3) molal representative of LO- and HI-sediment, respectively) as surrogates for (90)Sr and (137)Cs. The mineral composition of the weathered sediments showed that zeolite (chabazite-type) and feldspathoid (sodalite-type) were the major byproducts but different contents depending on the weathering conditions. Reactive transport modeling indicated that Cs leaching was controlled by ion-exchange, while Sr release was affected primarily by dissolution of the secondary minerals. The later release of K, Al, and Si from the HI-column indicated the additional dissolution of a more crystalline mineral (cancrinite-type). A two-site ion-exchange model successfully simulated the Cs release from the LO-column. However, a three-site ion-exchange model was needed for the HI-column. The study implied that the weathering conditions greatly impact the speciation of the secondary minerals and leaching behavior of sequestrated Sr and Cs.

Kinetics of Fe(II)-catalyzed transformation of 6-line ferrihydrite under anaerobic flow conditions.

The readsorption of ferrous ions produced by the abiotic and microbially mediated reductive dissolution of iron oxy-hydroxides drives a series of transformations of the host minerals. To further understand the mechanisms by which these transformations occur and their kinetics within a microporous flow environment, flow-through experiments were conducted in which capillary tubes packed with ferrihydrite-coated glass spheres were injected with inorganic Fe(II) solutions under circumneutral pH conditions at 25 degrees C. Synchrotron X-ray diffraction was used to identify the secondary phase(s) formed and to provide data for quantitative kinetic analysis. At concentrations at and above 1.8 mM Fe(II) in the injection solution, magnetite was the only secondary phase formed (no intermediates were detected), with complete transformation following a nonlinear rate law requiring 28 and 150 h of reaction at 18 and 1.8 mM Fe(II), respectively. However, when the injection solution consisted of 0.36 mM Fe(II), goethite was the predominant reaction product and formed much more slowly according to a linear rate law, while only minor magnetite was formed. When the rates are normalized based on the time to react half of the ferrihydrite on a reduced time plot, it is apparent that the 1.8 mM and 18 mM input Fe(II) experiments can be described by the same reaction mechanism, while the 0.36 input Fe(II) experiment is distinct. The analysis of the transformation kinetics suggests that the transformations involved an electron transfer reaction between the aqueous as well as sorbed Fe(II) and ferrihydrite acting as a semiconductor, rather than a simple dissolution and recrystallization mechanism. A transformation mechanism involving sorbed inner sphere Fe(II) alone is not supported, since the essentially equal coverage of sorption sites in the 18 mM and 1.8 mM Fe(II) injections cannot explain the difference in the transformation rates observed.

Contaminant desorption during long-term leaching of hydroxide-weathered Hanford sediments.

Mineral sorption/coprecipitation is thought to be a principal sequestration mechanism for radioactive (90)Sr and (137)Cs in sediments impacted by hyperalkaline, high-level radioactive waste (HLRW) at the DOE's Hanford site. However, the long-term persistence of neo-formed, contaminant bearing phases after removal of the HLRW source is unknown. We subjected pristine Hanford sediments to hyperalkaline Na-AI-NO(3)-OH solutions containing Sr, Cs, and I at 10(-5), 10(-5), and 10(-7) molal, respectively, for 182 days with either <10 ppmv or 385 ppmv pCO(2). This resulted in the formation of feldspathoid minerals. We leached these weathered sediments with dilute, neutral-pH solutions. After 500 pore volumes (PVs), effluent Sr, Cs, NO(3), Al, Si, and pH reached a steady-state with concentrations elevated above those of feedwater. Reactive transport modeling suggests that even after 500 PV, Cs desorption can be explained by ion exchange reactions, whereas Sr desorption is best described by dissolution of Sr-substituted, neo-formed minerals. While, pCO(2) had no effect on Sr or Cs sorption, sediments weathered at <10 ppmv pCO(2) did desorb more Sr (66% vs 28%) and Cs (13% vs 8%) during leaching than those weathered at 385 ppmv pCO(2). Thus, the dissolution of neo-formed aluminosilicates may represent a long-term, low-level supply of (90)Sr at the Hanford site.

Effects of physical and geochemical heterogeneities on mineral transformation and biomass accumulation during biostimulation experiments at Rifle, Colorado.

Electron donor amendment for bioremediation often results in precipitation of secondary minerals and the growth of biomass, both of which can potentially change flow paths and the efficacy of bioremediation. Quantitative estimation of precipitate and biomass distribution has remained challenging, partly due to the intrinsic heterogeneities of natural porous media and the scarcity of field data. In this work, we examine the effects of physical and geochemical heterogeneities on the spatial distributions of mineral precipitates and biomass accumulated during a biostimulation field experiment near Rifle, Colorado. Field bromide breakthrough data were used to infer a heterogeneous distribution of hydraulic conductivity through inverse transport modeling, while the solid phase Fe(III) content was determined by assuming a negative correlation with hydraulic conductivity. Validated by field aqueous geochemical data, reactive transport modeling was used to explicitly keep track of the growth of the biomass and to estimate the spatial distribution of precipitates and biomass. The results show that the maximum mineral precipitation and biomass accumulation occurs in the vicinity of the injection wells, occupying up to 5.4vol.% of the pore space, and is dominated by reaction products of sulfate reduction. Accumulation near the injection wells is not strongly affected by heterogeneities present in the system due to the ubiquitous presence of sulfate in the groundwater. However, accumulation in the down-gradient regions is dominated by the iron-reducing reaction products, whose spatial patterns are strongly controlled by both physical and geochemical heterogeneities. Heterogeneities can lead to localized large accumulation of mineral precipitates and biomass, increasing the possibility of pore clogging. Although ignoring the heterogeneities of the system can lead to adequate prediction of the average behavior of sulfate-reducing related products, it can also lead to an overestimation of the overall accumulation of iron-reducing bacteria, as well as the rate and extent of iron reduction. Surprisingly, the model predicts that the total amount of uranium being reduced in the heterogeneous 2D system was similar to that in the 1D homogeneous system, suggesting that the overall uranium bioremediation efficacy may not be significantly affected by the heterogeneities of Fe(III) content in the down-gradient regions. Rather, the characteristics close to the vicinity of the injection wells might be crucial in determining the overall efficacy of uranium bioremediation. These findings have important implications not only for uranium bioremediation at the Rifle site and for bioremediation of other redox sensitive contaminants at sites with similar characteristics, but also for the development of optimal amendment delivery strategies in other settings.

Mineral transformation and biomass accumulation associated with uranium bioremediation at Rifle, Colorado.

Injection of organic carbon into the subsurface as an electron donor for bioremediation of redox-sensitive contaminants like uranium often leads to mineral transformation and biomass accumulation, both of which can alter the flow field and potentially bioremediation efficacy. This work combines reactive transport modeling with a column experiment and field measurements to understand the biogeochemical processes and to quantify the biomass and mineral transformation/accumulation during a bioremediation experiment at a uranium contaminated site near Rifle, Colorado. We use the reactive transport model CrunchFlow to explicitly simulate microbial community dynamics of iron and sulfate reducers, and their impacts on reaction rates. The column experiment shows clear evidence of mineral precipitation, primarily in the form of calcite and iron monosulfide. At the field scale, reactive transport simulations suggest that the biogeochemical reactions occur mostly close to the injection wells where acetate concentrations are highest, with mineral precipitate and biomass accumulation reaching as high as 1.5% of the pore space. This work shows that reactive transport modeling coupled with field data can bean effective tool for quantitative estimation of mineral transformation and biomass accumulation, thus improving the design of bioremediation strategies.

Spatially resolved U(VI) partitioning and speciation: implications for plume scale behavior of contaminant U in the Hanford vadose zone.

A saline-alkaline brine containing high concentration of U(VI) was accidentally spilled at the Hanford Site in 1951, introducing 10 tons of U into sediments under storage tank BX-102. U concentrations in the deep vadose zone and groundwater plumes increase with time, yet how the U has been migrating is not fully understood. We simulated the spill event in laboratory soil columns, followed by aging, and obtained spatially resolved U partitioning and speciation along simulated plumes. We found after aging, at apparent steady state, that the pore aqueous phase U concentrations remained surprisingly high (up to 0.022 M), in close agreement with the recently reported high U concentrations (up to 0.027 M) in the vadose zone plume (1). The pH values of aged pore liquids varying from 10 to 7, consistent with the measured pH of the field borehole sediments varying from 9.5 to 7.4 (2), from near the plume source to the plume front. The direct measurements of aged pore liquids together with thermodynamic calculations using a Pitzer approach revealed that UO2(CO3)3(4-) is the dominant aqueous U species within the plume body (pH 8-10), whereas Ca2UO2(CO3)3 and CaUO2(CO3)32- are also significant in the plume frontvicinity (pH 7-8), consistent with that measured from field borehole pore-waters (3). U solid phase speciation varies at different locations along the plume flow path and even within single sediment grains, because of location dependent pore and micropore solution chemistry. Our results suggest that continuous gravity-driven migration of the highly stable U02(CO3)34 in the residual carbonate and sodium rich tank waste solution is likely responsible for the detected growing U concentrations in the vadose zone and groundwater.