An investigation of the complexity of the earthquake source time function using dynamic faulting models

By constructing more than 30 numerical models of spontaneous earthquake rupture on heterogeneous faults, the complexity of the far-field pulse shapes (the source time function), that is, the pulse consisting of several distinct subpulses, is investigated. The influence of multiple strength barriers and/or dry friction on the far-field waveforms has been studied using the simple model of a square fault with uniformly spaced square barriers. A central asperity is allowed to break to initiate the dynamic process. It is found that the pulse duration is shorter the greater the density of barriers (ratio of barrier size to barrier spacing). The presence of friction has been found to have three different consequences. First, the slipping area is reduced giving rise to a shorter pulse duration. Second, the rate of slip propagation is reduced which increases the pulse duration. Third, the frictional stress drop, if present, increases the amplitude of the pulse. The interaction of these three effects together with the effect of the barriers results in different pulses depending on the relative role of these effects. Breaking of multiple asperities with relatively short time delay between their fracturing produced a single pulse with ripples in its rising portion, the pulse duration depending on the total time required for the asperities to break. Breaking of, say, two asperities with sufficient time delay between them produced a double pulse which is in accord with the asperity model of complex earthquakes. In one case of a fault with barriers and friction, an additional peak was observed in some directions from the source, the peak form and amplitude being essentially the same as those of the peaks due to breaking of asperities. One may conclude that each subpulse observed on a complex far-field pulse is not necessarily interpretable as being due to the failure of an asperity on the fault. Furthermore, for the subpulses which can be related to asperity fracture, the duration is not determined by the size of the corresponding asperity but rather by the physical conditions on the fault surrounding the asperity.

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

[2]  Takashi Miyatake,et al.  NUMERICAL SIMULATIONS OF EARTHQUAKE SOURCE PROCESS BY A THREE-DIMENSIONAL CRACK MODEL. , 1980 .

[3]  Keiiti Aki,et al.  Fault plane with barriers: A versatile earthquake model , 1977 .

[4]  Hiroo Kanamori,et al.  Seismological aspects of the Guatemala Earthquake of February 4, 1976 , 1978 .

[5]  H. Kanamori,et al.  Effects of fault interaction on moment, stress drop, and strain energy release , 1981 .

[6]  Keiiti Aki,et al.  A numerical study of two-dimensional spontaneous rupture propagation , 1977 .

[7]  D. J. Andrews,et al.  Bulletin of the Seismological Society of America , 1985 .

[8]  Breaking of a single asperity: Rupture process and seismic radiation , 1983 .

[9]  J. Brune Tectonic stress and the spectra of seismic shear waves from earthquakes , 1970 .

[10]  Göran Ekström,et al.  A very broad band analysis of the Michoacan, Mexico, Earthquake of September 19, 1985 , 1986 .

[11]  T. Mikumo,et al.  Earthquake sequences on a frictional fault model with non-uniform strengths and relaxation times , 1979 .

[12]  S. Das,et al.  On the Numerical Boundary Integral Equation Method for Three-Dimensional Dynamic Shear Crack Problems , 1987 .

[13]  T. Miyatake Numerical simulations of earthquake source process by a three-dimensional crack model. Part II. Seismic waves and spectrum.:PART II. SEISMIC WAVES AND SPECTRUM , 1980 .

[14]  R. Madariaga Dynamics of an expanding circular fault , 1976, Bulletin of the Seismological Society of America.

[15]  Source complexity of small earthquakes near Matsushiro, Japan , 1989 .

[16]  S. Day Three-dimensional simulation of spontaneous rupture: The effect of nonuniform prestress , 1982, Bulletin of the Seismological Society of America.

[17]  J. Virieux,et al.  Dynamic faulting studied by a finite difference method : Bull seismol soc am, V72, N2, April 1982, P345–369 , 1982 .

[18]  Shamita Das,et al.  A numerical method for determination of source time functions for general three‐dimensional rupture propagation , 1980 .

[19]  Raul Madariaga,et al.  On the relation between seismic moment and stress drop in the presence of stress and strength heterogeneity , 1979 .