Results for “Earth and planetary sciences” — Page 2

DOE OSTI2024

Extensive Secondary Cratering From the InSight Sol 1034a Impact Event

Abstract Impact cratering is one of the fundamental processes throughout the history of the Solar System. The formation of new impact craters on planetary bodies has been observed with repeat images from orbiting satellites. However, the time gap between images is often large enough to preclude detailed analysis of smaller‐scale features such as secondary impact craters, which are often removed or buried over a short time period. Here we use a seismic event detected on Mars by the NASA InSight mission to investigate secondary cratering at a new impact crater. We strengthen the case that the seismic event that occurred on Sol 1034 (S1034a) is the result of a new impact cratering event. Using the exact timing of this event from InSight, we investigated the resulting new impact crater in orbital image data. The S1034a impact crater is approximately 9 m in diameter but is responsible for over 900 secondary impact events in the form of low albedo spots that are located at distances of up to almost 7 km from the primary crater. We suggest that the low albedo spots formed from relatively low energy ejecta, with individual ejecta block velocities less than 200 m s −1 . We estimate that the low albedo spots, the main evidence of secondary impact processes at this new impact event, fade within 200–300 days after formation.

Grindrod, P. M. [Natural History Museum London UK]↗

DOE OSTI2025

Recommendations for Comprehensive and Independent Evaluation of Machine Learning‐Based Earth System Models

Abstract Machine learning (ML) is a revolutionary technology with demonstrable applications across multiple disciplines. Within the Earth science community, ML has been most visible for weather forecasting, producing forecasts that rival modern physics‐based models. Given the importance of deepening our understanding and improving predictions of the Earth system on all time scales, efforts are now underway to develop Earth‐system models (ESMs) capable of representing all components of the coupled Earth system (or their aggregated behavior) and their response to external changes over long timescales. Building trust in ESMs is a much more difficult problem than for weather forecast models, not least because the model must represent the alternate (e.g., future or paleoclimatic) coupled states of the system for which there are no direct observations. Given that the physical principles that enable predictions about the response of the Earth system are often not explicitly coded in these ML‐based models, demonstrating the credibility of ML‐based ESMs thus requires us to build evidence of their consistency with the physical system. To this end, this paper puts forward five recommendations to enhance comprehensive, standardized, and independent evaluation of ML‐based ESMs to strengthen their credibility and promote their wider use.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Implementation of Sub‐Grid Scale Temperature Perturbations Induced by Non‐Orographic Gravity Waves in WACCM6

Abstract Atmospheric gravity waves can play a significant role on atmospheric chemistry through temperature fluctuations. A recent modeling study introduced a method to implement subgrid‐scale orographic gravity‐wave‐induced temperature perturbations in the Whole Atmosphere Community Climate Model (WACCM). The model with a wave‐induced temperature parameterization was able to reproduce for example, the influence of mountain wave events on atmospheric chemistry, as highlighted in previous literature. Here we extend the subgrid‐scale wave‐induced temperature parameterization to also include non‐orographic gravity waves arising from frontal activity and convection. We explore the impact of these waves on middle atmosphere chemistry, particularly focusing on reactions that are strongly sensitive to temperature. The non‐orographic gravity waves increase the variability of chemical reaction rates, especially in the lower mesosphere. As an example, we show that this, in turn, leads to increases in the daytime ozone variability. To demonstrate another impact, we briefly investigate the role of non‐orographic gravity waves in cirrus cloud formation in this model. Consistent with findings from the previous study focusing on orographic gravity waves, non‐orographic waves also enhance homogeneous nucleation and increase cirrus clouds. The updated method used enables the global chemistry‐climate model to account for both orographic and non‐orographic gravity‐wave‐induced subgrid‐scale dynamical perturbations in a consistent manner.

Yook, Simchan [Department of Earth Atmospheric and↗

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 OSTI2025

On the Dependence of Simulated Convection on Domain Size in CRMs

Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.

CRM↗

DOE OSTI2025

Navigating the Noise: Bringing Clarity to ML Parameterization Design With O $\boldsymbol{\mathcal{O}}$(100) Ensembles

Abstract Machine‐learning (ML) parameterizations of subgrid processes (here of turbulence, convection, and radiation) may one day replace conventional parameterizations by emulating high‐resolution physics without the cost of explicit simulation. However, uncertainty about the relationship between offline and online performance (i.e., when integrated with a large‐scale general circulation model) hinders their development. Much of this uncertainty stems from limited sampling of the noisy, emergent effects of upstream ML design decisions on downstream online hybrid simulation. Our work rectifies the sampling issue via the construction of a semi‐automated, end‐to‐end pipeline for size ensembles of hybrid simulations, revealing important nuances in how systematic reductions in offline error manifest in changes to online error and online stability. For example, removing dropout and switching from a Mean Squared Error to a Mean Absolute Error loss both reduce offline error, but they have opposite effects on online error and online stability. Other design decisions, like incorporating memory, converting moisture input from specific humidity to relative humidity, using batch normalization, and training on multiple climates do not come with any such compromises. Finally, we show that ensemble sizes of may be necessary to reliably detect causally relevant differences online. By enabling rapid online experimentation at scale, we can empirically settle debates regarding subgrid ML parameterization design that would have otherwise remained unresolved in the noise.

Lin, Jerry [Department of Earth System Sciences Un↗

DOE OSTI2024

MIBiG 4.0: advancing biosynthetic gene cluster curation through global collaboration

Specialized or secondary metabolites are small molecules of biological origin, often showing potent biological activities with applications in agriculture, engineering and medicine. Usually, the biosynthesis of these natural products is governed by sets of co-regulated and physically clustered genes known as biosynthetic gene clusters (BGCs). To share information about BGCs in a standardized and machine-readable way, the Minimum Information about a Biosynthetic Gene cluster (MIBiG) data standard and repository was initiated in 2015. Since its conception, MIBiG has been regularly updated to expand data coverage and remain up to date with innovations in natural product research. Here, we describe MIBiG version 4.0, an extensive update to the data repository and the underlying data standard. In a massive community annotation effort, 267 contributors performed 8304 edits, creating 557 new entries and modifying 590 existing entries, resulting in a new total of 3059 curated entries in MIBiG. Particular attention was paid to ensuring high data quality, with automated data validation using a newly developed custom submission portal prototype, paired with a novel peer-reviewing model. MIBiG 4.0 also takes steps towards a rolling release model and a broader involvement of the scientific community. MIBiG 4.0 is accessible online at https://mibig.secondarymetabolites.org/.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Studying Aerosol, Clouds, and Air Quality in the Coastal Urban Environment of Southeastern Texas

A multi-agency succession of field campaigns was conducted in southeastern Texas during July 2021 through October 2022 to study the complex interactions of aerosols, clouds and air pollution in the coastal urban environment. As part of the Tracking Aerosol Convection interactions Experiment (TRACER), the TRACER- Air Quality (TAQ) campaign the Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) and the Convective Cloud Urban Boundary Layer Experiment (CUBE), a combination of ground-based supersites and mobile laboratories, shipborne measurements and aircraft-based instrumentation were deployed. These diverse platforms collected high-resolution data to characterize the aerosol microphysics and chemistry, cloud and precipitation micro- and macro-physical properties, environmental thermodynamics and air quality-relevant constituents that are being used in follow-on analysis and modeling activities. We present the overall deployment setups, a summary of the campaign conditions and a sampling of early research results related to: (a) aerosol precursors in the urban environment, (b) influences of local meteorology on air pollution, (c) detailed observations of the sea breeze circulation, (d) retrieved supersaturation in convective updrafts, (e) characterizing the convective updraft lifecycle, (f) variability in lightning characteristics of convective storms and (g) urban influences on surface energy fluxes. The work concludes with discussion of future research activities highlighted by the TRACER model-intercomparison project to explore the representation of aerosol-convective interactions in high-resolution simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Using Radiogenic Noble Gas Nuclides to Identify and Characterize Rock Fracturing

Abstract Fracture‐released radiogenic noble gas nuclides are used to identify locations and constrain the volume of new fracture creation during subsurface detonations. Real‐time, in situ noble gases and reactive gases were monitored using a field‐deployed mass spectrometer and automated sampling system in a multilevel borehole array. Released gases were measured after two different detonations having distinct energy, pressure, and gas volume characteristics. Explosive‐derived gases (N 2 O, CO 2 ) and excess radiogenic 4 He and 40 Ar above atmospheric background are used to identify locations of gas transport and new fracture creation after each detonation. Fracture‐released radiogenic 4 He is used to constrain the volume of newly created fractures with a model of helium release from fracturing. Explosive by‐product gas was observed in multiple locations both near and distal to the shot locations for both detonations. Radiogenic 4 He and 40 Ar release from rock damage was observed in locations near the detonation after the second, more powerful detonation. Observed 4 He response is consistent with a model of diffusive release from newly created fractures. Volume of new fractures estimated from the 4 He release ranges from 1 to 5 m 2 with apertures ranging from 0.1 to 1 m. Our results provide evidence that radiogenic noble gases released during fracture creation can be identified at the field scale in real time and used to identify timing and location of fracture creation during deformation events. This technique could be useful in subsurface science and engineering problems where the location and amount of newly created rock fracturing is of interest including fault rupture, mine safety, subsurface detonation monitoring and reservoir stimulation.

58 GEOSCIENCES↗

DOE OSTI2022

A Non-perturbative Approach to Computing Seismic Normal Modes in Rotating Planets

In this work, a continuous Galerkin method based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial filtering eigensolver. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid regions. The relevant eigenvalue problem is solved by a combination of several highly parallel and computationally efficient methods. We validate our three-dimensional results in the non-rotating case using analytical results for constant elastic balls, as well as numerical results for an isotropic Earth model from standard “radial” algorithms. We also validate the computations in the rotating case, but only in the slowly-rotating regime where perturbation theory applies, because no other independent algorithms are available in the general case. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Earth and Mars models studying the effects of heterogeneity on a large range of scales.

58 GEOSCIENCES↗

DOE OSTI2025

Molecular insights into the composition, sources, and aging of atmospheric brown carbon

The light-absorbing chemical components of atmospheric organic aerosols are commonly referred to as Brown Carbon (BrC), reflecting the characteristic yellowish to brown appearance of aerosol.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Competition response of cloud supersaturation explains diminished Twomey effect for smoky aerosol in the tropical Atlantic

The Twomey effect brightens clouds by increasing aerosol concentrations, which activates more droplets and decreases cloud supersaturation in response to more competition for water vapor. To quantify this competition response, we used marine low cloud observations in clean and smoky conditions at Ascension Island in the tropical South Atlantic during the Layered Aerosol Smoke Interactions with Cloud (LASIC) campaign. These observations show similar increases in droplet number for increased accumulation-mode particles from surface-based and satellite cloud retrievals, demonstrating the importance of below-cloud aerosol measurements for retrieving aerosol–cloud interactions (ACI) in clean and smoky aerosol conditions. Four methods for estimating cloud supersaturation from aerosol–cloud measurements were compared, with cloud scene-based and parcel-based methods showing sufficient variability for a strong dependence on both aerosol accumulation number concentration and cloud-base updraft velocities. Decomposing aerosol-related changes in cloud albedo and optical depth shows the calculated competition response accounts for dampening the activation response by 12 to 35%, explaining the diminished Twomey effect at high aerosol concentrations observed for smoky conditions at LASIC and previously around the world. This result was consistent for independent supersaturation retrievals by cloud scene-based droplet number and cloud condensation nuclei and parcel-based multimode size-resolving Lagrangian methods. Translating aerosol effects to local radiative forcing with clean conditions as a proxy for preindustrial and smoky conditions for present-day showed that the competition response reduces cooling from the Twomey radiative forcing by 12 to 35%, providing an essential process-specific constraint for improving the representation of aerosol competition in climate model simulation of indirect aerosol forcing.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Planetary normal mode computation: Parallel algorithms, performance and reproducibility

This report is an extension of work entitled “Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.” by Shi et al., [1] originally presented at the SC18 conference. A highly parallel polynomial filtered eigensolver was developed and exploited to calculate the planetary normal modes. The proposed method is ideally suited for computing interior eigenpairs for large-scale eigenvalue problems as it greatly enhances memory and computational efficiency. In this article, the second-order finite element method is used to further improve the accuracy as only the first-order finite element method was deployed in the previous work. The parallel algorithm, its parallel performance up to 20k processors, and the great computational accuracy are illustrated. The reproducibility of the previous work was successfully performed on the Student Cluster Competition at the SC19 conference by several participant teams using a completely different Mars-model dataset on different clusters. Both weak and strong scaling performances of the reproducibility by the participant teams were impressive and encouraging. The analysis and reflection of their results are demonstrated and future direction is discussed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index

Despite its importance for determining global carbon fluxes, leaf respiration remains poorly constrained in land surface models (LSMs). We tested the sensitivity of the Energy Exascale Earth System Model Land Model – Functionally Assembled Terrestrial Ecosystem Simulator (ELM-FATES) to variation in the canopy gradients of leaf maintenance respiration (R dark ). We ran global and point simulations varying the canopy gradient of R dark to explore the impacts on forest structure, composition, and carbon cycling. In global simulations, steeper canopy gradients of R dark lead to increased understory survival and leaf biomass. Leaf area index (LAI) increased up to 77% in tropical regions compared with the default parameterization, improving alignment with remotely sensed benchmarks. Global vegetation carbon varied from 308 Pg C to 449 Pg C across the ensemble. In tropical forest simulations, steeper gradients of R dark had a large impact on successional dynamics. Results show the importance of canopy gradients in leaf traits and fluxes for determining plant carbon budgets and emergent ecosystem properties such as competitive dynamics, LAI, and vegetation carbon. The high-model sensitivity to canopy gradients in R dark highlights the need for more observations of how leaf traits and fluxes vary along light micro-environments to inform critical dynamics in LSMs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Achievement of the Planetary Defense Investigations of the Double Asteroid Redirection Test (DART) Mission

NASA's Double Asteroid Redirection Test (DART) mission was the first to demonstrate asteroid deflection, and the mission's Level 1 requirements guided its planetary defense investigations. Here, we summarize DART's achievement of those requirements. On 2022 September 26, the DART spacecraft impacted Dimorphos, the secondary member of the Didymos near-Earth asteroid binary system, demonstrating an autonomously navigated kinetic impact into an asteroid with limited prior knowledge for planetary defense. Months of subsequent Earth-based observations showed that the binary orbital period was changed by –33.24 minutes, with two independent analysis methods each reporting a 1 σ uncertainty of 1.4 s. Dynamical models determined that the momentum enhancement factor, β , resulting from DART's kinetic impact test is between 2.4 and 4.9, depending on the mass of Dimorphos, which remains the largest source of uncertainty. Over five dozen telescopes across the globe and in space, along with the Light Italian CubeSat for Imaging of Asteroids, have contributed to DART's investigations. These combined investigations have addressed topics related to the ejecta, dynamics, impact event, and properties of both asteroids in the binary system. A year following DART's successful impact into Dimorphos, the mission has achieved its planetary defense requirements, although work to further understand DART's kinetic impact test and the Didymos system will continue. In particular, ESA's Hera mission is planned to perform extensive measurements in 2027 during its rendezvous with the Didymos–Dimorphos system, building on DART to advance our knowledge and continue the ongoing international collaboration for planetary defense.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

DOE OSTI2025

Freeze-induced crystallization: An overlooked pathway for mineral genesis in natural waters

Natural ice plays salient roles in making the Earth habitable and sustainable. Previously overlooked, its role in chemical processes is now of emerging interest, particularly due to the freeze concentration effect, which can substantially promote chemical reactions during ice formation. We demonstrate here that ice formation can serve as a dynamic and unique pathway for mineral genesis. Freezing solutions containing dissolved manganese and carbonates produced rhodochrosite (Mn II CO 3 ) even under slightly undersaturated conditions. At room temperature, by contrast, this occurred when the solution saturation level was increased by ca. 30,000 times. The cryogenic rhodochrosite formed spherical aggregates of nano-polycrystallites, distinctly different from the cubic monocrystalline particles observed at room temperature. The distinct feature likely resulted from the combined effects of the intensified supersaturation induced by the freeze concentration effect and the low temperatures within liquid-like layers, conditions that make liquid-like layers an exceptional environment for mineral genesis, unlike typical natural water systems. The cryogenic rhodochrosite formation was successfully demonstrated using in situ, real-time X-ray absorption spectroscopy (XAS), enabling direct observation of freeze-induced solid formation. In conclusion, our findings reveal that freeze-induced crystallization may be an active mineralization pathway, potentially influencing elemental cycles within the cryosphere and contributing to minerals with distinguishing properties and reactivities in the environment.

Earth, Atmospheric, and Planetary Sciences↗