A Source Physics Interpretation of Nonself-Similar Double-Corner-Frequency Source Spectral Model JA19_2S

Author(s):  
Chen Ji ◽  
Ralph J. Archuleta

Abstract We investigate the relation between the kinematic double-corner-frequency source spectral model JA19_2S (Ji and Archuleta, 2020) and static fault geometry scaling relations proposed by Leonard (2010). We find that the nonself-similar low-corner-frequency scaling relation of JA19_2S model can be explained using the fault length scaling relation of Leonard’s model combined with an average rupture velocity ∼70% of shear-wave speed for earthquakes 5.3 < M< 6.9. Earthquakes consistent with both models have magnitude-independent average static stress drop and average dynamic stress drop around 3 MPa. Their scaled energy e˜ is not a constant. The decrease of e˜ with magnitude can be fully explained by the magnitude dependence of the fault aspect ratio. The high-frequency source radiation is generally controlled by seismic moment, static stress drop, and dynamic stress drop but is further modulated by the fault aspect ratio and the relative location of the hypocenter. Based on these two models, the commonly quoted average rupture velocity of 70%–80% of shear-wave speed implies predominantly unilateral rupture.

2021 ◽  
Author(s):  
Chen Ji ◽  
Ralph Archuleta

<p>Source spectral models developed for strong ground motion simulations are phenomenological models that represent the average effect that the source processes have on near fault ground motion. Their parameters are directly regressed from the observations and often do not have clear meaning for the physics of the source process. We investigate the relation between the kinematic double-corner frequency (DCF) source spectral model JA19_2S (Ji and Archuleta, BSSA, 2020) and static fault geometry scaling relations proposed by Leonard (2010). We derive scaling relations for the low and high corner frequency in terms of static stress drop, dynamic stress drop, fault rupture velocity, fault aspect ratio, and relative hypocenter location. We find that the non-self-similar low corner frequency  scaling relation of JA19_2S model for 5.3<<strong>M</strong><6.9 earthquakes is well explained using the fault length scaling relation of Leonard’s model combined with a constant rupture velocity. Earthquakes following both models have constant average static stress drop and constant average dynamic stress drop. The high frequency source radiation is controlled by seismic moment, static stress drop and dynamic stress drop but strongly modulated by the fault aspect ratio and the hypocenter’s relative location. The mean, scaled energy  (or apparent stress) decreases with magnitude due to the magnitude dependence of the fault aspect ratio. Based on these two models, the commonly quoted average rupture velocity of 70-80% of shear wave speed implies predominantly unilateral rupture.</p>


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.


2019 ◽  
Vol 38 (2) ◽  
pp. 123-129
Author(s):  
Takashi Mizuno ◽  
Joel Le Calvez ◽  
Jim Rutledge

We propose to utilize the corner frequency and seismic scalar moment relation as a new approach to monitor temporal changes of static stress drop as well as rupture velocity during development of a hydraulic fracture system. We introduce a single parameter M1 to describe a two-parameter relation (scalar moment and corner frequency relation) and analyze temporal variation of this two-parameter relation. Because M1 relates rupture velocity and static stress drop, we can infer temporal variation of rupture velocity and stress drop quantitatively. The parameter M1 is calculated in two case studies. We document that two types of fracturing processes exist: (1) stable rupture velocity and static stress drop during the development of rupture and (2) increase of rupture velocity and/or static stress drop while the fracture system develops. In the latter case, one possible scenario is increase of permeability at each fracture plane during development of the fracture system.


1982 ◽  
Vol 72 (6A) ◽  
pp. 2017-2036
Author(s):  
George L. Choy ◽  
John Boatwright

abstract The Miyagi-Oki earthquake of 12 June 1978, a large (Ms 7.8) interplate thrust event, occurred in a region which had not experienced earthquakes of magnitude greater than 7 since 1938. A sequence of four moderate-sized (5.4 < mb < 6.1) earthquakes encircled the rupture zone of the Miyagi-Oki earthquake over a period of 2 yr before the main shock. Broadband displacement and velocity records of body waves recorded digitally by stations of the Global Digital Seismograph Network are analyzed to determine the static and dynamic characteristics of the sequence. These characteristics include moment, radiated energy, dynamic and static stress drop, and apparent stress. Inversions of duration measurements made on the velocity waveforms permit quantifying the complexity of an event as well as constraining its rupture geometry. Intervals of 7 to 8 months separated the first three events; the main shock occurred 4 months after the third event. The rupture process of the third event was relatively complex; the event also had a substantially higher dynamic stress drop (175 bars) than did the stress drops of the first two events (9 and 10 bars, respectively). These observations suggest that the third event was an interme-diate-term precursor to the main shock. The fourth event, a short-term precursor to the main shock, occurred about 8 min before the main shock. Its dynamic stress drop (20 bars) was lower than that of the third event but higher than that of the first two events.


2020 ◽  
Vol 110 (1) ◽  
pp. 241-249
Author(s):  
Kazuhito Hikima ◽  
Akihiro Shimmura

ABSTRACT For crustal earthquakes, the scaling relationship between the seismic moment M0 and rupture area S varies with the size of the earthquake, due to the limited thickness of the seismogenic layer. In those M0–S scaling relations, in most cases, the calculated static stress drop is altered with the size of earthquake, although the change depends on the assumed fault model. However, it is not clear whether the dependence of the stress drop on M0 is physically reasonable. In this study, the scaling relation between M0 and S, which assumes a constant stress drop over a wide M0 range, is discussed based on the analytical stress drop formula of a rectangular strike-slip fault. In the proposed relation, M0 is proportional to S3/2 for small and medium faults and to S1 for long faults. In addition, the relation between M0 and S varies in the intermediate range, depending on the aspect ratio. The scaling relation showed good agreement with past event data when the saturated rupture width was set to around 15–20 km and the stress drop was set to about 3 MPa.


2021 ◽  
Author(s):  
Mathieu Causse

<p>Here, I use seismological observations (~70 broadband stations at distances between 100 and 400 km from the source) to characterize the rupture properties of the Petrinja mainshock (Mw 6.4). First, I perform a spectral analysis of the P-waves to compute the corner frequency. In order to remove the wave propagation effects and isolate the source properties, I use the largest foreshocks and aftershocks (Mw>4) as empirical Green’s functions (EGFs). Assuming a Brune’s source model, the obtained stress drop is ~20 MPa. This rather large value is in agreement with the short rupture length of ~8 km inferred by InSAR data (Ganas et al. 2021). In addition, the weak azimuthal variations of the corner frequencies indicates a bilateral rupture, that is a rupture nucleating close to the fault center. Second, I compute the apparent source time functions (i.e. the source time functions “seen” from any station) using an EGF deconvolution approach. The results indicate an average rupture duration of 5-6 s with weak azimuthal variation of the apparent rupture duration, in agreement with the spectral analysis. Finally, I perform a Bayesian inversion of the apparent source function, in order to obtain a kinematic model of the rupture propagation (slip distribution, rupture velocity). The preliminary results reveal a slow velocity of the rupture propagation. Such a slow rupture velocity associated with a large stress drop has been observed on other faults with slow slip rates (e.g. Causse et al. 2017). This work provides insight on the rupture process of this major event on a poorly documented fault. I am fully open for collaborations to further develop and enrich this study.</p><p><br><strong>References</strong><br>Causse, M., G. Cultrera, L. Moreau, A. Herrero, E. Schiapappietra and F. Courboulex. Bayesian rupture imaging in a complex medium. The 29 May 2012 Emilia, Northern Italy, earthquake (2017), Geophysical Research Letters, DOI : 10.1002/2017GL074698.<br>Ganas, A., Elias, P., Valkaniotis, S., Tsironi, V., Karasante, I., Briole, P., 2021, Petrinja earthquake moved crust 10 feet, Temblor, http://doi.org/10.32858/temblor.156</p>


1981 ◽  
Vol 71 (4) ◽  
pp. 939-957
Author(s):  
Ralph J. Archuleta ◽  
Stephen H. Hartzell

abstract Near-source ground motion at four azimuths but constant epicentral range is computed from a buried circular strike-slip fault in a half-space. Particle acceleration, velocity, and displacement at each station on the free surface is computed in the frequency band 0.0 to 5.0 Hz. The assumed dislocation is derived from the Kostrov (1964) displacement function for a continuously propagating stress relaxation. The azimuthal variations in the amplitudes and waveforms directly result from spatially varying slip on the fault, spatially varying radiation pattern over the fault, and the magnitude and direction of the rupture velocity. The near-source ground motions are dominated by the rupture in the direction of the receiver. Using a 100-bar effective stress (initial stress minus sliding friction) in a Poisson solid with β = 3.0 km/sec the shear wave speed, and shear modulus μ = 3.0 × 1011 dyne/cm2, the simulated earthquake has a moment Mo = 4.5 × 1025 dyne-cm. Using a rupture velocity of 0.9β, the peak acceleration is 1195 cm/sec2 and velocity 104 cm/sec for the receiver directly on strike. For a receiver 30° off strike, the maximum acceleration 236 cm/sec2 occurs on the vertical component.


Sign in / Sign up

Export Citation Format

Share Document