Virtual Geomagnetic Pole Calculator: From Paleomagnetic Directions to Pole Positions

simulator intermediate ~10 min
Loading simulation...
VGP = 85°N, 170°E — near the geographic pole

A site at 45°N, 10°E measuring D = 5°, I = 63° yields a VGP at 85°N, 170°E — consistent with a geocentric axial dipole field and confirming stable normal polarity magnetization.

Formula

p = arctan(2 / tan I)
φp = arcsin(sin φs cos p + cos φs sin p cos D)
tan I = 2 × tan λ (dipole formula)

From Compass Needles to Pole Positions

When a paleomagnetist measures the declination and inclination of remanent magnetization in a rock sample, the raw data is a direction in local coordinates. The Virtual Geomagnetic Pole (VGP) transformation converts this local direction into a global pole position — the point on the globe where a geocentric dipole would need to sit to produce the observed field. This transformation is the essential bridge between laboratory measurements and plate tectonic reconstructions.

The Dipole Formula

The calculation relies on the geocentric dipole field equations: tan I = 2 tan λ, where I is inclination and λ is magnetic latitude. This relationship gives the magnetic colatitude p = arctan(2/tan I), the angular distance from the site to the pole along the great circle defined by the declination. Spherical trigonometry then converts (site + colatitude + declination) into the VGP coordinates on the globe.

Testing the GAD Hypothesis

If Earth's time-averaged field truly behaves as a geocentric axial dipole, then VGPs from rocks of the same age on the same continent should cluster tightly around the geographic pole (for recent rocks) or around a consistent paleomagnetic pole (for ancient rocks). The remarkable success of this prediction — VGPs from globally distributed recent lavas cluster within ~5° of the geographic pole — validates the GAD hypothesis and underpins the entire paleomagnetic method.

VGPs During Reversals

During polarity reversals, the orderly dipole field breaks down and VGPs scatter wildly across the globe. Transitional VGP paths trace the reversal process — how the field unravels, reorganizes, and re-establishes in opposite polarity. Studies of these paths have revealed possible preferred longitudes (through the Americas and East Asia) that may reflect control by lower mantle heterogeneity on core flow, though this remains one of the most debated topics in paleomagnetism.

FAQ

What is a virtual geomagnetic pole?

A Virtual Geomagnetic Pole (VGP) is the position where a geocentric dipole would need to be located to produce the observed magnetic field direction (declination and inclination) at a given site. It is the fundamental data unit of paleomagnetism, converting local field measurements into global pole positions that can be compared across continents.

How is VGP calculated from declination and inclination?

First, compute the magnetic colatitude p = arctan(2/tan I) using the dipole formula. Then the pole latitude is φp = arcsin(sin φs cos p + cos φs sin p cos D), and pole longitude is λp = λs + arcsin(sin p sin D / cos φp), where φs, λs are the site coordinates and D, I are the measured declination and inclination.

What is the geocentric axial dipole hypothesis?

The GAD hypothesis states that when averaged over sufficient time (> 10,000 years), Earth's magnetic field is equivalent to a dipole at the center of the Earth aligned with the rotation axis. This is the foundational assumption of paleomagnetism — it allows paleomagnetic directions to be converted into paleolatitudes, enabling reconstruction of past continental positions.

Why do VGPs scatter during reversals?

During a polarity reversal, the axial dipole weakens and the field is dominated by non-dipole components. The VGPs calculated from these transitional directions scatter widely across the globe, creating VGP paths that trace the reversal. Preferential VGP paths through the Americas and East Asia have been reported but remain controversial.

Sources

Embed

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