skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, December 13 until 2:00 AM ET on Saturday, December 14 due to maintenance. We apologize for the inconvenience.


This content will become publicly available on March 1, 2025

Title: A new experimental monazite-xenotime thermometer: Application to metamorphic environments
Monazite and xenotime are common metamorphic phases that may be reliably U/Th-Pb dated to obtain absolute ages. The utility of these minerals is significantly enhanced if a crystallization thermometer can be developed and applied to better establish the temperature-time (T-t) paths of crustal terranes. Here we report experiments in which we have modeled the T-dependent Rare Earth element (REE) cationic exchange between coexisting monazite and xenotime to derive a new thermometer. We present a thermometer in which phosphates were cocrystallized from 1150 ◦C to 850 ◦C at 1 GPa in a Y-REE-P2O5-NaCl-H2O system, with oxygen fugacity buffered at the Ni-NiO equilibrium. The composition of the phosphates was quantified using a laser ablation inductively coupled plasma mass spectrometer (LA-ICP-MS). Results reveal strong correlations between log10 (Π XLREE Xnt ×XY Π Mzt XLREE Mzt ×XY Mzt ) (at. %) and 104/T(K) and our preferred calibration is: log (ΠXLREE Xnt × XY Π Mzt XLREE Mzt × XY Mzt ) = ( 􀀀 6714 ± 2264) T (K) 􀀀 (0.79±1.76) where LREE = La, Ce and Pr, α = activity of a cation in a phase, and ΠαY/REE Mzt/Xnt refers to the product of activities of Y and/or REEs in a phosphate phase. The errors are 2 s.e. The greatest strength of this thermometer is its versatility. One can obtain derivative thermometers based on select elements rather than the entire suite of REEs. We showcase our thermometer’s adaptability by applying it to two studies that have published REE data on monazite and xenotime from some quartzo-feldspathic psammites and garnet-bearing pelites that experienced amphibolite facies metamorphism from the Naver nappe in the Northern Highlands Terrane, Scotland. The main calibration shown above, as well as four derivative single-element thermometers (Y, La, Ce and Pr) were applied to the first study. Upon applying these thermometers, we find that the calculated metamorphic Ts agree well with the regional metamorphic facies. Thus, this versatile thermometer can be used in geologic environments where monazite and xenotime co-crystallized.  more » « less
Award ID(s):
1751903 2327940
PAR ID:
10511588
Author(s) / Creator(s):
;
Corporate Creator(s):
Editor(s):
Romano, C
Publisher / Repository:
Elsevier (Chemical Geology)
Date Published:
Journal Name:
Chemical Geology
Edition / Version:
1
Volume:
648
Issue:
C
ISSN:
0009-2541
Page Range / eLocation ID:
121939
Subject(s) / Keyword(s):
Rare-earth phosphates Monazite Xenotime Thermometry Metamorphism
Format(s):
Medium: X Size: 2MB Other: pdf
Size(s):
2MB
Sponsoring Org:
National Science Foundation
More Like this
  1. Charlier, B (Ed.)
    Quantifying the oxygen fugacity (fo2) of high temperature lithospheric fluids, including hydrothermal systems, presents a challenge because these fluids are difficult to capture and measure in the same manner as quenched glasses of silicate melts. The chemical properties of fluids can however be inferred through mineral proxies that interacted with the fluids through precipitation or recrystallization. Here, we present hydrothermal experiments to quantify the partition coefficients of rare earth elements (REEs) – including redox-sensitive Ce and Eu – between zircon and fluid. Experiments were conducted in a piston cylinder device at temperatures that range from 1200 to 800 ◦C under fo2-buffered conditions in a SiO2-ZrO2-NaCl-REE-oxide system, and similar experiments were performed in the absence of NaCl (31 total experiments). The fo2 was buffered to values that range from approximately 3 log units below to 7 log units above the fayalite magnetite quartz equilibrium. Zircon REE concentrations were quantified using laser ablation inductively coupled plasma mass spectrometry whereas the quenched fluids were extracted and measured by solution-based inductively coupled plasma mass spectrometry. Zircon Ce anomalies, quantified relative to La and Pr, exhibit sensitivity to oxygen fugacity and temperature and our preferred calibration is: log [􀀀 Ce Ce* ) D 􀀀 1 ] = (0.237 ± 0.040)× log(fo2) + 9437±640 T(K) 􀀀 5.02 ± 0.38 where the Ce anomalies are calculated from the partition coefficients for La, Ce, and Pr. Zircon Eu anomalies are also a function of oxygen fugacity though they exhibit no systematic dependence on T. Our preferred calibration is described by: 􀀀 Eu Eu* ) D = 1 1+100.30±0.04􀀀 [0.27±0.03]×ΔFMQ We performed additional calculations, in which lattice strain parabolas were fit to all non-redox sensitive rare earth elements that were added to the starting composition (i.e., La, Pr, Sm, Gd, Dy, Ho, Tm, Lu) as an alternate means to calculate anomalies. This method yields broadly similar results, though we prefer the La-Pr calibrations due to the non-systematic REE patterns frequently encountered with hydrothermal zircons; e.g., LREE zircon enrichment relative to other REEs. These experiments are applied to quantify the fo2 of fluids during mineralization of critical element-bearing systems, and separately to calculate the oxygen fugacity values of fluids formed during plate boundary processes. 
    more » « less
  2. Monazite-(Ce) and xenotime-(Y) occur as secondary minerals in iron-oxide-apatite (IOA) deposits, and their stability and composition are important indicators of timing and conditions of metasomatism. Both of these minerals occur as replacement of apatite and display slight but important variations in light (e.g. La, Ce, Pr, Nd, etc.) and heavy (e.g. Y, Er, Dy, Yb, etc.) REE concentrations [1,2]. The causes for these chemical variations can be quantified by combining thermodynamic modeling with field observations. Major challenges for determining the stability of these minerals in hydrothermal solutions are the underlying models for calculating the thermodynamic properties of REE-bearing mineral solid solutions and aqueous species as a function of temperature and pressure. The thermodynamic properties of monazite and xenotime have been determined using several calorimetric methods [3], but only a few hydrothermal solubility studies have been undertaken, which test the reliability and compatibility of both the calorimetric data and thermodynamic properties of associated REE aqueous species [4,5]. Here, we evaluate the conditions of REE metasomatism in the Pea Ridge IOA-REE deposit in Missouri, and combine newly available experimental solubility data to simulate the speciation of LREE vs. HREE, and the partitioning of REE as a function of varying fluid compositions and temperatures. Our new experimental data will be implemented in the MINES thermodynamic database (http:// tdb.mines.edu) for modeling the chemistry of crustal fluid-rock equilibria [6]. [1] Harlov et al. (2016), Econ. Geol. 111, 1963-1984;[2] Hofstra et al. (2016), Econ. Geol. 111, 1985-2016; [3] Navrotsky et al. (2015), J. Chem. Thermodyn. 88, 126-141; [4] Gysi et al. (2015), Chem. Geol. 83-95; [5] Gysi et al. (2018), Geochim. Cosmochim. Acta 242, 143-164; [6] Gysi (2017), Pure and Appl. Chem. 89, 581-596. 
    more » « less
  3. CePO4 and YPO4 are major components in monazite and xenotime, respectively, which are common hydrothermal phases in REE mineral deposits. Both minerals also occur as secondary minerals in iron-oxide-apatite deposits [1,2], and as accessory phases in high-grade metamorphic rocks where they display varying degrees of metasomatism. Studying the cause of their compositional variations using thermodynamic modeling may provide geochemical signals for interpreting P- T-x of crustal fluid-rock interaction. The thermodynamic properties of monazite and xenotime have been determined using several calorimetric methods [3], but only a few solubility studies have been undertaken, which test the reliability of both the calorimetric data and thermodynamic properties of associated REE aqueous species [4]. Combining available calorimetric data with the REE aqueous species from Haas et al. [5], implemented in the Supcrt92 database [6], yields several orders of magnitude differences when compared with our solubility measurements. To reconcile these discrepancies, we have used the GEMS code package [7,8] and GEMSFITS [9] for parameter optimization, and re- evaluated the standard Gibbs energies for aqueous REE species, while maintaining consistency with available calorimetric measurement of the REE phosphates. This study points to a need to revise the thermodynamic properties of the REE hydroxyl species, which will have an impact on the calculated solubilities of the REE phosphates and our understanding of the mobility of REE in hydrothermal fluids. Our new experimental data will be implemented in the MINES thermodynamic database (http:// tdb.mines.edu) [10] for modeling the chemistry of crustal fluid-rock equilibria. 
    more » « less
  4. Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (ΔfG°298) of REE aqueous species (REE3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of ΔfG°298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the ΔfG°298 values of REE3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the ΔfG°298 values of both REE3+ and REEOH2+ by ∼ 2–10 kJ/mol and ∼ 15–31 kJ/mol, respectively. For xenotime, adjustments of ΔfG°298 values by 1–26 kJ/mol are only necessary for the REE3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE3+ species at <150 °C and the LREEOH2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE3+, REEOH2+, and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. The optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems. 
    more » « less
  5. Monazite (CePO4) is a light rare earth element (REE) phosphate occurring as accessory mineral in metamorphic, igneous and sedimentary rocks, and is also a common mineral in REE mineral deposits. Metasomatism of monazite yields important clues about fluid-rock interaction in the crust, in particular, because its compositional variations may enable us to determine conditions of mineralization. The thermodynamic properties of monazite have been determined using several calorimetric methods, but up to the present time only a few solubility studies have been undertaken, which test the reliability of both, the thermodynamic properties of the REE phosphates and associated REE aqueous species. In this study, we have measured the solubility of the monoclinic REE phosphate end-members CePO4, SmPO4, and GdPO4 in aqueous perchloric acid solutions at temperatures from 100 to 250 °C at saturated water vapor pressure (swvp). The solubility products (Ks0) were determined according to the reaction: REEPO4 = REE3+ + PO43−. Combining available calorimetric data for the REE phosphates with the REE aqueous species from the Supcrt92 (slop98.dat) dataset, yields several orders of magnitude differences when compared with our solubility measurements. We have investigated ways to reconcile these discrepancies and propose a consistent set of provisional thermodynamic properties for REE aqueous species and REE phosphates that reproduce our measured solubility values. To reconcile these discrepancies, we have used the GEMS code package and GEMSFITS for parameter optimization by adjusting the standard Gibbs energy of REE3+ and REEOH2+ at 25 °C and 1 bar. An alternative optimization could involve adjustment of the standard Gibbs energy of REEPO4(s) and REEOH2+. Independently of the optimization method used, this study points to a need to revise the thermodynamic properties of REEOH2+ and possibly other REE hydroxyl species in future potentiometric studies. These revisions will have an impact on calculated solubilities of REE phosphates and our understanding of the mobility of REE in natural hydrothermal fluids. 
    more » « less