Spectroscopy and Diffraction

Despite the regular occurrence of both magnetite and iron-metabolizing bacteria in the same environments, it is currently unknown whether the iron(II) and iron(III) in magnetite can be cycled between different bacteria and whether or how magnetic properties are affected by this metabolic activity. We show through magnetic and spectroscopic measurements that the phototrophic Fe(II)-oxidizer Rhodopseudomonas palustris TIE-1 can oxidize solid-phase magnetite nanoparticles using light energy, leading to a decrease in the measured magnetic susceptibility (MS). This process likely occurs at the surface and is reversible in the dark by the Fe(III)-reducer Geobacter sulfurreducens resulting in an increase in MS. These results show that iron ions bound in highly crystalline mineral magnetite are bioavailable as electron stores and electron sinks under varying environmental conditions, making magnetite a potential “biogeobattery” during day/night cycles. These findings are relevant for environmental studies and reinforce the impact of microbial redox processes on the global iron cycle.

Subsurface injection of CO2 for enhanced hydrocarbon recovery, hydraulic fracturing of unconventional reservoirs, and geologic carbon sequestration produces a complex geochemical setting in which CO2-dominated fluids containing dissolved water and organic compounds interact with rocks and minerals. The details of these reactions are relatively unknown and benefit from additional experimentally derived data. In this study, we utilized an in situ X-ray diffraction technique to examine the carbonation reactions of forsterite (Mg2SiO4) during exposure to supercritical CO2 (scCO2) that had been equilibrated with aqueous solutions of acetate, oxalate, malonate, or citrate at 50 °C and 90 bar. The organics affected the relative abundances of the crystalline reaction products, nesquehonite (MgCO3·3H2O) and magnesite (MgCO3), likely due to enhanced dehydration of the Mg2+ cations by the organic ligands. These results also indicate that the scCO2 solvated and transported the organic ligands to the forsterite surface. This phenomenon has profound implications for mineral transformations and mass transfer in the upper crust.

Composite Portland cement-basalt caprock cores with fractures, as well as neat Portland cement columns, were prepared to understand the geochemical and geomechanical effects on the integrity of wellbores with defects during geologic carbon sequestration. The samples were reacted with CO2-saturated groundwater at 50 ºC and 10 MPa for 3 months under static conditions, while one cement-basalt core was subjected to mechanical stress at 2.7 MPa before the CO2 reaction. Micro-XRD and SEM-EDS data collected along the cement-basalt interface after 3-month reaction with CO2-saturated groundwater indicate that carbonation of cement matrix was extensive with the precipitation of calcite, aragonite, and vaterite, whereas the alteration of basalt caprock was minor. X-ray microtomography (XMT) provided three-dimensional (3-D) visualization of the opening and interconnection of cement fractures due to mechanical stress. Computational fluid dynamics (CFD) modeling further revealed that this stress led to the increase in fluid flow and hence permeability. After the CO2-reaction, XMT images displayed that calcium carbonate precipitation occurred extensively within the fractures in the cement matrix, but only partially along the fracture located at the cement-basalt interface. The 3-D visualization and CFD modeling also showed that the precipitation of calcium carbonate within the cement fractures after the CO2-reaction resulted in the disconnection of cement fractures and permeability decrease. The permeability calculated based on CFD modeling was in agreement with the experimentally determined permeability. This study demonstrates that XMT imaging coupled with CFD modeling represent a powerful tool to visualize and quantify fracture evolution and permeability change in geologic materials and to predict their behavior during geologic carbon sequestration or hydraulic fracturing for shale gas production and enhanced geothermal systems.

The garnet structure has been proposed as a potential crystalline nuclear waste form for accommodation of actinide elements, especially uranium (U). In this study, yttrium iron garnet (YIG) as a model garnet host was studied for the incorporation of U analogs, cerium (Ce), and thorium (Th), incorporated by a charge-coupled substitution with calci-um (Ca) for yttrium (Y) in YIG, namely 2Y3+ = Ca2+ + M4+, where M4+ = Ce4+ or Th4+. Single phase garnets Y3-xCa0.5xM0.5xFe5O12, synthesized by the citrate-nitrate combustion method, were obtained up to x = 0.7. Ce was confirmed to be tetravalent by X-ray absorption spectroscopy and X-ray photoelectron spectroscopy. X-ray diffraction and 57Fe-Mössbauer spectroscopy indicated that the samples are single phase, M4+ and Ca2+ cations are restricted to the c-site, the nature of M4+ has only a minor effect on the structure, and the local environments of both the tetrahedral and octahedral Fe3+ are systematically affected by the extent of substitution, especially on the tetrahedral sublattice. The charge coupled substitution has advantages in incorporating Ce/Th and in stabilizing the substituted phases, compared to a single substitution strategy. Enthalpies of formation of garnets were obtained by high temperature oxide melt solution calorimetry, and the enthalpies of substitution of Ce and Th were determined. The thermodynamic analysis demonstrates that the substituted garnets are entropically rather than energetically stabilized. This suggests that such garnets may form and persist in repositories at high temperature but might decompose near room temperature. These structural and thermodynamic findings shed light on possible incorporation of U in this garnet system.

Extensive research efforts have been invested in reducing model errors to improve the predictive ability of biogeochemical earth and environmental system simulators, with applications ranging from contaminant transport and remediation to impacts of biogeochemical elemental cycling (e.g., carbon and nitrogen) on local ecosystems and regional to global climate. While the bulk of this research has focused on improving model parameterizations in the face of observational limitations, the more challenging type of model error/uncertainty to identify and quantify is model structural error which arises from incorrect mathematical representations of (or failure to consider) important physical, chemical, or biological processes, properties, or system states in model formulations. While improved process understanding can be achieved through scientific study, such understanding is usually developed at small scales. Process-based numerical models are typically designed for a particular characteristic length and time scale. For application-relevant scales, it is generally necessary to introduce approximations and empirical parameterizations to describe complex systems or processes. This single-scale approach has been the best available to date because of limited understanding of process coupling combined with practical limitations on system characterization and computation. While computational power is increasing significantly and our understanding of biological and environmental processes at fundamental scales is accelerating, using this information to advance our knowledge of the larger system behavior requires the development of multiscale simulators. Accordingly there has been much recent interest in novel multiscale methods in which microscale and macroscale models are explicitly coupled in a single hybrid multiscale simulation. A limited number of hybrid multiscale simulations have been developed for biogeochemical earth systems, but they mostly utilize application-specific and sometimes ad-hoc approaches for model coupling. We are developing a generalized approach to hierarchical model coupling designed for high-performance computational systems, based on the Swift computing workflow framework. In this presentation we will describe the generalized approach and provide two use cases: 1) simulation of a mixing-controlled biogeochemical reaction coupling pore- and continuum-scale models, and 2) simulation of biogeochemical impacts of groundwater – river water interactions coupling fine- and coarse-grid model representations. This generalized framework can be customized for use with any pair of linked models (microscale and macroscale) with minimal intrusiveness to the at-scale simulators. It combines a set of python scripts with the Swift workflow environment to execute a complex multiscale simulation utilizing an approach similar to the well-known Heterogeneous Multiscale Method. User customization is facilitated through user-provided input and output file templates and processing function scripts, and execution within a high-performance computing environment is handled by Swift, such that minimal to no user modification of at-scale codes is required.