## Abstract

We present a quantitative model to describe the partitioning of rare earth elements (REE) and Y between clinopyroxene and anhydrous silicate melt as a function of pressure (P), temperature (T) and bulk composition (X). The model is based on the Brice (1975) equation, which relates the partition coefficient of element i (D_{i}) to that element of o (D_{o}) where the latter has the same ionic radius r_{o} as the crystallographic site of interest, in this case the clinopyroxene M2 site: [see original article for formula]. N_{A} is Avogadro's number, E_{M2} is the Young's Modulus of the site, R is the gas constant and T is in K. Values of E_{M2} obtained by fitting the Brice equation to experimental REE partition coefficient patterns are in good agreement with those obtained from the well-known correlation between bulk modulus, metal-oxygen distance and cation charge. Using this relationship to constrain E_{M2} for 3 + cations and then fitting the Brice equation to those experimental data where 3 or more REE partition coefficients had been simultaneously measured we obtained 82 values of D_{o} and r_{o}. The latter was found to be a simple and crystallochemically reasonable function of clinopyroxene composition. We show that for any clinopyroxene-melt pair if D for one middle REE (e.g. Sm or Gd) is known then the Brice equation can be used to predict Ds for all the other REE, with uncertainties similar to those involved in the actual measurements. The model was generalised using thermodynamic descriptions of REE components in crystal and melt phases to estimate the free energy of fusion (ΔG_{f}) of the fictive REE components REEMgAlSiO_{6} and Na_{0.5}REE_{0.5}MgSi_{2}O_{6}. For the melt we find that 6-oxygen melt components (CaMgSi_{2}O_{6}, NaAlSi_{2}O_{6}, Mg_{3}Si_{1.5}O_{6} etc) mix with constant activity coefficient over a wide range of natural compositions. Propagating ΔG_{f} into the Brice model we obtain an expression for D^{3+} _{o}in terms of the atomic fraction of Mg on the clinopyoroxene M1 site, the Mg-number of the melt, P and T. The D for any REE can be calculated from D^{3+} _{o} using the Brice equation. Over 92% of D_{REE} (454 points) calculated in this way lie within a factor 0.63-1.59 of the experimental value. The approach can be extended to calculate D for any REE at a given P (≤ 6GPa) and T (1200-2038 K) to within 0.60-1.66 times the true value given only the crystal and/or melt composition. The model has widespread applicability to geochemical modelling of all natural processes involving clinopyroxene, e.g. decompression mantle melting, enabling for the first time account to be taken of variations in partition coefficient in response to changing pressure, temperature and phase composition.

Original language | English |
---|---|

Pages (from-to) | 166-181 |

Number of pages | 16 |

Journal | Contributions to Mineralogy and Petrology |

Volume | 129 |

Issue number | 2-3 |

DOIs | |

Publication status | Published - 1997 |