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.     

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.     

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.     

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.     

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.     

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.     

Platinum solubility and partitioning in a felsic melt-vapor-brine assemblage

The solubility and partitioning of Pt in a S-free vapor - brine - rhyolite melt - Pt metal assemblage has been quantified at 800 °C, f_O2=NNO and pressures of 100 and 140 MPa. Vapor and brine were sampled at run conditions by trapping these phases as glass-hosted fluid inclusions as the melt cooled through the glass transition temperature. The vapor and brine were in equilibrium with the melt at the time of trapping and, thus, represent fluids which were sampled at the termination of each experimental run. The microthermometrically determined salinities of vapor and brine are ∼2 and ∼63 wt.% NaCl eq. and ∼9 and ∼43 wt.% NaCl eq. at 100 and 140 MPa, respectively. Platinum solubilities in vapor, brine and glass (i.e., quenched melt) were quantified by using laser ablation - inductively coupled plasma - mass spectrometry (LA-ICP-MS). Equilibrium is discussed with reference to the major and trace element concentrations of glass-hosted fluid inclusions as well as the silicate melt over run times that varied from 110 to 377 h at 140 MPa and 159 to 564 h at 100 MPa. Platinum solubility values (±1σ) in H₂O-saturated felsic melt are 0.28 ± 0.13 μg/g and 0.38 ± 0.06 μg/g at 140 and 100 MPa, respectively. Platinum solubility values () at 140 and 100 MPa, respectively, in aqueous vapor are 0.91 ± 0.29 μg/g and 0.37 ± 0.17 μg/g and in are brine 16 ± 10 μg/g and 3.3 ± 1.0 μg/g. The measured solubility data were used to calculate Nernst-type partition coefficients for Pt between vapor/melt, brine/melt and vapor/brine. The partition coefficient values () for vapor/melt, brine/melt and vapor/brine at 140 MPa are 2.9 ± 1.0, 67 ± 27, and 0.13 ± 0.05 and at 100 MPa are 1.0 ± 0.2, 6.8 ± 2.4, and 0.15 ± 0.05. The partitioning data were used to model the Pt-scavenging capacity of vapor and brine during the crystallization-driven degassing (i.e., second boiling) of a felsic silicate melt over a depth range (i.e., 3-6 km) consistent with the evolution of magmatic-hydrothermal ore deposits. Model calculations suggest that aqueous vapor and brine can scavenge sufficient quantities of Pt, and by analogy other platinum group elements (PGE), to produce economically important PGE-rich magmatic-hydrothermal ore deposits in Earth’s upper continental crust.     

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.     

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).     

Gold and copper partitioning in magmatic-hydrothermal systems at 800 °C and 100 MPa

Porphyry-type ore deposits sometimes contain fluid inclusion compositions consistent with the partitioning of copper and gold into vapor relative to coexisting brine at the depositional stage. However, this has not been reproduced experimentally at magmatic conditions. In an attempt to determine the conditions under which copper and gold may partition preferentially into vapor relative to brine at temperatures above the solidus of granitic magmas, we performed experiments at 800 °C, 100 MPa, oxygen fugacity () buffered by Ni-NiO, and fixed at either 3.5 × 10⁻² by using intermediate solid solution-pyrrhotite, or 1.2 × 10⁻⁴ by using intermediate solid solution-pyrrhotite-bornite. The coexisting vapor (∼3 wt.% NaCl eq.) and brine (∼68 wt.% NaCl eq.) were composed initially of NaCl + KCl + HCl + H₂O, with starting HCl set to <1000 μg/g in the aqueous mixture. Synthetic vapor and brine fluid inclusions were trapped at run conditions and subsequently analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Our experiments demonstrate that copper and gold partitioned strongly into the magmatic volatile phase(s) (MVP) (i.e., vapor or brine) relative to a silicate melt over the entire imposed range of . Nernst style partition coefficients between coexisting brine (b) and melt (m), D^b/m (±1σ), range from 3.6(±2.2) × 10¹ to 4(±2) × 10² for copper and from 1.2(±0.6) × 10² to 2.4(±2.4) × 10³ for gold. Partition coefficients between coexisting vapor (v) and melt, D^v/m range from 2.1 ± 0.7 to 18 ± 5 and 7(±3) × 10¹ to 1.6(±1.6) × 10² for copper and gold, respectively. Partition coefficients for all experiments between coexisting brine and vapor, D^b/v (±1σ), range from 7(±2) to 1.0(±0.4) × 10² and 1.7(±0.2) to 15(±2) for copper and gold, respectively. Observed average D^b/v at an of 1.2 × 10⁻⁴ were elevated, 95(±5) and 15 ± 1 for copper and gold, respectively, relative to those at the higher of 3.5 × 10⁻² where D^b/v were 10(±5) for copper and 7(±6) for gold. Thus, there is an inverse relationship between the and the D^b/v for both copper and gold with increasing resulting in a decrease in the D^b/v signifying increased importance of the vapor phase for copper and gold transport. This suggests that copper and gold may complex with volatile S-species as well as Cl-species at magmatic conditions, however, none of the experiments of our study at 800 °C and 100 MPa had a D^b/v ? 1. We did not directly determine speciation, but infer the existence of some metal-sulfur complexes based on the reported data. We suggest that copper and gold partition preferentially into the brine in most instances at or above the wet solidus. However, in most systems, the mass of vapor is greater than the mass of brine, and vapor transport of copper and gold may become more important in the magmatic environment at higher , lower , or near the critical point in a salt-water system. A D^b/v ? 1 at subsolidus hydrothermal conditions may also occur in response to changes in temperature, , , and/or acidity.Additionally, both copper and gold were observed to partition into intermediate solid solution and bornite much more strongly than into vapor, brine or silicate melt. This suggests that, although vapor and brine are both efficient at removing copper and gold from a silicate melt, the presence of Cu-Fe sulfides can sequester a substantial portion of the copper and gold contained within a silicate melt if the Cu-Fe sulfides are abundant.     

Highly siderophile element composition of the Earth’s primitive upper mantle: Constraints from new data on peridotite massifs and xenoliths

Osmium, Ru, Ir, Pt, Pd and Re abundances and ¹⁸⁷Os/¹⁸⁸Os data on peridotites were determined using improved analytical techniques in order to precisely constrain the highly siderophile element (HSE) composition of fertile lherzolites and to provide an updated estimate of HSE composition of the primitive upper mantle (PUM). The new data are used to better constrain the origin of the HSE excess in Earth’s mantle. Samples include lherzolite and harzburgite xenoliths from Archean and post-Archean continental lithosphere, peridotites from ultramafic massifs, ophiolites and other samples of oceanic mantle such as abyssal peridotites. Osmium, Ru and Ir abundances in the peridotite data set do not correlate with moderately incompatible melt extraction indicators such as Al₂O₃. Os/Ir is chondritic in most samples, while Ru/Ir, with few exceptions, is ca. 30% higher than in chondrites. Both ratios are constant over a wide range of Al₂O₃ contents, but show stronger scatter in depleted harzburgites. Platinum, Pd and Re abundances, their ratios with Ir, Os and Ru, and the ¹⁸⁷Os/¹⁸⁸Os ratio (a proxy for Re/Os) show positive correlations with Al₂O₃, indicating incompatible behavior of Pt, Pd and Re during mantle melting. The empirical sequence of peridotite-melt partition coefficients of Re, Pd and Pt as derived from peridotites () is consistent with previous data on natural samples. Some harzburgites and depleted lherzolites have been affected by secondary igneous processes such as silicate melt percolation, as indicated by U-shaped patterns of incompatible HSE, high ¹⁸⁷Os/¹⁸⁸Os, and scatter off the correlations defined by incompatible HSE and Al₂O₃. The bulk rock HSE content, chondritic Os/Ir, and chondritic to subchondritic Pt/Ir, Re/Os, Pt/Re and Re/Pd of many lherzolites of the present study are consistent with depletion by melting, and possibly solid state mixing processes in the convecting mantle, involving recycled oceanic lithosphere. Based on fertile lherzolite compositions, we infer that PUM is characterized by a mean Ir abundance of 3.5 ± 0.4 ng/g (or 0.0080 ± 0.0009*CI chondrites), chondritic ratios involving Os, Ir, Pt and Re (Os/Ir_PUM of 1.12 ± 0.09, Pt/Ir_PUM = 2.21 ± 0.21, Re/Os_PUM = 0.090 ± 0.002) and suprachondritic ratios involving Ru and Pd (Ru/Ir_PUM = 2.03 ± 0.12, Pd/Ir_PUM = 2.06 ± 0.31, uncertainties 1σ). The combination of chondritic and modestly suprachondritic HSE ratios of PUM cannot be explained by any single planetary fractionation process. Comparison with HSE patterns of chondrites shows that no known chondrite group perfectly matches the PUM composition. Similar HSE patterns, however, were found in Apollo 17 impact melt rocks from the Serenitatis impact basin [Norman M.D., Bennett V.C., Ryder G., 2002. Targeting the impactors: siderophile element signatures of lunar impact melts from Serenitatis. Earth Planet. Sci. Lett, 217-228.], which represent mixtures of chondritic material, and a component that may be either of meteoritic or indigenous origin. The similarities between the HSE composition of PUM and the bulk composition of lunar breccias establish a connection between the late accretion history of the lunar surface and the HSE composition of the Earth’s mantle. Although late accretion following core formation is still the most viable explanation for the HSE abundances in the Earth’s mantle, the “late veneer” hypothesis may require some modification in light of the unique PUM composition.     

Solubility of Pt in sulphide mattes: Implications for the genesis of PGE-rich horizons in layered intrusions

The partitioning of Pt in sulphide melt (matte) has been studied as a function of fS₂ and fO₂ at 1200 and 1300 °C. The results show that the solubility of Pt in mattes increases strongly with increasing fS₂ and decreases weakly with increasing fO₂. The increase in Pt solubility with increasing fS₂ is attributed to Pt dissolving in the melt as a sulphide species and the weak inverse dependence of Pt solubility on fO₂ to the diluting effect of increasing O in the melt at high fO₂. These results, coupled with measurements of Pt solubility in silicate melts taken from the literature, allow the calculation of Pt matte/silicate-melt partition coefficients () for a range of conditions pertinent to the formation of Pt-rich horizons in layered intrusions. The calculated values range between 10⁷ and 10¹¹, depending on fO₂ and fS₂, several orders of magnitude higher than previously published values. Our preferred value for for conditions appropriate to the Merensky Reef is 10⁷ and for the Stillwater Pt-rich horizon 10⁸. The new results are consistent with the magmatic hypothesis for Pt-rich horizons in layered intrusions.     

Temperature dependence of Pt and Rh solubilities in a haplobasaltic melt

The temperature dependence of the solubilities of Pt and Rh in a haplobasaltic (anorthite-diopside 1-bar eutectic composition) melt has been investigated at 1 bar and 1300 to 1550°C using the mechanically assisted equilibration technique (Dingwell et al., 1994). The experiments were performed at almost constant oxygen fugacity (log fO₂ = −2.5 ± 0.3) over the entire temperature range. Major element concentrations in the quenched glass samples were determined using an electron microprobe. Pt and Rh concentrations were obtained by laser ablation inductive coupled plasma mass spectrometry. From our data, we obtain the following expressions for the solubilities of pure Pt and pure Rh in anorthite-diopside eutectic melt at 1 bar and log fO₂ = −2.5:

     

Solubility and solution mechanisms of C-O-H volatiles in silicate melt with variable redox conditions and melt composition at upper mantle temperatures and pressures

Solubility and solution mechanisms in silicate melts of oxidized and reduced C-bearing species in the C-O-H system have been determined experimentally at 1.5 GPa and 1400 °C with mass spectrometric, NMR, and Raman spectroscopic methods. The hydrogen fugacity, f_H2, was controlled in the range between that of the iron-wüstite-H₂O (IW) and the magnetite-hematite-H₂O (MH) buffers. The melt polymerization varied between those typical of tholeiitic and andesitic melts.The solubility of oxidized (on the order of 1-2 wt% as C) and reduced carbon (on the order of 0.15-0.35 wt% as C) is positively correlated with the NBO/Si (nonbridging oxygen per silicon) of the melt. At given NBO/Si-value, the solubility of oxidized carbon is 2-4 times greater than under reducing conditions. Oxidized carbon dioxide is dissolved as complexes, whereas the dominant reduced species in melts are CH₃-groups forming bonds with Si⁴⁺ together with molecular CH₄. Formation of complexes results in silicate melt polymerization (decreasing NBO/Si), whereas solution of reduced carbon results in depolymerization of melts (increasing NBO/Si).Redox melting in the Earth’s interior has been explained with the aid of the different solution mechanisms of oxidized and reduced carbon in silicate melts. Further, effects of oxidized and reduced carbon on melt viscosity and on element partitioning between melts and minerals have been evaluated from relationships between melt polymerization and dissolved carbon combined with existing experimental data that link melt properties and melt polymerization. With total carbon contents in the melts on the order of several mol%, mineral/melt element partition coefficients and melt viscosity can change by several tens to several hundred percent with variable redox conditions in the range of the Earth’s deep crust and upper mantle.     

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.     

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

Experimental constraints on the mobility of Rhenium in silicate liquids

The volatization of Rhenium (Re) from melts of natural basalt, dacite and a synthetic composition in the CaO-MgO-Al₂O₃-SiO₂ system has been investigated at 0.1 MPa and 1250-1350 °C over a range of fO₂ conditions from log fO₂ = −10 to −0.68. Experiments were conducted using open top Pt crucibles doped with Re and Yb. Analysis of quenched glasses by laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) normal to the melt/gas interface showed concentration profiles for Re, to which a semi-infinite one-dimensional diffusion model could be applied to extract diffusion coefficients (D). The results show Re diffusivity in basalt at 1300 °C in air is log D_Re = −7.2 ± 0.3 cm²/s and increases to log D_Re = −6.6 ± 0.3 cm²/s when trace amounts of Cl were added to the starting material. At fO₂ conditions below the nickel-nickel oxide (NNO) buffer Re diffusivity decreases to and to in dacitic melt. In the CMAS composition, . The diffusivity of Re is comparable to Ar and CO₂ in basalt at 500 MPa favoring its release as a volatile. Our results support the contention that subaerial degassing is the cause of lower Re concentrations in arc-type and ocean island basalts compared to mid-ocean ridge basalts.     

Thermodynamic properties of H4SiO4 in the ideal gas state as evaluated from experimental data

Solid phases of silicon dioxide react with water vapor with the formation of hydroxides and oxyhydroxides of silica. Recent transpiration and mass-spectrometric studies convincingly demonstrate that H₄SiO₄ is the predominant form of silica in vapor phase at water pressure in excess of 10⁻² MPa. Available literature transpiration and solubility data for the reactions of solid SiO₂ phases and low-density water, extending from 424 to 1661 K, are employed for the determination of Δ_fG⁰, Δ_fH⁰ and S⁰ of H₄SiO₄ in the ideal gas state at 298.15 K, 0.1 MPa. In total, there are 102 data points from seven literature sources. The resulting values of the thermodynamic functions of H₄SiO₄(g) are: Δ_fG⁰ = −1238.51 ± 3.0 kJ mol⁻¹, Δ_fH⁰ = −1340.68 ± 3.5 kJ mol⁻¹ and S⁰ = 347.78 ± 6.2 J K⁻¹ mol⁻¹. These values agree quantitatively with one set of ab initio calculations. The relatively large uncertainties are mainly due to conflicting data for H₄SiO₄(g) from various sources, and new determinations of would be helpful. The thermodynamic properties of this species, H₄SiO₄(g), are necessary for realistic modeling of silica transport in a low-density water phase. Applications of this analysis may include the processes of silicates condensation in the primordial solar nebula, the precipitation of silica in steam-rich geothermal systems and the corrosion of SiO₂-containing alloys and ceramics in moist environments.     

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.     

Experimental partitioning of uranium between liquid iron sulfide and liquid silicate: Implications for radioactivity in the Earth’s core

Measurable uranium (U) is found in metal sulfide liquids in equilibrium with molten silicate at conditions appropriate for a planetary magma ocean: 1-10 GPa, 1750-2100 °C, 0-28 wt% S, and fO₂ 2 log units below IW. However, the transfer of U from metal sulfide to silicate under our experimental conditions is so complete that insufficient U would remain so as to be of any importance to the core’s heat budget. U content in the sulfide phase increases strongly with S content but shows no significant variability with either pressure or temperature. Maximum is 0.001 while most values are considerably lower.