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.

Can an Off-Nominal Landing by an MMRTG-Powered Spacecraft Induce a Special Region on Mars When No Ice Is Present?

This work aims at addressing whether a catastrophic failure of an entry, descent, and landing event of a Multimission Radioisotope Thermoelectric Generator-based lander could embed the heat sources into the martian subsurface and create a local environment that (1) would temporarily satisfy the conditions for a martian Special Region and (2) could establish a transport mechanism through which introduced terrestrial organisms could be mobilized to naturally occurring Special Regions elsewhere on Mars. Two models were run, a primary model by researchers at the Lawrence Berkeley National Laboratory and a secondary model by researchers at the Jet Propulsion Laboratory, both of which were based on selected starting conditions for various surface composition cases that establish the worst-case scenario, including geological data collected by the Mars Science Laboratory at Gale Crater. The summary outputs of both modeling efforts showed similar results: that the introduction of the modeled heat source could temporarily create the conditions established for a Special Region, but that there would be no transport mechanism by which an introduced terrestrial microbe, even if it was active during the temporarily induced Special Region conditions, could be transported to a naturally occurring Special Region of Mars.

Comparison of Electrostatic and Non-Electrostatic Models for U(VI) Sorption on Aquifer Sediments.

A non-electrostatic generalized composite surface complexation model (SCM) was developed for U(VI) sorption on contaminated F-Area sediments from the U.S. Department of Energy Savannah River Site, South Carolina. The objective of this study was to test if a simpler, semi-empirical, non-electrostatic U(VI) sorption model (NEM) could achieve the same predictive performance as a SCM with electrostatic correction terms in describing U(VI) plume evolution and long-term mobility. One-dimensional reactive transport simulations considering key hydrodynamic processes, Al and Fe minerals, as well as H+ and U surface complexation, with and without electrostatic correction terms, were conducted. The NEM was first calibrated with laboratory batch H+ and U(VI) sorption data on F-Area sediments, and then the surface area of the NEM was adjusted to match field observations of dissolved U(VI). Modeling results indicate that the calibrated NEM was able to perform as well as the previously developed electrostatic model in predicting the long-term evolution of H+ and U(VI) at the site, given the variability of field-site data. The electrostatic and NEM models yield somewhat different results for the time period when basin discharge was active; however, it is not clear which modeling approach may be better to model this early time period because groundwater quality data during this period were not available. A key finding of this study is that the applicability of NEM (and thus robustness of its predictions) to the field system evolves with time and is strongly dependent on the pH range that was used to develop the model.

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.

Impact of different environmental conditions on the aggregation of biogenic U(IV) nanoparticles synthesized by Desulfovibrio alaskensis G20.

This study investigates the impact of specific environmental conditions on the formation of colloidal U(IV) nanoparticles by the sulfate reducing bacteria (SRB, Desulfovibrio alaskensis G20). The reduction of soluble U(VI) to less soluble U(IV) was quantitatively investigated under growth and non-growth conditions in bicarbonate or 1,4-piperazinediethanesulfonic acid (PIPES) buffered environments. The results showed that under non-growth conditions, the majority of the reduced U nanoparticles aggregated and precipitated out of solution. High resolution transmission electron microscopy revealed that only a very small fraction of cells had reduced U precipitates in the periplasmic spaces in the presence of PIPES buffer, whereas in the presence of bicarbonate buffer, reduced U was also observed in the cytoplasm with greater aggregation of biogenic U(IV) particles at higher initial U(VI) concentrations. The same experiments were repeated under growth conditions using two different electron donors (lactate and pyruvate) and three electron acceptors (sulfate, fumarate, and thiosulfate). In contrast to the results of the non-growth experiments, even after 0.2 µm filtration, the majority of biogenic U(IV) remained in the aqueous phase resulting in potentially mobile biogenic U(IV) nanoparticles. Size fractionation results showed that U(IV) aggregates were between 18 and 200 nm in diameter, and thus could be very mobile. The findings of this study are helpful to assess the size and potential mobility of reduced U nanoparticles under different environmental conditions, and would provide insights on their potential impact affecting U(VI) bioremediation efforts at subsurface contaminated sites.

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

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

Effect of dissolved CO2 on a shallow groundwater system: a controlled release field experiment.

Capturing carbon dioxide (CO(2)) emissions from industrial sources and injecting the emissions deep underground in geologic formations is one method being considered to control CO(2) concentrations in the atmosphere. Sequestering CO(2) underground has its own set of environmental risks, including the potential migration of CO(2) out of the storage reservoir and resulting acidification and release of trace constituents in shallow groundwater. A field study involving the controlled release of groundwater containing dissolved CO(2) was initiated to investigate potential groundwater impacts. Dissolution of CO(2) in the groundwater resulted in a sustained and easily detected decrease of ~3 pH units. Several trace constituents, including As and Pb, remained below their respective detections limits and/or at background levels. Other constituents (Ba, Ca, Cr, Sr, Mg, Mn, and Fe) displayed a pulse response, consisting of an initial increase in concentration followed by either a return to background levels or slightly greater than background. This suggests a fast-release mechanism (desorption, exchange, and/or fast dissolution of small finite amounts of metals) concomitant in some cases with a slower release potentially involving different solid phases or mechanisms. Inorganic constituents regulated by the U.S. Environmental Protection Agency remained below their respective maximum contaminant levels throughout the experiment.

Influence of chelating agents on biogenic uraninite reoxidation by Fe(III) (Hydr)oxides.

Microbially mediated reduction of soluble U(VI) to U(IV) with subsequent precipitation of uraninite, UO(2(S)), has been proposed as a method for limiting uranium (U) migration. However, microbially reduced UO(2) may be susceptible to reoxidation by environmental factors, with Fe(III) (hydr)oxides playing a significant role. Little is known about the role that organic compounds such as Fe(III) chelators play in the stability of reduced U. Here, we investigate the impact of citrate, DFB, EDTA, and NTA on biogenic UO(2) reoxidation with ferrihydrite, goethite, and hematite. Experiments were conducted in anaerobic batch systems in PIPES buffer (10 mM, pH 7) with bicarbonate for approximately 80 days. Results showed EDTA accelerated UO(2) reoxidation the most at an initial rate of 9.5 µM day(-1) with ferrihydrite, 8.6 µM day(-1) with goethite, and 8.8 µM day(-1) with hematite. NTA accelerated UO(2) reoxidation with ferrihydrite at a rate of 4.8 µM day(-1); rates were less with goethite and hematite (0.66 and 0.71 µM day(-1), respectively). Citrate increased UO(2) reoxidation with ferrihydrite at a rate of 1.8 µM day(-1), but did not increase the extent of reaction with goethite or hematite, with no reoxidation in this case. In all cases, bicarbonate increased the rate and extent of UO(2) reoxidation with ferrihydrite in the presence and absence of chelators. The highest rate of UO(2) reoxidation occurred when the chelator promoted both UO(2) and Fe(III) (hydr)oxide dissolution as demonstrated with EDTA. When UO(2) dissolution did not occur, UO(2) reoxidation likely proceeded through an aqueous Fe(III) intermediate with lower reoxidation rates observed. Reaction modeling suggests that strong Fe(II) chelators promote reoxidation whereas strong Fe(III) chelators impede it. These results indicate that chelators found in U contaminated sites may play a significant role in mobilizing U, potentially affecting bioremediation efforts.

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation.

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4+ production during urea hydrolysis were incorporated in the model and captured critical changes in the major metal species. The electrical phase increases were potentially due to ion exchange processes that modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to understand complex biogeochemical processes in the subsurface.

The toxicity of lead to Desulfovibrio desulfuricans G20 in the presence of goethite and quartz.

An aqueous mixture of goethite, quartz, and lead chloride (PbCl(2)) was treated with the sulfate-reducing bacterium, Desulfovibrio desulfuricans G20 (D. desulfuricans G20), in a medium specifically designed to assess metal toxicity. In the presence of 26 muM of soluble Pb, together with the goethite and quartz, D. desulfuricans G20 grew after a lag time of 5 days compared to 2 days in Pb-, goethite-, and quartz-free treatments. In the absence of goethite and quartz, however, with 26 microM soluble Pb, no measurable growth was observed. Results showed that D. desulfuricans G20 first removed Pb from solutions then growth began resulting in black precipitates of Pb and iron sulfides. Transmission electron microscopic analyses of thin sections of D. desulfuricans G20 treated with 10 microM PbCl(2) in goethite- and quartz-free treatment showed the presence of a dense deposit of lead sulfide precipitates both in the periplasm and cytoplasm. However, thin sections of D. desulfuricans G20 treated with goethite, quartz, and PbCl(2) (26 microM soluble Pb) showed the presence of a dense deposit of iron sulfide precipitates both in the periplasm and cytoplasm. Energy-dispersive X-ray spectroscopy, selected area electron diffraction patterns, or X-ray diffraction analyses confirmed the structure of precipitated Pb inside the cell as galena (PbS) in goethite- and quartz-free treatments, and iron sulfides in treatments with goethite, quartz, and PbCl(2). Overall results suggest that even at the same soluble Pb concentration (26 microM), in the presence of goethite and quartz, apparent Pb toxicity to D. desulfuricans G20 decreased significantly. Further, accumulation of lead/iron sulfides inside D. desulfuricans G20 cells depended on the presence of goethite and quartz.

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).

Influence of heavy metals on microbial growth kinetics including lag time: mathematical modeling and experimental verification.

Heavy metals can significantly affect the kinetics of substrate biodegradation and microbial growth, including lag times and specific growth rates. A model to describe microbial metabolic lag as a function of the history of substrate concentration has been previously described by Wood et al. (Water Resour Res 31:553-563) and Ginn (Water Resour Res 35:1395-1408). In the present study, this model is extended by including the effect of heavy metals on metabolic lag by developing an inhibitor-dependent functional to account for the metabolic state of the microorganisms. The concentration of the inhibiting metal is explicitly incorporated into the functional. The validity of the model is tested against experimental data on the effects of zinc on Pseudomonas species isolated from Lake Coeur d'Alene sediments, Idaho, U.S.A., as well as the effects of nickel or cobalt on a mixed microbial culture collected from the aeration tank of a wastewater treatment plant in Athens, Greece. The simulations demonstrate the ability to incorporate the effect of metals on metabolism through lag, yield coefficient, and specific growth rates. The model includes growth limitation due to insufficient transfer of oxygen into the growth medium.

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.