Molecular Dynamics Simulations of Water Structure and Diffusion in Silica Nanopores

We present molecular dynamics (MD) simulations of water-filled silica nanopores such as those that occur in ordered oxide ceramics (MCM-41, SBA-15), controlled pore glasses (such as Vycor glass), mesoporous silica, bioglasses, and hydrous silica gel coatings of weathered minerals and glasses. Our simulations overlap the range of pore diameters (1–4 nm) where confinement causes the disappearance of bulk-liquid-like water. In<i> ≥</i>2 nm diameter pores, the silica surface carries three statistical monolayers of density-layered water, interfacial water structure is independent of confinement or surface curvature, and bulk-liquid-like water exists at the center of the pore (this last finding contradicts assumptions used in most previous neutron diffraction studies and in several MD simulation studies of silica nanopores). In 1 nm diameter pores, bulk-liquid-like water does not exist and the structural properties of interfacial water are influenced by confinement. Predicted water diffusion coefficients in 1–4 nm diameter pores agree with quasi-elastic neutron scattering (QENS) data and are roughly consistent with a very simple “core–shell” conceptual model whereupon the first statistical water monolayer is immobile and the rest of the pore water diffuses as rapidly as bulk liquid water

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 {101̅4} 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

Contaminant Desorption during Long-Term Leaching of Hydroxide-Weathered Hanford Sediments

Mineral sorption/coprecipitation is thought to be a principal sequestration mechanism for radioactive 90Sr and 137Cs 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−Al−NO3−OH solutions containing Sr, Cs, and I at 10−5, 10−5, and 10−7 molal, respectively, for 182 days with either pCO2. 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, NO3, 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, pCO2 had no effect on Sr or Cs sorption, sediments weathered at pCO2 did desorb more Sr (66% vs 28%) and Cs (13% vs 8%) during leaching than those weathered at 385 ppmv pCO2. Thus, the dissolution of neo-formed aluminosilicates may represent a long-term, low-level supply of 90Sr at the Hanford site

Contaminant Desorption During Long-Term Leaching of Hydroxide-Weathered Sediments

Contaminant Desorption During Long-Term Leaching of Hydroxide-Weathered Sediment

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 °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

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

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 be an effective tool for quantitative estimation of mineral transformation and biomass accumulation, thus improving the design of bioremediation strategies

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 δ³⁴S of sulfate enrichment in the first year, but in subsequent years elevated sulfide was detected up to 15 days after increased δ³⁴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 δ³⁴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

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 90Sr and 137Cs. 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

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 pCO₂ conditions. The laboratory experiment consisted of the injection of a solution at 4 bar pCO₂ 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