Deep Magma Ocean Research Articles

The niobium and tantalum concentration in the mantle constrains the composition of Earth’s primordial magma ocean

The bulk silicate Earth (BSE), and all its sampleable reservoirs, have a subchondritic niobium-to-tantalum ratio (Nb/Ta). Because both elements are refractory, and Nb/Ta is fairly constant across chondrite groups, this can only be explained by a preferential sequestration of Nb relative to Ta in a hidden (unsampled) reservoir. Experiments have shown that Nb becomes more siderophile than Ta under very reducing conditions, leading the way for the accepted hypothesis that Earth's core could have stripped sufficient amounts of Nb during its formation to account for the subchondritic signature of the BSE. Consequently, this suggestion has been used as an argument that Earth accreted and differentiated, for most of its history, under very reducing conditions. Here, we present a series of metal-silicate partitioning experiments of Nb and Ta in a laser-heated diamond anvil cell, at pressure and temperature conditions directly comparable to those of core formation; we find that Nb is more siderophile than Ta under any conditions relevant to a deep magma ocean, confirming that BSE's missing Nb is in the core. However, multistage core formation modeling only allows for moderately reducing or oxidizing accretionary conditions, ruling out the need for very reducing conditions, which lead to an overdepletion of Nb from the mantle (and a low Nb/Ta ratio) that is incompatible with geochemical observations. Earth's primordial magma ocean cannot have contained less than 2% or more than 18% FeO since the onset of core formation.

Metal–silicate partitioning of W and Mo and the role of carbon in controlling their abundances in the bulk silicate earth

The liquid metal–liquid silicate partitioning of molybdenum and tungsten during core formation must be well-constrained in order to understand the evolution of Earth and other planetary bodies, in particular because the Hf–W isotopic system is used to date early planetary evolution. The partition coefficients DMo and DW have been suggested to depend on pressure, temperature, silicate and metal compositions, although previous studies have produced varying and inconsistent models. Additionally, the high cationic charges of W and Mo in silicate melts make their partition coefficients particularly sensitive to oxygen fugacity. We combine 48 new high pressure and temperature experimental results with a comprehensive database of previous experiments to re-examine the systematics of Mo and W partitioning, and produce revised partitioning models from the large combined dataset. W partitioning is particularly sensitive to silicate and metallic melt compositions and becomes more siderophile with increasing temperature. We show that W has a 6+ oxidation state in silicate melts over the full experimental fO2 range of ΔIW −1.5 to −3.5. Mo has a 4+ oxidation state, and its partitioning is less sensitive to silicate melt composition but also depends on metallic melt composition. DMo stays approximately constant with increasing depth in Earth. Both W and Mo become more siderophile with increasing C content of the metal: we therefore performed experiments with varying C concentrations and fit epsilon interaction parameters: εCMo = −7.03 ± 0.30 and εCW = −7.38 ± 0.57.W and Mo along with C are incorporated into a combined N-body accretion and core–mantle differentiation model, which already includes the major rock-forming elements as well as S, and moderately and highly siderophile elements. In this model, oxidation and volatility gradients extend through the protoplanetary disk so that Earth accretes heterogeneously. These gradients, as well as the metal–silicate equilibration pressure, are fitted using a least squares optimisation so that the model Earth-like planet reproduces the composition of the bulk silicate Earth (BSE) in terms of 17 simulated element concentrations (Mg, Fe, Si, Ni, Co, Nb, Ta, V, Cr, S, Pt, Pd, Ru, Ir, W, Mo, and C). The effects of the interaction of W and Mo with Si, S, O, and C in metal are included. Using this model with six separate terrestrial planet accretion simulations, we show that W and Mo require the early accreting Earth to be sulfur-depleted and carbon-enriched so that W and Mo are efficiently partitioned into Earth’s core and do not accumulate in the mantle. When this is the case, the produced Earth-like planets possess mantle compositions matching the BSE for all simulated elements. However, there are two distinct groups of estimates of the bulk mantle’s C abundance in the literature: low (∼100 ppm) and high (∼800 ppm), and all six models are consistent with the higher estimated carbon abundance. The low BSE C abundance would be achievable when the effects of the segregation of dispersed metal droplets produced in deep magma oceans by the disproportionation of Fe2+ to Fe3+ plus metallic Fe is included.

Ag isotopic and chalcophile element evolution of the terrestrial and martian mantles during accretion: New constraints from Bi and Ag metal-silicate partitioning

The Earth's timing of accretion and acquisition of moderately volatile compounds is uncertain. Hafnium-W and Mn-Cr isotopic data can bracket the timing of early planetary differentiation and core formation. The Ag-Pd system has also been utilized but its application has been limited by a lack of high pressure and temperature metal-silicate partitioning for Pd and Ag. Because Ag (and Bi) are volatile chalcophile siderophile elements, understanding their early distribution can constrain the origin of volatile elements in differentiated bodies and planets. Unfortunately, neither Ag or Bi have been studied across the wide range of pressure and temperature conditions that are relevant to accretion and core-mantle differentiation. Here, new high-pressure and temperature multi-anvil metal-silicate equilibrium experiments for Bi and Ag have been carried out at conditions relevant to planetary accretion and metal silicate differentiation that allow a more refined and complete understanding of element partitioning during core formation. The new metal-silicate partitioning data are combined with previously reported data, and utilized to predict the distributions of Bi, Pd, and Ag at conditions of accretion relevant for Earth and Mars. Application of the new partitioning results to Earth shows that D(Bi) and D(Ag) (D = metal/silicate concentration ratio) are lowered due to the effect of pressure and Si alloyed in the metallic liquid, resulting in higher predicted mantle Bi and Ag abundances than in the bulk silicate Earth (BSE), as well as high and variable Pd/Ag. The unradiogenic Ag isotopic composition of the BSE could have been generated by early accretion of volatile-poor (high Pd/Ag) pre-cursors, followed by later accretion of volatile–rich (low Pd/Ag) material, in agreement with earlier studies of Pd-Ag and Mn-Cr (Schönbächler et al., 2010). However, these main accretion phases would have to be followed by segregation of a sulfide liquid (at least 1.5% of magma ocean) at high pressures (>30 GPa), to explain the primitive upper mantle (PUM) Bi, Pd, and Ag, as well as Au, Pt, Cu and Ni concentrations as proposed previously. If the early accreted bulk Earth was volatile depleted with high Pd/Ag ratios, portions of the mantle may contain ancient domains that developed positive 107Ag isotopic anomalies (as also argued by noble gases, Nd, W, and Os isotopes). In comparison, Bi, Pd, and Ag concentrations in the martian mantle could have been set by simple metal-silicate equilibrium. Mars accreted and differentiated relatively rapidly, while also developing a deep magma ocean with a high Pd/Ag ratio that could have evolved positive 107Ag anomalies, in contrast to Earth. Measurements on shergottites may reveal these predicted Ag isotopic anomalies.

The carbon content of Earth and its core

Earth's core is likely the largest reservoir of carbon (C) in the planet, but its C abundance has been poorly constrained because measurements of carbon's preference for core versus mantle materials at the pressures and temperatures of core formation are lacking. Using metal-silicate partitioning experiments in a laser-heated diamond anvil cell, we show that carbon becomes significantly less siderophile as pressures and temperatures increase to those expected in a deep magma ocean during formation of Earth's core. Based on a multistage model of core formation, the core likely contains a maximum of 0.09(4) to 0.20(10) wt% C, making carbon a negligible contributor to the core's composition and density. However, this accounts for ∼80 to 90% of Earth's overall carbon inventory, which totals 370(150) to 740(370) ppm. The bulk Earth's carbon/sulfur ratio is best explained by the delivery of most of Earth's volatiles from carbonaceous chondrite-like precursors.

Atmosphere Origins for Exoplanet Sub-Neptunes

Abstract Planets with 2 R ⊕ < R < 3 R ⊕ and orbital period <100 days are abundant; these sub-Neptune exoplanets are not well understood. For example, Kepler sub-Neptunes are likely to have deep magma oceans in contact with their atmospheres, but little is known about the effect of the magma on the atmosphere. Here we study this effect using a basic model, assuming that volatiles equilibrate with magma at T ∼ 3000 K. For our Fe–Mg–Si–O–H model system, we find that chemical reactions between the magma and the atmosphere and dissolution of volatiles into the magma are both important. Thus, magma matters. For H, most moles go into the magma, so the mass target for both H2 accretion and H2 loss models is weightier than is usually assumed. The known span of magma oxidation states can produce sub-Neptunes that have identical radius but with total volatile masses varying by 20-fold. Thus, planet radius is a proxy for atmospheric composition but not for total volatile content. This redox diversity degeneracy can be broken by measurements of atmosphere mean molecular weight. We emphasize H2 supply by nebula gas, but also consider solid-derived H2O. We find that adding H2O to Fe probably cannot make enough H2 to explain sub-Neptune radii because >103 km thick outgassed atmospheres have high mean molecular weight. The hypothesis of magma–atmosphere equilibration links observables such as atmosphere H2O/H2 ratio to magma FeO content and planet formation processes. Our model’s accuracy is limited by the lack of experiments (lab and/or numerical) that are specific to sub-Neptunes; we advocate for such experiments.

Deep magma ocean formation set the oxidation state of Earth's mantle.

The composition of Earth's atmosphere depends on the redox state of the mantle, which became more oxidizing at some stage after Earth's core started to form. Through high-pressure experiments, we found that Fe2+ in a deep magma ocean would disproportionate to Fe3+ plus metallic iron at high pressures. The separation of this metallic iron to the core raised the oxidation state of the upper mantle, changing the chemistry of degassing volatiles that formed the atmosphere to more oxidized species. Additionally, the resulting gradient in redox state of the magma ocean allowed dissolved CO2 from the atmosphere to precipitate as diamond at depth. This explains Earth's carbon-rich interior and suggests that redox evolution during accretion was an important variable in determining the composition of the terrestrial atmosphere.

In situ trace element and Sm-Nd isotope analysis of accessory minerals in an Eoarchean tonalitic gneiss from Greenland: Implications for Hf and Nd isotope decoupling in Earth's ancient rocks

Some of the Earth's oldest preserved continental crustal rocks, from southern West Greenland, contain contradictory radiogenic isotope signatures, whereby Hf isotope ratios are chondritic but Nd isotope ratios are distinctly super-chondritic. Models to explain this discrepancy are speculative and variously invoke deep magma ocean crystallisation, HfNd decoupling in subduction zones, or metamorphic disturbance during younger thermal events. Determining the cause of this discrepancy is essential for understanding Eoarchean crust-mantle differentiation. We employ, for the first time, micro-analysis of REE-rich accessory minerals to shed light on the Nd isotope evolution of a key tonalitic gneiss sample from southern West Greenland that displays the apparent HfNd isotope decoupling. The results show that the Sm-Nd isotope system was homogenized during a metamorphic event at ca. 2690 Ma. We suggest that metamorphic reactions involving consumption and re-crystallisation of REE-bearing phases were accompanied by LREE element mobility and the loss of unradiogenic Nd, shifting the bulk rock composition to a more radiogenic Nd isotope value. Our study provides the first direct evidence that the anomalous Nd isotope signatures in some Eoarchean gneisses are artefacts of the disturbance of the Sm-Nd isotope system, and not due to extensive differentiation of the bulk silicate Earth by magma ocean crystallisation or continental crust formation.

Dynamics and stability of an iron drop falling in a magma ocean

The latest stages of planetary accretion involved large impacts between differentiated bodies, hence large scale melting events. Consequently, the iron brought by the impactors sank within a deep magma ocean, before reaching the proto-core. Yet the fluid dynamics of this process remains poorly known. Here, we report numerical simulations of the sinking dynamics of an initially spherical liquid iron drop within a molten silicate phase, up to its possible fragmentation. We consider a 2D cylindrical axisymmetric geometry. We vary the viscosity of the molten silicates in the range of 0.05 Pa.s–100 Pa.s and the initial radius of the iron drop in the range of 1 mm–350 mm. Hence, we investigate Reynolds number in the range of [0.027–85600] and Weber number in the range of [0.073–7480]. Our numerical model constrains the morphology, dynamics and stability of the iron drop as a function of the dimensionless Weber and Reynolds numbers as well as of the viscosity ratio between the molten silicates and the liquid iron drop. In particular, we show that the maximal stable drop radius and the critical Weber number are monotonically increasing functions of the magma ocean viscosity. The momentum boundary layer thickness depends mainly on the drop radius and slightly on the magma ocean viscosity. Increasing the viscosity of the silicate phase prevents oscillations of the iron phase and limits the exchange surface. Oppositely, increasing the initial radius of the iron drop enhances its deformation and increases its relative exchange surface. Above the critical Weber number, we confirm that the fragmentation of the liquid iron occurs within a falling distance equal to 3.5–8 times the drop initial radius in the explored range of moderate Weber number, and we describe a variety of fragmentation regimes. Consequences for Earth’s formation models are briefly assessed.

An equation of state for high pressure-temperature liquids (RTpress) with application to MgSiO3 melt

The thermophysical properties of molten silicates at extreme conditions are crucial for understanding the early evolution of Earth and other massive rocky planets, which is marked by giant impacts capable of producing deep magma oceans. Cooling and crystallization of molten mantles are sensitive to the densities and adiabatic profiles of high-pressure molten silicates, demanding accurate Equation of State (EOS) models to predict the early evolution of planetary interiors. Unfortunately, EOS modeling for liquids at high P-T conditions is difficult due to constantly evolving liquid structure. The Rosenfeld-Tarazona (RT) model provides a physically sensible and accurate description of liquids but is limited to constant volume heating paths (Rosenfeld and Tarazona, 1998). We develop a high P-T EOS for liquids, called RTpress, which uses a generalized Rosenfeld-Tarazona model as a thermal perturbation to isothermal and adiabatic reference compression curves. This approach provides a thermodynamically consistent EOS which remains accurate over a large P-T range and depends on a limited number of physically meaningful parameters that can be determined empirically from either simulated or experimental datasets. As a first application, we model MgSiO3 melt representing a simplified rocky mantle chemistry. The model parameters are fitted to the MD simulations of both Spera et al. (2011) and de Koker and Stixrude (2009), recovering pressures, volumes, and internal energies to within 0.6 GPa, 0.1 Å3, and 6 meV per atom on average (for the higher resolution data set), as well as accurately predicting liquid densities and temperatures from shock-wave experiments on MgSiO3 glass. The fitted EOS is used to determine adiabatic thermal profiles, revealing the approximate thermal structure of a fully molten magma ocean like that of the early Earth. These adiabats, which are in strong agreement for both fitted models, are shown to be sufficiently steep to produce either a center-outwards or bottom-up style of crystallization, depending on the curvature of the mantle melting curve (liquidus), with a high-curvature model yielding crystallization at depths of roughly 80 GPa (Stixrude et al., 2009) whereas a nearly-flat experimentally determined liquidus implies bottom-up crystallization (Andrault et al., 2011).

Polyamorphic Transformations in Fe‐Ni‐C Liquids: Implications for Chemical Evolution of Terrestrial Planets

AbstractDuring the formation of the Earth's core, the segregation of metallic liquids from silicate mantle should have left behind evident geochemical imprints on both the mantle and the core. Some distinctive geochemical signatures of the mantle‐derived rocks likely own their origin to the metal‐silicate differentiation of the primitive Earth, setting our planet apart from undifferentiated meteorites as well as terrestrial planets or moons isotopically and compositionally. Understanding the chemical evolution of terrestrial planetary bodies requires knowledge on properties of both liquid iron alloys and silicates equilibrating under physicochemical conditions pertinent to the deep magma ocean. Here we report experimental and computational results on the pressure‐induced structural evolution of iron‐nickel liquids alloyed with carbon. Our X‐ray diffraction experiments up to 7.3 gigapascals (GPa) demonstrate that Fe‐Ni (Fe90Ni10) liquids alloyed with 3 and 5 wt % carbon undergo a polyamorphic liquid structure transition at approximately 5 GPa. Corroborating the experimental observations, our first‐principles molecular dynamic calculations reveal that the structural transitions result from the marked prevalence of three‐atom face‐sharing polyhedral connections in the liquids at >5 GPa. The structure and polyamorphic transitions of liquid iron‐nickel‐carbon alloys govern their physical and chemical properties and may thus cast fresh light on the chemical evolution of terrestrial planets and moons.

Melting phase relations in the MgSiO3\u2013CaSiO3 system at 24 GPa

The Earth’s lower mantle is composed of bridgmanite, ferropericlase, and CaSiO3-rich perovskite. The melting phase relations between each component are key to understanding the melting of the Earth’s lower mantle and the crystallization of the deep magma ocean. In this study, melting phase relations in the MgSiO3–CaSiO3 system were investigated at 24 GPa using a multi-anvil apparatus. The eutectic composition is (Mg,Ca)SiO3 with 81–86 mol% MgSiO3. The solidus temperature is 2600–2620 K. The solubility of CaSiO3 component into bridgmanite increases with temperature, reaching a maximum of 3–6 mol% at the solidus, and then decreases with temperature. The same trend was observed for the solubility of MgSiO3 component into CaSiO3-rich perovskite, with a maximum of 14–16 mol% at the solidus. The asymmetric regular solutions between bridgmanite and CaSiO3-rich perovskite and between MgSiO3 and CaSiO3 liquid components well reproduce the melting phase relations constrained experimentally.

Experimental constraints on the solidification of a nominally dry lunar magma ocean

The lunar magma ocean (LMO) concept has been used extensively for lunar evolution models for decades, but to date the full cooling and crystallization path of the LMO has not been studied experimentally. Here we present results of a high-pressure, high-temperature experimental study of the mineralogical and geochemical evolution accompanying the full solidification of a nominally dry LMO. Experiments used a bulk composition based on geophysical data, and assumed an initial LMO depth of 700 km. The effect of pressure within a deep magma ocean on solidification at different levels in the ocean was explicitly taken into account, by performing experiments at multiple pressures and constant temperature during each solidification step. Results show formation of a deep harzburgite (olivine + low-Ca pyroxene) layer in the first ∼50% of equilibrium crystallization. The crystallising mineral assemblage does not change until plagioclase and clinopyroxene appear at 68 PCS (per cent solid by volume), while low-Ca pyroxene stops forming. Olivine disappears at 83 PCS, and ilmenite and β-quartz start crystallizing at 91 and 96 PCS, respectively. At 99 PCS, we observe an extremely iron-rich (26.5 wt.% FeO) residual LMO liquid. Our results differ substantially from the oft-cited LMO solidification study of Snyder et al. (1992), which was based on a limited number of experiments at a single pressure. Differences include the mineralogy of the deepest sections of the solidified LMO (harzburgitic instead of dunitic), the formation of SiO2 in the lunar interior, and the development of extreme iron enrichment in the last remaining dregs of the LMO. Our findings shed new light on several aspects of lunar petrology, including the formation of felsic and iron-rich magmas in the Moon. Finally, based on our experiments the lunar crust, consisting of the light minerals plagioclase and quartz, would reach a thickness of ∼67.5 km. This is far greater than crustal thickness estimates from recent GRAIL mission gravitational data (34–43 km, Wieczorek et al., 2013). Although the initial depth of the LMO has an effect on the thickness of crust produced, this effect is not large enough to explain this discrepancy. Inefficient plagioclase segregation, trapping of magma in cumulate reservoirs, and Al sequestration in spinel cannot explain the discrepancy either. As plagioclase crystallization can be suppressed by the presence of H2O, this implies that the lunar magma ocean was water-bearing.

A sulfur-poor terrestrial core inferred from metal–silicate partitioning experiments

As a siderophile and a volatile element, sulfur's partitioning behavior allows constraints to be placed on processes in the primitive Earth. Sulfur's core–mantle distribution during Earth's accretion has consequences for core content and implications for volatile accretion. In this study, metal–silicate partitioning experiments of sulfur were conducted in a diamond anvil cell at pressures from 46 to 91 GPa and temperatures between 3100 and 4100 K, conditions that are directly relevant to core segregation in a deep magma ocean. The sulfur partition coefficients measured from these experiments are an order of magnitude less than those obtained from extrapolation of previous results to core formation conditions (e.g., Rose-Weston et al., 2009; Boujibar et al., 2014). These measurements challenge the idea that sulfur becomes a highly siderophile element at high pressures and temperatures. A relationship was derived that describes sulfur's partitioning behavior at the pressure–temperature range of core formation. This relationship combined with an accretion model was used to explore the effects of varying impactor sizes and volatile compositions on the sulfur contents of the Earth's core and mantle. The results show that homogeneous delivery of sulfur throughout accretion would overenrich the mantle in sulfur relative to the present day observations of 200 ± 80 ppm (Lorand et al., 2013) unless the bulk Earth sulfur content is lower than its cosmochemical estimate of ∼6400 ppm (e.g., McDonough, 2003). On the other hand, the mantle's sulfur content is matched if sulfur is delivered with large bodies (3 to 10% Earth mass) during the last 20% of Earth's accretion, combined with a chondritic late veneer of 0.5% Earth mass. These results are conditional on the lowered equilibration efficiency of large impactor cores in a terrestrial magma ocean. In each accretion scenario, the core sulfur content remains below ∼2 wt.% in close agreement with cosmochemical estimates and is a further indication that sulfur is not a dominant light element in the core.

Crystallization of silicon dioxide and compositional evolution of the Earth's core.

The Earth's core is about ten per cent less dense than pure iron (Fe), suggesting that it contains light elements as well as iron. Modelling of core formation at high pressure (around 40-60 gigapascals) and high temperature (about 3,500 kelvin) in a deep magma ocean predicts that both silicon (Si) and oxygen (O) are among the impurities in the liquid outer core. However, only the binary systems Fe-Si and Fe-O have been studied in detail at high pressures, and little is known about the compositional evolution of the Fe-Si-O ternary alloy under core conditions. Here we performed melting experiments on liquid Fe-Si-O alloy at core pressures in a laser-heated diamond-anvil cell. Our results demonstrate that the liquidus field of silicon dioxide (SiO2) is unexpectedly wide at the iron-rich portion of the Fe-Si-O ternary, such that an initial Fe-Si-O core crystallizes SiO2 as it cools. If crystallization proceeds on top of the core, the buoyancy released should have been more than sufficient to power core convection and a dynamo, in spite of high thermal conductivity, from as early on as the Hadean eon. SiO2 saturation also sets limits on silicon and oxygen concentrations in the present-day outer core.

Deformation-aided segregation of Fe-S liquid from olivine under deep Earth conditions: Implications for core formation in the early solar system

The planets and larger rocky bodies of the inner solar system are differentiated, and consist of metallic, iron-rich cores surrounded by thick shells of silicate. Core formation in these bodies, i.e. the segregation of metal from silicate, was a key process in the early solar system, and one which left a lasting geochemical signature. It is commonly assumed that extensive silicate melting and formation of deep magma oceans was required to initiate core formation, due to the inability of iron-rich melts to segregate from a solid silicate matrix. Here we assess the role of deformation in aiding segregation of core-forming melts from solid silicate under conditions of planetary deep interiors. Low-strain rate, high-pressure/temperature deformation experiments and high-resolution 2-D and 3-D textural analysis demonstrate that deformation fundamentally alters iron-rich melt geometry, promoting wetting of silicate grain boundaries and formation of extensive micron to sub-micron width Fe-rich melt bands. Deformation-aided Fe-S melt networks noted here contrast those observed in higher finite strain experiments conducted at lower pressure, and may reveal either an alternative mechanism for melt segregation at higher pressures, or an early stage process of melt segregation. Results suggest, however, that core-mantle chemical equilibration cannot be assumed in models of planetary formation, and that instead, the chemistry of rocky planets may record a complex, multi-stage process of core formation.

Silicate melts during Earth's core formation

Accretion from primordial material and its subsequent differentiation into a planet with core and mantle are fundamental problems in terrestrial and solar system. Many of the questions about the processes, although well developed as model scenarios over the last few decades, are still open and much debated. In the early Earth, during its formation and differentiation into rocky mantle and iron-rich core, it is likely that silicate melts played an important part in shaping the Earth's main reservoirs as we know them today.Here, we review several recent results in a deep magma ocean scenario that give tight constraints on the early evolution of our planet. These results include the behaviour of some siderophile elements (Ni and Fe), lithophile elements (Nb and Ta) and one volatile element (Helium) during Earth's core formation. We will also discuss the melting and crystallization of an early magma ocean, and the implications on the general feature of core-mantle separation and the depth of the magma ocean. The incorporation of Fe2+ and Fe3+ in bridgmanite during magma ocean crystallization is also discussed. All the examples presented here highlight the importance of the prevailing conditions during the earliest time of Earth's history in determining the composition and dynamic history of our planet.

Evidence for an early wet Moon from experimental crystallization of the lunar magma ocean

The Moon is thought to have been covered initially by a deep magma ocean, its gradual solidification leading to the formation of the plagioclase-rich highland crust. We performed a high-pressure, high-temperature experimental study of lunar mineralogical and geochemical evolution during magma ocean solidification that yields constraints on the presence of water in the earliest lunar interior. In the experiments, a deep layer containing both olivine and pyroxene is formed in the first ∼50% of crystallization, β-quartz forms towards the end of crystallization, and the last per cent of magma remaining is extremely iron rich. In dry experiments, plagioclase appears after 68 vol.% solidification and yields a floatation crust with a thickness of ∼68 km, far above the observed average of 34–43 km based on lunar gravity. The volume of plagioclase formed during crystallization is significantly less in water-bearing experiments. Using the relationship between magma water content and the resulting crustal thickness in the experiments, and considering uncertainties in initial lunar magma ocean depth, we estimate that the Moon may have contained at least 270 to 1,650 ppm water at the time of magma ocean crystallization, suggesting the Earth–Moon system was water-rich from the start. The Moon is thought to have initially had a magma ocean that gradually solidified. Crystallization experiments find that the resulting crustal thickness depends on water content and is consistent with significant water in the early Moon.

The effects of shear deformation on planetesimal core segregation: Results from in-situ X-ray micro-tomography

It is well accepted that the Earth formed by the accretion and collision of small (10–100 km), rocky bodies called planetesimals. W-Hf isotopic evidence from meteorites suggest that the cores of many planetesimals formed within a relatively short time frame of ~3 My. While a very hot, deep magma ocean is generally thought to have been the driving mechanism for core formation in large planetary bodies, it inadequately explains differentiation and core formation in small planetesimals due to temperatures potentially being insufficient for wide-scale silicate melting to occur. In order for these planetesimals to differentiate within such a relatively short time without a magma ocean, a critical melt volume of the metallic (core-forming) phase and sufficient melt connectivity and grain size must have existed to attain the required permeability and lead to efficient core formation. Shear deformation may increase the connectedness of melt and the permeability, and thus could have been a major contributing factor in the formation of planetesimal cores. This deformation may have been caused by large impacts and collisions experienced by the planetesimals in the early solar system. The purpose of this work is to test the hypothesis that shear deformation enhances the connectivity and permeability of Fe-S melt within a solid silicate (olivine) matrix, such that rapid core formation is plausible. A rotational Drickamer apparatus (RDA) was used to heat and torsionally deform a sample of solid olivine + FeS liquid through six steps of large-strain shear deformation. After each deformation step, X-ray microtomographs were collected in the RDA to obtain in situ three-dimensional images of the sample. The resulting digital volumes were processed and permeability simulations utilizing the lattice Boltzmann method were performed to determine the effect of shear deformation on connectivity and permeability within the sample. The resulting permeabilities of the sample at various steps of deformation are the same within uncertainty and do not exhibit a change with increasing deformation. Additionally, the migration velocity calculated from the permeability of the sample is not high enough for segregation to take place within the time frame of ~3 My. In addition to further constraining the mechanism of core formation in planetesimals, the image processing techniques developed in this study will be of great benefit to future studies utilizing similar methods.

Differentiated planetesimal impacts into a terrestrial magma ocean: Fate of the iron core

The abundance of moderately siderophile elements (“iron-loving”; e.g. Co, Ni) in the Earth's mantle is 10 to 100 times larger than predicted by chemical equilibrium between silicate melt and iron at low pressure, but it does match expectation for equilibrium at high pressure and temperature. Recent studies of differentiated planetesimal impacts assume that planetesimal cores survive the impact intact as concentrated masses that passively settle from a zero initial velocity and undergo turbulent entrainment in a global magma ocean; under these conditions, cores greater than 10 km in diameter do not fully mix without a sufficiently deep magma ocean. We have performed hydrocode simulations that revise this assumption and yield a clearer picture of the impact process for differentiated planetesimals possessing iron cores with radius = 100 km that impact into magma oceans. The impact process strips away the silicate mantle of the planetesimal and then stretches the iron core, dispersing the liquid iron into a much larger volume of the underlying liquid silicate mantle. Lagrangian tracer particles track the initially intact iron core as the impact stretches and disperses the core. The final displacement distance of initially closest tracer pairs gives a metric of core stretching. The statistics of stretching imply mixing that separates the iron core into sheets, ligaments, and smaller fragments, on a scale of 10 km or less. The impact dispersed core fragments undergo further mixing through turbulent entrainment as the molten iron fragments rain through the magma ocean and settle deeper into the planet. Our results thus support the idea that iron in the cores of even large differentiated planetesimals can chemically equilibrate deep in a terrestrial magma ocean.

Constraints on core formation in Vesta from metal–silicate partitioning of siderophile elements

It is now widely accepted that the asteroid 4-Vesta has an Fe-rich metallic core, but the composition of the core and the conditions prevailing during core–mantle differentiation are poorly constrained. In light of new constraints on Vesta’s geophysical and geochemical properties obtained by the DAWN mission, we have re-examined the conditions at which core–mantle differentiation in Vesta may have occurred by linking the estimated mantle depletions of siderophile elements P, Co, Ni, Cu, Ga, Ge, Mo and W in the vestan mantle to newly derived predictive equations for core–mantle partitioning of these elements. We extend the number of elements previously considered in geochemical modeling of vestan core formation and use published metal–silicate partitioning data obtained at low pressures to characterize the dependence of metal/silicate partition coefficients (D) on pressure, temperature, oxygen fugacity and composition of the silicate and metallic melt. In our modeling we implement newly derived mantle depletions of P, Co, Ni and Ga through analysis of published HED meteorite analyses and assess two contrasting bulk compositional models for Vesta.Modeling results using Monte Carlo simulations constrain vestan core formation to have occurred at mildly reducing conditions of approximately 2log units below the iron–wüstite (IW) buffer (ΔIW=−2.05±0.20) if the two most likely bulk compositions (binary mixtures of H+CM or H+CV chondritic meteorites) are considered, assuming a temperature range between 1725 and 1850K and a sulfur-free pure Fe core. If the core is assumed to be sulfur-rich (15wt.% S) as predicted by the latter bulk compositional models, observed depletions for all eight siderophile elements can be simultaneously satisfied at ΔIW=−2.35±0.10 and 1725–1850K for the H+CV bulk composition and ΔIW=−2.30±0.15 and 1725–1850K for the H+CM bulk composition. More reducing conditions are not consistent with the observed siderophile element depletions in the vestan mantle, independent of the sulfur content of the vestan core and of the bulk compositional model chosen.Our analysis shows that a previously proposed shallow magma ocean on Vesta during core formation is not consistent with the observed mantle depletion of Ga for the two considered bulk compositions irrespective of core composition. Instead, our results are consistent with the existence of a deep magma ocean during core formation on Vesta, requiring >50–100% mantle melting.