Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geosciences”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of Selected Publicly Available Geothermal Exploration Data Gaps

As part of a United States Department of Energy (DOE) supported retrospective analysis of DOE's Play Fairway Analysis (PFA) projects, the National Renewable Energy Laboratory (NREL) compiled and analyzed publicly available geothermal exploration datasets to identify and highlight data gaps in areas prospective for hosting geothermal resources. The analysis was intended to understand the existing geographic coverage of selected datasets commonly utilized both by the PFA projects and geothermal developers during resource assessments including geologic mapping, temperature gradient drilling, and aeromagnetic, gravimetric, and lidar surveys. Results indicate that broad areas of the western United States estimated to have geothermal potential lack sufficient geologic and geophysical coverage necessary for even regional resource exploration. The study directly informed the recent Geoscience Data Acquisition for Western Nevada, or GeoDAWN - which united DOE's Geothermal Technologies Office (GTO) with the U.S. Geological Survey (USGS) of the U.S. Department of the Interior to assist U.S. needs for energy and critical minerals. The study also has the potential to inform public investment in further data acquisition for characterization of the Earth both for geothermal and other natural resource assessments.

data↗
DOE OSTI2020

Grimsel Test Site - A Successful International Underground Research Laboratory for Many Decades - 20429

For more than 35 years, Nagra and its partners from around the world have been conducting underground research projects at the Grimsel Test Site (GTS, www.grimsel.com) to contribute to the development and confirmation of safe geological disposal concepts and for the characterization of suitable host rock formations. Over the years, the results of this internationally recognized research program have been, and continue to be, incorporated directly into exploration programs, modelling, safety, and engineering feasibility studies on options for deep geological repositories. Each project of the GTS program involves field-testing, laboratory studies, design and modelling tasks, and integrates all scientific and technical aspects. Each project phase is planned with a duration of three to five years to facilitate practical and administrative aspects and allow flexibility for updating the overall project plans with the latest findings. Scientific and engineering interaction among the different projects is ensured via an international steering committee meeting. Hosting an IAEA level C radiation- controlled zone, which allows use of radionuclides, including actinides such as thorium, uranium, neptunium, plutonium and americium, in in-situ experiments is one of the reasons why GTS also developed as a center of excellence for work with radioactive tracers under realistic in-situ boundary conditions. Last year, a new five-year program (2019 to 2023) started which includes projects with a planning horizon of decades. The new five-year program includes a new phase of in-situ experiments using radionuclides such as migration experiments in the Colloid Formation and Migration project (CFM), the Long-Term Diffusion experiment (LTD) and the newly established C-14 and I-129 Migration in cement project (CIM). The 'High Temperature effects on Bentonite' (HotBENT) project is starting in the current phase and is studying the effects of elevated temperatures (>175 deg. C) on bentonite materials. As a generic underground research laboratory (URL) it is expected that the GTS will provide in the coming years a platform for international collaboration, knowledge development and knowledge transfer for the next generation of scientists and engineers in the area of radioactive waste disposal and geosciences. A key role regarding knowledge transfer and training is provided by the well-established Grimsel Training Center (GTC), which (beside many URL related issues) also covers many general aspects of radioactive waste management. In this paper we provide an overview of the current program at the GTS, focusing on the experiments that study the migration of radionuclides through engineered barrier materials and the geosphere. (authors)

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

Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects

This is the conference paper accompanying an oral presentation “Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects” at the 17th International Conference on Greenhouse Gas Control Technologies GHGT-17 held in Calgary, Canada, October 20-24 , 2024. Carbon capture and storage (CCS) technology is critical for mitigating climate change but requires effective subsurface reservoir management to ensure safe containment of injected CO2. Accurate predictions of reservoir pressure and saturation are essential for assessing long-term CCS performance. Traditional numerical simulations, while effective, are computationally intensive, time-consuming, and constrained by data discretization. Previous work has shown the effectiveness of MeshGraphNets (MGN), a graph-based machine learning framework, as an innovative alternative for predicting reservoir behavior. MGN leverages graph neural networks (GNNs) and mesh representations to model complex geological formations, offering superior adaptability across different discretizations and reservoir configurations. Classic MGN implementations utilize an autoregressive technique to predict future behavior based on current predictions, but this technique is hampered by error accumulation over time. To enhance the model accuracy in time-series predictions, this study implemented a multi-step rollout strategy that integrates autoregressive predictions during training to stabilize prediction of saturation over time. Using the Illinois Basin – Decatur Project (IBDP) dataset, comprising 100 simulations of CO2 injection, pressure, and saturation changes, the framework demonstrated its ability to learn spatial dependencies and temporal dynamics. With inputs including permeabilities, porosities, and injection rates, MGN accurately predicted CO2 plume evolution over time, even with limited training data. Moreover, the addition of a multi-step rollout procedure during training improved the ability of MGN to predict stably over time by ~15%. This research positions MGN, enhanced with multi-step rollout capabilities, as a robust and efficient tool for CCS applications. It advances the field by enabling precise, computationally efficient predictions of reservoir behavior, providing a foundation for the broader adoption of machine learning frameworks in CCS and other geoscience domains.

Holcomb, Paul↗
DOE OSTI2025

Hydromechanical Modeling of Fault Rupture in Geologic CO 2 Sequestration: A Comparison of Two Failure Criteria

Abstract The large‐scale implementation of geological carbon sequestration has raised concerns about potential fault activation and induced seismicity, which could compromise storage integrity and pose seismic risks. We theoretically compared two failure criteria, Mohr‐Coulomb (MC) and Modified Cam‐Clay criteria (MCC), to assess fault rupture during CO 2 storage. Both criteria characterize fault behavior in a specific stress regime but differ in reducing the complexity of fault rupture to a few key mechanisms. Using a coupled hydromechanical model, we demonstrate that the choice of a failure criterion and the physics of fault weakening associated with these criteria strongly condition fault response to a given fluid injection. MC mainly relates rupture to friction, while MCC relates rupture to fault poro‐plasticity. Our findings highlight that the selection of a failure criterion, being inherently subjective, can significantly alter the predicted fault behavior during CO 2 storage, thereby impacting the reliability of geomechanical risk assessments.

Cao, M. [Energy Geosciences Division Lawrence Berk↗
DOE OSTI2024

Interfacial Engineering of Metal Chalcogenides‐based Heterostructures for Advanced Sodium‐Ion Batteries

Sodium-ion batteries (SIBs) have become one of the most promising candidates for large-scale energy storage applications. Metal chalcogenides anode materials based on alloying or conversion reactions are widely studied because of their high theoretical capacities and rich redox reactions. However, their intrinsic limitations such as high voltage hysteresis and large volume expansion hinder their further applications. The construction of heterostructures has become an attractive strategy to alleviate the above issues. Here, the formation of built in electric fields (BIEFs) at the heterointerfaces will accelerate the migration of Na + and electrons. Moreover, heterostructures can also enhance the structural stability, generate more active sites and provide additional capacity. It is worth noting that heterointerfacial properties play a significant role in promoting the overall electrochemical performance of the heterostructures. However, a systematic understanding of their interfacial engineering is currently lacking. This article reviews the research progress of metal chalcogenides-based heterostructure anode materials in the near term. First, the definition, classification and the roles of heterostructures are introduced. Second, the detailed research progress of the metal chalcogenide-based heterostructures anodes in SIBs is discussed. Finally, the future prospects and potential research directions of the heterostructures for batteries are discussed.

anode↗
DOE OSTI2025

Effect of CO 2 Concentration on the Microbial Activity of Orenia metallireducens (Strain Z6) in Surface Inert Materials

Carbon dioxide (CO 2 ) sequestration has garnered widespread attention as a key strategy for mitigating CO 2 emissions and combating the greenhouse effect. However, the mechanisms underlying the interactions between CO 2 , widespread siliceous minerals and biological processes remain unclear. The present study explored the potential impacts of different CO 2 concentrations on microbial activity, environmental conditions and their feedback on the fate of CO 2 . A total of 20 experimental conditions was created, with the variables including different natural and synthetic siliceous minerals (e.g., quartz sand and a type of commercial glass beads), the presence or absence of the iron-reducing microorganism Orenia metallireducens (strain Z6) and varying CO 2 concentrations (0%, 20%, 50%, 100%) in the presence of ferrihydrite and pyruvate. Geochemical, microbial and mineralogical analyses revealed that elevated CO 2 concentrations significantly inhibited microbial Fe(III) reduction and pyruvate metabolism. Interestingly, compared to cultures without mineral amendments or those with glass beads alone, the addition of quartz sand enabled strain Z6 to better withstand the environmental stress caused by elevated CO 2 , promoting pyruvate fermentation and iron reduction. In addition to an increased pH, the formation of siderite, hematite and vivianite was also observed in the bioactive systems. Although both glass beads and quartz sand were primarily composed of silica, differences in the mineral structure, elemental composition and acid neutralization capacity rendered quartz sand more chemically active and unexpectedly led to greater CO 2 sequestration.

CO2 stress↗
DOE OSTI2025

Hydrometallurgical Recycling of Black Mass of Spent Lithium-Ion Batteries Using Methanesulfonic Acid: Leaching, Kinetic Studies, and Potential for Total Recovery of Valuable Components

Methanesulfonic acid (MSA) exhibits several advantageous properties rendering it a promising candidate for circular hydrometallurgical processes. These properties include a high acidity (pKa = - 1.9) comparable to that of classical mineral acids as well as biodegradability, high stability, and high solubility of metal-MSA complexes in aqueous solutions. In this study, MSA was employed as a lixiviant for the leaching of metals (lithium, nickel, cobalt, and manganese) from the black mass of spent lithium-ion batteries (LIBs). The effect of various parameters, including MSA concentration, H 2 O 2 concentration, temperature, and pulp density, was systematically investigated. Under the optimized conditions (1.5 M MSA, 0.2 M H 2 O 2 , 60 °C, and 50 g/L pulp density), quantitative leaching of lithium was achieved within 30 min, while for nickel and cobalt it was after 2 h, and 4 h for manganese leaching. The leaching kinetics of Li, Ni, Co, and Mn were studies using the shrinking particle models (SPM) and the Avrami model. The results indicated that the Avrami model provided the best fit to the kinetic data, with apparent activation energies of 46.81 kJ/mol for Li, 58.61 kJ/mol for Ni, 59.69 kJ/mol for Co, and 58.86 kJ/mol for Mn, within the temperature range of 25-70 °C (except for Li, which was analyzed in the range of 25-60 °C), consistent with chemical reaction control. Subsequently, residual contaminants in the leaching residue were eliminated through pyrolysis. The quantitative leaching of metals in MSA solution (a green lixiviant), combined with the pyrolytic treatment of leaching residues, represents a circular strategy for the total recovery of valuable components from the black mass of spent LIBs.

25 ENERGY STORAGE↗
DOE OSTI2025

Construction of MoS 2 /MoO 3 Heterostructure with Ultrafast-Charged and Superior Low-Temperature Sodium Storage Properties

The electrochemical reaction kinetics of sodium-ion batteries (SIBs) become sluggish at low temperatures, resulting in significant reductions in energy density and power density. Rational design of anode materials with excellent low-temperature performance is of great significance for promoting the application of SIBs under extreme conditions. Here, in this study, the spontaneous hydrolysis and oxidation reactions of MoS 2 in aqueous solution are used to successfully construct the MoS 2 /MoO 3 heterostructure. The formation of built-in electric fields at the MoS 2 /MoO 3 heterointerfaces improves the electrochemical reaction kinetics, thereby enhancing the rate performance. In addition, the dual-phase material can effectively buffer the volume strain during the cycle process, thereby improving the cycle stability. Thus, the MoS 2 /MoO 3 displays ultrafast charging properties at room temperature (up to 244.6 mAh g –1 at 40 A g –1 , discharge/charge in 22 s). Even at −40 °C, it also exhibits a high capacity of 303.7 mAh g –1 and superior cycling performance (capacity retention rate up to 92.9% after 900 cycles at 2 A g –1 ).

MoS2↗
DOE OSTI2025

An international study on THM modelling of the full-scale heater experiment at Mont Terri laboratory

We present results from an international model comparison study of the Full-Scale Emplacement (FE) experiment in Opalinus Clay at the Mont Terri Laboratory, Switzerland. Based on a provided parameter set the teams decided which parameters they adopted for their models, whether they considered the excavation and the ventilation phase in addition to the heating phase and if they included technical features like the shotcrete or the EDZ. The teams were able to reproduce the measured parameters temperature, relative humidity and pore pressure. The modelled results for temperature agree very closely between the teams especially in the sensors in Opalinus Clay. All teams were able to reproduce the redistribution of water in the bentonite backfill due to heating. The evolution of the relative humidity showed similar trends with differences in the intensity of the dry out effect. To model the pore pressure evolution is more complex because it comprises the full interaction of the coupled THM processes. The spread between the pore pressure modelled by the teams was larger, with some teams overestimating the pressure increase due to heating and some teams overestimating the extent of drainage. The agreement of modelled results with measurements improves with larger distance to the heater. We conclude that the EDZ and the shotcrete potentially influence the behaviour of the rock causing higher differences closer to the heater. Further research is needed to better implement those influences into the models. Based on the calibrated models, the future evolution of temperature, relative humidity and pore pressure was predicted over the next 10 years following a change of the heat power applied in 2023 and 2024. Again, the predicted temperatures agree very closely between the teams. Most teams do not expect an increase in relative humidity during the next 10 years after the initial dry-out.

58 GEOSCIENCES↗
DOE OSTI2025

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

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

58 GEOSCIENCES↗
DOE OSTI2025

GNSS Installation at the ARM Southern Great Plains (SGP) Atmospheric Observatory Field Campaign Report

GRUAN, the GCOS Reference Upper Air Network, is an international reference observing network of sites measuring essential climate variables above Earth's surface, designed to fill an important gap in the current global observing system. GRUAN’s objectives are to (1) provide long-term, high-quality climate records, (2) constrain and calibrate data from more spatially comprehensive global observing systems (including satellites and current radiosonde networks), and (3) fully characterize the properties of the atmospheric column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Assessing Seismic Risk for CO2 Geologic Storage: Comparative Analysis of the Delaware Basin and Basin and Range Province Projects

ABSTRACT: Effective management of induced seismicity is critical for safe and sustainable CO2 storage. This study evaluates fault slippage risks in the Delaware Basin (Texas) and Basin and Range Province (Utah), integrating geological, operational, and geomechanical parameters to assess fault stability and seismic hazard mitigation. In the Delaware Basin, two sites were analyzed under an injection rate of 20,000 bbl/day over 25 years. One site showed low fault slip risk, while the other exhibited higher reactivation potential due to proximity to critically stressed faults. Sensitivity analysis revealed that increased pore pressure significantly heightened slip potential, highlighting the necessity of precise pressure control and real-time monitoring. In the Basin and Range Province, fault stability was evaluated at Neck of the Desert, Escalante Desert, Parowan, and Beaver sites under injection rates of 8,750 bbl/day per site over 30 years. Minimal fault slip risk was observed at Neck of the Desert and Escalante Desert sites, whereas Parowan and Beaver sites exhibited elevated slip potential due to semi-critically stressed faults sensitive to modest pore pressure increases. The findings demonstrate that fault slippage analysis, combined with sensitivity analysis of pore pressure and friction coefficients, is essential for understanding seismic risks. Continuous monitoring, adaptive injection management, and rigorous geomechanical analysis are key strategies for minimizing induced seismicity in CO2 sequestration projects.

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 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

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗