Gravitational Torque & Multidecadal Day Length Shifts
Gravitational Torque and Multidecadal Variations in Length of Day
1. Introduction to Earth’s Variable Rotation
1.1 Definition of Length of Day (LOD) and Rotational Fluctuations
Length of Day (LOD) measures the duration required for Earth to complete one full rotation on its axis relative to celestial reference points. Astronomically, LOD is defined relative to the International Celestial Reference System and differs from the standardized metric:
$$\text{Standard Day} = 86,400 \text{ SI seconds}$$
Rotational variations express as deviations from this nominal baseline:
$$\Delta\text{LOD} = \text{True LOD} - 86,400\text{ s}$$
Positive values of $\Delta\text{LOD}$ denote an increase in rotation period (planetary deceleration). Negative values indicate accelerated planetary spin.
Earth’s rotational velocity exhibits variations across distinct temporal domains:
- Short-Term Variations (Sub-daily to Interannual): Driven by atmospheric angular momentum (AAM), oceanic circulation, zonal wind patterns, and seasonal mass shifts. Amplitudes range from tenths of a millisecond to approximately 1 millisecond.
- Decadal Variations (10–30 Years): Associated with torsional oscillations in the fluid outer core, large-scale geomagnetic jerks, and intermediate climate oscillations.
- Multidecadal Variations (50–70 Years): Dominated by a pronounced ~60-year periodic fluctuation with amplitudes reaching 3 to 5 milliseconds.
- Secular Trends (Centennial to Deep Geological Time): Characterized by a steady slowing of approximately 1.5 to 2.3 milliseconds per century, driven by lunar and solar tidal dissipation.
+-------------------------------------------------------------------------------+
| Temporal Spectrum of Earth Rotation Variations |
+-------------------+-----------------------------+-----------------------------+
| Scale | Dominant Period | Primary Physical Driver |
+-------------------+-----------------------------+-----------------------------+
| Short-term | Days to 1–2 years | Atmosphere & Ocean Dynamics |
| Decadal | 10 to 30 years | Outer Core Torsional Waves |
| Multidecadal | 50 to 70 years (~60 years) | Gravitational Torque (IC-M) |
| Secular | Centennial / Millennial | Tidal Friction & GIA |
+-------------------+-----------------------------+-----------------------------+
The ~60-year multidecadal LOD oscillation represents the largest amplitude signal in non-secular Earth rotation changes.
1.2 The Scope of Deep Earth Dynamics vs. Surface Drivers
Rotational variations stem from two primary categories: external surface processes and internal planetary dynamics.
Surface processes operate via mass redistribution and angular momentum conservation within the surficial fluid envelopes. Atmospheric winds, Antarctic circumpolar current shifts, global precipitation cycles, and post-glacial rebound (Glacial Isostatic Adjustment, or GIA) continuously alter the planet’s principal moment of inertia ($C$) or exchange angular momentum with the solid lithosphere. These surface mechanisms cannot account for the amplitude or the persistent ~60-year periodicity observed in historical geodetic records. The total angular momentum capacity of the atmosphere and ocean is insufficient to sustain millisecond-scale variations over multiple decades.
Total Angular Momentum (L_total) = L_mantle + L_outer_core + L_inner_core + L_atmosphere + L_oceans = Constant
Deep Earth processes account for this deficit. The liquid outer core and solid inner core contain approximately 32% of Earth’s total mass and possess a substantial fraction of its rotational inertia. The Core-Mantle Boundary (CMB), situated at a depth of roughly 2,890 kilometers, and the Inner Core Boundary (ICB), at 5,150 kilometers, act as dynamic mechanical interfaces. Exchanges of angular momentum across these internal boundaries control long-period planetary rotational variations.
2. Theoretical Framework of Gravitational Torque
2.1 Core-Mantle Coupling Mechanisms
Angular momentum exchange across the CMB requires physical coupling mechanisms capable of exerting net torque between the liquid/solid core and the silicate mantle:
- Electromagnetic Coupling: Flowing iron-nickel fluid across the CMB interacts with the geomagnetic field, inducing electrical currents in the weakly conducting lowermost mantle (the $D’’$ layer). Lorentz forces ($\mathbf{J} \times \mathbf{B}$) generate shear stresses that transfer torque. While effective for sub-decadal and decadal fluctuations, electromagnetic coupling alone lacks the magnitude to sustain the ~60-year multidecadal amplitude without requiring unphysically high mantle electrical conductivity.
- Topographic Coupling: Non-hydrostatic boundary variations (topographic undulations of 1–3 km amplitude on the CMB) obstruct azimuthal fluid flows. Dynamic core pressure acting against these slopes generates pressure torques. Theoretical models show that while topographic torque contributes to broadband core-mantle exchange, it is prone to self-damping and turbulent dissipation.
- Viscous Coupling: Molecular and turbulent viscosity in the thin boundary layer at the CMB transfers momentum via laminar shear. The viscosity of liquid iron under CMB pressures is too low ($\sim 10^{-3} \text{ Pa}\cdot\text{s}$) to provide significant global torque.
- Gravitational Coupling: Asymmetric internal density structures within the solid mantle generate a non-axisymmetric gravitational potential. This potential exerts a non-contact, long-range volume torque on matching non-hydrostatic density anomalies within the solid inner core. Gravitational coupling provides the primary restoring force responsible for driving multidecadal rotational exchanges.
+-----------------------------------+
| Solid Silicate Mantle |
| (Density Anomalies / LLVPs) |
+-----------------+-----------------+
|
Gravitational Torque
|
+-----------------v-----------------+
| Solid Inner Core |
| (Triaxial Mass Anomalies) |
+-----------------+-----------------+
|
Viscous Relaxation /
MHD Convective Coupling
|
+-----------------v-----------------+
| Liquid Outer Core |
| (Geodynamo Advection Flow) |
+-----------------------------------+
2.2 Gravitational Interaction Between Mantle and Inner Core
The interior of the Earth is not spherically symmetric. Seismic tomography reveals large-scale heterogeneous structures in the deep mantle, termed Large Low Velocity Provinces (LLVPs), located beneath Africa and the Pacific Ocean. These features, along with subducted slab remnants, represent significant density variations ($\delta \rho_m$).
Similarly, the solid inner core exhibits non-hydrostatic triaxiality and density variations ($\delta \rho_{ic}$) frozen into its crystalline iron-nickel lattice or induced by dynamic core solidification patterns.
The gravitational potential field $V_m(\mathbf{r})$ generated by the lower mantle’s mass anomalies penetrates the fluid core and interacts directly with the density field $\rho_{ic}(\mathbf{r})$ of the inner core. The resulting gravitational torque $\boldsymbol{\Gamma}_g$ acting on the inner core is defined by:
$$\boldsymbol{\Gamma}g = \int{V_{ic}} \rho_{ic}(\mathbf{r}) \left( \mathbf{r} \times \nabla V_m(\mathbf{r}) \right) dV$$
When the inner core rotates relative to the mantle, its non-axisymmetric density anomalies misalign with the corresponding mantle gravitational potential wells. This misalignment creates a restoring torque proportional to the angle of misorientation $\theta$:
$$\Gamma_z = -\kappa (\theta_{ic} - \theta_m)$$
Where:
- $\Gamma_z$ is the axial gravitational torque component.
- $\kappa$ is the gravitational coupling constant (typically estimated between $10^{19} \text{ and } 10^{21} \text{ N}\cdot\text{m/rad}$).
- $\theta_{ic}$ and $\theta_m$ are the rotational angles of the inner core and mantle, respectively.
This restoring torque acts as a torsional spring. When the inner core advances relative to the mantle, $\Gamma_z$ decelerates the inner core and accelerates the mantle, conserving the total angular momentum of the solid Earth system. The interaction results in an oscillatory exchange of momentum with a fundamental period controlled by the coupling constant $\kappa$, the moment of inertia of the inner core ($C_{ic} \approx 2.5 \times 10^{34} \text{ kg}\cdot\text{m}^2$), and mantle rheological dissipation. This theoretical framework accounts for the observed ~60-year period and the 3–5 millisecond amplitude of multidecadal $\Delta\text{LOD}$ variations.
3. Observational Evidence and Geodetic Measurement
3.1 Historical and Modern Geodetic Techniques
Quantifying multidecadal LOD fluctuations requires continuous, multi-century observational baselines:
+----------------------------------------------------------------------------+
| Evolution of Rotational Tracking Techniques |
+-----------------------+-----------------------+----------------------------+
| Era | Primary Technique | Precision Level |
+-----------------------+-----------------------+----------------------------+
| Antiquity to ~1600 | Solar/Lunar Eclipses | ~100 to 500 ms |
| 1600 to 1950 | Star Occultations | ~10 to 50 ms |
| 1970 to Present | VLBI & SLR | < 0.01 ms (microsecond) |
| 1990 to Present | GNSS Constellations | Continuous Daily Baselines |
+-----------------------+-----------------------+----------------------------+
- Historical Eclipse Records: Solar and lunar eclipse timings documented in Babylonian, Chinese, Greek, and Arab annals dating from 700 BCE provide the baseline for secular deceleration and reveal historical long-period deviations from uniform spin.
- Telescopic Star Occultations (1620–1950s): Lunar occultations of background stars established a coherent, continuous record of multidecadal rotation anomalies spanning more than three centuries, identifying the cyclical ~60-year wave across the 18th, 19th, and 20th centuries.
- Very Long Baseline Interferometry (VLBI): By recording signals from extragalactic quasars, VLBI establishes the celestial reference frame and tracks Earth orientation parameters (EOP) with sub-millisecond precision.
- Satellite Laser Ranging (SLR) & GNSS: SLR tracks retroreflector-equipped satellites to determine low-degree gravitational field coefficients ($J_2$) and planetary rotational states. High-frequency Global Navigation Satellite System (GNSS) constellations provide continuous daily updates to Earth rotation rate models.
Multidecadal LOD Anomaly Pattern (Reconstructed Baseline)
ΔLOD (ms)
+4 | * * * *
+2 | * * * *
0 |--------*-------*---------------------*-------*-------- Nominal Baseline
-2 | * * * *
-4 | * * * *
+--------+--------+--------+--------+--------+--------+
1900 1920 1940 1960 1980 2000 2020
3.2 Corroborating Geomagnetic and Gravity Data
Multidecadal LOD cycles correlate with independent geophysical observables:
- Geomagnetic Secular Variation: Azimuthal core flows near the CMB advect magnetic field lines. Rapid changes in the second-time derivative of the geomagnetic field (geomagnetic jerks) align with phase shifts in the ~60-year $\Delta\text{LOD}$ cycle. Inversions of core surface flow patterns show that fluid outer core angular momentum inversely mirrors mantle angular momentum, verifying the internal transfer mechanism.
- Satellite Gravity Data (GRACE & GRACE-FO): The Gravity Recovery and Climate Experiment (GRACE) and its Follow-On mission measure global mass migration by tracking micrometer-scale changes in inter-satellite distance. Isolating surface hydrologic signals reveals underlying degree-2 and degree-3 gravitational variations that correspond to deep mantle-core mass redistributions.
- ESA Swarm Mission: The multi-satellite Swarm constellation isolates core-generated magnetic fields from lithospheric and ionospheric signals. High-resolution magnetic maps track core-flow oscillations and verify the localized torque balance driving mantle acceleration cycles.
4. Deep Earth Dynamics and Inner Core Oscillation
4.1 Inner Core Differential Rotation and Wobble
The solid inner core rotates semi-independently of the mantle because it is suspended within the low-viscosity liquid outer core. Seismic body waves ($PKIKP$ and $PKP$ branches) traversing the inner core display changing travel times and anisotropic waveforms over decades:
Seismic Ray Paths Probing the Inner Core:
Mantle / Surface
|
| (PKP path through outer core)
v
CMB (2890 km)
|
| (PKIKP path through solid inner core)
v
ICB (5150 km) --> [ Anisotropic Inner Core Lattice ]
These seismic shifts indicate differential rotation regimes:
- Super-Rotation: The inner core rotates faster than the mantle.
- Sub-Rotation: The inner core rotates slower than the mantle.
- Oscillatory Wobble: The inner core oscillates around an equilibrium position relative to the mantle.
Dynamic simulations and seismic doublet analyses confirm that the inner core does not continuously spin in a single direction relative to the mantle. Instead, it oscillates with an approximate 60- to 70-year periodicity.
This behavior is governed by a mechanical feedback loop:
- Convective flows in the liquid outer core exert hydromagnetic and viscous forces on the inner core.
- As the inner core departs from alignment with the mantle’s gravitational field, gravitational torque opposes the displacement.
- The gravitational torque slows inner core differential motion and reverses its relative direction, transferring angular momentum back into the mantle and modulating $\Delta\text{LOD}$.
4.2 Constraints on Mantle Viscosity and Core Density
The mechanics of gravitational coupling place bounds on lower mantle rheology and the structure of the $D’’$ layer:
- Lower Mantle Viscosity ($\eta_m$): Gravitational torque displaces non-hydrostatic anomalies, inducing viscoelastic stress in the lower mantle. The persistence of the ~60-year oscillation indicates that lower mantle viscosity must remain sufficiently high ($\eta \ge 10^{21} \text{–} 10^{22} \text{ Pa}\cdot\text{s}$) to prevent rapid viscous relaxation of gravitational potential wells. Rapid relaxation on multidecadal timescales would dissipate the gravitational spring and suppress sustained LOD oscillations.
- Inner Core Boundary (ICB) Viscoelasticity: Viscous deformation in the inner core permits localized relaxation under mantle gravitational pull. The observed phase lag between geomagnetic flow shifts and LOD inflection points provides operational bounds on inner core viscosity ($\eta_{ic} \approx 10^{16} \text{–} 10^{20} \text{ Pa}\cdot\text{s}$).
- Density Contrast across Core Interfaces: Non-hydrostatic density variations ($\Delta\rho$) along the core-mantle boundary are constrained to $\Delta\rho/\rho \approx 0.1% \text{ to } 1%$. These values align with seismic tomographic estimates for the compositional density of LLVPs.
5. Broader Geophysical and Technological Implications
5.1 Disentangling Climate Signatures from Core Dynamics
Separating deep Earth signals from surface environmental indicators is essential for long-term climate modeling:
- Atlantic Multidecadal Oscillation (AMO) & Ocean Heat Content: The Atlantic Multidecadal Oscillation and Pacific Decadal Oscillation exhibit cyclicities of 50–70 years, overlapping the core-driven ~60-year $\Delta\text{LOD}$ signal. Precise rotational decomposition prevents misidentifying internal rotational dynamics as externally forced climate variability.
- True Polar Wander (TPW) & Mean Sea-Level Adjustments: Long-term mass redistributions between core and mantle induce low-frequency polar motion (wobble of the rotational axis relative to the crust). Removing core-mantle torque signals enables accurate geoid deformation modeling, isolating ice-sheet melt and thermal expansion contributions to absolute sea-level rise.
+--------------------------------------------------------------------------------+
| Disentangling Multi-Scale Earth Dynamic Signals |
+---------------------------+------------------------+---------------------------+
| Signal Component | Primary Source | Impact on Geodesy/Climate |
+---------------------------+------------------------+---------------------------+
| 60-Year Oscillatory Cycle | Core Gravitational | True baseline for LOD |
| | Torque | filtering and TPW models |
+---------------------------+------------------------+---------------------------+
| Secular Deceleration | Lunar Tidal Drag | Long-term clock drift |
+---------------------------+------------------------+---------------------------+
| Multi-Year Trends | Ice Melt / GIA / Ocean | Polar motion, Geoid |
| | Mass Shifts | shape, Sea-level trends |
+---------------------------+------------------------+---------------------------+
5.2 Precision Timekeeping and Global Navigation Systems
Variations in LOD affect civil, astronomical, and space-based timing infrastructures:
- Coordinated Universal Time (UTC) and Leap Seconds: UTC is pegged to International Atomic Time (TAI), governed by atomic clocks. Universal Time (UT1) is based on Earth’s variable rotation. To maintain $| \text{UTC} - \text{UT1} | < 0.9 \text{ seconds}$, leap seconds are periodically inserted or omitted. Multidecadal acceleration phases reduce the frequency of positive leap seconds and can necessitate negative leap seconds.
- Global Navigation Satellite Systems (GNSS): Systems including GPS, GLONASS, Galileo, and BeiDou depend on sub-nanosecond signal synchronization. Precise orbital propagation and coordinate frame transformations (from the International Terrestrial Reference Frame to the International Celestial Reference Frame) require continuous compensation for LOD fluctuations driven by core-mantle torque.
- Deep Space Telemetry: Interplanetary navigation requires sub-milliradian antenna pointing precision. Unmodeled rotational shifts in Earth’s crust produce geometric projection errors over astronomical baselines.
6. Future Directions in Geodynamic Modeling
6.1 Advanced Magnetohydrodynamic (MHD) Simulations
Computational geodynamics continues to refine the resolution of coupled boundary processes:
- Resolving Low Ekman Numbers: Earth’s outer core exhibits an Ekman number ($E = \nu / 2\Omega L^2$) on the order of $10^{-15}$, reflecting low fluid viscosity relative to planetary rotation. Supercomputing simulations have reached $E \sim 10^{-6} \text{–} 10^{-7}$, narrowing the gap toward realistic turbulence modeling.
- Coupled Mantle-Core Systems: Modern models integrate 3D spherical mantle convection grids with convective dynamos, using real seismic tomography data directly at the CMB boundary rather than idealized spherical surfaces.
Coupled Geodynamic Model Architecture
+-------------------------------------------------------+
| Mantle Rheology Module (Viscoelasticity, LLVP Tomography) |
+---------------------------+---------------------------+
| Boundary Flux & Grav Torque
+---------------------------v---------------------------+
| High-Resolution MHD Dynamo Module (Low Ekman Number) |
+---------------------------+---------------------------+
| Mechanical / Field Coupling
+---------------------------v---------------------------+
| Inner Core Kinematics Module (Anisotropy & Deformation)|
+-------------------------------------------------------+
6.2 Next-Generation Gravitational and Seismic Probing
- Dense Seismic Arrays and Reprocessed Multiplets: Deployments of dense, broadband seismic networks (e.g., USArray, China Array) yield refined datasets of core-phase waveforms. Standardizing the processing of repeating doublet earthquakes provides clear views of ICB boundary displacement and inner core rotational drift.
- Quantum Gravimetry and Next-Gen Satellite Missions: Emerging quantum cold-atom gravimeters and planned successor missions to GRACE-FO aim to measure higher-degree gravitational harmonics with minimal instrument drift. These measurements will isolate localized core-mantle mass anomalies, providing direct real-time tracking of internal gravitational torque.
Frequently Asked Questions (FAQ)
What causes multidecadal variations in the length of day?
Multidecadal LOD variations are driven by angular momentum exchange between Earth’s solid mantle and core. Gravitational torque between non-axisymmetric density anomalies in the mantle (such as Large Low Velocity Provinces) and the solid inner core acts as an internal restoring mechanism, generating ~60-year periodic cycles in planetary spin.
How does gravitational torque differ from tidal friction?
Tidal friction occurs through external gravitational interactions with the Moon and Sun, converting rotational kinetic energy into heat through oceanic and solid-Earth tides to permanently slow Earth’s spin over geological timescales. Internal gravitational torque occurs within the planet, transferring momentum between internal layers in an oscillatory manner without dissipating the planet’s total angular momentum.
How do scientists measure changes in the length of day?
Rotational measurements rely on space geodetic methods: Very Long Baseline Interferometry (VLBI) measuring extragalactic quasars, Satellite Laser Ranging (SLR) tracking geodetic satellites, and continuous tracking through Global Navigation Satellite Systems (GNSS). Historical variations are constrained by star occultation records and ancient solar/lunar eclipse observations.
Does the changing length of day affect daily life or technology?
Because multidecadal variations modify the length of day by only a few milliseconds, they are imperceptible to human senses. Modern technologies depend on precise synchronization: satellite navigation systems (GPS), satellite orbit tracking, astronomical observations, and telecommunication protocols require continuous LOD tracking to align Earth-fixed frames with celestial coordinates.
Can climate change alter Earth’s rotation in the same way as core gravitational torque?
Climate change influences Earth rotation by transferring mass from melting glaciers and ice sheets into the oceans, altering Earth’s moment of inertia. This process induces secular and decadal adjustments in Earth rotation and polar motion, which differ in physical mechanism, spatial footprint, and temporal behavior from the ~60-year cyclical signals produced by deep Earth gravitational torque.
Artikel Terkait
-
Inside Reentry: Real-Time Heat Shield Dynamics
-
First Direct Radio Signal Detected from an Exoplanet
-
Was the Body Conscious Before Language Evolved?
-
CERN Disconnects the LHC: Roadmap and HL-LHC Upgrades
-
How Dogs Process Human Speech: Neuroscience Insights
-
Flexible Boron Breakthrough: 10M-Fold Conductivity Jump