{ "id": "epabh7bjfw", "slug": "quasistatic-approximation", "title": "Quasistatic approximation", "updated": "2026-09-30", "topic_path": [ { "id": "physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/topics/physical" }, { "id": "physical.earth", "label": "Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.earth" }, { "id": "physical.earth.planet", "label": "Earth systems and geophysics", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.earth.planet" } ], "geo": [ { "id": "geo.nongeo.t1800.physical.earth", "label": "Non-geographic · 1800 to 1945: Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1800.physical.earth", "path": [ { "id": "geo.nongeo", "label": "Non-geographic", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo" }, { "id": "geo.nongeo.t1800", "label": "Non-geographic · 1800 to 1945", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1800" }, { "id": "geo.nongeo.t1800.physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1800.physical" }, { "id": "geo.nongeo.t1800.physical.earth", "label": "Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1800.physical.earth" } ] }, { "id": "geo.nongeo.t1946.physical.earth.planet", "label": "Non-geographic · 1946 to 2000: Earth systems and geophysics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1946.physical.earth.planet", "path": [ { "id": "geo.nongeo", "label": "Non-geographic", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo" }, { "id": "geo.nongeo.t1946", "label": "Non-geographic · 1946 to 2000", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1946" }, { "id": "geo.nongeo.t1946.physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1946.physical" }, { "id": "geo.nongeo.t1946.physical.earth", "label": "Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1946.physical.earth" }, { "id": "geo.nongeo.t1946.physical.earth.planet", "label": "Earth systems and geophysics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1946.physical.earth.planet" } ] } ], "excerpt": "The quasistatic approximation is a geophysics modeling simplification that neglects inertial terms in the momentum balance, underpinning poroelastic consolidation, crustal deformation, ice-sheet flow, and atmospheric dynamics.", "snippet": "The quasistatic approximation is a geophysics modeling simplification that neglects inertial terms in the momentum balance, underpinning poroelastic consolidation, crustal deformation, ice-sheet flow, and atmospheric dynamics.", "node": "physical.earth.planet", "markdown": "# Quasistatic approximation\n\nThe quasistatic approximation is a modeling simplification in geophysics and related earth sciences that neglects inertial terms in the momentum balance when a process evolves slowly compared with the elastic wave timescale, reducing dynamic momentum balance to a quasi-equilibrium, or elliptic, problem while time dependence may remain in loading, constitutive behavior, and other governing equations. It underlies standard treatments of poroelastic consolidation, crustal deformation over geologic time, ice-sheet and ice-shelf flow, and quasi-geostrophic atmospheric dynamics.[1](https://www.sandia.gov/files/sierra/SM_Theory_5_30/main/quasistatics.html) • [2](https://onlinelibrary.wiley.com/doi/10.1002/nme.6671) • [3](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C)\n\n| Key fact | Detail |\n|---|---|\n| What is neglected | Inertial (acceleration) terms in the momentum balance, justified when inertial forces are negligible compared with internal and applied forces[1](https://www.sandia.gov/files/sierra/SM_Theory_5_30/main/quasistatics.html) |\n| Reduced problem | A static balance \\( [F]^{int}(d(t)) = [F]^{ext} \\) solved at each time; only one initial condition remains and time becomes a load parameter[1](https://www.sandia.gov/files/sierra/SM_Theory_5_30/main/quasistatics.html) |\n| Small parameters | Mach number (ice velocity over elastic wave speed) in glaciology; Rossby number in quasi-geostrophic theory; aspect ratio in hydrostatic primitive equations[3](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C) • [4](https://pordlabs.ucsd.edu/wryoung/theorySeminar/pdf14/VallisQGderivation.pdf) |\n| Slip-rate validity | Aseismic slip below about 0.1 mm/s makes accelerations negligible; slip above about 1 cm/s is seismic[5](https://doi.org/10.31224/osf.io/ps2ve) |\n| Computational gain | No CFL time-step limit; implicit quasistatic solvers take large steps limited only by accuracy |\n| Main failure mode | Quasistatic slip rates diverge as a seismic event approaches; quasi-dynamic models underpredict rupture speed and can mispredict earthquake size and recurrence[6](https://tectonics.caltech.edu/publications/pdf/Thomas-Dynamic_vsQuasidynamic-JGR-2014.pdf) |\n\n## How it works\n\nThe approximation rests on a comparison of forces. It is appropriate when inertial forces are negligible compared with the internal and applied forces in a system; what counts as negligible relies on intuition, and numerical experimentation is one way to gain that intuition.[1](https://www.sandia.gov/files/sierra/SM_Theory_5_30/main/quasistatics.html) Formally, the justification is a timescale separation. In ice-sheet modeling the requirement is that the ratio of ice velocity to elastic wave speed, the [Mach number](https://www.edgepedia.org/mach-number), remain small; over sub-daily timescales the quasi-static-creep approximation fails and elastic effects become important.[3](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C) In atmospheric dynamics, quasi-geostrophy is the lowest-order model in an asymptotic expansion in the small [Rossby number](https://www.edgepedia.org/rossby-number), and the hydrostatic primitive equations follow from the [Navier–Stokes equations](https://www.edgepedia.org/navier-stokes-equations) with the aspect ratio as the small parameter.[4](https://pordlabs.ucsd.edu/wryoung/theorySeminar/pdf14/VallisQGderivation.pdf)\n\nDropping the acceleration term changes the character of the equations. The momentum balance becomes elliptic, so information propagates infinitely fast within the model; reintroducing acceleration and visco-elasticity in ice yields a hyperbolic system in which information travels at a finite elastic wave speed.[3](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C) Because the lowest-order asymptotic model is a limit of the original equation set, it preserves invariants such as energy and potential vorticity.[4](https://pordlabs.ucsd.edu/wryoung/theorySeminar/pdf14/VallisQGderivation.pdf)\n\n## How it is done\n\nThe practitioner workflow has three steps. First, nondimensionalize the governing equations to expose the small parameter (Mach number, Rossby number, aspect ratio, or a slip-rate ratio). Second, drop the inertial term. In a discrete formulation, omitting the inertial term in the equations of motion yields the quasistatic problem \\( [F]^{int}(d(t)) = [F]^{ext} \\), a static balance evaluated at each time; only the single initial condition \\( d(0) = d_0 \\) remains, and time may be a generic parameterization of the loads rather than a physical clock.[1](https://www.sandia.gov/files/sierra/SM_Theory_5_30/main/quasistatics.html) Equivalently, the inertial term is ignored and time dependence enters only through constitutive models and loading conditions, so a quasistatic simulation is a series of static problems with time-varying properties and boundary conditions. Third, solve the reduced boundary-value problem and check consistency.\n\nConsistency checks are problem-specific. For consolidation, analytical one-dimensional examples have been used to show that for low-frequency soil-mechanics applications the complete Biot theory does not significantly differ from a simplified form, but also that for harmonic loadings inertia cannot be neglected even in the long-time behavior.[7](https://www.tugraz.at/fileadmin/user_upload/Institute/AM-BM/Files/Forschung/Preprints/preprint_2008_03.pdf)\n\n## Origin\n\nEarly theories of porous media saturated by viscous fluid were established as quasi-static theories, with no inertia effects taken into account; building on Terzaghi's work, a theoretical description of such porous materials was presented by Maurice A. Biot in his 1941 paper *General Theory of Three-Dimensional Consolidation* in the [Journal of Applied Physics](https://www.edgepedia.org/journal-of-applied-physics).[7](https://www.tugraz.at/fileadmin/user_upload/Institute/AM-BM/Files/Forschung/Preprints/preprint_2008_03.pdf) • [8](https://doi.org/10.1063/1.1712886) Biot's quasi-static consolidation equations are obtained by neglecting body forces and inertial terms, and they form the u-p formulation used as the basis of quasi-static finite-element modeling of wave-induced fluid flow.[9](https://marcelfrehner.ch/resourcen/Quintal_et_al_JGR_2011.pdf) Dynamic poroelastodynamics developed later.[7](https://www.tugraz.at/fileadmin/user_upload/Institute/AM-BM/Files/Forschung/Preprints/preprint_2008_03.pdf) In earthquake-cycle modeling, a fully dynamic spectral formulation for tectonic loading with spontaneous rupture on faults with rate- and state-dependent friction was reported by Nadia Lapusta and colleagues in 2000 in the [Journal of Geophysical Research: Solid Earth](https://www.edgepedia.org/journal-of-geophysical-research-solid-earth).[10](https://doi.org/10.1029/2000jb900250)\n\n## Variants\n\n**Quasistatic poroelasticity** retains elastic equilibrium coupled to pore-pressure diffusion; Biot's 1941 equations give three equations in four unknowns (the displacements u, v, w, and the stress \\( \\sigma \\)), closed by the continuity or storage equation.[8](https://doi.org/10.1063/1.1712886) Associated quasi-static poroelastic parameters include Skempton's coefficient, the ratio of pore-pressure increment to mean-stress increment under undrained conditions, and the Biot-Willis parameter, which serve equally well for the undrained bulk modulus and [Poisson's ratio](https://www.edgepedia.org/poissons-ratio).[11](https://link.springer.com/article/10.1007/BF00998332)\n\n**Quasi-dynamic rupture** prescribes stress and frictional conditions on the fault but approximates wave propagation by ignoring inertia and adding a radiation damping term to the equation of motion; fully dynamic rupture models instead generate the whole wavefield by including inertia.[12](https://par.nsf.gov/servlets/purl/10380873) Without the damping term \\( V/(2c_{s}) \\), the quasi-dynamic procedure would turn into a quasi-static one and would not allow solutions during fast, inertially controlled slip.[6](https://tectonics.caltech.edu/publications/pdf/Thomas-Dynamic_vsQuasidynamic-JGR-2014.pdf)\n\n**The shallow shelf approximation (SSA)** in glaciology results from neglecting the material time derivative in the ice-shelf momentum balance, yielding a non-linear elliptic equation.[3](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C) **Quasi-geostrophic balance** in atmospheric dynamics requires quasi-geostrophic horizontal velocity and quasi-hydrostatic pressure, and filters out unwanted motions such as sound waves that would otherwise obstruct prediction.[13](https://atmos.uw.edu/academics/classes/2010Q4/441/ch6.pdf)\n\n## Applications\n\n**Seismic-cycle and fault modeling.** Quasi-dynamic simulators resolve interseismic, nucleation, post-seismic, and dynamic rupture phases of earthquake cycles.[14](https://github.com/dunyuliu/EQquasi) Benchmark exercises such as the SEAS code comparison compare codes that incorporate full elastodynamic effects with quasi-static treatments.[15](https://par.nsf.gov/servlets/purl/10417587)\n\n**Crustal and salt tectonics.** By multiscale asymptotics, the inertia term in the momentum balance can be safely neglected when tracking crustal deformation over long geologic timescales; a blended transient/quasistatic scheme enforcing stress static equilibrium improves efficiency over transient dynamics algorithms, which are forced to resolve seismic events over geologic timescales.[2](https://onlinelibrary.wiley.com/doi/10.1002/nme.6671)\n\n**Rock physics and ice.** Biot's quasi-static consolidation equations underpin quasi-static finite-element modeling of seismic attenuation due to wave-induced fluid flow.[9](https://marcelfrehner.ch/resourcen/Quintal_et_al_JGR_2011.pdf) The SSA is the workhorse elliptic model for ice-shelf flow.[3](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C) [Quasi-geostrophic theory](https://www.edgepedia.org/quasi-geostrophic-theory) is standard for synoptic-scale atmospheric disturbances.[13](https://atmos.uw.edu/academics/classes/2010Q4/441/ch6.pdf)\n\n## Limitations and alternatives\n\nThe approximation fails whenever accelerations matter. In the quasi-static formulation of fault slip, slip rates become infinite as a seismic event approaches, so the method cannot describe fast seismic slip.[6](https://tectonics.caltech.edu/publications/pdf/Thomas-Dynamic_vsQuasidynamic-JGR-2014.pdf) Quantitatively, one review places the validity boundary for aseismic, creeping slip at about 0.1 mm/s, with slip faster than about 1 cm/s classified as seismic.\n\n**Quasi-dynamic errors are systematic.** Without additional coseismic weakening, quasi-dynamic (QD) and fully dynamic (FD) simulations give qualitatively similar slip patterns, but QD produces slower slip velocities and rupture speeds and more rupture arrest at velocity-strengthening patches; one comparison found an average rupture speed of 3.56 km/s in FD versus 0.98 km/s in QD.[6](https://tectonics.caltech.edu/publications/pdf/Thomas-Dynamic_vsQuasidynamic-JGR-2014.pdf) With additional coseismic weakening, the two approaches diverge qualitatively: near-periodic pulse-like FD events versus much larger crack-like QD events. Ignoring transient wave-mediated stress transfers may mispredict earthquake size and recurrence, average fault stress levels, and postseismic slip; seismic waves can also promote local fault weakening and modify rupture speed.[6](https://tectonics.caltech.edu/publications/pdf/Thomas-Dynamic_vsQuasidynamic-JGR-2014.pdf) • [12](https://par.nsf.gov/servlets/purl/10380873)\n\n**Switching thresholds are not standardized.** One spectral-element study switches from quasi-static to dynamic at 0.5 mm/s and back at 0.2 mm/s,[16](https://tectonics.caltech.edu/publications/pdf/Kaneko_JGR2011.pdf) while the SEAS benchmark uses a scheme switching to dynamic when \\( \\max(V) > 10 \\) mm/s and back when \\( \\max(V) < 1 \\) mm/s.[15](https://par.nsf.gov/servlets/purl/10417587) Published sources do not reconcile these choices.\n\n**Alternatives.** Fully dynamic rupture models retain inertia and the whole wavefield.[12](https://par.nsf.gov/servlets/purl/10380873) Blended transient/quasistatic schemes enforce static equilibrium while resolving transients where needed.[2](https://onlinelibrary.wiley.com/doi/10.1002/nme.6671) Pseudo-transient ([dynamic relaxation](https://www.edgepedia.org/dynamic-relaxation)) solvers iterate until the residual drops below tolerance; because the residual may include physical transient terms, the method is not limited to quasi-static problems.[17](https://gmd.copernicus.org/articles/19/5343/2026/gmd-19-5343-2026.html) The computational trade-off is explicit: in poroelastodynamics the CFL condition makes the critical time step the time for a [P wave](https://www.edgepedia.org/p-wave) to cross the smallest cell dimension, whereas implicit poroelastostatics has no such restriction and can take large time steps limited only by accuracy.\n\n## References\n\n1. [Quasistatics, Sierra/SM Theory Manual (Sandia National Laboratories)](https://www.sandia.gov/files/sierra/SM_Theory_5_30/main/quasistatics.html)\n2. [A blended transient/quasistatic Lagrangian framework for salt tectonics simulations (IJNME, 2021)](https://onlinelibrary.wiley.com/doi/10.1002/nme.6671)\n3. [Beyond the Stokes approximation: shallow visco-elastic ice-sheet models (Journal of Glaciology)](https://www.cambridge.org/core/journals/journal-of-glaciology/article/beyond-the-stokes-approximation-shallow-viscoelastic-icesheet-models/9FB2F9C588EB9084DD0731DAD308831C)\n4. [The Continuously Stratified Quasi-Geostrophic System (Vallis, textbook derivation)](https://pordlabs.ucsd.edu/wryoung/theorySeminar/pdf14/VallisQGderivation.pdf)\n5. [Comparing poroelastostatics and poroelastodynamics: Numerics, solvers and algorithms](https://doi.org/10.31224/osf.io/ps2ve)\n6. [Quasidynamic versus fully dynamic simulations of earthquakes and aseismic slip (Thomas et al., JGR 2014)](https://tectonics.caltech.edu/publications/pdf/Thomas-Dynamic_vsQuasidynamic-JGR-2014.pdf)\n7. [Historical review of poroelastodynamics (TU Graz preprint)](https://www.tugraz.at/fileadmin/user_upload/Institute/AM-BM/Files/Forschung/Preprints/preprint_2008_03.pdf)\n8. [Maurice A. Biot (1941). General Theory of Three-Dimensional Consolidation. Journal of Applied Physics.](https://doi.org/10.1063/1.1712886)\n9. [Quasi-static finite element modeling of seismic attenuation due to wave-induced fluid flow (Quintal et al., JGR 2011; author-hosted copy)](https://marcelfrehner.ch/resourcen/Quintal_et_al_JGR_2011.pdf)\n10. [Nadia Lapusta and colleagues (2000). Elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate‐ and state‐dependent friction. Journal of Geophysical Research Atmospheres.](https://doi.org/10.1029/2000jb900250)\n11. [Quasi-static poroelastic parameters in rock and their geophysical applications (Pure and Applied Geophysics)](https://link.springer.com/article/10.1007/BF00998332)\n12. [Review of dynamic earthquake rupture modeling approaches (SRL)](https://par.nsf.gov/servlets/purl/10380873)\n13. [Quasi-Geostrophic Analysis (Univ. of Washington course text, Ch. 6)](https://atmos.uw.edu/academics/classes/2010Q4/441/ch6.pdf)\n14. [EQquasi GitHub repository (EQsimu quasi-dynamic earthquake cycle simulator)](https://github.com/dunyuliu/EQquasi)\n15. [Incorporating Full Elastodynamic Effects into earthquake-cycle simulations (SEAS code comparison)](https://par.nsf.gov/servlets/purl/10417587)\n16. [Spectral-element simulations of long-term fault slip (Kaneko et al., JGR 2011)](https://tectonics.caltech.edu/publications/pdf/Kaneko_JGR2011.pdf)\n17. [Automatic tuning of iterative pseudo-transient solvers for modeling deformation of heterogeneous media (GMD, 2026)](https://gmd.copernicus.org/articles/19/5343/2026/gmd-19-5343-2026.html)\n\n---\n*Topic: Encyclopedia › Physical world and mathematics › Earth sciences › Earth systems and geophysics*\n\n*Initially written Sep 29, 2026 · Reviewed: Sep 30, 2026 · Edited: Sep 30, 2026 · Last review: Sep 30, 2026*\n\n*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*\n\nLicense: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license\n", "same_as": [], "url": "https://www.edgepedia.org/quasistatic-approximation", "markdown_url": "https://www.edgepedia.org/quasistatic-approximation.md", "license": { "name": "Edgepedia Community License 1.0", "url": "https://www.edgepedia.org/edgepedia/license", "summary": "Free with credit, commercial use included. AI training is open to everyone. For other uses, organizations over USD 100M in revenue or 100M monthly users license separately.", "spdx": "LicenseRef-Edgepedia-Community-1.0" }, "credit": "\"Quasistatic approximation\", Edgepedia (EdgeChat), https://www.edgepedia.org/quasistatic-approximation. Edgepedia Community License 1.0.", "credit_md": "\"[Quasistatic approximation](https://www.edgepedia.org/quasistatic-approximation)\", Edgepedia (EdgeChat), [https://www.edgepedia.org/quasistatic-approximation](https://www.edgepedia.org/quasistatic-approximation). [Edgepedia Community License 1.0](https://www.edgepedia.org/edgepedia/license).", "credit_html": "\"Quasistatic approximation\", Edgepedia (EdgeChat), https://www.edgepedia.org/quasistatic-approximation. Edgepedia Community License 1.0.", "speakable": "The quasistatic approximation is a geophysics modeling simplification that neglects inertial terms in the momentum balance, underpinning poroelastic consolidation, crustal deformation, ice-sheet flow, and atmospheric dynamics." }