One-atmosphere experiments in the system diopside-albite-anorthite (Di-Ab-An) have been used to explore the controls on plagioclase-melt partitioning of trace elements (Li, K, Rb, Cs, Mg, Zn, Sr, Ba, La, Sm, Y). By performing experiments along isotherms it was possible to isolate compositional controls from those of temperature. The lattice strain model accounts for variations in partition coefficients (<i>K</i><sub><i>p</i>(i)</sub>) for isovalent series of cations (<i>i</i><sup><i>n</i>+</sup>) and the experimental data are used to optimize the lattice strain parameters so that partition coefficients for <i>n</i> = 1 to 3 can be derived from <i>K</i><sub><i>p</i>(Na)</sub>, <i>K</i><sub><i>p</i>(Ca)</sub>, and <i>K</i><sub><i>p</i>(La)</sub>, respectively, at the conditions of interest. The optimized lattice strain parameters have magnitudes and compositional and thermal dependences that are consistent with elastic data for plagioclases. <i>K</i><sub><i>p</i>(La)</sub> is parameterized via an exchange reaction involving CaAl<sub>2</sub>Si<sub>2</sub>O<sub>8</sub> and the La-feldspar component La<sub>0.5</sub>Na<sub>0.5</sub>Al<sub>2</sub>Si<sub>2</sub>O<sub>8</sub>. The free energy change of this reaction is small, but well constrained by the data, such that for a given <i>P-T-X K</i><sub><i>p</i>(La)</sub> can be derived from a combination of <i>K</i><sub><i>p</i>(Na)</sub> and <i>K</i><sub><i>p</i>(Ca)</sub>. In the case where <i>K</i><sub><i>p</i>(Na)</sub> and <i>K</i><sub><i>p</i>(Ca)</sub> are independently constrained, the models reproduce the experimental values of 749 <i>K<sub>p</sub></i>s for Li, Mg, K, Rb, Sr, Y, Cs, Ba, and REE (obtained from 91 individual experiments varying over five orders of magnitude) within a factor of 1.5 for 53 percent of the data and a factor of 2 for 74 percent of the data. The isothermal experimental data reveal that additional compositional effects, arising from non-ideal mixing in plagioclase and in melt, also influence <i>K</i><sub><i>p</i>(<i>i</i>)</sub>. A thermodynamic model based on the free energy of fusion of albite (NaAlSi<sub>3</sub>O<sub>8</sub>) and anorthite (CaAl<sub>2</sub>Si<sub>2</sub>O<sub>8</sub>), and taking into account plagioclase and melt non-ideality in Di-Ab-An, is developed to express <i>K</i><sub><i>p</i>(Na)</sub>, <i>K</i><sub><i>p</i>(Ca)</sub> as functions of pressure (<i>P</i>), temperature (<i>T</i>) and composition (<i>X</i>). The thermodynamic models are extended to a wider compositional range by means of an empirical parameterization of 919 <i>K</i><sub><i>p</i>(Na)</sub> and 920 <i>K</i><sub><i>p</i>(Ca)</sub> from a database of 991 published experiments. These models are able to reproduce 919 measured partitioning values for Li, Na, Mg, Ca, K, Rb, Sr, Y, Cs, Ba, and REE (from 91 experiments) with 47 percent of the model values within a factor of 1.5 of the measured values and 69 percent within a factor of 2 for any <i>P-T-X</i> condition. The models can also be used to constrain the dependence of the chemical potential of trace cations in plagioclase solid solutions (μ<sub>i</sub><sup>pl</sup>) for use in diffusion chronometry where the fully equilibrated trace element concentration profile must be known in order to calculate timescales. It is shown that μ<sub>i</sub><sup>pl</sup> is not adequately described by existing empirical models used to predict <i>K<sub>p</sub></i> as a function of anorthite because of the conflating effects of <i>T</i> and <i>X</i> in the polythermal experimental datasets from which those models were derived.
Read full abstract