High-pressure melting relations in Fe-C-S systems: Implications for formation, evolution, and structure of metallic cores in planetary bodies

We present new high-pressure temperature experiments on melting phase relations of Fe-C-S systems with applications to metallic core formation in planetary interiors. Experiments were performed on Fe-5 wt% C-5 wt% S and Fe-5 wt% C-15 wt% S at 2-6 GPa and 1050-2000 °C in MgO capsules and on Fe-13 wt% S, Fe-5 wt% S, and Fe-1.4 wt% S at 2 GPa and 1600 °C in graphite capsules. Our experiments show that: (a) At a given P-T, the solubility of carbon in iron-rich metallic melt decreases modestly with increasing sulfur content and at sufficiently high concentration, the interaction between carbon and sulfur can cause formation of two immiscible melts, one rich in Fe-carbide and the other rich in Fe-sulfide. (b) The mutual solubility of carbon and sulfur increases with increasing pressure and no super-liquidus immiscibility in Fe-rich compositions is likely expected at pressures greater than 5-6 GPa even for bulk compositions that are volatile-rich. (c) The liquidus temperature in the Fe-C-S ternary is significantly different compared to the binary liquidus in the Fe-C and Fe-S systems. At 6 GPa, the liquidus of Fe-5 wt% C-5 wt% S is 150-200 °C lower than the Fe-5 wt% S. (d) For Fe-C-S bulk compositions with modest concentration of carbon, the sole liquidus phase is iron carbide, Fe₃C at 2 GPa and Fe₇C₃ at 6 GPa and metallic iron crystallizes only with further cooling as sulfur is concentrated in the late crystallizing liquid. Our results suggest that for carbon and sulfur-rich core compositions, immiscibility induced core stratification can be expected for planets with core pressure less than ∼6 GPa. Thus planetary bodies in the outer solar system such as Ganymede, Europa, and Io with present day core-mantle boundary (CMB) pressures of ∼8, ∼5, and 7 GPa, respectively, if sufficiently volatile-rich, may either have a stratified core or may have experienced core stratification owing to liquid immiscibility at some stage of their accretion. A similar argument can be made for terrestrial planetary bodies such as Mercury and Earth’s Moon, but no such stratification is predicted for cores of terrestrial planets such as Earth, Venus, and Mars with the present day core pressure in the order ?136 GPa, ?100 GPa, and ?23 GPa. (e) Owing to different expected densities of Fe-rich (and carbon-bearing) and sulfur-rich metallic melts, their settling velocities are likely different; thus core formation in terrestrial planets may involve rain of more than one metallic melt through silicate magma ocean. (f) For small planetary bodies that have core pressures <6 GPa and have a molten core or outer core, settling of denser carbide-rich liquid or flotation of lighter, sulfide-rich melt may contribute to an early, short-lived geodynamo.     

Trace element partitioning between apatite and silicate melts

We present new experimental apatite/melt trace element partition coefficients for a large number of trace elements (Cs, Rb, Ba, La, Ce, Pr, Sm, Gd, Lu, Y, Sr, Zr, Hf, Nb, Ta, U, Pb, and Th). The experiments were conducted at pressures of 1.0 GPa and temperatures of 1250 °C. The rare earth elements (La, Ce, Pr, Sm, Gd, and Lu), Y, and Sr are compatible in apatite, whereas the larger lithophile elements (Cs, Rb, and Ba) are strongly incompatible. Other trace elements such as U, Th, and Pb have partition coefficients close to unity. In all experiments we found D_Hf > D_Zr, D_Ta ≈ D_Nb, and D_Ba > D_Rb > D_Cs. The experiments reveal a strong influence of melt composition on REE partition coefficients. With increasing polymerisation of the melt, apatite/melt partition coefficients for the rare earth elements increase for about an order of magnitude. We also present some results in fluorine-rich and water-rich systems, respectively, but no significant influence of either H₂O or F on the partitioning was found. Furthermore, we also present experimentally determined partition coefficients in close-to natural compositions which should be directly applicable to magmatic processes.     

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

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

The partitioning of Pt-Re-Os between solid and liquid metal in the Fe-Ni-S system at high pressure: Implications for inner core fractionation

Coupled ¹⁸⁶Os/¹⁸⁸Os and ¹⁸⁷Os/¹⁸⁸Os enrichments of plume-derived lavas have been suggested to reflect contributions of materials from the outer core (Brandon et al., 1998). This hypothesis is based on the assumption that the Earth’s liquid outer core has high Pt/Os and slightly high Re/Os ratios as a result of the crystallization of the solid inner core, and shows coupled enrichments in the ¹⁸⁶Os/¹⁸⁸Os and ¹⁸⁷Os/¹⁸⁸Os ratios, reflecting the decay of ¹⁹⁰Pt and ¹⁸⁷Re to ¹⁸⁶Os and ¹⁸⁷Os, respectively. Partitioning experiments of Pt-Re-Os between solid and liquid metal were performed at 5-20 GPa and 1250-1400 °C, to examine the effects of pressure in the Fe-Ni-S system. The ratios (D_Os/D_Pt, D_Os/D_Re) of measured partition coefficients of Pt, Re and Os are almost constant with increasing pressure. D_Os/D_Pt increases significantly, whereas D_Os/D_Re decreases, with increasing sulphur content in the liquid metal. On the basis of the present experimental results, it is unlikely that the required Pt-Re-Os fractionation is generated during inner core crystallization, assuming that the light element in the Earth’s core is sulphur.     

Trace element partitioning in the Fe-S-C system and its implications for planetary differentiation and the thermal history of ureilites

Element partitioning in metal-light element systems is important to our understanding of planetary differentiation processes. In this study, solid-metal/liquid-sulfide, liquid-metal/liquid-sulfide and solid-metal/troilite partition coefficients (D) were determined for 18 elements (Ag, As, Au, Co, Cr, Cu, Ge, Ir, Ni, Os, Pd, Pt, Mo, Mn, Re, Ru, Se and W) in the graphite-saturated Fe-S-C system at 1 atm. Compared at the same liquid S concentration, the solid/liquid partition coefficients are similar to those in the Fe-S system, but there are systematic differences that appear to be related to interactions with carbon dissolved in the solid metal. Elements previously shown to be “anthracophile” generally have larger solid/liquid partition coefficients in the Fe-S-C system, whereas those that are not have similar or smaller partition coefficients in the Fe-S-C system. The partitioning of trace elements between C-rich and S-rich liquids is, in most cases, broadly similar to the partitioning between solid metal and S-rich liquid. The highly siderophile elements Os, Re, Ir and W are partitioned strongly into the C-rich liquid, with D ? 100. The partition coefficients for Pt, Ge and W decrease significantly at the transition to liquid immiscibility, while the partition coefficient for Mo increases sharply. The bulk siderophile element patterns of ureilite meteorities appear to be better explained by separation of S-rich liquid from residual C-rich metallic liquid at temperatures above the silicate solidus, rather than by separation of S-rich liquid from residual solid metal at lower temperatures.     

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

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

The influence of carbon on trace element partitioning behavior

Carbon has been proposed as a potential light element in planetary cores, included in models of planetary core formation, and found in meteoritic samples and minerals. To better understand the effect of C on the partitioning behavior of elements, solid/liquid partition coefficients (D = (solid metal)/(liquid metal)) were determined for 17 elements (As, Au, Co, Cr, Cu, Ga, Ge, Ir, Ni, Os, Pd, Pt, Re, Ru, Sb, Sn, and W) over a range of C contents in the Fe-Ni-C system at 1 atm. The partition coefficients for the majority of the elements increased as the C content of the liquid increased, an effect analogous to that of S for many of the elements. In contrast, three of the elements, Cr, Re, and W, were found to have anthracophile (C-loving) preferences, partitioning more strongly into the metallic liquid as the C content increased, resulting in decreases to their partition coefficients. For half of the elements examined, the prediction that partitioning in the Fe-Ni-S and Fe-Ni-C systems could be parameterized using a single set of variables was not supported. The effects of S and C on elemental partitioning behavior can be quite different; consequently, the presence of different non-metals can result in different fractionation patterns, and that uniqueness offers the opportunity to gain insight into the evolution of planetary bodies.     

The Fe-rich liquidus in the Fe-FeS system from 1 bar to 10 GPa

The composition and evolution of a metallic planetary core is determined by the behavior with pressure of the eutectic and the liquidus on the Fe-rich side of the Fe-FeS eutectic. New experiments at 6 GPa presented here, along with existing experimental data, inform a thermodynamic model for this liquidus from 1 bar to at least 10 GPa. Fe-FeS has a eutectic that becomes more Fe-rich but remains constant in T up to 6 GPa. The 1 bar, 3 GPa, and 6 GPa liquidi all cross at a pivot point at 1640 ± 5 K and FeS_37 ± 0.5. This liquid/crystalline metal equilibrium is T-x-fixed and pressure independent through 6 GPa. Models of the 1 bar through 10 GPa experimental liquidi show that with increasing P there is an increase in the T separation between the liquidus and the crest of the metastable two-liquid solvus. The solvus crest decreases in T with increasing P. The model accurately reproduces all the experimental liquidi from 1 bar to 10 GPa, as well as reproducing the 0-6 GPa pivot point. The 14 GPa experimental liquidus ( [Chen et al., 2008a] and Chen et al., 2008b) deviates sharply from the lower pressure trends indicating that the 0-10 GPa model no longer applies to this 14 GPa data.     

Experimental comparison of trace element partitioning between clinopyroxene and melt in carbonate and silicate systems, and implications for mantle metasomatism

Experiments in the systems diopside-albite (Di-Ab) and diopside-albite-dolomite (Di-Ab-Dmt), doped with a wide range of trace elements, have been used to characterise the difference between clinopyroxene-silicate melt and clinopyroxene-carbonate melt partitioning. Experiments in Di-Ab-Dmt yielded clinopyroxene and olivine in equilibrium with CO₂-saturated dolomitic carbonate melt at 3 GPa, 1375 °C. The experiments in Di-Ab were designed to bracket those conditions (3 GPa, 1640 °C and 0.8 GPa, 1375 °C), and so minimise the contribution of differential temperature and pressure to partitioning. Partition coefficients, determined by SIMS analysis of run products, differ markedly for some elements between Di-Ab and Di-Ab-Dmt systems. Notably, in the carbonate system clinopyroxene-melt partition coefficients for Si, Al, Ga, heavy REE, Ti and Zr are higher by factors of 5 to 200 than in the silicate system. Conversely, partition coefficients for Nb, light REE, alkali metals and alkaline earths show much less fractionation (<3). The observed differences compare quantitatively with experimental data on partitioning between immiscible carbonate and silicate melts, indicating that changes in melt chemistry provide the dominant control on variation in partition coefficients in this case. The importance of melt chemistry in controlling several aspects of element partitioning is discussed in light of the energetics of the partitioning process. The compositions of clinopyroxene and carbonate melt in our experiments closely match those of near-solidus melts and crystals in CMAS-CO₂ at 3 GPa, suggesting that our partition coefficients have direct relevance to melting of carbonated mantle lherzolite. Melts so produced will be characterised by elevated incompatible trace element concentrations, due to the low degrees of melting involved, but marked depletions of Ti and Zr, and fractionated REE patterns. These are common features of natural carbonatites. The different behaviour of trace elements in carbonate and silicate systems will lead to contrasted styles of trace element metasomatism in the mantle. Received: 15 July 1999 / Accepted: 18 February 2000     

Metal-silicate partitioning of cesium: Implications for core formation

Here we present the first set of metal-silicate partitioning data for Cs, which we use to examine whether the primitive mantle depletion of Cs can be attributed to core segregation. Our experiments independently varied pressure from 5 to 15 GPa, temperature from 1900 to 2400 °C, metallic sulfur content from pure Fe to pure FeS, silicate melt polymerization, expressed as a ratio of non-bridging oxygens to tetrahedrally coordinated cations (nbo/t) from 1.26 to 3.1, and fO₂ from two to four log units below the iron-wüstite buffer. The most important controls on the partitioning behavior of alkalis were the metallic sulfur content, expressed as X_S, and the nbo/t of the silicate liquid. Normalization of X_S to 0.5 yielded the following expressions for D-values as a function of nbo/t: log D_Na = −2.0 + 0.44 × (nbo/t), log D_K = −2.4 + 0.67 × ( nbo/t), and log D_Cs = −3.2 + 1.17 × (nbo/t). Normalization of nbo/t to 2.7 resulted in the following equations for D-values as a function of S content: log D_Na = −4.1 + 6.4 × X_S, log D_K = −7.7 + 13.9 × X_S, and log D_Cs = −12.1 + 23.3 × X_S.There appears to be a negative pressure effect up to 15 GPa, but it should be noted that this trend was not present before normalization, and is based on only two measurements. There is a positive trend in cesium’s metal-silicate partition coefficient with increasing temperature. D_Cs exhibits the largest change and increased by a factor of three over 500 °C. The effect of oxygen fugacity has not been precisely determined but in general, lowering fO₂ by two log units resulted in a rise in all D-values of approximately an order of magnitude. In general, the sensitivity of partition coefficients to changing parameters increased with atomic number.The highest D-value for Cs observed in this study is 0.345, which was obtained at nbo/t of 2.7 and a metal phase of pure FeS. This metallic composition has far more S than has been suggested for any credible core-forming metal. We therefore conclude that the depletion of Cs in Earth’s mantle is either caused by radically different behavior of Cs at pressures higher than 15 GPa or is not related to core formation. Even so, we have shown that a planet with a sufficient S inventory may incorporate significant amounts of alkali elements into its core.     

Evidence for high-pressure core-mantle differentiation from the metal-silicate partitioning of lithophile and weakly-siderophile elements

Liquid Fe metal-liquid silicate partition coefficients for the lithophile and weakly-siderophile elements Ta, Nb, V, Cr, Si, Mn, Ga, In and Zn have been measured in multianvil experiments performed from 2 to 24 GPa, 2023-2873 K and at oxygen fugacities of −1.3 to −4.2 log units relative to the iron-wüstite buffer. Compositional effects of light elements dissolved in the metal liquid (S, C) have been examined and experiments were performed in both graphite and MgO capsules, specifically to address the effect of C solubility in Fe-metal on siderophile element partitioning. The results were used to examine whether there is categorical evidence that a significant portion of metal-silicate equilibration occurred under very high pressures during core-mantle fractionation on Earth. Although the depletion of V from the mantle due to core formation is significantly greater than that of Nb, our results indicate that both elements have similar siderophile tendencies under reducing conditions at low pressures. With increasing pressure, however, Nb becomes less siderophile than V, implying that average metal-silicate equilibration pressures of at least 10-40 GPa are required to explain the Nb/V ratio of the mantle. Similarly the moderately-siderophile, volatile element ratios Ga/Mn and In/Zn are chondritic in the mantle but both volatility and core-mantle equilibration at low pressure would render these ratios strongly sub-chondritic. Our results indicate that pressures of metal-silicate partitioning exceeding 30-60 GPa would be required to render these element ratios chondritic in the mantle. These observations strongly indicate that metal-silicate equilibration must have occurred at high pressures, and therefore support core-formation models that involve deep magma oceans. Moreover, our results allow us to exclude models that envisage primarily low-pressure (<1 GPa) equilibration in relatively small planetary bodies. We also argue that the core cannot contain significant U as this would require metal-silicate equilibration at oxygen fugacities low enough for significant amounts of Ta to have also been extracted from the mantle. Likewise, as In is more siderophile than Pb but similarly volatile and also quite chalcophile it would have been difficult for Pb to enter the core without reversing the relative depletions of these elements in the mantle unless metal-silicate equilibration occurred at high pressures >20 GPa.     

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

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

Diffusion of Au, Pd, Re, and P in FeNi alloys at High Pressure

Multi-anvil and piston cylinder experiments were performed to determine the effect of both pressure and temperature on the diffusivities of several siderophile elements in a Fe₉₀Ni₁₀ core-analog alloy. Activation energies were calculated to be between 244 and 257 kJ/mol for Re, Pd, and Au at 10 GPa, and 264 kJ/mol for P at 1 GPa. It was found that pressure has a marked negative effect on the diffusivities of Au, Re, and Pd, and activation volumes for these elements were calculated to be between 3 and 6 cm³/mol at 10 GPa. The effect of both temperature and pressure on P diffusion is noticeably less, and the absolute diffusivity of phosphorus is consistently higher than that of the other elements. It is inferred that the reason for this difference is because P is diffusing via an interstitial mechanism as opposed to Re, Pd, and Au which occupy regular lattice sites in the crystal. The effect of pressure and temperature together with depth in the Earth suggests that these elements may continue to exhibit different behavior at more extreme conditions. The significance of these new results lies in the ability to place constraints on many time-dependent processes that pertain to the formation and evolution of planetary cores, as well as the formation and cooling histories of other metal-rich bodies in the solar system, such as meteorites.     

Mineral/melt partitioning of trace elements during hydrous peridotite partial melting

This experimental study examines the mineral/melt partitioning of incompatible trace elements among high-Ca clinopyroxene, garnet, and hydrous silicate melt at upper mantle pressure and temperature conditions. Experiments were performed at pressures of 1.2 and 1.6 GPa and temperatures of 1,185 to 1,370 °C. Experimentally produced silicate melts contain up to 6.3 wt% dissolved H ₂O, and are saturated with an upper mantle peridotite mineral assemblage of olivine+orthopyroxene+clinopyroxene+spinel or garnet. Clinopyroxene/melt and garnet/melt partition coefficients were measured for Li, B, K, Sr, Y, Zr, Nb, and select rare earth elements by secondary ion mass spectrometry. A comparison of our experimental results for trivalent cations (REEs and Y) with the results from calculations carried out using the Wood-Blundy partitioning model indicates that H ₂O dissolved in the silicate melt has a discernible effect on trace element partitioning. Experiments carried out at 1.2 GPa, 1,315 °C and 1.6 GPa, 1,370 °C produced clinopyroxene containing 15.0 and 13.9 wt% CaO, respectively, coexisting with silicate melts containing ~1–2 wt% H ₂O. Partition coefficients measured in these experiments are consistent with the Wood-Blundy model. However, partition coefficients determined in an experiment carried out at 1.2 GPa and 1,185 °C, which produced clinopyroxene containing 19.3 wt% CaO coexisting with a high-H ₂O (6.26±0.10 wt%) silicate melt, are significantly smaller than predicted by the Wood-Blundy model. Accounting for the depolymerized structure of the H ₂O-rich melt eliminates the mismatch between experimental result and model prediction. Therefore, the increased Ca ²⁺ content of clinopyroxene at low-temperature, hydrous conditions does not enhance compatibility to the extent indicated by results from anhydrous experiments, and models used to predict mineral/melt partition coefficients during hydrous peridotite partial melting in the sub-arc mantle must take into account the effects of H ₂O on the structure of silicate melts.     

Systematics of metal-silicate partitioning for many siderophile elements applied to Earth’s core formation

Superliquidus metal-silicate partitioning was investigated for a number of moderately siderophile (Mo, As, Ge, W, P, Ni, Co), slightly siderophile (Zn, Ga, Mn, V, Cr) and refractory lithophile (Nb, Ta) elements. To provide independent constrains on the effects of temperature, oxygen fugacity and silicate melt composition, isobaric (3 GPa) experiments were conducted in piston cylinder apparatus at temperature between 1600 and 2600 °C, relative oxygen fugacities of IW−1.5 to IW−3.5, and for silicate melt compositions ranging from basalt to peridotite. The effect of pressure was investigated through a combination of piston cylinder and multi-anvil isothermal experiments between 0.5 and 18 GPa at 1900 °C. Oxidation states of siderophile elements in the silicate melt as well as effect of carbon saturation on partitioning are also derived from these results. For some elements (e.g. Ga, Ge, W, V, Zn) the observed temperature dependence does not define trends parallel to those modeled using metal-metal oxide free energy data. We correct partitioning data for solute interactions in the metallic liquid and provide a parameterization utilized in extrapolating these results to the P-T-X conditions proposed by various core formation models. A single-stage core formation model reproduces the mantle abundances of several siderophile elements (Ni, Co, Cr, Mn, Mo, W, Zn) for core-mantle equilibration at pressures from 32 to 42 GPa along the solidus of a deep peridotitic magma ocean (∼3000 K for this pressure range) and oxygen fugacities relevant to the FeO content of the present-day mantle. However, these P-T-fO₂ conditions cannot produce the observed concentrations of Ga, Ge, V, Nb, As and P. For more reducing conditions, the P-T solution domain for single stage core formation occurs at subsolidus conditions and still cannot account for the abundances of Ge, Nb and P. Continuous core formation at the base of a magma ocean at P-T conditions constrained by the peridotite liquidus and fixed fO₂ yields concentrations matching observed values for Ni, Co, Cr, Zn, Mn and W but underestimates the core/mantle partitioning observed for other elements, notably V, which can be reconciled if accretion began under reducing conditions with progressive oxidation to fO₂ conditions consistent with the current concentration of FeO in the mantle as proposed by Wade and Wood (2005). However, neither oxygen fugacity path is capable of accounting for the depletions of Ga and Ge in the Earth’s mantle. To better understand core formation, we need further tests integrating the currently poorly-known effects of light elements and more complex conditions of accretion and differentiation such as giant impacts and incomplete equilibration.     

Eutectic temperatures and melting relations in the Fe–O–S system at high pressures and temperatures

We have investigated melting relations in the Fe–O–S ternary system in the pressure range of 15–27 GPa and 1873 K. Subsolidus phase relations are Fe, Fe₃S₂, and FeO up to 17 GPa and Fe, Fe₃S, and FeO above this pressure. The eutectic temperature slightly decreases from ambient pressure to 17 GPa, whereas increases above this pressure. The eutectic temperature in this study is 100 K lower than that in the Fe–S binary system. The oxygen content in the Fe–O–S eutectic liquid drops when the coexisting solid phases changes from FeS to Fe₃S₂. The cotectic lines in the ternary phase diagram lie close to the Fe–FeS binary axis. The isothermal sections indicate that oxygen solubility in the Fe–O–S liquid increases with increasing temperature, and with increasing sulfur content. The solubility of sulfur in the solid Fe has a maximum value at the eutectic temperature, and decreases with increasing temperature. Our results could have important implications for formation and composition of the Martian core.     

Melting relationships in the system Fe-Feo at high pressures: Implications for the composition and formation of the earth's core

A reconnaissance investigation has been carried out on melting relationships in the system Fe-FeO at pressures up to 25 GPa and temperatures up to 2200° C using an MA-8 apparatus. Limited studies were also made of the Co-CoO and Ni-NiO systems. In the system FeFeO, the rapid exsolution of FeO from liquids during quenching causes some difficulties in interpretation of textures and phase relationships. The Co-CoO and Ni-NiO systems are more tractable experimentally and provide useful analogues to the Fe-FeO system. It was found that the broad field of liquid immiscibility present at ambient pressure in the Co-CoO system had disappeared at 18 GPa, 2200° C and that the system displayed complete miscibility between molten Co and CoO, analogous to the behaviour of the Ni-NiO system at ambient pressure. The phase diagram of the system Fe-FeO at 16 GPa and from 1600–2200° C was constructed from interpretations based on the textures of quenched run products. The solubility of FeO in molten iron is considerably enhanced by high pressures. At 16 GPa, the Fe-FeO eutectic contains about 10–15 mol percent FeO and the eutectic temperature in this iron-rich region of the system occurs at 1700±25° C, some 350° C below the melting point of pure iron at the same pressure. The solubility of FeO in molten Fe increases rapidly as temperature increases from 1700 to 2200° C. A relatively small liquid immiscibility field is present above 1900° C but is believed to be eliminated above 2200° C. This inference is supported by thermodynamic calculations on the positions of key phase boundaries. A single run carried out on an Fe₅₀ FeO₅₀ composition at 25 GPa and 2200° C demonstrated extensive and probably complete miscibility between Fe and FeO liquids under these conditions. The melting point of iron is decreased considerably by solution of FeO at high pressures; moreover, the melting point gradient (dP/dT) of the Fe-FeO eutectic is much smaller than that of pure iron and is also smaller than that of mantle pyrolite under the P, T conditions studied. These characteristics make it possible for melting of metal phase and segregation of the core to proceed within the Earth under conditions where most of the mantle remains below solidus temperatures. Under these conditions, the core would inevitably contain a large proportion of dissolved FeO. It is concluded therefore, that oxygen is likely to be the principal light element in the core. The inner core may not be composed of pure iron, as often proposed. Instead, it may consist of a crystalline oxide solid solution (Ni, Fe)₂O.     

The effect of sodium and titanium on crystal-melt partitioning of trace elements

Eclogite of recycled slab origin has often been invoked in the source region of mid-ocean ridge and ocean-island basalts. Melts of this eclogitic material are expected to be enriched in incompatible elements including major elements such as Na and Ti. In order to investigate the controls on trace element chemistry of a melt from such a recycled component we have performed trace element partitioning experiments in the simple systems CMAS (CaO-MgO-Al₂O₃-SiO₂), NCMAS (Na₂O-CMAS) and Ti-CMAS (TiO₂-CMAS) at 3 GPa and 1298-1500°C using analogue eclogitic compositions.We show that sodium has a profound effect on clinopyroxene-melt partition coefficients. NCMAS is characterized by elevated partition coefficients, relative to CMAS for all elements except Li. The increase is more pronounced for more highly charged cations, resulting in negative partitioning anomalies for Li, Sr, Ba and Pb. In contrast to sodium, titanium has very little effect on trace element partitioning for all elements except Nb and Ta, which are retained by the rutile that is saturated in this run.     

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:

     

Trace element partitioning between amphibole and hydrous silicate glasses at 0.6–2.6 GPa

Partitioning behavior between amphibole and silicate glass of thirty-three minor and trace elements (Sc, Ti, V, Cr, Co, Rb, Sr, P, Y, Zr, Nb, Cs, Ba, K, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta, Pb, Th, and U) have been determined experimentally. Products of crystallization of hydrous basalt melts from 0.6 GPa/860 °C up to 2.6 GPa/970 °C were obtained in a multi-anvil apparatus. Major and trace element compositions of amphibole and glass were determined with a combination of electron microprobe and laser ablation inductively coupled plasma mass spectrometry. The main mineral phase is calcic amphibole, and the coexisting glass compositions are tonalite, granodiorite, and granite. The compatibility of rare earth elements increase at 915 °C and then decrease at 970 °C, but the compatibility of most of these elements shows a continued, significant increase with increasing pressure. For high-field strength elements, large ion lithophile elements, actinide compatibility decrease with increasing temperature or pressure, but transition metals show a continued increase in compatibility within the temperature–pressure conditions. From mathematical and graphical fitting, we determined best-fit values for the ideal ionic radius (r₀, 1.01–1.04 Å), the strain-free partition coefficient (D₀, 1.18–1.58), and apparent Young’s modulus (E, 142–370 GPa) for the M4 site in amphibole according to the lattice strain model. The \(D_{0}^{M4}\) for rare earth elements rises at 915 °C and then drops at 970 °C at 0.6 GPa. However, the \(D_{0}^{M4}\) values are positively proportional to the pressure for rare earth elements in the amphibole-glass pairs at 0.6–2.6 GPa and 970 °C. Furthermore, the derived best-fit values for \(r_{0}^{M4}\) and \(E^{M4}\) are almost constant and trend to increase with rising temperature and pressure, respectively. The partition coefficient is distinctly different for different melt compositions. The rare earth elements become more enriched in amphibole if the quenched glass is granodiorite or granite compared to the tonalitic glasses.