Lithostatic Pressure Calculator
Estimate overburden pressure from rock density, depth, and gravity using the standard geomechanics equation P = rho x g x h.
Results
Enter your values and click calculate.
Expert Guide: How to Calculate Lithostatic Pressure Accurately
Lithostatic pressure is the pressure exerted by the weight of overlying rock at depth. In structural geology, geomechanics, basin modeling, drilling, and metamorphic petrology, this pressure is a foundational variable because it controls stress state, compaction, porosity evolution, and many deformation processes. If you work in subsurface engineering, geothermal exploration, petroleum systems, carbon storage, mining, or tectonics, knowing how to calculate lithostatic pressure correctly is essential for reliable interpretation and safe operations.
The standard calculation is straightforward in form but can become sensitive to assumptions in real settings. The simplest equation is:
P = rho x g x h, where P is pressure (Pa), rho is average overburden density (kg/m3), g is gravitational acceleration (m/s2), and h is depth (m).
This formula assumes a vertical, static overburden with constant average density. In nature, density changes with lithology, compaction, fluid content, and temperature. Yet even with these complexities, the equation remains the core framework. Most advanced workflows still begin from this relation and then integrate depth dependent density functions, layered stratigraphy, or geophysical constraints.
Why lithostatic pressure matters in practice
- Well design and drilling safety: overburden stress helps define mud window limits and geomechanical stability.
- Reservoir geomechanics: stress state influences fault reactivation risk and fracture behavior.
- Metamorphic geology: pressure estimates constrain burial history and P-T paths.
- Basin analysis: compaction and porosity predictions require realistic overburden pressure profiles.
- Planetary geology: lower gravity bodies like Mars and the Moon produce significantly lower pressure at equal depth.
Step by step method for accurate calculations
- Choose depth reference carefully. Decide whether depth is measured from ground surface, seafloor, or another geodetic datum.
- Select appropriate density. Use bulk density, not grain density. For mixed lithologies, use a weighted average or layered model.
- Convert units to SI. Depth in meters and density in kg/m3 reduce conversion errors.
- Use the correct gravity value. Earth standard is 9.80665 m/s2; planetary applications require site specific values.
- Compute pressure in pascals. Convert to MPa or psi as needed for engineering use.
- Sanity check against expected gradients. On Earth, common lithostatic gradients are often around 20 to 30 MPa/km depending on density.
Typical rock densities and pressure gradients
Because the gradient dP/dz equals rho x g, density directly scales pressure increase per unit depth. The table below uses Earth gravity and representative bulk densities commonly used in early stage models.
| Material or setting | Typical bulk density (kg/m3) | Lithostatic gradient (MPa/km, Earth) | Notes |
|---|---|---|---|
| Unconsolidated sediment | 2200 | 21.6 | Shallow basins, high porosity intervals |
| Sandstone and shale average | 2500 | 24.5 | Common basin scale first pass estimate |
| Granitic crust average | 2700 | 26.5 | Typical continental upper crust estimate |
| Basaltic crust | 3000 | 29.4 | Dense mafic crust and volcanic provinces |
| Upper mantle peridotite | 3300 | 32.4 | Common mantle scale approximation |
Worked examples
Suppose average bulk density is 2700 kg/m3 and depth is 3 km on Earth. Convert depth to meters: 3000 m. Then:
P = 2700 x 9.80665 x 3000 = 79,433,865 Pa = 79.43 MPa
If you keep the same depth and density on Mars (g = 3.71 m/s2), pressure is only:
P = 2700 x 3.71 x 3000 = 30,051,000 Pa = 30.05 MPa
This large difference illustrates why gravity must be explicit in planetary or comparative geology calculations.
Depth pressure benchmarks for crustal interpretation
| Depth (km) | Pressure at 2700 kg/m3 (MPa) | Pressure (psi) | Common geological context |
|---|---|---|---|
| 1 | 26.5 | 3,843 | Shallow burial, early compaction |
| 5 | 132.4 | 19,204 | Many hydrocarbon source intervals |
| 10 | 264.8 | 38,409 | Deep sedimentary basins, brittle-ductile transition vicinity in some terrains |
| 20 | 529.6 | 76,819 | Lower crustal levels in many continental settings |
| 35 | 926.8 | 134,433 | Approximate pressure near average continental Moho depth in many regions |
Common mistakes and how to avoid them
- Mixing true vertical depth and measured depth: use vertical depth for stress calculations unless you intentionally model trajectory effects.
- Using grain density instead of bulk density: grain density can overestimate stress in porous rocks.
- Ignoring density increase with depth: shallow averages may underpredict deep pressures if compaction is significant.
- Unit conversion errors: ft to m and g/cm3 to kg/m3 are frequent error points.
- Confusing lithostatic and pore pressure: they are related but not identical; effective stress depends on both.
Layered calculation for higher fidelity models
In real basins and crustal profiles, density often varies with depth. Instead of one average rho, divide the column into layers and sum each layer contribution:
P_total = sum(rho_i x g x h_i)
This is simple to implement in spreadsheets or code and often yields much better agreement with log calibrated density profiles. If density logs or seismic inversion outputs are available, you can sample rho(z) and numerically integrate pressure. That is standard in many advanced geomechanics workflows.
Lithostatic pressure versus hydrostatic and overpressure systems
Hydrostatic pressure is fluid pressure from a connected water column, while lithostatic pressure is total vertical stress from rock overburden. In compacting sedimentary basins, pore pressure can deviate from hydrostatic and approach lithostatic in extreme overpressure zones, but these are conceptually different quantities. Engineers often compare them using pressure gradients. Hydrostatic gradients in freshwater settings are near 9.8 MPa/km, whereas lithostatic gradients for typical rocks are often around 22 to 30 MPa/km on Earth.
Effective stress, central to deformation and failure, is total stress minus pore pressure (in simplified scalar form). Therefore, a robust lithostatic estimate is one half of the stress interpretation; accurate pore pressure estimation is the other half.
Data quality and uncertainty management
Any pressure estimate is only as good as the density and depth model behind it. A practical way to handle uncertainty is to run low, base, and high scenarios. For example, if density uncertainty is plus or minus 100 kg/m3 at 4 km depth, pressure uncertainty can be on the order of plus or minus 3.9 MPa on Earth. Scenario ranges improve risk communication for drilling windows, seal integrity assessment, and fault stability studies.
- Use density logs where available and calibrate with core measurements.
- Check stratigraphic tops and depth datum consistency across datasets.
- Document assumptions for gravity, lithology, and compaction state.
- Track conversion precision, especially when reporting MPa and psi together.
Authoritative references for deeper study
For trusted background data and unit standards, review resources from government and university institutions:
- USGS density fundamentals and measurement context
- NIST SI unit guidance and consistent unit usage
- MIT OpenCourseWare geophysics material for stress and subsurface interpretation
Final takeaways
Calculating lithostatic pressure is conceptually simple, but professional grade results require disciplined input handling. Always align depth reference, density definition, gravity, and units before calculating. Use a single layer estimate for quick screening, then move to layered integration when project decisions depend on tighter accuracy. With these practices, lithostatic pressure becomes a reliable building block for geomechanical interpretation, basin modeling, and subsurface risk management.