10.3 Isostasy

Archimedes' Principle for the Earth

Isostasy is the condition of gravitational equilibrium in which the Earth's lithosphere “floats” on the denser, fluid-like asthenosphere, much as a block of wood floats in water. The concept follows directly from Archimedes' principle: a floating body displaces a mass of fluid equal to its own mass. Mountains, with their thick, low-density crustal roots, ride higher because they displace more mantle, while ocean basins, underlain by thin, dense oceanic crust, sit lower.

The idea emerged in the 19th century from geodetic surveys in India. George Everest (1847) and John Henry Pratt noticed that the gravitational deflection of a plumb line near the Himalaya was much less than predicted from the visible mass of the mountains alone. George Biddell Airy (1855) and Pratt (1855) independently proposed models to explain this deficit — both invoking subsurface mass compensation, but through different mechanisms. Together with the later flexural isostasy model of Vening Meinesz (1931), these form the three classical models of isostatic compensation.

Airy Isostasy: Variable Thickness, Constant Density

In the Airy model, isostatic compensation is achieved by varying the thickness of a constant-density crust. Mountains have deep crustal roots that extend downward into the mantle, displacing denser mantle material. The higher the mountain, the deeper the root. The crust is treated as a series of columns of uniform density $\rho_c$ floating on a denser mantle of density $\rho_m$.

The Airy model successfully explains the first-order observation from seismology that mountain belts have thick crust and ocean basins have thin crust. Receiver function studies beneath the Himalaya show crustal thicknesses of 65–75 km, consistent with a root supporting the 5 km mean elevation of Tibet. Similarly, the Andes show crustal thicknesses up to 70 km beneath the Altiplano.

Pratt Isostasy: Variable Density, Constant Depth

In the Pratt model, isostatic compensation is achieved by lateral variations in the density of crustal columns that all extend to a common compensation depth $D$. Elevated regions have lower-density crust, while depressed regions (like ocean basins) have higher-density crust. All columns have the same total mass per unit area above the compensation depth.

The Pratt model is most applicable to oceanic lithosphere, where lateral density variations arise naturally from differences in thermal state. Young oceanic lithosphere at mid-ocean ridges is hot and therefore less dense, standing higher (the ridge itself). As it cools and moves away from the ridge, it densifies and subsides — precisely the pattern predicted by the Pratt model. The depth–age relationship of the ocean floor ($d \propto \sqrt{t}$) is essentially a Pratt isostatic response to lateral density variations caused by conductive cooling.

Vening Meinesz (Flexural) Isostasy

The Vening Meinesz model (also called flexural isostasy) recognizes that the lithosphere has finite elastic strength and distributes applied loads over an area wider than the load itself, rather than compensating each column independently as in the Airy and Pratt models. A point load on an elastic plate produces a broad depression surrounded by a flexural forebulge, as described by the flexure equation in section 10.2.

The isostatic response depends on the wavelength of the load relative to the flexural wavelength of the plate. This is captured by the isostatic response function in the spectral domain:

Flexural isostasy is the most physically realistic of the three models. It naturally transitions between rigid support of small loads (like individual volcanoes) and full Airy compensation of very large loads (like entire mountain belts). The single free parameter — the effective elastic thickness $T_e$ — controls the transition wavelength, making flexural isostasy a powerful tool for probing the mechanical properties of the lithosphere.

Gravity Anomalies & Isostatic Compensation

Gravity anomalies provide the primary observational test of isostatic models. Different corrections to the observed gravity field isolate different aspects of mass distribution:

Isostatic Rebound & Mantle Viscosity

When a load is removed from the lithosphere (e.g., by deglaciation), the depressed surface rebounds toward its equilibrium position. This post-glacial rebound (or glacial isostatic adjustment, GIA) is one of the most important geophysical observables for constraining mantle viscosity, because the timescale of rebound depends directly on how fast the mantle flows to fill the space vacated by the rebounding lithosphere.

Observations of post-glacial rebound have yielded among the best estimates of mantle viscosity:

Compensation Time & Dynamic Equilibrium

The compensation time is the time required for the lithosphere to reach isostatic equilibrium after a loading or unloading event. It depends on the mantle viscosity and the spatial scale of the load. For the upper mantle viscosity of$\sim 10^{21}$ Pa·s and load wavelengths of ~1,000 km, the compensation time is on the order of 10,000–30,000 years.

This means that loads changing faster than the compensation time produce transient isostatic disequilibrium. During glacial cycles (period ~100,000 years), the ice sheets grow and retreat faster than the mantle can respond, creating persistent departures from equilibrium. The current ongoing rebound of Fennoscandia and Hudson Bay is a direct consequence of this lag.

For orogenic processes, which act over millions to tens of millions of years, the lithosphere is generally close to isostatic equilibrium. The crust beneath the Himalaya and Tibet is approximately in Airy isostatic balance, as confirmed by the near-zero free-air gravity anomaly averaged over the orogen. However, at shorter wavelengths and near active thrusts, departures from equilibrium are common, reflecting the finite elastic strength of the plate (flexural isostasy) and the time lag in the isostatic response to rapidly advancing thrust loads.