Probabilistic Earthquake Location in 2D

Demonstrates the effect of data and model errors on the estimation of epicenter location.

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.
INPUTS
eqX, eqX: true epicenter location (km)
staX, staY: station coordinates (km)
params: parameter object (see description below)
OUTPUT
- 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.

Examples:

1) Event at [10,10] using the default station distribution
probEqLoc2D(10, 10);

2) Same event with 3 stations
probEqLoc2D(10, 10, [10 20 30], [20 15 25]);