Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrothermal fluids”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Selective transport of light vs. heavy rare earth elements by sulfate/bisulfate complexes in hydrothermal fluids

Here, this study explores the transport of rare earth elements (REE) in acidic sulfate-bearing hydrothermal fluids and the implications for the fractionation of light/heavy REE in critical mineral deposits. The speciation of Nd (light REE) and Yb (heavy REE) sulfate complexes were determined via in situ Raman spectroscopy using fused SiO 2 capillary cells up to 300 ºC at saturated water vapor pressure and in a hydrothermal diamond anvil cell up to 500 ºC and 540 MPa. The REE monosulfate (REESO 4 + ) and REE disulfate (REE(SO 4 ) 2 - ) species are stable below 150 to 250 ºC but become less stable at higher temperatures, particularly the light REE, due to the decreased solubility of REE sulfate solids. At higher pressure, these REE sulfate complexes display an increased stability field up to 400 ℃. The REE bisulfate complex (REEHSO 4 2+ ) was identified with a wide stability field below 400 ºC for the HREE in acidic Yb 2 (SO 4 ) 3 -bearing solutions, whereas in Nd 2 (SO 4 ) 3 -bearing solutions, the LREE bisulfate complex is restricted to below 100 ℃. These results suggest that bisulfate is a previously unrecognized selective ligand for heavy REE transport at low temperature in acidic oxidized hydrothermal fluids. Such fluids are responsible for hydrothermal alteration in many REE deposits. Prediction of phase stabilities across pressure, temperature, and composition space (P-T- x ) is crucial for predicting the role of aqueous REE sulfate complexes in the mobilization of REE in crustal fluids.

58 GEOSCIENCES↗
DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗
DOE OSTI2024

Accurate and Data‐Efficient Micro X‐ray Diffraction Phase Identification Using Multitask Learning: Application to Hydrothermal Fluids

Traditional analysis of highly distorted micro X‐ray diffraction (μ‐XRD) patterns from hydrothermal fluid environments is a time‐consuming process, often requiring substantial data preprocessing and labeled experimental data. Herein, the potential of deep learning with a multitask learning (MTL) architecture to overcome these limitations is demonstrated. MTL models are trained to identify phase information in μ‐XRD patterns, minimizing the need for labeled experimental data and masking preprocessing steps. Notably, MTL models show superior accuracy compared to binary classification convolutional neural networks. Additionally, introducing a tailored cross‐entropy loss function improves MTL model performance. Most significantly, MTL models tuned to analyze raw and unmasked XRD patterns achieve close performance to models analyzing preprocessed data, with minimal accuracy differences. This work indicates that advanced deep learning architectures like MTL can automate arduous data handling tasks, streamline the analysis of distorted XRD patterns, and reduce the reliance on labor‐intensive experimental datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

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

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2025

In situ Raman investigation of Dy complexation in Cl-bearing aqueous solutions at 20–300 °C

Raman spectroscopy provides a versatile tool for in situ characterization of aqueous rare earth elements (REE) speciation at the molecular level. Complexation of REE with ligands such as Cl – and OH – is of particular interest for understanding the mobility of REE in NaCl-bearing hydrothermal fluids responsible for enriching REE to economic levels in nature. Raman spectroscopic studies of REE speciation in Cl-bearing aqueous fluids are primarily conducted at ambient temperature, whereas natural systems indicate temperatures of >100–600 °C. In this study, the speciation of Dy in acidic chloride-bearing hydrothermal solutions was investigated using confocal Raman spectroscopy with a new capillary Raman heating stage at 20–300 °C. Background solutions (pure water, NaCl-solutions) and solutions with 0.14–1.8 mol kg –1 dissolved DyCl 3 were sealed in quartz capillary cells. Comparison of the spectra for Dy chloride solutions with those for background solutions and the spectra for reference Dy-bearing solids was used to identify Raman bands specific to Dy–O and Dy–Cl bonds. The Raman band for the Dy–O stretching mode of hydrated Dy 3+ aqua ions was measured at 365–384 cm –1 and a Raman band for the Dy–Cl stretching modes of Dy chloride complexes was measured near 240 cm –1 . The Dy–O band decreases systematically with temperature, whereas the Dy–Cl band systematically increases, indicating a systematic increase in the stability of Dy chloride complexes with temperature. Here, this study provides the framework for expanding the use of in situ Raman spectroscopy to investigate the speciation of REE in aqueous solutions to hydrothermal conditions.

58 GEOSCIENCES↗
DOE OSTI2025

Stability of aqueous neodymium complexes in carbonate-bearing solutions from 100–600 °C

Rare earth element exploration requires a quantitative understanding of factors governing their mobilization and economic concentration. However, the behavior of rare earth elements in carbonate- bearing hydrothermal fluids associated with carbonatite-hosted deposits is poorly understood, and conflicting mechanisms of rare earth transport by anionic ligands and alkali behavior have been described. Here, we report quantitative data to characterize the role of carbonate-bearing solutions in the hydrothermal mobilization of neodymium. Solubility studies of neodymium phosphate were performed at temperatures ranging from 100 to 600 °C in carbonate-bearing solutions. The thermodynamic data determined for the predominant complex were used to model the separation of neodymium from thorium in a simple flow-through system based on fluid and mineral compositions characteristic of carbonatite deposits. Our data suggest that neodymium transport is controlled by the stability of the carbonate species NdCO 3 OH o , and at temperatures of 500–600 °C, the concentrations of neodymium in solutions can reach ~1000 ppm.

58 GEOSCIENCES↗
DOE OSTI2024

Subsurface microbial community structure shifts along the geological features of the Central American Volcanic Arc

Subduction of the Cocos and Nazca oceanic plates beneath the Caribbean plate drives the upward movement of deep fluids enriched in carbon, nitrogen, sulfur, and iron along the Central American Volcanic Arc (CAVA). These compounds fuel diverse subsurface microbial communities that in turn alter the distribution, redox state, and isotopic composition of these compounds. Microbial community structure and functions vary according to deep fluid delivery across the arc, but less is known about how microbial communities differ along the axis of a convergent margin as geological features (e.g., extent of volcanism and subduction geometry) shift. Here, we investigate changes in bacterial 16S rRNA gene amplicons and geochemical analysis of deeply-sourced seeps along the southern CAVA, where subduction of the Cocos Ridge alters the geological setting. We find shifts in community composition along the convergent margin, with communities in similar geological settings clustering together independently of the proximity of sample sites. Microbial community composition correlates with geological variables such as host rock type, maturity of hydrothermal fluid and slab depth along different segments of the CAVA. This reveals tight coupling between deep Earth processes and subsurface microbial activity, controlling community distribution, structure and composition along a convergent margin.

Science & Technology - Other Topics↗
DOE OSTI2025

"Hidden" hydrothermal technical potential & technoeconomics: Revealing permeability & fluids with more data

Historical hydrothermal estimates have largely relied on temperature or heat flow estimates ignoring the need for natural flowing fluids. More accurate hydrothermal estimates require some indication of permeability and fluids that naturally exist in the subsurface. This paper describes a novel approach that includes proxies of permeability and fluids in hydrothermal estimates by leveraging the relatively data-rich Great Basin. Specifically, nameplate capacities (megawatts) of operating geothermal plants, negative (0 megawatt) locations and 48 geophysical and geologic features are used to used in eXtreme Gradient Boosting (XGBoost) regression to make hydrothermal capacity predictions. Additionally, this work inputs the XGBoost-based hydrothermal predictions into the Renewable Energy Potential (reV) model to quantify technical capacity, its uncertainty and techno-economics. Compared to historical hydrothermal estimates, these predictions adhere to the 37 operating geothermal plants and negative locations. We present a method for subsampling the negative sites to bring the labels into balance that uses the geologic domain knowledge to proportionally represent negatives. Overall, the distributions of the hydrothermal technical capacity and the site levelized cost of energy are respectively much tighter, lower and more accurate than the previous estimates for the Great Basin, as they include geological and geophysical surrogates for permeability and fluids. Percentile (50th and 90th, median and high estimate, respectively) models provide bookends for these metrics.

13 HYDRO ENERGY↗
DOE OSTI2024

Co-Processing Fast Pyrolysis Bio-Oils and Hydrothermal Liquefaction Biocrudes in Fluid Catalytic Cracking and Hydroprocessing in Refineries

Co-processing of biogenic feedstocks within the existing petroleum infrastructure represents an opportunity to incorporate large volumes of bio-carbon into transportation fuels. However, upgrading intermediate liquids from biomass and waste to final products efficiently and economically, while minimizing impacts to existing refinery infrastructure, remains a notable barrier to commercialization. The co-processing project, Bio-oil Co-processing with Refinery Streams, aimed to generate foundational knowledge for processing renewable intermediates in petroleum refineries and offer co-processing strategies, as well as develop new methods for biogenic carbon tracking and measurement.

bio-oil↗
DOE OSTI2025

Spectrophotometric determination of the stability of La hydroxyl complexes at near neutral to alkaline pH from 25 to 75 °C

The hydrolysis of rare earth elements (REE) potentially controls their mobility during fluid-rock interaction in a broad range of pH and temperature conditions. However, there is still a lack of thermodynamic data for modeling accurately the stability of REE hydroxyl complexes in hydrothermal aqueous fluids. Here, in this study, UV–Vis spectrophotometric experiments were conducted from 25 to 75 °C in near-neutral to alkaline NaOH-bearing aqueous solutions with varying lanthanum (La) concentrations (0 to ∼0.23 mmol/kg). The color indicator m-cresol purple was used to determine in situ pH and derive the average OH− ligand number ($\overrightarrow{n}$) and formation constants for the La hydroxyl complexes (LaOH 2+ , La(OH) 2 + , and La(OH) 3 0 ). From 25 to 50 °C, $\overrightarrow{n}$ ranges between ∼1 and 2 at pH from 7.0 to 9.3. At 75 °C, $\overrightarrow{n}$ ranges between ∼1.5 and 3 at pH from 6.3 to 8.8. These results suggest the predominance of LaOH 2+ and La(OH) 2 + complexes from 25 to 50 °C, and an increased predominance of La(OH) 3 0 at 75 °C. The cumulative formation constants (β n °, n = 1 to 3) are derived for the reaction La 3+ + nOH − = La(OH) n 3-n , and fitted between 25 and 250 °C by combining the UV–Vis and literature solubility data. The resulting logβ n ° are expressed as function of temperature (T in Kelvin): logβ 1 ° = −1.786 + 0.0133 T + 1.049·10 3 /T; logβ 2 ° = −5.797 + 0.0267 T + 2.713·10 3 /T; logβ 3 ° = 6.435 + 0.0223 T + 512.7/T. A comparison between these new fits and existing extrapolations using the Helgeson-Kirkham-Flowers equation of state indicates significant differences in the predicted hydrolysis of La. The latter extrapolations should therefore be updated for the hydrolysis of REE.

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

Evidence for metal sources, fluid-mixing processes, and S isotope recycling within the feeder zone of an Irish type Zn-Pb deposit

The origin and evolution of fluids in Irish-type Zn-Pb deposits remains debated, particularly regarding the mobility of metals such as Cu and Ni, sources of sulphur, and the role of fluid mixing and replacement. The Lisheen Zn-Pb deposit, Ireland, offers a well-defined natural laboratory to investigate these questions. While most studies have focused on the Waulsortian Limestone Formation, the primary sulphide host, less is known about mineralisation in underlying units, such as the Lisduff Oolite Member (LOM). The LOM displays enrichment in Cu and Ni and displays intense replacement textures compared to other hosts at Lisheen, making it an ideal target for studying metal mobility and sulphur recycling in carbonate-hosted systems. Through characterising and studying LOM-hosted sulphides, valuable insights into mineralisation processes, especially related to Cu-Ni metals, can be defined. This study integrates petrography, EMPA, and in situ sulphur isotope (δ 34 S) analysis to investigate sulphide paragenesis, mineral chemistry, and fluid evolution across LOM ore zones. Results reveal a multistage mineralising system involving extensive replacement of early pyrite (Py0, δ 34 S = −28.4 to −21.9 ‰) by sphalerite and galena, with zoned pyrite (Py1) enriched in As-Cu-Ni-Tl. The δ 34 S values and trace element trends indicate mixing between hydrothermal and bacteriogenic sulphur-rich fluids, with evidence for sulphur recycling during replacement. Pyrite textures and compositions capture this evolving fluid regime, with trace element enrichment linked to paragenetic stage. The steel ore region, adjacent to major fault intersections, records intense hydrothermal fluid interaction, hosting Ni- and As-rich phases such as nickeline, gersdorffite, and arsenopyrite. These findings highlight the importance of structural controls and fluid mixing in metal transport and deposition, positioning the LOM as a key stratigraphic unit for understanding ore-forming processes in Irish-type systems. These results have implications for targeting similar carbonate-hosted systems globally, especially where deeper or structurally complex ore zones remain underexplored.

58 - GEOSCIENCES↗
DOE OSTI2026

Chromium-for-Aluminum Substitution in Synthetic Serpentine

Cr-bearing clay minerals are products of hydrothermal alteration and fluid–rock interactions of ultramafic rocks that form serpentine minerals. Cr is typically observed to substitute for Al in serpentine minerals, but the crystal chemistry and environmental constraints on this substitution are unknown. Here, we synthesized endmember and Cr-substituted amesite, a typical Al-serpentine mineral, via the hydrothermal method. We found that the phase purity highly depends on the pH of the hydrothermal solution, which should be controlled at ~12.7 to avoid the formation of impurity phases. Additionally, amesite can incorporate Cr at a concentration equivalent to ~39.5% substitution of Al. The Cr-free and Cr-substituted amesite are highly defective and contain multiple polytypes, including 6R 2 , 2M 1 , and possibly 2H 2 . However, the relative proportions of these polytypes do not change with increasing chromium substitution.

Al-lizardite↗
DOE OSTI2026

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2025

The Effect of Hydrothermal Alteration and Microcracks on Hydraulic Properties and Poroelastic Deformation: A Case Study of the Blue Mountain Geothermal Field

Abstract Geothermal energy plays a vital role in decarbonizing electricity and heat supply. Effective utilization of geothermal resources hinges on identifying or generating permeable reservoir zones and understanding how effective pressure variations affect fluid circulation and reservoir properties by poroelastic deformation. Hydrothermal alteration can modify the petrophysical properties of geothermal reservoir rocks, which may increase or decrease its productivity. Understanding these alteration effects is essential to predict and optimize long‐term sustainable geothermal operations. Here, we investigate the impact of hydrothermal alteration on poroelastic and hydraulic properties of diverse lithologies in a series of deformation tests performed at several confining (0–80 MPa) and pore pressure (10–30 MPa) levels. Experimental results of hydrothermally altered dikes and phyllites obtained from the Blue Mountain geothermal field (Nevada, USA) are compared to thermally cracked La Peyratte granite (France) and correlated with petrophysical properties, mineral composition, and microstructures. Argillic alteration of dikes increases porosity and storage capacity but lowers thermal conductivity and increases pore compressibility. Conversely, silicate precipitation in phyllites increases stiffness and thermal conductivity but also reduces porosity and permeability. Experimentally determined effective pressure coefficients range from 0.1 to 0.9, differ for permeability and volumetric strain and decrease with increasing effective pressure. The presence of compliant microcracks and crack‐like pores significantly increases the stress sensitivity of La Peyratte granite and silicified phyllites. This study demonstrates how thermal and chemical alteration impacts poromechanical and petrophysical characteristics of geothermal targets, which ultimately govern reservoir stability and subsidence, induced seismicity as well as fluid and heat extraction efficiency during geothermal operations.

Schuster, Valerian [Helmholtz Centre Potsdam GFZ G↗