Source rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data

The spatio-temporal slip distribution of the 2003 Tokachi-oki, Japan, earthquake was estimated from teleseismic body wave and strong ground motion data. To perform stable inversion, we applied smoothing constraints to the slip distribution with respect to time and space, and determined the optimal weights of constraints using an optimized Akaike’s Bayesian Information Criterion (ABIC). We found that the rupture propagates mainly along the dip direction, and the length of the rupture area is shorter than its width. The mean rise time in the shallow asperity is significantly longer than that in the deep asperity, which might be attributed to variable frictional properties or lower strength of the plate interface at shallower depths. The average rupture velocity of deep asperity extends to the shear-wave velocity. The derived source parameters are as follows: seismic moment Mo = 1.7×1021 Nm (Mw 8.0); source duration = 50 sec. We also estimated the shear stress change due to the mainshock on and around the major fault zone. It appears that many aftershocks on the plate boundary took place in and adjacent to the zones of stress increase due to the rupture of the mainshock.

[1]  Douglas S. Dreger,et al.  Seismic remote sensing for the earthquake source process and near‐source strong shaking: A case study of the October 16, 1999 Hector Mine earthquake , 2000 .

[2]  Y. Yagi,et al.  Comparison of the coseismic rupture with the aftershock distribution in the Hyuga‐nada Earthquakes of 1996 , 1999 .

[3]  H. Shiobara,et al.  A detailed subduction structure in the Kuril trench deduced from ocean bottom seismographic refraction studies , 1989 .

[4]  H. Kanamori,et al.  Rupture Process of the Kobe, Japan, Earthquake of Jan. 17, 1995, Determined from Teleseismic Body Waves. , 1996 .

[5]  Masayuki Kikuchi,et al.  Inversion of complex body waves , 1982 .

[6]  Stephen H. Hartzell,et al.  Aftershock patterns and main shock faulting , 1988 .

[7]  T. Lay,et al.  Tsunami earthquakes possibly widespread manifestations of frictional conditional stability , 2002 .

[8]  S. Yoshida.  Waveform inversion for rupture process using a non-flat seafloor model: Application to 1986 Andreanof Islands and 1985 Chile earthquakes , 1992 .

[9]  Masayuki Kikuchi,et al.  Inversion of complex body waves—III , 1991, Bulletin of the Seismological Society of America.

[10]  Y. Okada.  Internal deformation due to shear and tensile faults in a half-space , 1992, Bulletin of the Seismological Society of America.

[11]  Richard G. Gordon,et al.  Current plate motions , 1990 .

[12]  M. Takeo,et al.  Fault Heterogeneity of Inland Earthquakes in Japan , 1990 .

[13]  SPACE-TIME PATTERN OF LARGE EARTHQUAKES OCCURRING OFF THE PACIFIC COAST OF THE JAPANESE ISLANDS , 1974 .

[14]  Hirotugu Akaike,et al.  Likelihood and the Bayes procedure , 1980 .

[15]  Yuji Yagi,et al.  Waveform inversion for seismic source processes using ABIC with two sorts of prior constraints: Comparison between proper and improper formulations , 2003 .

[16]  Shigeo Kinoshita,et al.  Kyoshin Net (K-NET) , 1998 .

[17]  Spatial distribution for moment tensor solutions of the 2003 Tokachi-oki earthquake (MJMA = 8.0) and aftershocks , 2004 .

[18]  Thomas H. Heaton,et al.  Inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 Imperial Valley, California, earthquake , 1983 .

[19]  K. Kohketsu.  The extended reflectivity method for synthetic near-field seismograms , 1985 .

[20]  C. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[21]  Akira Nishitani,et al.  Geodetic data inversion using ABIC to estimate slip history during one earthquake cycle with viscoelastic slip-response functions , 2004 .