Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of Scalable Reactive Transport Framework for PFAS

Research on PFAS chemicals is extensive and covers all aspects, including analytical quantification, determination of properties, toxicology, ex situ treatment, and in situ remediation. PFAS chemicals are very stable and thus persistent/recalcitrant in the environment. Although there are many unknowns about PFAS chemicals, degradation pathways, reaction rates, etc., many different sorption, oxidation, reduction, biological, and innovative treatment approaches are being developed. Sorption with activated carbon is currently the only fully available in situ treatment technology for PFAS-impacted groundwater. Given the wide array of PFAS chemicals and transformation products, remediation may need multi-step treatment trains to fully address the PFAS contamination. The work here provides kinetic reaction modules that represent an initial cut at functionality representing PFAS migration and reaction in groundwater aquifer flow and transport models. One reaction kinetics module provides a method to model kinetically limited adsorption using a mass transfer model. The second reaction module represents biological transformation of 8:2 FTOH and daughter species, illustrating how a complex reaction pathway network can be represented. Both reaction modules allow for spatially variable parameter values so that a variety of remediation approaches (e.g., a PRB or volumetric treatment or variations in geochemical conditions) can be simulated. The intent with these PFAS reaction modules is to provide tools for practitioners to aid in the selection, design, and assessment of potential in situ PFAS remediation strategies. It is anticipated that, as PFAS remediation technologies and scientific understanding advances, these modules would be refined or replaced to match new knowledge.

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

Can we safely go to 200 °C? An integrated approach to assessing impacts to the engineered barrier system in a high-temperature repository

Abstract. This presentation gives on overview of the complex thermo-hydro-mechanical-chemical (THMC) processes occurring during the disposal of heat-producing high-level radioactive waste in geologic repositories. A specific focus is on the role of compacted bentonite, which is commonly used as an engineered backfill material for emplacement tunnels because of its low permeability, high swelling pressure, and radionuclide retention capacity. Laboratory and field tests integrated with THMC modeling have provided an effective way to deepen our understanding of temperature-related perturbations in the engineered barrier system; however, most of this work has been conducted for maximum temperatures around 100 ∘C. In contrast, some international disposal programs have recently started investigations to understand whether local temperatures in the bentonite of up to 200 ∘C could be tolerated with no significant changes to safety relevant properties. Raising the maximum temperature is attractive for economical and safety reasons but faces the challenge of exposing the bentonite to significant temperature increases. Strong thermal gradients may induce complex moisture transport processes while geochemical processes, such as cementation and perhaps also illitization effects may occur, all of which could strongly affect the bentonite and near-field rock properties. Here, we present initial investigations of repository behavior exposed to strongly elevated temperatures. We will start discussing our current knowledge base for temperature effects in repositories exposed to a maximum temperature of 100 ∘C, based on data and related modeling analysis from a large heater experiment conducted for over 18 years in the Grimsel Test Site in Switzerland. We then show results from coupled THMC simulations of a nuclear waste repository in a clay formation exposed to a maximum temperature of 200 ∘C. We also explore preliminary data from a bench-scale laboratory mock-up experiment, which was designed to represent the strong THMC gradients occurring in a “hot” repository, and we finally touch on a full-scale field heater test to be conducted soon in the Grimsel Test Site underground research laboratory in Switzerland (referred to as HotBENT).

Birkholzer, Jens T. (ORCID:0000000279891912)↗
DOE OSTI2023

Noble gas insights into early impact delivery and volcanic outgassing to Earth's atmosphere: A limited role for the continental crust

Earth's atmosphere, crust and mantle have evolved together through continuous geochemical exchange throughout Earth's history. Constraints on the transport of volatile elements and compounds between these reservoirs are crucial for understanding how Earth could have stayed habitable for extended periods of time. In this work, we present a new forward model of He, Ne and Ar in the mantle, crust and atmosphere. We explore concentrations of noble gases at the end of accretion, bulk silicate Earth K/U ratios, crustal growth scenarios, and upper mantle processing rates throughout Earth's history. We search for parameter combinations that simultaneously satisfy observational constraints on present-day mantle 4 He/ 3 He and 40 Ar/ 36 Ar (sensitive to mantle outgassing and continental crust growth, which depletes the mantle of U, Th and K), atmospheric 20 Ne/ 22 Ne (which tracks the mix of outgassed vs. delivered Ne in the atmosphere), and atmospheric 40 Ar/ 36 Ar in the past and today (sensitive to volatile delivery, mantle outgassing and crustal growth). Leveraging this intertwined set of noble gas abundances and isotopic compositions yields new constraints on initial noble gas abundances in the mantle and on the proportions of atmospheric volatiles originating from delivery by impact degassing, mantle outgassing, and degassing of the continental crust. We find that atmospheric Ar isotopic evolution is primarily sensitive to the mantle processing rate history; the atmospheric Ar isotopic record should therefore not be used to reconstruct continental crust growth, but instead provides valuable insights into mantle processing rates. Our model predicts a measurably low 20 Ne/ 22 Ne ratio of ~9.7 in Archean atmospheric samples. Most of atmospheric primordial 36 Ar was directly delivered by chondritic bodies and not transferred to the atmosphere during an intense early episode of mantle outgassing. Nitrogen delivered by impact degassing could account for the present-day atmospheric nitrogen inventory.

58 GEOSCIENCES↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2025

Building confidence in models for complex barrier systems for radionuclides

The modeling and simulation of the Cement-clay Interaction-Diffusion field (CI-D) experiment at the Mont Terri site in Switzerland presented here demonstrates that it is possible to capture the multiscale physical and chemical features of natural and engineered barrier systems for radionuclides. The simulations are successfully carried out with the newly developed CrunchODiTi high-performance computing software that accounts for multiple continua, including a continuum representing the electrical double layer (EDL) developed along negatively charged clay particles in clay rock. The simulation also accounts for both the complex three-dimensional (3D) geometry, expected as the norm in a geological waste repository, and the anisotropy of the geological formation. In addition, the high resolution of the model makes it possible to include "skin effects" developed at the interface between highly reactive materials, in this case between the high pH cement and the circumneutral but electrostatic Opalinus Clay. The successful history matching with the field experiment demonstrates that the distinct geochemical and physical properties of the cement and the Opalinus Clay in the CI-D experiment can be accounted for. Such analyses are essential for developing a defensible safety case for the underground storage of radioactive waste.

Sarsenbayev, Dauren↗
DOE OSTI2024

Potential applications of microbial genomics in nuclear non-proliferation

As nuclear technology evolves in response to increased demand for diversification and decarbonization of the energy sector, new and innovative approaches are needed to effectively identify and deter the proliferation of nuclear arms, while ensuring safe development of global nuclear energy resources. Preventing the use of nuclear material and technology for unsanctioned development of nuclear weapons has been a long-standing challenge for the International Atomic Energy Agency and signatories of the Treaty on the Non-Proliferation of Nuclear Weapons. Environmental swipe sampling has proven to be an effective technique for characterizing clandestine proliferation activities within and around known locations of nuclear facilities and sites. However, limited tools and techniques exist for detecting nuclear proliferation in unknown locations beyond the boundaries of declared nuclear fuel cycle facilities, representing a critical gap in non-proliferation safeguards. Microbiomes, defined as “characteristic communities of microorganisms” found in specific habitats with distinct physical and chemical properties, can provide valuable information about the conditions and activities occurring in the surrounding environment. Microorganisms are known to inhabit radionuclide-contaminated sites, spent nuclear fuel storage pools, and cooling systems of water-cooled nuclear reactors, where they can cause radionuclide migration and corrosion of critical structures. Microbial transformation of radionuclides is a well-established process that has been documented in numerous field and laboratory studies. These studies helped to identify key bacterial taxa and microbially-mediated processes that directly and indirectly control the transformation, mobility, and fate of radionuclides in the environment. Expanding on this work, other studies have used microbial genomics integrated with machine learning models to successfully monitor and predict the occurrence of heavy metals, radionuclides, and other process wastes in the environment, indicating the potential role of nuclear activities in shaping microbial community structure and function. Results of this previous body of work suggest fundamental geochemical-microbial interactions occurring at nuclear fuel cycle facilities could give rise to microbiomes that are characteristic of nuclear activities. These microbiomes could provide valuable information for monitoring nuclear fuel cycle facilities, planning environmental sampling campaigns, and developing biosensor technology for the detection of undisclosed fuel cycle activities and proliferation concerns.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Bed form-induced hyporheic exchange and geochemical hotspots

Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hotspots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. The results in this research suggest that out-of-phase bed forms generate more complex hyporheic flow patterns which reduce the efficiency of solute transformations. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. Here this study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes including the rate and spatio-temporal distribution of reaction hotspots.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Data and scripts associated with a manuscript investigating dissolved organic matter and microbial community linkages across seven globally distributed rivers

This data package is associated with the publication “Meta-metabolome ecology reveals that geochemistry and microbial functional potential are linked to organic matter development across seven rivers” submitted to Science of the Total Environment. This data package includes the data necessary to replicate the analyses presented within the manuscript to investigate dissolved organic matter (DOM) development across broad spatial distances and within divergent biomes. Specifically, we included the Fourier transform ion cyclotron mass spectrometry (FTICR-MS) data, geochemistry data, annotated metagenomic data, and results from ecological null modeling analyses in this data package. Additionally, we included the scripts necessary to generate the figures from the manuscript. Complete metagenomic data associated with this data package can be found at the National Center for Biotechnology (NCBI) under Bioproject PRJNA946291. This dataset consists of (1) four folders; (2) a file-level metadata (flmd) file; (3) a data dictionary (dd) file; (4) a factor sheet describing samples; and (5) a readme. The FTICR Data folder contains (1) the processed Fourier transform ion cyclotron mass spectrometry (FTICR-MS) data; (2) a transformation-weighted characteristics dendrogram generated from the FTICR-MS data; and (3) the script used to generate all FTICR-MS related figures. The Geochemical Data folder contains (1) the single geochemistry data file and (2) the R script responsible for generating associated figures. The Metagenomic Data folder contains (1) annotation information across different levels; (2) carbohydrate active enzyme (CAZyme) information from the dbCAN database (Yin et al., 2012); (3) phylogenetic tree data (FASTAs, alignments, and tree file); and (4) the scripts necessary to analyze all of these data and generate figures. The Null Modeling Data folder contains (1) data generated during null modeling for each river and all rivers combined and (2) the R scripts necessary to process the data. All files are .csv, .pdf, .tsv, .tre, .faa, .afa, .tree, or .R.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Roughrider Carbon Storage Hub (Final Report)

The Roughrider Carbon Storage Hub was a 2-year project (October 2023 – September 2025) conducted by the Energy & Environmental Research Center (EERC) focused on advancing the feasibility of a commercial-scale carbon dioxide (CO 2 ) geologic storage hub in McKenzie County, North Dakota. The project’s objective was to investigate the potential that stacked storage complexes (multiple deep saline formations) can safely and economically store at least 50 million tonnes of CO 2 within 30 years. The captured CO 2 would be sourced from industrial emitters including project partner ONEOK, Inc.’s gas-processing plants and a planned gas-to-liquids facility. Drilling of the Roughrider 1 stratigraphic test well (14,979-ft total depth) was completed in November 2024. The wellbore intersected four candidate storage formations: Inyan Kara, Broom Creek, Mission Canyon, and Black Island–Deadwood. Operational challenges, including a stuck drill string, were resolved without long-term impact. A comprehensive logging and coring program was conducted, followed by successful well abandonment and site reclamation. Over 660 ft of 4-in. whole core was retrieved. Core plug samples were processed and analyzed for petrophysical and geochemical properties. Results confirmed promising porosity and permeability in the Inyan Kara and Broom Creek Formations and removal of the Mission Canyon and Black Island–Deadwood horizons from further investigation. Data derived from the logging and coring program were used to improve initial geologic models built from legacy data. CO 2 injection simulations showed that the Inyan Kara alone can feasibly store the target mass of CO 2 . Because of subtle differences in geologic structure and porosity trends between the formations, a stacked storage scenario using the Broom Creek and Inyan Kara Formations resulted in a larger overall plume area than using the Inyan Kara alone. Preliminary CO 2 pipeline routes from the industrial sources were mapped utilizing existing rights of way and evaluated for capacity and cost using U.S. Department of Energy Office of Fossil Energy and Carbon Management/National Energy Technology Laboratory models and U.S. Environmental Protection Agency emissions data. Integrating capture, transport, and storage cost estimates with policy incentives (e.g., 45Q credits) provided a total cost-per-ton analysis. Results indicate that the small scale of the volumes to be transported over the cumulative large distances does not support the project’s financial viability. However, the groundwork laid during this project from geological, regulatory, and social perspectives positions the Roughrider hub site as a promising candidate for commercial carbon storage in North Dakota, especially if the economy of scale is introduced for CO 2 transportation to the hub site.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Impact of Mineral Reactive Surface Area on Forecasting Geological Carbon Sequestration in a CO2-EOR Field

Mineral reactive surface area (RSA) is one of the key factors that control mineral reactions, as it describes how much mineral is accessible and can participate in reactions. This work aims to evaluate the impact of mineral RSA on numerical simulations for CO2 storage at depleted oil fields. The Farnsworth Unit (FWU) in northern Texas was chosen as a case study. A simplified model was used to screen representative cases from 87 RSA combinations to reduce the computational cost. Three selected cases with low, mid, and high RSA values were used for the FWU model. Results suggest that the impact of RSA values on CO2 mineral trapping is more complex than it is on individual reactions. While the low RSA case predicted negligible porosity change and an insignificant amount of CO2 mineral trapping for the FWU model, the mid and high RSA cases forecasted up to 1.19% and 5.04% of porosity reduction due to mineral reactions, and 2.46% and 9.44% of total CO2 trapped in minerals by the end of the 600-year simulation, respectively. The presence of hydrocarbons affects geochemical reactions and can lead to net CO2 mineral trapping, whereas mineral dissolution is forecasted when hydrocarbons are removed from the system.

58 GEOSCIENCES↗
DOE OSTI2022

NEWTS USGS Produced Waters Database

Data from the USGS Produced Waters Database on produced water compositions from across the US. Data has been filtered for completeness of data attributes. Data has been transposed and converted for ease of input into aqueous chemistry modeling software including OLI Systems and Geochemist's Workbench. Original data from: Blondes, M. S., Gans, K. D., Engle, M. A., Kharaka, Y. K., Reidy, M. E., Saraswathula, V., Thordsen, J. J., Rowan, E. L., & Morrissey, E. A. (2019). U.S. Geological Survey National Produced Waters Geochemical Database v2.3 [Data set]. U.S. Geological Survey. https://doi.org/10.5066/F7J964W8

Geochemists Workbench↗
DOE OSTI2022

Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation

Shallow subsurface microbial and geochemical processes in watersheds are dynamic and responsive to seasonal and long-term environmental change. At high elevation or latitude, such changes may occur during winter months when normal sampling is impossible, impractical, unsafe, or may actually alter the processes being studied. Yet accurate modeling of biogeochemical parameters in these environments requires sampling that can capture events and persistent cold climate trends. Gaps in current models exist where data have not been collected during extended snow- or ice-covered periods. The goals of our exploratory study were to examine the biogeochemical processes occurring in a well-studied high altitude watershed during winter when the system is largely inaccessible. We tested the hypothesis that during snow cover in the East River (ER) watershed, episodic excursions of microbial community structure and biogeochemical processes and concentrations fluctuated from values extrapolated from pre- and post-snow time periods. Through our research we identified previously unresolved subsurface-surface transport phenomena and a heretofore under-appreciated methane cycle occurring during wintertime months, thus providing valuable insight into watershed processes. Our research showed that autonomous sampling techniques, first used in marine systems, is extremely versatile and can be adapted for use in harsh continental settings to examine temporal changes in freshwater chemistry and microbial community assemblage. Further, we found that by sampling during periods when systems are hard to access, will lend insight and better accountability for global budgets for methane, an important greenhouse gas. By strengthening our knowledge of wintertime biogeochemical relationships in high altitude watersheds our research leads to more complete descriptions of cryptic processes and ultimately will inform and improve the reactive transport models that simulate processes in these systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗
DOE OSTI2023

Endothermic physiology of extinct megatooth sharks

The evolution of the extinct megatooth shark, Otodus megalodon , and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark’s gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Impact of reservoir parameters and wellbore permeability uncertainties on CO 2 and brine leakage potential at the Shenhua CO 2 Storage Site, China

Here, this paper describes a system-level risk assessment for the Shenhua CO 2 storage site, China, using the National Risk Assessment Partnership Integrated Assessment Model for Carbon Storage (NRAP-IAM-CS). Here, we begin by determining the optimal number of Monte Carlo (MC) simulations to achieve CO 2 and brine leakage result convergence. Then, we calculate mass CO 2 and brine leakage to the atmosphere and a hypothetical shallow aquifer. Finally, we assess the geochemical impacts in the event of leakage as if there were a shallow freshwater aquifer at the Shenhua site. Simulation results show that leakage results tend to stabilize after 300 Monte Carlo simulations. When the three wells on site are assigned a permeability of 10 -11 m 2 (representing significantly leaking wells), moderate CO 2 and brine leakages occur, and the percentage of CO 2 leakage exceeds the threshold value we set based on the Intergovernmental Panel on Climate Change (IPCC). This is, however, unlikely to be the case for the Shenhua site. For all the other scenarios, the CO 2 leakage is trivial although there is still the possibility of CO 2 leakage into the groundwater aquifer exceeding the 1% threshold over 1000 years, assuming constant legacy wellbore permeability. For the significantly leaking scenario, there is a 10% probability to have a moderate (2.7 × 10 7 m 3 ) leakage-affected volume in the shallow aquifer with the pH below 6.5, and a large (1.4 × 10 8 m 3 ) volume with the total dissolved solids (TDS) above 500 ppm, hence pH and TDS may be considered for site monitoring plans. Based on the simulation results, there is a very low probability of significant CO 2 /brine leakage through the existing wells at the Shenhua CO 2 storage site.

58 GEOSCIENCES↗
DOE OSTI2025

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Data‐Driven Insights into Rare Earth Mineralization: Machine Learning Applications Using Functional Material Synthesis Data

Understanding rare‐earth element (REE) mineralization mechanisms is essential for developing efficient separation strategies. Although the geochemical pathways that generate REE deposits are qualitatively known, quantitative links between specific conditions and mineralization outcomes remain limited. Herein, the repurpose laboratory REE hydrothermal synthesis data—originally collected for functional‐materials fabrication—as a surrogate for studying mineralization with data‐driven methods. The compiled 1,200+ hydrothermal reaction records and trained three machine‐learning models—K‐nearest neighbors (KNN), random forest (RF), and extreme gradient boosting (XGB)—to predict product elements and phases from precursors, additives, reaction conditions, and engineered features. Validation shows XGB achieves the highest accuracy. Feature importance indicates thermodynamic properties of cations and anions dominate model decisions. Correlations reveal positive relationships among precursor concentration, reaction time, pH, and temperature, consistent with classical crystallization behavior. XGB‐based regressors are built to predict crystallization temperature and pH from precursor/product attributes. Performance is strongest when similar training examples exist, while accuracy declines for underrepresented reactions, notably REE carbonates and heavy‐REE systems. Overall, the study shows that functional‐materials datasets can illuminate REE mineralization and provide priors for exploration and processing. Expanding datasets with less‐studied chemistries and conditions will improve generality and support deposit discovery and more efficient REE recovery.

feature importance analysis↗