1.3 Evidence for Plate Tectonics
Paleomagnetism & Apparent Polar Wander
Paleomagnetism — the study of the ancient geomagnetic field recorded in rocks — provided the first quantitative, physics-based evidence for continental motion. When igneous rocks cool through the Curie temperature or sedimentary grains settle in a magnetic field, they acquire a natural remanent magnetization (NRM) that records the direction of the ambient field at the time of formation.
The direction of the paleomagnetic field is described by two angles: the declination $D$ (the horizontal angle from geographic north) and the inclination $I$ (the dip below horizontal). For a geocentric axial dipole (GAD) field — the time-averaged geomagnetic field — the inclination is related to the geographic (paleo)latitude by the fundamental equation:
The dipole equation (magnetic inclination–latitude relation):
\[ \tan(I) = 2\tan(\lambda) \]
where $I$ is the magnetic inclination and $\lambda$ is the geographic latitude. At the equator ($\lambda = 0°$), $I = 0°$ (horizontal field); at the poles ($\lambda = 90°$), $I = 90°$ (vertical field).
By measuring the paleomagnetic direction in rocks of known age and solving the dipole equation, one can calculate the position of the paleomagnetic pole (the "virtual geomagnetic pole" or VGP) at the time the rock formed. When poles of successive ages from a single continent are plotted, they define a path called the Apparent Polar Wander Path (APWP).
Paleomagnetic Pole Calculation
Given a site at coordinates $(\lambda_s, \phi_s)$ with measured declination $D$ and inclination $I$, the paleomagnetic colatitude $p$ and pole position $(\lambda_p, \phi_p)$ are:
\[ p = \cot^{-1}\!\left(\frac{\tan I}{2}\right) \]
\[ \lambda_p = \sin^{-1}\!\left(\sin\lambda_s \cos p + \cos\lambda_s \sin p \cos D\right) \]
\[ \phi_p = \phi_s + \sin^{-1}\!\left(\frac{\sin p \sin D}{\cos \lambda_p}\right) \]
where $p$ is the magnetic colatitude (angular distance from site to pole), and $(\lambda_p, \phi_p)$ is the computed paleomagnetic pole position.
The critical observation by Runcorn, Creer, and Irving in the 1950s was that APWPs constructed from different continents were discordant — they did not coincide. Since there can be only one geomagnetic pole at any given time (the GAD hypothesis), the discrepancy must mean that the continents themselves have moved relative to each other. When the continents are reassembled in their Pangaean configuration, the APWPs converge onto a single path, providing powerful quantitative support for continental drift.
GPS Geodesy & Modern Plate Velocity Measurements
Since the 1980s, space-geodetic techniques have allowed direct measurement of plate motions in real time, providing the most compelling modern evidence for plate tectonics. Three principal techniques are used:
GPS / GNSS
The Global Positioning System provides position measurements with millimeter-level precision. Continuous GPS stations (e.g., IGS network) track station velocities over years to decades, resolving plate motions of 1–100+ mm/yr.
VLBI
Very Long Baseline Interferometry measures the time delay of quasar signals arriving at radio telescopes on different continents, yielding intercontinental baseline changes with sub-millimeter precision.
SLR
Satellite Laser Ranging measures the round-trip travel time of laser pulses to retroreflector-equipped satellites (e.g., LAGEOS), providing station positions and plate motions with centimeter-level accuracy.
All geodetic measurements are expressed in the International Terrestrial Reference Frame (ITRF), currently ITRF2020. The ITRF is a no-net-rotation frame that defines positions and velocities of hundreds of geodetic stations worldwide. GPS velocities from the ITRF confirm the plate velocity predictions of geological models (NUVEL-1A, MORVEL) to within 1–2 mm/yr for most plate pairs, demonstrating that plate motions have been approximately steady over the past few million years.
Notable GPS Results
North America–Europe separation: ~23 mm/yr (consistent with Mid-Atlantic Ridge spreading)
Pacific–North America (San Andreas): ~46 mm/yr right-lateral
India–Eurasia convergence: ~35–40 mm/yr (driving Himalayan collision)
Australia northward motion: ~65–70 mm/yr (fastest continental plate)
Global Earthquake Distribution
The global distribution of earthquakes is one of the most visually compelling lines of evidence for plate tectonics. Approximately 95% of the world's seismic energy is released along narrow belts that define the boundaries between tectonic plates. The remaining ~5% occurs within plate interiors (intraplate seismicity), often associated with ancient zones of weakness or stress concentrations.
~95%
Seismic Energy at Plate Boundaries
670 km
Maximum Earthquake Depth
~80%
Energy in Subduction Zones
The depth distribution of earthquakes varies systematically with tectonic setting:
Divergent Boundaries (Mid-Ocean Ridges)
Shallow earthquakes only (0–30 km depth), confined to the thin, hot lithosphere at the ridge axis. Magnitudes rarely exceed M 6. Focal mechanisms show normal faulting (extension) perpendicular to the ridge strike, with strike-slip mechanisms on transform offsets.
Transform Boundaries
Shallow to intermediate depth earthquakes (0–30 km), generally restricted to the active transform segment between ridge offsets. The San Andreas Fault system, a continental transform, generates earthquakes to ~20 km depth (the seismogenic zone) with predominantly strike-slip focal mechanisms.
Convergent Boundaries (Subduction Zones)
Earthquakes span the full range from shallow (0–70 km) through intermediate (70–300 km) to deep focus (300–670 km). The inclined zone of seismicity — the Wadati–Benioff zone — delineates the geometry of the descending slab. The deepest earthquakes occur at ~670 km, near the 660-km phase transition (ringwoodite to bridgmanite + ferropericlase), where slabs may stall or penetrate into the lower mantle.
Volcanic Distribution & the Ring of Fire
Like earthquakes, volcanic activity is concentrated along plate boundaries, particularly at convergent margins. The Ring of Fire — a horseshoe-shaped belt of volcanoes and seismic activity encircling the Pacific Basin — contains approximately 75% of the world's active and dormant volcanoes and accounts for about 90% of the world's earthquakes.
Subduction Zone Volcanism
Volcanic arcs form 100–200 km above the descending slab, where released fluids (primarily H2O from dehydration of serpentine, chlorite, and amphibole) lower the solidus of the overlying mantle wedge, triggering partial melting. The resulting magmas are characteristically calc-alkaline, enriched in large-ion lithophile elements (LILE) and depleted in high-field-strength elements (HFSE).
Mid-Ocean Ridge Volcanism
The mid-ocean ridge system is the world's most volcanically active feature, producing ~21 km³ of new oceanic crust per year through decompression melting of upwelling asthenosphere. Ridge basalts (MORB) are tholeiitic, depleted in incompatible elements, and remarkably uniform in composition worldwide, reflecting a well-mixed depleted mantle source.
Intraplate volcanism (hotspots such as Hawaii, Iceland, and Réunion) is not directly associated with plate boundaries but is attributed to mantle plumes — buoyant upwellings from the deep mantle that produce age-progressive volcanic chains as plates move over relatively stationary plume sources.
Gravity Anomalies & Satellite Geodesy
Gravity measurements provide information about the subsurface density distribution and the state of isostatic compensation. At convergent plate boundaries, deep-sea trenches are associated with large negative Bouguer gravity anomalies (as negative as −300 mGal), indicating a mass deficit beneath the trench. This deficit arises from the depression of the ocean floor, the low-density sedimentary fill, and the dynamic suction of the subducting slab.
Bouguer anomaly (simplified, onshore):
\[ \Delta g_B = g_{\text{obs}} - g_{\text{ref}} + \delta g_{\text{FA}} - \delta g_{\text{Boug}} + \delta g_{\text{terr}} \]
where $\delta g_{\text{FA}} = 0.3086 \, h$ mGal is the free-air correction (h in meters),$\delta g_{\text{Boug}} = 2\pi G \rho h = 0.04193 \, \rho \, h$ mGal is the Bouguer slab correction, and $\delta g_{\text{terr}}$ is the terrain correction.
Modern satellite missions have revolutionized our ability to map Earth's gravity field:
GRACE / GRACE-FO
The Gravity Recovery and Climate Experiment (2002–2017) and its follow-on (2018–present) use twin satellites to measure time-variable gravity. GRACE detects mass redistribution from ice sheet melting, groundwater changes, and post-glacial rebound with microgal sensitivity, resolving features at ~300 km spatial resolution.
InSAR
Interferometric Synthetic Aperture Radar measures ground deformation with millimeter precision by comparing the phase of radar signals from repeat satellite passes. InSAR has mapped co-seismic and inter-seismic deformation at plate boundaries worldwide, providing detailed strain maps that constrain fault locking and slip rates.
Gravity Signatures of Plate Boundaries
The gravity field across convergent margins shows a distinctive pattern:
Trench (Negative Anomaly)
A pronounced free-air gravity low (−100 to −200 mGal) over the trench axis, reflecting the deep bathymetric depression and the dynamic deflection of the plate. The Bouguer anomaly is even more negative (−200 to −350 mGal) because the water column represents a large mass deficit relative to rock.
Outer Rise (Slight Positive)
A modest free-air gravity high (+20 to +50 mGal) on the oceanic plate seaward of the trench, corresponding to the flexural bulge where the plate bends before subducting. The wavelength of this bulge (~150–250 km) constrains the effective elastic thickness of the lithosphere.
Volcanic Arc (Positive Anomaly)
A positive free-air anomaly (+50 to +100 mGal) over the volcanic arc, reflecting the topographic mass of the arc and the dense roots of the arc crust. The island arc or continental margin is approximately isostatically compensated.
Vening Meinesz's Contribution
Felix Vening Meinesz conducted pioneering submarine gravity measurements in the 1920s–1930s across the Indonesian arc, discovering the enormous negative gravity anomalies associated with ocean trenches. He correctly interpreted these as zones of crustal downbuckling, anticipating the concept of subduction by several decades.
Summary of Evidence
| Evidence Type | Key Observation | Tectonic Implication |
|---|---|---|
| Paleomagnetism | Discordant APWPs from different continents | Continents have moved relative to each other |
| GPS Geodesy | Station velocities of 1–100+ mm/yr | Plates are moving now at measurable rates |
| Seismicity | Earthquakes concentrated in narrow belts | Active deformation at plate boundaries |
| Volcanism | Ring of Fire, mid-ocean ridge system | Magma generation linked to plate dynamics |
| Gravity | Negative anomalies at trenches | Dynamic processes at convergent margins |
| Satellite Geodesy | InSAR deformation maps, GRACE mass changes | Ongoing strain accumulation and release |