A finite-difference algorithm for full waveform teleseismic tomography

Geophysical Journal International

Volumn 181, Issue 2, 2010, Pages 1017-1040

  Roecker, S ^a   Baker, B ^a   McLaughlin, J ^a  

^a RENSSELAER POLYTECHNIC INSTITUTE   (United States)

Author keywords

Body waves; Computational seismology; Seismic tomography; Wave scattering and diffraction

Indexed keywords

ALGORITHM; BODY WAVE; DATA INVERSION; FINITE DIFFERENCE METHOD; SEISMIC DATA; SEISMIC MIGRATION; SEISMIC SOURCE; SEISMIC TOMOGRAPHY; SPECTRAL ANALYSIS; TELESEISMIC WAVE; WAVE DIFFRACTION; WAVE FIELD; WAVE PROPAGATION; WAVE SCATTERING; WAVEFORM ANALYSIS;

EID: 77953103150     PISSN: 0956540X     EISSN: 1365246X     Source Type: Journal    
DOI: 10.1111/j.1365-246X.2010.04553.x     Document Type: Article

References (58)

1. Akçelik V, Biros G, Ghattas O. Parallel multiscale Gauss-Newton-Krylov methods for inverse wave propagation. 2002, 10.1109/SC.2002.10002, http://www.sc-conference.org/sc2002, doi, published on CD-ROM and at
2. Akçelik V, Biros G, Ghattas O. High resolution forward and inverse earthquake modeling on terascale computers. 2003, 10.1109/SC.2003.10056, http://www.sc-conference.org/sc2003, doi, published on CD-ROM and at
3. Aki K, Christoffersson A, Husebye ES. Determination of the three-dimensional seismic structure of the lithosphere. J. geophys. Res. 1977, 82:277-296.
4. Bleibinhaus F, Hole JA, Ryberg T, Fuis GS. Structure of the California Coast Ranges and San Andreas Fault at SAFOD from seismic waveform inversion and reflection imaging. J. geophys. Res. 2007, 112:B06315. 10.1029/2006JB004611, doi
5. Bostock MG, Rondenay S, Shragge J. Multiparameter 2-D inversion of scattered teleseismic body-waves. I. Theory for oblique incidence. J. geophys. Res. 2001, 160:30 771-30 782.
6. Brenders A, Pratt RG. Full waveform tomography for lithospheric imaging; results from a blind test in a realistic crustal model. Geophys. J. Int. 2007, 168:133-151.
7. Brittan J, Forgues E, Pratt R, Morgan J, Warner M, Macintyre H, Marin L. Wavefield inversion across the edge of the Chicxulub impact structure. J. Conf. Abstr. 1997, 1:65-72.
8. Chen P, Zhao L, Jordan TH. Full 3D tomography for the crustal structure of the Los Angeles region. Bull. seism. Soc. Am. 2007, 97(4):1094-1120.
9. Collino F, Tsogka C. Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics 2001, 66:294-307.
10. Dahlen FA. Finite-frequency sensitivity kernels for boundary topography perturbations. Geophys. J. Int. 2005, 162:525540.
11. George A, Lui JWH. Computer Solutions of Large Positive Definite Systems 1981, Prentice-Hall Inc., Englewood Cliffs, NJ
12. Haskell NA. The dispersion of surface waves on multilayered media. Bull. seism. Soc. Am. 1953, 43(1):17-34.
13. Haskell NA. Crustal reflection of plane SH waves. J. geophys. Res. 1960, 65(12):4147-4150.
14. Haskell NA. Crustal reflection of plane P and SV waves. J. geophys. Res. 1962, 67(12):4751-4767.
15. Herrmann RB, Ammon CJ, Julia J. Application of joint receiver-function surface-wave dispersion for local structure in Eurasia. Proceedings of the 23rd Seismic Research Review; Worldwide Monitoring of Nuclear Explosions, I 2001, 46-54. in
16. Hole JA, Zelt CA, Pratt RG. Advances in controlled-source seismic imaging. EOS, Trans. Am. geophys. Un. 2005, 86:177-181.
17. Jo CH, Shin CS, Suh JH. An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator. Geophysics 1996, 61(2):529-537.
18. Karypis G, Kumar V. METIS: A Software Package for Partitioning Unstructured Graphs, Partitioning Meshes, and Computing Fill-Reducing Orderings of Sparse Matrices, Version 4.0 1998, http://www.cs.umn.edu/karypis, (Manual), University of Minnesota, Department of Computer Science/Army HPC Research Center
19. Komatitsch D, Tromp J. A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation. Geophys. J. Int. 2003, 154(1):146-153. 10.1046/j.1365-246X.2003.01950.x, doi
20. Kosarev GL, Petersen NV, Vinnik LP, Roecker SW. Receiver functions for the Tien Shan analog broadband network; contrasts in the evolution of structure across the Talasso-Fergano Fault. J. geophys. Res. 1993, 98:4437-4448.
21. Kumar P. The lithosphere-asthenosphere boundary in the North-West Atlantic region. Earth planet. Sci. Lett. 2005, 236:249-257.
22. Küpper FJ. Theoretische Untersuchung über die Mehrfachaufstellung von Geophonen. Geophys. Prospect. 1958, 10.1111/j.1365-2478.1958.tb01652.x, doi
23. Lailly P. The seismic inverse problem as a sequence of before stack migrations. Conference on Inverse Scattering: Theory and Application 1983, 206-220. Bednar J BRedner RRobinson EWeglein A. pp, eds, SIAM, Philadelphia
24. Langston CA. Structure under Mount Rainier, Washington, inferred from teleseismic body waves. J. geophys. Res. 1979, 84:4749-4762.
25. Levander A, Niu F, Symes WW. Imaging teleseismic P to S scattered waves using the Kirchhoff integral. Geophysical Monograph 157 2005, 149-169. American Geophysical Union, Washington, DC, United States
26. Levin V, Park J. P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation. Geophys. J. Int. 1997, 131:253-266.
27. Li Z, Roecker S, Li Z, Wei B, Wang H, Schelochkov G, Bragin V. Tomographic Image of the crust and upper mantle beneath the western Tien Shan from the MANAS broadband deployment: Possible evidence for lithospheric delamination. Tectonophysics 2009, 477(1):49-57.
28. Maggi A, Tape C, Chen M, Chao D, Tromp J. An automated time-window selection algorithm for seismic tomography. Geophys. J. Int. 2009, 178:257-281.
29. Marfurt KJ. Accuracy of finite difference and finite element modeling of the scalar and elastic wave equations. Geophysics 1984, 49:533-549.
30. Mathews JH, Fink KD. Numerical Methods: Using Matlab 2004, 4th edn, Prentice-Hall, Upper Saddle River, NJ
31. Mora PR. Nonlinear two-dimensional elastic inversion of multi-offset seismic data. Geophysics 1987, 52:1211-1228.
32. Oreshin S, Vinnik L, Peregoudov D, Roecker S. Lithosphere and asthenosphere of the Tien Shan imaged by S receiver functions. Geophys. Res. Lett. 2002, 29(8). 1029/2001GLO14441, doi
33. Poppeliers C, Pavlis GL. Three-dimensional, prestack, plane wave migration of teleseismic P-to-S converted phases. 1: theory. J. geophys. Res. 2003a, 108(B2). 10.1029/2001JB000216, doi
34. Poppeliers C, Pavlis GL. Three-dimensional, prestack, plane wave migration of teleseismic P-to-S converted phases; 2, Stacking multiple events. J. geophys. Res. 2003b, 108(B5). 10.1029/2001JB001583, doi
35. Pratt RG. Frequency domain elastic wave modeling by finite differences: a tool for cross-hole seismic imaging. Geophysics 1990, 55:626-632.
36. Pratt RG. Seismic waveform inversion in the frequency domain. Part 1: theory and vertification in a physical scale model. Geophysics 1999, 64:888-901.
37. Pratt RG. Waveform tomography: theory and practice. 2003, http://www.geophys.geos.vt.edu/hole/ccss/prattCCSS.pdf, available at
38. Pratt RG, Worthington MH. Inverse theory applied to multi-source crosshole tomography. Part I: acoustic wave-equation method. Geophys. Prospect. 1990, 38:287-310.
39. Pratt RG, Shin C, Hicks G. Gauss-newton and full newton methods in frequency-space seismic waveform inversion. Geophys. J. Int. 1998, 133:341-362.
40. Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in FORTRAN: The Art of Scientific Computing 1992, 994. 2nd edn, Cambridge Univ. Press, Cambridge, pp
41. Roecker SW, Sabitova TM, Vinnik LP, Burmakov YA, Golvanov MI, Mamatkanova R, Munirova L. Three-dimensional elastic wave velocity structure of the western and central Tien Shan. J. geophys. Res. 1993, 98:15 779-15 795.
42. Ruud B, Hestholm S. 2D surface topography boundary conditions in seismic wave modelling. Geophys. Prospect. 2001, 49:445-460.
43. Schneider JB. Plane Wave in FDTD Simulations and a nearly perfect total-field/scattered-field boundary. IEEE Trans. Antennas Propagat. 2004, 48:494-501.
44. Song ZM, Williamson PR. Frequency-domain acoustic wave modeling and inversion of cross-hole data. Part I: 2.5-D modeling method. Geophysics 1995a, 60:784-795.
45. Song ZM, Williamson PR, Pratt RG. Frequency-domain acoustic wave modeling and inversion of cross-hole data. Part II: inversion method, synthetic experiments and real data results. Geophysics 1995b, 60:796-809.
46. Štekl I, Pratt RG. Accurate visco-elastic modeling by frequency-domain finite differences using rotated operators. Geophysics 1998, 63:1179-1794.
47. Takenaka H, Kennett BLN. A 2.5-D time-domain elastodynamic equation for plane-wave incidence. Geophys. J. Int. 1996, 125:F5-F9.
48. Tan T, Potter M. Perfect Plane wave source for total-field/scattered-field formulation in FDTD using time domain multipoint 1D auxiliary grid. Antennas and Propagation Society International Symposium 2006, IEEE 2006, 3777-3780.
49. Tape C, Liu Q, Maggi A, Tromp J. Adjoint tomography of the southern California crust. Science 2009, 325:988992.
50. Tarantola A. Inversion of seismic reflection data in the acoustic approximation. Geophysics 1984, 49:1259-1266.
51. Vinnik L. Detection of waves converted from P to SV in the mantle. Phys. Earth planet. Inter. 1977, 15:39-45.
52. Vinnik LP, Peregoudov D, Makeyeva L, Oreshin S, Roecker S. Towards 3-D fabric in the continental lithosphere and asthenosphere: the Tien Shan. Geophys. Res. Lett. 2002a, 10.1029/2001GL014588
53. Vinnik LP, Roecker SW, Kosarev GL, Oreshin SI, Koulakov IYu. Crustal structure and dynamics of the Tien Shan. Geophys. Res. Lett. 2002b, 29:41-44.
54. Vinnik LP, Reigber Ch, Aleshin IM, Kosarev GL, Kaban MK, Oreshin SI, Roecker S. Receiver function tomography of the central Tien Shan. Earth planet. Sci. Lett. 2004, 225:131-146.
55. Williamson PW, Pratt RG. A critical review of acoustic wave modelling procedures in 2.5 D. Geophysics 1995, 60:591-595.
56. Wilson D, Aster R. Seismic structure of the lithosphere in the Southwestern United States using receiver functions. Leading Edge 2003, 22:234-237.
57. Woodward MJ. Wave-equation tomography. Geophysics 1992, 57:15-26.
58. Zheng Y, Huang X. Anisotropic perfectly matched layers for elastic waves in Cartesian and curvilinear coordinates. Earth Resources Laboratory 2002 Industry Consortium Meeting 2002, in, Dept. of Earth, Atmospheric, and Planetary sciences, Massachusetts Institute of Technology, Cambridge, MA, USA