Probabilistic Earthquake Location in 2D
Updated 8 Feb 2017
Demonstrates the effect of observation and velocity model errors on the estimation of epicenter location. Assumed data are single P-phase arrivals at each station. The method is based on forward modelling of hyperbolas and their joint likelihood.
eqX, eqX: true epicenter location (km)
staX, staY: station coordinates (km)
params: parameter object (see description below)
- figure of the used parabolas and the resulting spatial PDF of the epicenter
- error PDFs of the estimated magnitude, horizontal and vertical coordinates
This code was introduced in the following paper to demonstrate the effect of station distribution on magnitude errors.
Kamer & Hiemer (2015), Data-driven spatial b value estimation with applications to California seismicity: To b or not to b, JGR Solid Earth, doi:10.1002/2014JB011510.
1) Event at [10,10] using the default station distribution
2) Same event with 3 stations
probEqLoc2D(10, 10, [10 20 30], [20 15 25]);
Yavor Kamer (2023). Probabilistic Earthquake Location in 2D (https://www.mathworks.com/matlabcentral/fileexchange/61280-probabilistic-earthquake-location-in-2d), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
- Sciences > Earth, Ocean, and Atmospheric Sciences > Geology >
- Sciences > Earth, Ocean, and Atmospheric Sciences > Seismology >
- Engineering > Civil and Environmental Engineering > Earthquake Engineering >
Inspired by: interparc, Generate Random Numbers from a 2D Discrete Distribution
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.