Bulletin of the Seismological Society of America
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Bulletin of the Seismological Society of America; August 1980; v. 70; no. 4; p. 1137-1148
© 1980 Seismological Society of America
This Article
Right arrow Full Text (PDF)
Right arrow References
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by THURBER, C. H.
Right arrow Articles by ELLSWORTH, W. L.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation

Rapid solution of ray tracing problems in heterogeneous media

C. H. THURBER* and W. L. ELLSWORTH

U.S. GEOLOGICAL SURVEY, 345 MIDDLEFIELD ROAD, MENLO PARK, CALIFORNIA 94025

Abstract

The determination of local earthquake hypocenters and orgin times from first-P-arrival times by Geiger's method requires a technique for finding the minimum travel time (and derivatives) between the source and the station. Sophisticated ray tracing techniques have been developed for this purpose for use in complex velocity structures. Unfortunately, the two common techniques, shooting and bending, are generally prohibitively expensive for routine use in data analysis. The bending method is also particularly vulnerable to the problem of local minima in travel time. A method has been developed known as the ray initializer, which can be used to circumvent these problems in many cases. First, the technique can find a reasonable estimate of the minimum-time ray path in a quick and efficient manner. The velocity in a region local to the source and receiver is laterally averaged to yield an approximate layered velocity model. One-dimensional ray tracing techniques are used to find the minimum-time path for this layered structure. The ray path estimate can then be used as the starting path in a bending routine, a procedure resulting in more rapid convergence and the avoidance of local minima. Second, the travel time found by numerical integration along the estimated ray path is an excellent approximation to the actual travel time. Thus, in many cases, the ray initializer can be substituted for a three-dimensional ray tracing routine with a tremendous increase in efficiency and only a small loss in accuracy. It is found that the location of an explosion, derived using the ray initializer, is nearly identical to a complete ray tracing solution, even for a highly complex velocity structure.

Footnotes

* Present address: Department of Earth and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139.




This article has been cited by other articles:


Home page
Bulletin of the Seismological Society of AmericaHome page
W. Kim and C.-E. Baag
Rapid and Accurate Two-Point Ray Tracing Based on a Quadratic Equation of Takeoff Angle in Layered Media with Constant or Linearly Varying Velocity Functions
Bulletin of the Seismological Society of America, August 1, 2002; 92(6): 2251 - 2263.
[Abstract] [Full Text] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
H. MAGISTRALE and T. Rockwell
The central and southern Elsinore fault zone, southern California
Bulletin of the Seismological Society of America, December 1, 1996; 86(6): 1793 - 1803.
[Abstract] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
T. J. MOSER, G. NOLET, and R. SNIEDER
Ray bending revisited
Bulletin of the Seismological Society of America, February 1, 1992; 82(1): 259 - 288.
[Abstract] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
H. J. AL-SHUKRI and B. J. MITCHELL
Three-dimensional attenuation structure in and around the New Madrid seismic zone
Bulletin of the Seismological Society of America, June 1, 1990; 80(3): 615 - 632.
[Abstract] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
W. A. PROTHERO, W. J. TAYLOR, and J. A. EICKEMEYER
A fast, two-point, three-dimensional raytracing algorithm using a simple step search method
Bulletin of the Seismological Society of America, June 1, 1988; 78(3): 1190 - 1198.
[Abstract] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
A. J. MICHAEL
Effects of three-dimensional velocity structure on the seismicity of the 1984 Morgan Hill, California, aftershock sequence
Bulletin of the Seismological Society of America, June 1, 1988; 78(3): 1199 - 1221.
[Abstract] [PDF]


Home page
ScienceHome page
C. H. Thurber and C. H. THURBER
Seismic Detection of the Summit Magma Complex of Kilauea Volcano, Hawaii
Science, January 13, 1984; 223(4632): 165 - 167.
[Abstract] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
S. R. TAYLOR and J. F. SCHEIMER
P-velocity models and earthquake locations in the Livermore Valley Region, California
Bulletin of the Seismological Society of America, August 1, 1982; 72(4): 1255 - 1275.
[Abstract] [PDF]




HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 1980 by the Seismological Society of America.