Isostasy Calculator - Airy Pratt
Calculate isostatic equilibrium using Airy and Pratt models for crustal root depth, density variations, and compensation depth.
Isostasy Calculator - Airy Pratt
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: Airy: Root = h x (rho_c / (rho_m - rho_c)); Pratt: rho = rho_ref x D / (D + h)
Worked example โ Root: 22.50 km | Total Crust: 57.50 km | Moho: 57.50 km depth
Formula
Airy: Root = h x (rho_c / (rho_m - rho_c)); Pratt: rho = rho_ref x D / (D + h)
In the Airy model, root depth depends on elevation and the density contrast between crust and mantle. In the Pratt model, crustal density varies so that all columns have equal mass per unit area down to a uniform compensation depth D.
Worked Examples
Example 1: Airy Model: Himalayan Mountain Root
Problem:Calculate the crustal root depth for a mountain with 5 km elevation using the Airy model. Crustal density = 2.7 g/cm3, mantle density = 3.3 g/cm3, normal crustal thickness = 35 km.
Solution:Density contrast = 3.3 - 2.7 = 0.6 g/cm3 Root depth = elevation x (crustal density / density contrast) Root = 5 x (2.7 / 0.6) = 5 x 4.5 = 22.5 km Total crustal thickness = 35 + 22.5 = 57.5 km Moho depth = 57.5 km below surface
Result:Root: 22.50 km | Total Crust: 57.50 km | Moho: 57.50 km depth
Example 2: Pratt Model: Mid-Ocean Ridge Density
Problem:Calculate the crustal density beneath a 2 km elevated mid-ocean ridge using the Pratt model. Reference density = 2.8 g/cm3, compensation depth = 100 km.
Solution:Column density = reference density x D / (D + h) Column density = 2.8 x 100 / (100 + 2) = 280 / 102 = 2.7451 g/cm3 Density reduction = 2.8 - 2.7451 = 0.0549 g/cm3 Mass balance check: 2.8 x 100 = 280; 2.7451 x 102 = 280.0 (balanced)
Result:Column Density: 2.7451 g/cm3 | Density Reduction: 0.0549 g/cm3
Frequently Asked Questions
What is isostasy and why is it important in geology?
Isostasy is the gravitational equilibrium between the Earth's lithosphere and asthenosphere, where lighter crustal blocks float on the denser mantle material below, similar to how icebergs float in water. This concept is fundamental to understanding why mountains have deep crustal roots, why continents stand higher than ocean floors, and how the Earth's surface responds to loading and unloading events. When weight is added to the crust through ice sheets, sediment deposition, or volcanic buildup, the crust sinks into the mantle. When weight is removed through erosion or ice sheet melting, the crust rebounds upward, a process called isostatic adjustment. This explains phenomena like post-glacial rebound in Scandinavia and Canada, where land is still rising centuries after ice age glaciers melted.
How does the Airy isostasy model work?
The Airy model of isostasy, proposed by George Biddell Airy in 1855, assumes that the Earth's crust has a uniform density but varies in thickness. Mountains are supported by deep crustal roots that extend into the denser mantle, much like icebergs with deeper keels floating higher above water. The root depth is calculated as the elevation multiplied by the crustal density divided by the density contrast between mantle and crust. For typical values of 2.7 grams per cubic centimeter for crust and 3.3 for mantle, the density contrast is 0.6, meaning each kilometer of elevation requires approximately 4.5 kilometers of crustal root. This model successfully explains observations from seismic studies showing that mountain ranges like the Himalayas and Andes have significantly thicker crust beneath them.
How does the Pratt isostasy model differ from Airy?
The Pratt model, proposed by John Henry Pratt in 1855, takes a fundamentally different approach from Airy. Instead of varying crustal thickness with uniform density, Pratt assumes a uniform compensation depth but varying crustal density. Higher topography is underlain by less dense rock, while lower topography has denser rock. All crustal columns extend to the same compensation depth, typically around 100 kilometers. The density of each column is calculated so that the mass per unit area is equal everywhere, achieving isostatic equilibrium. The Pratt model works well for explaining mid-ocean ridges where thermal expansion reduces rock density and creates elevated seafloor. In practice, both Airy and Pratt mechanisms operate simultaneously in different geological settings.
What is the Moho discontinuity and how does isostasy affect it?
The Mohorovicic discontinuity, commonly called the Moho, is the boundary between the Earth's crust and the underlying mantle, defined by a sharp increase in seismic wave velocities. Isostasy directly controls the depth of the Moho because crustal roots extend the crust deeper into the mantle beneath elevated terrain. Under continental plains, the Moho typically lies at about 30 to 40 kilometers depth. Beneath major mountain ranges, isostatic compensation pushes the Moho to 60 to 80 kilometers deep, as observed beneath the Tibetan Plateau where the crust is approximately 70 kilometers thick. Under oceanic crust, the Moho is much shallower at only 5 to 10 kilometers depth. Seismic refraction surveys confirm these depth variations and provide evidence supporting isostatic models.
How is isostatic rebound observed in the real world?
Isostatic rebound, or glacial isostatic adjustment, is one of the most dramatic observable effects of isostasy occurring over human timescales. During the last ice age, massive ice sheets up to 3 kilometers thick depressed the crust beneath Scandinavia, Canada, and other glaciated regions. Since the ice melted approximately 10,000 to 20,000 years ago, these regions have been slowly rebounding. Scandinavia is currently rising at rates up to 10 millimeters per year near the Gulf of Bothnia. Hudson Bay in Canada is still rising at about 12 millimeters per year. This rebound has measurable effects: ancient shorelines are now found hundreds of meters above current sea level, and sea level gauge records show apparent sea level fall in rebounding regions. The ongoing adjustment will continue for thousands of years until full equilibrium is reached.
What is the difference between local and regional isostatic compensation?
Airy and Pratt are both local models: each vertical column is treated as floating independently, with no strength coupling it to its neighbours. Real lithosphere behaves as an elastic plate over the ductile asthenosphere, so a load is partly supported by flexure of the surrounding region. Vening Meinesz formalised this as regional compensation in 1931, characterised by an effective elastic thickness, typically a few kilometres for young hot lithosphere and several tens of kilometres for old cratonic lithosphere. Narrow loads such as individual seamounts are held up mainly by plate bending, while broad loads hundreds of kilometres across approach local Airy behaviour.
How do free-air and Bouguer gravity anomalies show whether a region is compensated?
The free-air anomaly corrects observed gravity only for elevation, using a gradient of about 0.3086 milligal per metre. Over terrain that is fully compensated at depth, the excess mass of the topography and the mass deficiency of the root or low-density column largely cancel, so the long-wavelength free-air anomaly stays near zero. The Bouguer anomaly additionally removes the attraction of the rock between the station and sea level, roughly 0.0419 times density in grams per cubic centimetre, in milligal per metre. Because that correction strips out the topography but not the compensating root, compensated mountain belts show strongly negative Bouguer anomalies, often exceeding minus 300 milligal.
Why do ocean basins need an anti-root rather than a crustal root?
In the Airy formulation the root scales with elevation, so a negative elevation gives a negative root, meaning the crust thins rather than thickens. Physically, a basin has less crustal mass above the compensation level, and the deficit is made up by mantle rising closer to the surface, which is why the Moho sits at only about 5 to 10 kilometres beneath oceanic crust. A submarine basin is also loaded by seawater rather than air, so a rigorous ocean calculation replaces the missing rock with water at roughly 1.03 grams per cubic centimetre, which reduces the required crustal thinning compared with a dry basin of the same depth.
How quickly does the crust return to isostatic equilibrium after a load is removed?
Isostatic models describe the final equilibrium state, not the path to it. The rate is set by how fast mantle material can flow laterally beneath the load, and it depends on both mantle viscosity and the wavelength of the load, with the sense of that wavelength dependence set by the viscosity structure. For a uniform viscous half-space the classical relaxation time varies as viscosity divided by the product of density, gravity and load wavelength, so broad loads relax faster than narrow ones. If flow is instead channelled through a low-viscosity asthenospheric layer beneath stiffer mantle, the dependence reverses and long wavelengths relax more slowly. That is precisely why rebound histories are inverted to recover the viscosity profile rather than assumed from it. Sediment loading in a subsiding basin is slow enough that the crust stays close to equilibrium throughout.
Which isostasy model should I use for a given geological setting?
Choose Airy where topography is supported by crustal thickening, which covers collisional orogens and continental margins: seismic refraction there shows the Moho deepening in step with surface elevation. Choose Pratt where elevation reflects lateral density change at nearly constant thickness, most clearly at mid-ocean ridges, where hot young lithosphere is thermally expanded and therefore lighter. Continental hotspot swells and thermally uplifted plateaus also fit Pratt better. In practice the two mechanisms overlap, so the useful test is whether independent seismic data show a thickened crust or a normal crust of anomalous density.
References
Background & Theory
History
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer ยท Editorial policy
Related Calculators
๐งฎEarth Energy Balance Calculator
Calculate earth energy balance with inputs, formulas, and instant results.
๐งฎBouguer Correction Calculator
Calculate bouguer correction with inputs, formulas, and instant results.
๐งฎCrustal Density Calculator
Calculate crustal density with inputs, formulas, and instant results.
๐งฎCrustal Thickness Calculator (H-k)
Calculate crustal thickness from receiver functions with inputs, formulas, and instant results.
๐งฎEarthquake Magnitude to Energy Calculator
Calculate earthquake magnitude to energy with inputs, formulas, and instant results.
๐งฎGutenberg-Richter Recurrence Calculator
Calculate earthquake recurrence (gutenbergโrichter) with inputs, formulas, and instant results.
๐งฎEarthโs Rotation Period Variation Calculator
Calculate earthโs rotation period variation with inputs, formulas, and instant results.
๐งฎElastic Moduli Converter (e, G, K, ฮ)
Calculate elastic moduli converter (e, g, k, ฮฝ) with inputs, formulas, and instant results.