Fact-checked by Grok 2 weeks ago

Earth's inner core

Earth's inner core is the innermost geologic layer of the planet, consisting of a dense, sphere primarily composed of iron and nickel alloys, with a radius of approximately 1,220 kilometers (759 miles) and temperatures reaching up to 5,400°C (9,800°F). Despite the extreme heat, immense pressure from the overlying layers keeps it , distinguishing it from the surrounding liquid outer core. The inner core constitutes about 15% of Earth's total volume when combined with the outer core, contributing significantly to the planet's overall density and mass. Seismological studies reveal that the inner core exhibits seismic , with compressional wave speeds varying directionally, particularly aligned with axis, and features like an innermost inner core region spanning a few hundred kilometers where crystalline iron structures may differ. It grows slowly at a rate of about 1 millimeter per year as material from the outer core solidifies at the inner core , which has subtle topographic variations of a few kilometers. Lighter elements, such as oxygen and (along with ~5-10% ), make up around 10% of its composition, influencing its physical properties and the broader core-mantle dynamics. The inner core's existence and properties were first inferred in through analysis of seismic waves, which show that shear waves can propagate through it—unlike the liquid outer core—providing key evidence for its solid state. This central layer plays a crucial role in generating , as the convective motion in the adjacent liquid outer core interacts with the solid inner core to sustain the geodynamo process that protects the planet from solar radiation. As of , research highlights hemispherical asymmetries, with the of the upper inner core exhibiting slower compressional wave speeds by about 1%, potential relative to , and structural changes at the inner core's surface over the past two decades, suggesting a more dynamic and potentially less uniformly solid structure. Ongoing seismic observations continue to refine our understanding of its structure, including shear wave velocities and , underscoring the inner core's dynamic over geological time.

Discovery and Historical Context

Initial Hypotheses and Seismic Evidence

In the early , geophysicists began hypothesizing about Earth's deep interior using seismic data from earthquakes. Richard Dixon Oldham, analyzing seismograms from distant earthquakes, identified distinct wave arrivals: primary (P) waves, secondary (S) waves, and surface waves, which allowed him to infer a central core region where S-waves were absent, suggesting a liquid outer core at depths around 2,900 km. This 1906 hypothesis marked the first seismic evidence for a dense, fluid core comprising about one-third of Earth's radius, based on wave transmission patterns and a beyond 103° epicentral distance. Building on this, Danish seismologist Inge Lehmann provided the pivotal evidence for a solid inner core in , through meticulous analysis of P-wave records from northern hemisphere earthquakes, including those in recorded at Scandinavian stations like Godthaab and Umanak. She observed anomalous later-arriving P-waves, termed P' (now known as the PKIKP phase), in the core between 112° and 154° epicentral distance, which could not be explained by a uniform liquid core. Lehmann proposed that these waves resulted from transmission through a solid inner region with increased P-wave velocity (approximately 8.6 km/s compared to 8 km/s in the outer core), acting as a refracting . Further supporting this, Lehmann identified short-period precursors to P', interpreted as reflections (PKiKP phase) from the inner core boundary (ICB), indicating a sharp velocity contrast at the solid-liquid . These PKIKP and PKiKP phases, observed at distances up to 143°, confirmed the existence of a distinct solid inner core boundary at approximately 5,150 km depth. Initial travel-time analyses in yielded an inner core radius estimate of about 1,400 km, refined in the 1940s by Beno Gutenberg and using expanded seismic datasets to approximately 1,220 km.

Evolution of Models (1900s–Present)

Following Inge Lehmann's 1936 discovery of the inner boundary, seismic studies in the and began to probe its physical state, with early indications of low suggesting a rather than liquid composition. By the early , detailed analysis of wave propagation and rigidity moduli provided definitive evidence for the inner 's solidity, overturning prior assumptions of a fully fluid and establishing it as a distinct within the liquid outer . In the 1980s, observations of differential travel times for PKP waves—specifically PKIKP phases traversing the inner core—revealed systematic deviations that could not be explained by isotropic models, leading to the confirmation of large-scale seismic aligned with Earth's rotational . This breakthrough, based on global data, implied directional variations in wave speeds, with faster propagation along polar paths than equatorial ones, prompting the development of cylindrically symmetric models for the inner core. The 1990s and 2000s saw significant advancements through the integration of expansive global seismic networks, such as those coordinated by the , which deployed dense arrays of broadband stations to capture high-resolution inner core signals from distant earthquakes. These datasets facilitated refinements to reference Earth models like PREM, including updates to inner core boundary radii and velocity gradients in models such as iasp91 (1991) and ak135 (1995), enhancing accuracy in delineating the core's interfaces and internal structure. Entering the 21st century, seismic analyses uncovered evidence of , with the inner core exhibiting super-rotation relative to at rates of about 0.3 to 0.5 degrees per year from 1996 to around 2009, as inferred from temporal changes in PKP travel times and patterns. This super-rotation, first proposed by Song and Richards in 1996, was linked to electromagnetic coupling at the core-mantle boundary. However, post-2010 observations indicated a deceleration and potential reversal, with the inner core's slowing to match or lag behind the surface by 2010, as confirmed by a 2024 () study analyzing repeating doublets and misfits. By 2025, models incorporating reversals and patterns have revealed dynamic shape changes at the inner 's surface, including localized undulations and structural transformations possibly driven by outer , challenging the notion of a static spherical and suggesting ongoing deformation over decadal timescales. These findings, derived from high-fidelity seismic arrays monitoring South Sandwich Islands events, imply that the inner 's boundary may exhibit topographic variations of several kilometers, prompting revised geodynamical simulations of evolution. Additionally, September 2025 research provided experimental evidence that the inner may exist in a superionic state, where light elements like carbon are mobile within a solid iron , explaining observed softening and ultralow wave velocities.

Methods of Study

Seismic Wave Analysis

Seismic waves generated by earthquakes provide the primary means to probe the Earth's inner core, as these waves traverse the planet and are recorded by global networks. Compressional P-waves, which propagate through both solids and liquids by alternating and , and shear S-waves, which involve transverse particle motion and thus only travel through solids, both traverse the inner core. The propagation of S-waves through the inner core, in contrast to their absence in the liquid outer core, confirms the inner core's solidity, a conclusion solidified by observations in the through detection of inner core phases. Key observations of inner core structure derive from travel-time anomalies in PKP phases, which are P-waves that traverse the outer and refract through the inner . These phases exhibit triplications—multiple arrivals at certain epicentral distances—due to the velocity contrast at the inner boundary (ICB), allowing estimation of the inner radius at approximately 1,220 km. For instance, differential travel times of PKP branches (e.g., DF, bc, and ab) reveal anomalies of up to several seconds, particularly for paths sampling the uppermost inner , indicating lateral variations in . Additionally, splitting of J-waves (SKJ phases), which are shear waves grazing the ICB, provides evidence of , as the waves arrive with split polarizations and time delays, suggesting directional dependence in wave speeds aligned with the axis. In the Preliminary Reference Earth Model (PREM), the average P-wave velocity in the inner core is approximately 11.0 km/s, increasing radially inward due to compression. The inner core radius is inferred from the triplication in PKP waveforms, where ray paths refract at the ICB and produce overlapping arrivals for epicentral distances around 140°–180°. Ray path curvature in PREM is described by the relation for travel-time derivative with respect to epicentral distance: \frac{dt}{d\Delta} = \frac{r \sin i}{v} where r is the radius, v is the wave speed, and i is the ray incidence angle at the boundary; this equation aids in modeling how waves bend through the core, constraining boundary sharpness and velocity gradients. Modern techniques enhance detection of low-amplitude core phases, such as through , which uses seismic arrays to beamform signals and suppress noise, enabling observation of faint inner core shear waves that individual stations cannot resolve. For example, on array data extracts PKP precursors by correlating ambient noise or waves, revealing fine-scale heterogeneity near the ICB. Recent studies leveraging repeating earthquakes—events with similar source locations—have detected temporal changes in inner core wave speeds, with waveform mismatches indicating sub-rotation and structural evolution over decades, as documented in analyses up to 2023 and extended in 2025 observations of shape variations.

Laboratory Simulations and Other Techniques

Laboratory simulations of Earth's inner core conditions primarily rely on high-pressure experiments using diamond anvil cells (DACs) to replicate the extreme pressures of up to 360 GPa at the inner core boundary. These experiments compress iron samples to investigate phase transitions, such as the shift from body-centered cubic (bcc) to hexagonal close-packed (hcp) structures, which is observed above approximately 15 GPa under ambient temperatures. techniques integrated with DACs have mapped the γ (face-centered cubic) to ε (hcp) transition boundary up to 69 GPa, providing insights into the stability of hcp iron dominant in the inner core. Recent DAC studies at inner core pressures confirm that while bcc iron remains mechanically stable, hcp is thermodynamically favored, influencing seismic anisotropy interpretations. Additionally, DAC experiments measuring sound velocities in hcp iron up to core conditions reveal collective atomic motions that affect elasticity. Neutrino-based techniques offer an independent probe of the inner core's electron density profile, unaffected by seismic wave assumptions. Large detectors like Borexino and KamLAND, designed for low-energy detection, enable through effects influenced by matter density. A seminal method uses oscillations to remotely measure electron densities, potentially resolving the inner core's by distinguishing element contributions in the core-mantle boundary region. Geoneutrino detections from and decays, observed at 4σ significance by both experiments, constrain the overall heat budget and indirectly limit radioactive elements in the core, supporting iron-dominated models with minimal alloys. Geomagnetic and paleomagnetic records of field reversals serve as indirect indicators of inner dynamics, revealing interactions between the solid inner and the fluid outer dynamo. Paleomagnetic data from reversals, documented in the Geomagnetic Polarity Time Scale spanning 170 million years, show asymmetric field behaviors linked to inner growth and convection patterns. Numerical models incorporating heterogeneity demonstrate that regional variations trigger reversals, with inner anisotropy modulating the time-averaged paleofield. These reversals occur on timescales of centuries to millennia, constraining flow velocities and electromagnetic coupling. Advancements in 2025 have refined understandings of inner formation and structure through computational and methods. -informed simulations indicate that carbon concentrations of about 10 mol% (reducing to ~481 K) stabilize during solidification, resolving the paradox of inner onset under realistic cooling rates. These models, validated against data, suggest up to 15 mol% carbon enables at 330–360 GPa, implying higher carbon abundance than previously estimated. Complementing this, ultrasonic experiments on high-pressure hcp-structured samples mimic inner anisotropy by measuring premelting effects. Pulse-echo ultrasonics at 3 GPa near melting temperatures show velocity drops of 14.5% and Poisson's ratio increases to ~0.444, aligning with observed low inner velocities and confirming premelting's role in elastic properties.

Physical Properties

Size, Shape, and Volume

The Earth's inner core is a solid primarily composed of iron and nickel, with a of approximately 1,220 km, extending from the center of the planet outward to the inner core boundary (ICB) at a depth of about 5,150 km from the surface. This places the inner core at roughly 19% of Earth's total of 6,371 km. Early geophysical models, based on propagation and the (PREM), assumed a highly due to the symmetric nature of observed P-wave and S-wave velocities across the ICB. The of the inner core can be calculated using the formula for the volume of a , V = \frac{4}{3} \pi r^3, where r is the . Substituting r \approx 1,220 km yields a of approximately $7.6 \times 10^9 km³, representing about 0.7% of Earth's total of $1.083 \times 10^{12} km³. This modest volumetric fraction underscores the inner core's concentrated in planetary , with its implications depending on densities around 12–13 g/cm³ that amplify its gravitational influence despite the small size. While traditionally modeled as perfectly spherical, recent seismic analyses indicate deviations from ideal , including potential ellipticity with less than 0.5% deviation from a perfect , possibly arising from rotational forces and gravitational coupling with the overlying . Furthermore, from waveform comparisons of earthquakes recorded between 2004 and 2024 reveals structural changes at the ICB, suggesting the inner core's surface has deformed by up to a few kilometers in height over the past two decades, challenging assumptions of static geometry. These shifts, detected through variations in seismic signal shapes propagating through , imply dynamic interactions with the liquid outer core, though the overall form remains dominantly spherical.

Pressure, Density, and Mass

The within Earth's inner arises from the , governed by the equation \nabla P = \rho \mathbf{g}, where the balances the gravitational force per unit volume, and \mathbf{g} decreases toward the center due to the diminishing in a spherically symmetric distribution. At the inner boundary (ICB), pressure reaches approximately 330 GPa, increasing slightly to about 365 GPa at the center, reflecting the integrated weight of overlying material across the relatively thin inner core layer. The density profile of the inner core, as modeled in the (PREM), shows an average value of 12.8–13.1 g/cm³, with a gradual increase to roughly 13.1 g/cm³ at the center due to compressional effects under extreme pressure. A notable discontinuity occurs at the ICB, where density jumps by about 0.6 g/cm³ (from ~12.2 g/cm³ in the outer core to ~12.8 g/cm³ in the inner core), attributed to the depletion of light elements during solidification that preferentially remain in the liquid outer core. The total mass of the inner core is approximately $9.9 \times 10^{22} kg, constituting about 1.7% of Earth's overall mass, consistent with its high average density and spherical volume of radius ~1,220 km. This mass distribution supports the light element explanation for the ICB density contrast, as the solidification process enriches the solid inner core in heavier iron-nickel alloys while excluding lighter components like silicon or hydrogen. Recent 2025 seismic analyses indicate structural transformations near the inner core's surface, suggesting a less rigidly solid state with potential local variations due to deformation under gravitational and electromagnetic forces, which may refine models of its overall profile.

Temperature and Thermal State

The thermal regime of Earth's inner core is dominated by extreme temperatures that increase gradually from the inner core boundary (ICB) toward , reflecting its nearly adiabatic structure. Estimates place the at between 5,400 and 6,000 , while at the ICB it is approximately 5,700 , comparable to the surface of . The profile follows an adiabatic gradient, described by the \frac{dT}{dr} = \frac{\alpha T g}{C_p}, where \alpha is the thermal expansivity, T is , g is , and C_p is the at constant pressure; this results in a small increase of only tens of across the inner core's radius. The inner core remains solid despite these high temperatures due to the elevated pressures, with its boundary at the temperature of the iron alloy comprising the core. The curve of iron under core pressures (330–360 GPa) is extrapolated using approximations from the Lindemann law, which posits when vibrations reach a critical ; these models indicate a of approximately 6,000 for pure iron at inner core conditions. Recent experimental measurements using laser-driven shock compression and confirm this range, with values up to 6,200 as an upper limit. Heat contributing to the inner core's state originates primarily from released during ongoing solidification at the ICB, which powers in the overlying outer core, and from potential radiogenic decay of trace incorporated into the core during Earth's formation. A 2025 study highlights the role of carbon as an impurity, demonstrating that even small concentrations (around 2–3 wt%) lower the by facilitating under realistic , thus influencing the core's evolution without requiring excessive cooling. At the ICB, a thin boundary layer in the outer core exhibits a temperature jump of 100–200 relative to the adiabatic profile, accommodating the transition to the solid inner core.

Viscosity and Mechanical Behavior

The viscosity of Earth's inner core is estimated to range from $10^{18} to $10^{24} Pa·s, reflecting its behavior as a highly viscous solid that acts elastically on short seismic timescales (seconds to minutes) but can exhibit plastic deformation over longer geological periods (millions of years). This extreme viscosity arises from the material's composition, primarily hexagonal close-packed (hcp) iron under immense pressure and temperature, which resists flow but allows gradual creep under sustained stress. Such properties imply that the inner core maintains structural integrity during rapid seismic wave propagation while permitting slow adjustments to align with external forces, contributing to observed seismic anisotropy through mechanical reorientation of crystals. Creep in the inner core is dominated by dislocation glide within the hcp-iron lattice, a mechanism where defects in the crystal structure enable deformation at low strain rates typical of the core's dynamics ($10^{-14} to $10^{-18} s^{-1}). This process, akin to power-law dislocation creep, allows the material to flow viscously under deviatoric stresses on the order of tens of pascals, contrasting with diffusion creep that may prevail near the inner core boundary. The viscoelastic response can be characterized by the Maxwell relaxation time \tau = \eta / \mu, where \eta is the viscosity and \mu is the shear modulus, approximately 300 GPa for hcp iron at inner core conditions; for viscosities in the $10^{18}–$10^{24} Pa·s range, \tau spans approximately $10^6 to $10^{13} seconds, far exceeding seismic periods but permitting deformation over core ages. Post-seismic deformation models, incorporating viscoelastic relaxation of the , support these low effective rigidity estimates by simulating how the responds to or electromagnetic stresses, with below $10^{18} Pa·s enabling observable alignments and attenuations in seismic data. Recent 2025 experimental evidence indicates a less rigidly solid state than previously assumed, attributed to superionic phases in Fe-C alloys where carbon ions diffuse freely within the iron , implying an even lower effective and enhanced softening at the .

Composition and Phase

Primary Constituents

The Earth's inner core is primarily composed of an , with iron constituting approximately 85–90% by weight and making up 5–10% by weight. This dominant metallic is inferred from cosmochemical models of , which draw on the siderophile element partitioning observed in iron meteorites, where the Ni/Fe ratio closely matches that expected for the core. In addition to these major elements, the inner core incorporates light elements such as carbon (C), (S), oxygen (O), and (Si) at concentrations approximately 10% by weight in total, which account for the observed density deficit relative to pure iron under core pressures. At the extreme pressures exceeding 330 GPa in the inner core, the predominantly adopts the hexagonal close-packed (hcp) ε-Fe phase, which is stable for iron-nickel mixtures under these conditions. However, trace inclusions of the face-centered cubic (fcc) phase may exist locally, potentially stabilized by or other impurities at high temperatures near 5000–6000 . A significant advancement in 2025 revealed that carbon serves as a key in the inner core, present at concentrations of approximately 3.8% by weight, which lowers the supercooling required for and enhances the stability of solidification processes. This finding, derived from simulations of Fe-C alloys, indicates that even this carbon content enables homogeneous at the inner core boundary, resolving longstanding paradoxes in core formation models.

Solidification Processes and Carbon's Role

The inner core of grows primarily through fractional crystallization at the inner core boundary (ICB), where iron-rich alloys solidify from the overlying liquid outer core. This process involves the progressive freezing of the core melt as the planet cools, with lighter elements partitioning into the remaining , thereby driving compositional in the outer core. The solidification releases significant , estimated to contribute substantially to the core's thermal budget and the powering of the geodynamo. Over geological timescales, this growth induces a temperature decrease at the ICB of approximately 100 per gigayear (Gyr), reflecting the slow cooling of the core-mantle system. Dendritic growth characterizes the solidification at the ICB, arising from the morphological of the solid-liquid under constitutional conditions. As the interface advances, solute rejection ahead of the advancing front creates a solute-enriched in the melt, leading to local undercooling that destabilizes the planar front and promotes branching dendritic structures. This results in the formation of a thin mushy zone at the top of the inner , consisting of interconnected solid dendrites and interdendritic pockets. Models indicate that the mush fraction in this zone is on the order of 20–30%, influencing permeability and fluid flow at the boundary. Recent research highlights carbon's critical role in facilitating inner core solidification, as revealed in 2025 studies of Fe-C nucleation. Carbon acts as a alloying element that depresses the of iron by 500–1,000 K depending on concentration, enabling the core to reach the nucleation threshold with modest of 266–481 K. This depression stabilizes the initial solid against disruptive convective flows in the outer core, preventing premature and allowing sustained growth. Without sufficient carbon (estimated at ~3.8 wt% in the core), the required would exceed geophysical limits, potentially inhibiting inner core formation altogether. In dendritic growth models incorporating constitutional , the interface advance rate v can be approximated by the relation derived from solute and undercooling : v = \frac{\Delta T}{m C_0 K} where \Delta T is the undercooling, m is the liquidus , C_0 is the initial solute concentration in the melt, and K is the . This equation underscores how solute effects control the growth kinetics at the ICB, with carbon influencing m and C_0 to promote solidification.

Internal Structure

Seismic Anisotropy

Seismic in Earth's inner core manifests as directional variations in the propagation speeds of seismic , primarily compressional P-waves, due to the ordered alignment of crystalline structures within the solid iron-nickel . This results in faster wave travel along certain paths relative to others, providing key insights into the inner core's fabric and deformation history. The dominant form is axial , characterized by cylindrical aligned with Earth's axis, where V_p is approximately 3–4% faster parallel to the polar axis than in equatorial directions. Non-axial variations, including hemispherical differences, further complicate this pattern, with stronger in the compared to the eastern, reflecting regional asymmetries in alignment. Measurements of this rely on analyzing travel times and splitting of seismic phases that traverse the inner . For instance, wave splitting in SKJKS phases, which involve waves (J phases) propagating through the inner , shows delays (Δt) of about 0.5–1 s, indicative of the velocity contrast between quasi-longitudinal and quasi-transverse components. Transverse wave (V_s) are estimated at around 3.5 km/s in equatorial directions, with faster propagation at oblique angles to the rotation axis. These observations are derived from seismic arrays and modeling of repeating earthquakes. The underlying cause of inner core seismic anisotropy is the development of crystallographic preferred orientation (CPO) in hexagonal close-packed (hcp) iron , which align preferentially during deformation. This alignment occurs primarily through , a deformation mechanism driven by differential stresses from , inner core rotation, or interactions at the core-mantle boundary. in the inner core, estimated at around 10^{21} ·s, facilitates this process, enabling crystals to reorient over geological timescales without fracturing. Recent observations as of 2025 indicate temporal changes in inner core , linked to disturbances that alter the core's shape, such as viscous deformations at its surface induced by outer core . These shifts, detected via variations in seismic waveform residuals over decades, suggest dynamic adjustments in fabric that could influence global and generation.

Layering and Substructures

The Earth's inner core displays radial stratification, consisting of multiple distinct substructures inferred from seismic wave analyses. The outermost portion features an isotropic layer approximately 60–100 km thick, characterized by relatively uniform seismic properties and potentially mushy textures due to or compositional heterogeneity at the inner core boundary. This region, sometimes termed the F-layer, transitions into a broader anisotropic zone where crystalline alignment imparts directional variations in wave propagation. Further inward lies an anisotropic transition layer, referred to as the E-layer, marking the shift to more pronounced structural differences toward the core's center. At the innermost region, the innermost inner core (IMIC) forms a compact with a of roughly 300–750 km, exhibiting distinct seismic signatures including elevated P-wave velocities (V_p) by 1–2% relative to the overlying inner core material. This substructure arises from differences in iron crystallinity, potentially involving a face-centered cubic (fcc) of iron under extreme pressure-temperature conditions. A 2025 study analyzing seismic data provided confirmatory evidence for this innermost layer, revealing reverberating seismic waves that highlight its anisotropic distinction and support the presence of phase-related boundaries within the inner core, including and confirming a transition in elastic structure. These substructures underscore the complex solidification history of the inner core, as captured in high-resolution models from recent thermodynamic simulations.

Lateral and Temporal Variations

The Earth's inner core exhibits notable hemispherical asymmetry in its seismic properties, with the displaying stronger P-wave compared to the . Seismic observations indicate that in the reaches values of 3–4% or higher, increasing with depth, while the shows weaker of approximately 0.5–1.5%. This difference, characterized by a ΔV_p variation of about 1–2% between hemispheres, is attributed to differential growth rates influenced by outer core patterns that favor faster solidification in the eastern sector. Recent studies have revealed temporal variations in the inner core's structure, including shape transformations over decades driven by interactions with outer core flows. A 2025 analysis of seismic data from repeating earthquakes documented non-uniform changes in the inner core's surface morphology, with evidence of viscous deformation and topographic undulations forming over periods of years to decades. These alterations are linked to convective disturbances in the outer core, which exert traction on the inner core , causing localized bulging and flattening. Additionally, seismic waveform comparisons suggest a slowdown in the inner core's rotation relative to since around 2010, manifesting as a backward of approximately 0.1° per year. Evidence for these temporal changes comes primarily from arrays of repeating earthquakes, where seismic traversing the inner core show systematic shifts in travel times. For instance, differential traveltimes for PKIKP have increased by about 0.3 seconds over intervals spanning 20–30 years, corresponding to a relative wave speed decrease of roughly 0.3% per in certain paths. These observations, derived from global doublets recorded between 1991 and 2023, highlight ongoing structural rather than static heterogeneity. Such lateral and temporal variations imply a highly dynamic inner core surface, featuring "hills and valleys" with amplitudes up to several kilometers, sculpted by convective disturbances from the outer core. This influences and material exchange at the inner , potentially affecting the geodynamo's and long-term core . While the overall radial remains intact, these surface features underscore the inner core's responsiveness to in the surrounding outer core.

Dynamics and Motion

Rotational Dynamics

The Earth's inner core exhibits relative to the overlying , characterized by variations in its over decadal timescales. Seismic observations using doublets—pairs of similar earthquakes recorded at the same station—have revealed that from the to around 2009, the inner core underwent a phase of super-rotation, spinning eastward faster than at rates of approximately 0.3° to 0.5° per year. This super-rotation was inferred from systematic changes in the travel times of seismic waves traversing the inner core, particularly PKP waves, which showed temporal misalignments consistent with the inner core's eastward drift. Post-2010, seismic analyses indicate a marked deceleration of this super-rotation, transitioning to a slower westward sub-rotation relative to at rates on the order of 0.05° to 0.15° per year. Studies published in 2023 and 2024 using repeating multiplets documented this shift, with evidence of reversals suggesting the inner began backtracking around 2008–2009, leading to a slow westward sub-rotation as observed through 2023, accompanied by annual-scale variations. These variations are part of a broader multidecadal , potentially with a period of about 70 years, driven by interactions at the -mantle boundary. A 2025 study using seismic data up to 2023 has further revealed annual-scale variability in the rotation rate and non-rotational changes near the inner surface, likely due to viscous deformation influenced by outer flows and interactions. The is governed by the relation \omega_{ic} = \omega_m + \Delta \omega, where \omega_{ic} is the inner core's , \omega_m is 's , and \Delta \omega represents the relative differential component arising from and magnetic alignments. Key coupling mechanisms include electromagnetic torques generated by fluid flows in the outer core, which interact with the geomagnetic field to drive the inner core's motion, and gravitational interactions with that exert a stabilizing pull through heterogeneities at the core-mantle boundary. These torques balance to produce the observed oscillatory behavior, with electromagnetic effects dominating short-term variations and gravitational coupling influencing longer-term alignments. Shape changes at the inner core's surface may contribute to these rotational oscillations by altering gravitational and viscous interactions with the outer core.

Growth and Boundary Interactions

The Earth's inner core grows through the gradual solidification of molten iron from the surrounding outer core, with an estimated radial rate of approximately 0.5–1 mm per year. This slow accretion has resulted in the inner core accumulating about 1% of Earth's total mass over the planet's 4.5 billion-year history. Growth is asymmetric, occurring faster in the due to variations in influenced by outer core dynamics and hemispherical differences in seismic , while maintains an overall spherical shape. The rate of mass increase for the inner core can be described by the formula \frac{\partial M}{\partial t} = 4\pi r^2 \rho v, where M is the mass, r is the radius, \rho is the density of the solidifying material, and v is the crystallization velocity at the inner core boundary (ICB). This expression quantifies the flux of material freezing onto the inner core surface, balancing heat loss from the core with latent heat release during solidification. At the ICB, interactions with the outer core involve convective plumes that can erode the boundary through localized melting, potentially destabilizing a mushy zone—a hybrid layer of solid crystals and residual liquid—estimated to be 4–8 km thick in certain regions. Seismic evidence indicates that this mushy zone's stability is influenced by outer core convection, which drives phase changes and maintains a dynamic interface despite the overall growth. Recent observations from 2025 reveal that disturbances in the outer core, such as turbulent flows, induce structural transformations in the inner , rendering it less rigidly than previously assumed and altering its surface shape on observable timescales. These changes, detected via seismic of repeating earthquakes, highlight the inner 's viscous response to outer core activity, with implications for its overall and geodynamic evolution.

Influence on Earth's Magnetic Field

The Earth's inner core plays a crucial role in the geodynamo process that generates the , primarily through its high electrical conductivity and interaction with the surrounding fluid outer core. With an electrical conductivity of approximately $10^6 S/m, the solid inner core enables the frozen-flux approximation in numerical geodynamo models, where magnetic field lines are effectively "frozen" into the conducting material and advected by the inner core's motion relative to the outer core flow. This mechanism allows the inner core to contribute to the field's large-scale structure by resisting rapid diffusion of , thereby supporting the maintenance of a coherent dipole-dominated field observed at Earth's surface. Seismic in the inner , characterized by faster wave along polar directions compared to equatorial ones, aligns with the geomagnetic field's axial geometry and helps stabilize its . This arises from the preferred of iron crystals in the inner , which couples electromagnetically with the outer to dampen instabilities that could lead to field reversals. The ongoing growth of the inner further enhances this stability by releasing and compositional buoyancy—primarily from light elements like or oxygen rejected during solidification—which drives vigorous in the outer and sustains the geodynamo's power. Paleointensity records link the inner 's to the of the geomagnetic , with evidence indicating a as early as 3.4–3.45 billion years ago in rocks. This persistence suggests that inner formation provided a critical boost to the geodynamo by introducing compositional , enabling the field to endure beyond the initial thermal cooling phase of Earth's . Recent observations as of 2025 reveal a slowdown in the inner 's relative to , decelerating since around 2010 and potentially reversing direction in a ~70-year . This deceleration may influence the geodynamo by altering electromagnetic coupling with the outer , potentially affecting the geomagnetic field's strength over decadal timescales.

Age and Formation

Thermodynamic Models

Thermodynamic models of Earth's inner formation rely on and balances within to estimate its , integrating heat loss from core cooling with equilibria at the inner (ICB). The solidification process releases as liquid iron alloys crystallize, quantified by the budget Q_{\text{lat}} = L \Delta M, where L is the of fusion (approximately 300 kJ/kg for iron under core conditions) and \Delta M is the solidified since . This heat contributes to powering outer core and the geodynamo, with models indicating that the inner likely nucleated less than 1 billion years after core formation to balance the total released, assuming no significant radioactive heating in the . equilibria further constrain the process, as the inner 's growth depends on the temperature drop below the iron alloy's at ICB pressures (around 330–360 GPa), driving gradual solidification over time. Adiabatic cooling in the outer core imposes a lower bound on the inner core's age through entropy balance considerations. The outer core's , primarily from ohmic in the (350–700 MW K⁻¹), limits how rapidly the core can cool without violating superadiabatic requirements. Models incorporating these balances predict a minimum inner core age of approximately 500–600 million years, as faster cooling would exceed the entropy available to sustain observed geomagnetic strengths. This constraint arises from integrating the core's equation, where the at the ICB must align with adiabatic gradients to prevent stable that could inhibit . Recent 2025 models incorporating carbon as a light-element in highlight how such effects can delay but support an inner core age of approximately 1 billion years. Carbon reduces the needed for (to ~200–400 K at 1–15 mol%) by depressing the melting temperature of iron alloys, enabling realistic thermal histories consistent with paleomagnetic records of an active . Recent 2025 modeling further constrains the age to 0.5–2 Gyr, with carbon enabling with minimal (~266 K at ~15 mol% C), resolving paradoxes for ages within the last ~1 Gyr without excessive delays. Core-mantle boundary (CMB) heat flux provides additional constraints, requiring 5–15 TW to match observed inner core growth rates and sustain . This flux drives core cooling, with inner core solidification accounting for 20–50% of it through and compositional , ensuring the thermal budget aligns with the core's present (about 1220 km). Variations in flux, influenced by , directly impact growth models, reinforcing age estimates from 0.5–2 Gyr.

Paleomagnetic and Geochemical Evidence

Paleomagnetic records indicate that virtual geomagnetic pole (VGP) paths from volcanic and sedimentary rocks show increased clustering near the geographic poles beginning approximately 1 billion years ago, signifying the onset of a more stable axial dipole field driven by the solid inner 's influence on outer convection. This enhanced dipolarity reduced field variability and secular variation compared to earlier periods, where VGPs exhibited greater scatter indicative of a less organized regime. Numerical geodynamo models link this stabilization to inner nucleation around 1 Ga, as the growing solid provided electromagnetic that damped nondipolar components and promoted axial symmetry in the magnetic field. Supporting evidence against an early inner core formation comes from high paleointensity measurements in rocks, such as those from late (ca. 2.77 Ga) volcanics in the , , which record field strengths of approximately 85 μT—comparable to or stronger than modern values. These robust intensities imply a vigorous geodynamo powered solely by thermal in a fully liquid core, without the additional compositional buoyancy from inner core crystallization. Similarly, single-zircon paleointensities from 3.4–4.2 Ga detrital grains confirm sustained high field strengths over 1 billion years, challenging models of inner core nucleation before 3 Ga. Geochemical analyses indicate Os isotope heterogeneity in the mantle consistent with late siderophile element exchange, though timing remains debated with major events around 4.5–4.4 Ga and potential prolonged interactions. Peridotite xenoliths from Archean cratons display unradiogenic ¹⁸⁷Os/¹⁸⁸Os values, reflecting Re/Os fractionation in the mantle, but direct ties to inner core growth are uncertain. A 2025 study of ruthenium (Ru) isotopes in ocean island basalts from Hawaii identifies ε¹⁰⁰Ru anomalies (0.09 ± 0.03) in the mantle source, interpreted as signatures of core material leakage through the core-mantle boundary. Combined with anomalous tungsten isotopes, these data indicate ongoing core-mantle exchange, highlighting dynamic boundary processes that postdate early core formation. This empirical evidence complements thermodynamic models by underscoring active interactions at the core-mantle interface.