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.     

Redox equilibria of iron and silicate melt structure: Implications for olivine/melt element partitioning

Olivine/melt partitioning of ΣFe, Fe²⁺, Mg²⁺, Ca²⁺, Mn²⁺, Co²⁺, and Ni²⁺ has been determined in the systems CaO-MgO-FeO-Fe₂O₃-SiO₂ (FD) and CaO-MgO-FeO-Fe₂O₃-Al₂O₃-SiO₂ (FDA3) as a function of oxygen fugacity (f_O2) at 0.1 MPa pressure. Total iron oxide content of the starting materials was ∼20 wt%. The f_O2 was to used to control the Fe³⁺/ΣFe (ΣFe: total iron) of the melts. The Fe³⁺/ΣFe and structural roles of Fe²⁺ and Fe³⁺ were determined with ⁵⁷Fe resonant absorption Mössbauer spectroscopy. Changes in melt polymerization, NBO/T, as a function of f_O2 was estimated from the Mössbauer data and existing melt structure information. It varies by ∼100% in melts coexisting with olivine in the FDA3 system and by about 300% in the FD system in the Fe³⁺/ΣFe range of the experiments (0.805-0.092). The partition coefficients ( in olivine/wt% in melt) are systematic functions of f_O2 and, therefore, NBO/T of the melt. There is a -minimum in the FDA3 system at NBO/T-values corresponding to intermediate Fe³⁺/ΣFe (0.34-0.44). In the Al-free system, FD, where the NBO/T values of melts range between ∼1 and ∼2.9, the partition coefficients are positively correlated with NBO/T (decreasing Fe³⁺/ΣFe). These relationships are explained by consideration of solution behavior in the melts governed by Qⁿ-unit distribution and structural changes of the divalent cations in the melts (coordination number, complexing with Fe³⁺, and distortion of the polyhedra).     

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.     

Molybdenum, tungsten and manganese partitioning in the system pyrrhotite-Fe-S-O melt-rhyolite melt: Impact of sulfide segregation on arc magma evolution

Experiments were performed to determine the partitioning of molybdenum, tungsten and manganese among a rhyolitic melt (melt), pyrrhotite (po), and an immiscible Fe-S-O melt (Fe-S-O). Sulfide phases such as these may be isolated from a silicate melt along with other crystallizing phases during the evolution of arc magma, and partition coefficients are required to model the effect of this process on molybdenum and tungsten budgets.We developed an experimental design to take advantage of properties of the phases under study. Careful control of temperature allowed pyrrhotite and magnetite to be stable along with an Fe-S-O melt, and this phase assemblage allowed the composition of run-product pyrrhotite to be used to calculate both fS₂ and fO₂ for the experiments. At run temperature, (1042 ± 2 °C), a rhyolitic melt can be formed at low pressure, under nominally dry conditions, which removed the need for confining pressure as well as externally imposed fugacities. The silica-saturated melt allowed the charges to be contained in sealed evacuated silica tubes without danger of reaction, and with closed system behavior for molybdenum and tungsten.Experiments were run for durations up to 2000 min. Molybdenite (mb) and wolframite (wo) were added to the experiments as sources for molybdenum and tungsten, respectively. Manganese was added to the system as both a component of the starting rhyolitic pumice, and of Mn-bearing wolframite. Oxygen fugacity in these experiments was fixed at the Ni-NiO oxygen fugacity buffer. Sulfur fugacity was 10⁻¹ bar. Run products were analyzed by EPMA and LA-ICP-MS. Analysis of the run products yielded ( standard deviation of the mean): , , , and . The partition coefficients for manganese in this system are and .Simple Rayleigh fractionation modeling suggests that oxidized felsic melts produced through fractional crystallization may have lost as much as 14% of their initial molybdenum, but only 2% of their initial tungsten, through the removal of an Fe-S-O melt along with crystalline phases. Modeling consistent with conditions of oxygen and sulfur fugacity influenced by assimilation of sulfide (with low concentrations of molybdenum and tungsten) from, for example, sedimentary rock, results in evolved magmas significantly depleted in molybdenum, but only moderately depleted in tungsten. The molybdenum:tungsten ratio can vary by two orders of magnitude. These systematics may help to explain some of the variability in metal ratios of intrusion-related hydrothermal ore deposits.     

Oxygen fugacity dependence of Os solubility in haplobasaltic melt

Os equilibrium solubilities were determined at 1350 °C over a wide range of oxygen fugacities (−12 < log fO₂ < −7) applying the mechanically assisted equilibration technique (MAE) at 10⁵ Pa (= 1 bar). Os concentrations in the glass samples were analysed using ID-NTIMS. Additional LA-ICP-MS and SEM analyses were performed to detect, visualize and analyse the nature and chemistry of “nanonuggets.” Os solubilities determined range at a constant temperature of 1350 °C from 0.63 ± 0.04 to 37.4 ± 1.16 ppb depending on oxygen fugacity. At the highest oxygen fugacities, Os³⁺ can be confirmed as the main oxidation state of Os. At low oxygen fugacities (below log fO₂ = −8), samples are contaminated by nanonuggets which, despite the MAE technique, were still not removed entirely from the melt. However, the present results indicate that applying MAE technology does reduce the amount of nanonuggets present significantly, resulting in the lowest Os solubility results reported to date under these experimental conditions, and extending the experimentally accessible range of fO₂ for these studies to lower values. Calculated metal/silicate melt partition coefficients are therefore higher compared to previous studies, making Os more siderophile. Neglecting the as yet unknown temperature dependence of the Os metal/silicate melt partition coefficient, extrapolation of the obtained Os solubilities to conditions for core-mantle equilibrium, results in a , while metallic alloy/silicate melt partition coefficients range from 1.4 × 10⁶ to 8.6 × 10⁷, in agreement with earlier findings. Therefore remains too high by 2-4 orders of magnitude to explain the Os abundance in the Earth’s mantle as result of core-mantle equilibrium during core formation.     

Trace element partitioning between ilmenite, armalcolite and anhydrous silicate melt: Implications for the formation of lunar high-Ti mare basalts

We performed a series of experiments at high pressures and temperatures to determine the partitioning of a wide range of trace elements between ilmenite (Ilm), armalcolite (Arm) and anhydrous lunar silicate melt, to constrain geochemical models of the formation of titanium-rich melts in the Moon. Experiments were performed in graphite-lined platinum capsules at pressures and temperatures ranging from 1.1 to 2.3 GPa and 1300-1400 °C using a synthetic Ti-enriched Apollo ‘black glass’ composition in the CaO-FeO-MgO-Al₂O₃-TiO₂-SiO₂ system. Ilmenite-melt and armalcolite-melt partition coefficients (D) show highly incompatible values for the rare earth elements (REE) with the light REE more incompatible compared to the heavy REE ( 0.0020 ± 0.0010 to 0.069 ± 0.010 for ilmenite; 0.0048 ± 0.0023 to 0.041 ± 0.008 for armalcolite). D values for the high field strength elements vary from highly incompatible for Th, U and to a lesser extent W (for ilmenite: 0.0013 ± 0.0008, 0.0035 ± 0.0015 and 0.039 ± 0.005, and for armalcolite 0.008 ± 0.003, 0.0048 ± 0.0022 and 0.062 ± 0.03), to mildly incompatible for Nb, Ta, Zr, and Hf (e.g. 0.28 ± 0.05 and : 0.76 ± 0.07). Both minerals fractionate the high field strength elements with D_Ta/D_Nb and D_Hf/D_Zr between 1.3 and 1.6 for ilmenite and 1.3 and 1.4 for armalcolite. Armalcolite is slightly more efficient at fractionating Hf from W during lunar magma ocean crystallisation, with D_Hf/D_W = 12-13 compared to 6.7-7.5 for ilmenite. The transition metals vary from mildly incompatible to compatible, with the highest compatibilities for Cr in ilmenite (D ∼ 7.5) and V in armalcolite (D ∼ 8.1). D values show no clear variation with pressure in the small range covered.Crystal lattice strain modelling of D values for di-, tri- and tetravalent trace elements shows that in ilmenite, divalent elements prefer to substitute for Fe while armalcolite data suggest REE replacing Mg. Tetravalent cations appear to preferentially substitute for Ti in both minerals, with the exception of Th and U that likely substitute for the larger Fe or Mg cations. Crystal lattice strain modelling is also used to identify and correct for very small (∼0.3 wt.%) melt contamination of trace element concentration determinations in crystals.Our results are used to model the Lu-Hf-Ti concentrations of lunar high-Ti mare basalts. The combination of their subchondritic Lu/Hf ratios and high TiO₂ contents requires preferential dissolution of ilmenite or armalcolite from late-stage, lunar magma ocean cumulates into low-Ti partial melts of deeper pyroxene-rich cumulates.     

The partitioning behavior of silver in a vapor-brine-rhyolite melt assemblage

The partitioning of silver in a sulfur-free rhyolite melt-vapor-brine assemblage has been quantified at 800 °C, pressures of 100 and 140 MPa and f_O2≈NNO (nickel-nickel oxide). Silver solubility (±2σ) in rhyolite increases 5-fold from 105 ± 21 to 675 ± 98 μg/g as pressure increases from 100 to 140 MPa. Nernst-type partition coefficients describing the mass transfer of silver at 100 MPa between vapor and melt, brine and melt and vapor and brine are 32 ± 30, 1151 ± 238 and 0.026 ± 0.004, respectively. At 140 MPa, values for for vapor and melt, brine and melt, and vapor and brine are 32 ± 10, 413 ± 172 and 0.06 ± 0.03, respectively. Apparent equilibrium constant values (±2σ) describing the exchange of silver and sodium between vapor and melt, , at 100 and 140 MPa are 105 ± 68 and 14 ± 6. The average values (±2σ) for silver and sodium exchange between brine and melt, , at 100 and 140 MPa are 313 ± 288 and 65 ± 12. These data indicate that the mass transfer of silver from rhyolite melt to an exsolved volatile phase(s) is enhanced at 100 MPa relative to 140 MPa, suggesting that decompression increases the silver ore-generative potential of an evolving silicate magma. Model calculations using the new data suggest that the evolution of low-density, aqueous fluid (i.e., vapor) may be responsible for the the silver tonnage of many porphyry-type and perhaps epithermal-type ore deposits. For example, Halter et al. (Halter W. E., Pettke T. and Heinrich C. A. (2002) The origin of Cu/Au ratios in porphyry-type ore deposits. Science296, 1842-1844) used detailed silicate and sulfide melt inclusion and vapor and brine fluid inclusions analyses to estimate a melt volume on the order of 15 km³ to satisfy the copper budget at the Bajo de la Alumbrera copper-, gold-, silver-ore deposit. Using their melt volume estimate with the data presented here, model calculations for a 15-km³ felsic melt, saturated with pyrrhotite and magnetite, suggest that a low-salinity magmatic vapor may scavenge on the order of 7 × 10¹² g of silver from the melt. This quantity of silver exceeds the discovered 2 × 10⁹ g of Ag at Alumbrera. Calculated tonnages for numerous other deposits yield similar results. The excess silver in the vapor, remaining after porphyry formation, is then available to precipitate at lower PTconditions in the stratigraphically higher epithermal environment. These data suggest that silver, and perhaps other ore metals, in the porphyry-epithermal continuum may be derived solely from the time-integrated flux of dominantly low-salinity vapor exsolved from a series of sequential magma batches.     

Alkali exchange equilibria between a silicate melt and coexisting magmatic volatile phase: an experimental study at 800°C and 100 MPa

Many experimental studies have been performed to evaluate the composition of coexisting silicate melts and magmatic volatile phases (MVP). However, few studies have attempted to define the relationship between melt chemistry and the acidity of a chloride-bearing fluid. Here we report data on melt composition as a function of the HCl concentration of coexisting brines. We performed 35 experimental runs with a NaCl-KCl-HCl-H₂O brine (70 wt% NaCl [equivalent])-silicate melt (starting composition of Qtz_0.38Ab_0.33Or_0.29, anhydrous) assemblage at 800°C and 100 MPa. We determined an apparent equilibrium constant

     

Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle

The solubility of carbon in Fe and Fe-5.2 wt.% Ni melts, saturated with graphite, determined by electron microprobe analysis of quenched metal melts was 5.8 ± 0.1 wt.% at 2000 °C, 6.7 ± 0.2 wt.% at 2200 °C, and 7.4 ± 0.2 wt.% at 2410 °C at 2 GPa, conditions relevant for core/mantle differentiation in a shallow magma ocean. These solubilities are slightly lower than low-pressure literature values and significantly beneath calculated values for even higher pressures [e.g., Wood B. J. (1993) Carbon in the core. Earth Planet. Sci. Lett.117, 593-607]. The trend of C solubility versus temperature for Fe-5.2 wt.% Ni melt, within analytical uncertainties, is similar to or slightly lower (∼0.2-0.4 wt.%) than that of pure Fe. Carbon content of core melts and residual mantle silicates derived from equilibrium batch or fractional segregation of core liquids and their comparison with our solubility data and carbon content estimate of the present day mantle, respectively, constrain the partition coefficient of carbon between silicate and metallic melts, in a magma ocean. For the entire range of possible bulk Earth carbon content from chondritic to subchondritic values, of 10⁻⁴ to 1 is derived. But for ∼1000 ppm bulk Earth carbon, is between 10⁻² and 1. Using the complete range of possible for a magma ocean at ∼2200 °C, we predict maximum carbon content of the Earth’s core to be ∼6-7 wt.% and a preferred value of 0.25 ± 0.15 wt.% for a bulk Earth carbon concentration of ∼1000 ppm.     

Apatite solubility in carbonatitic liquids and trace element partitioning between apatite and carbonatite at high pressure

We have measured apatite solubility in calcic carbonatitic liquids and determined apatite/melt partition coefficients for a series of trace elements, including the rare earth elements (REE), high field strength elements (HFSE), Rb, Sr, U-Th-Pb. Experiments were performed between 4 and 6 GPa, from 1200 to 1380 °C, using the multianvil apparatus. Trace element concentrations were determined by laser ablation ICP-MS and electron microprobe. In addition, a specific protocol was designed to measure carbon concentration in the apatites, using the electron microprobe. Two starting apatite samples were used in order to test for the effect of apatite chemistry on partitioning behavior.Apatite solubility is lower in calcitic melts by a factor 3-5 compared to dolomitic melts (3-5.5 vs. 10-18 wt.% P₂O₅ in melt). We interpret this difference in terms of solubility product in the liquid and propose an empirical model for apatite saturation that takes into account melt calcium content. We conclude that calcitic melts that may form by melting of carbonated eclogites could be saturated with residual apatite, contrary to dolomitic melts formed in carbonated peridotites.Compatibility behavior of the REE depends on apatite silica content: REE are compatible in apatites containing 3.5-5 wt.% SiO₂, with values between 1.5 and 4, whereas REE are incompatible in apatites containing 0.2 wt.% SiO₂. HFSE, U, Th, and Y are compatible in silica-rich apatite, with while . Strontium is always retained in the melt, with of the order of 0.5. Lead appears to be incompatible in apatite, although this finding is weakened by almost complete Pb loss to sample container. High silica concentration favors REE incorporation in apatite by allowing for charged balanced coupled substitution. Sulfur and carbonate may also favor REE incorporation in apatite. Our results allow to reconcile previously published experimental determinations of REE partitioning. We use our experimentally determined partition coefficients to investigate the impact of residual apatite during partial melting of recycled carbonated material (eclogite + sediments) and discuss how the chemical characteristics of the produced liquids can be affected by residual apatite.     

The solution behavior of H2O in peralkaline aluminosilicate melts at high pressure with implications for properties of hydrous melts

Solubility and solution mechanisms of H₂O in depolymerized melts in the system Na₂O-Al₂O₃-SiO₂ were deduced from spectroscopic data of glasses quenched from melts at 1100 °C at 0.8-2.0 GPa. Data were obtained along a join with fixed nominal NBO/T = 0.5 of the anhydrous materials [Na₂Si₄O₉-Na₂(NaAl)₄O₉] with Al/(Al+Si) = 0.00-0.25. The H₂O solubility was fitted to the expression, X_H2O=0.20+0.0020f_H2O-0.7X_Al+0.9(X_Al)², where X_H2O is the mole fraction of H₂O (calculated with O = 1), f_H2O the fugacity of H₂O, and X_Al = Al/(Al+Si). Partial molar volume of H₂O in the melts, , calculated from the H₂O-solulbility data assuming ideal mixing of melt-H₂O solutions, is 12.5 cm³/mol for Al-free melts and decreases linearly to 8.9 cm³/mol for melts with Al/(Al+Si) ∼ 0.25. However, if recent suggestion that is composition-independent is applied to constrain activity-composition relations of the hydrous melts, the activity coefficient of H₂O, , increases with Al/(Al+Si).Solution mechanisms of H₂O were obtained by combining Raman and ²⁹Si NMR spectroscopic data. Degree of melt depolymerization, NBO/T, increases with H₂O content. The rate of NBO/T-change with H₂O is negatively correlated with H₂O and positively correlated with Al/(Al+Si). The main depolymerization reaction involves breakage of oxygen bridges in Q⁴-species to form Q² species. Steric hindrance appears to restrict bonding of H⁺ with nonbridging oxygen in Q³ species. The presence of Al³⁺ does not affect the water solution mechanisms significantly.     

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

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

Clinopyroxene dissolution in basaltic melt

The history of magmatic systems may be inferred from reactions between mantle xenoliths and host basalt if the thermodynamics and kinetics of the reactions are quantified. To study diffusive and convective clinopyroxene dissolution in silicate melts, diffusive clinopyroxene dissolution experiments were conducted at 0.47–1.90 GPa and 1509–1790 K in a piston-cylinder apparatus. Clinopyroxene saturation is found to be roughly determined by MgO and CaO content. The effective binary diffusivities, D_MgO and D_CaO, and the interface melt saturation condition, , are extracted from the experiments. D_MgO and D_CaO show Arrhenian dependence on temperature. The pressure dependence is small and not resolved within 0.47–1.90 GPa. in the interface melt increases with increasing temperature, but decreases with increasing pressure. Convective clinopyroxene dissolution, where the convection is driven by the density difference between the crystal and melt, is modeled using the diffusivities and interface melt saturation condition. Previous studies showed that the convective dissolution rate depends on the thermodynamics, kinetics and fluid dynamics of the system. Comparing our results for clinopyroxene dissolution to results from a previous study on convective olivine dissolution shows that the kinetic and fluid dynamic aspects of the two minerals are quite similar. However, the thermodynamics of clinopyroxene dissolution depends more strongly on the degree of superheating and composition of the host melt than that of olivine dissolution. The models for clinopyroxene and olivine dissolution are tested against literature experiments on mineral–melt interaction. They are then applied to previously proposed reactions between Hawaii basalts and mantle minerals, mid-ocean ridge basalts and mantle minerals, and xenoliths digestion in a basalt at Kuandian, Northeast China.     

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.     

Carbon dioxide in silicate melts: A molecular dynamics simulation study

We have performed a series of molecular dynamics simulations aimed at the evaluation of the solubility of CO₂ in silicate melts of natural composition (from felsic to ultramafic). In making in contact within the simulation cell a supercritical CO₂ phase with a silicate melt of a given composition, we have been able to evaluate (i) the solubility of CO₂ in the P-T range 1473-2273 K and 20-150 kbar, (ii) the density change experienced by the CO₂-bearing melt, (iii) the respective concentrations of CO₂ and species in the melt, (iv) the lifetime and the diffusivity of these species and (v) the structure of the melt around the carbonate groups. The main results are the following:(1) The solubility of CO₂ increases markedly with the pressure in the three investigated melts (a rhyolite, a mid-ocean ridge basalt and a kimberlite) from about ∼2 wt% CO₂ at 20 kbar to ∼25 wt% at 100 kbar and 2273 K. The solubility is found to be weakly dependent on the melt composition (as far as the present compositions are concerned) and it is only at very high pressure (above ∼100 kbar) that a clear hierarchy between solubilities occurs (rhyolite < MORB < kimberlite). Furthermore at a given pressure the calculated solubility is negatively correlated with the temperature.(2) In CO₂-saturated melts, the proportion of carbonate ions is positively correlated with the pressure at isothermal condition and is negatively correlated with the temperature at isobaric condition (and vice versa for molecular CO₂). Furthermore, at fixed (P, T) conditions the proportion of carbonate ions is higher in CO₂-undersaturated melts than in the CO₂-saturated melt. Although the proportion of molecular CO₂ decreases when the degree of depolymerization of the melt increases, it is still significant in CO₂-saturated basic and ultrabasic compositions at high temperatures. This finding is at variance with experimental data on CO₂-bearing glasses which show no evidence of molecular CO₂ as soon as the degree of depolymerization of the melt is high (e.g. basalt). These conflicting results can be reconciled with each other by noticing that a simple low temperature extrapolation of the simulation data predicts that the proportion of molecular CO₂ in basaltic melts might be negligible in the glass at room temperature.(3) The carbonate ions are found to be transient species in the liquid phase, with a lifetime increasing exponentially with the inverse of the temperature. Contrarily to a usual assumption, the diffusivity of carbonate ions into the liquid silicate is not vanishingly small with respect to that of CO₂ molecules: in MORB they differ from each other by a factor of ∼6 at 1473 K and only a factor of ∼2 at 2273 K. Although the bulk diffusivity of CO₂ is governed primarily by the diffusivity of CO₂ molecules, the carbonate ions contribute significantly to the diffusivity of CO₂ in depolymerized melts.(4) Concerning the structure of the CO₂-bearing silicate melt, the carbonate ions are found to be preferentially associated with NBO’s of the melt, with an affinity for NBOs which exceeds that for BOs by almost one order of magnitude. This result explains why the concentration in carbonate ions is positively correlated with the degree of depolymerization of the melt and diminishes drastically in fully polymerized melts where the number of NBO’s is close to zero. Furthermore, the network modifier cations are not randomly distributed in the close vicinity of carbonate groups but exhibit a preferential ordering which depends at once on the nature of the cation and on the melt composition. However at the high temperatures investigated here, there is no evidence of long lived complexes between carbonate groups and metal cations.     

The partitioning behavior of As and Au in S-free and S-bearing magmatic assemblages

The partitioning of As and Au between rhyolite melt and low-salinity vapor (2 wt% NaCl eq.) in a melt-vapor-Au metal ± magnetite ± pyrrhotite assemblage has been quantified at 800 °C, 120 MPa and f_O2=NNO. The S-bearing runs have calculated values for the fugacities of H₂S, SO₂ and S₂ of logf_H2S=1.1, logf_SO2=-1.5, and logf_S2=-3.0. The ratio of H₂S to SO₂ is on the order of 400. The experiments constrain the effect of S on the partitioning behavior of As and Au at magmatic conditions. Calculated average Nernst-type partition coefficients (±1σ) for As between vapor and melt, , are 1.0 ± 0.1 and 2.5 ± 0.3 in the S-free and S-bearing assemblages, respectively. These results suggest that sulfur has a small, but statistically meaningful, effect on the mass transfer of As between silicate melt and low-salinity vapor at the experimental conditions. Efficiencies of removal, calculated following Candela and Holland (1986), suggest that the S-free and S-bearing low-salinity vapor can scavenge approximately 41% and 63% As from water-saturated rhyolite melt, respectively, during devolatilization assuming that As is partitioned into magnetite and pyrrhotite during second boiling. The S-free data are consistent with the presence of arsenous acid, As(OH)₃ in the vapor phase. However, the S-bearing data suggest the presence of both arsenous acid and a As-S complex in S-bearing magmatic vapor. Apparent equilibrium constants, , describing the partitioning of As between melt and vapor are −1.3 (0.1) and −1.1 (0.1) for the S-free and S-bearing runs, respectively. The increase in the value of with the addition of S suggests a role for S in complexing and scavenging As from the melt during degassing.The calculated vapor/melt partition coefficients (±1σ) for Au between vapor and melt, , in S-free and S-bearing assemblages are 15 ± 2.5 and 12 ± 0.3, respectively. Efficiencies of removal (Candela and Holland, 1986) for the S-free melt, calculated assuming that magnetite is the dominant Au-sequestering solid phase during crystallization (Simon et al., 2003), suggest that magmatic vapor may scavenge on the order of 72% Au from a water-saturated melt. Efficiencies of removal calculated for the S-bearing assemblage, assuming pyrrhotite and magnetite are the dominant Au-sequestering solid phases, indicate that vapor may scavenge on the order of 60% Au from the melt. These model calculations suggest that the loss of pyrrhotite and magnetite from a melt, owing to punctuated differentiation during ascent and emplacement, does not prohibit the ability of a rhyolite melt to generate a large-tonnage Au deposit. Apparent equilibrium constants describing the partitioning of Au between melt and vapor were calculated using the mean values for the S-free and S-bearing assemblages; only S-bearing data from runs longer than 400 h were used as shorter runs may not have reached equilibrium with respect only to vapor/melt partitioning of Au. The values for are −4.4 (0.1) and −4.2 (0.2) for the S-free and S-bearing runs, respectively. These data suggest that the presence of S does not affect the mass transfer of Au from degassing silicate melt to an exsolved, low-salinity vapor in a low-f_S2 assemblage (i.e., pyrrhotite-magnetite at NNO) at the experimental conditions reported here. Efficiencies of removal are calculated and used to model the mass transfer of Au from a crystallizing silicate melt to an exsolved, low-salinity vapor phase. The calculations suggest that the model, absolute tonnage of Au scavenged and transported by S-free and S-bearing vapors, from a crystallizing melt, would be comparable and that the time-integrated flux of low-salinity vapor could be responsible for a significant quantity of the Au in magmatic-hydrothermal ore deposits.     

Copper partitioning in a melt-vapor-brine-magnetite-pyrrhotite assemblage

The effect of sulfur on the partitioning of Cu in a melt-vapor-brine ± magnetite ± pyrrhotite assemblage has been quantified at 800 °C, 140 MPa, f_O2 = nickel-nickel oxide (NNO), logf_S2=-3.0 (i.e., on the magnetite-pyrrhotite curve at NNO), logf_H2S=-1.3 and logf_SO2=-1. All experiments were vapor + brine saturated. Vapor and brine fluid inclusions were trapped in silicate glass and self-healed quartz fractures. Vapor and brine are dominated by NaCl, KCl and HCl in the S-free runs and NaCl, KCl and FeCl₂ in S-bearing runs. Pyrrhotite served as the source of sulfur in S-bearing experiments. The composition of fluid inclusions, glass and crystals were quantified by laser-ablation inductively coupled plasma mass spectrometry. Major element, chlorine and sulfur concentrations in glass were quantified by using electron probe microanalysis. Calculated Nernst-type partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-free system. The partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-bearing system. Apparent equilibrium constants (±1σ) describing Cu and Na exchange between vapor and melt and brine and melt were also calculated. The values of are 34 ± 21 and 128 ± 29 in the S-free and S-bearing runs, respectively. The values of are 33 ± 22 and60 ± 5 in the S-free and S-bearing runs, respectively. The data presented here indicate that the presence of sulfur increases the mass transfer of Cu into vapor from silicate melt. Further, the nearly threefold increase in suggests that Cu may be transported as both a chloride and sulfide complex in magmatic vapor, in agreement with hypotheses based on data from natural systems. Most significantly, the data demonstrate that the presence of sulfur enhances the partitioning of Cu from melt into magmatic volatile phases.     

Solution behavior of reduced COH volatiles in silicate melts at high pressure and temperature

The solubility and solution mechanisms of reduced COH volatiles in Na₂OSiO₂ melts in equilibrium with a (H₂ + CH₄) fluid at the hydrogen fugacity defined by the iron-wüstite + H₂O buffer [f_H2(IW)] have been determined as a function of pressure (1-2.5 GPa) and silicate melt polymerization (NBO/Si: nonbridging oxygen per silicon) at 1400 °C. The solubility, calculated as CH₄, increases from ∼0.2 wt% to ∼0.5 wt% in the melt NBO/Si-range ∼0.4 to ∼1.0. The solubility is not significantly pressure-dependent, probably because f_H2(IW) in the 1-2.5 GPa range does not vary greatly with pressure. Carbon isotope fractionation between methane-saturated melts and (H₂ + CH₄) fluid varied by ∼14‰ in the NBO/Si-range of these melts.The (C..H) and (O..H) speciation in the quenched melts was determined with Raman and ¹H MAS NMR spectroscopy. The dominant (C..H)-bearing complexes are molecular methane, CH₄, and a complex or functional group that includes entities with CCH bonding. Minor abundance of complexes that include SiOCH₃ bonding is tentatively identified in some melts. There is no spectroscopic evidence for SiC or SiCH₃. Raman spectra indicate silicate melt depolymerization (increasing NBO/Si). The [CH₄/CCH]^melt abundance ratio is positively correlated with NBO/Si, which is interpreted to suggest that the (CCH)-containing structural entity is bonded to the silicate melt network structure via its nonbridging oxygen. The ∼14‰ carbon isotope fractionation change between fluid and melt is because of the speciation changes of carbon in the melt.     

Experimental studies of metal-silicate partitioning of Sb: Implications for the terrestrial and lunar mantles

The terrestrial mantle has a well defined Sb depletion of ∼7 ± 1 (Jochum and Hofmann, 1997), and the lunar mantle is depleted relative to the Earth by a factor of ∼50 ± 5 (Wolf and Anders, 1980). Despite these well defined depletions, there are few data upon which to evaluate their origin—whether due to volatility or core formation. We have carried out a series of experiments to isolate several variables such as oxygen fugacity, temperature, pressure, and silicate and metallic melt compositions, on the magnitude of . The activity of Sb in FeNi metal is strongly composition dependent such that solubility of Sb as a function of fO₂ must be corrected for the metal composition. When the correction is applied, Sb solubility is consistent with 3+ valence. Temperature series (at 1.5 GPa) shows that decreases by a factor of 100 over 400 °C, and a pressure series exhibits an additional decrease between ambient pressure (100 MPa) and 13 GPa. A strong dependence upon silicate melt composition is evident from a factor of 100 decrease in between nbo/t values of 0.3 and 1.7. Consideration of all these variables indicates that the small Sb depletion for the Earth’s mantle can be explained by high PT equilibrium partitioning between metal and silicate melt . The relatively large lunar Sb depletion can also be explained by segregation of a small metallic core, at lower pressure conditions where is much higher (2500).     

Core formation and the oxidation state of the Earth: Additional constraints from Nb, V and Cr partitioning

We have combined metal-silicate partitioning data from the literature with new experimental results at 1.5-8 GPa and 1480-2000 °C to parameterize the effects of pressure, temperature and composition on the partitioning of V, Cr and Nb between liquid Fe metal (with low S and C content) and silicate melt.Using information from the steelmaking literature to correct for interactions in the metal phase, we find that, for peridotitic silicate melts, metal-silicate partition coefficients are given by: