Seismic radiation from circular cracks growing at variable rupture velocity

1994 ◽  
Vol 84 (4) ◽  
pp. 1199-1215
Author(s):  
Tamao Sato

Abstract We have investigated the far-field body waves emitted by a circular fault growing at variable rupture velocity. The slip motion on the fault is constructed by assuming that the self-similar slip distribution holds at every successive instant of rupture formation for a circular crack. The present method permits us to evaluate not only the high-frequency but also the low-frequency radiation. It is also possible to handle both abrupt and continuous change in rupture velocity. The resultant expression for the radiation does not involve an integral. It is expressed in a closed form that contains information about the isochrone properties of the two extreme points having the minimum and maximum distances from the observer. Because of its simplicity and high computational speed, the present method provides us with a basic tool for simulating the radiation from seismic sources with a multitude of irregular rupture growth. The property of acceleration pulse radiated by continuous change in rupture velocity has been investigated. The time duration of the velocity change is assumed to be short but to take a finite time. The shape of the rupture front that is effective in radiating the high-frequency radiation is presumed to be semicircular. The acceleration pulse demonstrates a directivity with respect to both the pulse width and amplitude. The pulse width is, on average, given by the duration of the velocity change, and it is modulated by the directivity factor, which depends on the rupture velocity averaged over the change in rupture velocity. The pulse width radiated toward the growth direction of the semi-circular rupture front is shorter than that radiated toward the opposite side. The spectral amplitude of the acceleration pulse depends linearly on the strain, the radius of the rupture front at which the rupture velocity starts to change, the magnitude of the change in rupture velocity, and the generalized radiation pattern coefficient. The directivity of the radiation pattern coefficient is stronger than that for the case of an abrupt change in rupture velocity. In the case where the rupture stops completely, the radiation pattern coefficient is simply controlled by the magnitude of the change in rupture velocity. If we deal with a partial drop of the rupture velocity, however, we must consider the average rupture velocity as well to fully describe the radiation pattern. A procedure is presented for retrieving the model parameters from acceleration pulse data. The present results with regard to the acceleration pulse do not require the coherent movement of the rupture front over the entire circle. The results are applicable to more general cases where the rupture front moves coherently over a certain restricted segment, though the range of the segment has to become wider as the duration of the change in rupture velocity increases.

1982 ◽  
Vol 72 (4) ◽  
pp. 1049-1068
Author(s):  
John Boatwright

abstract A model for the far-field acceleration radiated by an incoherent rupture is constructed by combining Madariaga's (1977) theory for the high-frequency radiation from crack models of faulting with a simple statistical source model. By extending Madariaga's results to acceleration pulses with finite durations, the peak acceleration of a pulse radiated by a single stop or start of a crack tip is shown to depend on the dynamic stress drop of the subevent, the total change in rupture velocity, and the ratio of the subevent radius to the acceleration pulse width. An incoherent rupture is approximated by a sample from a self-similar distribution of coherent subevents. Assuming the subevents fit together without overlapping, the high-frequency level of the acceleration spectra depends linearly on the rms dynamic stress drop, the average change in rupture velocity, and the square root of the overall rupture area. The high-frequency level is independent, to first order, of the rupture complexity. Following Hanks (1979), simple approximations are derived for the relation between the rms dynamic stress drop and the rms acceleration, averaged over the pulse duration. This relation necessarily depends on the shape of the body-wave spectra. The body waves radiated by 10 small earthquakes near Monticello Dam, South Carolina, are analyzed to test these results. The average change of rupture velocity of Δv = 0.8β associated with the radiation of the acceleration pulses is estimated by comparing the rms acceleration contained in the P waves to that in the S waves. The rms dynamic stress drops of the 10 events, estimated from the rms accelerations, range from 0.4 to 1.9 bars and are strongly correlated with estimates of the apparent stress.


1984 ◽  
Vol 74 (6) ◽  
pp. 2061-2082
Author(s):  
Paul Spudich ◽  
L. Neil Frazer

Abstract We analyze the problem of calculating high-frequency ground motions (>1 Hz) caused by earthquakes having arbitrary spatial variations of rupture velocity and slip velocity (or stress drop) over the fault. We approximate the elastic wave Green's functions by far-field body waves, which we calculate using geometric ray theory. However, we do not make the traditional Fraunhofer approximation, so our method may be used close to large faults. The method is confined to high frequencies (greater than about 1 Hz) due to the omisson of near-field terms, and must be used at source-observer distances less than a few source depths, due to the omission of surface waves. It is easily used in laterally varying velocity structures. Assuming a simple parameterization of the slip function, the computational problem collapses to the evaluation of a series of line integrals over the fault, with one line integral per each time ti in the observer seismogram. The path of integration corresponding to observation time ti consists of only those points on the fault which radiate body waves arriving at the observer at exactly time ti. This path is an isochron of the arrival time function. An isochron velocity may be defined that depends on rupture velocity and resembles the usual directivity function. Observed ground motions are directly dependent upon this isochron velocity. Ground velocity is proportional to isochron velocity and ground acceleration is proportional to isochron acceleration in dislocation models of rupture. Ground acceleration may also be related to spatial variations of slip velocity on the fault, using the square of isochron velocity as a constant of proportionality. We show two rupture models, one with variable slip velocity and the other with variable rupture velocity, that cause the same ground acceleration at a single observer. The computational method is shown to produce reasonably accurate synthetic seismograms, compared to a method using complete Green's functions, and requires about 0.5 per cent of the computer time. It may be very effective in calculating ground motions in the frequency band 1 to 10 Hz at observers within a few source depths of large earthquakes, where most of the high-frequency motions may be caused by direct P and S waves. We suggest a possible method for inverting ground motions for both slip velocity and rupture velocity over the fault.


1978 ◽  
Vol 68 (5) ◽  
pp. 1253-1263
Author(s):  
W. H. Bakun ◽  
R. M. Stewart ◽  
C. G. Bufe

abstract On December 12, 1972 at 0351 and 0355 GMT, two earthquakes with magnitudes equal to 3.0 and 2.8, respectively, occurred on the Cienega Road section of the San Andreas fault in central California. The two events have the same hypocenter location and fault-plane soultion. Observed seismograms for these two events at 28 stations within about 65 km of and surrounding the epicenters are systematically different in a pattern that is consistent with different directions of rupture expansion for the two events. The 0351 GMT event preferentially radiated high-frequency (f ⪚ 10 Hz) body waves to the southeast consistent with unilateral rupture propagation toward the southeast while the 0355 GMT event rupture expanded more toward the northwest.


1983 ◽  
Vol 73 (3) ◽  
pp. 723-734 ◽  
Author(s):  
Michel Campillo

abstract We compute the near-field, high-frequency radiation from a circular crack expanding with constant rupture velocity and discuss the characteristics of the stopping phases. We then introduce rupture velocity jumps in the fracture process. The computed accelerations show the dominant role played by the rupture front kinematics. The high acceleration pulses are associated with sudden changes of the rupture velocity. For a sudden jump (or a sudden stop), there is no theoretical high-frequency limit to the spectral density of acceleration. In order to account for fmax, we introduce a smooth deceleration of the rupture front over a time t′ in place of a sudden stop. This results in a spectral fall-off for frequencies greater than 1/t′ and supports the interpretation of fmax as a source effect.


Sign in / Sign up

Export Citation Format

Share Document