The lithosphere bends elastically under loads (mountains, sediment basins) with a characteristic bending length (related to elastic thickness). Flexural models, validated by gravity and topography, estimate effective elastic thickness as a function of age and temperature.
Use forward modeling to simulate basin geometry under varying loads. Compare predictions to observed topography and gravity to invert elastic thickness.
From isostasy, you know that the lithosphere floats on the denser asthenosphere, and that loads on the surface — mountains, ice sheets, sediment piles — must be compensated by displacement of mantle material below. But the simple Airy model treats the lithosphere as if it has no strength: each column sinks independently, like blocks of wood floating in water. Real lithosphere is not that weak. It has rigidity, and it bends as a coherent plate rather than sinking in disconnected columns. This bending behavior is lithospheric flexure, and it changes the geometry of isostatic compensation in important ways.
Think of the lithosphere as an elastic beam resting on a fluid foundation (the asthenosphere). When you place a point load on a beam — say, a volcanic island — the beam does not just sink directly beneath the load. It bends over a broad region: it deflects downward under the island, creating a surrounding moat (a flexural depression), and bows slightly upward farther away, forming a flexural bulge. The Hawaiian Islands are a textbook example: the seafloor is depressed in an arc around each island and slightly elevated in a ring beyond. The width and amplitude of this deflection pattern depend on a single key parameter: the flexural rigidity of the plate, which is controlled by its effective elastic thickness (Tₑ).
Effective elastic thickness is not the same as the total thickness of the lithosphere — it represents the thickness of an idealized perfectly elastic plate that would produce the same bending. Young, hot oceanic lithosphere near a mid-ocean ridge might have Tₑ of only 5–10 km because the rock is warm and weak. Old, cold oceanic lithosphere can have Tₑ of 30–40 km. Continental lithosphere varies widely (10–100+ km) depending on thermal state and composition. The governing equation is the flexural equation: D∇⁴w + (ρ_m − ρ_fill)gw = q(x), where D is flexural rigidity (proportional to Tₑ³), w is deflection, ρ_m and ρ_fill are mantle and infill densities, g is gravity, and q is the applied load. Larger D means the plate distributes loads over a wider area; smaller D means the deflection is narrow and deep, approaching the Airy limit.
In practice, Tₑ is estimated by comparing observed topography and gravity anomalies to predictions from flexural models. A sedimentary basin next to a mountain belt, for instance, has a shape controlled by the flexural response to the mountain load. If you forward-model the basin geometry for different values of Tₑ and find the one that best matches the observed basin width and depth, you have constrained the plate's strength. This approach links surface observables — topography, gravity, basin stratigraphy — directly to the mechanical and thermal properties of the lithosphere, making flexural analysis one of the most powerful tools in geodynamics.