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.     

Effect of melt composition on the partitioning of trace elements between titanite and silicate melt

Isobaric and isothermal experiments were performed to investigate the effect of melt composition on the partitioning of trace elements between titanite (CaTiSiO₅) and a range of different silicate melts. Titanite-melt partition coefficients for 18 trace elements were determined by secondary ion mass spectrometry (SIMS) analyses of experimental run products. The partition coefficients for the rare earth elements and for Th, Nb, and Ta reveal a strong influence of melt composition on partition coefficients, whereas partition coefficients for other studied monovalent, divalent and most quadrivalent (i.e., Zr, Hf) cations are not significantly affected by melt composition. The present data show that the influence of melt composition may not be neglected when modelling trace element partitioning.It is argued that it is mainly the change of coordination number and the regularity of the coordination space of trace elements in the melt structure that controls partition coefficients in our experiments. Furthermore, our data also show that the substitution mechanism by which trace elements are incorporated into titanite crystals may be of additional importance in this context.     

The effect of liquid composition on the partitioning of Ni between olivine and silicate melt

We report the results of experiments designed to separate the effects of temperature and pressure from liquid composition on the partitioning of Ni between olivine and liquid, \(D_{\text{Ni}}^{\text{ol/liq}}\). Experiments were performed from 1300 to 1600 °C and 1 atm to 3.0 GPa, using mid-ocean ridge basalt (MORB) glass surrounded by powdered olivine in graphite–Pt double capsules at high pressure and powdered MORB in crucibles fabricated from single crystals of San Carlos olivine at one atmosphere. In these experiments, pressure and temperature were varied in such a way that we produced a series of liquids, each with an approximately constant composition (~12, ~15, and ~21 wt% MgO). Previously, we used a similar approach to show that \(D_{\text{Ni}}^{\text{ol/liq}}\) for a liquid with ~18 wt% MgO is a strong function of temperature. Combining the new data presented here with our previous results allows us to separate the effects of temperature from composition. We fit our data based on a Ni–Mg exchange reaction, which yields \(\ln \left( {D_{\text{Ni}}^{\text{molar}} } \right) = \frac{{ -\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{RT} + \frac{{\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{R} - \ln \left( {\frac{{X_{\text{MgO}}^{\text{liq}} }}{{X_{{{\text{MgSi}}_{ 0. 5} {\text{O}}_{ 2} }}^{\text{ol}} }}} \right).\) Each subset of constant composition experiments displays roughly the same temperature dependence of \(D_{\text{Ni}}^{\text{ol/liq}}\) (i.e.,\(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\)) as previously reported for liquids with ~18 wt% MgO. Fitting new data presented here (15 experiments) in conjunction with our 13 previously published experiments (those with ~18 wt% MgO in the silicate liquid) to the above expression gives \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 3641 ± 396 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 1.597 ± 0.229. Adding data from the literature yields \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 4505 ± 196 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 2.075 ± 0.120, a set of coefficients that leads to a predictive equation for \(D_{\text{Ni}}^{\text{ol/liq}}\) applicable to a wide range of melt compositions. We use the results of our work to model the melting of peridotite beneath lithosphere of varying thickness and show that: (1) a positive correlation between NiO in magnesian olivine phenocrysts and lithospheric thickness is expected given a temperature-dependent \(D_{\text{Ni}}^{\text{ol/liq}} ,\) and (2) the magnitude of the slope for natural samples is consistent with our experimentally determined temperature dependence. Alternative processes to generate the positive correlation between NiO in magnesian olivines and lithospheric thickness, such as the melting of olivine-free pyroxenite, are possible, but they are not required to explain the observed correlation of NiO concentration in initially crystallizing olivine with lithospheric thickness.     

The univariant curve liquid = forsterite + anorthite + diopside in the system CMAS at 1 bar: solid solutions and melt structure

One atmosphere liquid-present experiments were carried out in the CMAS system using an ordinary quench furnace apparatus. The runs, including reversal and duplicate experiments, describe the univariant curve l=fo+an+ di between the invariant points Q: l+fo+an=di+sp and F: l+fo=an+di+oen, located respectively at 1245±1° C and 1244±1° C. The thermal divide on this curve M₃: l= fo+an+di is located at 1275±1° C and plots well within the silica-saturated field, in agreement with Longhi's (1987) experiments. Along the univariant curve l=fo+an+di, liquid composition evolves away from the thermal divide either toward invariant points Q or F and pierces the silica saturation plane, i.e., the join Di-An-En, in the silica saturated field. In this compositional range, the Al solubility in clinopyroxene changes drastically from one side of the thermal divide to the other, with great increase of Al solubility in the silica-undersaturated field. Four endmembers must be used to describe the complex solid solution of anorthite: CaAl₂Si₂O₈, CaMgSi₃O₈, MgAl₂Si₂O₈ and [] Si₄O₈. The last two of these are present only within the silica-saturated field. Unlike clinopyroxene, the Mg content of anorthite is insensitive to the thermal barrier but is only sensitive to silica-saturation plane. Olivine composition can be described by a binary solid solution of forsterite and monticellite with no Ca in the M₂ site. As with anorthite, olivine compositions exhibit a marked change with crossing of the silica-saturation plane. The above features imply that the solubility of minor elements in crystalline phases (Al in clinopyroxene, Ca in olivine and Mg in anorthite) selectively respond to only one or another of these particular plane. Results have many important consequences. One is the likelihood of changes in melt speciation depending on position with respect to the thermal divide and the silica-saturation plane.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Influence of Al for Si substitution in the tetrahedral network of silicate melts

The influence on olivine/melt transition metal (Mn, Co, Ni) partitioning of substitution in the tetrahedral network of silicate melt structure has been examined at ambient pressure in the 1450-1550 °C temperature range. Experiments were conducted in the systems NaAlSiO₄-Mg₂SiO₄- SiO₂ and CaAl₂Si₂O₈-Mg₂SiO₄-SiO₂ with about 1 wt% each of MnO, CoO, and NiO added. These compositions were used to evaluate how, in silicate melts, substitution and ionization potential of charge-balancing cations affect activity-composition relations in silicate melts and mineral/melt partitioning.The exchange equilibrium coefficient, , is a positive and linear function of melt Al/(Al + Si) at constant degree of melt polymerization, NBO/T. The is negatively correlated with the ionic radius, r, of the M-cation and also with the ionization potential (Z/r², Z = electrical charge) of the cation that serves to charge-balance Al³⁺ in tetrahedral coordination in the melts. The activity coefficient ratio, (γ_M/γ_Mg)^melt, is therefore similarly correlated.These melt composition relationships are governed by the distribution of Al³⁺ among coexisting Q-species in the peralkaline (depolymerized) melts coexisting with olivine. This distribution controls Q-speciation abundance, which, in turn, controls (γ_M/γ_Mg)^melt and . The relations between melt structure and olivine/melt partitioning behavior lead to the suggestion that in natural magmatic systems mineral/melt partition coefficients are more dependent on melt composition and, therefore, melt structure the more alkali-rich and the more felsic the melt. Moreover, mineral/melt partition coefficients are more sensitive to melt composition the more highly charged or the smaller the ionic radius of the cation of interest.     

The partitioning of Zr and Nb between diopside and melts in the system diopside-albite-anorthite

The crystal liquid partitioning of Zr and Nb has been measured experimentally between diopsidic clinopyroxene and melts in the system Di-Ab-An. Nb was found to be excluded from diopside (D(Nb) is always less than 0.02). D(Zr) is quite variable, ranging from 0.05 to 0.45. D(Zr) is positively correlated with the Al content of both the melt and the pyroxene and is negatively correlated with temperature. Both D(Zr) and D(Nb) were found to be independent of oxygen fugacity. This implies that neither Zr or Nb suffer valence changes over a range of oxygen fugacities spanning both lunar and terrestrial conditions.     

The influence of melt structure on trace element partitioning near the peridotite solidus

This experimental study examines the mineral/melt partitioning of Na, Ti, La, Sm, Ho, and Lu among high-Ca clinopyroxene, plagioclase, and silicate melts analogous to varying degrees of peridotite partial melting. Experiments performed at a pressure of 1.5 GPa and temperatures of 1,285 to 1,345 °C produced silicate melts saturated with high-Ca clinopyroxene, plagioclase and/or spinel, and, in one case, orthopyroxene and garnet. Partition coefficients measured in experiments in which clinopyroxene coexists with basaltic melt containing ~18 to 19 wt% Al₂O₃ and up to ~3 wt% Na₂O are consistent with those determined experimentally in a majority of the previous studies, with values of ~0.05 for the light rare earths and of ~0.70 for the heavy rare earths. The magnitudes of clinopyroxene/melt partition coefficients for the rare earth elements correlate with pyroxene composition in these experiments, and relative compatibilities are consistent with the effects of lattice strain energy. Clinopyroxene/melt partition coefficients measured in experiments in which the melt contains ~20 wt% Al₂O₃ and ~4 to 8 wt% Na₂O are unusually large (e.g., values for Lu of up to 1.33±0.05) and are not consistent with the dependence on pyroxene composition found in previous studies. The magnitudes of the partition coefficients measured in these experiments correlate with the degree of polymerization of the melt, rather than with crystal composition, indicating a significant melt structural influence on trace element partitioning. The ratio of non-bridging oxygens to tetrahedrally coordinated cations (NBO/T) in the melt provides a measure of this effect; melt structure has a significant influence on trace element compatibility only for values of NBO/T less than ~0.49. This result suggests that when ascending peridotite intersects the solidus at relatively low pressures (~1.5 GPa or less), the compatibility of trace elements in the residual solid varies significantly during the initial stages of partial melting in response to the changing liquid composition. It is unlikely that this effect is important at higher pressures due to the increased compatibility of SiO₂, Na₂O, and Al₂O₃ in the residual peridotite, and correspondingly larger values of NBO/T for small degree partial melts.Editorial responsibility: T.L. Grove     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

A possible effect of melt structure on the Mg-Fe partitioning between olivine and melt

Partitioning of Mg and Fe²⁺ between olivine and mafic melts has been determined experimentally for eight different synthetic compositions in the temperature range between 1335 and 1425°C at 0.1 MPa pressure and at fo₂ ∼1 log unit below the quartz-fayalite-magnetite buffer. The partition coefficient [K_D = (Fe²⁺/Mg)^ol/(Fe²⁺/Mg)^melt] increases from 0.25 to 0.34 with increasing depolymerization of melt (NBO/T of melt from 0.25-1.2), and then decreases with further depolymerization of melt (NBO/T from 1.2-2.8). These variations are similar to those observed in natural basalt-peridotite systems. In particular, the variation in NBO/T ranges for basaltic-picritic melts (0.4-1.5) is nearly identical to that obtained in the present experiments. Because the present experiments were carried out at constant pressure (0.1 MPa) and in a relatively small temperature range (90°C), the observed variations of Mg and Fe²⁺ partitioning between olivine and melt must depend primarily on the composition or structure of melt. Such variations of K_D may depend on the relative proportions of four-, five-, and six-coordinated Mg²⁺ and Fe²⁺ in melt as a function of degree of NBO/T.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Melt polymerization and Q-speciation in alkaline earth silicate systems

The two most abundant network-modifying cations in magmatic liquids are Ca²⁺ and Mg²⁺. To evaluate the influence of melt structure on exchange of Ca²⁺ and Mg²⁺ with other geochemically important divalent cations (m-cations) between coexisting minerals and melts, high-temperature (1470-1650 °C), ambient-pressure (0.1 MPa) forsterite/melt partitioning experiments were carried out in the system Mg₂SiO₄-CaMgSi₂O₆-SiO₂ with ?1 wt% m-cations (Mn²⁺, Co²⁺, and Ni²⁺) substituting for Ca²⁺ and Mg²⁺. The bulk melt NBO/Si-range (NBO/Si: nonbridging oxygen per silicon) of melt in equilibrium with forsterite was between 1.89 and 2.74. In this NBO/Si-range, the NBO/Si(Ca) (fraction of nonbridging oxygens, NBO, that form bonds with Ca²⁺, Ca²⁺-NBO) is linearly related to NBO/Si, whereas fraction of Mg²⁺-NBO bonds is essentially independent of NBO/Si. For individual m-cations, rate of change of K_D(m−Mg) with NBO/Si(Ca) for the exchange equilibrium, m^melt + Mg^olivine ? m^olivine + Mg^melt, is linear. K_D(m−Mg) decreases as an exponential function of increasing ionic potential, Z/r² (Z: formal electrical charge, r: ionic radius—here calculated with oxygen in sixfold coordination around the divalent cations) of the m-cation. The enthalpy change of the exchange equilibrium, ΔH, decreases linearly with increasing Z/r² [ΔH = 261(9)-81(3)·Z/r² (Å⁻²)]. From existing information on (Ca,Mg)O-SiO₂ melt structure at ambient pressure, these relationships are understood by considering the exchange of divalent cations that form bonds with nonbridging oxygen in individual Qⁿ-species in the melts. The negative ∂K_D(m−Mg)/∂(Z/r²) and ∂(ΔH)/∂(Z/r²) is because increasing Z/r² is because the cations forming bonds with nonbridging oxygen in increasingly depolymerized Qⁿ-species where steric hindrance is decreasingly important. In other words, principles of ionic size/site mismatch commonly observed for trace and minor elements in crystals, also govern their solubility behavior in silicate melts.     

The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

The effect of melt composition on the wetting angle between silicate melts and olivine

The wetting angle between silicate melts containing Ca, Li, Na, or K and olivine single crystals have been measured as part of an investigation of the dependence of the solid-liquid interfacial energy on melt composition and olivine orientation. The wetting angle increases with increasing silica content of the melt on (100) surfaces, but decreases with increasing silica content on (010) and (001) surfaces. For a given silica content, the wetting angle on (100) decreases in going from Ca to Li to Na to K, while the wetting angle on (010) and (001) increases in going from Ca to K-bearing melts. Based on published values for liquid-vapor interfacial energies, the observed changes in wetting angle with changes in melt composition indicate that the solid-liquid interfacial energy increases with increasing silica content of the melt for the (100) surface. However, for (010) and (001) surfaces, the variation of the solid-liquid interfacial energy with silica content depends upon whether Ca or K is present in the melt. In addition, the solid-liquid interfacial energy depends upon the orientation of the olivine in the following manner: _sl ⁽⁰¹⁰⁾ _sl ⁽⁰⁰¹⁾ _sl ⁽¹⁰⁰⁾ .     

The dependence of the partitioning of iron and europium between plagioclase and hydrous tonalitic melt on oxygen fugacity

The dependence of iron and europium partitioning between plagioclase and melt on oxygen fugacity was studied in the system SiO₂(Qz)—NaAlSi₃O₈(Ab)—CaAl₂Si₂O₈(An)—H₂O. Experiments were performed at 500 MPa and 850 °C/750 °C under water saturated conditions. The oxygen fugacity was varied in the log f _O2-range from −7.27 to −15.78. To work at the most reducing conditions the classical double-capsule technique was modified. The sample and a C—O—H bearing sensor capsule were placed next to each other within a BN jacket to minimise loss of hydrogen to the vessel atmosphere. By this setup redox conditions slightly more reducing than the FeO—Fe₃O₄ buffer could be maintained even in 96 h runs. Raman spectra showed that the BN was modified by reaction with hydrogen resulting in a low hydrogen permeability. The partition coefficients determined for Eu at 850 °C and 500 MPa vary from 0.095 at conditions of the Cu—Cu₂O buffer to 1.81 at the most reducing conditions (C—O—H sensor). In the same f _O2 interval the partition coefficient for Fe varies from 0.55 at oxidising conditions to 0.08 at the most reducing conditions. The partitioning of Sm, which was added as a reference for a trivalent REE, does not vary with the oxygen fugacity, yielding an average value for D = 0.07. Lowering the temperature to 750 °C for a given f _O2 decreases the partition coefficient of Eu and increases that of Fe. Comparison with published data at 1 atm and at higher temperatures shows that both temperature and composition of the melt have strong effects on the partitioning behaviour. As the change of the partition coefficients in the geologically relevant f _O2 range is quite strong, element partitioning of Eu and Fe might be used to estimate redox conditions for the genesis of igneous rocks. Furthermore, by modelling the partitioning data it is possible to extract information about the redox state of the melt. Resulting ferric-ferrous ratios show significant differences from those predicted by empirical models. Received: 14 October 1998 / Received: 5 March 1999     

Effect of water on the fluorine and chlorine partitioning behavior between olivine and silicate melt

Halogens show a range from moderate (F) to highly (Cl, Br, I) volatile and incompatible behavior, which makes them excellent tracers for volatile transport processes in the Earth’s mantle. Experimentally determined fluorine and chlorine partitioning data between mantle minerals and silicate melt enable us to estimate Mid Ocean Ridge Basalt (MORB) and Ocean Island Basalt (OIB) source region concentrations for these elements. This study investigates the effect of varying small amounts of water on the fluorine and chlorine partitioning behavior at 1280?°C and 0.3 GPa between olivine and silicate melt in the Fe-free CMAS+F–Cl–Br–I–H₂O model system. Results show that, within the uncertainty of the analyses, water has no effect on the chlorine partitioning behavior for bulk water contents ranging from 0.03 (2) wt% H₂O (D_Cl ^ol/melt = 1.6?±?0.9 × 10^?4) to 0.33 (6) wt% H₂O (D_Cl ^ol/melt = 2.2?±?1.1 × 10^?4). Consequently, with the effect of pressure being negligible in the uppermost mantle (Joachim et al. Chem Geol 416:65–78, 2015), temperature is the only parameter that needs to be considered for the determination of chlorine partition coefficients between olivine and melt at least in the simplified iron-free CMAS+F–Cl–Br–I–H₂O system. In contrast, the fluorine partition coefficient increases linearly in this range and may be described at 1280?°C and 0.3 GPa with (R ²?=?0.99): \(D_{F}^{\text{ol/melt}}\ =\ 3.6\pm 0.4\ \times \ {{10}^{-3}}\ \times \ {{X}_{{{\text{H}}_{\text{2}}}\text{O}}}\left( \text{wt }\!\!\%\!\!\text{ } \right)\ +\ 6\ \pm \ 0.4\times \,{{10}^{-4}}\). The observed fluorine partitioning behavior supports the theory suggested by Crépisson et al. (Earth Planet Sci Lett 390:287–295, 2014) that fluorine and water are incorporated as clumped OH/F defects in the olivine structure. Results of this study further suggest that fluorine concentration estimates in OIB source regions are at least 10% lower than previously expected (Joachim et al. Chem Geol 416:65–78, 2015), implying that consideration of the effect of water on the fluorine partitioning behavior between Earth’s mantle minerals and silicate melt is vital for a correct estimation of fluorine abundances in OIB source regions. Estimates for MORB source fluorine concentrations as well as chlorine abundances in both mantle source regions are within uncertainty not affected by the presence of water.     

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.     

The effect of oxygen fugacity on the partitioning of Re between crystals and silicate melt during mantle melting

Interpretation of Re-Os isotopic systematics applied to mantle and mantle-derived rocks is currently hindered by the poorly understood behaviour of Re and Os during partial melting. Of particular interest is the incompatibility of Re and how it partitions between melt and the different mantle phases. Here, we study the partitioning behaviour of Re between the common upper mantle minerals (garnet, spinel, clinopyroxene, orthopyroxene, and olivine) and silicate melt under temperature (1275-1450 °C) and pressure (1.5-3.2 GPa) conditions relevant for basaltic magma genesis, over a range of oxygen fugacity (?O₂) large enough (QFM+5.6 to QFM−2.9) to demonstrate the effects of changing the oxidation state of Re from 4+ to 6+. Rhenium crystal/silicate-melt partition coefficients vary by 4-5 orders of magnitude, from moderately compatible to highly incompatible, for pyroxenes, garnet, and spinel as the oxidation state of Re changes from 4+ to 6+, but Re in either oxidation state is incompatible in olivine. Because the changeover from the one Re oxidation state to the other occurs over the range of ?O₂s pertinent to partial melting in the Earth’s mantle, bulk Re crystal/silicate-melt partition coefficients during mantle melting are also expected to vary significantly according to the oxidation state of the system. For instance, assuming QFM−0.7 and QFM+1.6 as average ?O₂ for mid-ocean ridge (MORBs) and island arc (IABs) basalts, respectively, a difference of at least one order of magnitude for bulk Re partition coefficients is expected (excluding any influence from a sulphide phase). Hence, Re is probably much more incompatible during the genesis of IABs compared to MORBs. Our results also demonstrate that Re⁴⁺ has a partitioning behaviour similar to Ti⁴⁺ rather than Yb, and is accordingly not a sensitive indicator of garnet in the source. The lower concentrations of Re observed in ocean island basalts (OIBs) compared to MORBs are therefore not a result of being generated deeper in the mantle where garnet is stable, leaving the hypothesis of late-stage loss of Re from OIB lavas by degassing as the most plausible explanation.     

Efremovka 101.1: a CAI with ultrarefractory REE patterns and enormous enrichments of Sc, Zr, and Y in Fassaite and Perovskite

Inclusion 101.1 from the CV3 carbonaceous chondrite Efremovka is a compact Type A Ca-Al-rich inclusion (CAI) highly enriched in ultrarefractory (UR) oxides. It is the first complete CAI with a UR rare earth element (REE) pattern found in a CV3 chondrite. The inclusion is petrographically complex and was formed in a multistage process. It consists of several lithologically unrelated units.The core contains abundant Y- and Zr-perovskite, Sc- and Zr-rich fassaite, and metallic FeNi enclosed in melilite. All mineral species (except spinel) in all lithological units exhibit the same basic UR REE pattern. Four different populations of perovskites are distinguished by different Y/Zr ratios. A few of the perovskites have Y/Zr ratios similar to those obtained from crystal/liquid fractionation experiments. Perovskites from the other three populations have either chondritic, lower than chondritic Y/Zr ratios or extremely low Zr contents. Ca isotopic ratios differ among three perovskites from different populations, demonstrating a variety of sources and formational processes. Most fassaites crystallized in situ through reaction between the CAI liquid and preexisting perovskites. This process induced redistribution of Zr, Y, Sc, and V between perovskite and fassaite, thus overprinting the original abundances in perovskite. Fassaite reaction rims around FeNi metals are also encountered. They are enriched in V, which was gained from the metal through oxidation of V in metal during fassaite crystallization. The relative abundances of Zr, Y, and Sc in perovskites are complementary to the abundances of these elements in Sc- and Zr-fassaite, indicating subsolidus partitioning of these elements between the two phases. Perovskites are enriched in Y and depleted in Sc and Zr in comparison to fassaites.The core contains two complete captured CAIs, several sinuous fragments, and fine-grained polygonal refractory fragments. An assemblage of andradite-wollastonite-hedenbergite and pure metallic iron is encountered as enclaves in the interior of some sinuous fragments. Metallic Fe and wollastonite formed by reduction of preexisting andradite and hedenbergite nebular alteration products upon inclusion in the highly reduced CAI melt. Numerous spinel clusters and framboids with varying V₂O₃ and Cr₂O₃ concentrations are enclosed in individual melilite crystals in the host CAI and captured CAIs. The rim sequence of the host consists of six layers (from the inside outward): (a) FeO-poor spinel, (b) Sc-bearing fassaite, (c) Al-diopside, (d) Al- and Ca-bearing olivine, (e) pure diopside, and (f) Ca-poor olivine. Like the constituents of the CAI core, all mineral layers of the rim sequence, except spinel, have the same UR REE pattern. However, the total REE abundances decrease systematically by 1 order of magnitude from layer 2 to layer 6. This feature strongly suggests formation of the rim sequence by successive condensation from a unique reservoir enriched in UR elements and excludes formation by flash heating. Petrography, mineral chemistry, REE, refractory lithophile element abundances, and Ca isotopic compositions demonstrate the complex multistage formation history of a CAI that on the surface looks like a regular Type A inclusion.     

The effect of sulfur on the partitioning of Ni and other first-row transition elements between olivine and silicate melt

The effect of sulfur dissolved as sulfide (S²⁻) in silicate melts on the activity coefficients of NiO and some other oxides of divalent cations (Ca, Cr, Mn, Fe and Co) has been determined from olivine/melt partitioning experiments at 1400 °C in six melt compositions in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS), and in derivatives of these compositions at 1370 °C, obtained from the six CMAS compositions by substituting Fe for Mg (FeCMAS). Amounts of S²⁻ were varied from zero to sulfide saturation, reaching 4100 μg g⁻¹ S in the most sulfur-rich silicate melt. The sulfide solubilities compare reasonably well with those predicted from the parameterization of the sulfide capacity of silicate melts at 1400 °C of O’Neill and Mavrogenes (2002), although in detail systematic deviations indicate that a more sophisticated model may improve the prediction of sulfide capacities.The results show a barely discernible effect of S²⁻ in the silicate melt on Fe, Co and Ni partition coefficients, and also surprisingly, a tiny but resolvable effect on Ca partitioning, but no detectable effect on Cr, Mn or some other lithophile incompatible elements (Sc, Ti, V, Y, Zr and Hf). Decreasing Mg# of olivine (reflecting increasing FeO in the system) has a significant influence on the partitioning of several of the divalent cations, particularly Ca and Ni. We find a remarkably systematic correlation between and the ionic radius of M²⁺, where M = Ca, Cr, Mn, Fe, Co or Ni, which is attributable to a simple relationship between size mismatch and excess free energies of mixing in Mg-rich olivine solid solutions.Neither the effect of S²⁻ nor of Mg#^ol is large enough by an order of magnitude to account for the reported variations of obtained from electron microprobe analyses of olivine/glass pairs from mid-ocean ridge basalts (MORBs). Comparing these MORB glass analyses with the Ni-MgO systematics of MORB from other studies in the literature, which were obtained using a variety of analytical techniques, shows that these electron microprobe analyses are anomalous. We suggest that the reported variation of with S content in MORB is an analytical artifact.Mass balance of melt and olivine compositions with the starting compositions shows that dissolved S²⁻ depresses the olivine liquidus of haplobasaltic silicate melts by 5.8 × 10⁻³ (±1.3 × 10⁻³) K per μg g⁻¹ of S²⁻, which is negligible in most contexts. We also present data for the partitioning of some incompatible trace elements (Sc, Ti, Y, Zr and Hf) between olivine and melt. The data for Sc and Y confirm previous results showing that and decrease with increasing SiO₂ content of the melt. Values of average 0.01 with most falling in the range 0.005-0.015. Zr and Hf are considerably more incompatible than Ti in olivine, with and about 10⁻³. The ratio / is well constrained at 0.611 ± 0.016.     

Experimental determination of crystal/melt partitioning of Ga and Ge in the system forsterite-anorthite-diopside

The crystal/liquid partitioning of Ga and Ge has been measured experimentally between forsterite, diopside, anorthite and spinel and melts in the pseudoternary system forsterite-anorthite-diopside at one atmosphere pressure and 1300°C. Gallium is incompatible with forsterite and diopside [D (Ga) = 0.024 and 0.19, respectively], is only slightly incompatible in anorthite [D (Ga) = 0.86] and is highly compatible in spinel [D (Ga) = 4.6]. The partition coefficient for Ge is within a factor of two of unity for forsterite, diopside and anorthite [D (Ge) = 0.62, 1.4 and 0.51, respectively], but Ge is incompatible in spinel [D (Ge) = 0.1]. The coefficients for the exchange of Ga and Al and the exchange of Ge and Si between minerals and melts generally are within a factor of two of unity, as is expected from the geochemical coherence of these element pairs in natural samples. The application of our results to the interpretation of natural basaltic and mantle samples from the Earth and basalts from the Moon and the Shergottite Parent Body demonstrates that it is possible to discriminate between different mantle source compositions using Ga/Al and Ge/Si ratios. The Ge variation among lunar mare basalts may be indicative of a heterogeneous lunar mantle. The substantial depletion of Ge in Chassigny relative to the other SNC meteorites may be evidence of either a heterogeneous Shergottite Parent Body (SPB) mantle, or of different geochemical behavior for Ge in the SPB.