Solid phase evolution in the Biosphere 2 hillslope experiment as predicted by modeling of hydrologic and geochemical fluxes

A reactive transport geochemical modeling study was conducted to help predict the mineral transformations occurring over a ten year time-scale that are expected to impact soil hydraulic properties in the Biosphere 2 (B2) synthetic hillslope experiment. The modeling sought to predict the rate and extent of weathering of a granular basalt (selected for hillslope construction) as a function of climatic drivers, and to assess the feedback effects of such weathering processes on the hydraulic properties of the hillslope. Flow vectors were imported from HYDRUS into a reactive transport code, CrunchFlow2007, which was then used to model mineral weathering coupled to reactive solute transport. Associated particle size evolution was translated into changes in saturated hydraulic conductivity using Rosetta software. We found that flow characteristics, including velocity and saturation, strongly influenced the predicted extent of incongruent mineral weathering and neo-phase precipitation on the hillslope. Results were also highly sensitive to specific surface areas of the soil media, consistent with surface reaction controls on dissolution. Effects of fluid flow on weathering resulted in significant differences in the prediction of soil particle size distributions, which should feedback to alter hillslope hydraulic conductivities

Biogeochemical Coupling of Fe and Tc Speciation in Subsurface Sediments: Implications to Long-Term Tc Immobilization

The project has been focused on biochemical processes in subsurface sediments involving Fe that control the valence state, solubility, and effective mobility of 99Tc. Our goal has been to understand the Tc biogeochemistry as it may occur in suboxic and biostimulated subsurface environments. Two objectives have been pursued: (1) To determine the relative reaction rates of 99Tc(VII)O2(aq) with metal reducing bacteria and biogenic Fe(II); and to characterize the identity, structure, and molecular speciation of Tc(IV) products formed through reaction with both biotic and abiotic reductants. (2) To quantify the biogeochemical factors controlling the reaction rate of O2 with Tc(IV)O2?nH2O in sediment resulting from the direct enzymatic reduction of Tc(VII) by DIRB and/or the reaction of Tc(VII) with the various types of biogenic Fe(II) produced by DIRB

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

Toward catchment hydro-biogeochemical theories

Headwater catchments are the fundamental units that connect the land to the ocean. Hydrological flow and biogeochemical processes are intricately coupled, yet their respective sciences have progressed without much integration. Reaction kinetic theories that prescribe rate dependence on environmental variables (e.g., temperature and water content) have advanced substantially, mostly in well-mixed reactors, columns, and warming experiments without considering the characteristics of hydrological flow at the catchment scale. These theories have shown significant divergence from observations in natural systems. On the other hand, hydrological theories, including transit time theory, have progressed substantially yet have not been incorporated into understanding reactions at the catchment scale. Here we advocate for the development of integrated hydro-biogeochemical theories across gradients of climate, vegetation, and geology conditions. The lack of such theories presents barriers for understanding mechanisms and forecasting the future of the Critical Zone under human- and climate-induced perturbations. Although integration has started and co-located measurements are well under way, tremendous challenges remain. In particular, even in this era of "big data," we are still limited by data and will need to (1) intensify measurements beyond river channels and characterize the vertical connectivity and broadly the shallow and deep subsurface; (2) expand to older water dating beyond the time scales reflected in stable water isotopes; (3) combine the use of reactive solutes, nonreactive tracers, and isotopes; and (4) augment measurements in environments that are undergoing rapid changes. To develop integrated theories, it is essential to (1) engage models at all stages to develop model-informed data collection strategies and to maximize data usage; (2) adopt a "simple but not simplistic," or fit-for-purpose approach to include essential processes in process-based models; (3) blend the use of process-based and data-driven models in the framework of "theory-guided data science." Within the framework of hypothesis testing, model-data fusion can advance integrated theories that mechanistically link catchments' internal structures and external drivers to their functioning. It can not only advance the field of hydro-biogeochemistry, but also enable hind- and fore-casting and serve the society at large. Broadly, future education will need to cultivate thinkers at the intersections of traditional disciplines with hollistic approaches for understanding interacting processes in complex earth systems.This article is categorized under:Science of Water > Method

Modeling acid-gas generation from boiling chloride brines

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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°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.</p> <p>Conclusion</p> <p>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°C).</p

The Center for Environmental Kinetics Analysis: an NSF- and DOE-funded Environmental Molecular Science Institute (EMSI) at Penn State

Physicochemical and microbiological processes taking place at environmental interfaces influence natural processes as well as the transport and fate of environmental contaminants, the remediation of toxic chemicals, and the sequestration of anthropogenic CO2. A team of scientists and engineers has been assembled to develop and apply new experimental and computational techniques to expand our knowledge of environmental kinetics. We are also training a cohort of talented and diverse students to work on these complex problems at multiple length scales and to compile and synthesize the kinetic data. Development of the human resources capable of translating molecular-scale information into parameters that are applicable in real world, field-scale problems of environmental kinetics is a major and relatively unique objective of the Institute's efforts. The EMSI team is a partnership among 10 faculty at The Pennsylvania State University (funded by the National Science Foundation Divisions of Chemistry and Earth Sciences), one faculty member at Juniata College, one faculty member at the University of Florida, and four researchers drawn from Los Alamos National Laboratory, Pacific Northwest National Laboratory, and Lawrence Berkeley National Laboratory (funded by the Department of Energy Division of Environmental Remediation Sciences). Interactions among the applied and academic scientists drives research approaches aimed toward solving important problems of national interest. The Institute is organized into three interest groups (IGs) focusing on the processes of dissolution (DIG), precipitation (PIG), and microbial reactions at surfaces (BIG). Some of the research activity from each IG is highlighted to the right. The IGs interact with each other as each interest group studies reactions across the molecular, microscopic, mesoscopic and, in most cases, field scales. For example, abiotic dissolution and precipitation reactions of Fe oxides as studied in the Dissolution IG provides the baseline for kinetic behavior as the BIG researches the interaction of microorganisms with these same minerals. The attachment of bacteria and redox chemistry that occurs between microorganisms and minerals are critical factors in maintaining groundwater quality and remediation of many toxic waste sites and is one of the main thrusts of research within our EMSI. The IGs also participate in using visualization tools to promote greater understanding of complex environmental data. As a whole, CEKA is also working to compile environmental kinetics data into a cyberinfrastructure and database. The database can be accessed at: http://keystone.ist.psu.edu/

A Unified Approach for Process-Based Hydrologic Modeling: 2. Model Implementation and Case Studies

This work advances a unified approach to process-based hydrologic modeling, which we term the “Structure for Unifying Multiple Modeling Alternatives (SUMMA).” The modeling framework, introduced in the companion paper, uses a general set of conservation equations with flexibility in the choice of process parameterizations (closure relationships) and spatial architecture. This second paper specifies the model equations and their spatial approximations, describes the hydrologic and biophysical process parameterizations currently supported within the framework, and illustrates how the framework can be used in conjunction with multivariate observations to identify model improvements and future research and data needs. The case studies illustrate the use of SUMMA to select among competing modeling approaches based on both observed data and theoretical considerations. Specific examples of preferable modeling approaches include the use of physiological methods to estimate stomatal resistance, careful specification of the shape of the within-canopy and below-canopy wind profile, explicitly accounting for dust concentrations within the snowpack, and explicitly representing distributed lateral flow processes. Results also demonstrate that changes in parameter values can make as much or more difference to the model predictions than changes in the process representation. This emphasizes that improvements in model fidelity require a sagacious choice of both process parameterizations and model parameters. In conclusion, we envisage that SUMMA can facilitate ongoing model development efforts, the diagnosis and correction of model structural errors, and improved characterization of model uncertainty