skip to main content

Title: High Pressure Melting Curve of Fe Determined by Inter‐Metallic Fast Diffusion Technique

The heat extracted from the core by the overlying mantle across the core‐mantle boundary controls the thermal evolution of the core. This in turn leads to the solidification of the inner core in association with the exsolution of light alloying elements into the liquid outer core. Although the temperature (T) at the inner core boundary (ICB) would be adjusted to account for the effects of the light elements, the melting T of Fe places an upper bound at the ICB and it is a vital point in the thermal profile of the core. Here, we determine the melting T of Fe in the multi‐anvil press by characterizing the interface of Fe‐W interaction. Our data place a tighter constraint on the melting curve of Fe between 8 and 21 GPa, that is directly applicable to small planetary bodies and serves as an anchor for melting curve of Fe at higher pressure.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Experiments in laser‐heated diamond anvil cells (LH DACs) are conducted to assess phase diagrams of planetary materials at high pressure‐temperature (P‐T) conditions; thus, reliable determination of temperature in LH DAC experiments is essential. Radiometric temperature determination in LH DACs relies on the assumption of sample's wavelength‐independent optical properties (graybody assumption), which is not justified for major lower mantle materials. The result is that experimental phase diagrams contain systematic unconstrained errors. Here we estimate the systematic error in radiometric temperature of nongray polycrystalline bridgmanite (Bgm; Mg0.96Fe2+0.036Fe3+0.014Si0.99O3) in a LH DAC by modeling emission and absorption of thermal radiation in a sample with experimentally‐constrained optical properties. A comparison to experimental data validates the models and reveals that thermal spectra measured in LH DAC experiments record the interaction of radiation with the hot nongray sample. The graybody assumption in the experiments on translucent Bgm (light extinction coefficient,k < 250 cm‐1at 500–900 nm) yields temperatures ∼5% higher than the maximum temperature in the sample heated to ∼1900 K. In contrast, the graybody temperature of dark Bgm (k > ∼1500 cm−1), such as that produced upon melt quenching in LH DACs, underestimates the maximum temperature by ∼10%. Our experimental results pose quantitative constraints on the effect of nongray optical properties on the uncertainty of radiometric temperature determination in Bgm in the LH DACs. Evaluating nongray temperature in the future would enable a revision of the Bgm to post‐perovskite phase transition and the high‐pressure melting curve of Bgm.

    more » « less
  2. Light elements in Earth’s core play a key role in driving convection and influencing geodynamics, both of which are crucial to the geodynamo. However, the thermal transport properties of iron alloys at high-pressure and -temperature conditions remain uncertain. Here we investigate the transport properties of solid hexagonal close-packed and liquid Fe-Si alloys with 4.3 and 9.0 wt % Si at high pressure and temperature using laser-heated diamond anvil cell experiments and first-principles molecular dynamics and dynamical mean field theory calculations. In contrast to the case of Fe, Si impurity scattering gradually dominates the total scattering in Fe-Si alloys with increasing Si concentration, leading to temperature independence of the resistivity and less electron–electron contribution to the conductivity in Fe-9Si. Our results show a thermal conductivity of ∼100 to 110 W⋅m −1 ⋅K −1 for liquid Fe-9Si near the topmost outer core. If Earth’s core consists of a large amount of silicon (e.g., > 4.3 wt %) with such a high thermal conductivity, a subadiabatic heat flow across the core–mantle boundary is likely, leaving a 400- to 500-km-deep thermally stratified layer below the core–mantle boundary, and challenges proposed thermal convection in Fe-Si liquid outer core. 
    more » « less
  3. Abstract

    An internally generated magnetic field once existed on the Moon. This field reached high intensities (∼10–100μT, perhaps intermittently) from ∼4.3 to 3.6 Gyr ago and then weakened to ≲5μT before dissipating by ∼1.9–0.8 Gyr ago. While the Moon’s metallic core could have generated a magnetic field via a dynamo powered by vigorous convection, models of a core dynamo often fail to explain the observed characteristics of the lunar magnetic field. In particular, the core alone may not contain sufficient thermal, chemical, or radiogenic energy to sustain the high-intensity fields for >100 Myr. A recent study by Scheinberg et al. suggested that a dynamo hosted in electrically conductive, molten silicates in a basal magma ocean (BMO) may have produced a strong early field. However, that study did not fully explore the BMO’s coupled evolution with the core. Here we show that a coupled BMO–core dynamo driven primarily by inner core growth can explain the timing and staged decline of the lunar magnetic field. We compute the thermochemical evolution of the lunar core with a 1D parameterized model tied to extant simulations of mantle evolution and BMO solidification. Our models are most sensitive to four parameters: the abundances of sulfur and potassium in the core, the core’s thermal conductivity, and the present-day heat flow across the core–mantle boundary. Our models best match the Moon’s magnetic history if the bulk core contains ∼6.5–8.5 wt% sulfur, in agreement with seismic structure models.

    more » « less
  4. Abstract

    A hemispherical pattern in inner core attenuation (Q) has been observed from the apparentQof inner core transmitted waves (PKPDF). How the elastic scattering (QSC) originating from kilometer‐scale inner core heterogeneities relates to large‐scale apparentQvariations remains elusive. Such inner core scattered (ICS) energy is characterized by emergent, long‐lasting, high‐frequency (1–4 Hz) coda immediately following pre‐critical reflections (PKiKP) from the inner core boundary (ICB). In this study, we develop a framework to systematically investigate ICS and examine its hemispherical pattern using data from two arrays that sample opposing sides of the Pacific quasi‐hemispherical boundary. We use all the viable data from earthquakes (Mw ≥ 5.8) within the past 2–3 decades recorded at distances of 50°–75° by YKA and ILAR—small‐aperture (∼10–20 km) seismic arrays situated in northwestern North America. Two independent waveform stripping methods are applied to extract ICS from the background energy, and various stacks are performed to identify ICS and its spatial pattern. We found that individual beams lacking clear evidence ofPKiKPcoda waves reveal ICS characteristics when stacked together, implying that ICS is ubiquitous. We used a modified phonon‐based simulation to reinforce the idea that ICS is primarily created by volumetric heterogeneity within the inner core as opposed to ICB topography. With simplified two‐layer ICS models, our results suggest that ICS within the eastern quasi‐hemisphere is slightly stronger than in the western quasi‐hemisphere, although intra‐hemispherical variations are as significant, and our sampling is limited to patches of the northern hemisphere.

    more » « less
  5. Abstract

    At nearly 2,900‐km depth, the core‐mantle boundary (CMB) represents the largest density increase within the Earth going from a rocky mantle into an iron‐alloy core. This compositional change sets up steep temperature gradients, which in turn influences mantle flow, structure, and seismic velocities. Here we resolve the thermodynamic parameters of (Mg,Fe)O and compute the melting phase relations of the MgO‐FeO binary system at CMB conditions. Based on this phase diagram, we revisit iron infiltration into solid ferropericlase along the CMB by morphological instability and find that the length scale of infiltration is comparable with the high electrical conductivity layer inferred from core nutations. We also compute the (Mg,Fe)O‐SiO2pseudo‐binary system and find that the solidus melting temperatures near the CMB decrease with FeO and SiO2content, becoming potentially important for ultralow velocity zones. Therefore, an ultralow velocity zone composed of solid‐state bridgmanite and ferropericlase may be relatively enriched in MgO and depleted in SiO2and FeO along a hot CMB.

    more » « less