Earth's Core Research Articles

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. The potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: Tm=1969.23+19.15P−0.012P2. In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

The fastest growing modes in rotating magnetoconvection with anisotropic diffusivities and their links to the Earth's core conditions

The fastest growing modes in rotating magnetoconvection with anisotropic diffusivities and their links to the Earth's core conditions

Condition for metal fragmentation during Earth-forming collisions

The long-term evolution of the deep Earth depends on its initial temperature and composition. These were set by the large planetary collisions that formed the Earth. After each collision, the metallic core of the impactor fell into a molten silicate magma ocean. Previous investigations showed that, as it sank, the impactor core fragmented into drops. The overall fragmentation of the core controlled the efficiency of chemical transfers between the impactor metal and the magma ocean, and, as a consequence, the composition of the Earth's core and mantle. However, because previous studies lack an impact stage, it is unclear whether the projectile core fragmented during the impact at Earth's surface, or deeper in the magma ocean.To answer this question, we conduct laboratory experiments modeling the collision of single-phase and two-phase impactors. In a first series of experiments, we investigate the impact of a single-phase centimetric liquid volume, representing the impactor core, onto a lighter immiscible liquid, representing the magma ocean. Our experiments approach the dynamical regime of planetary collisions for which inertia is large compared to surface tension. Varying the velocity and size of the impactor, we determine the conditions under which the impactor fragments into drops. We find that fragmentation occurs when the Froude number, which measures the relative importance of inertia to gravity, is larger than 40, regardless of surface tension. This fragmentation results from the growth of a turbulent Rayleigh-Taylor instability at the interface between the impacting liquid and the target pool. In contrast, when Fr<10, the impactor remains coherent. In a second series of experiments, we use two-phase impactors to show that these results hold for impactors that are differentiated into a core and a mantle.Applied to planet formation, our results suggest that the core of impactors less than 330 km in radius impacting at the escape velocity onto an Earth-sized planet fully fragments into droplets during the impact process, whereas the core of a giant Mars-sized impactor remains coherent. We derive a model for the depth at which the impactor core fragments in the magma ocean as a function of the impactor size and velocity. This model predicts that impactors with a radius less than 800 km fully fragment before reaching the bottom of the magma ocean. For velocities higher than twice the escape speed, some degree of fragmentation is unavoidable for any impactor size.

Reassessment of Birch's Law on hcp‐Fe From Ultrasonic Sound Velocity Measurement and Implications on the Velocity Profiles of Earth's Inner Core

AbstractWe performed in situ X‐ray diffraction and ultrasonic sound velocity measurements on hcp‐Fe up to 15 GPa, 873 K in a multi‐anvil apparatus. The elastic moduli and their pressure and temperature derivatives were determined by fitting the velocity and density data to the third‐order finite strain equations, yielding KS0 = 169.0(57) GPa, KS0′ = 5.4(6), (∂KS/∂T)P = −0.031(3) GPa/K, G0 = 104.5(27) GPa, G0′ = 1.7(2), (∂G/∂T)P = −0.060(2) GPa/K. Within the experimental P‐T range, we find significant temperature effect on the density‐velocity (ρ − VP or S) relations and caution the use of the temperature‐independent Birch's law for extrapolation to core conditions. Furthermore, temperature‐induced velocity decrease is more significant in VS than in VP, offering a possible explanation for the high Poisson's ratio in the core. Extrapolations based on our results together with previous experimental data suggest that VP of hcp‐Fe aligns with PREM at Earth's core conditions, while Vs and density are approximately 10% and 2.7% higher than PREM, respectively.

A Scaling Relation for Core Heating by Giant Impacts and Implications for Dynamo Onset

AbstractAccretional heating of Earth's interior during formation is pivotal to its subsequent thermal and chemical evolution. In particular, impact heating of Earth's core is expected, but its amplitude and radial distribution within the core is unknown and could influence the onset of the geodynamo. The uncertainty is due, in part, to the lack of constraints on the temperature of the interior following formation due to the difficulty of preserving a record of such a high energy environment, and the assertion that super‐heating during formation would be rapidly lost through magma ocean cooling. Here we systematically investigate core heating due to giant impacts using a Smoothed Particle Hydrodynamics (SPH) code with simulations spanning a range of impact angles, velocities, and masses. From these simulations we derive a scaling relation for core heating that depends on the impact parameters and predicts the radial core temperature profile following the impact. Our findings show that a significant amount of heat is deposited into the core, with a canonical impact scenario resulting in an average core temperature increase of about 3000 K, approximately 500 K higher than that of the overlying mantle. In this case the heat distribution within the core produces a strong thermal stratification. We use a parameterized cooling model to estimate that the core could have cooled to an adiabatic state ∼290 Myr after a canonical impact, which is consistent with the observed time span between the age of the Moon and evidence for an active geodynamo.

Iron alloys of volatile elements in the deep Earth's interior.

Investigations into the compositional model of the Earth, particularly the atypical concentrations of volatile elements within the silicate portion of the early Earth, have attracted significant interest due to their pivotal role in elucidating the planet's evolution and dynamics. To understand the behavior of such volatile elements, an established 'volatility trend' has been used to explain the observed depletion of certain volatile elements. However, elements such as Se and Br remain notably over-depleted in the silicate Earth. Here we show the results from first-principles simulations that explore the potential for these elements to integrate into hcp-Fe through the formation of substitutional alloys, long presumed to be predominant constituents of the Earth's core. Based on our findings, the thermodynamic stability of these alloys suggests that these volatile elements might indeed be partially sequestered within the Earth's core. We suggest potential reservoirs for volatile elements within the deep Earth, augmenting our understanding of the deep Earth's composition.

Imaging the top of the Earth's inner core: a present-day flow model.

Despite considerable progress in seismology, mineral physics, geodynamics, paleomagnetism, and mathematical geophysics, Earth's inner core structure and evolution remain enigmatic. One of the most significant issues is its thermal history and the current thermal state. Several hypotheses involving a thermally-convecting inner core have been proposed: a simple, high-viscosity, translational mode, or a classical, lower-viscosity, plume-style convection. Here, we use state-of-the-art seismic imaging to probe the outermost shell of the inner core for its isotropic compressional speed and compare it with recently developed attenuation maps. The pattern emerging in the resulting tomograms is interpreted with recent data on the viscosity of iron as the inner core surface manifestation of a thermally-driven flow, with a positive correlation among compressional speed and attenuation and temperature. Although the outer-core convection controls the heat flux across the inner core boundary, the internally driven inner-core convection is a plausible model that explains a range of observations for the inner core, including distinct anisotropy in the innermost inner core.

Assessing the Availability of Spatial Geothermal Information Using Business Tools and Exploratory Data Analysis

Geothermal energy is a renewable energy source that utilizes heat from the earth's core. Because of that, it is a reliable source of energy demonstrating great potential, which is still used to a small extent. Barriers to the wider use of geothermal energy include high investment costs, location constraints, and the quality of resources at various depths. Opposition from local communities is also a barrier. It results from low environmental awareness of geothermal opportunities, exacerbated by the lack of free access to good-quality geothermal spatial information. Therefore, the purpose of this study is to assess geothermal spatial information availability for selected geoportals using selected business tools and exploratory data analysis. The authors examined geoportals in terms of the information provided, the way it is presented, and the features available. They performed the assessment of their degree of similarity using competitive profile assessment, image graphs, cluster analysis, and affinity analysis. The research confirmed that the available spatial information is characterized by high variability, which indicates the lack of uniform rules for collecting, gathering, storing, and sharing geothermal data. The foregoing leads to information chaos, which can impede investors' investment decisions regarding implementing a geothermal installation on a property.

Electride Formation of HCP-Iron at High Pressure: Unraveling the Origin of the Superionic State of Iron-Rich Compounds in Rocky Planets.

Electride possesses electrons localized at interstitial sites without attracting nuclei. It brings outstanding material properties not only originating from its own loosely bounded characteristics but also serving as a quasiatom, which even chemically interacts with other elemental ions. In elemental metals, electride transitions have been reported in alkali metals where valence electrons can easily gain enough kinetic energy to escape nuclei. However, there are few studies on transition metals. Especially iron, the key element of human technology and geophysics, has not been studied in respect of electride formation. In this study, it is demonstrated that electride formation drives the superionic state in iron hydride under high-pressure conditions of the earth's inner core. The electride stabilizes the iron lattice and provides a pathway for hydrogen diffusion by severing the direct interaction between the metal and the volatile element. The coupling between lattice stability and superionicity is triggered near 100GPa and enhanced at higher pressures. It is shown that the electride-driven superionicity can also be generalized for metal electrides and other rocky planetary cores by providing a fundamental interaction between the electride of the parent metal and doped light elements.

Role of partial stable stratification on fluid flow and heat transfer in rotating thermal convection

The liquid iron core of the Earth undergoes vigorous convection driven by thermal and compositional buoyancy. The dynamics of convective fluid motions and heat transfer in such conditions are determined by background rotation, geometrical symmetry, and thermal interactions across the boundaries. In this study, rotating thermal convection in a horizontal fluid layer is considered to understand the fluid flow characteristics in the Earth's outer core focusing on the regions close to the rotational axis. The effects of a partial stable stratification on fluid flow and heat transfer are investigated to ascertain the physical significance of thermal core–mantle interaction on geomagnetic field generation driven by core fluid motion. It is found that even with non-linear evolution, convective instabilities retain the fundamental characteristics of linear onset modes. Mildly supercritical regimes lead to near laminar flows with the transition to turbulent convection occurring for strongly driven convection around 50–100 times enhanced buoyancy. Axial symmetry breaking and preferential damping of small-scale vortical structures are the hallmark of penetrative convection. Rapid rotation sustains small-scale helical flows in stable regions, a necessary ingredient for the sustenance of Earthlike dipolar magnetic fields. Coherent flow structures for strongly turbulent convection are obtained using reduced-order modeling. The overall total heat transfer is suppressed (up to 25%) due to the stable stratification although convective efficiency is enhanced (up to 30%) in the unstable regions favored by rapid rotation. Flow suppression is overcome under strong buoyancy forces, a relevant dynamical regime for deep-seated dynamo action in the Earth's core.

Thermoelastic Properties of B2‐Type FeSi Under Deep Earth Conditions: Implications for the Compositions of the Ultralow‐Velocity Zones and the Inner Core

AbstractThe CsCl‐type (B2) phase of FeSi (B2‐FeSi) has been proposed as a candidate phase in the ultralow‐velocity zones (ULVZs) at the base of the lower mantle and in the Earth's inner core. However, the elastic properties of B2‐FeSi under relevant conditions remain unclear. Here we determine the density, elastic constants, and velocities of B2‐FeSi at high pressures (90–390 GPa) and temperatures (3,000–6,000 K) relevant to the Earth's lower most mantle and the inner core, using first‐principles molecular dynamics simulations. At the base of the lower mantle, B2‐FeSi shows significantly lower velocities and a higher density than those of the ambient mantle. Mechanical mixing models suggest the presence of ∼27–39 vol% B2‐FeSi in the silicate mantle is consistent with the reduced velocities and the elevated density of ULVZs observed seismically. On the other hand, the hcp‐Fe and B2‐FeSi mixture exhibits higher bulk sound velocity compared to the PREM under inner core conditions. Adding superionic H in the interstitial sites of B2‐FeSi lowers its density but has little effect on the bulk sound velocity of B2‐FeSi, precluding H‐bearing B2‐FeSi as a major component in the Earth's inner core.

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. In this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

Reviews

Earth's Core: Geophysics of a Planet's Deepest Interior, by Vernon F. Cormier, Michael I. Bergman, and Peter L. Olson, 2021 Remote Sensing Digital Image Analysis, by John A. Richards, 2022

High thermal conductivity of the Earth's core and geodynamo

The geomagnetic field is generated by dynamo processes in the Earth's core. This process is supported by the cooling of the planet and the growth of the solid core. It is known that the magnetic field existed long before the appearance of the solid core, although according to available estimates and calculations, this could not have happened. Moreover, according to some models, thermal conductivity in the Earth's core may be three times higher than generally accepted values. In this case, the magnitude of the convective heat flow decreases and the generation of the magnetic field stops. In the above-mentioned models it was assumed that the magnitude of the heat flow at the core – mantle boundary over the entire existence of the Earth linearly decreased by only 15–20 %. The latter is a rough guess. Significantly large changes in heat flow are predicted by models of mantle cooling, in which the growth of convective heat flow in ancient times occurs due to a decrease in the viscosity of the mantle with an increase in the temperature of the substance. Below we consider a model of the cooling of the Earth with a threefold value of the thermal conductivity of the core. It is shown that the use of a combined model of core and mantle cooling makes it possible to significantly increase the heat flow at the core – mantle boundary in ancient times. To slow down the growth of the solid core, a subadiabatic layer is included in the model. As a result, the size of the solid core in the model satisfies seismological observations. The model allows us to obtain a sufficient amount of energy to generate a geomagnetic field, starting from the moment the liquid core appeared and to the present. The appearance of a solid core 2.4 billion years ago, which does not lead to sharp changes in heat flow in the model, is consistent with paleomagnetic observations that do not record changes in the behavior of the magnetic field. The model does not exclude the existence of a multipole magnetic field at the initial stage.

A global timekeeping problem postponed by global warming.

The historical association of time with the rotation of Earth has meant that Coordinated Universal Time (UTC) closely follows this rotation1. Because the rotation rate is not constant, UTC contains discontinuities (leap seconds), which complicates its use in computer networks2. Since 1972, all UTC discontinuities have required that a leap second be added3. Here we show that increased melting of ice in Greenland and Antarctica, measured by satellite gravity4,5, has decreased the angular velocity of Earth more rapidly than before. Removing this effect from the observed angular velocity shows that since 1972, the angular velocity of the liquid core of Earth has been decreasing at a constant rate that has steadily increased the angular velocity of the rest of the Earth. Extrapolating the trends for the core and other relevant phenomena to predict future Earth orientation shows that UTC as now defined will require a negative discontinuity by 2029. This will pose an unprecedented problem for computer network timing and may require changes in UTC to be made earlier than is planned. If polar ice melting had not recently accelerated, this problem would occur 3 years earlier: global warming is already affecting global timekeeping.

On the Inner‐Core Differential‐Rotation (Un)Resolvability From Earthquake Doublets: The Traps of Data Selection

AbstractThe phenomenon of differential rotation of the Earth's inner core relative to the mantle is a subject of interest in geodynamo modeling that has been validated by seismological observations, mainly via the earthquake‐doublets method. Although recent studies converge on the time‐varying differential rotation of the inner core relative to the mantle, favoring a decadal variation, the inferred models significantly differ. Here, considering the data selection, the observed data structure, and the subjective model parameterizations, which we avoid by employing a Bayesian transdimensional approach, we show that the recent best‐fit model by Yang and Song (2023, https://doi.org/10.1038/s41561‐022‐01112‐z) featuring the 70‐year decadal variation is not obtained when all available data are considered. Namely, including only a small number of discarded earthquake doublets (&lt;10%) changes the period of the inner‐core differential rotation fluctuation to 20–30 years. More earthquake‐doublet data are required to address the non‐uniqueness of the inversion problem.

Equation of State of Liquid Fe7C3 and Thermodynamic Modeling of the Liquidus Phase Relations in the Fe‐C System

AbstractWe calculated the pressure‐volume relations of liquids Fe and Fe7C3 up to ∼360 GPa and 3000–8000 K based on the first‐principles molecular dynamics simulations. The results demonstrate that liquid Fe7C3 is similar in compressibility and thermal expansivity to liquid Fe at the Earth's core pressure range. We then obtained a thermodynamic model of the liquidus phase relations in the Fe‐C system at high pressures by using the thermal equation of state (EoS) of liquid Fe7C3 as well as that of liquid Fe3C estimated by interpolation between the volumes of Fe and Fe7C3 calculated in this study. The previously reported eutectic points in the Fe‐Fe3C system are reproduced with interaction parameters that are dependent on temperature but independent on pressure at &gt;∼50 GPa. The melting curve of Fe7C3 should be close to that of Fe, which leads to a change in the eutectic system from Fe‐Fe3C to Fe‐Fe7C3 above 250 GPa as observed by previous experiments. The thermodynamic model also suggests that the solidus and liquidus of Fe3C are relatively low. The Fe‐Fe7C3 liquidus phase relations at the Earth's inner core boundary (ICB) pressure indicate the Fe7C3 inner core is unlikely. While the carbon content in the Earth's outer core may be small, the addition of 1 wt% C to liquid Fe drops the liquidus temperature to crystallize hcp‐Fe by ∼450 K, suggesting that the presence of carbon in the outer core is preferable for a relatively low core‐mantle boundary (CMB) temperature.

Lorenz number and transport properties of Fe: Implications to the thermal conductivity at Earth's core-mantle boundary

Abstract The electrical resistivity (ρ) and thermal conductivity (κ) of the Earth's core compositions are essential parameters for constraining the core's thermal state, the inner core age, and the evolutionary history of the geodynamo. However, controversies persist between experimental and computational results regarding the electronic transport properties (ρ and κ) of the Earth's core. Iron is the major element in the core, and its transport properties under high pressure and high temperature conditions are crucial for understanding the core's thermal state. We measured the ρ values of solid iron using the four-wire van der Pauw method at 300 K and pressures ranging from 3 to 26 GPa within a multi-anvil press. For comparison, we calculated the ρ and κ values of hexagonal close-packed (hcp) iron at conditions of 300-4100 K and 22-136 GPa using the first-principles molecular dynamics (FPMD) method. Our calculations generally align with prior studies, indicating that the electrical resistivity of solid hcp iron at Earth's core-mantle boundary (CMB) conditions is ~76-83 μΩ∙cm. The resistivity of hcp iron changes small as it melts from solid to liquid at pressures from 98 to 134 GPa. The impact of temperature and pressure on the Lorenz numbers of solid hcp iron are investigated according to our calculation results and previous studies. Under the CMB's pressure conditions, the κ of hcp iron initially decreases with increasing temperature and subsequently increases. The electron-electron scattering plays a dominant role at low temperatures and causes the decrease in κ. At high temperatures, the increase of electronic specific heat significantly increases the Lorentz number and κ. Overall, we estimate the κ of solid hcp iron at CMB's condition to be 114 ± 6 W/m/K, slightly lower than the room temperature value of 129 ± 9 W/m/K at the same pressure. Our model shows that a 0-525 km thickness of a thermally stratified layer may exist beneath the Earth's CMB depending on the core's heat flow and thermal conductivity.

Decadal length-of-day and geomagnetic changes imply more complex Earth's core motions

Decadal length-of-day and geomagnetic changes imply more complex Earth's core motions

Earth's “Missing” Chlorine May Be in the Core

AbstractThe budget of volatile halogens in the bulk silicate Earth, fluorine and chlorine, differs distinctly from that of chondritic meteorites. Arguably, the bulk silicate Earth shows a low abundance of Cl, while the F budget is in line with the expected volatility trend. One hypothesis for the missing Cl is its sequestration into Earth's core during planetary segregation, but experimental data on partitioning between silicate and metallic melts are limited in pressure and remain inconclusive. Here we use computational quantum mechanical methods to study F and Cl geochemistry during core–mantle differentiation over a wide pressure and temperature range. Our calculations reveal that with increasing pressure and temperature, chlorine shows an enhanced affinity for iron metal. The Cl metal–silicate partition coefficient increase from −1.89 ± 0.84 at 10 GPa and 3000 K to 1.62 ± 0.80 at 130 GPa and 5000 K, while corresponding calculations for F show minimal variations in metal–silicate partitioning across the pressure–temperature conditions investigated, yielding values between −3.61 ± 0.81 and −3.37 ± 0.80. The shift in shows a transition from lithophile to siderophile behavior. Further calculations on isotopic partitioning show negligible fractionation of Cl isotopes (37Cl/35Cl) between the mantle and core. Our results suggest that metallic iron within Earth's mantle and core may serve as an important Cl reservoir, potentially accounting for up to 40% of Earth's overall Cl inventory.