Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Comparison of mercury (Hg) bioaccumulation with mono- and mixed Lemna minor and Spirodela polyrhiza cultures

Mercury (Hg) is a prevalent and harmful contaminant that persists in the environment. For phytoremediation, it is important to discover which plants can bioaccumulate meaningful amounts of Hg while also tolerating its toxicity. Additionally, increasing biodiversity could create a more resilient and self-sustaining system for remediation. This study explores whether mixed populations of Lemna minor and Spirodela polyrhiza can better bioaccumulate and tolerate Hg than monocultures. Mono- and mixed cultures of L. minor and S. polyrhiza were grown in mesocosms of 0.5 μg/L or 100 μg/L Hg ( HgCl2) spiked water for 96 h. Change in weight of duckweed was used to assess Hg tolerance. Diffusive gradients in thin-films (DGTs) were used as surrogate monitoring devices for bioavailable levels of Hg. For biomass growth, the mixed culture of the L. minor was greater than the monoculture at the high dose. The L. minor accumulated more Hg in the mixed culture at the low dose while the S. polyrhiza was higher in the mixed at the high dose. Hg speciation in water was modeled using Windermere Humic Aqueous Model 7 (WHAM7) to compare the bioavailable species indicated by the DGTs. Potentially due to the controlled conditions, the WHAM7 output of bioavailable Hg was almost 1:1 to that estimated by the DGTs, indicating good predictive capability of geochemical modeling and passive sampler DGT on metal bioavailability. Altogether, the mixed cultures statistically performed as well as or better than the monocultures when tolerating and bioaccumulating Hg. However, there needs to be further work to see if the significant differences translate into practical differences worth the extra resources to maintain multiple species.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
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↗
DOE OSTI2024

Jackson, L., Johnson, M.B., Latrach, A., Grimes, D., Martinez, C., and Mclaughlin, J.F., 2024, Multidisciplinary geotechnical data collection, curation, and analysis for conformity with the regulatory framework for geologic carbon storage in Wyoming, USA: Geological Society of America Abstracts with Programs. Vol. 56, No. 5, 2024, doi: 10.1130/abs/2024AM-405024

Title: Multidisciplinary Geotechnical Data Collection, Curation, and Analysis for Conformity with the Regulatory Framework for Geologic Carbon Storage in Wyoming, USA. Text: Construction and operation of wells for geologic sequestration of carbon dioxide necessitate that they are permitted under the Environmental Protection Agency’s Underground Injection Control Class VI requirements. Class VI wells conform to stringent requirements to ensure long-term safety and integrity of the storage site and the protection of Underground Sources of Drinking Water. Entities pursuing Class VI permitting must provide comprehensive geologic site characterization, including regional geologic structure and stratigraphy, aquifer information, reservoir and confining unit geomechanical properties, geochemical analyses, assessment of trapping capacity and mechanisms, and a variety of other of multidisciplinary geotechnical data. The Wyoming Class VI Site Characterization Database Project is focused on developing a geologic site characterization database of geotechnical information, which has been compiled and verified from established, public databases/entities and scientific literature to expedite Class VI permitting in Sweetwater County within the Greater Green River Basin of southern Wyoming. The preliminary suite of compiled data from 14,000 wells includes 8,000 wells with logs and 7,250 wells with formation tops, ~70 wells with core data (e.g., X-Ray diffraction, petrographic, and petrophysical data), ~2,500 water analyses, ~740 seismic events data, and ~520 bottom-hole temperature measurements. Future work on—and stemming from—this project will include new core analyses, calculation and interpolation of subsurface temperature gradients, mechanical earth models, geochemical simulations, storage capacity estimation, stratigraphic column generation and correlation, and construction of subsurface maps. Finally, this work will help to inspire and facilitate subsurface data compilation and curation beyond Sweetwater County, Wyoming.

42 ENGINEERING↗
DOE OSTI2023

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

The Stability of Nd Hydroxyl Complexes at Near Neutral to Alkaline pH from 25 to 75 °C: Implications for Monazite Solubility in Hydrothermal Aqueous Fluids

Speciation of rare earth elements (REE) in aqueous fluids controls their mobilization during fluid-rock interaction. Thermodynamic modeling provides important insights into the factors controlling REE mobility and fractionation in aqueous fluids, in which aqueous REE hydroxyl complexes are particularly important across a broad pH range. However, the stability of these REE hydroxyl complexes remains poorly determined due to limited experimental data above ambient temperature. In this work, UV–vis spectrophotometry experiments were conducted to determine the cumulative hydrolysis constants for neodymium (Nd) from 25 to 75 °C in near-neutral to alkaline solutions. The color indicator m -cresol purple was used to determine in situ pH. Alkaline NaOH-bearing aqueous solutions were doped with varying initial Nd concentrations (0 to ∼0.155 mmol/kg), which resulted in the release of protons (H + ) and pH decrease during hydrolysis of Nd 3+ according to Nd 3+ + n H 2 O = Nd(OH) n 3– n + H + . The average OH – ligand number coordinated to Nd 3+ was found to increase from 1.0 to 1.6 at 25 °C and from 1.2 to 2.8 at 75 °C over a pH range from 6.3 to 9.0. The measured speciation shows an increased predominance of Nd(OH) 3 0 over the Nd(OH) 2 + and Nd(OH) 2+ species with increased temperature and pH. The increasing cumulative formation constants (log*β° n , n = 1 to 3) retrieved from 25 to 75 °C differ by 0.1 to 1.9 logarithmic units in comparison to existing thermodynamic databases. These updated thermodynamic data have important implications for geochemical modeling of the speciation of Nd hydroxyl complexes and the solubility of monazite as a function of pH.

58 GEOSCIENCES↗
DOE OSTI2025

Selective recovery of native copper from basalt tailings using alkaline glycinate solution

Mafic and ultramafic rocks present an intriguing pathway for CO 2 capture through strategic enhancements to the natural silicate weathering cycle. Simultaneously, these rock types are often hosts to appreciable amounts of metals critical to U.S. energy independence, particularly in the context of alkaline mine tailings from historical metal mining. The research and scientific prerogative, then, is to identify promising mafic and ultramafic feedstocks and conduct exploratory studies to effectively recover these critical minerals while preserving the potential for mineral carbonation. The Keweenaw Peninsula of Michigan hosts the largest native Cu reservoir within basalt in the world and experienced a rich history of Cu production spanning from pre-historic times to the mid-1900s. Today, much of the Cu mining legacy remains as mine waste tailings. In this study, we examined the Cu extractability from Keweenaw Basalt tailings using a sodium glycinate solution, in comparison with acid and sodium hydroxide leaching. Experimental results showed an 85 % Cu extraction rate using sodium glycinate as the extraction solution with negligible release of other cations from the basalt. The kinetic and extraction mechanisms of Cu selective recovery using glycinate solution were discussed using time-resolved experimental data and kinetic geochemical modeling. Theoretical estimation of carbon mineralization potential of all the existing basalt waste tailings (∼500 million tons) can reach 85.5 MMT CO 2 . A total of 0.786 MMT Cu can be recovered with sodium glycinate, with a value of 7.7 billion USD. In conclusion, this novel application of alkaline glycinate for selective Cu recovery from basalt mine tailings demonstrates the viability of selective metal recovery using a non-hazardous chemical while preserving CO 2 capture potential and presents a potential pathway toward reducing energy-related emissions and providing an unconventional domestic source of critical minerals.

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

Chemical structure and complex growth modes of magnesium silicate hydrate: Nanoparticle orientation, aggregation, and fusion

The extent of utilization of magnesium silicate hydrate (MSH) in construction is limited partly because of insufficient data ascertaining the kinetics of its precipitation. In this report MSH grown homogeneously or heterogeneously in the range of magnesium to silicon solution molar concentration ratios, [Mg]/[Si] = 0.5-1.5, and temperature = 25-80 °C was analyzed for its chemical structure and morphology. Infrared spectroscopy and geochemical modeling show that increasing the [Mg]/[Si] ratio results in less silicate polymerization in MSH. Atomic force microscopy reveals the presence of nanoparticles that aggregate and fuse to form films. Oriented attachment of nanoparticles and features indicative of enhanced crystallinity were observed at higher temperatures and longer reaction times. These data provide direct evidence for the persistence of hierarchical structures (i.e., nanoparticles forming 3D microparticles and 2D film layers) at the nanoscale to mesoscale. These findings offer insights into the precise chemical synthesis of MSH and its widespread use as a binder for construction purposes.

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

Mechanisms of dissolution from gibbsite step edges elucidated by ab initio molecular dynamics with enhanced sampling

Surface reactions underpin our collective understanding of mineral dissolution. This impacts a range of predictive geochemical models (for weathering, the fate and transport of metals) as well as industrial utilization of minerals. Yet atomistic details are rarely known due to the complex mineral/fluid interfacial environment. There is significant need to understand the mechanistic details of dissolution reactions and how they depend on surface morphology and solution conditions. This work utilizes surface pit models in conjunction with changes to solution composition that mimic pH to explore surface detachment of aluminate from gibbsite, which is a primary source of Al in soils and within the industrial processing of aluminum. Ab initio molecular dynamics simulations with enhanced sampling has been used to explore the detailed process of the detachment, the results of which indicate two potential pathways that are differentiated based upon the extent of water hydration. Here, the heights of the energy barriers depend upon the local morphology which influence the number of bridges (quasi-)simultaneously broken (1 or 2) or the Al-O coordination of the neighboring aluminum atoms (5 or 6) at the armchair edge. pH effects are significant, with a nearly 50% reduction in barrier height under alkaline conditions that are relevant to geothermal fluids and Al extraction from minerals.

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

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗
DOE OSTI2024

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗
DOE OSTI2023

Role of pH and Eh in geothermal systems: Thermodynamic examples and impacts on scaling and corrosion

The global demand for sustainable and clean renewable energy sources is rapidly increasing to meet growing energy demands while curtailing the CO2 emissions from fossil fuels. Harnessing the natural heat or geothermal energy of the Earth's interior can simultaneously address both issues by offering a clean and renewable energy source. Enhanced geothermal systems require the movement and cycling of fluid through a porous rock media. Further, this process requires understanding the chemistry occurring at the fluid-rock and fluid-metal interfaces both in the deep subsurface and surface, where fluid is circulated for energy transfer or storage. In this review paper, we evaluate the current state of research on how two key variables, pH and Eh, control the fluid-rock reactions controlling mineral scaling and the steel corrosion reactions associated with geothermal systems. These variables are frequently used in geochemical modeling to predict chemical reactions in geothermal systems. Prediction of the fluid-rock behavior can be utilized to formulate new strategies for extending the lifespan and efficiency of heat extraction and utilization and maximizing the safety of geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗
DOE OSTI2025

Complex carbonate phases drive geologic CO2 mineralization

Abstract Geologic carbon sequestration in mafic and ultramafic reservoirs is a scalable strategy for carbon dioxide removal, offering permanent storage via mineralization as stable carbonates. However, there is limited information on the structure and composition of key mineralization endpoints during sequestration. Here, we unravel the atomic structure, composition, and nanoscale morphology of carbonates recovered from a field-scale demonstration of CO 2 mineralization in basalt. Using transmission electron microscopy, we mapped mineralogical variations from the initial to later stages of subsurface carbonate growth and identified a previously unknown cation-ordered ankerite phase that exerts a primary control over carbonation processes. This study has provided a new understanding of subsurface carbonation pathways which will impact the parameterization of predictive geochemical models for future sequestration efforts in basalt formations.

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

Aperture and roughness govern iron oxide passivation in olivine fractures during carbon mineralization

Efficient carbon dioxide mineralization in silicate-rich formations depends on reaction rates within fracture networks, yet the role of fracture microstructure remains poorly constrained. We investigated how fracture aperture and surface roughness influence mineralization by reacting synthetic forsteritic olivine fractures with controlled apertures (0.5–1.6 millimeters) and roughness (6–16 micrometers) with supercritical carbon dioxide at 13.7 megapascals and 185 degrees Celsius. Optical microscopy and Raman spectroscopy showed that rougher fracture surfaces promoted iron oxide precipitation, which passivated olivine surfaces, limiting magnesium availability for carbonate formation while removing iron as a reactant. These effects were more pronounced in smaller-aperture fractures. A geochemical model incorporating surface topography and passivation reproduced the observed trends, confirming that increased roughness can reduce reactivity despite greater surface area. Our results demonstrate that fracture microstructure strongly influences iron oxide passivation and carbon dioxide mineralization, with important implications for predicting long-term performance of engineered carbon dioxide storage systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part II: hydrolysis of Er from 25 to 75 °C

Aqueous speciation of rare earth elements (REE) controls their mobilization, fractionation, and enrichment in the natural waters. Geochemical modeling of their speciation is key to improve our understanding of the formation of economic mineral deposits, for developing mineral separation and mine tailing recovery technologies, and for characterizing the geochemistry of thermal water. However, our ability to predict the fate of REE in a wide pH and temperature range is limited by the scarcity of thermodynamic data for the REE hydroxyl complexes. In part I of this study (H. J. Han and A. P. Gysi, Dalton Trans., 2024, 53, 13129–13141), the optical properties of m-cresol purple (mCP) were determined using UV-Vis spectrophotometry between 25 and 75 °C in order to develop a method for deriving the hydrolysis constants of erbium (Er). Here, UV-Vis spectrophotometry experiments were conducted as a function of temperature between 35 and 75 °C to determine the hydrolysis of Er in near-neutral to alkaline solutions using mCP as an in situ pH color indicator. Here, the experiments were conducted with Er concentrations from 0 to ~0.253 mmol kg –1 in low ionic strength solutions (≤0.001 mol kg –1 ).

58 GEOSCIENCES↗
DOE OSTI2023

Inorganic Hydrogeochemistry in the 21st Century

Abstract Chemical and isotopic processes occur in every segment of the hydrological cycle. Hydrogeochemistry—the subdiscipline that studies these processes—has seen a transformation from “witch's brew” to credible science since 2000. Going forward, hydrogeochemical research and applications are critical to meeting urgent societal needs of climate change mitigation and clean energy, such as (1) removing CO 2 from the atmosphere and storing gigatons of CO 2 in soils and aquifers to achieve net‐zero emissions, (2) securing critical minerals in support of the transition from fossil fuels to renewable energies, and (3) protecting water resources by adapting to a warming climate. In the last two decades, we have seen extensive activity and progress in four research areas of hydrogeochemistry related to water‐rock interactions: arsenic contamination of groundwater; the use of isotopic and chemical tracers to quantify groundwater recharge and submarine groundwater discharge; the kinetics of chemical reactions and the mineral‐water interface's control of contaminant fate and transport; and the transformation of geochemical modeling from an expert‐only exercise to a widely accessible tool. In the future, embracing technological advances in machine learning, cyberinfrastructure, and isotope analytical tools will allow breakthrough research and expand the role of hydrogeochemistry in meeting society's needs for climate change mitigation and the transition from fossil fuels to renewable energies.

Zhu, Chen (ORCID:0000000153746787)↗