10.2 Lithospheric Flexure
The Bending Lithosphere
The lithosphere possesses long-term mechanical strength and responds to applied loads by bending, much like an elastic beam on a fluid foundation. When a mountain belt grows, an ice sheet advances, or a volcanic island chain forms, the lithosphere flexes downward beneath the load and upward in adjacent regions, creating characteristic patterns of subsidence and uplift that profoundly shape sedimentary basins, coastlines, and topography.
The theory of lithospheric flexure, rooted in the mechanics of thin elastic plates, is one of the most successful quantitative frameworks in tectonics. It explains features ranging from the deep trenches flanking subduction zones to the broad foreland basins ahead of mountain belts, from the moats surrounding oceanic islands to the patterns of post-glacial rebound.
The Flexure Equation
The lithosphere is modeled as a thin elastic plate floating on a denser fluid substratum (the asthenosphere). When a vertical load $q(x)$ is applied to the plate, it deflects by an amount $w(x)$. The restoring force comes from the buoyancy of the displaced asthenospheric fluid. The governing equation in one dimension is:
General flexure equation (1D):
\[ D \frac{d^4 w}{dx^4} + \Delta\rho \, g \, w = q(x) \]
where $D$ is the flexural rigidity of the plate, $w(x)$ is the vertical deflection, $\Delta\rho = \rho_m - \rho_{\text{infill}}$ is the density contrast between the mantle and the material filling the deflection (water, sediment, or air), and $q(x)$ is the applied load per unit area.
This is a fourth-order ordinary differential equation, reflecting the fact that the plate's resistance to bending involves the fourth derivative of the deflection. The equation is the geodynamic analog of the Euler–Bernoulli beam equation and was first applied to geological problems by Vening Meinesz (1931) and Gunn (1943).
Flexural Rigidity & Elastic Thickness
The flexural rigidity $D$ quantifies the plate's resistance to bending. It depends on the elastic properties of the lithosphere and, critically, on the cube of the effective elastic thickness $T_e$:
Flexural rigidity:
\[ D = \frac{E \, T_e^3}{12(1 - \nu^2)} \]
where $E \approx 70\text{--}100$ GPa is Young's modulus,$\nu \approx 0.25$ is Poisson's ratio, and $T_e$ is the effective elastic thickness. The cubic dependence on $T_e$ means that even modest changes in elastic thickness produce dramatic changes in flexural behavior.
The effective elastic thickness $T_e$ is not the total thickness of the lithosphere but rather the thickness of the layer that behaves elastically on geological timescales. It is always less than the total mechanical thickness because the uppermost crust deforms brittlely and the lowermost lithosphere creeps ductilely:
Young Oceanic Lithosphere (<25 Ma)
$T_e \approx 5\text{--}10$ km. The lithosphere is thin and warm, with a shallow brittle-ductile transition. The plate is weak and flexes easily, producing small-wavelength deformation features.
Old Oceanic Lithosphere (>80 Ma)
$T_e \approx 20\text{--}40$ km. Cooling thickens the elastic core. For oceanic lithosphere, $T_e$ closely tracks the depth to the 300–600°C isotherm, which deepens as the plate ages according to the half-space cooling model:$T_e \propto \sqrt{t}$.
Continental Lithosphere
$T_e \approx 20\text{--}100$ km. Highly variable due to differences in thermal state, composition, and tectonic history. Cratons (old, cold, thick lithosphere) have the largest$T_e$ values (~60–100 km), while young or thermally perturbed continental lithosphere can have $T_e < 20$ km.
Flexural Parameter & Analytical Solutions
The characteristic length scale of flexural deformation is the flexural parameter (or flexural wavelength parameter) $\alpha$:
Flexural parameter:
\[ \alpha = \left[ \frac{4D}{\Delta\rho \, g} \right]^{1/4} \]
Typical values: $\alpha \approx 50\text{--}70$ km for young oceanic lithosphere,$\alpha \approx 100\text{--}250$ km for old oceanic and continental lithosphere. The flexural wavelength is $\lambda \approx 2\pi\alpha$.
For a line load $V_0$ (force per unit length) applied to a continuous elastic plate, the analytical solution for the deflection is:
Line load deflection (continuous plate):
\[ w(x) = \frac{V_0 \, \alpha^3}{8D} \, e^{-|x|/\alpha} \left[ \cos\!\left(\frac{|x|}{\alpha}\right) + \sin\!\left(\frac{|x|}{\alpha}\right) \right] \]
The deflection decays exponentially with distance from the load, modulated by oscillatory (cosine and sine) terms that produce the characteristic forebulge.
Moat (Deflection Maximum)
The maximum deflection occurs directly beneath the load at $x = 0$, where$w_0 = V_0 \alpha^3 / (8D)$. This is the deepest part of the flexural depression — corresponding to the trench at a subduction zone or the deepest part of a foreland basin.
Forebulge
At a distance of approximately $x \approx \pi\alpha$ from the load, the deflection becomes negative (upward), creating a low-amplitude topographic rise called the forebulge. The amplitude of the forebulge is only about 4% of the maximum deflection, but it is a diagnostic feature of flexural loading. In foreland basins, the forebulge often forms a topographic high that separates the foredeep from the back-bulge basin.
Broken Plate vs. Continuous Plate
The flexure equation has different solutions depending on the boundary conditions at the load. Two end-member cases are commonly considered:
Continuous Plate
The plate extends beneath and beyond the load without a break. The deflection and slope are continuous everywhere. Produces a symmetric deflection profile about the load axis. Appropriate for loads applied to the interior of a plate (e.g., seamounts, ice sheets).
Broken Plate
The plate is broken (free edge) at the point of loading. Only one side of the plate is loaded, and the boundary conditions specify zero bending moment at the free end. The broken plate produces a deeper, narrower deflection and a more prominent forebulge. This is the appropriate model for subduction zones and foreland basins where the plate is broken at a thrust front.
The broken plate solution predicts a maximum deflection that is twice that of the continuous plate for the same applied load, because the load is supported by flexure on only one side. In practice, the real lithosphere often lies between these two end members, with partial mechanical coupling across faults.
Gravity Admittance & Constraining Te
The gravity admittance (or transfer function) describes the relationship between topography and gravity in the spectral domain. It is the primary observational tool for determining the effective elastic thickness $T_e$ of the lithosphere.
For a flexurally supported topographic load with wavenumber $k$, the theoretical free-air admittance is:
Flexural admittance (simplified):
\[ Z(k) = 2\pi G \rho_c \left(1 - e^{-k t_c} \cdot \frac{\Delta\rho}{\rho_c} \cdot W(k) \right) \]
where $W(k) = 1/(1 + Dk^4/(\Delta\rho \, g))$ is the isostatic response function and $t_c$ is the crustal thickness. By comparing observed admittance from satellite gravity and topography data with theoretical predictions for different $T_e$ values, the best-fitting elastic thickness can be determined.
Short-wavelength loads (high $k$) are supported by the plate's elastic strength and are uncompensated (admittance approaches the Bouguer slab value $2\pi G \rho_c$). Long-wavelength loads (low $k$) exceed the plate's capacity to support them and are isostatically compensated (admittance approaches zero). The transition wavelength between these regimes is controlled by $T_e$: a stronger plate (larger $T_e$) supports loads to longer wavelengths before compensation becomes significant.
Foreland Basins
Foreland basins are sedimentary basins that form on the undeformed plate adjacent to a mountain belt, created by the flexural subsidence of the lithosphere under the weight of the orogenic load. They are among the most important hydrocarbon provinces and sediment repositories on Earth, and their geometry provides direct constraints on the flexural properties of the underlying lithosphere.
Foredeep
The deepest part of the foreland basin, immediately adjacent to the thrust front. Sediment thickness can exceed 5–8 km. The foredeep migrates toward the foreland as the thrust front advances and the load progresses onto previously undeformed lithosphere.
Forebulge
The flexural forebulge is the low-amplitude topographic high cratonward of the foredeep, located approximately $\pi\alpha$ from the thrust front. Forebulge erosion and unconformities are common features. In the Alps, the Jura Mountains have been interpreted as sitting on the forebulge of the Alpine foreland.
Back-Bulge Basin
A shallow, broad depression cratonward of the forebulge. It is a region of subtle subsidence caused by the oscillatory decay of the flexural deflection. Sediment thicknesses are typically modest (hundreds of meters) but can cover very large areas.
Peripheral Foreland Basins
Form on the subducting or underthrusting plate at a continent–continent collision zone. The orogenic load is the colliding mountain belt itself. Examples: Indo-Gangetic plain (Himalayan foredeep), Molasse Basin (Alpine foredeep), Persian Gulf (Zagros foredeep).
Retroarc Foreland Basins
Form on the overriding plate behind a volcanic arc at an ocean–continent subduction zone. The orogenic load is the thickened crust of the magmatic arc and fold-thrust belt. Examples: Andean foreland basins (Chaco, Gran Chaco), Western Canada Sedimentary Basin (Cretaceous Seaway).
The Himalayan foreland basin is the type example of a peripheral foreland. The Indo-Gangetic Plain represents the foredeep, filled with up to 6 km of Cenozoic sediment shed from the Himalaya. The Ganges and Indus river systems transport this enormous sediment flux to the sea, building the Bengal and Indus submarine fans — the two largest sediment accumulations on Earth. The depocenter has migrated southward over time as the Himalayan thrust front has advanced onto the Indian plate.
Key Numbers
5–100 km
Range of Te
10²¹–10²&sup5; N·m
Range of Flexural Rigidity
50–250 km
Flexural Parameter α
~6 km
Himalayan Foredeep Sediment