Isovalent trace element partitioning between minerals and melts: A computer simulation study

We present a new approach for the rationalisation of trace element partitioning between silicate melts and minerals, which is not based on the empirical, parameterised continuum models in common use. We calculate the energetics of ion substitution using atomistic simulation techniques, which include an explicit evaluation of the relaxation energy (strain energy) contribution to this process. Solution energies are estimated for isovalent impurities in CaO, diopside, orthoenstatite, and forsterite. These show a parabolic dependence on ionic radius, similar to the variation of mineral-melt partition coefficients with ionic radius. The success of the empirical models, which often include only the strain energy, appear to have been due to the partial cancellation of energy terms, and to the empirical fitting of the parameters included in these models. Our approach can be readily extended to aliovalent substitution.     

Laser Ablation ICPMS study of trace element partitioning between plagioclase and basaltic melts: an experimental approach

Plagioclase-melt partition coefficients (D) for 34 trace elements at natural concentration levels were determined experimentally in a natural MORB composition at atmospheric pressure using thin Pt-wire loops. Experiments were carried out at three temperatures (1,220, 1,200, and 1,180°C), and at three different oxygen fugacities (fO₂ = IW, QFM, air) in order to assess the effect of fO₂ on the partitioning of elements with multiple valence (Fe, Eu, Cr). Run products were analyzed by laser-ablation ICP-MS. Most trace element Ds increase slightly as temperature decreases, except for D _Zr, D _Fe, D _Eu and D _Cr that vary systematically with fO₂. Applying the Lattice Strain Model to our data suggests the presence of Fe²⁺ entirely in the octahedral site at highly to moderate reducing conditions, while Fe³⁺ was assigned wholly to the tetrahedral site of the plagioclase structure. Furthermore, we provide a new quantitative framework for understanding the partitioning behaviour of Eu, which occurs as both 2+ and 3+ cations, depending on fO₂and confirm the greater compatibility of Eu²⁺, which has an ionic radius similar to Sr, relative to Eu³⁺ in plagioclase and the higher Eu²⁺/ Eu³⁺ under reducing conditions. For petrogenetic basaltic processes, a combined fractionation of Eu²⁺–Sr and Fe–Mg by plagioclase has considerable potential as an oxybarometer for natural magmatic rocks.     

Rare earth element partitioning between minerals from anhydrous spinel peridotite xenoliths

REE abundances in minerals from spinel peridotite xenoliths from West Germany, the south-western U.S. and Mongolia decrease in the order clinopyroxene > orthopyroxene > olivine > spinel. While clinopyroxenes are similar in absolute chondrite-normalized concentrations to those known from other studies, orthopyroxenes and olivines are significantly lower in LREE although comparable in HREE. Spinels are much lower in all REE than any previously reported values and are completely negligible for the REE budget of peridotites.Partition coefficients for most orthopyroxene/clinopyroxene pairs increase systematically from La to Lu. Olivine/clinopyroxene and spinel/clinopyroxene partition coefficients increase from the intermediate rare earth elements to Lu and normally are higher for La compared to Sm.The application of Nagasawa's (1966) elastic lattice model suggests that all heavy but only minor amounts of the light REE substitute into structural positions of orthopyroxene and olivine.Significant differences between orthopyroxene/clinopyroxene partition coefficients for various xenoliths may be assigned to dependences upon equilibration temperature and bulk chemistry.Apart from grain surface contaminations, fluid inclusions which are practically always present in mantle minerals, can highly concentrate light rare earth elements and thus may be responsible for unexpectedly high concentrations of incompatible elements frequently reported for mantle olivines or orthopyroxenes.     

Rare earth element partitioning between titanite and silicate melts: Henry’s law revisited

We present detailed experimental results on the partitioning of rare earth elements (REE) between titanite and a range of different silicate melts. Our results show that Henry’s law of trace element partitioning depends on bulk composition, the available partners for heterovalent substitution, crystal composition, and melt composition. We illustrate that the partition coefficients for Sm depend very strongly on the bulk concentration of Sm in the system. The substitution mechanism, by which rare earth elements are incorporated into the crystal structure, plays an important role for trace element partitioning and also for the onset of Henry’s law. Our data show that there are clear differences between substitution mechanisms of major elements compared to elements which are present only as traces. Our experiments also clearly show that the onset of Henry’s law depends on the concentrations of the sum of all trace elements which are incorporated into the crystal by the same substitution mechanism. For geochemical modelling of magmatic processes involving titanite, and indeed other accessory phases, it is of crucial importance to first evaluate whether the REE, and other trace elements, are present as traces or as major elements, only then appropriate D values may be chosen.     

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.     

Trace element partitioning between carbonatitic melts and mantle transition zone minerals: Implications for the source of carbonatites

The occurrence of CO₂-rich lavas (carbonatites, kimberlites) and carbonate-rich xenoliths provide evidence for the existence of carbonatitic melts in the mantle. To model the chemical composition of such melts in the deep mantle, we experimentally determined partition coefficients for 23 trace elements (including REE, U-Th, HFSE, LILE) between deep mantle minerals and carbonatite liquids at 20 and 25 GPa and 1600 °C. Under these conditions, majoritic garnet and CaSiO₃ perovskite are the main reservoirs for trace elements. This study used both femtosecond LA-ICP-MS and SIMS techniques to measure reliable trace element concentrations. Comparison of the two techniques shows a general agreement, except for Sc and Ba. Our experimentally determined partition coefficients are consistent with the lattice strain model. The data suggest an effect of melt structure on partition coefficients in this pressure range. For instance, strain-free partition coefficient (D₀) for majorite-carbonatite melts do not follow the order of cation valence, , observed for majorite-CO₂-free silicate melts. The newly determined partition coefficients were combined with trace element composition of majoritic garnets found as inclusions in diamond to model trace element patterns of deep-seated carbonatites. The result compares favorably with natural carbonatites. This suggests that carbonatites can originate from the mantle transition zone.     

An experimental study of rare earth element partitioning between a shergottite melt and pigeonite: implications for the oxygen fugacity of the mart ian interior

A series of experiments on a synthetic, pigeonite-saturated, basaltic shergottite were carried out to constrain the variation of D(Eu/Gd)_{pigeonite/melt} and D(Eu/Sm)_{pigeonite/melt} with oxygen fugacity (f_O2). The experiments have been run under both dry and hydrous conditions. The shergottite was doped with 0.1, 0.5, and 1.0 wt.% Eu, Gd, and Sm oxides in different experiments and was equilibrated at liquidus conditions for 24 hours. D(Eu/Gd)_{pigeonite/melt} in dry melts ranges from 0.156 ± 0.014 (f_O2 = IW − 1) to 0.630 ± 0.102 (IW + 3.5). D(Eu/Sm)_{pigeonite/melt} in dry melts ranges from 0.279 ± 0.021 (IW − 1) to 1.114 ± 0.072 (IW + 3.5). Due to difficulties with low-f_O2 experiments, hydrous distribution coefficients were measured, but were not used in the calibration of the Eu-oxybarometers. These two Eu-oxybarometers provide an accurate way to measure f_O2 recorded during pigeonite crystallization, thereby yielding a record of f_O2 during the earliest period of Martian meteorite parent magma crystallization history.Using this new calibration, Martian meteorite pigeonite cores record f_O2 values of IW − 0.6 (±0.3) (QUE94201) to IW + 1.9 (±0.6) (Shergotty). These new values differ in magnitude, but not trend, from previously published data. The pigeonite Eu-oxybarometer yields an f_O2 range in the basaltic shergottites of 2 to 3 orders of magnitude. Several processes have been proposed to explain the origin of this f_O2 range, the majority of which rely on assimilation of an oxidized source. A potential correlation between this new pigeonite data and recent Fe-Ti oxide data, however, is consistent with intrinsic f_O2 differences in the magma source region being responsible for the measured f_O2 variations. This implies that the Martian meteorite source region, the mantle or lithosphere, may be heterogeneous in nature. However, the process of assimilation cannot be completely ruled out in that an assimilation event that took place before crystallization commenced would result in the overprinting of the source region f_O2 signature.     

Experimental determination of rare earth element partitioning between hydrous silicate melt,amphibole and garnet peridotite minerals at upper mantle pressures and temperatures

Partition coefficients of Ce, Sm and Tm involving garnet peridotite minerals, amphibole and hydrous silicate melt have been determined experimentally in the temperature and pressure ranges 950–1075°C and 10–25 kbar.Only several parts per million to several tens of parts per million of rare earth element (REE) can dissolve in the minerals before the crystal-liquid partition coefficients begin to vary as a function of REE content. The concentration ranges of constant partition coefficient increase with increasing temperature and are also positively correlated with the magnitude of the crystal-liquid partition coefficients. The upper concentration limits of constant partition coefficient and the value of the crystal-liquid partition coefficient for REE decrease in the order garnet > clinopyroxene > amphibole > orthopyroxene > olivine.Partition coefficients may vary by at least an order of magnitude as a function of bulk composition of the liquid phase (e.g. changing from basaltic to andesitic). The approximate ranges of the values of the partition coefficients as a function of bulk liquid composition are as follows:

$\text{Ce}\text{Sm}\text{Tm}\text{K}{}^{\mathrm{<mtext>ga-liq</mtext>}}\text{0.01–0.1}\text{0.3–3.4}\text{1–10}\text{K}{}^{\mathrm{<mtext>cpx-liq</mtext>}}\text{0.05–0.4}\text{0.09–0.7}\text{0.04–0.4}\text{K}{}^{\mathrm{<mtext>amph-liq</mtext>}}\text{0.04–0.4}\text{0.08–0.8}\text{0.07–0.7}\text{K}{}^{\mathrm{<mtext>opx-liq</mtext>}}\text{0.04–0.1}\text{0.05–0.1}\text{0.08–0.1}\text{K}{}^{\mathrm{<mtext>ol-liq</mtext>}}\text{0.01–0.02}\text{0.01–0.02}\text{0.01–0.02}$

where the values increase with increasing acidity of the melt.     

Compositional effects on element partitioning between Mg-silicate perovskite and silicate melts

High-pressure melting experiments were performed at ~26 GPa and ~2,200–2,400°C on synthetic peridotite compositions with varying FeO and Al₂O₃ contents and on a synthetic CI chondrite analogue composition. Peridotite liquids show a crystallisation sequence of ferropericlase (Fp) followed down temperature by Mg-silicate perovskite (MgPv) + Fp, which contrasts a sequence of MgPv followed by MgPv + Fp observed in the chondritic composition. The difference in crystallisation sequence is a consequence of the different bulk Mg/Si ratios. MgPv/melt partition coefficients for major, minor and trace elements were determined by electron microprobe and secondary ion mass spectrometry. Partition coefficients of tri- and tetravalent elements increase with increasing Al concentration in MgPv. A lattice strain model indicates that Al³⁺ substitutes predominantly onto the Si-site in MgPv, whereas most elements substitute onto the Mg-site, which is consistent with a charge-compensating coupled substitution mechanism. MgPv/melt partition coefficients for Mg (D_Mg) and Si (D_Si) are related to the melt Mg/Si ratio such that D_Si becomes lower than D_Mg at low Mg/Si melt ratios. We use a crystal fractionation model, based on upper mantle refractory lithophile element ratios, to constrain the amount of MgPv and Ca-silicate perovskite (CaPv) that could have fractionated during a Hadean magma ocean event and could still be present as a chemically distinct heterogeneity in the lower mantle today. We show that a fractionated crystal pile composed of 96% MgPv and 4% CaPv could comprise up to 13 wt% of the entire mantle.     

Determination of melt influence on divalent element partitioning between anorthite and CMAS melts

We propose a theory for crystal-melt trace element partitioning that considers the energetic consequences of crystal-lattice strain, of multi-component major-element silicate liquid mixing, and of trace-element activity coefficients in melts. We demonstrate application of the theory using newly determined partition coefficients for Ca, Mg, Sr, and Ba between pure anorthite and seven CMAS liquid compositions at 1330 °C and 1 atm. By selecting a range of melt compositions in equilibrium with a common crystal composition at equal liquidus temperature and pressure, we have isolated the contribution of melt composition to divalent trace element partitioning in this simple system. The partitioning data are fit to Onuma curves with parameterizations that can be thermodynamically rationalized in terms of the melt major element activity product (a_Al2O3)(a_SiO2)² and lattice strain theory modeling. Residuals between observed partition coefficients and the lattice strain plus major oxide melt activity model are then attributed to non-ideality of trace constituents in the liquids. The activity coefficients of the trace species in the melt are found to vary systematically with composition. Accounting for the major and trace element thermodynamics in the melt allows a good fit in which the parameters of the crystal-lattice strain model are independent of melt composition.     

An experimental determination of rare earth partition coefficients between a chloride containing vapor phase and silicate melts

The partitioning behavior of cerium, europium, gadolinium and ytterbium between an aqueous “vapor” phase and water saturated silicate melt have been experimentally examined using a new experimental approach employing radioactive tracers and a double-capsule technique. Equilibrium was established by reversing the partition coefficient¹ and by betatrack autoradiography. Aqueous solution compositions were varied by adding different amounts of chloride and in some cases fluoride or carbon dioxide. The H₂O contents of the Spruce Pine pegmatite melts were varied by conducting experiments at 4.0 kb, 800°C and at 1.25 kb, 800°C. A jadeite-nepheline composition (75 wt% Jadeite) also was employed at 4.0 kb, 800°C.The chloride experiments (Spruce Pine 4 kb, 800°C) show a linear relationship between the cube of the chloride molality and the partition coefficients of the trivalent rare earths. Europium, under the experimental fO₂ conditions (quartz-fayalite-magnetite buffer), varied linearly as the fifth power of the chloride molality. At the chloride molalities examined (<1.1 mC1), all the rare earths partitioned preferentially into the melt phase (K_P^RE <1). Relative to pure water, the presence of chloride and fluoride fon increased the partitioning of the individual rare earths into the vapor phase, while carbon dioxide did not. Europium anomalies were recorded 1n all experiments, particularly those involving the Spruce P1ne melt at 4.0 kb and 800°C which displayed a large positive europium anomaly at all chloride molalities. Furthermore, a relative fractionation of the trivalent rare earths was also observed in these experiments, such that K_P^Ce>K_P^Gd>K_P^Yb. The smaller ytterbium ion was consistently concentrated in the melt phase relative to the other rare earths in all experiments on the Spruce Pine composition. Experiments on the jadeite-nepheline composition showed no relative fractionation and a positive europium anomaly. The 1.25 kb experiment on the Spruce Pine composition showed a negative europium anomaly in plots of K_p^RE vs. REE.The overall rare earth partitioning at a constant chloride molality (mCl = .914) was such that K_P^SP(1.25 kb) > K_P^SP(4.0 kb) > K_P^Jd-Ne(4.0 kb), where SP = Spruce Pine, Jd-Ne = jadeitenepheli Using the model of Burnnam (1975), It is suggested that the trivalent rare earth partitioning is related to the cube of the melt octahedral site concentration; a property which 1n hydrous melts 1s dependent on melt composition and hydroxyl molality. Excellent agreement was found for the Spruce Pine melt, whereas the jadeite-nepheline melt gave apparent hydroxyl molalities which were too high for the measured partition coefficient. Additional octahedral sites are proposed for this unusual composition perhaps due to some aluminum in 6-fold coordination. The apparent compositional variation of europium partitioning at a constant oxygen fugacity is believed to be related to both the octahedral melt site concentration for trlvalent europium and an 8-coordinated site concentration for divalent europium. Any parameter which affects the numbers of these sites (_PH₂O, melt composition) will affect the rare earth partitioning. The observed dependency of the partition coefficient on the structural state of the melt could be as significant as its dependency on crystalline structural constraints. Furthermore, since _PH₂O can drastically effect the melt structural state, its effects could be reflected in melt/crystal partition coefficients.     

Distribution of titanium and the rare earth elements between peridotitic minerals

The concentrations of titanium and rare earth elements (REE) in olivines, orthopyroxenes, clinopyroxenes and spinels from four anhydrous, spinel-bearing peridotite xenoliths have been determined. The distribution of titanium (used as an analogue for the high field strength elements: HFSE) relative to the REE between clinopyroxenes and orthopyroxenes varies as a function of the whole rock composition and modal mineralogy. The distribution coefficients for titanium and the REE in these peridotites do not reflect mineral-melt equilibria. It is believed that subsolidus distribution coefficients for HFSE relative to REE vary with temperature. Ratios of various incompatible elements (e.g., Ti/Eu, Zr/Sm, Hf/Sm and P/Nd) in peridotite minerals differ from those in most primary basalts. However, the abundance ratios of incompatible elements in the bulk peridotite are comparable to those found in modern basalts. Given this and the differing contribution of melt from each phase during melting, near constant ratios of such incompatible elements in primary and primitive basalts and komatiites reflect the buffering of the melt by its residue. These ratios are fixed in the magma during the initial stages of melting because of similar and low distribution coefficients between melt and bulk residue for these element pairs. Differences in the relative abundances of titanium and REE in clinopyroxenes and orthopyroxenes demonstrate that mantle normalized abundance patterns for clinopyroxene are not equivalent to those of the whole rock. Therefore, claims of a widespread HFSE-depleted reservoir in the upper mantle base solely on the relative abundances of incompatible elements in peridotitic clinopyroxenes are invalid.     

Partitioning of rare earth elements between phosphate-rich carbonatite melts and mantle peridotites

Summary Near solidus equilibria in the system mantle peridotite-carbonate-phosphate doped with Ce and Yb have been studied at 20 kbar and 950°C. Carbonatitic melts in this system may be quenched into homogeneous glasses. Such melts intensely extract REE from rock-forming mantle minerals, and their migration may cause processes of mantle metasomatism.

---

Verteilung von Seltenen Erden zwischen phosphatreichen karbonatitischen Schmelzen und Mantel-Peridotiten

Zusammenfassung Gleichgewichte nahe dem Solidus im System Mantel-Peridotit-Karbonat-Phosphat, das mit Ce und Yb dotiert wurde, wurden bei 20 kbar und 950°C untersucht. Karbonatitische Schmelzen in diesem System können zu homogenen Gläsern abgeschreckt werden. Solche Schmelzen extrahieren SEE aus gesteinsbildenden Mantelmineralen und ihre Migration könnte für Vorgänge der Mantel-Metasomatose verantwortlich sein.

---

---

With 2 figures     

Trace element partitioning between type B CAI melts and melilite and spinel: Implications for trace element distribution during CAI formation

Calcium- and aluminum-rich inclusions (CAIs), occurring in chondritic meteorites and considered the oldest materials in the solar system, can provide critical information about the environment and time scale of creation of planetary materials. However, interpretation of the trace element and isotope compositions of CAIs, particularly the light elements Li, Be, and B, is hampered by the lack of constraint on melilite-melt and spinel-melt partition coefficients. We determined melilite-melt and spinel-melt partition coefficients for 21 elements by performing controlled cooling rate (2 °C/h) experiments at 1 atmosphere pressure in sealed platinum capsules using a synthetic type B CAI melt. Trace element concentrations were measured by secondary ion mass spectrometry (SIMS) and/or laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Melilites vary only slightly in composition, ranging from Åk_31-43. Results for the partitioning of trace elements between melilite and melt in three experiments and between spinel and melt in two experiments show that partition coefficients are independent of trace element concentration, are in good agreement for different analytical techniques (SIMS and LA-ICP-MS), and are in agreement with previous measurements in the literature. Partition coefficients between intermediate composition melilites and CAI melt are the following: Li, 0.5; Be, 1.0; B, 0.22; Rb, 0.012; Sr, 0.68; Zr, 0.004; Nb, 0.003; Cs, 0.002; Ba, 0.018; La, 0.056; Nd, 0.065; Sm, 0.073; Eu, 0.67; Er, 0.037; Yb, 0.018; Hf, 0.001; Ta, 0.003; Pb, 0.15; U, 0.001; Th, 0.002. Site size energetics analysis is used to assess isovalent partitioning into the different cation sites. The Young’s modulus deduced from +2 cations partitioning into the melilite X site agrees well with the bulk modulus of melilite based on X-ray diffraction methods. The changes in light element partitioning as melilite composition varies are predicted and used in several models of fractional crystallization to evaluate if the observed Li, Be, and B systematics in Allende CAI 3529-41 are consistent with crystallization from a melt. Models of crystallization agree reasonably well with observed light element variations in areas previously interpreted to be unperturbed by secondary processes [Chaussidon, M., Robert, F., McKeegan, K.D., 2006. Li and B isotopic variations in an Allende CAI: Evidence for the in situ decay of short-lived ¹⁰Be and for the possible presence of the short-lived nuclide ⁷Be in the early solar system. Geochim. Cosmochim. Acta70, 224-245], indicating that the trends of light elements could reflect fractional crystallization of a melt. In contrast, areas interpreted to have been affected by alteration processes are not consistent with crystallization models.     

Rare earth and high field strength element partitioning between iron-rich clinopyroxenes and felsic liquids

Rare earth elements are commonly assumed to substitute only for Ca in clinopyroxene because of the similarity of ionic radii for REE³⁺ and Ca²⁺ in eightfold coordination. The assumption is valid for Mg-rich clinopyroxenes for which observed mineral/melt partition coefficients are readily predicted by the lattice strain model for substitution onto a single site (e.g., Wood and Blundy 1997). We show that natural Fe-rich pyroxenes in both silica-undersaturated and silica-oversaturated magmatic systems deviate from this behavior. Salites (Mg# 48–59) in phonolites from Tenerife, ferrohedenbergites (Mg# 14.2–16.2) from the rhyolitic Bandelier Tuff, and ferroaugites (Mg# 9.6–32) from the rhyolitic Rattlesnake Tuff have higher heavy REE contents than predicted by single-site substitution. The ionic radius of Fe²⁺ in sixfold coordination is substantially greater than that of Mg²⁺; hence, we propose that, in Fe-rich clinopyroxenes, heavy REE are significantly partitioned between eightfold Ca sites and sixfold Fe and Mg sites such that Yb and Lu exist dominantly in sixfold coordination. We also outline a REE-based method of identifying pyroxene/melt pairs in systems with multiple liquid and crystal populations, based upon the assumption that LREE and MREE reside exclusively in eightfold coordination in pyroxene. Contrary to expectations, interpolation of mineral/melt partition coefficient data for heavy REE does not predict the behavior of Y. We speculate that mass fractionation effects play a role in mineral/melt lithophile trace element partitioning that is detectable among pairs of isovalent elements with near-identical radii, such as Y and Ho, Zr and Hf, and Nb and Ta.     

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.     

Trace element partitioning and melt structure: An experimental study at 1 atm pressure

Diopside-melt and forsterite-melt rare earth (REE) and Ni partition coefficients have been determined as a function of bulk compositions of the melt. Available Raman spectroscopic data have been used to determine the structures of the melts coexisting with diopside and forsterite. The compositional dependence of the partition coefficients is then related to the structural changes of the melt.The melts in all experiments have a ratio of nonbridging oxygens to tetrahedral cations ($\text{NBO}\text{T}$) between 1 and 0. The quenched melts consist of structural units that have, on the average, 2 (chain), 1 (sheet) and 0 (three-dimensional network) nonbridging oxygens per tetrahedral cation. The proportions of these structural units in the melts, as well as the overall $\text{NBO}\text{T}$, change as a function of the bulk composition of the melt.It has been found that Ce, Sm, Tm and Ni crystal-liquid partition coefficients ($\text{K}{}^{\mathrm{crystal?liq}}{}_{\mathrm{i}}\text{=}\text{C}{}^{\mathrm{crystal}}{}_{\mathrm{i}}\text{C}{}^{\mathrm{liq}}{}_{\mathrm{i}}$) decrease linearly with increasing $\text{NBO}\text{T}$. The values of the individual REE crystal-liquid trace element partition coefficients have different functional relations to $\text{NBO}\text{T}$, so that the degree of light REE enrichment of the melts would depend on their $\text{NBO}\text{T}$.The solution mechanisms of minor oxides such as CO₂, H₂O, TiO₂, P₂O₅ and Fe₂O₃ in silicate melts are known. These data have been recast as changes of $\text{NBO}\text{T}$ of the melts with regard to the type of oxide and its concentration in the melt. From such data the dependence of crystal-liquid partition coefficients on concentration and type of minor oxide in melt solution has been calculated.