Miroslav Hallo - Software
Bayesian earthquake hypocenter location from direct P and S-waves arrival times
- Open source Matlab functions for assessment of the earthquake hypocenter location by grid search method in the homogeneous or 1D layered medium.
The posterior probability density function (PDF) of the earthquake hypocenter location is following theory by Tarantola (2005, Chapter 7.1).
Hence, the solution contains also information about hypocenter location uncertainty.
Codes are published under the GNU General Public License (GNU GPL) for non-commercial use.
Methodology described in:
- Hallo, M., Opršal, I., Asano, K., Gallovič, F. (2019): Seismotectonics of the 2018 Northern Osaka M6.1 earthquake and its aftershocks: joint movements on strike-slip and reverse faults in inland Japan, Earth, Planets and Space, 71:34.
-
List of the programs:
- File crustal.dat contains 1D velocity model of Earth's crust.
- The code a1_Pert_Model.m perturb velocity model and perform Monte-Carlo simulations to get uncertainties of P and S wave arrival times.
- The code a2_LatLon2xyz.m converts coordinates of stations from WGS84 Lat/Lon into local coordinates Easting/Northing
- The code a3_hypo_homo.m compute earthquake hypocenter location by grid search method in the homogeneous medium
- The code a4_hypo_layer.m compute earthquake hypocenter location by grid search method in the 1D layered medium
- The a3_hypo_homo.m is much faster than a4_hypo_layer.m, but less precise as the homogeneous medium is used. We recommend to use a3_hypo_homo.m to rough search for hypocentre and the a4_hypo_layer.m for precise adjustments
Download
MATLAB (revision 12/2018):
-
Requirements
MATLAB:
-
The codes do not require any additional MatLab Toolboxes
The codes were tested on various versions of MatLab software, and it should work without any problems even on older versions. Tested versions: Matlab 2012b 32-bit (Windows 7), Matlab 2018b 64-bit (Windows 7).
