Mogi: point source of dilatation in elastic half-space

Computes surface displacements, tilt and strain due to point source of dilatation at depth
2.5K Downloads
Updated 15 Sep 2019

View License

The Mogi [1958] model allows to calculate analytic solution for surface deformation due to a point source in an elastic half-space. This model is widely used to simulate ground deformation produced by local perturbation like volcanic magma chamber. It computes displacements, tilt and strain in a polar space, due to a volume variation at depth [Anderson, 1936] or an isotropic pressure variation in a spherical source [Mogi, 1958].

The proposed Matlab script is a literal transcription of Mogi's simple equations, extended to non-isotropic medium (Poisson's ratio different from 0.25). All parameters can be vectorized. See help for syntax, and script comments for details.

Cite As

François Beauducel (2024). Mogi: point source of dilatation in elastic half-space (https://www.mathworks.com/matlabcentral/fileexchange/25943-mogi-point-source-of-dilatation-in-elastic-half-space), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R11.1
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Seismology in Help Center and MATLAB Answers

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes
1.3.0.1

- improves documentation
- comment input arguments checking

1.3.0.0

Bug correction: radial strain value (optional er output argument) was wrong... Now corrected.

1.2.0.0

IMPORTANT: corrects an error when specifying elastic parameters E and nu.
Adds some checks on input arguments.
Adds a scheme explaining the model parameters.

1.1.0.0

Possibility to specify a Poisson's ratio different from 0.25 for non-isotropic medium. Previous input parameters are still allowed.

1.0.0.0