principles and applications of microearthquake networks
principles and applications of microearthquake networks
principles and applications of microearthquake networks
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
% 4. Seismic Ray Tracing for Minimum Time Path<br />
4.4.4. One-Dimensional Multilayer Velocity Models<br />
The most commonly used velocity model in <strong>microearthquake</strong> studies<br />
consists <strong>of</strong> a sequence <strong>of</strong> horizontal layers <strong>of</strong> constant velocity, each layer<br />
having a higher velocity than the layer above it. For the source <strong>and</strong> the<br />
station at the same elevation, the ray paths <strong>and</strong> travel times are well<br />
known (e.g., Officer, 1958, pp. 239-242). For an earthquake source at<br />
depth, the specific formulas involved have been given by Eaton (1969, pp.<br />
2638). We now derive an equivalent set <strong>of</strong> formulas by a systematic<br />
approach using results from Section 4.4.1.<br />
Let us first consider the simple case <strong>of</strong> a single layer <strong>of</strong> thickness h, <strong>and</strong><br />
velocity v1 over a half-space <strong>of</strong> velocity vz (Fig. 20). For a surface source<br />
at point A, we consider two possible ray paths to a surface station at point<br />
B. The travel time <strong>of</strong> the direct wave, Td, along p atha is<br />
(4.88) Td = A/u1<br />
where A is the epicentral distance between the source <strong>and</strong> the station, i.e.,<br />
A = AT. The travel time <strong>of</strong> the refracted wave, T,, along pathACDB is<br />
(4.89) T, = (AT/vl) + (m/u,) + (m/v,)<br />
Using Snell's law, we have<br />
(4.90) sin 8/v, = sin 90°/u2 = sin 4/v1<br />
so that 8 = qb (Fig. 20), <strong>and</strong> AT = m. Consequently, Eq. (4.89) can be<br />
written in terms <strong>of</strong> 8, hl, <strong>and</strong> A as (see Officer, 1958, p. 240),<br />
(4.91)<br />
A 2h,<br />
T, = - + - cos 8<br />
v2 1' 1<br />
Using Eqs. (4.88) <strong>and</strong> (4.91), we can plot travel time versus epicentral<br />
A<br />
I<br />
A(source)<br />
-<br />
B(station)<br />
h, Velocity=V,<br />
-<br />
r<br />
A' C D B'<br />
Fig. 20.<br />
half-space.<br />
Velocity = V2<br />
Diagram <strong>of</strong> direct <strong>and</strong> refracted paths in a velocity model <strong>of</strong> one layer over a