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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2020

Paleotopography controls weathering of Cambrian-age profiles beneath the Great Unconformity, St. Francois Mountains, SE Missouri, U.S.A.

ABSTRACT Precambrian (1.4 to 1.5 Ga) granite and rhyolite in the St. Francois Mountains at the northeast corner of the Ozark Plateau in Missouri has been altered down to a depth up to 8 meters below the Great Unconformity (the contact between Paleozoic sedimentary rock and underlying Precambrian). Petrographic, geochemical, and mineralogic data indicate that at least two events generated this alteration: 1) surficial weathering due to subaerial exposure of the granite before Cambrian burial—this material is preserved as a paleosol; and 2) alteration due to reaction with basinal fluids channeled along the unconformity from nearby sedimentary basins long after burial by Paleozoic strata. To assess the variation between surficial weathering and basinal fluid alteration, we measured and sampled for petrologic, geochemical, and mineralogic data in the rock at and just below the Great Unconformity at three paleoelevations. Whole-rock geochemical oxide and X-ray diffraction data indicate that K-metasomatism and highly crystalline illite occurred in each profile. The K increase reflects crystallization of authigenic feldspar and illite from basinal fluids that were channeled along the Great Unconformity during younger Paleozoic fluid-flow events. Each profile also exhibits an upward increase in altered feldspar crystals and highly crystalline kaolinite, and an upward decrease in Ca and Na. Such changes reflect soil formation due to reaction with meteoric water before Cambrian burial, indicating that the altered granite was a paleosol before Paleozoic basinal fluid-flow events. Notably, the paleosol at the highest paleoelevation displays the least amount of paleoweathering and the paleosol at the lowest displays the greatest amount of paleoweathering. These results demonstrate that not only can characteristics of the paleosol just below the Great Unconformity be recognized in the St. Francois Mountains, despite subsequent alteration, but also it is possible to detect variations in soil thickness that were controlled by slope steepness and, therefore, water availability and/or soil creep or failure. This spatial relationship is compatible with studies of modern soils which indicate that soil character varies with position on a slope.

Geology↗
DOE OSTI2023

RECLAIM: Electrochemical Lithium and Nickel Extraction with Concurrent Carbon Dioxide Mineralization

This project aims to develop technologies toward the production of Li and Ni from CO 2 -reactive minerals and rocks that contain Ca and Mg, while sequestering CO 2 . The technologies include acoustic stimulation and electrolytic proton production, which use electricity to liberate valuable metal ions from the surrounding mineral matrix at sub-boiling temperatures (~20-80°C). Feedstocks include Li/Ni/Ca/Mg-rich igneous and sedimentary rocks and minerals such as holmquistite, hectorite, lateritic Ni ores enriched in Mg-rich serpentine. Other examples of feedstocks that could be used include coarse grained orthocumulate, mine tailings, etc. In the first step, the electrolytic breakdown of water to H + and OH – (as well as O 2 and H 2 ) will be induced by applying an electric potential to simultaneously produce acidity and alkalinity. The solid feedstocks will be dissolved in the acidic anolyte under acoustic stimulation. Solubilized species will then undergo separation using our unique Li + -selective membranes driven by electrodialysis, leading to concentration of Li as soluble lithium silicate, from which Li 2 CO 3 can be derived by reaction with CO 2(g) while supplying cathode alkalinity. Also at the cathode, stepwise electrolytic precipitation of metal hydroxides, including Ni(OH) 2 , Ca(OH) 2 , and Mg(OH) 2 is performed along a cascade reactor in which both the applied potential and pH increase downstream, thereby simultaneously extracting and separating the metals. Finally, exposure of Ca(OH) 2 and Mg(OH) 2 to dilute CO 2 will result in CO 2 uptake in the form of Ca and Mg carbonates. The primary accomplishments include sample procurement and characterization. Progress has been made on precursor identification and mapping, establishing protocols for assessing metal extractability under sonication, and preliminary technoeconomic analysis (TEA) and life cycle analysis (LCA).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fluid Induced Elemental and Mineralogy Alterations of Caney Shale

ABSTRACT: Studies were conducted on cores and rock cuttings recovered from two wells drilled in Caney shale in the Southern Oklahoma. Mineralogical compositions of these samples were obtained by X-Ray Diffraction (XRD) measurements using two methods of sample preparation: conventional powder (by crushing and grinding of rock fragments) and spray-drying of micronized rock powder method. Microstructure and microanalysis were acquired by Scanning Electron Microscopy/Energy Dispersive Spectroscopy (SEM/EDS). Batch experiments were then conducted using selected rock samples and fracturing fluid whilst deionized water was used as control sample for fluid/rock interaction. These were conducted at 95°C and ambient pressure for various time periods to assess the geochemical reactivity of the Caney shale. The rock samples were subsequently recovered from fluid, dried and analysis repeated using identical XRD and SEM/EDS measurements, following the same conditions. Preliminary experimental results show rock mineralogical compositions are predominantly quartz with relatively high clay content. The clay minerals are mainly illite. Carbonate content is less than 10% by weight for most samples, with other non-clay minerals pyrite, albite, dolomite, and some apatite. In general, sprayed samples showed more intense peaks and less full width at half maxima (FWHM) relative to conventional rock-powder samples. Post-experimental mineralogical changes were observed in samples with clay portions breaking down to amorphous entities. As expected, dissolution of carbonates (calcite and dolomite) was observed and newly detected dissolution of biogenic micro quartz resulted in formation of amorphous silica. 1.0 INTRODUCTION Formation of sedimentary rocks involves processes where rock sediments are emplaced with fluids starting interactions that persist over long periods. During the lithification process, sediments and fluids react and attain equilibrium state. This equilibrium is disturbed during drilling, completion, hydraulic fracturing and other operations that involve the introduction of engineered fluids into the subsurface. To return to equilibrium, dissolution and precipitation of minerals may occur. These reactions have implications on the petrophysical properties of the reservoir. Whiles dissolution and precipitation of minerals can enhance porosity and permeability, the reverse can also be true. For example, when acidic fluids are introduced into formations with high carbonate contents, carbonate minerals dissolve into solution thus increasing the porosity and permeability of the formation. In contrast, when alkaline fluids with pH above 11 are introduced into an illite-rich shale formation, the rock-fluid interaction leads to de-flocculation of illite platelets causing substantial increase in clay fines migration, consequently significant reduction in permeability (Salles et al., 2008).

Awejori, G.↗
DOE OSTI2023

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗
DOE OSTI2020

Shifting modes of iron sulfidization at the onset of OAE-2 drive regional shifts in pyrite δ 34 S records

We report total reduced inorganic S isotope ratios (δ 34 S CRS ) shift toward more negative values across much of the southern North Atlantic just before the onset of the Cenomanian-Turonian Ocean Anoxic Event (OAE-2). At the same time, there is no parallel isotopic change in the significantly larger pool of kerogen (organic) S, which indicates that the distribution and S-isotope composition of sulfide in the environment likely did not drive the change in δ 34 S CRS . Here, we investigate possible explanations for the negative shift in δ 34 S CRS values and their divergence from organic S by isolating iron sulfides for morphological identification and grain-specific isotopic analysis using secondary ion mass spectrometry (SIMS). In pre- and syn-OAE-2 sedimentary rocks from Demerara Rise, we find four distinct morphologies of iron sulfides: pyrite framboids (1–20 μm diameter), irregular pyrite aggregates (1–38 μm diameter), large cemented pyrite aggregates (~60 μm diameter), and irregular and cemented aggregates of the pyrite polymorph marcasite (1–45 μm diameter). These different textural groups have distinct S-isotopic compositions that are largely consistent through the onset of OAE-2. As such, the secular change in bulk δ 34 S CRS values likely reflects the changing proportions of these phases stratigraphically across OAE-2. All textural groups feature resolvable intra-grain δ 34 S variability, suggesting that the environments in which they formed were characterized by dynamic sulfide δ 34 S values and/or by partial closed-system distillation. We use grain-specific δ 34 S distributions to rule out shoaling of the chemocline within the sediments as a mechanism for the observed decrease in δ 34 S CRS . Instead, we propose that changes in the reactivity of the iron species delivered to Demerara Rise over the ~200 kyr leading up to the onset of OAE-2 impacted the relative contributions of pyrite with S-isotope signatures reflecting the water column, shallow sediments, and deeper sediments to the bulk sedimentary δ 34 S CRS value. Specifically, the change in iron reactivity at the onset of OAE-2 favored the production of 34 S-depleted large, cemented aggregates and framboids at the expense of more 34 S-enriched irregular aggregates. Our results underscore that bulk δ 34 S CRS measurements integrate multiple reduced phases that form via distinct reaction mechanisms and potentially in different parts of the depositional environment. Grain-specific SIMS analyses dramatically enrich our ability to interpret pyrite isotopic patterns in the geologic record.

58 GEOSCIENCES↗
DOE OSTI2023

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗
DOE OSTI2025

High-resolution analysis of clay minerals and amorphous materials in martian analog environments

Mudrocks and mud-sized sediments (i.e., silt to clay) dominate the surface of Earth and Mars. These fine-grained sediments preserve a rich history of sedimentary processes from source to sink and shed light on ancient climates. However, both the physical and chemical nature of these materials make them difficult to fully characterize with traditional laboratory techniques. Here, we explore a cross-disciplinary and high-resolution approach using synchrotron radiation for X-ray diffraction, pair distribution function analysis, and submicrometer-scale X-ray fluorescence, combined with transmission electron microscopy, to better understand the nanostructure and composition of mud-sized sediments from a glacio-fluvial watershed in southwest Iceland. Our results demonstrate that sediments in the cold and wet climate of Iceland are more altered than previously thought, as evidenced by the identification of kaolinite and mixed-layer kaolinite-smectite. Additionally, sediments are enriched in amorphous materials and nanocrystalline phases, as determined from grain morphologies and compositions consistent with allophane, hisingerite, ferrihydrite, and halloysite. These alteration products are present as intimate mixtures that vary across depositional sites, demonstrating the dynamic nature of the secondary assemblage from source to sink. This work has implications for Mars, where, for example, basalt-sourced sedimentary rocks from Gale crater are abundant in clay minerals and amorphous materials. Finally, this work underpins the importance of using high-resolution techniques, a coordinated methodology, and developing innovative approaches for future planetary sample return missions (e.g., Mars sample return).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental evaluation of the role of redox during glauconite-CO 2 -brine interactions

Greensands formations are globally abundant sedimentary rocks rich in Fe clays (typically glauconite) that commonly contain natural hydrocarbon accumulations and may be important reservoirs for geologic storage of anthropogenic CO 2 . Diagenesis in greensands is commonly accompanied by the conversion of primary glauconite to siderite (FeCO 3 ), a process that could be exploited for the permanent trapping of CO 2 . Importantly, siderite formation after glauconite requires that the mostly oxidized Fe in the primary Fe clay minerals is reduced during diagenetic interactions. Here, we explore the effect of solution redox state on the stability of glauconite in sandstones with implications for the diagenetic and/or engineered formation of siderite. We performed two flow-through experiments on intact, glauconite-rich sandstone cores at 150 °C and 150 bar. Both experiments employed a 1 mol NaCl/kg, 0.1 mol NaHCO 3 /kg solution charged with ~0.58 mol CO 2 /kg solution, but the redox state of the injected fluid was manipulated between experiments in order to compare glauconite reactivity and siderite saturation state at oxidizing and reducing end-member conditions. After reaction with the oxidizing (O 2 (aq) ≈ 6 μmol/kg) fluid, chemical and Mӧssbauer spectroscopic analyses indicate the production of Fe(III)-oxy/hydroxide minerals from glauconite, whereas, in the reducing (H 2 (aq) ≈ 5–40 mmol/kg) experiment, thermodynamic calculations and coupled chemical, mineralogical, and Mӧssbauer analyses suggest glauconite dissolution and precipitation of an Fe(II) mineral, likely siderite, and minor magnetite formation. These experimental results, along with thermodynamic calculations, confirm that solution redox state is the master variable dictating siderite formation in greensands.

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

An inelastic homogenization framework for layered materials with planes of weakness

Many geologic materials have a composite structure, in which macroscopic mechanical behavior is determined by the properties, shape, and heterogeneous distribution of individual constituents. In particular, sedimentary rocks commonly exhibit a layered microstructure, with distinct bedding planes that can also form planes of weakness. In this work, we present a homogenization framework for modeling inelastic layered media. The proposed constitutive model allows for distinct micro-constitutive laws for each layer, explicit representation of layer distributions, as well as incorporation of imperfect bonding at the interface between adjacent layers. No a priori assumptions are needed regarding the specific constitutive models used for the layers and interfaces, providing significant modeling flexibility. The overall framework provides a simple and physically-motivated way of defining anisotropic material behavior as an emergent property of the layered microstructure. The model is calibrated using triaxial and true-triaxial experimental data to demonstrate its ability to describe anisotropic deformation and multiple modes of failure.

58 GEOSCIENCES↗
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 OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2020

U-Pb zircon age constraints on the earliest eruptions of the Deccan Large Igneous Province, Malwa Plateau, India

Climate instability driven by emission of volatiles during emplacement of large igneous provinces (LIPs) is frequently invoked as a potential cause of mass extinctions. However, documenting this process in the geologic record requires a holistic understanding of eruption rates, the location of eruptive centers, and potential sources of climate-changing volatiles. We present new chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology from Malwa Plateau basalts on the northern margin of the Deccan LIP, India. These basalts have been previously interpreted as either an extension of the province's main volcanic stratigraphy or as an independent eruptive center active up to millions of years prior to the main eruptive phase. Our data instead demonstrate that the lower Malwa Plateau basalts are temporally correlative with the first pulse of Deccan volcanism and provide new constraints on its initiation and duration. Paleomagnetic data further indicate that upper Malwa Plateau basalts may be age-equivalent to the second, third, and fourth pulses of Deccan volcanism. The relative thicknesses of age-equivalent packages of basalt are consistent with eruption of the Deccan LIP from a southward-migrating eruptive center. The first eruptive pulse is coeval with a ~200 kyr Late Maastrichtian warming event preserved globally in contemporaneous stratigraphic sections. We propose that the first pulse of Deccan magmatism was more voluminous in the north, where it erupted through organic-rich sedimentary rocks of the Narmada-Tapti rift basin. Thermal metamorphism of these sediments could have been a source of sufficient CO 2 to drive the Late Maastrichtian warming event, which if true would reconcile the apparent dampened warming signals associated with later Deccan eruptive pulses.

58 GEOSCIENCES↗
DOE OSTI2025

Elevated-Mn ChemCam targets illuminating Mn redox cycling and diagenesis in the Bradbury Rise, Gale Crater, Mars

Manganese plays a crucial role as a paleo-environmental and geological indicator due to its sensitivity to redox potential and pH variations in the environment. On Earth, the association between the rise of atmospheric oxygen during the Great Oxidation Event and the presence of Mn in the sedimentary rock record underscores its significance. Here, in this study, we reexamined ChemCam targets from the first 600 sols of the Mars Science Laboratory mission, focusing on identifying instances of above-average Mn within these targets. These elevated-Mn targets were categorized into distinct geologic classes, revealing a pattern linking heightened Mn levels with diagenetically altered materials, such as calcium-sulfate veins and concretions, as well as clay minerals within the same targets, indicating a compelling relationship between Mn enrichment and diagenetic processes. High concentrations of Mn were observed in chemically altered targets, suggesting the occurrence of multiple fluid events: the first to alter the material and the second to deposit Mn. The observed patterns suggest multiple diagenetic events and redox cycling that facilitated the deposition and transport of Mn subsequent to the initial dissolution of basaltic materials. This research sheds light on the complexity of martian diagenetic processes and their implications for the planet’s environmental evolution.

58 GEOSCIENCES↗
DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2023

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗