Determination of zircon/melt trace element partition coefficients from SIMS analysis of melt inclusions in zircon

Partition coefficients (^zircon/meltD_M) for rare earth elements (REE) (La, Ce, Nd, Sm, Dy, Er and Yb) and other trace elements (Ba, Rb, B, Sr, Ti, Y and Nb) between zircon and melt have been calculated from secondary ion mass spectrometric (SIMS) analyses of zircon/melt inclusion pairs. The melt inclusion-mineral (MIM) technique shows that D_REE increase in compatibility with increasing atomic number, similar to results of previous studies. However, D_REE determined using the MIM technique are, in general, lower than previously reported values. Calculated D_REE indicate that light REE with atomic numbers less than Sm are incompatible in zircon and become more incompatible with decreasing atomic number. This behavior is in contrast to most previously published results which indicate D > 1 and define a flat partitioning pattern for elements from La through Sm. The partition coefficients for the heavy REE determined using the MIM technique are lower than previously published results by factors of ≈15 to 20 but follow a similar trend. These differences are thought to reflect the effects of mineral and/or glass contaminants in samples from earlier studies which employed bulk analysis techniques.D_REE determined using the MIM technique agree well with values predicted using the equations of Brice (1975), which are based on the size and elasticity of crystallographic sites. The presence of Ce⁴⁺ in the melt results in elevated D_Ce compared to neighboring REE due to the similar valence and size of Ce⁴⁺ and Zr⁴⁺. Predicted ^zircon/meltD values for Ce⁴⁺ and Ce³⁺ indicate that the Ce⁴⁺/Ce³⁺ ratios of the melt ranged from about 10⁻³ to 10⁻². Partition coefficients for other trace elements determined in this study increase in compatibility in the order Ba < Rb < B < Sr < Ti < Y < Nb, with Ba, Rb, B and Sr showing incompatible behavior (D_M < 1.0), and Ti, Y and Nb showing compatible behavior (D_M > 1.0).The effect of partition coefficients on melt evolution during petrogenetic modeling was examined using partition coefficients determined in this study and compared to trends obtained using published partition coefficients. The lower D_REE determined in this study result in smaller REE bulk distribution coefficients, for a given mineral assemblage, compared to those calculated using previously reported values. As an example, fractional crystallization of an assemblage composed of 35% hornblende, 64.5% plagioclase and 0.5% zircon produces a melt that becomes increasingly more enriched in Yb using the D_Yb from this study. Using D_Yb from Fujimaki (1986) results in a melt that becomes progressively depleted in Yb during crystallization.     

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.     

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.     

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 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 investigation of large-ion-lithophile-element-, high-field-strength-element- and rare-earth-element-partitioning between calcic amphibole and basaltic melt: the effects of pressure and oxygen fugacity

The effects of pressure and oxygen fugacity (fO₂) on trace element partitioning between pargasitic amphibole and alkali-basaltic melts have been determined at pressures from 1.5 to 2.5 GPa and oxygen fugacities at 2 log units above and below the nickel–nickel oxide buffer. Amphibole crystallization experiments were performed in a piston cylinder apparatus and partition coefficients between amphibole and quenched melt of large-ion-lithophile elements (LILE: Rb, Sr, Ba), high-field-strength elements (HFSE: Zr, Nb, Ta, Hf, U, Th) and rare-earth elements (REE: La to Lu; +Y) were measured with a LASER ablation inductively coupled plasma – mass spectrometer. Increasing pressure from 1.5 to 2.5 GPa at similar temperatures and approximately constant fO₂ increases D _Rb but decreases D _Zr and D _Hf and D _REE (D _La, D _Ce, D _Pr). An empirical relationship was observed between D _Zr and (Ti/Al)_M2 in the amphibole, which can be described by:

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     

Trace element partitioning in plagioclase feldspar

Compilation and interpretation of experimental and natural Nernst partition coefficient (^{plagioclase/melt}D) data show that, with a few exceptions, increases in ^{plagioclase/melt}D correlate with decreasing anorthite-content of plagioclase. In contrast, increases of ^{plagioclase/melt}D for Ga, Sc, Cu, Zn, Zr, Hf and Ti, are better correlated against decreasing melt MgO or increasing melt SiO₂ contents. ^{plagioclase/melt}D for Ti and the rare earth elements (REE) show little dependence on temperature, but increase as the melt water content increases. ^{plagioclase/melt}D for K and Sr are sensitive to pressure. Variations of D₀ (the strain compensated partition coefficient), r₀ (the size of the site into which REE substitute), and E (Young’s Modulus of this site) were parameterized against variations of melt SiO₂, the An-content of plagioclase, and other combinations of variables, allowing ^{plagioclase/melt}D_REE-Y to be calculated from a variety of input parameters. The interrelations of temperature, melt MgO and SiO₂ content, and plagioclase anorthite-content for wet and dry systems were also parameterized to facilitate interpolation where such data are lacking. When combined, these semi-empirical parameterizations yield ^{plagioclase/melt}D results comparable to available experimental and natural data.     

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.     

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.     

The influence of melt composition on the partitioning of REEs, Y, Sc, Zr and Al between forsterite and melt in the system CMAS

Partition coefficients for a range of Rare Earth Elements (REEs), Y, Sc, Al and Zr were determined between forsteritic olivine (nearly end-member Mg₂SiO₄) and ten melt compositions in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS) at 1 bar and 1400 °C, with concentrations of the trace elements in the olivine and the melt measured by laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). The REEs and Sc were added at levels sufficient to ensure that concentrations in the olivine were well above the detection limits. The REE partition coefficients decrease with increasing silica in the melt, indicating strong bonding between REEO_1.5 and SiO₂ in the melt. The variation of as a function of ionic radius is well described by the Brice equation for each composition, although a small proportion of this variation is due to the increase in the strength of the REEO_1.5-SiO₂ interactions in the melt with ionic radius. Scandium behaves very similarly to the REEs, but a global fit of the data from all ten melt compositions suggests that deviates somewhat from the parabolas established by the REE and Y, implying that Sc may substitute into olivine differently to that of the REEs. In contrast to the behaviour of the large trivalent cations, the concentration of Al in olivine is proportional to the square root of its concentration in the melt, indicating a coupled substitution in olivine with a high degree of short-range order. The lack of any correlation of REE partition coefficients with Al in olivine or melt suggests that the REE substitution in olivine is charge-balanced by cation vacancies. The partition coefficient of the tetravalent trace element Zr, which is highly incompatible in olivine, depends on the CaO content of the melt.     

Rare earth element partitioning between clinopyroxene and silicate liquid at moderate to high pressure

Experimental determination of over seventy sets of clinopyroxene/silicate liquid (glass) partitition coefficients (D) for four rare earth elements (REE — La, Sm, Ho, Lu) in a range of REE-enriched natural rock compositions (basalt, basaltic andesite, andesite and rhyodacite) demonstrate a convex upward pattern, favouring the heavy REE (Ho, Lu) and markedly discriminating against the light REE (La). These patterns are consistent with previously documented clinopyroxene D values reported from natural phenocryst/matrix pairs and from experimental work using either REE-enriched compositions and electron microprobe analytical techniques (as in the present study) or natural or synthetic undoped compositions and mass spectrometric, ion probe or X-ray autoradiographic analytical techniques. However, the large data base in the present study allows evaluation of the effect of compositional and physical parameters on REE partitioning relationships. Considering D_Ho, it is shown that (1) D increases 6-fold with increasing SiO₂ content of the coexisting liquid from 50 to 70 wt% SiO₂ (2) D increases 4-fold with decreasing temperature from 1,120°C to 900° C (3) D increases 2-fold with increasing pressure from 2.5 to 20 kb. (4) D increases 2-fold fO₂ increases from approximately that of the MW buffer to the HM buffer (5) D remains unchanged within experimental error as the water content of the melt changes from 0.3 to 10% by weight H₂O.The absolute REE content of the clinopyroxene shows no consistent trend with temperature, but decreases slightly with increasing pressure, paralleling an increase in the jadeite component of the pyroxene. Thus the increase in D with increasing pressure is attributed to changes in the silicate liquid structure, which discriminate against accommodation of REE with increasing pressure. The clinopyroxene REE content increases with increasing fO₂, and in this case the increase in D with increasing fO₂ may be attributed mainly to this change in the clinopyroxene composition. Application of the present results to geochemical modelling allows a more appropriate choice of D values, according to the liquid composition and physical conditions applicable in the modelled system. They may also be used to evaluate cognate or xenocrystic relationships between clinopyroxene megacrysts and their host matrix.     

Partial melting of secondary pyroxenite at 1 and 1.5 GPa,and its role in upwelling heterogeneous mantle

We performed partial melting experiments at 1 and 1.5 GPa, and 1180–1400 °C, to investigate the melting under mantle conditions of an olivine-websterite (GV10), which represents a natural proxy of secondary (or stage 2) pyroxenite. Its subsolidus mineralogy consists of clinopyroxene, orthopyroxene, olivine and spinel (+garnet at 1.5 GPa). Solidus temperature is located between 1180 and 1200 °C at 1 GPa, and between 1230 and 1250 °C at 1.5 GPa. Orthopyroxene (±garnet), spinel and clinopyroxene are progressively consumed by melting reactions to produce olivine and melt. High coefficient of orthopyroxene in the melting reaction results in relatively high SiO₂ content of low melt fractions. After orthopyroxene exhaustion, melt composition is controlled by the composition of coexisting clinopyroxene. At increasing melt fraction, CaO content of melt increases, whereas Na₂O, Al₂O₃ and TiO₂ behave as incompatible elements. Low Na₂O contents reflect high partition coefficient of Na between clinopyroxene and melt (\(D_{{{\text{Na}}_{ 2} {\text{O}}}}^{{{\text{cpx}}/{\text{liquid}}}}\)). Melting of GV10 produces Quartz- to Hyperstene-normative basaltic melts that differ from peridotitic melts only in terms of lower Na₂O and higher CaO contents. We model the partial melting of mantle sources made of different mixing of secondary pyroxenite and fertile lherzolite in the context of adiabatic oceanic mantle upwelling. At low potential temperatures (T _P < 1310 °C), low-degree melt fractions from secondary pyroxenite react with surrounding peridotite producing orthopyroxene-rich reaction zones (or refertilized peridotite) and refractory clinopyroxene-rich residues. At higher T _P (1310–1430 °C), simultaneous melting of pyroxenite and peridotite produces mixed melts with major element compositions matching those of primitive MORBs. This reinforces the notion that secondary pyroxenite may be potential hidden components in MORB mantle source.     

Clinopyroxene/liquid trace element partitioning in natural trachyte–trachyphonolite systems: insights from Campi Flegrei (southern Italy)

Trace element partition coefficients between clinopyroxenes and associated glassy matrix (^Cpx/L D) have been determined for 13 REE, HFSE^4+,5+, U, Th, Sr, Pb, Sc and V from combined LA-ICP-MS/EMP analyses in selected trachytes and trachyphonolites from Campi Flegrei. Composition of clinopyroxene and glass is pretty homogeneous in the trachyphonolites, pointing to an overall attainment of the equilibrium conditions. In trachytes, conversely, phases show some compositional heterogeneity (due to the presence of clinopyroxene xenocrysts) that requested a more careful petrographic and geochemical inspection of the samples to assess the equilibrium clinopyroxene composition. In the trachyte clinopyroxenes, REE are compatible from Nd to Lu (^Cpx/L D up to 2.9), like Y, Ti, Sc and V. The ^Cpx/L D for Eu is lower than those of the adjacent REE, highlighting Eu²⁺ contribution. High D values are also shown by U, Th, Pb, Zr, Hf, Nb and Ta relatively to basaltic and andesitic systems, whereas the D _Sr is roughly similar to that found for less evolved magmas. Trachyphonolites are characterized by an overall decrease of the ^Cpx/L D for highly-charged cations (with the exception of V), and by a slight increase of D _Sr. REE are still compatible from Nd to Lu (^Cpx/L D up to 2.1), like Ti, Y, Sc and V. This variation is also predicted for REE and Y by models based on the elastic strain theory, being consistent with the slightly lower polymerization degree estimated for the trachyphonolites. However, the observed ^Cpx/L D _(REE,Y) are matched by the modelled ones only considering very low T (≤825°C), which are believed unlikely. This mismatch cannot be attributed to effects induced by the water-rich composition of the trachyte–trachyphonolite suite, since they would lower the observed ^Cpx/L D _(REE,Y). Moreover, the anomalous inflections of measured ^Cpx/L D for HREE suggests some crystal-chemical control, such as the entrance of these elements in a site distinct from M2. It is concluded that the large ^Cpx/L D determined for trachytes and trachyphonolites are likely induced by hitherto unconstrained changes of the Z³⁺ activities related to the composition of melt and/or solid. All these considerations strongly highlight the importance of a direct characterization of trace element partitioning in natural samples from magmatic systems poorly characterized by experimental studies. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Rare earth elements in phoscorites and carbonatites of the Devonian Kola Alkaline Province,Russia: Examples from Kovdor,Khibina, Vuoriyarvi and Turiy Mys complexes

The Devonian (ca. 385–360 Ma) Kola Alkaline Province includes 22 plutonic ultrabasic–alkaline complexes, some of which also contain carbonatites and rarely phoscorites. The latter are composite silicate–oxide–phosphate–carbonate rocks, occurring in close space-time genetic relations with various carbonatites. Several carbonatites types are recognized at Kola, including abundant calcite carbonatites (early- and late-stage), with subordinate amounts of late-stage dolomite carbonatites, and rarely magnesite, siderite and rhodochrosite carbonatites. In phoscorites and early-stage carbonatites the rare earth elements (REE) are distributed among the major minerals including calcite (up to 490 ppm), apatite (up to 4400 ppm in Kovdor and 3.5 wt.% REE₂O₃ in Khibina), and dolomite (up to 77 ppm), as well as accessory pyrochlore (up to 9.1 wt.% REE₂O₃) and zirconolite (up to 17.8 wt.% REE₂O₃). Late-stage carbonatites, at some localities, are strongly enriched in REE (up to 5.2 wt.% REE₂O₃ in Khibina) and the REE are major components in diverse major and minor minerals such as burbankite, carbocernaite, Ca- and Ba-fluocarbonates, ancylite and others. The rare earth minerals form two distinct mineral assemblages: primary (crystallized from a melt or carbohydrothermal fluid) and secondary (formed during metasomatic replacement). Stable (C–O) and radiogenic (Sr–Nd) isotopes data indicate that the REE minerals and their host calcite and/or dolomite have crystallized from a melt derived from the same mantle source and are co-genetic.     

Kimberlite petrogenesis: Insights from clinopyroxene-melt partitioning experiments at 6 GPa in the CaO-MgO-Al2O3-SiO2-CO2 system

In this experimental study, we examine the mineral-melt partitioning of major and trace elements between clinopyroxene and CO₂-rich kimberlitic melts at a pressure of 6 GPa and temperatures of 1410°C and 1430°C. The melts produced contain ∼ 28 wt% dissolved CO₂, and are saturated with olivine and clinopyroxene. To assess the effects of temperature, crystal and melt compositions on trace element partitioning, experiments were performed in the model CaO-MgO-Al₂O₃-SiO₂-CO₂ system. Our results reveal that all the elements studied, except Al, Mg, Si, and Ga, are incompatible in clinopyroxene. Partition coefficients show a considerable range in magnitude, from ∼ 10⁻³ for D_U and D_Ba to ∼ 2.5 for D_Si. The two experimental runs show similar overall partitioning patterns with the D values being lower at 1430°C. Rare earth elements display a wide range of partition coefficients, D_La (0.012-0.026) being approximately one order of magnitude lower than D_Lu (0.18-0.23). Partition coefficients for the 2+ and 3+ cations entering the M2-site exhibit a near-parabolic dependence on radius of the incorporated cations as predicted from the lattice strain model. This underlines the contribution made by the crystal structure toward controlling the distribution of trace elements. Using data obtained in this study combined with that in the published literature, we also discuss the effects that other important parameters, namely, melt composition, pressure, and temperature, could have on partitioning.Our partition coefficients have been used to model the generation of the Group I (GI) kimberlites from South Africa. The numerical modeling shows that kimberlitic melts can be produced by ∼0.5% melting of a MORB-type depleted source that has been enriched by small-degree melts originating from a similar depleted source. This result suggests that the source of GI kimberlites may be located at the lithosphere-asthenosphere transition. Percolation of small degree melts from the asthenosphere would essentially create a metasomatic horizon near the bottom of the non-convecting sublithospheric mantle. Accumulation of such small degree melts together with the presence of volatiles and conductive heating would trigger melting of the ambient mantle and subsequently lead to eruption of kimberlitic melts. Additionally, our model shows that the GI source can be generated by metasomatism of a 2 Ga old MORB source ca. 1 Ga ago. Assuming that MORB-type mantle is the most depleted source of magmas on earth, then this is the oldest age at which the GI source could have existed. However, this age most likely reflects the average age of a series of metasomatic events than that of a single event.