Search NASA⌕ Search

SEARCH · Search NASA

Results for “sedimentary 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2023

The role of stress and fluid saturation on the acoustic response of fractured rock

Standard rock physics models are formulated to describe the behavior of porous sedimentary reservoirs, with clean sandstones being the archetypal system; however, many situations demand geophysical monitoring of rocks with significantly different structures, such as low porosity, fractured reservoirs. Conventional models also suggest that these “stiff” reservoirs can be challenging to monitor seismically due to small fluid substitution effects, but the presence of fractures leads to stress dependence which may be leveraged for remote monitoring purposes. Using samples from the Duperow Formation (dolostone) obtained from the Danielson test well in Kevin Dome, MT, we conducted ultrasonic and multi-scale structural (profilometry, synchrotron micro-tomography, pressure sensitive film) measurements on naturally fractured core in order to characterize the effects of fluid substitution and effective stress on the acoustic response of fractured reservoir rock with a focus in particular on the textural and seismic characteristics of natural fractures. We find that changes in effective stress can yield changes in velocity of up to 20% and changes in attenuation up to 200%. Measured fluid substitution effects are resolvable, but stress effects dominate. These measurements provide insight into the physical processes controlling acoustic response of fractured rocks in general and can also be used to inform monitoring efforts in fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2023

A Geo-Data Science Method for Assessing Unconventional Rare-Earth Element Resources in Sedimentary Systems

Abstract Rare-earth elements (REEs) supply raw materials that constitute many of our modern critical infrastructure, defense, technology, and electrification needs. Despite REE accumulations occurring in conventional bedrock and ion-adsorption deposits sourced from weathering of igneous rocks, unconventional host materials such as coal and related sedimentary strata have been identified as promising sources of REEs to meet growing demand. To maximize the potential of unconventional resources such as REE-coal systems, new approaches are needed overcome challenges from mineral systems with no known deposits and areas with sparse geochemical data. This article presents a systematic knowledge-data resource assessment method for predicting and identifying REE resource potential and occurrence in these unconventional systems. The method utilizes a geologic and geospatial knowledge-data approach informed and guided by REE accumulation mechanisms to systematically assess and identify areas of higher enrichment. An assessment of the Powder River Basin is presented as a test case to demonstrate the method workflow and results. The key output is a potential enrichment score map reported with varying confidence levels based on the amount of supporting evidence. Results from the test case indicate several locations with promising potential for different types of coal-REE deposits, demonstrating the viability of the method for exploration and assessment of unconventional REE resources. The method is flexible by design and, with sufficient applicable knowledge and data, can be adapted for assessing critical mineral systems in other sedimentary systems as well.

58 GEOSCIENCES↗
DOE OSTI2020

Characterizing the uppermost 100 m structure of the San Jacinto fault zone southeast of Anza, California, through joint analysis of geological, topographic, seismic and resistivity data

We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (<100 m) fault zone at this structurally complex site. Geological mapping at the surface shows the CF has three main subparallel strands within a <100 m zone with varying degrees of rock damage. These strands intersect units of banded gneiss and tonalite, and various sedimentary units. Near the surface, the weathered but more intact tonalite and gneiss to the southwest have relatively high VP. The low-lying flat sedimentary basins around the two southwestern-most CF strands and elevated damaged gneiss to the northeast have lowest VP <500 m s–1. The high relief of the northeast gneiss unit may in part be explained by its extensive damage and inferred increased relative rock uplift. Resistivity imaging shows the unconsolidated dry basin sediments (maximum >1300 Ohm.m) contrasted against the compacted fine-grained (potentially wet) materials within the CF core and the Bautista Formation (minimum <40 Ohm.m), which is slightly elevated above the flat basins. The inverse relationship between VP (increases) and resistivity (decreases) in the uppermost ~15 m can be characterized as log–log linear with slopes of –2.6 to –4. At depths >30 m, the velocity heterogeneity near the surface merges into larger-scale structures that are generally slower on the northeast side of the CF core compared to the southwest side (as much as ~40 per cent reduction in average VP). A previous study revealed a 20–37 per cent variability in peak ground velocities across the SGB site from local earthquakes. The upper end of that range is associated with the near-surface unconsolidated sedimentary basins and northeast damaged gneiss unit. Preliminary analysis of time-dependent topography mostly shows effects of changing vegetation and anthropogenic activity.

Geochemistry & Geophysics↗
DOE OSTI2024

Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide

The project Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. Sandia Technologies, LLC, and co-investigator Conrad Geoscience Corporation, examined the potential for large-scale, permanent CO2 storage in sedimentary strata within the Newark Rift Basin. The Newark Rift Basin underlies an industrialized, developed region comprising parts of New York, New Jersey, and Pennsylvania. The project characterized and investigated the suitability of Triassic age sedimentary formations for potential geologic CO2 storage. The project team drilled and cored two test wells to define the sedimentary geologic formations underlying the basin and to document or reach basement rock. With this geologic characterization phase, an integration of seismic, geologic, borehole, and formation core results provided a higher resolution assessment of CO2 storage potential. The Stockton Formation is known to be a potentially favorable geologic storage formation in the basin. In 2011, the 1-NYSTA Tandem Lot stratigraphic test well was drilled to a depth of 6,855 feet in the northern portion of the Newark Basin in southern New York State. Approximately 9 miles south-southeast on the Lamont Doherty Campus, TW-4 was drilled and cored in 2013 to a depth of 1,802 feet and contacted apparent igneous basement at a depth of 1,712 feet. Both wells penetrated the Palisades Sill ranging from 800 feet thick in the eastern well to approximately 1,800 feet in thickness at the 1-NYSTA Tandem Lot deep drill site. A diabase sill can provide an excellent seal and dense confining layer for potential CO2 storage reservoirs and flow layers that are situated beneath it within the Stockton Sandstone. The Stockton Sandstone was encountered beneath the sill in the TW-4 well on the Lamont campus, and data integration suggests that it was likely observed near total depth in the deep 1-NYSTA Tandem Lot well. The test wells confirm and define reservoirs are present beneath the sill and offer CO2 storage potential. The integration of geologic and reservoir characterization of well logs, formation cores, and formation fluids indicated Triassic age lacustrine playa lake and mudbank shales of the Upper Passaic Group can provide an effective seal for the porous and permeable underlying sandstone reservoir layers. This project acquired seismic data, drilled borehole well logs, acquired core samples, and integrated these findings to provide a better understanding of the subsurface geologic formations in the Newark Rift Basin. These findings have contributed to a higher degree of accuracy in predicting potential geologic storage opportunities, while refining geologic storage capacity estimates for the indicated reservoirs and flow units.

.las↗
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2020

Microanalytical Approaches to Characterizing REE in Appalachian Basin Underclays

The search for a reliable U.S. domestic source of rare earth elements (REE) is necessary to support the demand of advanced energy applications (e.g., catalysts, electronics, magnets). Sedimentary deposits may be sources for selectively recovering REE and critical metals—specifically the interbedded seat rock, or underclay, that underlies or forms the floor of a coal seam. This material is often a major component of coal waste fines and refuse and thus readily available. This study examines several Appalachian Basin underclays associated with actively mined coal seams as potential feedstocks for the REE. Multimodal microanalytical electron microscopy (SEM, FIB-SEM, EMPA) synchrotron-based µXRF, and image processing techniques are coupled with detailed elemental and mineral data to classify the 2D and 3D petrophysical properties of the materials. The REE contents of Appalachian Basin underclays were measured from 235–399 ppm and predominantly observed as discrete REE-bearing minerals such as monazite and xenotime on the order of 10–100 µm in size. These REE-bearing minerals typically accounted for less than 1% of the scanned areas and volumes under SEM and FIB-SEM analysis, with the exception of regions enriched in crandallite. Synchrotron-based µXRF elemental maps further identified several REE deposition environments in different underclays, including micro-scale (10–100 µm) light REEs co-localizing with Ca and P, micro-scale heavy REEs with Fe, and large-scale light REEs (>200 µm) co-localizing with Sr, Ba, Ca and P.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spatial Seal Database for Prospective Storage Resources in the USA

The goal of the Spatial Seal Database for Prospective Storage Resources in the USA is to provide relevant information and spatial extents of caprock and seal rocks. A lack of aggregated information is readily available that focuses on the caprock and seal units within sedimentary basins. The EPA class VI permit requires an assessment of the confining zone as part of submitting a permit. The data catalog of seal unit names and relevant properties with the seal spatial extent database aims to help provided important data for carbon storage based assessments. The data catalog and database are designed to show what seal data is available in a sedimentary basin and guide stakeholders to the original data source for those datasets.

Pantaleone, Scott↗
DOE OSTI2024

Prospective Seal Unit Spatial Extent Database for U.S. Sedimentary Basins

The Prospective Seal Unit Spatial Extent Database for U.S. Sedimentary Basins contains a series of spatial datasets representing spatial extents of publicly available data for caprock and seal rock units within the Appalachian Basin, Denver-Julesburg Basin, Great Valley Basin (Sacramento and San Joaquin Basins), Illinois Basin, Michigan Basin, San Juan Basin, U.S. Gulf Coast Basin, and Williston Basin. The database is designed to support carbon storage feasibility and resources assessment for carbon transport and storage (CTS) projects while displaying the spatial extent of prospective seal units and provide a guide to the original data source. This database leverages publicly available data resources from authoritative sources (e.g. U.S. Geological Survey, State Geologic Surveys, and published reports), and aims to help guide users to understand the seal unit's spatial coverage and data gaps from the regional to sub-basin/field scale. The database is organized by seal unit/formation, including the spatial extent for data found to be available for the seal unit. The various datasets represented include spatial extents of the lithologic formation, depth to top structural contour maps, and thickness/isopach maps. Included in this submission are the following resources: 1. Geodatabase/Dataset: “prospective-seal-unit-extents-2025.gdb” 2. ReadMe: “readme-prospective-seal-unit-spatial-extent-dataset-2025.pdf” 3. Data Catalog: “prospective-seal-unit-spatial-extents-data-catalog-2025.xlsx” 4. Data Sources Key: “data-source.csv” Please see NETL disclaimers here: https://netl.doe.gov/home/disclaimer

Basin↗
DOE OSTI2020

Assessing the thermal maturity of black shales using vitrinite reflectance: Insights from Devonian black shales in the eastern United States

Thermal maturity of source rocks indicated by vitrinite reflectance (R o ) is an important and reliable parameter to determine the petroleum potential of sedimentary basins. Originally, R o was used to characterize the degree of coalification of coals and thus works well for coal seams, but it has limitations when used to determine the thermal maturity of black shale successions. A suite of six coal-shale pairs from the Upper Devonian black shale formations in the eastern United States (New Albany Shale, Ohio Shale, Dunkirk Shale, and Rhinestreet Shale) was selected to study the applicability of R o in assessing the thermal maturity of black shale successions. The results show that vitrinite in the studied coal samples is dominated by collotelinite, whereas vitrinite in black shales occurs as small dispersed particles (~5 μm) in the mineral matrix. When comparing the size and morphology of dispersed vitrinite particles in shales and collotelinite fragments in coals, vitrinite in shales and adjacent coals should have the same origin. The measured mean random R o of vitrinite in coals ranges from 0.51–0.68%, and is 0.04–0.11% (average 0.07%) lower than that of dispersed vitrinite particles in enclosing shales. This observation contrasts with previously reported R o suppression in liptinite-rich black shales. A high standard deviation of R o measurements indicates a highly heterogeneous nature of vitrinite in black shales. Three mechanisms can be envisioned to contribute to the higher R o of dispersed vitrinite particles in black shales relative to enclosed coals. First, small vitrinite particles may become more oxidized because small particles are more frequently suspended during transport to the site of deposition and reworking of bottom sediments after deposition than large pieces of driftwood that will turn into enclosed coal lenses during burial diagenesis. Second, misidentification of zooclast (e.g., chitinozoan) fragments as vitrinite due to loss of diagnostic morphology. Third, various macerals in the vitrinite group (e.g., corpogelinite and collotelinite) may have had different original reflectance values. Based on the petrographic characteristics of vitrinite and vitrinite-like particles, caution should be applied when using the R o values of dispersed vitrinite and vitrinite-like particles in black shales as indicators of thermal maturity.

58 GEOSCIENCES↗
DOE OSTI2023

Geocatalytically generated methane from low-maturity coal and shale source rocks at low temperatures (80–120 °C) over 52 months

Geocatalytic methanogenesis has been proposed to contribute to methane generation from low-maturity coal and shale source rocks. This study contributes further evidence for geocatalytic methanogenesis from low-maturity source rocks based on long-term experiments lasting up to 52 months. Samples from the Upper Devonian New Albany Shale (R o 0.54 %) and Springfield Coal No. 2 (Pennsylvanian; R o 0.54 %) were heated in glass tubes at 80, 100, and 120 °C for 52 months. Sample aliquots from the Upper Cretaceous Second White Specks Formation (R o 0.42 %) were heated in gold tubes at 80 and 100 °C for 42 months at elevated hydrostatic pressures of 100 to 300 MPa. The product gases — methane (CH 4 ) and carbon dioxide (CO 2 ) — were collected and quantified, and gas yields were corrected for leakage from imperfectly closed pores in samples during heating. The results show that longer heating produced more CH 4 . The average CH 4 yields from New Albany Shale and Springfield Coal No. 2 are 0.47 and 3.0 μmol CH 4 per gram of total organic carbon (TOC) over 52 months of heating. Elevated hydrostatic pressure caused lower CH 4 yields from the Second White Specks Formation (3.82 to 1.20 μmol g –1 TOC), suggesting that pressure can retard methanogenesis. Maceral type critically controls the methanogenesis potential of low-maturity coal and shale source rocks. Finally, results of this study provide important insights to the origin of natural gas in low-maturity sedimentary basins.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Geology and Hydrothermal Alteration of the Raft River Geothermal System, Idaho

The Raft River geothermal system is located in southern Idaho, near the Utah-Idaho state boarder in the Raft River Valley. The field, which is owned and operated by U.S. Geothermal, has been selected as an EGS demonstration site by the U. S. Department of Energy. This paper summarizes ongoing geologic and petrologic investigations being conducted in support of this project. The reservoir is developed in fractured Proterozoic schist and quartzite, and Archean quartz monzonite cut by younger diabase intrusions. The basement complex was deformed during the mid Tertiary and covered by approximately 5000 ft of late Tertiary sedimentary and volcanic deposits. Listric normal faults of Cenozoic age disrupt the Tertiary deposits but do not offset the basement rocks. RRG-9, the target well, was drilled southwest of the main well field to a measured depth (MD) of 6089 ft. The well is deviated to the west and cased to a depth of 2316 ft MD. It penetrated the Proterozoic reservoir rocks at a depth of 5286 ft MD and encountered a maximum temperature of 139 degrees C. Bottomhole temperatures in other deep wells range from 133 to 149 degrees C. X-ray diffraction and thin section analyses are being conducted on samples from 5 deep wells, RRG- 1, 2, 3, 7 and 9, to determine the characteristics of the rock types and hydrothermal alteration within the geothermal system. Thin section analyses of samples from RRG-9 document the presence of strong alteration and brecciation at the contact between the Tertiary and basement complex. The Tertiary rocks consist of ash-flow tuffs, lava flows, tuffaceous siltstone, greywacke, and sandstone. No core is available from RRG-9 but core was obtained from RRG-3C. The sample is a brecciated and altered siltstone from the base of the Tertiary sequence and is similar to rocks at the base of the Tertiary deposits in RRG-9. Here, the results of thermal and quasi-static mechanical property measurements that were conducted on the core sample are presented.

Jones, Clay↗
DOE OSTI2025

The influence of kinematics of deformation on polycrystalline halite dynamic recrystallization: Full-field simulation of simple shear versus pure shear

Rock salt, composed mainly of halite, flows viscoplastically over a wide range of geological conditions, strongly impacting the dynamic evolution of sedimentary basins and orogens. Understanding how dislocation creep, which involves dislocation glide, intracrystalline recovery and dynamic recrystallization, influences the microstructure and rheology of halite under various deformation kinematics and temperatures is crucial for enhancing knowledge of salt flow dynamics. Here, this study employs a full-field numerical simulation method to compare the viscoplastic deformation of polycrystalline halite under simple shear and pure shear conditions up to a natural strain of ε = 1.5 at temperatures ranging from 100 °C to 300 °C. The results are presented in terms of crystallographic preferred orientation (CPO), grain shape preferred orientation (SPO), subgrain boundary direction, grain size and strain rate distribution. The results indicate that the crystallographic anisotropy of individual halite crystals is transferred to the polycrystalline scale, resulting in strain localization, particularly in simple shear simulations. The kinematics of deformation affect the evolution and distribution of high strain-rate bands, determining the direction of intragranular substructures and the morphology of strain-induced grain boundaries, with minimal impact on grain size. The intensity of grain boundary migration increases with temperature, significantly influencing grain morphology and size, thereby obscuring strain localization, while having little effect on CPOs. At low strain (ε < 1.0), CPOs relative to both the maximum shortening direction and the grain SPO are similar regardless of the deformation kinematics. At high strain (ε > 1.0), simple shear CPOs exhibit three stronger {100} maxima with a monoclinic symmetry relative to the grain SPO compared to the six {100} maxima with an orthotropic symmetry relative to the grain SPO generated under pure shear. Therefore, microstructures and CPOs can serve as indicators of the strain path in polycrystalline halite under various conditions, aiding in determining the shear sense and elucidating the deformation kinematics of salt structures.

58 GEOSCIENCES↗
DOE OSTI2023

Organic matter sulfurization and organic carbon burial in the Mesoproterozoic

Throughout the Proterozoic Era, sedimentary organic carbon burial helped set the pace of global oxygenation and acted as a major modulator of atmospheric CO 2 and climate. Although Proterozoic rocks generally contain low concentrations of organic matter (OM), there are key exceptions to this rule, including the relatively OM-rich Arctic Bay shales from Baffin Island, Canada (Bylot Supergroup, Borden Basin, ~1.05 Ga). The mechanisms driving elevated OM concentrations in these and other Proterozoic shales remain poorly understood. In the Mesozoic and Cenozoic, organic matter sulfurization can be a major driver of enhanced OM burial across a range of redox conditions comparable to those inferred for many Proterozoic environments. Therefore, in this study, we evaluate the role of sulfurization in driving OM preservation in the Mesoproterozoic Borden Basin and discuss its relevance to Proterozoic systems in general. We present the first evidence for syngenetic-to-early-diagenetic OM sulfurization in a Proterozoic basin, which begins to fill a several-billion-year gap in our record of organic S across Earth history. We find that OM sulfurization was particularly extensive in shales from a relatively shallow-water section (Alpha River) but less extensive in shales deposited in deeper water (Shale Valley), which is consistent with models that infer sulfidic ‘wedges’ or O 2 -minimum-zone-type structures on shelf margins at least intermittently at this time. At the shallower site, organic S and pyrite are similarly 34 S-depleted and thus likely formed at roughly the same time near the sediment–water interface under conditions previously interpreted to have been ferruginous to intermittently sulfidic. In contrast, at the deeper-water site, large S-isotope differences between pyrite and organic S along with low apparent OM sulfurization intensities indicate that pyrite formation was favored over OM sulfurization during early sedimentation under variable but primarily ferruginous conditions. Although Mesoproterozoic biomass can be substantially sulfurized, indicators of sulfurization intensity are not correlated with OM concentrations, and therefore sulfurization does not appear to have been the primary driver of enhanced OM concentrations in Arctic Bay Formation shales. The link between sulfurization and total OM preservation may have been modulated during the deposition of Arctic Bay Formation shales by differences in iron availability, nutrient cycling, and particle dynamics in the Mesoproterozoic.

58 GEOSCIENCES↗
DOE OSTI2025

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES↗
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2026

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗