Results for “STRATIGRAPHY”

DOE OSTI2023

Potential for constraining sequence stratigraphy and cycle stratigraphy with U-Pb dating of carbonates

Studies of sequence stratigraphy and cycle stratigraphy in carbonate formations increasingly improve our understanding of depositional systems and basinal to global correlations. However, there is still a scarcity of dated sequences to determine the time length of sequences and cycles. Furthermore, the durations of sequences and cycles provide tests of the processes that produce them. Poor age resolution mostly stems from the lack of recognizable volcanic ash deposits in many sections. Through continuing advances in analytical techniques and approaches to sampling, the current state of Usingle bondPb dating of carbonate sedimentary rocks provides opportunities to fill in age gaps. The approach should build on well-constrained 3-dimensional stratigraphic architectures and diagenetic stratigraphies.

58 GEOSCIENCES↗

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 OSTI2021

The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings

The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.

58 GEOSCIENCES↗

DOE OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗

DOE OSTI2024

Sensitivity of modelled passive margin stratigraphy to variations in sea level, sediment supply and subsidence

Abstract We produced a 10 Myr synthetic stratigraphic section using a forward stratigraphic model that generates marine deltaic stratigraphy over geological timescales. We recursively fit the model using a Bayesian inversion algorithm to test: (1) if it could be accurately reconstructed; (2) if the parameters used to create it could be recovered; and (3) the sensitivity of the model output to given model parameters and the attendant physical processes. The original synthetic stratigraphic section was produced with cyclical sea‐level variations of 40 and 30 m with 2.4 and 10 Myr periods respectively. Sediment was also supplied cyclically, in 2.4 and 10 Myr cycles with amplitudes of 30 and 80 tons/100 kyr, respectively, varying from a mean of 232 tons/100 kyr. Parameter values were sampled to fit the model using a Markov chain Monte Carlo algorithm, resulting in a ±5 m (1σ) variation between the experimental output and the original. Sea level varied by ±7 m (1σ) within the posterior distribution of parameters. As a result, both the 10 Myr and 2.4 Myr sea‐level cycles could be extracted from the original output. The variation in sediment supply was approximately ±38 tons/100 kyr (1σ) and, as a result, only the larger long‐term supply variations could be accurately recovered in refitting the model. The variation in thermal, flexural and total subsidence across those parameter sets is less than ±10 m (1σ). The original section experienced 150 m of total subsidence at the depocentre. Our results demonstrate the distinct and interpretable imprint of sea level and subsidence on continental margin stratigraphy can be quantified. Moreover, we conclude that sea‐level change produces a defined effect on the geometries of stratigraphic architecture, and that techniques applied for the purpose of delineating sea‐level variation from continental margin strata have a well‐founded conceptual basis.

Schmelz, William J.↗

DOE OSTI2025

Snow ALbedo eVOlution (SALVO) Campaign Qualitative Characterization of Snowpack Stratigraphy from April - June, 2024 in Utqiagivk, AK

The physical structure of the snowpack is impacted by and impacts the exchange of energy with the atmosphere. Qualitative characterization of snowpack stratigraphy in snow pits provides context for what physical processes are changing the snow (e.g., melt, wind-packing, hoar formation). 36 snowpits were dug as part of the 2024 Snow Albedo Evolution (SALVO) field campaign. Snow grain stratigraphy was recorded using standard International Association of Cryospheric Sciences (IACS) snow grain classification ( https://cryosphericsciences.org/publications/snow-classification/) and hand hardness tests (in order of least hard to hardest: F (fist), 4F (four fingers), 1F (one finger), P (pencil), K (knife), I (ice)). These snowpits were generally co-located with library site albedo measurements and are given the same index (Date, Site, Location, Position) when appropriate (see broadband and spectral albedo data archived from the SALVO II project).

Depth to bottom of snow layer↗

DOE OSTI2022

Supplemental Material: Cretaceous sequence stratigraphy of the northern Baltimore Canyon Trough: Implications for basin evolution and carbon storage

We focus on potential reservoirs in the Logan Canyon Formation. Based on lithostratigraphy, Libby-French originally divided the Logan Canyon sandstones into lower and upper Logan Canyon sandstone units separated by the Sable Shale equivalent. Miller et al. reinterpreted the well-log correlations with biostratigraphy using a sequence stratigraphic approach and found three distinct Logan Canyon highstand sandstones associated with three Aptian to lower Cenomanian stratigraphic sequences, the LC3, LC2, and LC1, respectively. In addition, Miller et al. found that the Sable Shale Formation identified by Libby-French was placed above and below the LC1 basal sequence boundary in some wells, highlighting issues with correlating units based only on lithology. The Upper Cretaceous strata in the BCT are termed the Dawson Canyon Equivalent Formation, correlate to the DCx composite sequence, and cap the Logan Canyon sequences below made up of a very thick, calcareous Cretaceous sequence stratigraphy mudstone, divided by a single, thick, laterally continuous Coniacian to Campanian sandstone called the Middle Sandstone.

58 GEOSCIENCES↗

DOE OSTI2022

Cretaceous sequence stratigraphy of the northern Baltimore Canyon Trough: Implications for basin evolution and carbon storage

We evaluate the Cretaceous stratigraphy and carbon sequestration potential of the northern Baltimore Canyon Trough (NBCT) using >10,000 km of multi-channel seismic profiles integrated with geophysical logs, biostratigraphy, and lithology from 29 offshore wells. We identify and map six sequences resolved primarily at the stage level. Accommodation was dominated by thermal and non-thermal subsidence, though sequence boundaries correlate with regional and global sea-level changes, and the record is modified by igneous intrusion, active faulting, and changes in sediment supply and sources. Our stratigraphic maps illustrate a primary southern (central Appalachian) Early Cretaceous source that migrated northward during the Aptian and Albian. During the Cenomanian, sedimentation rates in the NBCT increased and depocenters shifted northward and landward. We show that deposition occurred in three phases: (1) earlier Cretaceous paleoenvironments were primarily terrestrial indicated by variable amplitude, chaotic seismic facies, serrated gamma logs, and heterolithic sandstones and mudstones with terrestrial microfossils; (2) the Albian to Cenomanian was dominated by deltaic paleoenvironments indicated by blocky, funnel-shaped, gamma-ray logs and clinoforms characterized by continuous high-amplitude seismic reflections with well-defined terminations; and (3) the Cenomanian and younger was marine shelf, inferred from mudstoneprone lithologies, peak gamma-ray values in well logs, and foraminiferal evidence. Long-term transgression and maximum water depths at the Cenomanian/Turonian boundary correlative with Ocean Anoxic Event 2 were followed by a regression and relative sea-level fall. We show that porous and permeable sandstones of three Aptian to Cenomanian highstand systems tracts are high-volume reservoirs for supercritical CO 2 storage that are confined by overlying deep water mudstones.

58 GEOSCIENCES↗

DOE OSTI2021

Evaluating the Thickness and Stratigraphy of Ejecta Materials at the Chang’e-4 Landing Site

The landing site of China’s Chang’e-4 (CE-4) probe is located on the mare basalts on the floor of the Von Kármán crater on the lunar far side. The Von Kármán crater is inside the ancient and highly cratered South Pole–Aitken basin, which has experienced complex emplacement sequences of both near and distant ejecta materials. These issues complicated the interpretation of the CE-4 surface in situ measurements of the visible and near-infrared spectrometer and the lunar penetrating radar (LPR) onboard Yutu-2 rover. To evaluate the sources and amounts of all principle foreign materials at the CE-4 landing site, we thoroughly examine the ejecta delivered by crater-forming events that occurred later than the formation of the mare basalts at the CE-4 landing site. We found a total of 16 craters that may have delivered ejecta thicker than 10 cm level superposed on the mare basalts at the CE-4 landing site. Crater Finsen, Von Kármán L, Von Kármán L′, and Maksutov are the top four major foreign material sources, and each of them contributed ejecta thicker than 1 m. Our surveys confirm that the ejecta from Finsen crater are the most dominant foreign materials in the uppermost few meters at the CE-4 landing site and the total impact ejecta deposited upon the mare basalts at the landing site is estimated to be thinner than 30 m. We found that the estimations from Pike’s model are the most consistent with the Yutu-2 LPR observations.

58 GEOSCIENCES↗

DOE OSTI2020

Sedimentology and stratigraphy of a deep-water gas hydrate reservoir in the northern Gulf of Mexico

Not provided.

Geology↗

DOE OSTI2021

Deposition, Diagenesis, and Sequence Stratigraphy of the Pennsylvanian Morrowan and Atokan Intervals at Farnsworth Unit

Farnsworth Field Unit (FWU), a mature oilfield currently undergoing CO2-enhanced oil recovery (EOR) in the northeastern Texas panhandle, is the study area for an extensive project undertaken by the Southwest Regional Partnership on Carbon Sequestration (SWP). SWP is characterizing the field and monitoring and modeling injection and fluid flow processes with the intent of verifying storage of CO2 in a timeframe of 100–1000 years. Collection of a large set of data including logs, core, and 3D geophysical data has allowed us to build a detailed reservoir model that is well-grounded in observations from the field. This paper presents a geological description of the rocks comprising the reservoir that is a target for both oil production and CO2 storage, as well as the overlying units that make up the primary and secondary seals. Core descriptions and petrographic analyses were used to determine depositional setting, general lithofacies, and a diagenetic sequence for reservoir and caprock at FWU. The reservoir is in the Pennsylvanian-aged Morrow B sandstone, an incised valley fluvial deposit that is encased within marine shales. The Morrow B exhibits several lithofacies with distinct appearance as well as petrophysical characteristics. The lithofacies are typical of incised valley fluvial sequences and vary from a relatively coarse conglomerate base to an upper fine sandstone that grades into the overlying marine-dominated shales and mudstone/limestone cyclical sequences of the Thirteen Finger limestone. Observations ranging from field scale (seismic surveys, well logs) to microscopic (mercury porosimetry, petrographic microscopy, microprobe and isotope data) provide a rich set of data on which we have built our geological and reservoir models.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2024

The dynamics and stratigraphy of distributary channel networks (Final Technical Report)

Sedimentary reservoirs are frequently composed of the deposits of distributary channel networks. These networks are associated with coarse sediments (and high permeability) within the channel network, and relatively fine sediments (and low permeability) outside the channel network. The shape of the network is complex, and a product of the environmental boundary conditions under which it grows (Figure 1). These systems control heterogeneity at small and large scales within sedimentary basins, and understanding this heterogeneity and the links between scales represents the important scientific frontier addressed by this project.

58 GEOSCIENCES↗

DOE OSTI2025

Revealing the Moon's Taurus‐Littrow Landslide via Integrated Analysis of Pristine Apollo 17 Soil Core 73001/2

The “light mantle” deposit at the base of South Massif in the Moon's Taurus‐Littrow Valley was a primary science target for the Apollo 17 exploration. The possibility that it was a landslide triggered by ejecta from Tycho Crater is critical for establishing the age of Tycho and constraining recent lunar impact chronology; however, the mechanism of emplacement of the deposit has recently been questioned. The newly opened 73001/73002 double drive tube from Station 3 sampled 70.6 cm deep into the regolith and represents the first stratigraphic section of an extraterrestrial landslide deposit returned to Earth. Here we provide an overview of the stratigraphy of the 73001/73002 core based on top to bottom variations revealed by coordinated laboratory analyses and explore constraints on the emplacement of the light mantle deposit. Briefly, the upper ∼10 cm of 73002 contains a disturbed zone from space weathering and emplacement of ejecta from a nearby crater that excavated and ejected basaltic material. Below 10 cm is a nearly uniform unit of immature regolith. These data support a single event for the emplacement of the deposit at this location, followed by weathering and mixing of materials from nearby crater ejecta in the upper 10 cm. Slight variations in chemistry and clast components may reflect the relative stratigraphy of the South Massif slope, with material toward the bottom of 73001 originating from lower slopes and material from higher up in the core representing regolith from higher up the South Massif slopes.

58 GEOSCIENCES↗

DOE OSTI2024

INTERDUNE DEPOSITS AND STRATIGRAPHIC FRAMEWORK OF THE NAVAJO SANDSTONE IN THE LAKE POWELL AREA, UTAH, USA

The Western Interior of the USA was dominated by desert environments for much of the Jurassic. Of these desert environments, the most extensive by far was the Navajo Sandstone erg, which, in conjunction with its correlative Nugget and Aztec Sandstones, is thought to have covered well over 500,000 km2, and is one of the largest ergs preserved in the rock record. The Navajo Sandstone is also one of a few geologic units which show promise as a potential CO2 sequestration reservoir in Utah. It outcrops along the eastern and southern margin of the state, and subcrops throughout most of the central area. The Navajo Sandstone consists of stacked eolian deposits, composed of well sorted, clean sandstone, deposited in commonly tabular, high angle cross bedded paleodunes, which are on average 15 m thick. Interspersed between these paleodunes are interdune deposits. Interdune deposits consist of sandy to calcareous, predominantly planar bedded, lensoid units, ranging from <1 m up to 7 m thick. In some cases, these interdune deposits occur along more extensive stratigraphic bounding surfaces. Dating within the Navajo Sandstone has historically been very difficult, as a result of its predominantly clean, aeolian sandstone composition. In addition, complexity of internal stratigraphic surfaces coupled with a low structural dip and generally flat topography has resulted in a poor understanding of the internal stratigraphy of the Navajo Sandstone. Exceptional three-dimensional exposures of the Navajo Sandstone in the upper portions of Lake Powell permit us to relate interdune carbonates to a hierarchical bounding surfaces model. Here we use this combined approach to interpret the internal stratigraphy of the Navajo Sandstone in this region. This work forms part of a larger study of the Navajo Sandstone across Utah and its potential as a CO2 reservoir, by developing a more comprehensive stratigraphic framework, and increasing understanding of stratigraphic complexity within one of the most significant eolian systems in the world.

Mahon, Elizabeth↗

DOE OSTI2024

Timescales of mafic magmatic fractionation documented by paleosecular variation in basaltic drill core, Snake River Plain volcanic province, Idaho, USA

Abstract The timescales over which fractional crystallization and recharge work in mafic volcano-plutonic provinces is subject to great uncertainty. Currently modeled processes are subject to the scale of measurement: monogenetic basaltic fields accumulate over hundreds of thousands of years, consistent with U-Th-Ra isotopic variations that imply 50% crystallization of basic magmas on timescales of 100,000 years or more, whereas crystal diffusion modeling implies phenocryst residence times of ~1–1000 years. Monogenetic basalts of the Snake River Plain in southern Idaho, USA, are up to 2 km thick and postdate passage over the Yellowstone–Snake River Plain hotspot. Detailed lithologic and geophysical logging of core from deep drill holes, along with chemical stratigraphy and high-resolution paleomagnetic inclination measurements, document individual eruptive units, compound lava flows, and basaltic flow groups that accumulated over 1–6 m.y. Hiatuses are commonly marked by loess or fluvial interbeds that vary from ~0.1 m thick to 20 m thick. Radiometric (40Ar-39Ar, detrital zircon U-Pb) and paleomagnetic timescale ages show that the deepest hole (Kimama drill hole, 1912 m total depth) accumulated over ~6 m.y. Cycles of fractional crystallization and recharge are recognized in the chemical stratigraphy as up-section shifts in major and trace elements; these fractionation cycles commonly represent 40%–50% fractionation. Individual fractionation cycles may comprise 20–40 eruptive units (8–17 lava flows) with little to no change in paleomagnetic inclination (0°–1°), whereas adjacent cycles may differ by several degrees from one another or reflect changes in polarity. Rates of paleosecular variation in Holocene lavas and sediments dated using 14C document significant shifts in magnetic inclination over short timescales, ranging from ~0.05° to 2°/decade, with an average of ~0.5°/decade and a minimum rate of 0.05°/decade. This implies that fractionation cycles with ≤1° variation in magnetic inclination formed on timescales of a few decades up to a few centuries (20–200 years). Thus, the lavas collectively represent only a few thousand years of eruptive activity, with major flow groups separated in time by tens to hundreds of thousands of years. We suggest that the rates defined by paleosecular variation capture the timescales of magmatic chamber evolution (fractionation/recharge) in the seismically imaged mid-crustal sill complex; in contrast, we suggest that crystal diffusion modeling captures the residence times in shallow subvolcanic magmatic chambers that underlie individual monogenetic volcanoes.

Geology↗

DOE OSTI2022

Organic carbon burial by river meandering partially offsets bank erosion carbon fluxes in a discontinuous permafrost floodplain

Abstract. Arctic river systems erode permafrost in their banks and mobilize particulate organic carbon (OC). Meandering rivers can entrain particulate OC from permafrost many meters below the depth of annual thaw, potentially enabling the production of greenhouse gases. However, the amount and fate of permafrost OC that is mobilized by river erosion is uncertain. To constrain OC fluxes due to riverbank erosion and deposition, we collected riverbank and floodplain sediment samples along the Koyukuk River, which meanders through discontinuous permafrost in the Yukon River watershed, Alaska, USA, with an average migration rate of 0.52 m yr−1. We measured sediment total OC (TOC) content, radiocarbon activity, water content, bulk density, grain size, and floodplain stratigraphy. Radiocarbon activity and TOC content were higher in samples dominated by silt as compared to sand, which we used to map OC content onto floodplain stratigraphy and estimate carbon fluxes due to river meandering. Results showed that the Koyukuk River erodes and re-deposits a substantial flux of OC each year due to its depth and high migration rate, generating a combined OC flux of a similar magnitude to the floodplain net ecological productivity. However, sediment being eroded from cutbanks and deposited as point bars had similar OC stocks (mean ± 1 SD of 125.3±13.1 kg OC m−2 in cutbanks versus 114.0±15.7 kg OC m−2 in point bars) whether or not the banks contained permafrost. We also observed radiocarbon-depleted biospheric OC in both cutbanks and permafrost-free point bars. These results indicate that a substantial fraction of aged biospheric OC that is liberated from floodplains by bank erosion is subsequently re-deposited in point bars rather than being oxidized. The process of aging, erosion, and re-deposition of floodplain organic material may be intrinsic to river–floodplain dynamics, regardless of permafrost content.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗