Phosphorous-solubility in carbonatite melts: Apatite crystallization modeled via its solubility product
Gino Sartori, Max W. Schmidt
Abstract
We model apatite-saturation in carbonatite melts based on a compilation of experimental data ranging from 650 to 1430 °C and 1 to 60 kbar. The data show a very strong correlation of inverse temperature with the apatite solubility product, a relation expressed by the equation. ln(CaO5·PO2.53)=-27450/T+5.79 (T in Kelvin). Within the available dataset, F and Cl do not play a discernable role. Application of the solubility product to natural Ca-carbonatites indicates that a few rocks with >8 wt% P2O5 have cumulative apatite while most Ca-carbonatites (with typically <5 wt% P2O5) are apatite undersaturated at their liquidus temperatures, defined by calcite crystallization. To address true carbonatite liquids, we model calcite fractionation and melt evolution for natural rock compositions with 5, 10 and 20 mol% H2O and/or (Na,K)2CO3 added, 5% representing the lower bound for any carbonatite formation model. Both H2O or (Na,K)2CO3 cause very similar liquidus depressions of ∼10 °C/mol%. The model result is that saturation of apatite occurs in most natural carbonatite melts only after >45, 30–55, and 10–30 mol% calcite-fractionation for 5, 10, and 20 mol% fluxing components added, respectively. We further estimate the melt fractions necessary to dissolve all apatite in carbonatite melts generated from carbonated MORB and pelites, opening the discussion on an unlikely restitic nature of subducted apatites. In both the crystallization and forward melting cases, apatite crystallization or dissolution is mostly governed by temperature, surprisingly, carbonatite melt evolution through calcite-fractionation has a minor influence on the solubility product.