Okada: Surface deformation due to a finite rectangular source

Computes Okada's 1985 solution for displacements, tilts and strains due to fault dislocation.
5.6K 次下载
更新时间 2014/5/28

查看许可证

The Okada [1985] model calculates analytical solution for surface deformation due to shear and tensile faults in an elastic half-space. This model is widely used to simulate ground deformation produced by local perturbation like tectonic faults (earthquakes) or volcanic dykes (magmatic intrusion). Given rectangular fault geometry (length, width, depth, strike, dip) and 3-component dislocation amplitude (rake, slip and open), it computes the displacements, tilts and strains at the free-surface.
The proposed Matlab script is a literal transcription of the Okada's equations, except that it is transposed in a geographical referential (East, North, Up), where the fault is defined by a strike angle relative to the North, and dislocation parameters are given by: rake, slip and opening (instead of U1, U2, U3), following Aki & Richards [1980] definition. All coordinates and depth are relative to fault centroid. Lamé's constants λ and μ are replaced by Poisson's ratio ν (with a default value of 0.25 for isotropic medium), since the equations are independent of other elastic parameters. The equations are also vectorized for (x,y) coordinates and all input parameters except dip angle.

To check the consistency of numerical calculations, run the script okada85_checklist.m, a transcription of table 2 cases 2, 3, and 4 checklist from [Okada, 1985] paper (needs also the roundsd.m function).

See help for further details, syntax, example, and script comments for technical details.

引用格式

François Beauducel (2024). Okada: Surface deformation due to a finite rectangular source (https://www.mathworks.com/matlabcentral/fileexchange/25982-okada-surface-deformation-due-to-a-finite-rectangular-source), MATLAB Central File Exchange. 检索来源 .

MATLAB 版本兼容性
创建方式 R13
兼容任何版本
平台兼容性
Windows macOS Linux

Community Treasure Hunt

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

Start Hunting!
版本 已发布 发行说明
1.7.0.0

fixes a problem with K1 function (tilt) when DIP=90. Thanks to Halldor Geirsson.

1.6.0.0

- allows vectorization of RAKE, SLIP and OPEN
- checklist function to validate numerical consistancy

1.5.0.0

Corrects another bug in the plot figure.
Minor modifications in the code and help text.

1.4.0.0

Bug correction on fault centroid exact position. Add optional plot and example.

1.3.0.0

Coordinates are now relative to fault centroid (instead of middle top edge). This is more convenient for earthquake hypocenter use.

1.2.0.0

ATTENTION: correction of 3 errors in some equations (functions I1, K2 and uyy_tf) affecting some component values. Now the function fits all the numerical values of Okada's checklist cases 2, 3 and 4. Thanks to Dmitry Nikolski for his contribution.

1.1.0.0

Corrects some output arguments description in the help (inversion of uN and uE), and better explanation of input arguments.
Adds a 3-D graphical output the of fault geometry when using without output argument.

1.0.0.0