3.3 Isostasy and its applications

Isostasy explains how Earth's crust floats on the denser mantle, balancing weight and buoyancy. This concept is central to understanding gravity anomalies, crustal thickness variations, and why the Earth's surface moves vertically over time.

Geophysicists use several models to describe isostasy, including the Airy, Pratt, and flexural approaches. These models help interpret gravity data, estimate crustal properties, and explain dynamic processes like mountain building and basin formation.

Isostasy: Concept and Principles

Isostatic Equilibrium and Compensation

Isostasy is the state of gravitational equilibrium between Earth's crust and mantle, where the lighter crust "floats" on the denser mantle much like an iceberg floats on water. Thicker or less dense crust rides higher, while thinner or denser crust sits lower.

Isostatic equilibrium is reached when the weight of a crustal column is exactly balanced by the buoyancy force from the mantle, so there's no net vertical force. The depth of compensation is the depth below which pressure from overlying columns becomes equal everywhere. This depth roughly corresponds to the base of the lithosphere, though it's often discussed in relation to the Moho (crust-mantle boundary), which varies with crustal thickness and density.

Isostatic Adjustments and Disturbances

When something disturbs the equilibrium, the crust slowly rises or sinks to restore balance. These are isostatic adjustments, and they're driven by processes that add or remove mass at the surface:

• Glaciation/deglaciation: Ice sheets load the crust downward; when they melt, the crust rebounds. Scandinavia is still rising at up to ~10 mm/yr after the last ice age, and Hudson Bay in Canada shows similar post-glacial rebound.
• Erosion and deposition: Removing material from mountains causes uplift, while depositing sediment in basins causes subsidence (e.g., the Gulf Coast and North Sea basins).
• Tectonic loading: Mountain building (Himalayas, Andes) or emplacement of dense igneous intrusions (Bushveld Complex, South Africa) can push the crust out of equilibrium.

The timescale of adjustment depends on mantle viscosity and the spatial scale of the disturbance. Post-glacial rebound operates over thousands to tens of thousands of years, while responses to large-scale tectonic loading can take millions of years.

Models of Isostasy

Airy and Pratt Isostatic Models

The Airy model assumes the crust has a uniform density but varies in thickness. Topographic highs are compensated by deep crustal roots extending into the mantle, while topographic lows (like ocean basins) have thinner crust. Think of wooden blocks of different heights floating in water: taller blocks stick up higher but also extend deeper below the waterline.

• The compensation condition can be written as: $\rho_c \cdot t = \rho_m \cdot r$, where $\rho_c$ is crustal density, $t$ is the elevation above a reference, $\rho_m$ is mantle density, and $r$ is the root depth below the reference Moho.
• This model works well in regions with large crustal roots, such as the Himalayas and the Andes, where seismic data confirm thickened crust beneath high topography.

The Pratt model assumes the crust has a uniform thickness but varies in density. Mountains are underlain by less dense crust, while lowlands and ocean floors are underlain by denser crust. The compensation depth (base of all columns) is the same everywhere.

• Mathematically: $\rho_1 \cdot h_1 = \rho_2 \cdot h_2 = \text{constant}$ for all columns extending to the compensation depth.
• This model applies well in regions with lateral density variations, such as the Basin and Range Province (where thermal expansion reduces crustal density beneath elevated terrain) and the East African Rift System.

Flexural Isostatic Model and Combined Approaches

The Vening Meinesz (flexural) model treats the lithosphere as an elastic plate that bends under applied loads rather than responding purely locally. A point load doesn't just push the crust down directly beneath it; the rigidity of the plate spreads the deflection over a wider area, creating a depression under the load flanked by a subtle peripheral bulge.

• The key parameter is flexural rigidity ($D$), which depends on the elastic thickness ($T_e$) of the lithosphere: $D = \frac{E T_e^3}{12(1 - \nu^2)}$, where $E$ is Young's modulus and $\nu$ is Poisson's ratio.
• This model explains features like the moat and arch around the Hawaiian Islands (where the volcanic load bends the Pacific plate) and the broad rebound pattern of the Fennoscandian Shield.

In practice, Earth's isostatic behavior combines elements of all three models. The relative importance of each depends on the geological setting and spatial scale. At very large scales (hundreds of km), local Airy or Pratt compensation dominates. At intermediate scales, flexural effects become significant. Combined approaches like the Airy-Heiskanen and Pratt-Hayford models incorporate elements of both end-member models for more realistic representations. Choosing the right model depends on the tectonic setting, crustal age and composition, and available geophysical constraints.

Isostasy and Gravity Anomalies

Gravity Anomalies and Their Relationship to Isostasy

A gravity anomaly is the difference between observed gravity at a point and the theoretical gravity predicted by a reference model (typically the reference ellipsoid). Because isostatic compensation redistributes mass at depth, it directly shapes the gravity field measured at the surface.

• A positive gravity anomaly indicates excess mass beneath the observation point, which could mean dense rocks near the surface or topography that isn't fully compensated by a crustal root.
• A negative gravity anomaly indicates a mass deficit, such as a thick low-density crustal root or overcompensated topography.

If a mountain range is perfectly isostatically compensated, the extra mass of the topography above is offset by the low-density root below, and the gravity anomaly (after appropriate corrections) should be near zero. Departures from zero tell you the region is out of isostatic equilibrium.

Types of Gravity Anomalies and Their Calculation

Each type of gravity anomaly removes different effects, progressively isolating subsurface structure:

1. Free-air anomaly: Start with observed gravity, then apply the free-air correction to account for the station's elevation above the reference ellipsoid. This anomaly reflects the total mass distribution beneath the observation point but doesn't remove the gravitational effect of topographic mass itself.

2. Bouguer anomaly: Take the free-air anomaly and apply the Bouguer correction, which removes the gravitational attraction of the rock mass between the station and the reference level (plus terrain corrections for irregular topography). The result highlights subsurface density variations. Over continents, Bouguer anomalies are typically negative because the correction removes topographic mass but the compensating root remains as a mass deficit.

3. Isostatic anomaly: Take the Bouguer anomaly and remove the calculated gravitational effect of the assumed isostatic compensation (using an Airy, Pratt, or flexural model). What's left reveals departures from isostatic equilibrium. A near-zero isostatic anomaly means the chosen model fits well.

Notable examples: the Tibetan Plateau shows a strongly positive free-air anomaly (high elevation with large total mass) but a strongly negative Bouguer anomaly (thick, low-density crustal root). The Mid-Atlantic Ridge shows a negative Bouguer anomaly due to hot, low-density mantle material beneath the ridge.

Applying Isostatic Principles

Crustal Thickness and Density Estimation

By assuming an isostatic model and using observed topography and gravity data, you can estimate crustal thickness and density variations across a region. For example, under the Airy model, you can invert topographic data to predict Moho depth. These gravity-derived estimates are then validated against independent seismic data (from seismic refraction surveys or receiver function analysis). Where the two agree, confidence in the crustal model is high. Where they disagree, it points to regions that may be out of isostatic equilibrium or where the chosen model is too simple.

Vertical Motions and Sedimentary Basin Evolution

Isostatic models predict how the crust responds vertically to changes in surface loading:

• Post-glacial rebound: After ice sheets melt, the unloaded crust rises. GPS measurements in Scandinavia, Canada, and Antarctica track this ongoing uplift, and the rate of rebound provides constraints on mantle viscosity (typically $10^{20}$ to $10^{21}$ Pa·s for the upper mantle).
• Sedimentary isostasy: As sediments accumulate in a basin, their weight causes further subsidence, which creates more accommodation space for additional sediment. This positive feedback loop explains why basins like the Gulf Coast and North Sea can accumulate sedimentary sequences many kilometers thick. The relationship is roughly that for every 1 km of sediment added, the basin floor subsides by an additional amount depending on the density contrast: $w \approx h \cdot \frac{\rho_s}{\rho_m - \rho_s}$, where $\rho_s$ is sediment density and $\rho_m$ is mantle density.

Mountain Building and Geodynamic Modeling

During orogenesis (mountain building), crustal thickening from collision or compression triggers isostatic uplift. The Tibetan Plateau (~5 km average elevation) is supported by a crustal root extending to ~70 km depth, roughly double the global continental average of ~35 km. The Altiplano-Puna Plateau in the Central Andes shows a similar pattern.

Isostatic compensation is built into numerical geodynamic models that simulate long-term lithospheric deformation. Thin-sheet models and finite element models incorporate isostatic response to tectonic loading, allowing researchers to test hypotheses about the forces driving crustal evolution. Without isostasy, these models would predict unrealistic topography and gravity fields.