Tomography of the 1995 Kobe earthquake area: comparison of finite‐frequency and ray approaches

SUMMARY We determined a detailed 3-D crustal model in the 1995 Kobe earthquake (M 7.2) area in southwest Japan using both finite-frequency and ray tomography methods. Our finite-frequency tomography technique is based on the single-scattering theory. The finite-frequency sensitivity kernel derived in this study reflects correctly the sensitivity of the heterogeneity off the geometrical ray path and the existence of Fresnel volume, and the kernel depends on the dominant frequency of the observed wave. The dominant frequency is estimated directly from the earthquake magnitude based on a relation that is obtained by regressively analyzing the displacement spectra of 20 earthquakes in the study area. We used 141 118 P-wave and 133 648 S-wave high-quality arrival-time data from 2813 Kobe aftershocks and 3140 other local earthquakes during 1995–2010. The tomographic images obtained with the finite-frequency and ray tomography methods show a high level of similarity, which is verified quantitatively by adopting the structural similarity index. Our results show that the Kobe main shock hypocentre is located in a distinctive zone characterized by a high Poisson's ratio and a low product VP×VS of P- and S-wave velocities, which is interpreted as a fluid-filled, fractured rock matrix that may have triggered the 1995 Kobe earthquake.

[1]  F. Tajima,et al.  Seismic scaling relations and aftershock activity from the sequences of the 2004 mid Niigata and the 2005 west off Fukuoka earthquakes (MW 6.6) in Japan , 2007 .

[2]  H. Kanamori,et al.  Tomography of the Source Area of the 1995 Kobe Earthquake: Evidence for Fluids at the Hypocenter? , 1996, Science.

[3]  R. Hilst,et al.  Reply to comment by F.A. Dahlen and G. Nolet on , 2005 .

[4]  Walter H. F. Smith,et al.  New, improved version of generic mapping tools released , 1998 .

[5]  T. Mizuno,et al.  Crack density and saturation rate in the 1995 Kobe Earthquake Region , 1999 .

[6]  Kojiro Irikura,et al.  Source model of the 1995 Hyogo-Ken Nanbu earthquake and simulation of near-source ground motion , 1998, Bulletin of the Seismological Society of America.

[7]  Vlastislav Cerveny,et al.  Fresnel volume ray tracing , 1992 .

[8]  E. Humphreys,et al.  Seismic heterogeneity and small‐scale convection in the southern California upper mantle , 2010 .

[9]  Maarten V. de Hoop,et al.  On sensitivity kernels for ‘wave-equation’ transmission tomography , 2005 .

[10]  Guust Nolet,et al.  Computing traveltime and amplitude sensitivity kernels in finite-frequency tomography , 2007, J. Comput. Phys..

[11]  Hiroo Kanamori,et al.  Effect of directivity on estimates of radiated seismic energy , 2004 .

[12]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[13]  S. Ueki,et al.  Crustal tomography under the Median Tectonic Line in Southwest Japan using P and PmP data , 2009 .

[14]  Hiroo Kanamori,et al.  Static and Dynamic Scaling Relations for Earthquakes and Their Implications for Rupture Speed and Stress Drop , 2004 .

[15]  Guust Nolet,et al.  Fréchet kernels for finite-frequency traveltimes—I. Theory , 2000 .

[16]  Yoji Kobayashi,et al.  Dehydration of serpentinized slab mantle: Seismic evidence from southwest Japan , 2001 .

[17]  Bo Holm Jacobsen,et al.  Sensitivity Kernels for Time-Distance Inversion , 2000 .

[18]  M. Santosh,et al.  Dissecting large earthquakes in Japan: Role of arc magma and fluids , 2010 .

[19]  O. P. Mishra,et al.  Influence of fluids and magma on earthquakes: seismological evidence , 2002 .

[20]  Nicholas Rawlinson,et al.  Seismic tomography: a window into deep Earth , 2010 .

[21]  Jean Virieux,et al.  Three‐dimensional seismic tomography from P wave and S wave microearthquake travel times and rock physics characterization of the Campi Flegrei Caldera , 2005 .

[22]  T. Shimamoto,et al.  Internal structure and permeability of the Nojima fault, southwest Japan , 2008 .

[23]  J. Virieux,et al.  Finite‐frequency tomography in a crustal environment: Application to the western part of the Gulf of Corinth , 2008 .

[24]  M. Salah,et al.  3-D seismic structure of Kii Peninsula in southwest Japan: evidence for slab dehydration in the forearc , 2003 .

[25]  Y. Ohno,et al.  Crustal structure, fault segmentation, and activity of the Median Tectonic Line in Shikoku, Japan , 2006 .

[26]  W. Ellsworth,et al.  Apparent break in earthquake scaling due to path and site effects on deep borehole recordings , 2003 .

[27]  J. Lees,et al.  Poisson's ratio and porosity at Coso geothermal area, California , 2000 .

[28]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

[29]  Liangguo Dong,et al.  Sensitivity kernels for seismic Fresnel volume tomography , 2009 .

[30]  Guust Nolet,et al.  Comment on ‘On sensitivity kernels for ‘wave-equation’ transmission tomography’ by de Hoop and van der Hilst , 2005 .

[31]  K. Aki.  Scaling law of seismic spectrum , 1967 .

[32]  Akira Hasegawa,et al.  Tomographic imaging of P and S wave velocity structure beneath northeastern Japan , 1992 .

[33]  Charles S. Mueller,et al.  Source pulse enhancement by deconvolution of an empirical Green's function , 1985 .

[34]  Hiroo Kanamori,et al.  The physics of earthquakes , 2004 .

[35]  J. Virieux,et al.  A New Passive Tomography of the Aigion Area (Gulf of Corinth, Greece) from the 2002 Data Set , 2006 .

[36]  L. Chiao,et al.  Imaging seismic velocity structure beneath the Iceland hot spot: A finite frequency approach , 2004 .

[37]  Z. Suzuki,et al.  Objective estimation of source parameters and local Q values by simultaneous inversion method , 1982 .

[38]  Richard M. Allen,et al.  Automatic detection and rapid determination of earthquake magnitude by wavelet multiscale analysis of the primary arrival , 2006 .

[39]  Susan E. Hough,et al.  Empirical Green's function analysis: Taking the next step , 1997 .

[40]  D. Komatitsch,et al.  Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method , 2004 .

[41]  W. Friederich.  The S-velocity structure of the East Asian mantle from inversion of shear and surface waveforms , 2003 .

[42]  S. Hough.  Empirical Green's Function Analysis of Recent Moderate Events in California , 2001 .

[43]  Kim B. Olsen,et al.  Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models , 2005 .

[44]  A fault model of the 1995 Kobe earthquake derived from the GPS data on the Akashi Kaikyo Bridge and other datasets , 1998 .

[45]  M. Ritzwoller,et al.  Empirically determined finite frequency sensitivity kernels for surface waves , 2010 .

[46]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[47]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[48]  H. Kanamori,et al.  Source Properties of Earthquakes near the Salton Sea Triggered by the 16 October 1999 M 7.1 Hector Mine, California, Earthquake , 2002 .

[49]  L. Boschi.  Global multiresolution models of surface wave propagation: comparing equivalently regularized Born and ray theoretical solutions , 2006 .

[50]  Roel Snieder,et al.  Wavefield smoothing and the effect of rough velocity perturbations on arrival times and amplitudes , 1996 .

[51]  Guust Nolet,et al.  Three-dimensional sensitivity kernels for finite-frequency traveltimes: the banana–doughnut paradox , 1999 .

[52]  R. Snieder,et al.  The effect of small-scale heterogeneity on the arrival time of waves , 2001 .