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.