Sulfur partitioning between apatite and melt and effect of sulfur on apatite solubility at oxidizing conditions

The effect of sulfur on phosphorus solubility in rhyolitic melt and the sulfur distribution between apatite, ±anhydrite, melt and fluid have been determined at 200 MPa and 800–1,100 °C via apatite crystallization and dissolution experiments. The presence of a small amount of sulfur in the system (0.5 wt.% S) under oxidizing conditions increases the solubility of phosphorus in the melt, probably due to changing calcium activity in the melt as a result of the formation of Ca-S complexing cations. Apatite solubility geothermometers tend to overestimate temperature in Ca-poor, S-bearing system at oxidizing conditions. In crystallization experiments, the sulfur content in apatite decreases with decreasing temperature and also with decreasing sulfur content of the melt. The sulfur partition coefficient between apatite and rhyolitic melt increases with decreasing temperature (Kd_S^apatite/melt=4.5–14.2 at T=1,100–900 °C) under sulfur-undersaturated conditions (no anhydrite). The sulfur partition coefficient is lower in anhydrite-saturated melt (~8 at 800 °C) than in anhydrite-undersaturated melt, suggesting that Kd_S^apatite/melt depends not only on the temperature but also on the sulfur content of the melt. These first results indicate that the sulfur content in apatite can be used to track the evolution of sulfur content in a magmatic system at oxidizing conditions.Editorial responsibility: J. Hoefs     

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.     

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 between silicate minerals and carbonatite at 25 kbar and application to mantle metasomatism

Summary Experiments at 25 kbar and 1000°C, on a model trace element-enriched carbonatite-eridotite mix, produced augite + pargasite ± garnet ± dolomite coexisting with a carbonatite melt. Proton microprobe analysis of the phases showed that key trace elements (Rb, Ba, Sr, Nb, Ta, Zr, Y and REE) all partitioned strongly into the melt (with the exception of Y, Ho and Lu in garnet), verifying that carbonatite is potentially a highly effective metasomatizing agent. The data also indicate that carbonatitic metasomatism will impart higher Ba/Rb, Ba/Nb, Nb/Ta, Sr/Ta, La/Ta, and lower Zr/Y, with little change to Sr/Nb, in affected mantle.

---

Spurenelementverteilung zwischen Silikatmineralen und Karbonatit bei 25 kbar: Anwendung für die Mantel-Metasomatose

Zusammenfassung Experimente mit einer Modell-mischung von Karbonatit-Peridotit, angereichert mit Spurenelementen, produzierten bei 25 kbar und 1000°C Augit + Pargasit ± Granat ±Dolomit coexistierend mit einer Karbonatitschmelze. Protonmikrosonden-Analyse der Phasen zeigte, dass alle Schlüsselspurenelemente (Rb, Ba, Sr, Nb, Ta, Zr, Y and REE) stark in der Schmelze angereichert werden (mit der Ausnahme von Y, Ho und Lu in Granat), was beweist, dass Karbonatit potentiell ein sehr effektives Agens für Metasomatose ist. Die Daten zeigen weiterhin, dass karbonatitische Metasomatose in betroffenen Mantel höhere Ba/Rb, Ba/Nb, Nb/Ta, Sr/Ta, La/Ta und niedrigere Zr/Y produziert, mit geringen Äderungen für Sr/Nb.

---

---

With 1 Figure     

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.     

Aluminum-dependent trace element partitioning in clinopyroxene

As technical advances have dramatically increased our ability to analyze trace elements, the need for more reliable data on the compositional dependence of trace element partitioning between minerals and melt has become increasingly important. The late-Cretaceous Carmacks Group of south central Yukon comprises a succession of primitive high-Mg ankaramitic lavas characterized by shoshonitic chemical affinities and containing large complexly zoned clinopyroxene phenocrysts. The compositional zonation of the clinopyroxene phenocrysts is characterized by relatively Fe-rich (Mg^# = Mg/(Mg + Fe) = 0.85), but mottled, cores surrounded by mantles of cyclically-zoned clinopyroxene whose Mg^# varies repeatedly between 0.9 and 0.80. These cyclically zoned clinopyroxene mantles appear to record the repeated influx and mixing of batches of primitive with more evolved magma in a deep sub-crustal (∼1.2 GPa) magma chamber(s). Laser ablation ICP-MS was used to analyze the trace element variation in these zoned clinopyroxenes. The results indicate more than a threefold variation in the absolute concentrations of Th, Zr, rare earth elements (REE), and Y within individual clinopyroxene phenocrysts, with no apparent change in the degree of REE or high field strength element (HFSE) fractionation. The variation in absolute abundances of trace elements correlates closely with the major element composition of the clinopyroxene, with the most enriched clinopyroxene having the lowest Mg^# and highest Al contents. The problem is that the amount of crystal fractionation required to explain the major element variation (∼20%) in these clinopyroxene phenocrysts cannot explain the increase in the abundance of the incompatible trace elements, which would require more than 70% crystal fractionation, if constant partition coefficients are assumed. The anomalous increase in incompatible trace elements appears to reflect an increase in their partition coefficients with increasing Al^IV in the clinopyroxene; with an increase in Al₂O₃ from 1.5 to 4.0 wt.% during ∼20% crystal fractionation over a temperature decrease of ∼100°C being associated with more that a threefold increase in the partition coefficients of Th, Zr, REE, and Y. The magnitude of these increases may indicate that the substitution of these trace elements into clinopyroxene is better modeled in some natural systems by a local charge balance model, rather than the distributed charge model that better replicates the results of annealed experiments. These findings indicate that the effect of Al on the partition coefficients of incompatible trace elements in clinopyroxene may be under appreciated in natural magmatic systems and that the application of experimentally determined clinopyroxene partition coefficients to natural systems must be done with caution.     

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.     

Barium partitioning between alkali feldspar and silicate liquid at high temperature and pressure

Barium partitioning between alkali feldspar and a natural trachyte liquid, enriched with barium, has been determined as a function of pressure and temperature from 10 to 25 kb and 900°–1100° C. Both long duration experiments and a re-equilibration experiment suggest close approach to equilibrium. Partition coefficients (D _Ba) decrease as both temperature and pressure increase (e.g., D _Ba changes from 8.71 at 10 kb, 900° C to 1.48 at 25 kb, 1100° C). Water activity also controls the barium partitioning with a marked decrease in D _Ba ^af/liq for addition of less than 0.8 wt% H₂O, but with no apparent additional effect for higher water contents in the bulk composition (e.g., from 0.8–4.2 wt% H₂O). The composition of alkali feldspar also has a significant effect on D _Ba ^af/liq , but the data obtained do not allow derivation of a complete D-Or relationship. These new data suggest that Henry's Law is obeyed for most of the barium concentrations examined, and the limit of Henry's Law behaviour for barium in alkali feldspar is as high as 6 wt% BaO in alkali feldspar and 1.2 wt% BaO in the melt, similar to the results of Long (1978). The experimental results broadly overlap with natural data for D _Ba, determined from coexisting alkali feldspar phenocrysts and glass (or groundmass).     

Chromium solubility in MgSiO3 ilmenite at high pressure

The crystal structure and chemical composition of crystals of (Mg_1?x Cr _x )(Si_1?x Cr _x )O₃ ilmenite (with x = 0.015, 0.023 and 0.038) synthesized in the model system Mg₃Cr₂Si₃O₁₂–Mg₄Si₄O₁₂ at 18–19 GPa and 1,600 °C have been investigated. Chromium was found as substitute for both Mg at the octahedral X site and Si at the octahedral Y site, according to the reaction Mg²⁺ + Si⁴⁺ = 2Cr³⁺. Such substitutions cause a shortening of the <X–O> and a lengthening of the <Y–O> distances with respect to the values typically observed for pure MgSiO₃ ilmenite and eskolaite Cr₂O₃. Although no high Cr contents are considered in the pyrolite model, Cr-bearing ilmenite may be the host for chromium in the Earth’s transition zone. The successful synthesis of ilmenite with high Cr contents and its structural characterization are of key importance because the study of its thermodynamic constants combined with the data on phase relations in the lower-mantle systems can help in the understanding of the seismic velocity and density profiles of the transition zone and the constraining composition and mineralogy of pyrolite in this area of the Earth.     

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.     

Low solubility of He and Ar in carbonatitic liquids: Implications for decoupling noble gas and lithophile isotope systems

In order to asses the importance of carbonatitic liquids in transporting noble gases in the mantle, the solubilities of He and Ar in carbonatitic liquids were estimated from analyses of calcium-potassium carbonate glasses that had been synthesized at 1 bar and temperatures between 850 and 950 °C under He or Ar enriched atmospheres. Despite poor reproducibility related to difficulties synthesizing carbonatite glass, we have been able to estimate He and Ar solubilities in carbonatite liquids to be 1 × 10⁻⁸ and 2 × 10⁻⁹ mol g⁻¹ at 1 bar respectively (with ?50% uncertainty). Despite the significant uncertainties on these estimates, it is clear that the noble gases are not massively soluble in carbonatite liquids (within error, these solubilities are identical to their equivalent solubilities in tholeiitic melts). Assuming the results of these low pressure experiments can be applied to mantle conditions, it seems unlikely that carbonatite metasomatism per se transports noble gases within the mantle. It is nevertheless possible that partitioning of lithophile trace elements (including the important radioelements, U, K and Th) and noble gases between a carbonatitic melt and a silicate melt could effectively decouple lithophile and noble gas isotope systematics because the carbonatitic melt expressedly does not transport noble gases, yet is known to efficiently transport incompatible trace elements.     

Oxygen isotope partitioning between phosphate and carbonate in mammalian apatite

The oxygen isotope compositions of phosphate and structural carbonate in mammalian enamel and bone apatite are linked to that of body water at constant body temperature near 37°C, but the isotope systematics of oxygen in structural carbonate are not well understood. Using coupled measurements of the oxygen isotope composition of structural carbonate and phosphate from horse tooth enamel, the apparent oxygen isotope fractionation factor between structural carbonate and body water is estimated to be 1.0263 ± 0.0014. These estimates provide a quantitative basis for using the oxygen isotope composition of structural carbonate in mammalian biogenic apatite for ecological, climatological, and physiological reconstruction.     

Element partitioning between olivine and silicate melt under high pressure

Element partitioning between olivine and silicate melt has been investigated at pressures 1–14 GPa, by using a 6–8 type multi-anvil high pressure apparatus. In order to observe systematics in the partitioning of trivalent ions, Li was added to the starting materials in order to increase the concentration of trivalent ions in olivine. With increasing pressure, it was found that partition coefficients of most of the elements gradually decreased. Trivalent ions generally showed parabolic pattern on partition coefficient — ionic radius diagram. When pyrolite-like material was used as the starting material, partition coefficient of Al, D_Al, gradually increased with increase in pressure while the partition coefficients of the other elements decreased, and the D_Al deviated from the parabolic pattern of other trivalent ions. The deviation of D_Al from the D pattern of the other trivalent ions was also found when olivine was employed as main component of the starting material. This result may be ascribed to the compositional change of coexisting silicate melt with increase in pressure.     

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.     

Constraints on core formation from Pt partitioning in mafic silicate liquids at high temperatures

A new high temperature piston cylinder design has enabled the measurement of platinum solubility in mafic melts at temperatures up to 2500 °C, 2.2 GPa pressure, and under reducing conditions for 1-10 h. These high temperature and low f_O2 conditions may mimic a magma ocean during planetary core formation. Under these conditions, we measured tens to hundreds of ppm Pt in the quenched silicate glass corresponding to , 4-12 orders of magnitude lower than extrapolations from high f_O2 experiments at 1 bar and at temperatures no higher than 1550 °C. Moreover, the new experiments provide coupled textural and compositional evidence that noble metal micro-nuggets, ubiquitous in experimental studies of the highly siderophile elements, can be produced on quench: we measure equally high Pt concentrations in the rapidly quenched nugget-free peripheral margin of the silicate as we do in the more slowly quenched nugget-bearing interior region. We find that both temperature and melt composition exercise strong control on and that Pt⁰ and Pt¹⁺ may contribute significantly to the total dissolved Pt such that low f_O2 does not imply low Pt solubility. Equilibration of metal alloy with liquid silicate in a hot primitive magma might not have depleted platinum to the extent previously believed.     

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.     

Fe---Mg partitioning between coexisting garnet and phengite at high pressure, and comments on a garnet-phengite geothermometer

The Fe---Mg exchange reaction between coexisting garnet and phengite has been studied by reacting a natural phengite (mg = 67) in the presence of quartz and water at pressures of 20–35 kb and temperatures of 800–1000°C. Compositions of coexisting garnet and mica indicate a linear relation between both the InK_D((Fe/Mg) garnet/(Fe/Mg) phengite) and temperature, and InK_D and pressure in the above P.T range. This Fe---Mg exchange reaction between garnet and phengite is shown to be dependent on the Ca-content of the garnet, and on the mg number of the bulk composition. These two composition effects have been studied by usin phengitic mica mixes with mg numbers of 23 and 46, and by using a synthetic basaltic composition. The overall results allow broad empirical calibration of separate geothermometers for pelitic and basaltic systems, respectively. However, because of non-ideality in the exchange reaction, this geothermometer should not be used in any practical application outside the composition ranges studied. Also, careful consideration of the presence of Fe³⁺ in phengite must be made. If the Fe³⁺ content of the natural phengite is unknown, then the temperatures obtained will be maximum temperatures only.     

Mg-Fe partitioning between olivine and ultramafic melts at high pressures

Precise determination of the partitioning of Mg and Fe²⁺ between olivine and ultramafic melt has been made at pressures from 5 to 13 GPa using a MA-8 type multi-anvil high-pressure apparatus (PREM) installed at Earthquake Research Institute, University of Tokyo. A very short rhenium capsule (<100 μm sample thickness) was adopted to minimize temperature variation within the sample container. Synthetic gels with the composition of the upper mantle peridotite were used as starting materials to promote the homogeneity. Analyses of quenched melts and coexisting olivines were made with an electron probe microanalyzer. The obtained partition coefficient, K_D [=(FeO/MgO)^ol/(FeO/MgO)^melt], decreases from 0.35 to 0.25 with increasing pressure from 5 to 13 GPa, suggesting a negative correlation between pressure and K_D above 5 GPa. Our result is consistent with a parabolic relationship between K_D and degree of polymerization (NBO/T) of melts reported by previous studies at lower pressures. The negative correlation between pressure and K_D suggests that olivine crystallizing in a magma ocean becomes more Mg-rich with depth and that primary magmas generated in the upper mantle become more Fe-rich with depth than previously estimated.