The environmental controls on efficiency of enhanced rock weathering in soils.

Enhanced rock weathering (ERW) in soils is a promising carbon removal technology, but the realistically achievable efficiency, controlled primarily by in situ weathering rates of the applied rocks, is highly uncertain. Here we explored the impacts of coupled biogeochemical and transport processes and a set of primary environmental and operational controls, using forsterite as a proxy mineral in soils and a multiphase multi-component reactive transport model considering microbe-mediated reactions. For a onetime forsterite application of ~ 16 kg/m2, complete weathering within five years can be achieved, giving an equivalent carbon removal rate of ~ 2.3 kgCO2/m2/yr. However, the rate is highly variable based on site-specific conditions. We showed that the in situ weathering rate can be enhanced by conditions and operations that maintain high CO2 availability via effective transport of atmospheric CO2 (e.g. in well-drained soils) and/or sufficient biogenic CO2 supply (e.g. stimulated plant-microbe processes). Our results further highlight that the effect of increasing surface area on weathering rate can be significant-so that the energy penalty of reducing the grain size may be justified-only when CO2 supply is nonlimiting. Therefore, for ERW practices to be effective, siting and engineering design (e.g. optimal grain size) need to be co-optimized.

Microbial U Isotope Fractionation Depends on the U(VI) Reduction Rate.

U isotope fractionation may serve as an accurate proxy for U(VI) reduction in both modern and ancient environments, if the systematic controls on the magnitude of fractionation (ε) are known. We model the effect of U(VI) reduction kinetics on U isotopic fractionation during U(VI) reduction by a novel Shewanella isolate, Shewanella sp. (NR), in batch incubations. The measured ε values range from 0.96 ± 0.16 to 0.36 ± 0.07‰ and are strongly dependent on the U(VI) reduction rate. The ε decreases with increasing reduction rate constants normalized by cell density and initial U(VI). Reactive transport simulations suggest that the rate dependence of ε is due to a two-step process, where diffusive transport of U(VI) from the bulk solution across a boundary layer is followed by enzymatic reduction. Our results imply that the spatial decoupling of bulk U(VI) solution and enzymatic reduction should be taken into account for interpreting U isotope data from the environment.

Reactive Transport Model of Sulfur Cycling as Impacted by Perchlorate and Nitrate Treatments.

Microbial souring in oil reservoirs produces toxic, corrosive hydrogen sulfide through microbial sulfate reduction, often accompanying (sea)water flooding during secondary oil recovery. With data from column experiments as constraints, we developed the first reactive-transport model of a new candidate inhibitor, perchlorate, and compared it with the commonly used inhibitor, nitrate. Our model provided a good fit to the data, which suggest that perchlorate is more effective than nitrate on a per mole of inhibitor basis. Critically, we used our model to gain insight into the underlying competing mechanisms controlling the action of each inhibitor. This analysis suggested that competition by heterotrophic perchlorate reducers and direct inhibition by nitrite produced from heterotrophic nitrate reduction were the most important mechanisms for the perchlorate and nitrate treatments, respectively, in the modeled column experiments. This work demonstrates modeling to be a powerful tool for increasing and testing our understanding of reservoir-souring generation, prevention, and remediation processes, allowing us to incorporate insights derived from laboratory experiments into a framework that can potentially be used to assess risk and design optimal treatment schemes.

Deep-sea oil plume enriches indigenous oil-degrading bacteria.

The biological effects and expected fate of the vast amount of oil in the Gulf of Mexico from the Deepwater Horizon blowout are unknown owing to the depth and magnitude of this event. Here, we report that the dispersed hydrocarbon plume stimulated deep-sea indigenous γ-Proteobacteria that are closely related to known petroleum degraders. Hydrocarbon-degrading genes coincided with the concentration of various oil contaminants. Changes in hydrocarbon composition with distance from the source and incubation experiments with environmental isolates demonstrated faster-than-expected hydrocarbon biodegradation rates at 5°C. Based on these results, the potential exists for intrinsic bioremediation of the oil plume in the deep-water column without substantial oxygen drawdown.

Modeling acid-gas generation from boiling chloride brines.

BACKGROUND: This study investigates the generation of HCl and other acid gases from boiling calcium chloride dominated waters at atmospheric pressure, primarily using numerical modeling. The main focus of this investigation relates to the long-term geologic disposal of nuclear waste at Yucca Mountain, Nevada, where pore waters around waste-emplacement tunnels are expected to undergo boiling and evaporative concentration as a result of the heat released by spent nuclear fuel. Processes that are modeled include boiling of highly concentrated solutions, gas transport, and gas condensation accompanied by the dissociation of acid gases, causing low-pH condensate. RESULTS: Simple calculations are first carried out to evaluate condensate pH as a function of HCl gas fugacity and condensed water fraction for a vapor equilibrated with saturated calcium chloride brine at 50-150 degrees C and 1 bar. The distillation of a calcium-chloride-dominated brine is then simulated with a reactive transport model using a brine composition representative of partially evaporated calcium-rich pore waters at Yucca Mountain. Results show a significant increase in boiling temperature from evaporative concentration, as well as low pH in condensates, particularly for dynamic systems where partial condensation takes place, which result in enrichment of HCl in condensates. These results are in qualitative agreement with experimental data from other studies. CONCLUSION: The combination of reactive transport with multicomponent brine chemistry to study evaporation, boiling, and the potential for acid gas generation at the proposed Yucca Mountain repository is seen as an improvement relative to previously applied simpler batch evaporation models. This approach allows the evaluation of thermal, hydrological, and chemical (THC) processes in a coupled manner, and modeling of settings much more relevant to actual field conditions than the distillation experiment considered. The actual and modeled distillation experiments do not represent expected conditions in an emplacement drift, but nevertheless illustrate the potential for acid-gas generation at moderate temperatures (<150 degrees C).

A reaction-transport model for calcite precipitation and evaluation of infiltration fluxes in unsaturated fractured rock.

The percolation flux in the unsaturated zone (UZ) is an important parameter addressed in site characterization and flow and transport modeling of the potential nuclear-waste repository at Yucca Mountain, NV, USA. The US Geological Survey (USGS) has documented hydrogenic calcite abundances in fractures and lithophysal cavities at Yucca Mountain to provide constraints on percolation fluxes in the UZ. The purpose of this study was to investigate the relationship between percolation flux and measured calcite abundances using reactive transport modeling. Our model considers the following essential factors affecting calcite precipitation: (1) infiltration, (2) the ambient geothermal gradient, (3) gaseous CO(2) diffusive transport and partitioning in liquid and gas phases, (4) fracture-matrix interaction for water flow and chemical constituents, and (5) water-rock interaction. Over a bounding range of 2-20 mm/year infiltration rate, the simulated calcite distributions capture the trend in calcite abundances measured in a deep borehole (WT-24) by the USGS. The calcite is found predominantly in fractures in the welded tuffs, which is also captured by the model simulations. Simulations showed that from about 2 to 6 mm/year, the amount of calcite precipitated in the welded Topopah Spring tuff is sensitive to the infiltration rate. This dependence decreases at higher infiltration rates owing to a modification of the geothermal gradient from the increased percolation flux. The model also confirms the conceptual model for higher percolation fluxes in the fractures compared to the matrix in the welded units, and the significant contribution of Ca from water-rock interaction. This study indicates that reactive transport modeling of calcite deposition can yield important constraints on the unsaturated zone infiltration-percolation flux and provide useful insight into processes such as fracture-matrix interaction as well as conditions and parameters controlling calcite deposition.

Calibration of Yucca Mountain unsaturated zone flow and transport model using porewater chloride data.

In this study, porewater chloride data from Yucca Mountain, NV are analyzed and modeled by three-dimensional chemical transport simulation and analytical methods. The simulation modeling approach is based on a continuum formulation of coupled multiphase fluid flow and tracer transport processes through fractured porous rock using a dual-continuum concept. Infiltration rate calibrations were performed using the porewater chloride data. Model results of chloride distributions were improved in matching the observed data with the calibrated infiltration rates. Statistical analyses of the frequency distribution for overall percolation fluxes and chloride concentration in the unsaturated zone system demonstrate that the use of the calibrated infiltration rates had an insignificant effect on the distribution of simulated percolation fluxes but significantly changed the predicted distribution of simulated chloride concentrations. An analytical method was also applied to model transient chloride transport. The method was verified by three-dimensional simulation results to be capable of capturing major chemical transient behavior and trends. Effects of lateral flow in the Paintbrush nonwelded unit on percolation fluxes and chloride distribution were studied by three-dimensional simulations with increased horizontal permeability. The combined results from these model calibrations furnish important information for the UZ model studies, contributing to performance assessment of the potential repository.

Experimental and numerical simulation of dissolution and precipitation: implications for fracture sealing at Yucca Mountain, Nevada.

Plugging of flow paths caused by mineral precipitation in fractures above the potential repository at Yucca Mountain, Nevada could reduce the probability of water seeping into the repository. As part of an ongoing effort to evaluate thermal-hydrological-chemical (THC) effects on flow in fractured media, we performed a laboratory experiment and numerical simulations to investigate mineral dissolution and precipitation under anticipated temperature and pressure conditions in the repository. To replicate mineral dissolution by vapor condensate in fractured tuff, water was flowed through crushed Yucca Mountain tuff at 94 degrees C. The resulting steady-state fluid composition had a total dissolved solids content of about 140 mg/l; silica was the dominant dissolved constituent. A portion of the steady-state mineralized water was flowed into a vertically oriented planar fracture in a block of welded Topopah Spring Tuff that was maintained at 80 degrees C at the top and 130 degrees C at the bottom. The fracture began to seal with amorphous silica within 5 days.A 1-D plug-flow numerical model was used to simulate mineral dissolution, and a similar model was developed to simulate the flow of mineralized water through a planar fracture, where boiling conditions led to mineral precipitation. Predicted concentrations of the major dissolved constituents for the tuff dissolution were within a factor of 2 of the measured average steady-state compositions. The mineral precipitation simulations predicted the precipitation of amorphous silica at the base of the boiling front, leading to a greater than 50-fold decrease in fracture permeability in 5 days, consistent with the laboratory experiment.These results help validate the use of a numerical model to simulate THC processes at Yucca Mountain. The experiment and simulations indicated that boiling and concomitant precipitation of amorphous silica could cause significant reductions in fracture porosity and permeability on a local scale. However, differences in fluid flow rates and thermal gradients between the experimental setup and anticipated conditions at Yucca Mountain need to be factored into scaling the results of the dissolution/precipitation experiments and associated simulations to THC models for the potential Yucca Mountain repository.