This page is subdivided into the following sections:
The determine if a candidate hypocenter is falsified by observations or compatible with these observations, the Degree of Compatibility Function DOC(x,y,z,t) is defined.
First, we have to ask how to falsify a single 4D point h(x,y,z,t) as hypocenter:
In the following, I will focus on the treatment of arrival times. Directional data are much easier to handle and I think the concept is easy to transfer.
Observation data is usually represented by a single value and an uncertainty measure, inspired by the two parameters needed to define a normal (gaussian) distribution. The mechanism of least squares optimization does indeed imply normal distributions.
Considering arrival times of seismic phases, it is often unclear how to determine, estimate or guess the standard deviation, since phase picking is still a very subjective process. In the case of having only few seismograms of possibly poor quality, the usual data model therefore seems to be ill suited.
Much easier than determining the exact onset time of a seismic phase is to define a time interval which only contains the onset. It is in fact always possible to define a finite length interval which contains, say, the P arrival with 100% certainty. In some cases, this interval would be very long (up to several minutes), but with more than one such interval, it can still be used to constrain the hypocenter.
LOCSMITH therefore defines all data, arrival times as well as angles and other direction data, as intervals: a scalar datum is given by the closed interval [xmin, xmax], such that the true value is between these bounds with 100% certainty. Computer scientists have developed a whole algebra for the so-called interval arithmetic based on this data model.
A strict implementation of interval arithmetic, however, would make the program extremely slow. LOCSMITH only uses the data model, but carries out its computations with normal airthmetics.
The Degree of Compatibility function DOC(x,y,z,t) simply gives the number of observations with which a candidate solution h(x,y,z,t) is compatible. It is therefore an integer valued function of the hypocenter coordinates.
To determine the DOC for arrival times, follow the following list of steps:
Note that the degree of compatibility function DOC(x,y,z,t) is single-valued, integer-valued, positive and non-continuous (stair-like).