# Apatite saturation revisited: new model formulations and applications to igneous rocks

**Authors:** Benjamin Z. Klein, Othmar Müntener, Jack Gillespie, Felix Marxer

PMC · DOI: 10.1007/s00410-026-02300-5 · Contributions to Mineralogy and Petrology. Beitrage Zur Mineralogie Und Petrologie · 2026-02-19

## TL;DR

This paper introduces new models to predict when apatite forms in igneous rocks, improving accuracy and understanding of its behavior under various conditions.

## Contribution

The paper presents two new, accurate models for apatite saturation in igneous melts, outperforming existing methods.

## Key findings

- The first model predicts apatite saturation temperature with ~32°C accuracy, surpassing existing models.
- The second model estimates melt P2O5 content at apatite saturation with better than a factor of two accuracy.
- Apatite stability is not significantly affected by volatile species or pressure in the tested ranges.

## Abstract

Apatite is the primary phosphate mineral in the Earth’s crust and is present in a range of intrusive and extrusive igneous rocks, typically at minor to trace quantities. Similarly, apatite is a stable phase in many equilibrium crystallization experiments conducted with phosphorus-bearing starting compositions. We leverage this experimental stability to produce a large compilation of apatite-saturated liquid compositions, supplemented by additional apatite trace element and volatile partitioning and explicit apatite solubility experiments, as well as analyses of natural rhyolitic glasses. Using this compilation, we calibrate two new, independent models: for apatite saturation temperature as a function of melt P2O5 and SiO2 contents and aluminum saturation index (ASI; molar Al/(2Ca + Na + K)); and for melt P2O5 contents at apatite saturation as a function of temperature, melt SiO2 contents and ASI. The first model reproduces apatite saturation temperatures with an accuracy of ~ 32 °C, significantly outperforming existing apatite saturation models as well as recent zircon saturation thermometers. The second model reproduces melt P2O5 contents at apatite saturation by better than a factor of two across four orders of magnitude. Our new calibrations show that, within uncertainty, apatite stability does not differ regardless of the specific volatile species present (H2O, F or Cl) or on the quantity of H2O dissolved in the melt. Further, we find that apatite stability is not sensitive to pressures at least in the range of 1 atm to 2 GPa. This model accurately describes the saturation of apatite in a wide range of liquids, including metaluminous liquids, moderately alkaline liquids, and most peraluminous liquids. Experimental and natural peraluminous liquids with unusually high P2O5 contents that are not well described by our model contain CaO:P2O5 ratios below apatite stoichiometry (‘perphosphorous melts’), indicating that apatite saturation in these liquids is at minimum jointly controlled by CaO and P2O5 contents.

The online version contains supplementary material available at 10.1007/s00410-026-02300-5.

## Linked entities

- **Chemicals:** P2O5 (PubChem CID 14812), SiO2 (PubChem CID 24261), H2O (PubChem CID 962), F (PubChem CID 24524), Cl (PubChem CID 312)

## Full-text entities

- **Chemicals:** amphibole (MESH:D017636), P2O5 (MESH:C012500), C (MESH:D002244), Lu (MESH:D008187), carbonate (MESH:D002254), silicate (MESH:D017640), basalt (MESH:C060346), P (MESH:D010758), phosphate (MESH:D010710), Granite (MESH:C007886), oxygens (MESH:D010100), Nd (MESH:D009354), H2O (MESH:D014867), Apatite (MESH:D001031), Sm (MESH:D012493), F (MESH:D005461), FeO (MESH:C034236), Li (MESH:D008094), Fe (MESH:D007501), Al2O3 (MESH:D000537), andalusite (MESH:C121090), Cl (MESH:D002713), whitlockite (MESH:C021767), ASI (-), Al (MESH:D000535), sulfur (MESH:D013455), Zr (MESH:D015040), He (MESH:D006371), SiO2 (MESH:D012822), K (MESH:D011188), Na (MESH:D012964), Pb (MESH:D007854), halogen (MESH:D006219), Mn (MESH:D008345), K2O (MESH:C068440), Mg (MESH:D008274), CaO (MESH:C016538), Ca (MESH:D002118), Y (MESH:D015019), Th (MESH:D013910), zircon (MESH:C003784), CO2 (MESH:D002245), rutile (MESH:C009495), Na2O (MESH:C096707), alkali (MESH:D000468), MgO (MESH:D008277), U (MESH:D014501), Hf (MESH:D006195), monazite (MESH:C015370), OH (MESH:C031356)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12920737/full.md

## Figures

11 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12920737/full.md

---
Source: https://tomesphere.com/paper/PMC12920737