SEARCH · Engineering Papers

Results for “lithology”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Long‐term growth and persistence of granitic inselbergs in a semi‐arid cratonic landscape

Granitic inselbergs are enduring landforms that punctuate planation surfaces in tectonically stable settings, yet the mechanisms and timescales of their persistence remain debated. Here, in this study, we present a multinuclide cosmogenic dataset ( 10 Be, 26 Al and 21 Ne) from bedrock, fluvial sediments and planation surfaces in the semi-arid Central Ceará Domain, north-eastern Brazil, to quantify denudation rates and reconstruct the long-term dynamics of residual relief formation. This triple-nuclide approach enables independent validation of denudation contrasts and long-term exposure histories, providing quantitative constraints on inselberg evolution that a single-nuclide approach could not resolve. Bedrock denudation rates (1–4 m Myr −1 ) are systematically lower than basin-wide rates (8–14 m Myr −1 ), revealing a regional pattern of vertical differential denudation. When combined with summit elevations and surface exposure ages of up to 5 Myr, these contrasts suggest that cumulative relief growth has been operating over timescales of tens of millions of years (1–57 Myr). Morphological and structural observations further indicate localized mass wasting and heterogeneous erosion styles, modulated by rock fabric and joint density. Cosmogenic nuclide inventories in detrital sediments reveal partial reworking from older sediment stocks, including surface pebbles with exposure ages up to 5 Myr, which record multistage of burial and re-exposure histories and reflect episodic remobilization in low-connectivity, semi-arid catchments. Rather than static relicts of ancient surfaces, inselbergs in this cratonic setting emerge as dynamic landforms actively sustained by slow but persistent differential denudation. These results provide new constraints on long-term cratonic landscape evolution and underscore the interplay between lithological resistance, structural inheritance and sediment transfer processes in sustaining topographic relief. They also call for a critical reassessment of classical inselberg models, suggesting that hybrid mechanisms, rather than single-process paradigms, more accurately capture the complexity of residual landform development in cratonic interiors.

Cosmogenic nuclides↗

DOE OSTI2023

A phased workflow to define permit‐ready locations for large volume CO 2 injection and storage

Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

58 GEOSCIENCES↗

DOE OSTI2020

Relating land surface, water table, and weathering fronts with a conceptual valve model for headwater catchments

Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the earth-air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re-orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes toward an “ideal hill”. In this steady state, reaction-front valves partition water volumes into shallow and deep flowpaths. Furthermore, these flows dissolve low- and high-solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Effects of structure and volcanic stratigraphy on groundwater and surface water flow: Hat Creek basin, California, USA

Abstract Hydrogeologic systems in the southern Cascade Range in California (USA) develop in volcanic rocks where morphology, stratigraphy, extensional structures, and attendant basin geometry play a central role in groundwater flow paths, groundwater/surface-water interactions, and spring discharge locations. High-volume springs (greater than 3 m 3 /s) flow from basin-filling (<800 ka) volcanic rocks in the Hat Creek and Fall River tributaries and contribute approximately half of the average annual flow of the Pit River, the largest tributary to Shasta Lake. A hydrogeologic conceptual framework is constructed for the Hat Creek basin combining new geologic mapping, water-well lithologic logs, a database of active faults, LiDAR mapping of faults and volcanic landforms, streamflow measurements and airborne thermal infrared remote sensing of stream temperature. These data are used to integrate the geologic structure and the volcanic and volcaniclastic stratigraphy to create a three-dimensional interpretation of the hydrogeology in the basin. Two large streamflow gains from focused groundwater discharge near Big Spring and north of Sugarloaf Peak result from geologic barriers that restrict lateral groundwater flow and force water into Hat Creek. The inferred groundwater-flow barriers divide the aquifer system into at least three leaky compartments. The two downstream compartments lose streamflow in the upstream reaches (immediately downstream of the groundwater-flow barriers) and gain in downstream reaches with the greatest inflows immediately upstream of the barriers.

Geology↗

DOE OSTI2020

An Experimental Method to Measure Gaseous Diffusivity in Tight and Partially Saturated Porous Media via Continuously Monitored Mass Spectrometry

Detection of radioxenon and radioargon produced by underground nuclear explosions is one of the primary methods by which the Comprehensive Nuclear-Test–Ban Treaty (CTBT) monitors for nuclear activities. However, transport of these noble gases to the surface via barometric pumping is a complex process relying on advective and diffusive processes in a fractured porous medium to bring detectable levels to the surface. To better understand this process, experimental measurements of noble gas and chemical surrogate diffusivity in relevant lithologies are necessary. However, measurement of noble gas diffusivity in tight or partially saturated porous media is challenging due to the transparent nature of noble gases, the lengthy diffusion times, and difficulty maintaining consistent water saturation. Here, the quasi-steady-state Ney–Armistead method is modified to accommodate continuous gas sampling via effusive flow to a mass spectrometer. An analytical solution accounting for the cumulative sampling losses and induced advective flow is then derived. Experimental results appear in good agreement with the proposed theory, suggesting the presence of retained groundwater reduces the effective diffusivity of the gas tracers by 10–1000 times. Finally, by using a mass spectrometer, the method described herein is applicable to a broad range of gas species and porous media.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Geophysical methods reveal the soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie

Abstract Purpose The land-lake interface is a unique zone where terrestrial and aquatic ecosystems meet, forming part of the Earth’s most geochemically and biologically active zones. The unique characteristics of this interface are yet to be properly understood due to the inherently high spatiotemporal variability of subsurface properties, which are difficult to capture with the traditional soil sampling methods. Geophysical methods offer non-invasive techniques to capture variabilities in soil properties at a high resolution across various spatiotemporal scales. Methods We combined electromagnetic induction (EMI), electrical resistivity tomography (ERT), and ground penetrating radar (GPR) with data from soil cores and in situ sensors to investigate hydrostratigraphic heterogeneities across land-lake interfaces along the western basin of Lake Erie. Results EMI revealed high spatial heterogeneities in ECa distribution across the land-lake interfaces, with higher values in the wetland and transition zones compared to the upland zone. Soil ECa maps matched soil maps from a public database with the hydric soil units delineated as high conductivity zones (ECa > 40 mS/m). ERT and GPR showed vertical variation in soil properties with clear stratigraphic boundaries, and correlation of ERT profiles with lithologs from piezometers revealed the stratigraphic units of silt–clay and till sequence down to 3.5 m depth which are consistent with the surficial geology of the study area. Conclusions These results validate the use of multiple geophysical methods for extrapolating soil properties and mapping stratigraphic structures at land-lake interfaces, thereby providing the missing information required to improve the earth system model (ESM) of coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Formation and interior evolution of Mars

The formation and earliest evolution of Mars involved core formation followed by magma ocean solidification and primordial crust formation. Here, in this study, we review the formation and early interior evolution of Mars as understood through meteorite studies. Martian meteorites exhibit ages that span the entire history of the Solar System, but most derive from compositionally distinct sources that formed early in the history of Mars. Geochemical and isotopic studies of Martian meteorites have provided critical information on Mars’ building blocks, bulk composition, and its differentiation into a core, mantle and crust. Radiogenic isotope systems provide vital insights into the petrogenesis of Martian meteorite source regions and into the nature and timescales of primordial differentiation on Mars. Application of short-lived isotope chronometers to Martian meteorites indicates that the Martian core formed early within <10 million years (Ma), while magma ocean differentiation and crust formation occurred soon thereafter, within the first tens of Ma of Solar System history. Bulk silicate Mars is variably depleted in moderately volatile elements relative to chondrites, but to a lesser extent than the Earth. The nature and evolution of more volatile elements and volatile species on Mars (e.g., H 2 O, C, S, halogens) remains an active area of investigation, but overall suggest that Mars is more enriched in volatile elements than the Moon but less than the Earth. While geochemical exploration of Mars coupled with new discoveries of Martian meteorites has provided a wealth of information, the compositional and lithologic diversity on Mars significantly exceeds that seen in the suites of Martian meteorites currently studied. This implies that our current perspective on the planet's evolution is very likely biased, highlighting the need for future sample return missions to Mars.

Kruijer, Thomas S.↗

DOE OSTI2021

Rapid determination of supercritical CO 2 and brine relative permeability using an unsteady-state flow method

In this study, relative permeability of supercritical CO 2 (scCO 2 ) and brine was determined in reactive and non-reactive rock cores using a combination of unsteady-state methodology and computed tomography. Experiments were conducted using a medical grade CT scanner to determine saturation using a custom Python script. The saturation and differential pressure across the core were then used to derive four empirical constants to calculate relative permeability. This methodology increases temporal efficiency while reducing experimental complexity. Additionally, we show that the method can be used to determine scCO 2 relative permeability in a wide range of lithologies and flow rates, and with the ability to account for matrix dissolution during scCO 2 flooding.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Well–based monitoring of CO 2 geological sequestration operations in saline aquifers: Critical insights into key questions

Geological carbon sequestration in saline aquifers is one of the most promising strategies to help mitigate emissions of CO 2 to the atmosphere. Significant challenges in ensuring the security of the sequestration process rest in the evolution and expansion of the CO 2 plume in the subsurface. The ability to track the movement of the injected CO 2 poses another challenge. Critical questions related to the integrity of the sequestration operations in saline aquifers relate to plume characteristics that can we monitor using well-based variables. We addressed this and related questions using an integrated modeling framework through a numerical investigation of carbon sequestration in saline aquifers during long-term and post-injection periods. This modeling paradigm incorporates the effect of structural, geological, and petrophysical characteristics. That way, we can account for critical physicochemical processes, rock-fluid interactions, and lithology dependencies. The well fluid variables investigated include fluid composition, pH, fluid density, and ion activity. We learned that fluid property analytics and diagnostics can be powerful tools to estimate the movement of CO 2 and its storage in different trapping mechanisms. These analytics can help optimize operational aspects and simplify reservoir-scale models while still reflecting the complex nature of the CO 2 interactions underground and offering insights into plume evolution.

58 GEOSCIENCES↗

DOE OSTI2025

Geologic hydrogen: From natural occurrences to anthropogenic generation – A review of fundamentals, potential, challenges and prospects

Growing demand for hydrogen is exposing the environmental and economic limits of reforming-based and carbon-managed supply chains, while the scale-up of electrolytic capacity remains capital-constrained. Geologic hydrogen, defined as molecular H₂ generated and stored within the Earth's crust offers a complementary, potentially lower-cost resource, yet exploration is still ad hoc. This review (1) revisits a global inventory of confirmed hydrogen seeps and subsurface occurrences; (2) analyzes the controlling reactions, migration pathways, and trapping conditions governing these occurrences; (3) proposes a process-based geologic hydrogen system concept analogous to, yet distinct from, the petroleum system; and (4) evaluates potential geologic hydrogen systems within the United States as a representative case study. Here, we contrast natural systems powered by serpentinization, mantle degassing or radiolysis with anthropogenic systems that stimulate the same reactions or convert in-situ hydrocarbons. Stable hydrogen accumulations require generation rates that outpace combined physical, chemical and microbial losses; the Bourakébougou field (Mali) exemplifies a self-recharging, free-gas reservoir sustained by meteoric-water serpentinization beneath an efficient caprock. Prospective geologic hydrogen resources are likely to occur in regions where iron-rich lithologies, deep-seated faults, and low-permeability sealing formations coexist. Applying this principle, we highlight three promising hydrogen play types in U.S. geological terrains: ophiolite belts (Appalachian and Californian regions), the Midcontinent Rift and the Lake Superior banded‑iron formations. Multiphysics numerical models and positive-unlabeled machine-learning workflows help to accelerate play screening and de-risk future production; yet, reaction kinetics, stimulation strategies, and full techno-economic and life-cycle assessments remain pivotal knowledge gaps.

Anthropogenic hydrogen generation↗

DOE OSTI2021

Microfracture-pore structure characterization and water-rock interaction in three lithofacies of the Lower Eagle Ford Formation

The microfracture-pore structures and water-rock interactions of eight samples of three typical lithofacies (three limestone, three wackestone, and two mudstone) from the Lower Eagle Ford Formation have been investigated using an integrated methodology. The methods used were X-ray diffraction (XRD), total organic carbon (TOC) content, pyrolysis, thin-section petrography, scanning electron microscopy (SEM), mercury intrusion porosimetry (MIP), spontaneous water imbibition, and water vapor & nitrogen physisorption. Petrophysical properties such as organic richness, maturity, mineral composition, types and patterns of microfractures and pores, porosity, pore-throat size distribution, pore size distribution, pore connectivity, water adsorption and desorption behavior were determined from these methods, and the results compared among the three lithofacies. For example, the three limestone EAS samples and three wackestone DM samples show a low porosity range from 2.02–3.35%, whereas the two mudstone VRU samples exhibit a much large average porosity of 11.7% by MIP results. Nitrogen physisorption results show that pore size distributions have a good similarity in N2 adsorption volume for the samples with similar lithology. Microfractures are morphologically divided into seven types, with their definitions & pictorial examples presented and possible formation mechanisms proposed. In addition, the water vapor physisorption behavior and their controlling factors are discussed. Here, the water adsorption capacity shows a positive relationship (R 2 = 0.90) with clay content and a negative relationship (R 2 = 0.99) with carbonate content. In general, this integrated study of microfracture-pore structure characterization and water-rock interaction provides some insights of water vapor adsorption behavior in different lithofacies and whether microfractures could be an important pathway for fluid flow.

58 GEOSCIENCES↗

DOE OSTI2020

Onset of magma ocean solidification on Mars inferred from Mn-Cr chronometry

The mantle of Mars probably differentiated through the crystallization of a magma ocean during the first tens of million years (Ma) of Solar System evolution. However, the exact timescale of large-scale silicate differentiation of the martian mantle is debated, and in particular, it remains unclear when differentiation commenced. In this paper, we applied the short-lived 53 Mn- 53 Cr system to martian meteorites in order to date the onset of large-scale mantle differentiation on Mars. The new Cr isotope data demonstrate that martian meteorites exhibit no resolvable radiogenic 53 Cr variations, and instead have a uniform +20.2±1.2 (95% conf.) parts-per-million excess in 53 Cr/ 52 Cr relative to the terrestrial mantle. The investigated groups of martian meteorites are lithologically varied and derive from diverse mantle sources that probably had variable Mn/Cr. Hence, the lack of 53 Cr variability among martian meteorites demonstrates that silicate differentiation on Mars occurred after the extinction of 53 Mn. Provided that the sources of the martian meteorites have Mn/Cr variations that are typical of the terrestrial planets, this result implies that the onset of large-scale silicate differentiation must have occurred later than 20±5 Ma after Solar System formation. The onset of silicate differentiation on Mars inferred here is significantly later than time estimates for segregation of the martian core which conservatively occurred within <10 Ma after Solar System formation. Thus, the new Mn-Cr data imply that there was a small, but resolvable, time gap of at least 5 Ma between core formation and magma ocean solidification on Mars. If the age of core segregation is taken at face value, our results imply that the martian magma ocean remained mostly molten over several Ma. This inferred longevity of the magma ocean is inconsistent with thermal models predicting rapid (<1 Ma) solidification of the martian magma ocean. Although there is currently no unique solution to this conundrum, our results can potentially be explained by a protracted history of impact bombardment that delayed differentiation in a shallow magma ocean on Mars, or perhaps more readily, by the presence of an early and dense atmosphere that acted as an insulator and prevented the magma ocean from cooling quickly.

58 GEOSCIENCES↗

DOE OSTI2021

Organic sulfur fluxes and geomorphic control of sulfur isotope ratios in rivers

Pyrite oxidation plays a critical role in the relationship between weathering and climate, and its impact on the global carbon cycle has previously been constrained through inversion models utilizing observations of river sulfate ($SO^{2–}_{4}$) and its 34 S/ 32 S isotope ratio (δ 34 S SO4 ). However, measurements from some rivers have suggested that SSO4 can be substantially impacted by processes such as microbial sulfate reduction and/or sulfur assimilation and cycling, rather than simply reflecting a weighted mixture of lithologic sulfur sources. To study the prevalence and controls on $SO^{2–}_{4}$ transformations, in this study we measured dissolved major element concentrations and δ 34 S SO4 in river water samples from throughout western Iceland. Our analyses focused on samples from a small catchment hosting the Efri Haukadalsá river, a system with relatively uniform and isotopically constrained basaltic bedrock. We also measured sediment δ 34 S and sulfur speciation using sulfur K-edge X-ray absorption spectroscopy on sediment and vegetation samples from this catchment. Values of dissolved δ 34 S SO4 in the Efri Haukadalsá ranged from 2.5‰ to 23.7‰ and had a linear relationship with Cl – /$SO^{2–}_{4}$ ratios, indicating that $SO^{2–}_{4}$ predominantly derived from basalt weathering and meteoric precipitation. The lower δ 34 S SO4 values were found in fluvial valleys with V-shaped cross sections, while higher values of δ 34 S SO4 occurred in U-shaped, glacially eroded valleys with thick alluvial fills blanketing the valley floor. Spectroscopic observations identified organic sulfur phases in suspended river sediment, floodplain deposits, and vegetation. Mass balance calculations quantified the organic sulfur flux as less than 10% of $SO^{2–}_{4}$ export, and sediment δ 34 S values were comparable to river δ 34 S SO4 . We interpreted these isotopic and chemical patterns as reflecting differences in the availability of unweathered bedrock across the Efri Haukadalsá catchment, with V-shaped valleys having greater access to fresh sulfide-bearing minerals than alluviated U-shaped valleys; this interpretation is in contrast to one in which the elevated δ 34 S SO4 values reflect fractionation during sulfur transformations along alluvial reaches. These results validated the application of river inversion models for constraining weathering fluxes and affirmed that pyrite oxidation globally, even in the presence of river sulfur cycling, modulates the abundance of atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Gallium isotopic constraints for the origin of the Earth-Moon system

Comparatively heavy isotopic compositions of moderately volatile elements (MVE) in lunar rocks have been advocated to reflect the loss of light isotopes during devolatilization processes from the Moon. In this study we present new gallium (Ga) isotope data for lunar highland rocks, with a focus on the Ferroan Anorthosite Suite (FAS). These are commonly thought to be direct crystallization products from the late lunar magma ocean (LMO) and should contain the majority of the Ga inventory of the Moon. As such, FAS rocks are crucial for identifying the processes that drove Ga isotope fractionation as well as for inferring the Ga isotopic composition of the bulk Moon. Our data reveal that FAS samples have a range in δ 71 Ga from -0.27 to 0.22‰ and are generally isotopically light in Ga compared to other lunar lithologies, but straddle values typical of terrestrial rocks. Although Ga is defined as an MVE, these Ga isotope variations do not correspond with concentrations of more volatile elements, indicating that Ga isotope variations in the FAS are not primarily controlled by devolatilization processes. Instead, the Ga isotopic compositions of bulk FAS rocks broadly correlate with the composition of plagioclase, with the calcium content of plagioclase decreasing as Ga becomes isotopically heavier. This suggests that fractionation of Ga isotopes in FAS rocks was caused by the preferential incorporation of isotopically light Ga into plagioclase during the later solidification stages of the LMO. The progressive crystallization and extraction of plagioclase forces the residual melt towards increasingly heavier Ga isotope ratios, corroborating similar conclusions derived from correlations between δ 71 and Eu* in the mare basalt suite. Using Ga isotope partitioning calculations, we demonstrate that an isotope fractionation coefficient between plagioclase and coexisting melt of -0.3 to -0.4‰ could explain the observed range of δ 71 values in FAS, mare basalt suite rocks, and KREEP. These calculations allow for a first order estimate of the Ga isotopic composition of the bulk silicate Moon prior to plagioclase fractionation and suggest it was close to the composition of the bulk silicate Earth. This would imply that the Moon did not lose a substantial fraction of its Ga inventory during accretion, consistent with new constraints from Rb-Sr isotope systematics that indicate the Moon's volatile deficit was primarily inherited from Theia. In conjunction with the overlap in non-mass-dependent isotope ratios, these collective observations could be reconciled if Theia and the proto-Earth formed in similar regions of the inner Solar System that were already volatile-depleted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A microstructural signature of the coesite-quartz transformation: New insights from high-pressure experiments and EBSD

Ultra-high pressure (UHP) metamorphism is difficult to identify in continental crust as few petrological barometers are suitable for dominantly felsic lithologies. In such cases, burial to extreme depths is commonly identified through the preservation of coesite, a high-pressure polymorph of SiO 2 that typically forms at depths exceeding ∼ 100 km (i.e., > 2 GPa pressure). Unfortunately, coesite readily transforms to quartz upon exhumation, meaning that UHP terranes may often be overlooked. While some studies have suggested that quartz may inherit an orientation signature indicative of former coesite, both the specific nature of this signature and the conditions favouring its development remain uncertain. Here, to address this problem, we combine electron backscatter diffraction analysis of natural and experimental samples to explore microstructural evolution across the coesite-quartz phase transformation. We demonstrate that neighbouring domains of quartz commonly feature an 84 ± 4° rotation of [c] axes around the pole of a common {m} plane. This orientation relationship is a product of epitaxy, whereby the {$11\bar{2}2$} Japan twin plane in quartz nucleates on the (010) plane in coesite. In supercell simulations, the nucleation of Japan twins can be explained by the energetically favourable alignment of quartz tetrahedra on parental coesite tetrahedra. Through experiments, we demonstrate that this signature emerges over a broad range of conditions, regardless of the availability of nucleation sites (e.g., grain boundaries) or the density of crystal lattice defects (e.g., dislocations). Overall, our work provides a quantitative and unambiguous tool for identifying UHP terranes from quartz in isolation.

Coesite↗

DOE OSTI2023

Wettability variation and its impact on CO 2 storage capacity at the Wyoming CarbonSAFE storage hub: An experimental approach

Meeting global and national net zero carbon emission targets will require geologic carbon disposal. The U.S. Department of Energy (DOE) has accordingly funded significant research in this area, including the Wyoming CarbonSAFE project at Dry Fork Station (DFS) in Campbell County, Wyoming. This work studied wettability on micro- and macro-scales, CO 2 storage potential, and the correlation between the two to support the Wyoming CarbonSAFE project’s subsurface assessment. During the study, a target formation’s wettability was found to affect how much CO 2 can be stored in a given formation. Here, in this study, representative rock samples were selected from the target storage formations— Lakota, Hulett, and Minnelusa—based on the heterogeneity of the lithology, permeability, and porosity of the respective formations. The rock samples are all fine-grained sandstone with variable cementation and bedding structure, including different scales of laminated bedding. The porosity and permeability vary within the range of 9.0–14.3% and 0.1–28.9 mD, respectively. These rock samples were prepared for the micro-scale wettability (contact angle measurement), macro-scale wettability (wettability index derived from unsteady-state flow characterization for the core plugs), and CO 2 storage evaluation. The macro-scale experiments suggested that wettability appeared to dominate the CO 2 storage potential performance during the drainage process, where less water-wet behavior promoted higher CO 2 storage potential. The micro-scale wettability tests showed that the rock samples at the studied reservoir conditions behaved water-wet and became more water-wet as pressure increased. This kind of wettability change discourages further CO 2 storage potential yet benefits the CO 2 residual trapping as the CO 2 injection proceeds for the studied area. The results allow the recommendation of the best reservoir candidate for storage based on wettability that affects CO 2 storage. The work presented in this study provides valuable insights into wettability’s effect on the CO 2 storage capacity and wettability’s importance when identifying the optimal CO 2 storage formation to meet the project’s goals.

58 GEOSCIENCES↗

DOE OSTI2026

Integrated geological, economic, and risk assessment of underground hydrogen storage

Underground hydrogen storage (UHS) is a promising option to buffer variable renewable power and support the hydrogen economy. Yet this technology is early-stage, and key uncertainties remain about how coupled processes such as geochemical and geomechanical impacts affect long-term storage security. In this paper, we seek to integrate the results from a large, multi-scale research program to quantify core feasibility metrics and situate these studies within ongoing efforts. Through this work, we evaluate hydrogen recoverability during storage by assessing diffusive losses, losses to geochemical interactions and their subsequent impact on geomechanical properties, and losses to residual saturation during injection and withdrawal. Molecular and laboratory studies show that lithologic heterogeneity, pore geometry, and mineral surface chemistry govern hydrogen–rock interactions, controlling loss pathways. High-pressure coreflooding is used to estimate relative permeabilities and measure parameters needed for reservoir-scale sensitivity analyses. These core-scale experiments found early hydrogen breakthrough at low saturations driven by capillary and viscous fingering, which limits pore-space utilization at early times, while repeated injection and withdrawal cycles improve deliverability over time. Technoeconomic factors and purity requirements for end uses are also considered as part of our comprehensive feasibility analysis. Finally, we underscore the importance of geophysical monitoring and tailored injection strategies to maintain integrity and efficiency. Together, these results provide a quantitative foundation for safe, scalable deployment of UHS and highlight priorities for future integrated studies.

08 HYDROGEN↗

DOE OSTI2025

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗