Paleointensity Simulator: Measuring the Strength of Earth's Ancient Magnetic Field

simulator intermediate ~10 min
Loading simulation...
Banc = 40 μT — ancient field comparable to present day

An NRM/TRM slope of -0.8 with Blab = 50 μT yields Banc = 40 μT, corresponding to a virtual dipole moment of ~8×10²² A·m² — close to the present-day value of ~30 μT at the equator.

Formula

Banc = |slope| × Blab (Thellier equation)
VDM = (4π R³ / μ₀) × Banc / √(1 + 3sin²λ)
slope = -NRM_lost / TRM_gained

Measuring Ancient Field Strength

While paleomagnetic directions are recorded directly by the orientation of magnetic minerals, field intensity is encoded in the magnitude of magnetization — a much more difficult quantity to extract reliably. The Thellier double-heating method, developed in 1959, remains the gold standard: by systematically replacing the rock's natural magnetization with a laboratory-induced one at progressively higher temperatures, the ratio of ancient to laboratory field strength can be determined with high precision.

The Arai Plot

The Thellier experiment produces data plotted on an Arai diagram: NRM remaining versus TRM acquired at each temperature step. An ideal specimen yields a perfectly straight line whose slope directly gives the paleointensity ratio. In practice, thermal alteration (mineral changes during heating) causes the line to curve at high temperatures, and pTRM checks — repeating earlier steps to detect alteration — are essential quality controls. Only specimens passing strict linearity and alteration criteria yield reliable results.

Virtual Dipole Moments

To compare paleointensity measurements from different latitudes, the local field strength is converted to a Virtual Dipole Moment (VDM) using the dipole formula. This removes the latitude dependence (the field is twice as strong at the poles as at the equator) and provides a global measure of geodynamo strength. The present-day dipole moment is about 8×10²² A·m²; values during the Mesozoic Dipole Low dropped below 4×10²² A·m².

The Geodynamo Through Time

Paleointensity records spanning billions of years reveal a dynamic geodynamo. The field was relatively weak and volatile before the inner core nucleated (estimated at 1-1.5 Ga), strengthened dramatically afterward, experienced the prolonged Mesozoic Dipole Low during a period of infrequent reversals, and has fluctuated on million-year timescales ever since. Understanding these long-term intensity variations constrains models of core evolution, inner core growth, and the thermal history of the planet.

FAQ

What is paleointensity?

Paleointensity is the strength (intensity) of Earth's magnetic field at some time in the geological past, determined from the magnetization of rocks. Unlike paleomagnetic directions which are relatively straightforward to measure, paleointensity determination is technically challenging and requires specialized experiments like the Thellier double-heating method.

How does the Thellier method work?

The Thellier method progressively replaces the rock's natural remanent magnetization (NRM) with a laboratory thermoremanent magnetization (TRM) through stepwise heating in a known field. Plotting NRM lost vs TRM gained at each step creates an Arai plot — if the relationship is linear, the slope gives the ratio of ancient to laboratory field: Banc = |slope| × Blab.

What is an Arai plot?

An Arai plot graphs NRM remaining (y-axis) versus TRM gained (x-axis) at each temperature step in a Thellier experiment. An ideal sample produces a straight line whose slope equals -Banc/Blab. Curvature indicates thermal alteration (chemical change during heating) or multi-domain behavior, which invalidate the paleointensity estimate.

How has Earth's field strength changed over time?

The virtual axial dipole moment (VADM) has fluctuated between ~2×10²² and ~12×10²² A·m² over the past 3.5 billion years. Notable features include the Mesozoic Dipole Low (~80-120 Ma), the field strengthening associated with inner core nucleation (~1-1.5 Ga), and short-term variations linked to polarity reversals and excursions.

Sources

Embed

<iframe src="https://homo-deus.com/lab/paleomagnetism/paleointensity/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub