Rank: 446097 based on 0 downloads (last 30 days) and 0 file submitted
photo

Matang

E-mail

Personal Profile:

 

Watch this Author's files

 

Comments and Ratings by Matang
Updated File Comments Rating
17 Aug 2012 Okada: Surface deformation due to a finite rectangular source Computes Okada's 1985 solution for displacements, tilts and strains due to fault dislocation. Author: François Beauducel

Great code!
Try this:

Case of the 2011, Mw 9.0 Tohoku earthquake

Slip (D):
Mo = µ D S
µ = 30 GPa (elastic shear modulus)
S = 600000 m * 200000 m (fault surface)
log10(M0) = 1.5 * Mw + 9.1
Mw = 9.0
obtained: D = 11.0586 m

%----------------------
%Epicenter
LatEpic=38.297;
LonEpic=142.372;
% Observation coordinate at GEONET GPS station (#0175)
LatGPS_sta=38.68269882693;
LonGPS_sta=141.4493698139;

E=(LonGPS_sta-LonEpic)*111*1000;
N=(LatEpic-LatGPS_sta)*111*1000;

[uE,uN,uZ] = okada85(E,N,14000,193,24,600000,200000,81,11.0586,0)
%----------------------

got the the solution:
uE = 4.0607
uN = -1.8354
uZ = -0.4523

The GPS data from the ARIA team at JPL and Caltech (version 2.0) for the station:
uE = 4.0394
uN = -1.612
uZ = -0.6559

--Muzli
GFZ-Potsdam

Contact us