Forward modelling the multiply scattered 2.5-D teleseismic P waves accelerated by the multilevel fast multipole method

Geophysical Journal International

Volumn 176, Issue 2, 2009, Pages 505-517

  Çakir, Özcan ^a  

^a NEVSEHIR HACI BEKTAS VELI UNIVERSITY   (Turkey)

Author keywords

Body waves; Composition of the mantle; Computational seismology; Wave propagation; Wave scattering and diffraction

Indexed keywords

ALGORITHM; BODY WAVE; FORWARD MODELING; GREEN FUNCTION; HETEROGENEITY; LITHOSPHERIC STRUCTURE; P-WAVE; REFLECTIVITY; SEISMIC REFLECTION; TELESEISMIC WAVE; THREE-DIMENSIONAL MODELING; TWO-DIMENSIONAL MODELING; WAVE DIFFRACTION; WAVE PROPAGATION; WAVE SCATTERING;

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

References (83)

1. Aki, K. & Richards, P.G., 1980. Quantitative Seismology: Theory and Methods, ed. Freeman, W.H., San Francisco, CA.
2. Ammon, C.J., Randall, G.E. & Zandt, G., 1990. On the nonuniqueness of receiver function inversions, J. geophys. Res., 95, 15 303-15 318.
3. Angus, D.A., Wilson, D.C, Sandvol, E. & Ni, J.F., 2006. Lithospheric structure of the Arabian and Eurasian collision zone in eastern Turkey from S -wave receiver functions, Geophys. J. Int., 166, 1335-1346.
4. Assumpção, M., James, D. & Snoke, A., 2002. Crustal thicknesses in SE Brazilian Shield by receiver function analysis: Implications for isostatic compensation, J. geophys. Res., 107, doi:.
5. Berteussen, K.A., 1977. Moho depth determinations based on spectral analysis of NORSAR long period P waves, Phys. Earth. planet. Inter., 134, 13-27.
6. Bohlen, T., 2002. Parallel 3-D viscoelastic finite difference seismic modelling, Comp. Geosci., 28, 887-899.
7. Bostock, M.G., 2002. Kirchhoff-approximate inversion of teleseismic wavefields, Geophys. J. Int., 149, 787-795.
8. Bostock, M.G. & Rondenay, S., 1999. Migration of scattered teleseismic body waves, Geophys. J. Int., 137, 732-746.
9. Bostock, M.G., Rondenay, S. & Shragge, J., 2001. Multiparameter two-dimensional inversion of scattered teleseismic body waves, 1: Theory for oblique incidence, J. geophys. Res., 106, 30771-30782.
10. Bouchon, M., 2003. A review of the discrete wavenumber method, Pure appl. Geophys., 160, 445-465.
11. Bunse-Gerstner, A. & Gutiérrez-Cañas, I., 2006. A preconditioned GMRES for complex dense linear systems from electromagnetic wave scattering problems, Linear Algebra Appl., 416, 135-147.
12. Chen, L., Wen, L., Zheng, T., 2005. A wave equation migration method for receiver function imaging, 2: Application to the Japan subduction zone, J. geophys. Res., 110, B11310, doi:.
13. Cummins, P.R., Takeuchi, N. & Geller, R.J., 1997. Computation of complete synthetic seismograms for laterally heterogeneous models using the Direct Solution Method, Geophys. J. Int., 130, 1-16.
14. Çakir, Ö., 2006. The multilevel fast multipole method for forward modelling the multiply scattered seismic surface waves, Geophys. J. Int., 167, 663-678.
15. Çakir, Ö. & Erduran, M., 2004. Constraining crustal and uppermost mantle structure beneath station TBZ (Trabzon, Turkey) by receiver function and dispersion analyses, Geophys. J. Int., 158, 955-971.
16. Du, Z.J. & Foulger, G.R., 1999. The crustal structure beneath the northwest fjords, Iceland, from receiver functions and surface waves, Geophys. J. Int., 139, 419-432.
17. Dueker, K. & Yuan, H., 2004. Upper mantle P-wave velocity structure from PASSCAL teleseismic transects across Idaho, Wyoming and Colorado, Geophys. Res. Lett., 31, L08603, doi:.
18. Erduran, M. & Çakir, Ö., 2001. Non-uniqueness in the determination of crustal and upper-mantle velocity structure from receiver function (in Turkish), Yerbilimleri, 23, 99-112.
19. Faccioli, E., Maggio, F., Paolucci, R. & Quarteroni, A., 1997. 2D and 3D elastic wave propagation by a pseudo-spectral domain decomposition method, J. Seismol., 1, 237-251.
20. Foss, S.K., de Hoop, M.V. & Ursin, B., 2005. Linearized 2.5-dimensional parameter imaging inversion in anisotropic elastic media, Geophys. J. Int., 161, 722-738.
21. Frederiksen, A.W. & Revenaugh, J., 2004. Lithospheric imaging via teleseismic scattering tomography, Geophys. J. Int., 159, 978-990.
22. Fu, L.Y. & Bouchon, M., 2004. Discrete wavenumber solutions to numerical wave propagation in piecewise heterogeneous media, I: Theory of two-dimensional SH case, Geophys. J. Int., 157, 481-498.
23. Fujiwara, H., 1998. The fast multipole method for integral equations of seismic scattering problems, Geophys. J. Int., 133, 773-782.
24. Fujiwara, H., 2000. The fast multipole method for solving integral equations of three-dimensional topography and basin problems, Geophys. J. Int., 140, 198-210.
25. Gao, H. & Zhang, J., 2006. Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: A grid method approach, Geophys. J. Int., 165, 875-888.
26. Gilbert, H. J., Sheehan, A. F., Dueker, K. G. & Molnar, P., 2003. Receiver functions in the western United States, with implications for upper mantle structure and dynamics, J. Geophys. Res., 108, 2229, doi:.
27. Gök, R., Pasyanos, M. E. & Zor, E., 2006. Lithospheric structure of the continent-continent collision zone: Eastern Turkey, Geophys. J. Int., 169, 1079-1088.
28. Greengard, L. & Rokhlin, V., 1987. A Fast Algorithm for Particle Simulations, J. Comp. Phys., 73, 280-292.
29. Gumerov, N.A. & Duraiswami, R., 2004. Computation of scattering from clusters of spheres using the fast multipole method, J. Acoust. Soc. Am., 117, 1744-1761.
30. Gumerov, N.A., Duraiswami, R. & Borovikov, E.A., 2003. Data structures, optimal choice of parameters, and complexity results for generalized multilevel fast multipole methods in d dimensions Computer Science Technical Report No. 4458, Vol. 91, University of Maryland, College Park, Maryland.
31. Haines, A.J., Hulme, T. & Yu, J., 2004. General elastic wave scattering problems using an impedance operator approach, I: mathematical development, Geophys. J. Int., 159, 643-657.
32. Hara, T., 2004. Waveform inversion for 3-D earth structure using the direct solution method implemented on vector-parallel supercomputer, Phys. Earth. planet. Int., 146, 65-74.
33. Hertweck, T., Jäger, C., Goertz, A. & Schleicher, J., 2003. Aperture effects in 2.5D Kirchhoff migration: A geometrical explanation, Geophysics, 68, 1673-1684.
34. Hestholm, S., 1999. Three-dimensional finite difference viscoelastic wave modelling including surface topography, Geophys. J. Int., 139, 852-878.
35. Honda, R. & Yomogida, K., 2003. Contribution of vertically travelling plane S-waves to dynamic and static displacements near a finite fault, Geophys. J. Int., 152, 443-454.
36. Huang, L.J. & Fehler, M.C., 2000. Quasi-Born Fourier migration, Geophys. J. Int., 140, 521-534.
37. Iturrarán-Viveros, U., Sánchez-Sesma, F.J. & Janod, F., 2004. Scattering of elastic waves in heterogeneous media using the direct solution method, Geophys. J. Int., 156, 222-236.
38. Julià, J., Ammon, C.J., Herrmann, R.B. & Correig, A.M., 2000. Joint inversion of receiver function and surface wave dispersion observations, Geophys. J. Int., 143, 99-112.
39. Kawai, K., Takeuchi, N. & Geller, R.J., 2006. Complete synthetic seismograms up to 2 Hz for transversely isotropic spherically symmetric media, Geophys. J. Int., 164, 411-424.
40. Kenneth, B.L.N., 1983. Seismic Wave Propagation in Stratified Media, Cambridge University Press, New York.
41. Komatitsch, D. & Tromp, J., 1999. Introduction to the spectral element method for three-dimensional seismic wave propagation, Geophys. J. Int., 139, 806-822.
42. Komatitsch, D. & Tromp, J., 2002. Spectral-element simulations of global seismic wave propagation, I: Validation, Geophys. J. Int., 149, 390-412.
43. Langston, C.A., 1979. Structure under Mountain Rainer, Washington, inferred from teleseismic body waves, J. geophys. Res., 84, 4749-4762.
44. Lawrence, J.F. & Wiens, D.A., 2004. Combined receiver-function and surface wave phase-velocity inversion using a niching genetic algorithm: application to Patagonia, Bull. seism. Soc. Am., 94, 977-987.
45. Levander, A., 2003. USArray design implications for wavefield imaging in the lithosphere and upper mantle, Leading Edge, 22, 250-255.
46. Levander, A., Niu, F., Lee, C.A. & Cheng, X., 2006. Imag(in)ing the continental lithosphere, Tectonophysics, 416, 167-185.
47. Li, X., Kind, R. & Yuan, X., 2003. Seismic study of upper mantle and transition zone beneath hotspots, Phys. Earth planet. Inter., 136, 79-92.
48. Li, A., Fischer, K.M., Van Der Lee, S. & Wysession, M.E., 2002. Crust and upper mantle discontinuity structure beneath eastern North America, J. geophys. Res., 107, 2100,.
49. Lysmer, J. & Drake, L.A., 1972. A finite element method for seismology, in Methods of Computational Physics, Vol. 11, pp. 181-216, Academic Press, New York.
50. Maupin, V., 2001. A multiple-scattering scheme for modelling surface wave propagation in isotropic and anisotropic three-dimensional structures, Geophys. J. Int., 146, 332-348.
51. Moczo, P., Kristek, J. & Halada, L., 2000. 3D Fourth-Order Staggered-Grid Finite-Difference Schemes: Stability and Grid Dispersion, Bull. seism. Soc. Am., 90, 587-603.
52. Müller, G., 1969. Theoretical seismograms for some types of point-sources in layered media. Part III: Single force and dipole sources of arbitrary orientation, Z. Geophys., 35, 347-371.
53. Müller, G., 1985. The reflectivity method: A tutorial, J. Geophys., 58, 153-174.
54. Nair, S.K., Gao, S.S., Liu, K.H. & Silver, P.G., 2006. Southern African crustal evolution and composition: Constraints from receiver function studies, J. geophys. Res., 111, B02304, doi:.
55. Nihei, K.T. & Li, X., 2007. Frequency response modelling of seismic waves using finite difference time domain with phase sensitive detection (TD-PSD), Geophys. J. Int., 169, 1069-1078.
56. Nishimura, N., 2002. Fast multipole accelerated boundary integral equation methods, Appl. Mech. Rev., 55, 299-324.
57. Niu, F., Levander, A., Ham, S. & Obayashi, M., 2005. Mapping the subducting Pacific slab beneath southwest Japan with Hi-net receiver functions, Earth planet. Sci. Lett, 239, 9-17.
58. Novais, A. & Santos, L.T., 2005. 2.5D finite-difference solution of the acoustic wave equation, Geophys. Prospect., 53, 523-531.
59. Nowack, R.L., Dasgupta, S., Schuster, G.T. & Sheng, J.M., 2006. Correlation migration using Gaussian beams of scattered teleseismic body waves, Bull. seism. Soc. Am., 96, 1-10.
60. Özalaybey, S., Savage, M.K., Sheehan, A.F., Louie, J.N. & Brune, J.N., 1997. Shear-wave velocity structure in the northern Basin and Range province from the combined analysis of receiver functions and surface waves, Bull. seism. Soc. Am., 87, 183-199.
61. Pain, C.C., Saunders, J.H., Worthington, M.H., Singer, J.M., Stuart-Bruges, W., Mason, G. & Goddard, A., 2005. A mixed finite-element method for solving the poroelastic Biot equations with electrokinetic coupling, Geophys. J. Int., 160, 592-608.
62. Papageorgiou, A.S. & Pei, D., 1998. A Discrete Wavenumber boundary element method for study of the 3-D response 2-D scatterers, Earthq. Eng. Struct. Dyn., 27, 619-638.
63. Pollitz, F.F., 2002. Regional seismic wavefield computation on a 3-D heterogeneous earth model by means of coupled travelling wave synthesis, Pure appl. Geophys., 159, 2085-2112.
64. Rokhlin, V., 1985. Rapid solution of integral equations of classical potential theory, J. Comp. Phys., 60, 187-207.
65. Rondenay, S., Bostock, M.G. & Shragge, J., 2001. Multiparameter two-dimensional inversion of scattered teleseismic body waves, 3: application to the Cascadia 1993 data set, J. geophys. Res., 106, 30 795-30 807.
66. Rost, S. & Thomas, C., 2002. Array seismology: Methods and applications, Rev. Geophys., 40, 1008, doi:.
67. Ryberg, T & Weber, M., 2000. Receiver function arrays: A reflection seismic approach, Geophys. J. Int., 141, 1-11.
68. Saad, Y., 2003. Iterative Methods for Sparse Linear Systems, 2nd edn, Society for Industrial and Applied Mathematics, Philadelphia.
69. Schwab, F.A. & Knopoff, L., 1972. Fast surface wave and free mode computations, in Methods of Computational Physics, Vol. 11, pp. 87-180, Academic Press, New York.
70. Shragge, J., Bostock, M.G. & Rondenay, S., 2001. Multiparameter two-dimensional inversion of scattered teleseismic body waves, 2: numerical examples, J. geophys. Res., 106, 30783-30793.
71. Snieder, R., 2002. General theory of elastic wave scattering, Scattering and Inverse Scattering in Pure and Applied Science, pp. 528-542, eds Pike, R. & Sabatier, P., Academic Press, San Diego.
72. Sodoudi, F. et al., 2006. Lithospheric structure of the Aegean obtained from P and S receiver functions, J. geophys. Res., 111, B12307, doi:.
73. Tadeu, A., Mendes, P.A. & António, J., 2006. The simulation of 3D elastic scattering produced by thin rigid inclusions using the traction boundary element method, Comp. Struct., 84, 2244-2253.
74. Vai, R., Castillo-Covarrubias, J.M., Sánchez-Sesma, F.J., Komatitsch, D. & Vilotte, J.P., 1999. Elastic wave propagation in an irregularly layered medium, Soil Dyn. Earthq. Eng., 18, 11-18.
75. van der Meijde, M., van der Lee, S. & Giardini, D., 2003. Crustal structure beneath broad-band seismic stations in the Mediterranean region, Geophys. J. Int., 152, 729-739.
76. Vinnik, L. P., 1977. Detection of waves converted from P to SV in the mantle, Phys. Earth planet. Inter., 15, 39-45.
77. Vinnik, L.P., Reigber, C., Aleshin, I.M., Kosarev, G.L., Kaban, M.K., Oreshin, S.I. & Roecker, S.W., 2004. Receiver function tomography of the central Tien Shan, Earth Planet. Sci. Lett., 225, 131-146.
78. Vuik, K., Sevink, A.G.J. & Herman, G.C., 1996. A preconditioned Krylov subspace method for the solution of least squares problems in inverse scattering, J. Comp. Phys., 123, 330-340.
79. Wilson, D., Aster, R., Ni, J., Grand, S., West, W., Gao, W., Baldridge, W.S. & Semken, S., 2005. Imaging the seismic structure of the crust and upper mantle beneath the Great Plains, Rio Grande Rift, and Colorado Plateau using receiver functions, J. geophys. Res., 110, B05306, doi:.
80. Wu, R.S., Jin, S. & Xie, X.B., 2000. Seismic wave propagation and scattering in heterogeneous crustal waveguides using screen propagators, I: SH waves, Bull. seism. Soc. Am., 90, 401-413.
81. Xie, X.B. & Wu, R.S., 2001. Modeling elastic wave forward propagation and reflection using the complex screen method, J. Acoust. Soc. Am., 109, 2629-2635.
82. Yamauchi, M., Hirahara, K. & Shibutani, T., 2003. High resolution receiver function imaging of the seismic velocity discontinuities in the crust and the uppermost mantle beneath southwest Japan, Earth Planet. Space, 55, 59-64.
83. Zhu, L., 2000. Crustal structure across the San Andreas Fault, southern California from teleseismic converted waves, Earth planet. Sci. Lett., 179, 183-190.