Source process of the 1995 Kobe earthquake: Determination of spatio-temporal slip distribution by Bayesian modeling

The Kobe (Hyogo-ken Nanbu) earthquake ( MS 6.8) occurred on 16 January 1995, in western Japan, bringing about a tremendous disaster. We constructed its source model by determining a spatio-temporal slip distribution using strong-motion seismograms recorded at a number of stations near the hypocenter. The source area is expressed using both a rough and a detailed subfault discretization system and a source-time function expanded with a finite number of basis functions at each subfault is determined. To perform stable inversions with many parameters, we construct a Bayesian model with unknown hyperparameters, incorporating prior constraints on the smoothness of spatio-temporal slip distribution in the observation equation. Using Akaike's Bayesian Information Criterion (ABIC), the optimal set of hyperparameters are objectively determined from observed data, which enables us to uniquely determine a spatio-temporal slip distribution. After rough analysis in which the actual rupture area is estimated as 45 × 20 km, we carried out a detailed analysis with a finer subfault discretization system (5 × 5 km). In our final source model, most of the slip exists beneath Awaji Island, and the rupture below the city of Kobe, where the damage was most serious, is relatively smaller. The large slips of about 2 m in the shallowest subfaults located southwest of the hypocenter correspond to the surface offsets observed along the Nojima fault. The overall focal mechanism of the source is right-lateral strike slip with a small thrust component, while total seismic moment is 2.0 × 1019 Nm. The source-time functions at shallow subfaults show a longer duration than those at deeper ones, being hypothesized to represent a possible combined effect of low rupture propagation velocity and varying friction behavior in shallow regions of crust.

[1]  H. Quin Dynamic stress drop and rupture dynamics of the October 15, 1979 Imperial Valley, California, earthquake , 1990 .

[2]  Stephen H. Hartzell,et al.  Comparison of seismic waveform inversion results for the rupture history of a finite fault: Application to the 1986 North Palm Springs, California, earthquake , 1989 .

[3]  M. Takeo,et al.  The 1993 Kushiro-Oki, Japan, Earthquake: A high stress-drop event in a subducting slab , 1993 .

[4]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

[5]  M. Takeo The rupture process of the 1989 offshore ito earthquakes preceding a submarine volcanic eruption , 1992 .

[6]  C. Scholz The Mechanics of Earthquakes and Faulting , 1990 .

[7]  D. Wald,et al.  Spatial and temporal distribution of slip for the 1992 Landers, California, earthquake , 1994, Bulletin of the Seismological Society of America.

[8]  Christopher H. Scholz,et al.  The brittle-plastic transition and the depth of seismic faulting , 1988 .

[9]  Shingo Yoshida,et al.  Waveform inversion using ABIC for the rupture process of the 1983 Hindu Kush earthquake , 1989 .

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

[11]  T. Mikumo,et al.  Dynamic rupture analysis: Inversion for the source process of the 1990 Izu‐Oshima, Japan, earthquake (M = 6.5) , 1993 .

[12]  Allen H. Olson,et al.  Finite faults and inverse theory with applications to the 1979 Imperial Valley earthquake , 1982 .

[13]  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 .

[14]  Masaya Matsuura,et al.  Geodetic data inversion using a Bayesian information criterion for spatial distribution of fault slip , 1992 .

[15]  M. Bouchon A simple method to calculate Green's functions for elastic layered media , 1981 .

[16]  T. Heaton Evidence for and implications of self-healing pulses of slip in earthquake rupture , 1990 .

[17]  M. Takeo An inversion method to analyze the rupture processes of earthquakes using near-field seismograms , 1987 .

[18]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[19]  Shamita Das,et al.  Inversion for seismic slip rate history and distribution with stabilizing constraints: Application to the 1986 Andreanof Islands Earthquake , 1990 .

[20]  Reconstruction of dynamic rupture process of an earthquake with constraints of kinematic parameters , 1992 .

[21]  Gregory C. Beroza,et al.  Linearized inversion for fault rupture behavior: Application to the 1984 Morgan Hill, California, earthquake , 1988 .

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

[23]  K. Irikura,et al.  Rupture process of the 1983 Japan Sea (Akita-Oki) earthquake using a waveform inversion method , 1986 .

[24]  Estimation of High-frequency Wave Radiation Areas on the Fault Plane of the 1995 Hyogo-ken Nanbu Earthquake by the Envelope Inversion of Acceleration Seismograms , 1996 .

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

[26]  B. Kennett,et al.  Seismic waves in a stratified half space. , 1979 .

[27]  D. Wald Slip History of the 1995 Kobe, Japan, Earthquake Determined from Strong Motion, Teleseismic, and Geodetic Data , 1996 .

[28]  Shingo Yoshida,et al.  Joint Inversion of Near- and Far-field Waveforms and Geodetic Data for the Rupture Process of the 1995 Kobe Earthquake , 1996 .

[29]  David D. Jackson,et al.  A Bayesian approach to nonlinear inversion , 1985 .

[30]  Thomas H. Heaton,et al.  RUPTURE HISTORY OF THE 1984 MORGAN HILL, CALIFORNIA, EARTHQUAKE FROM THE INVERSION OF STRONG MOTION RECORDS , 1986 .

[31]  Thomas H. Heaton,et al.  Rupture model of the 1989 Loma Prieta earthquake from the inversion of strong-motion and broadband teleseismic data , 1991, Bulletin of the Seismological Society of America.

[32]  M. Takeo,et al.  The dynamic rupture process of the 1993 Kushiro-oki earthquake , 1996 .

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

[34]  S. Yoshida,et al.  Simultaneous inversion of waveform and geodetic data for the rupture process of the 1984 Naganoken–Seibu, Japan, earthquake , 1990 .