Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗
DOE OSTI2024

Wellbore cement alteration and roles of CO 2 and shale during underground hydrogen storage

To mitigate climate change and adopt renewable energy, energy storage is crucial and can be done in the form of hydrogen gas (H 2 ). Subsurface geologic reservoirs are positioned to store H 2 on the largest scales for the longest terms of all potential options. However, H 2 injection may boost reactions that consume hydrogen, generate undesired gases, and alter pore structures of geomedia. To explore the extent of H 2 -associated biotic reactions at a near wellbore location, four experiments were conducted under underground storage conditions with wellbore cement cores and, in most instances, shale samples submerged in synthetic formation brine. Post-reaction gas, aqueous, and solid phase samples were analyzed using olfactory screening and, later, gas chromatography (GC-MS), inductively coupled plasma optical emission spectroscopy (ICP-OES), scanning electron microscopy (SEM), and synchrotron micro-scale x-ray fluorescence (μ-XRF). Within a period of 16 weeks, hydrogen sulfide (H 2 S) was generated in systems containing both H 2 and shale. XRF mapping identified a zone enriched in iron(II) and reduced sulfur along the rim of cement cross sections that was largely associated with CO 2 -induced cement carbonation. Shale did not show noticeable alteration, but there is evidence it contributed to the initial inoculation of the system and provided nutrients for microbes via water-rock interactions. Here, this study considers both rock formations and wellbore cement not previously evaluated concurrently. Findings support understanding and modeling of H 2 -associated biogeochemical reactions during underground hydrogen storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗
DOE OSTI2025

Selective recovery of native copper from basalt tailings using alkaline glycinate solution

Mafic and ultramafic rocks present an intriguing pathway for CO 2 capture through strategic enhancements to the natural silicate weathering cycle. Simultaneously, these rock types are often hosts to appreciable amounts of metals critical to U.S. energy independence, particularly in the context of alkaline mine tailings from historical metal mining. The research and scientific prerogative, then, is to identify promising mafic and ultramafic feedstocks and conduct exploratory studies to effectively recover these critical minerals while preserving the potential for mineral carbonation. The Keweenaw Peninsula of Michigan hosts the largest native Cu reservoir within basalt in the world and experienced a rich history of Cu production spanning from pre-historic times to the mid-1900s. Today, much of the Cu mining legacy remains as mine waste tailings. In this study, we examined the Cu extractability from Keweenaw Basalt tailings using a sodium glycinate solution, in comparison with acid and sodium hydroxide leaching. Experimental results showed an 85 % Cu extraction rate using sodium glycinate as the extraction solution with negligible release of other cations from the basalt. The kinetic and extraction mechanisms of Cu selective recovery using glycinate solution were discussed using time-resolved experimental data and kinetic geochemical modeling. Theoretical estimation of carbon mineralization potential of all the existing basalt waste tailings (∼500 million tons) can reach 85.5 MMT CO 2 . A total of 0.786 MMT Cu can be recovered with sodium glycinate, with a value of 7.7 billion USD. In conclusion, this novel application of alkaline glycinate for selective Cu recovery from basalt mine tailings demonstrates the viability of selective metal recovery using a non-hazardous chemical while preserving CO 2 capture potential and presents a potential pathway toward reducing energy-related emissions and providing an unconventional domestic source of critical minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Reconstructing porous media using generative flow networks

One area of intense scientific interest for the study of sandstones, carbonates, and shale at the pore scale is the use of limited image and petrophysical data to generate multiple realizations of a rock’s pore structure. Such images aid efforts to quantify uncertainty in petrophysical properties, including porosity–permeability transforms. We develop and evaluate a deep learning-based method to synthesize porous media volumes using a so-called generative flow model trained on x-ray microscope images of rock texture and pore structure. These models are optimized on a log-likelihood objective and they synthesize large and realistic three-dimensional images. Further, we demonstrate the rapid generation of sandstone image volumes that display realism as gauged by quantitative comparison of topological features using Minkowski functionals of porosity, specific surface area, and the Euler–Poincaré characteristic (i.e., pore connectivity). We also evaluate the single-phase permeability using Navier–Stokes and lattice Boltzmann methods and show that transport properties of the generated samples match measured trends.

58 GEOSCIENCES↗
DOE OSTI2022

Chemical and isotopic evolution of flowback fluids from the Utica Gas Shale Play, Eastern Ohio USA

For this work, hydraulic fracturing flowback fluids were collected from two Utica/Point Pleasant well pads in eastern Ohio. One site was in the wet gas zone (UPPW), while the site ~50 km south consisted of four wells on the same pad in the dry gas zone (UPPS). Samples of input fluids also were collected before and during hydraulic fracturing. Flowback fluids are Na-Ca-Cl brines with total dissolved salt (TDS) concentrations that increase over several months from ~100 to 200 g/L. The slightly higher TDS of the dry gas fluids are in part due to recycled flowback used as input fluids, and in part due to lower volume of water used in hydraulic fracturing. Concentrations of most major ions (Ca 2+ , Mg 2+ , Na + , Sr2 + , Fe, Mn, Cl – , Br – ) are similar for the five wells sampled, although small but systematic changes occur in the major element ratios over time. Most notably, an increase in the Sr/Cl ratio corresponds to an increase in the 87 Sr/ 86 Sr ratio in the fluids, suggesting interactions with a more radiogenic Sr source in the subsurface. Dissolved Ba concentrations and Ra activities were different between the two sites, reflecting a high SO 4 2– fluid used at the UPPW site, and water-rock reactions occurring during hydraulic fracturing at the UPPS site. Water oxygen (δ 18 O) and hydrogen (δD) isotopes for input fluids used in the UPPW4 well and fresh water used for the UPPS wells fall on the Global Meteoric waterline (GMWL). Flowback fluids from both sites are relatively enriched in δ 18 O and δD compared to the input, but do not appear to follow a simple mixing trend, suggesting reaction and isotopic exchange with carbonates and fractionation due to imbibition in the rock. Cl isotopes, δ 37 Cl in the FP fluids varied from ~ –0.43 to +0.13‰, the largest variation was observed in the earlier stages of flowback, while in the later stages δ 37 Cl exhibited a small but systematic increase over time, suggesting diffusion control of isotopic composition. Some of the trace species measured (dissolved Fe, Mn, thiosulfate and organic acids) do not follow the same trends as the major ions, suggesting contributions from input fluids, or microbially mediated reactions are exerting control on their concentrations. Input water chemistry exerts an important control on the Sr and Ba concentrations in flowback water. High SO 4 2– in the input fluids used for hydraulically fracturing the wet gas well leads to precipitation of barite-celestite in flowback fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Evaluation of coal-associated sediments, wastes, and AMD sludge in the Southern Appalachian Basin as feedstock materials for REE and Li recovery

Critical minerals (CM) such as rare earth elements (REE+) and Lithium (Li) are essential to technological innovation, energy transitions, global economic and defense security, necessitating the search for unconventional resources and efficient recovery methods to avert supply chain disruptions. Here, this study evaluates coal-associated sediments (underclay and roof rock) and wastes from the Pennsylvanian Pottsville Formation of the Southern Appalachian Basin (SAB) as potential feedstocks for CM recovery. A total of 34 samples (15 underclays, 12 roof rocks, 5 Acid Mine Drainage (AMD) sludges, and 2 coal mining wastes) were characterized using XRD, XRF, ICP-MS, and μ-XRF analytical methods. The REE+ and Li concentrations of these materials ranged from 46.8 to 334.4 ppm and from 11.1 to 519 ppm, respectively, with one underclay sample (Hendrix 3456) yielding the highest values for both. Bulk mineralogy for all samples was dominated by aluminosilicate clay phases, particularly illite and kaolinite. All samples exhibited REY def, rel% values >26% and C outl indices that ranged from 0.69 to 0.94, classifying their REE ore potential as Category II (Promising) as defined by Seredin and Dai (2012). Extractability tests (EPA method 3051 A) yielded low REE+ and Li recoveries, with maximum values of 3.3% and 3.6%, respectively, suggesting associations with resistant minerals like clay and phosphates. Elemental mapping indicates that REE+ is associated with phosphate, whereas statistical analysis suggests that REE+ are associated with aluminosilicates, suggesting heterogeneous associations or minimal phosphate contribution. Li also correlated positively with Al 2 O 3 , indicating an aluminosilicate host. This study highlights the potential of coal-associated sediments in the SAB.

Clay minerals↗
DOE OSTI2020

Hydro-mechanical coupled modeling of hydraulic fracturing using the mesh fragmentation technique

The inclusion of special interface elements between standard finite elements is becoming a common technique to model discontinuities in porous media using the finite element method (FEM). Recently, it was proved that high aspect ratio (HAR) elements can reproduce very satisfactorily the main features of behavior observed in a single hydraulically-induced fracture. In this paper, the mesh fragmentation technique (MFT) is proposed to extend the application of HAR elements to simulate more general problems involving the formation of hydraulic fractures in rocks by introducing this type of element in-between the standard FE of a typical mesh. The proposed framework makes use of standard FE techniques. Therefore, it is relatively easy to upgrade an existing FE program for continuous porous media to deal with evolving discontinuities. The proposed approach is verified and validated against available analytical and numerical solutions. The crack patterns predicted by the MFT are also compared with the ones observed in the laboratory. The formation of multiple fractures is also analyzed. In all cases, the proposed technique was able to properly capture the main features of rock behavior subjected to hydraulic fracturing.

58 GEOSCIENCES↗
DOE OSTI2020

Estimation of stress and stress-induced permeability change in a geological nuclear waste repository in a thermo-hydrologically coupled simulation

Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed here. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Climatic controls on erosional efficiency vary with lithology across the Himalaya

Although climate can strongly influence erosional efficiency (i.e., erosion rate for a given topography), demonstrating its impact in tectonically active areas has been challenging due to other confounding controlling factors, such as lithology. Here, we show that 10 Be-derived erosion rates and efficiencies in the Himalayan orogen exhibit distinct relationships with climatic factors depending on lithology. We compile 173 10 Be-derived, basin-averaged erosion rates across the orogen, including 12 newly measured rates from the Dibang and Lohit valleys in the easternmost Himalaya, regions characterized by high precipitation magnitudes and variability. We group basins based on lithologies separated by orogen-scale thrust faults and quantify erosional efficiency coefficients based on the relationships between erosion rates and topographic metrics. Our results show that erosion rates and erosional efficiency from sedimentary and metasedimentary rocks along the Himalayan range front display a positive, nonlinear correlation with climatic factors, such as the number of extreme rainfall events and mean annual precipitation rates. In contrast, erosion rates from crystalline lithologies in the hanging wall of the Main Central thrust show a strong correlation with fluvial topography, whereas erosional efficiency shows no statistically significant correlation with climatic factors. Rapid erosion rates and high erosional efficiencies in the eastern Himalayan range front are likely driven by extreme precipitation on tectonically active, steep slopes composed of mechanically weak metasedimentary rocks. Our findings highlight the importance of the interplay between controlling factors, which include tectonics, lithology, and climate, that drive surface erosion and influence the topographic evolution of orogenic systems.

Channel steepness↗
DOE OSTI2021

Geochemical reactions and alteration of pore architecture in saturated shale after injection of stimulation fluid

Pore architecture regulates fluid flow between unconventional shale reservoir and hydraulically-induced fractures. Imbibition of stimulation fluid may change this architecture and alter hydrocarbon flow. Hydrothermal experiments were conducted at reservoir conditions (125 °C, 45 MPa) to test two hypotheses: 1) Shale, not stimulation fluid, dominates the geochemistry of an unconventional reservoir containing formation water; and 2) Mineral dissolution and precipitation induced by stimulation fluid is transient and manifests across micro-, meso-, and macro-scales of pore architecture. Experiments replicated a shut-in well in the Cretaceous Baxter Shale, Green River Basin, Wyoming USA. Stimulation fluid was injected into one experiment after formation water (I = 0.35 mol/kg, pH = 6.4) and core reacted for 48 days. This novel approach equilibrated rock and formation water and saturated pores with formation water before introducing stimulation fluid. The second experiment served as a control. Trends of aqueous calcium, silica and aluminum in the injection experiment suggest transient dissolution of calcite and feldspar and/or quartz as well as clay and barite precipitation; mineralogic evidence was limited to calcite dissolution and barite precipitation. The rock maintained reducing conditions (Eh = +0.08 to -0.16 V) despite injection of oxidizing stimulation fluid (Eh = 1.1 V). pH of the stimulation fluid-formation water mixture evolved from 2.3 to pre-injection values (~6) within 24 h. The results indicate that mineral dissolution and precipitation manifest in macropores with barely detectable alteration to micro- and mesopores. Formation water or organic matter in pores may have inhibited access of stimulation fluid to micro- and mesopores

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Distribution of rare earth elements and other critical elements in beneficiated Pennsylvania anthracites

The Pennsylvania Anthracite Fields are in a complex tectonic and metamorphic terrain, historically hosting one of the largest concentrations of coal mining in the USA. Anthracite mining now largely consists of the surface mining of the pillars remaining from the prime years of underground mining. The geochemical study of the sized coal products and the refuse (largely rock) from three preparation plants (breakers) demonstrates that Principal components analysis (PCA) of select major oxide, minor element, and rare earth elements illustrates some differentiation among the products from the individual plants. The rock samples, with abundant quartz and metamorphic Al-Si minerals and with a lower ash-basis REE concentration than the coals, were distinctly separated from the coal samples on the PCA plots. Plots of Gd N /Gd N * vs. Eu N /Eu N * and Ce N /Ce N * vs. Eu N / Eu N * showed differentiation between the individual suites of coals showed that the refuse samples had distinct REE distributions compared to the associated coals. Further, several minor and trace elements show enrichments in the coal samples. Lithium, with concentrations of up to 314 ppm on an ash basis, is among the most promising of the critical elements, exceeding the enrichment of the REY and Sc.

58 GEOSCIENCES↗
DOE OSTI2022

How to avoid multiple scattering in strongly scattering SANS and USANS samples

Small Angle Neutron Scattering (SANS) and Ultra Small Angle Neutron Scattering (USANS) are the only available experimental techniques to provide seamless non-destructive measurements of the geometry of the accessible and inaccessible pore structure of rocks from sub-nanopore size to the scale of macropores. They have therefore become the measurement of choice for tight reservoir rocks such as organic rich shales. A simplifying assumption in the analysis is, however, that during the path of neutrons through the sample each neutron is only scattered once. Shales are samples with a high scattering power and Multiple Scattering (MS) may occur which requires special modelling for deconvolution of the results. The approach to avoid MS is to simply reduce the sample thickness to <0.15–0.5 mm. Here, in this work, we present a systematic method on wavelength selection and preparation of samples to optimise extraction of microstructural data and minimise parasitic errors. Experimentally measured SAS transmission (TSAS) values are used as a practical criterion for estimation of the extent of MS. Generous beamtime allocations allowed robust testing revealing that sample thicknesses can be twice as thick as predicted using the standard protocol. Analysing thicker samples is particularly beneficial for statistically relevant characterisation of heterogeneous samples making the new protocol the method of choice for such samples.

(U)SANS↗
DOE OSTI2022

Effect of oxidative breakers on organic matter degradation, contaminant mobility and critical mineral release during shale-fracturing fluid interactions in the Marcellus Shale

Oil and gas production from organic-rich shale formations has become viable through advancements in multistage hydraulic fracturing. However, fluctuations in oil and gas prices, coupled with the sharp decline in gas production after the initial days of fracturing operations, drive operators to devise strategies for enhancing hydrocarbon production. The use of highly reactive fracturing fluids that include strong oxidizing agents, also known as breakers, can potentially increase well productivity. The oxidative breakers are used to reduce the viscosity of gel-based fluids after the proppant is transported into fracture zones of the target formation. These breakers can also degrade the organic matter, enhancing the release of hydrocarbons. However, chemical byproducts generated by the interaction of oxidative breakers with the shale matrix have not been extensively studied. Here, this study investigated the fluid-rock interactions between Marcellus Shale and synthetic hydraulic fracturing fluid (HFF) solutions comprising three different oxidative breakers, i.e., ammonium persulfate, sodium bromate, and sodium hypochlorite, commonly used in the Appalachian Basin, USA. Our results demonstrate that the type of oxidizing breakers used in the HFF controlled the type and amount of volatile organic compounds (VOCs) and semi-volatile organic compounds (SVOCs) produced. In all HFF reacted effluents, we observed the transformation of VOCs and the presence of organic acids in variable amounts. However, effluents from HFF containing sodium bromate and sodium hypochlorite breakers showed the presence of several halogenated organic compounds. Changes in major ions and mineralogy indicate that carbonate dissolution and barite precipitation were ubiquitous in all shale reacted effluents. Our results also demonstrate that the addition of oxidative breakers increased the concentration of several major and trace elements in the effluents. These elements fall under the critical mineral (CM) or critical element category due to their high demand in emerging technologies and susceptibility to supply chain disruption due to variety of factors. However, before oxidative breakers can be used at a larger field scale to enhance the release and recovery of CM and hydrocarbons, a better understanding is required of the potential environmental impacts associated with the generation and transformation of contaminants during breaker-fluid-rock interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Integrated experimental studies of pore structure and fluid uptake in the Bossier Shale in eastern Texas, USA

Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.

58 GEOSCIENCES↗
DOE OSTI2025

Geologic hydrogen as an emerging fuel: experimental insights, thermodynamics, kinetics, and reactive transport modeling

Geologic hydrogen (GeoH 2 ) is emerging as a viable clean energy source. It is largely produced through serpentinization, a geological process in which ultramafic rocks react with water under suitable temperature and pressure. Here, this review synthesizes the current understanding of H 2 generation by serpentinization, with an emphasis on reaction mechanisms, kinetics, and thermodynamics, as well as on modeling flow and transport of reacting fluids in geological formations. We describe the role of mineral assemblages, such as olivine and pyroxene, fluid-rock interactions, and catalytic surfaces, in influencing GeoH 2 yield and reaction rates. By integrating models of reaction kinetics, subsurface reactive flow and transport, and the serpentinization process, and by accounting for the thermodynamic state of the system, this review aims to guide future GeoH 2 research and to evaluate the potential of natural hydrogen as a sustainable clean energy source.

Moradi, Rasoul [Univ. of Southern California, Los ↗
DOE OSTI2020

Sulfur isotope fractionation between aqueous and carbonate-associated sulfate in abiotic calcite and aragonite

Sulfate (SO 4 2– ) incorporated into calcium carbonate minerals enables measurements of sulfur (S) isotope ratios in carbonate rocks. This Carbonate Associated Sulfate (CAS) in marine carbonate minerals is thought to faithfully represent the S isotope composition of the seawater sulfate incorporated into the mineral, with little or no S isotope fractionation in the process. However, comparison between different calcifying species reveals both positive and negative S isotope fractionation between CAS and seawater sulfate, and a large range of S isotope ratios can be found within a single rock sample, depending on the component measured. To better understand the isotopic effects associated with sulfate incorporation into carbonate minerals, we precipitated inorganic calcite and aragonite over a range covering more than two orders of magnitude of sulfate concentration and precipitation rate. Coupled measurements of CAS concentration, S isotope composition and X-ray absorption near-edge spectra (XANES) permit characterization and explanation of the observed dependence of S isotope fractionation between CAS and aqueous sulfate (CAS-SO 4 2– isotope fractionation) on sulfate concentration and precipitation rate. In aragonite, the CAS-SO 4 2– isotope fractionation is 1.0 ± 0.3‰ and independent of the sulfate (and CAS) concentration. In contrast, we note the CAS-SO 4 2– isotope fractionation in calcite covaries strongly with the sulfate concentration and weakly with the precipitation rate, between values of 1.3 ± 0.1 and 3.1 ± 0.6‰. We suggest that the correlation between aqueous sulfate concentration and CAS-SO 4 2– isotope fractionation in calcite reflects a dependence of the equilibrium S isotope fractionation on the concentration of CAS, through the effect of the sulfate impurity on the carbonate mineral’s energetic state.

58 GEOSCIENCES↗
DOE OSTI2021

Iron, magnesium, and titanium isotopic fractionations between garnet, ilmenite, f ayalite, biotite, and tourmaline: Results from NRIXS, ab initio , and study of mineral separates from the Moosilauke metapelite

Interpreting isotopic signatures documented in natural rocks requires knowledge of equilibrium isotopic fractionation factors. Here, we determine equilibrium Fe isotope fractionation factors between several common rock-forming minerals using a comparative approach involving three independent methods: (i) isotopic analyses of natural minerals from a metapelite from Mt. Moosilauke, New Hampshire, for which equilibration temperature and pressure are well constrained to be near the aluminosilicate triple point (T ≃ 500 °C, P ≃ 4 kbar), (ii) Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements of Fe force constants of minerals, and (iii) Density Functional Theory (DFT) ab initio calculations of Fe force constants of minerals. The minerals studied for Fe isotopes include, in increasing order of their β-factors: garnet < ilmenite ≈ fayalite < biotite < tourmaline < muscovite ≈ plagioclase. Some of this ordering is affected by the presence of Fe 3+ in the minerals, which tends to form stiffer bonds and be associated with heavy Fe isotope enrichments relative to Fe 2+ . We are, however, able to assess the magnitude of the effect of the ratio Fe 3+ /ΣFe on equilibrium fractionation factors, notably on the ilmenite-hematite solid solution. Equilibrium Fe isotopic fractionation factors between garnet, ilmenite, biotite, tourmaline and fayalite are determined. We also report Mg and Ti isotopic compositions of selected Moosilauke minerals that allow us to better constrain the equilibrium fractionation factors for garnet-biotite-tourmaline (Mg isotopes) and biotite-ilmenite (Ti isotopes). We show how the newly determined equilibrium fractionation factors can be used to address diverse problems in Earth and planetary sciences, notably (i) Fe and Mg isotopic fractionation during anatexis, (ii) Fe isotopic fractionation in lunar ilmenite, and (iii) Ti isotopic fractionation during fluvial transport of minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗