Experimental determination of REE partition coefficients between amphibole and basaltic to andesitic liquids at high pressure

Amphibole/liquid partition coefficients for the REE(D^amph/liq_REE) obtained from natural rocks increase systematically with bulk rock compositional change from basalt to rhyolite and are higher for the middle to heavy REE. Five new experimentally determined sets of D_REE (La, Sm, “Eu²⁺”, Ho, Lu)and 4 published sets are similar to these data and provide values for use in geochemical modelling of magmatic processes involving amphibole, over a range of temperature, pressure and oxygen fugacity. The partition coefficients increase significantly with decreasing temperature, and increase slightly with increasing oxygen fugacity. Pressure does not appear to have a major effect, although one data set suggests increased pressure lowers D^amph/liq_REE in a basaltic composition.     

Experimental rare earth element partition coefficients for garnet,clinopyroxene and amphibole coexisting with andesitic and basaltic liquids

The partitioning of La, Sm, Dy, Ho and Yb between garnet, calcic clinopyroxene, calcic amphibole and andesitic and basaltic liquids has been studied experimentally. Glasses containing one or more REE in concentrations of 500–2000 ppm were crystallized at pressures of 10–35 kbar, and temperatures of 900–1520°C. Water was added to stabilize amphibole and to allow study of partition coefficients over wide temperature ranges. Major element and REE contents of crystal rims and adjacent glass were determined by EPMA, with limits of detection for individual REE of 100–180 ppm. Measured partition coefficients, D_REE^Cryst-liq, are independent of REE concentration over the concentration ranges used.D-values show an inverse dependence on temperature, best illustrated for garnet. At a given temperature, they are almost always higher for equilibria involving andesitic liquid. Garnet shows by far the greatest range of D-values, with e.g. D_La < 0.05 and D_Yb ~ 44 for andesitic liquid at 940°C. D_Yb falls to ~ 12 at 1420°C. D_Sm^Ga-liq also correlates negatively with temperature and positively with the grossular content of garnet. Patterns of D_ree^Cryst-Liq for calcic clinopyroxenes and amphiboles are sub-parallel, with D-values for amphibole generally higher. Both individual D-values and patterns for the crystalline phases studied are comparable with those determined for phenocryst-matrix pairs in natural dacites, andesites and basalts.D-values and patterns are interpreted in terms of the entry of REE³⁺ cations into mineral structures and liquids of contrasted major element compositions. The significance of the partition coefficients for models of the genesis of andesitic and Hy-normative basaltic magmas is assessed. Most magmas of these types in island arcs are unlikely to be produced by melting of garnet-bearing sources such as eclogite or garnet lherzolite.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

Partition coefficients for REE between garnets and liquids: implications of non-Henry's Law behaviour for models of basalt origin and evolution

Partition coefficients for the rare earth elements (REE) Ce, Sm and Tm between coexisting garnets and hydrous liquids have been determined at high pressure and temperatures (30 kbar and 1300 and 1500°C). Two synthetic systems were studied, Mg₃Al₂Si₃O₁₂-H₂O and Ca₃Al₂Si₃O₁₂-H₂O, in addition to a natural pyrope-bearing system.Deviations from Henry's Law behaviour occur at geologically relevant REE concentrations. At concentrations < 3 ppm Ce, < 12 ppm Sm, < 80 ppm Tm in pyrope and < 100 ppm Ce, < 250 ppm Sm, < 1000 ppm Tm in grossular (at 30 kbar and 1300°C), D^{garnet liquid}_REE increases as the REE concentration in the garnet decreases. At higher concentrations, D_REE is constant. D^{grossular liquid}_REE also constant when the garnet contains less than about 2 ppm Sm or Tm. The REE concentration at which D_REE becomes constant increases with increasing temperature, decreasing REE ionic radius and increasing Ca content of the garnet.Partitioning behaviour of Ce, Sm and Tm between a natural pyrope-rich garnet and hydrous liquid is analogous to that in the synthetic systems and substantiates the substitution model proposed by Harrison and Wood (1980).Values of D_REE^{garnet/liquid} for which Henry's Law is obeyed are systematically higher for grossular than for pyrope (D^{pyrope/liquid} = 0.067(Ce), 0.108(Sm), 0.155(Tm) and D^{grossular/Liquid} = 0.65(Ce), 0.75(Sm), 4.55(Tm).The implications of non-Henry's Law partitioning of REE for models of basalt petrogenesis involving garnet are far-ranging. Deviations from Henry's Law permit refinements to be made to calculated REE abundances once basic model parameters have been defined.     

Experimental measurements of zircon/melt trace-element partition coefficients

Zircon was grown from trace-element doped hydrous peralkaline rhyolite melts with buffered oxygen fugacities in cold-seal experiments at 0.1 and 0.2 GPa and 800 °C and piston-cylinder experiments at 1.5 GPa and 900-1300 °C. Zircon and glass were present in all run products, and small monazite crystals were present in eight of the 12 experiments. Average diameters of zircon crystals ranged from 5 to 20 μm at 800 °C to 30-50 μm at 1300 °C. Zircon crystals have thin rims, and adjacent glass has a narrow (∼1 μm thick) compositional boundary layer. Concentrations obtained through in-situ analysis of cores of run product zircon crystals and melt pools were used to calculate trace-element partition coefficients D^zircon/melt for P, Sc, Ti, V, Y, La, Ce, Pr, Nd, Eu, Gd, Ho, Yb, Lu, Hf, Th, and U. In most cases Lu was the most (D 12-105) and La the least (0.06-0.95) compatible elements. D values from this study fall within the range of previously measured values for Rare Earth Elements (REE). However, D values measured experimentally show less fractionation than those recently measured using natural phenocryst/matrix pairs. For example, D_Lu/D_La measured experimentally in this study range between 27 and 206 compared to a value of 706,522 for a natural zircon/dacite pair [Sano, Y., Terada, K., and Fukuoka, T. 2002 High mass resolution ion microprobe analysis of rare earth elements in silicate glass, apatite and zircon: lack of matrix dependency. Chem. Geol.184, 217-230]. Although D values from this study show good agreement with the lattice strain model, D values from natural phenocryst/matrix pairs combined with measured zircon compositions better reproduce host-rock (magma) compositions of igneous rocks. They also yield more reasonable estimates of magma compositions when combined with compositions of ‘‘out-of-context” zircons. For example, compositions of the Hadean detrital zircons from Jack Hills, Australia yield LREE-enriched magmas when combined with D values from phenocryst/matrix pairs yields, but yield LREE-depleted magmas when experimentally determined D values are used. We infer that experimentally measured D^zircon/melt values represent disequilibrium partitioning resulting from rapid zircon growth during short laboratory timescales. Rapid growth causes development of observed diffusive boundary layers in the melt adjacent to zircon crystals. D values from phenocryst/matrix pairs are therefore recommended for petrogenetic modeling.     

Experimental study of polybaric REE partitioning between olivine, pyroxene and melt of the Yamato 980459 composition: Insights into the petrogenesis of depleted shergottites

A synthetic composition representing the Yamato 980459 martian basalt (shergottite) has been used to carry out phase relation, and rare earth element (REE) olivine and pyroxene partitioning experiments. Yamato 980459 is a sample of primitive basalt derived from a reduced end-member among martian mantle sources. Experiments carried out between 1-2 GPa and 1350-1650 °C simulate the estimated pressure-temperature conditions of basaltic melt generation in the martian mantle. Olivine-melt and orthopyroxene-melt partition coefficients for La, Nd, Sm, Eu, Gd and Yb (D_REE values) were determined by LA-ICPMS, and are similar to the published values for terrestrial basaltic systems. We have not detected significant variation in D-values with pressure over the range investigated, and by comparison with previous studies carried out at lower pressure.We apply the experimentally obtained olivine-melt and orthopyroxene-melt D_REE values to fractional crystallization and partial melting models to develop a three-stage geochemical model for the evolution of martian meteorites. In our model we propose two ancient (∼4.535 Ga) sources: the Nakhlite Source, located in the shallow mantle, and the Deep Mantle Source, located close to the martian core-mantle boundary. These two sources evolved distinctly on the ε¹⁴³Nd evolution curve due to their different Sm/Nd ratios. By partially melting the Nakhlite Source at ∼1.3 Ga, we are able to produce a slightly depleted residue (Nakhlite Residue). The Nakhlite Residue is left undisturbed until ∼500 Ma, at which point the depleted Deep Mantle Source is brought up by a plume mechanism carrying with it high heat flow, melts and isotopic signatures of the deep mantle (e.g., ε¹⁸²W, ε¹⁴²Nd, etc.). The plume-derived Deep Mantle Source combines with the Nakhlite Residue producing a mixture that becomes a mantle source (herein referred to as “the Y98 source”) for Yamato 980459 and the other depleted shergottites with the characteristic range of Sm/Nd ratios of these meteorites. The same hot plume provides a heat source for the formation of enriched and intermediate shergottites. Our model reproduces the REE patterns of nakhlites and depleted shergottites and can explain high ε¹⁴³Nd in depleted shergottites. Furthermore, the model results can be used to interpret whole rock Rb-Sr and Sm-Nd ages of shergottites.     

Partitioning of trace elements between rutile and silicate melts: Implications for subduction zones

Mineral/melt trace element partition coefficients were determined for rutile (TiO₂) for a large number of trace elements (Zr, Hf, Nb, Ta, V, Co, Cu, Zn, Sr, REE, Cr, Sb, W, U, Th). Whilst the high field strength elements (Zr, Hf, Nb, Ta) are compatible in rutile, other studied trace elements are incompatible (Sr, Th, REE). In all experiments we found D_Ta > D_Nb, D_Hf > D_Zr and D_U > D_Th. Partition coefficients for some polyvalent elements (Sb, W, and Co) were sensitive to oxygen fugacity. Melt composition exerts a strong influence on HFSE partition coefficients. With increasing polymerization of the melt, rutile/melt partition coefficients for the high field strength elements Zr, Hf, Nb and Ta increase about an order of magnitude. However, D_Nb/D_Ta and D_Hf/D_Zr are not significantly affected by melt composition. Because D_U ? D_Th, partial melting of rutile-bearing eclogite in subducted lithosphere may cause excesses of ²³⁰Th over ²³⁸U in some island arc lavas, whereas dehydration of subducted lithosphere may cause excesses of ²³⁸U over ²³⁰Th. From our partitioning results we infer partition coefficients for protactinium (Pa) which we predict to be much lower than previously anticipated. Contrary to previous studies, our data imply that rutile should not significantly influence observed ²³¹Pa-²³⁵U disequilibria in certain volcanic rocks.     

Trace element partitioning between apatite and silicate melts

We present new experimental apatite/melt trace element partition coefficients for a large number of trace elements (Cs, Rb, Ba, La, Ce, Pr, Sm, Gd, Lu, Y, Sr, Zr, Hf, Nb, Ta, U, Pb, and Th). The experiments were conducted at pressures of 1.0 GPa and temperatures of 1250 °C. The rare earth elements (La, Ce, Pr, Sm, Gd, and Lu), Y, and Sr are compatible in apatite, whereas the larger lithophile elements (Cs, Rb, and Ba) are strongly incompatible. Other trace elements such as U, Th, and Pb have partition coefficients close to unity. In all experiments we found D_Hf > D_Zr, D_Ta ≈ D_Nb, and D_Ba > D_Rb > D_Cs. The experiments reveal a strong influence of melt composition on REE partition coefficients. With increasing polymerisation of the melt, apatite/melt partition coefficients for the rare earth elements increase for about an order of magnitude. We also present some results in fluorine-rich and water-rich systems, respectively, but no significant influence of either H₂O or F on the partitioning was found. Furthermore, we also present experimentally determined partition coefficients in close-to natural compositions which should be directly applicable to magmatic processes.     

Distinct site preferences for heavy and light REE in amphibole and the prediction of Amph/LDREE

New experimental amphibole/melt partition coefficients from a variety of geologically relevant amphibole (pargasite, kaersutite, and K-richterite) and melt compositions obtained under conditions of interest to upper-mantle studies are combined with the results of X-ray single-crystal structure refinement. The ideal cation radii (r₀), calculated using the lattice-site elastic-strain model of Blundy and Wood (1994) under the hypothesis of complete REE (rare earth elements) ordering at ^[8]M4, mostly differ significantly from those obtained from both the structure refinement and the ionic radius of ^[8]Ca²⁺. Heavier REE may also strongly deviate from the parabolic trends defined by the other REE. On the basis of the crystal-chemical knowledge of major-element site-preference in amphibole and the occurrence of two sites with different co-ordination within the M4 cavity (M4 for Ca and Na, M4′ for Fe²⁺ and Mg), we propose a new model for REE incorporation. LREE order at the ^[8]M4 site, whereas HREE prefer the M4′ site with lower co-ordination in amphiboles with a significant cummingtonite component, and may also enter the M2 octahedron, at least in richterite. This more complex model is consistent with the observed ^Amph/LD, and drops the usual assumption that REE behave as a homogeneous group and order at the M4 site. The availability of multiple crystal-chemical mechanisms for REE³⁺ incorporation explains why measured and estimated ^Amph/LD_HREE may differ by up to one order of magnitude. When REE enter two different sites within the same cavity, a fit performed on the basis of a single curve may appear correct, but the values obtained for r₀ are biased towards those of the dominant site, and the Young's modulus is underestimated. When REE are incorporated in multiple sites in different cavities, the observed pattern cannot be reduced to a single curve, and the partition coefficients of heavy REE would be strongly underestimated by a single-site fit. The simplistic assumption that REE occupy a single site within the amphibole structure can thus substantially bias predictive models based on the elastic-strain theory. Our combined approach allows linkage between fine-scale site preference and the macroscopic properties of minerals and provides more reliable predictive models for mineral/melt partitioning. After the possible site-assignments have been identified, the shape of the Onuma curves constructed from accurately determined ^Amph/LD_REE now allows the active mechanisms for REE incorporation in amphiboles to be recognised even where site populations are not available. The REE preference for polyhedra with smaller size and lower co-ordination than those occupied by Ca invalidates the general idea that Ca acts as a “carrier” for REE. Received: 17 March 1999 / Accepted: 11 June 1999     

Experimental determination of F and Cl partitioning between lherzolite and basaltic melt

We experimentally determined F and Cl partition coefficients together with that of 19 trace elements (including REE, U-Th, HFSE and LILE) between basaltic melt and lherzolite minerals: olivine, orthopyroxene, clinopyroxene, plagioclase and garnet. Under conditions from 8 to 25 kbars and from 1,265 to 1,430°C, compatibilities of F and Cl are globally ordered as D ^Cpx/melt > D ^Opx/melt > D ^Grt/melt > D ^Ol/melt > D ^Plag/melt, and D _F ^mineral/melt is larger than D _Cl^mineral/melt. Four other major results were brought to light. (1) Chlorine partition coefficients positively correlate with the jadeite component in orthopyroxene, and increase of the CaTs component promotes Cl incorporation in clinopyroxene. (2) Variations of fluorine partition coefficients correlate strongly with melt viscosity. (3) F and Cl partition coefficients correlate with the Young’s modulus (E ₀) of pyroxene octahedral sites (M sites) and with Raman vibrational modes of pyroxenes. This demonstrates a fundamental interaction between the M site of pyroxenes and the incorporation of F and Cl. (4) We also determined the parameters of the lattice-strain model applied to 3+ cation trace elements for the two M sites in orthopyroxene and clinopyroxene: D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1 and E ₀^M2.     

Distribution of REE between clinopyroxene and basaltic melt along a mantle adiabat: effects of major element composition, water, and temperature

The distribution of rare earth elements (REE) between clinopyroxene (cpx) and basaltic melt is important in deciphering the processes of mantle melting. REE and Y partition coefficients from a given cpx-melt partitioning experiment can be quantitatively described by the lattice strain model. We analyzed published REE and Y partitioning data between cpx and basaltic melts using the nonlinear regression method and parameterized key partitioning parameters in the lattice strain model (D ₀, r ₀ and E) as functions of pressure, temperature, and compositions of cpx and melt. D ₀ is found to positively correlate with Al in tetrahedral site (Al ^T ) and Mg in the M2 site (Mg^M2) of cpx and negatively correlate with temperature and water content in the melt. r ₀ is negatively correlated with Al in M1 site (Al^M1) and Mg^M2 in cpx. And E is positively correlated with r ₀. During adiabatic melting of spinel lherzolite, temperature, Al ^T , and Mg^M2 in cpx all decrease systematically as a function of pressure or degree of melting. The competing effects between temperature and cpx composition result in very small variations in REE partition coefficients along a mantle adiabat. A higher potential temperature (1,400°C) gives rise to REE partition coefficients slightly lower than those at a lower potential temperature (1,300°C) because the temperature effect overwhelms the compositional effect. A set of constant REE partition coefficients therefore may be used to accurately model REE fractionation during partial melting of spinel lherzolite along a mantle adiabat. As cpx has low Al and Mg abundances at high temperature during melting in the garnet stability field, REE are more incompatible in cpx. Heavy REE depletion in the melt may imply deep melting of a hydrous garnet lherzolite. Water-dependent cpx partition coefficients need to be considered for modeling low-degree hydrous melting.     

Partitioning behavior of trace elements between dacitic melt and plagioclase, orthopyroxene, and clinopyroxene based on laser ablation ICPMS analysis of silicate melt inclusions

Partitioning behavior of Sc, Ti, V, Mn, Sr, Y, Zr, Nb, Ba, La, Nd, Sm, Eu, Gd, Dy, Ho, Yb, Hf, and Pb between dacitic silicate melt and clinopyroxene, orthopyroxene, and plagioclase has been determined based on laser ablation-inductively coupled plasma mass spectrometric (LA-ICPMS) analysis of melt inclusions and the immediately adjacent host mineral. Samples from the 1988 eruption of White Island, New Zealand were selected because petrographic evidence suggests that all three mineral phases are in equilibrium with each other and with the melt inclusions. All three phenocryst types are found as mineral inclusions within each of the other phases, and mineral inclusions often coexist with melt inclusions in growth-zone assemblages. Compositions of melt inclusions do not vary between the different host minerals, suggesting that boundary layer processes did not affect compositions of melt inclusions and that post-trapping modifications have not occurred.Partition coefficients were calculated from the host and melt inclusion compositions and results were compared to published values. All trace elements examined in this study except Sr are incompatible in plagioclase, and all measured trace elements except for Mn are incompatible in orthopyroxene. In clinopyroxene, Sc, V, and Mn are compatible, and Y, Ti, HREE, and the MREE are only slightly incompatible. Most partition coefficients overlap the wide range of values reported in the literature, but the White Island data are consistently at the lower end of the range in published values. Results from the literature obtained using modern microanalytical techniques such as secondary ion mass spectrometry (SIMS) or proton induced X-ray emission spectroscopy (PIXE) also fall at the lower end of the published values, whereas partition coefficients determined from bulk analysis of glass and crystals separated from volcanic rocks typically extend to higher values. Rapid crystal growth-rates, crystal zonation, or the presence of accessory mineral inclusions in phenocrysts likely accounts for the wide range and generally higher partition coefficients obtained using bulk sampling techniques. The results for 3+ cations from this study are consistent with theoretical predictions based on a lattice strain model for site occupancy. The results also confirm that the melt inclusion-mineral (MIM) technique is a reliable method for determining partition coefficients, as long as the melt inclusions have not experienced post-entrapment reequilibration.     

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

S-rich apatite-hosted glass inclusions in xenoliths from La Palma: constraints on the volatile partitioning in evolved alkaline magmas

The composition of S-rich apatite, of volatile-rich glass inclusions in apatite, and of interstitial glasses in alkaline xenoliths from the 1949 basanite eruption in La Palma has been investigated to constrain the partitioning of volatiles between apatite and alkali-rich melts. The xenoliths are interpreted as cumulates from alkaline La Palma magmas. Apatite contains up to 0.89 wt% SO₃ (3560 ppm S), 0.31 wt% Cl, and 0.66 wt% Ce₂O₃. Sulfur is incorporated in apatite via several independent exchange reactions involving (P⁵⁺, Ca²⁺) vs. (S⁶⁺, Si⁴⁺, Na⁺, and Ce³⁺). The concentration of halogens in phonolitic to trachytic glasses ranges from 0.15 to 0.44 wt% for Cl and from <0.07 to 0.65 wt% for F. The sulfur concentration in the glasses ranges from 0.06 to 0.23 wt% SO₃ (sulfate-saturated systems). The chlorine partition coefficients (D_Cl^{apatite/glass}) range from 0.4 to 1.3 (average D_Cl^{apatite/glass} = 0.8), in good agreement with the results of experimental data in mafic and rhyolitic system with low Cl concentrations. With increasing F in glass inclusions D_F^{apatite/glass} decreases from 35 to 3. However, most of our data display a high partition coefficient (~30) close to D_F^{apatite/glass} determined experimentally in felsic rock. D_S^{apatite/glass} decreases from 9.1 to 2.9 with increasing SO₃ in glass inclusions. The combination of natural and experimental data reveals that the S partition coefficient tends toward a value of 2 for high S content in the glass (>0.2 wt% SO₃). D_S^{apatite/glass} is only slightly dependent on the melt composition and can be expressed as: SO_{3 apatite} (wt%) = 0.157 * ln SO_{3 glass} (wt%) + 0.9834. The phonolitic compositions of glass inclusions in amphibole and haüyne are very similar to evolved melts erupted on La Palma. The lower sulfur content and the higher Cl content in the phonolitic melt compared to basaltic magmas erupted in La Palma suggest that during magma evolution the crystallization of haüyne and pyrrhotite probably buffered the sulfur content of the melt, whereas the evolution of Cl concentration reflects an incompatible behavior. Trachytic compositions similar to those of the (water-rich) glass inclusions analyzed in apatite and clinopyroxene are not found as erupted products. These compositions are interpreted to be formed by the reaction between water-rich phonolitic melt and peridotite wall-rock.     

An experimental investigation of the partitioning of REE between garnet and liquid with reference to the role of defect equilibria

The partitioning of samarium and thulium between garnets and melts in the systems Mg₃Al₂-Si₃O₁₂-H₂O and Ca₃Al₂Si₃O₁₂-H₂O has been studied as a function of REE concentration in the garnets at 30 kbar pressure. Synthesis experiments of variable time under constant P, T conditions indicate that garnet initially crystallizes rapidly to produce apparent values of D _Sm (D _Sm=concentration of Sm in garnet/concentration of Sm in liquid) which are too large in the case of pyrope and too small in the case of grossular. As the experiment proceeds, Sm diffuses out of or into the garnet and the equilibrium value of D _Sm is approached. Approximate values of diffusion coefficients for Sm in pyrope garnet obtained by this method are 6 × 10^–13 cm² s^–1 at 1,300 ° C and 2 × 10^–12 cm² s^–1 at 1,500 ° C, and for grossular, 8.3 × 10^–12 cm² s^–1 at 1,200 ° C and 4.6 × 10^–11 cm² s^–1 at 1,300 ° C. The equilibrium values of D _Sm have been reversed by experiments with Sm-free pyrope and Sm-bearing glass, and with Sm-bearing grossular and Sm-free glass.Between 12 ppm and 1,000 ppm Sm in pyrope at 1,300 ° C and between 80 ppm and >2 wt.% Tm in pyrope at 1,500 ° C, partition coefficients are constant and independent of REE concentration. Above 100 ppm of Sm in garnet at 1,500 ° C, partition coefficients are independent of Sm concentration. At lower concentrations, however, D _Sm is dependent upon the Sm content of the garnet. The two regions may be interpreted in terms of charge-balanced substitution of Sm₃Al₅O₁₂ in the garnet at high Sm concentrations and defect equilibria involving cation vacancies at low concentrations. At very low REE concentrations (< 1 ppm Tm in grossular at 1,300 ° C) D_REE garnet/liquid again becomes constant with an apparent Henry's Law value greater than that at high concentrations. This may be interpreted in terms of a large abundance of cation vacancies relative to the number of REE ions.The importance of defects in the low concentration region has been confirmed by adding other REE (at 80 ppm level) to the system Mg₃Al₂Si₃O₁₂-H₂O at low Sm concentrations. These change D _Sm in the defect region, demonstrating their role in the production of vacancies.Experiments on a natural pyropic garnet indicate that defect equilibria are of importance to REE partitioning within the concentration ranges found in nature.     

Divergent T–ƒO<Subscript>2</Subscript> paths during crystallisation of H<Subscript>2</Subscript>O-rich and H<Subscript>2</Subscript>O-poor magmas as recorded by Ce and U in zircon,with implications for TitaniQ and TitaniZ geothermometry

During solidification of magma chambers as systems closed to chemical exchange with environs, the residual siliceous melt may follow a trend of rising, constant, or decreasing oxidation state, relative to reference buffers such as nickel?+?nickel oxide (NNO) or fayalite?+?magnetite?+?quartz. Titanomagnetite–hemoilmenite thermometry and oxybarometry on quenched volcanic suites yield temperature versus oxygen fugacity arrays of varied positive and negative slopes, the validity of which has been disputed for several years. We resolve the controversy by introducing a new recorder of magmatic redox evolution employing temperature- and redox-sensitive trace-element abundances in zircon. The zircon/melt partition coefficients of cerium and uranium vary oppositely in response to variation of magma redox state, but vary in tandem as temperature varies. Plots of U/Pr versus Ce⁴⁺/Ce³⁺ in zircon provide a robust test for change in oxidation state of the melt during zircon crystallisation from cooling magma, and the plots discriminate thermally induced from redox-induced variation of Ce⁴⁺/Ce³⁺ in zircon. Temperature-dependent lattice strain causes Ce⁴⁺/Ce³⁺ in zircon to increase strongly as zircon crystallises from cooling magma at constant Ce⁴⁺/Ce³⁺ ratio in the melt. We examine 19 zircon populations from igneous complexes in varied tectonic settings. Variation of zircon Ce⁴⁺/Ce³⁺ due to minor variation in melt oxidation state during crystallisation is resolvable in 11 cases but very subordinate to temperature dependence. In many zircon populations described in published literature, there is no resolvable change in redox state of the melt during tenfold variation of Ce⁴⁺/Ce³⁺ in zircons. Varied magmatic redox trends indicated by different slopes on plots of zircon U/Pr versus Ce⁴⁺/Ce³⁺ are corroborated by Fe–Ti-oxide-based T–?O₂ trends of correspondingly varied slopes. Zircon and Fe–Ti-oxide compositions agree that exceptionally, H₂O-rich arc magmas tend to follow a trend of rising oxidation state of the melt during late stages of fluid-saturated magmatic differentiation at upper-crustal pressures. We suggest that H₂ and/or SO₃ and/or Fe²⁺ loss from the melt to segregating fluid is largely responsible. Conversely, zircon and Fe–Ti-oxide compositions agree in indicating that H₂O-poor magmas tend to follow a T–?O₂ trend of decreasing oxidation state of the melt during late stages of magmatic differentiation at upper-crustal pressures, because the precipitating mineral assemblage has higher Fe³⁺/Fe²⁺ than coexisting rhyolitic melt. We present new evidence showing that the Fe–Ti-oxide oxybarometer calibration by Ghiorso and Evans (Am J Sci 308(9):957–1039, 2008) retrieves experimentally imposed values of ?O₂ in laboratory syntheses of Fe–Ti-oxide pairs to a precision of ±?0.2 log unit, over a large experimental temperature range, without systematic bias up to at least log ?O₂?≈?NNO?+?4.4. Their titanomagnetite–hemoilmenite geothermometer calibration has large systematic errors in application to Ti-poor oxides that precipitate from very oxidised magmas. A key outcome is validation of Fe–Ti-oxide-based values of melt TiO₂ activity for use in Ti-in-zircon thermometry and Ti-in-quartz thermobarometry.     

Some experimentally determined zircon/liquid partition coefficients for the rare earth elements

Partition coefficients for La, Sm, Ho and Lu (REE) between synthetic zircon and felsic, peralkaline liquid were determined at 800°C and 2 kbar water pressure by adding small amounts of REE to experimental charges and analyzing zircons in the quenched runs with an electron microprobe. The preferred zircon/liquid partition coefficients obtained by this method are: La, 1.4–2.1; Sm, 26–40; Ho, 340+; Lu, 72–126. These results confirm the strong heavy rare earth enrichment discovered by Nagasawa (1970) in zircon separates from dacites and granites, but they also introduce a modification to the supposed zircon/liquid partition coefficient pattern. The heavy REE end of the pattern is concave downward, in qualitative resemblance to some garnet/liquid and hornblende/liquid REE partitioning patterns.     

Kiglapait geochemistry VII: Yttrium and the rare earth elements

Based on 51 wholerock analyses by XRF and summation over the layered group, the Kiglapait Intrusion contains 4.7_?1.6^+1.2 ppm Y, which resides principally in augite and apatite. Using liquid compositions calculated by summation, the partition coefficient D^AUG/L_Y is 0.95 ± 0.12 from 84 to 97 PCS (percent solidified) and 1.5 ± 0.4 above 97 PCS. For feldspar, the most likely value for D is 0.028 ± 0.02 (N = 6).REE analyses for 13 whole rocks were interpreted with the aid of yttrium models to yield trends for wholerocks and liquids vs PCS. Summations over the rocks of the layered group gave La = 2.5, Ce = 5.8, Nd = 3.9, Sm = 1.0, Eu = 0.8, Tb = 0.17, Yb = 0.37, and Lu = 0.06 ppm, with 2 s.d. errors near ± 30%. All these elements are highly incompatible until the arrival of augite, which affects chiefly the HREE, and apatite, which affects all (but more strongly, the LREE). The net result is that after apatite arrival at 94 PCS, the liquid compositions are nearly constant, hence D^WR/L_REE ≈ 1.0. These results are compatible with the mineralogy of the intrusion and the estimated partition coefficients for feldspar, olivine, augite, apatite, and Fe-Ti oxide minerals. For pre-apatite liquids, D^FSP/L_REE vary regularly with the normative di content of the liquid and change by an order of magnitude, hence the bulk liquid composition must be considered in any attempt to invert the compositions of feldspars to parent liquids.The Eu anomaly at first decreases in Kiglapait liquids due to plagioclase fractionation, but then increases due to removal of augite and apatite with negative Eu anomalies. The features dominantly responsible for Eu partitioning are liquid structure and, for monoclinic ternary feldspars, crystal structure. The former is best monitored by the augite or diopside content of the liquid and the latter, by the K content of the feldspar.The chondrite-normalized REE pattern for the intrusion has La_N = 7.4, Lu_N = 1.6, (Ce/Yb)_N = 3.6, and Eu/Eu* = 2.4, indicating its feldspar-rich nature. The chilled margin of the nearby Hettasch Intrusion has a similar but more evolved pattern, corresponding roughly to the Kiglapait liquid at 70 PCS. As with other data, those for the REE suggest source differences for the two intrusions rather than a relationship due to fractionation.     

Magmatic evolution of the differentiated ultramafic, alkaline and carbonatite intrusion of Vuoriyarvi (Kola Peninsula, Russia). A LA-ICP-MS study of apatite

The nature of the petrogenetic links between carbonatites and associated silicate rocks is still under discussion (i.e., [Gittins J., Harmer R.E., 2003. Myth and reality of the carbonatite–silicate rock “association”. Period di Mineral. 72, 19–26.]). In the Paleozoic Kola alkaline province (NW Russia), the carbonatites are spatially and temporally associated to ultramafic cumulates (clinopyroxenite, wehrlite and dunite) and alkaline silicate rocks of the ijolite–melteigite series [(Kogarko, 1987), (Kogarko et al., 1995), (Verhulst et al., 2000), (Dunworth and Bell, 2001) and (Woolley, 2003)]. In the small (≈ 20 km²) Vuoriyarvi massif, apatite is typically a liquidus phase during the magmatic evolution and so it can be used to test genetic relationships. Trace elements contents have been obtained for both whole rocks and apatite (by LA-ICP-MS). The apatites define a single continuous chemical evolution marked by an increase in REE and Na (belovite-type of substitution, i.e., 2Ca²⁺ = Na⁺ + REE³⁺). This evolution possibly reflects a fractional crystallisation process of a single batch of isotopically homogeneous, mantle-derived magma.The distribution of REE between apatite and their host carbonatite have been estimated from the apatite composition of a carbonatite vein, belonging to the Neskevara conical-ring-like vein system. This carbonatite vein is tentatively interpreted as a melt. So, the calculated distribution coefficients are close to partition coefficients. Rare earth elements are compatible in apatite (D > 1) with a higher compatibility for the middle REE (D_Sm : 6.1) than for the light (D_La : 4.1) and the heavy (D_Yb : 1) REE.