Seismic interferometry of teleseicmic S‐wave coda for retrieval of body waves: an application to the Philippine Sea slab underneath the Japanese Islands

SUMMARY The reconstruction of surface waves from the cross-correlation of random wavefields has recently been extensively inspected by theoretical and experimental approaches. However, the reconstruction of body waves has not been extensively studied. In this study, we present a method for extracting body waves, that is, direct P and S waves, and reflected waves from the Philippine Sea slab. We use the cross-correlation of the wavefield generated by the teleseismic S waves observed by Hi-net tiltmeters with a passband of 0.07–0.5 Hz. To enhance the contribution of the S coda, we reduced the source-time function and the deterministic phases sS, ScS and SS, and the surface wave in S coda. The records of the processed S coda observed at every station-pair are cross-correlated for each teleseismic event. The cross-correlation functions (CCFs) of different earthquakes are stacked to retrieve wavefields propagating between the station-pairs. As a result, direct P, S and the reflected waves could be extracted when the stacked CCFs are aligned as a function of distance of separation between two stations. The gradients of the traveltime curves for the direct P and S waves are approximately 5–7 and 3k m s −1 , respectively. These velocities correspond to the seismic velocities of the crust or the P-wave velocity of the uppermost mantle. In order to enhance the reflected waves, we searched for the reflection points by assuming that the later phases in the CCFs are SS reflections and map the amplitudes onto the depth sections. As a result, the negative phases dipping to the north can be traced right below the region of hypocentre distribution. These negative phases probably correspond to the oceanic Moho within the Philippine Sea slab. We also show that the oceanic Moho can also be traced by assuming PP reflections. These results indicate that the CCFs plausibly contain information regarding both P and S waves propagating between the two receivers, and are capable of detecting reflected phases in addition to the direct waves.

[1]  Roel Snieder,et al.  Application of seismic interferometry to extract P- and S-wave propagation and observation of shear-wave splitting from noise data at Cold Lake, Alberta, Canada , 2008 .

[2]  Charles A. Langston,et al.  Ambient seismic noise tomography and structure of eastern North America , 2008 .

[3]  Peter Gerstoft,et al.  Extracting time‐domain Green's function estimates from ambient seismic noise , 2005 .

[4]  R. Snieder.  Extracting the Green's function from the correlation of coda waves: a derivation based on stationary phase. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  R. Weaver,et al.  On the emergence of the Green's function in the correlations of a diffuse field: pulse-echo using thermal phonons. , 2001, Ultrasonics.

[6]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[7]  C. Wilson,et al.  Teleseismic shot-profile migration , 2006 .

[8]  Hiroshi P. Sato,et al.  Interferometric seismic imaging of crustal structure using scattered teleseismic waves , 2007 .

[9]  Roel Snieder,et al.  Interferometry by deconvolution: Part 2 — Theory for elastic waves and application to drill-bit seismic imaging , 2008 .

[10]  Michael G. Bostock,et al.  Multiparameter two-dimensional inversion of scattered teleseismic body waves 1. Theory for oblique incidence , 2001 .

[11]  Michel Campillo,et al.  Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise , 2004 .

[12]  Keiiti Aki,et al.  Space and Time Spectra of Stationary Stochastic Waves, with Special Reference to Microtremors , 1957 .

[13]  M. Ishida.  Geometry and relative motion of the Philippine Sea plate and Pacific plate beneath the Kanto-Tokai district, Japan , 1992 .

[14]  Gerard T. Schuster,et al.  Crosscorrelogram migration of inverse vertical seismic profile data , 2006 .

[15]  F. Yamazaki,et al.  Subduction of the Philippine Sea Plate beneath the Tokai Area, Central Japan , 1989 .

[16]  Michael G. Bostock,et al.  Multiparameter two‐dimensional inversion of scattered teleseismic body waves 2. Numerical examples , 2001 .

[17]  Walter H. F. Smith,et al.  New version of the generic mapping tools , 1995 .

[18]  Öz Yilmaz,et al.  Seismic data processing , 1987 .

[19]  Hiroyuki Fujiwara,et al.  Recent Progress of Seismic Observation Networks in Japan , 2004 .

[20]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[21]  K. Obara,et al.  Simple relationship between seismic activity along Philippine Sea slab and geometry of oceanic Moho beneath southwest Japan , 2008 .

[22]  T. Iwasaki,et al.  Configuration of subducting Philippine Sea plate and crustal structure in the central Japan region , 2003 .

[23]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[24]  Hiroo Kanamori,et al.  Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events , 1994 .

[25]  Peter Gerstoft,et al.  Surface wave tomography from microseisms in Southern California , 2005 .

[26]  K. Obara,et al.  Three‐dimensional crustal S wave velocity structure in Japan using microseismic data recorded by Hi‐net tiltmeters , 2008 .

[27]  K. Hirahara,et al.  Upper mantle imaging beneath the Japan Islands by Hi-net tiltmeter recordings , 2006 .

[28]  H. Nakahara.  Theoretical background of retrieving Green's function by cross-correlation: one-dimensional case , 2006 .

[29]  J. Nakajima,et al.  Three-dimensional seismic velocity structure and configuration of the Philippine Sea slab in southwestern Japan estimated by double-difference tomography , 2008 .

[30]  Moho and Philippine Sea plate structure beneath central Honshu Island, Japan, from teleseismic receiver functions , 2004 .

[31]  Kees Wapenaar,et al.  Retrieval of reflections from seismic background‐noise measurements , 2007 .

[32]  Gerard T. Schuster,et al.  Autocorrelogram migration: Theory , 2003 .

[33]  H. Nakahara.  A systematic study of theoretical relations between spatial correlation and Green's function in one-, two- and three-dimensional random scalar wavefields , 2006 .

[34]  Michael G. Bostock,et al.  Multiparameter two-dimensional inversion of scattered teleseismic body waves 3. Application to the Cascadia 1993 data set , 2001 .

[35]  Richard L Weaver,et al.  Fluctuations in diffuse field-field correlations and the emergence of the Green's function in open systems. , 2005, The Journal of the Acoustical Society of America.

[36]  J. Nakajima,et al.  Subduction of the Philippine Sea plate beneath southwestern Japan: Slab geometry and its relationship to arc magmatism , 2007 .

[37]  Kees Wapenaar,et al.  Relations between reflection and transmission responses of three-dimensional inhomogeneous media , 2004 .

[38]  R. Snieder,et al.  Correcting for bias due to noise in coda wave interferometry , 2006 .

[39]  勝彦 汐見,et al.  Hi-net/KiK-net 観測点における地中地震計設置方位の推定 , 2003 .

[40]  A. Paul,et al.  Long-Range Correlations in the Diffuse Seismic Coda , 2003, Science.

[41]  Kees Wapenaar,et al.  Spurious multiples in seismic interferometry of primaries , 2006 .

[42]  A. Kato,et al.  High Pore Fluid Pressure May Cause Silent Slip in the Nankai Trough , 2004, Science.

[43]  Michel Campillo,et al.  High-Resolution Surface-Wave Tomography from Ambient Seismic Noise , 2005, Science.

[44]  Jon F. Claerbout,et al.  Synthesis of a layered medium from its acoustic transmission response , 1968 .

[45]  K. Wapenaar.  Retrieving the elastodynamic Green's function of an arbitrary inhomogeneous medium by cross correlation. , 2004, Physical review letters.

[46]  T E Browder,et al.  Search for the lepton-flavor-violating decay tau- -->micro-eta at Belle. , 2004, Physical review letters.

[47]  Peter Gerstoft,et al.  P‐waves from cross‐correlation of seismic noise , 2005 .

[48]  Configuration of subducting Philippine Sea plate and crustal structure in the central Japan region , 2003 .

[49]  T. Kunugi,et al.  Mapping of the 410‐ and 660‐km discontinuities beneath the Japanese islands , 2005 .

[50]  G. Schuster,et al.  Autocorrelogram migration: IVSPWD test , 2003 .