On the heterogeneity of the earthquake rupture

Author(s):  
Luca Malagnini ◽  
Douglas S Dreger ◽  
Robert M Nadeau ◽  
Irene Munafò ◽  
Massimo Cocco

Summary The scaling of earthquake parameters with seismic moment and its interpretation in terms of self-similarity is still debated in the literature. We address this question by examining a worldwide compilation of corner frequency-based and elastic rebound theory (ERT)-based fault slip, area and stress drop values for earthquakes ranging in magnitude from -0.7 to 7.8. We find that corner frequency estimates of slip (and stress drop) scale differently than those inferred from the ERT approach, where the latter deviates from the generally accepted constant stress drop behavior of so-called self-similar scaling models. We also find that average slips from finite-source models are consistent with corner frequency scaling, whereas peak slip values are more consistent with the ERT scaling. The different scaling of corner frequency- and ERT-based estimates of slip and stress drop with earthquake size is interpreted in terms of heterogeneity of the rupture process. ERT-based estimates of stress drop decrease with seismic moment suggesting a self-affine behavior. Despite the inferred heterogeneity at all scales, we do not observe a clear effect on the Brune stress drop scaling with earthquake size.

2021 ◽  
Vol 228 (1) ◽  
pp. 134-146
Author(s):  
Jian Wen ◽  
Jiankuan Xu ◽  
Xiaofei Chen

SUMMARY The stress drop is an important dynamic source parameter for understanding the physics of source processes. The estimation of stress drops for moderate and small earthquakes is based on measurements of the corner frequency ${f_c}$, the seismic moment ${M_0}$ and a specific theoretical model of rupture behaviour. To date, several theoretical rupture models have been used. However, different models cause considerable differences in the estimated stress drop, even in an idealized scenario of circular earthquake rupture. Moreover, most of these models are either kinematic or quasi-dynamic models. Compared with previous models, we use the boundary integral equation method to simulate spontaneous dynamic rupture in a homogeneous elastic full space and then investigate the relations between the corner frequency, seismic moment and source dynamic parameters. Spontaneous ruptures include two states: runaway ruptures, in which the rupture does not stop without a barrier, and self-arresting ruptures, in which the rupture can stop itself after nucleation. The scaling relationships between ${f_c}$, ${M_0}$ and the dynamic parameters for runaway ruptures are different from those for self-arresting ruptures. There are obvious boundaries in those scaling relations that distinguish runaway ruptures from self-arresting ruptures. Because the stress drop varies during the rupture and the rupture shape is not circular, Eshelby's analytical solution may be inaccurate for spontaneous dynamic ruptures. For runaway ruptures, the relations between the corner frequency and dynamic parameters coincide with those in the previous kinematic or quasi-dynamic models. For self-arresting ruptures, the scaling relationships are opposite to those for runaway ruptures. Moreover, the relation between ${f_c}$ and ${M_0}$ for a spontaneous dynamic rupture depends on three factors: the dynamic rupture state, the background stress and the nucleation zone size. The scaling between ${f_c}$ and ${M_0}$ is ${f_c} \propto {M_0^{ - n}}$, where n is larger than 0. Earthquakes with the same dimensionless dynamic parameters but different nucleation zone sizes are self-similar and follow a ${f_c} \propto {M_0^{ - 1/3}}$ scaling law. However, if the nucleation zone size does not change, the relation between ${f_c}$ and ${M_0}$ shows a clear departure from self-similarity due to the rupture state or background stress.


1995 ◽  
Vol 38 (2) ◽  
Author(s):  
M. Di Bona ◽  
M. Cocco ◽  
A. Rovelli ◽  
R. Berardi ◽  
E. Boschi

The strong motion accelerograms recorded during the 1990 Eastern Sicily earthquake have been analyzed to investigate source and attenuation parameters. Peak ground motions (peak acceleration, velocity and displacement) overestimate the values predicted by the empirical scaling law proposed for other Italian earthquakes, suggesting that local site response and propagation path effects play an important role in interpreting the observed time histories. The local magnitude, computed from the strong motion accelerograms by synthesizing the Wood-Anderson response, is ML = 5.9, that is sensibly larger than the local magnitude estimated at regional distances from broad-band seismograms (ML = 5.4). The standard omega-square source spectral model seems to be inadequate to describe the observed spectra over the entire frequency band from 0.2 to 20 Hz. The seismic moment estimated from the strong motion accelerogram recorded at the closest rock site (Sortino) is Mo = 0.8 x 1024 dyne.cm, that is roughly 4.5 times lower than the value estimated at regional distances (Mo = 3.7 x 1024 dyne.cm) from broad-band seismograms. The corner frequency estimated from the accelera- tion spectra i.5 J; = 1.3 Hz, that is close to the inverse of the dUl.ation of displacement pulses at the two closest recording sites. This value of corner tì.equency and the two values of seismic moment yield a Brune stress drop larger than 500 bars. However, a corner frequency value off; = 0.6 Hz and the seismic moment resulting from regional data allows the acceleration spectra to be reproduced on the entire available frequency band yielding to a Brune stress drop of 210 bars. The ambiguity on the corner frequency value associated to this earthquake is due to the limited frequency bandwidth available on the strong motion recordil1gs. Assuming the seismic moment estimated at regional distances from broad-band data, the moment magnitude for this earthquake is 5.7. The higher local magnitude (5.9) compared with the moment magnitude (5.7) is due to the weak regional attenuation. Beside this, site amplifications due to surface geology have produced the highest peak ground motions among those observed at the strong motion sites.


2019 ◽  
Vol 5 (1) ◽  
pp. 18-23
Author(s):  
Tri Kusmita ◽  
Kirbani Brotopuspito ◽  
Hetty Triastuty

The source parameters describe the different physical properties of seismic volumes under the volcanoes. Source parameters that can be used to distinguish seismic events that are generated by different types of volcanoes activities. Temporary changes of the spectral source parameters provided a description of the main events during the eruption process.  Source parameters are calculated by correlating the relationship between source frequency at spectral displacement (corner frequency) and source parameters based on spectral sources of the Brune model (1970). The angular frequency obtained by applying the FFT algorithm to the VTA spectral displacement. The source parameters analyzed from this VTA earthquake are the spectral slope, seismic moment, stress drop, length of rupture, moment magnitude and radiation energy. Based on the obtained corner frequency (12 Hz-13 Hz), seismic moment, moment magnitude and energy radiation respectively were at 0.2 -1.9 x 1012 Nm, 0.7 - 2 Mw, and 0.1 - 9.5 x 1015 erg. The length of rupture were from 144.2 to 243.1 m, the spectra slope has 2.1 - 7.8 dB/cm, and stress drop are 0.1 - 7,6 bar. From the results of this study, it can be concluded that the changes of spectra characteristic and fluctuate of source patrameters value of VTA earthquakes was asosiated with the different  volcanic activity of Sinabung. Keywords: spectral, VTA, source parameter, volcanic earthquake


2021 ◽  
Author(s):  
Prosanta Kumar Khan ◽  
Bandana Baruah

<p>We investigate the source parameters of 87 local earthquakes (3.5 ≤ M<sub>L</sub> ≤ 5.0) that occurred in West Brahmaputra basin and its neighbouring area, using body wave displacement spectra. Seismic moment, corner frequency, source dimension and static stress drop are estimated using a grid search method based on the model of circular source. The measured seismic moments, corner frequency and moment magnitude ranges from   to  N-m, 0.7 to 12.1 and 3.0 to 4.8, respectively. The average ratio of corner frequency of P - and S - waves is 2.21.<strong> </strong>The scaling relationship of seismic moment against corner frequency is also studied for various tectonics regimes separately. Median stress drop values of individual earthquake vary from ~ 0.1 to 38.5 MPa, with an average value of about ~ 6 MPa. Spatial variation of stress drop observed for different tectonic unit reveals a higher stress drop values associated with West Brahmaputra basin, Shillong-Mikir plateau and Indo-Myanmar subduction zone suggesting a higher stress accumulation that may increase the probability of higher magnitude earthquake. The empirical relationship between M<sub>L</sub> and M<sub>W</sub> scale is also derived for hazard assessment.</p>


2021 ◽  
Vol 254 ◽  
pp. 02016
Author(s):  
Naylya Sycheva

Based on the method of polarity of signs of P-waves, the focal mechanisms of 1674 earthquakes with M ≥ 1.6, which occurred on the territory of the Bishkek geodynamic proving ground (BGPG) from 1994 to 2020, were determined. Some characteristics of the complete catalogue are presented. Quantitative distributions by the type of mechanisms and diagrams of azimuths of the main stress axes are constructed. A variety of focal mechanism of earthquake is observed, most of them are reverse fault, oblique reverse fault, and horizontal strike-slip fault. The compression axis for most of the events has a north-northwest direction and a sub-horizontal position. For 183, dynamic parameters (DP, source parameters) were obtained: spectral density Ω0, corner frequency f0, scalar seismic moment M0, source radius (Brune radius) r, and stress drop Δσ. The correlations between DP and energy characteristic (magnitude) and scalar seismic moment are investigated. The smallest correlation coefficient was obtained for stress drop.


1975 ◽  
Vol 65 (4) ◽  
pp. 899-913
Author(s):  
Robert B. Herrmann

Abstract The observed relationship between magnitude and duration is shown to be a result of the particular shape of the signal coda as a function of time. If the envelope of the coda follows a t−q relationship with increasing time, then the magnitude, mτ, based on a duration τ is consequently of the form m τ = q log ⁡ 10 τ + r . A study of the duration-magnitude and duration-moment relationships for a set of central United States earthquakes indicates that the linear relationship between mτ and log10τ is valid only over a limited range. The departure from the simple linear dependence is explained in terms of instrumental response and the shift of the source-spectrum corner frequency with increasing event size.


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.


1987 ◽  
Vol 77 (4) ◽  
pp. 1127-1146
Author(s):  
Giuseppe De Natale ◽  
Raul Madariaga ◽  
Roberto Scarpa ◽  
Aldo Zollo

Abstract Time and frequency domain analyses are applied to strong motion data recorded in Friuli, Italy, during 1976 to 1977. An inversion procedure to estimate spectral parameters (low frequency level, corner frequency, and high frequency decay) has been applied to displacement spectra using a simple earthquake source model with a single corner frequency. The data were digitized accelerograms from ENEA-ENEL portable and permanent networks. Instrument-corrected SH waves were selected from a set of 138 three-component, hand-digitized records and 28 automatically digitized records. Thirty-eight events with stations having 8 to 32 km epicentral distance were studied. Different stress drop estimates were performed showing high values (200 to 300 bars, on the average) with seismic moments ranging from 2.8 × 1022 to 8.0 × 1024 dyne-cm. The observation of systematic higher values of Brune stress drop (obtained from corner frequencies) with respect to other time and frequency domain estimates of stress release, and the evidence on time series of multiple rupture episodes suggest that the observed corner frequencies are most probably related to subevent ruptures rather than the overall fault size. Seven events recorded at more than one station show a good correlation between rms, Brune, and dynamic stress drops, and a constant scaling of this parameter as a function of the seismic moment. When single station events are also considered, a slight moment dependence of these three stress drop estimates is observed differently. This may be explained by an inadequacy of the ω−2 high-frequency decay of the source model or by high-frequency attenuation due to propagation effects. The high-frequency cutoff of acceleration spectra indicates the presence of an Fmax in the range of 5 to 14 Hz, except for the stations where local site effects produce spectral peaks.


2021 ◽  
Author(s):  
Meng Li ◽  
Casper Pranger ◽  
Ylona van Dinther

<p>Numerical models are well-suited to overcome limited spatial-temporal observations to understand earthquake sequences, which is fundamental to ultimately better assess seismic hazard. However, high-resolution numerical models in 3D are computationally time and memory consuming. This is not optimal if the aspects of lateral or depth variations within the results are not needed to answer a particular objective. In this study we quantify and summarize the limitations and advantages for simulating earthquake sequences in all spatial dimensions.</p><p> </p><p>We simulate earthquake sequences on a strike-slip fault with rate-and-state friction from 0D to 3D using both quasi-dynamic and fully dynamic approaches. This cross-dimensional comparison is facilitated by our newly developed, flexible code library <em>Garnet</em>, which adopts a finite difference method with a fully staggered grid. We have validated our models using problems BP1-QD & FD and BP4-QD & FD of the SEAS (Sequences of Earthquakes and Aseismic Slip) benchmarks from the Southern California Earthquake Center.</p><p> </p><p>Our results demonstrate that lower-dimensional/quasi-dynamic models are qualitatively similar in terms of earthquake cycle characteristics to their higher-dimensional/fully-dynamic counterparts, while they could be hundreds to millions times faster at the same time. Quantitatively, we observe that certain earthquake parameters, such as stress drop and fracture energy release, can be accurately reproduced in each of these simpler models as well. However, higher dimensional models generally produce lower maximum slip velocities and hence longer coseismic durations. This is mainly due to lower rupture speeds, which result from increased energy consumption along added rupture front directions. In the long term, higher dimensional models produce shorter recurrence interval and hence smaller total slip, which is mainly caused by the higher interseismic stress loading rate inside the nucleation zone. The same trend is also observed when comparing quasi-dynamic models to fully dynamic ones. We extend a theoretical calculation that to first order approximates the aforementioned physical observables in 3D to all other dimensions. These theoretical considerations confirm the same trend as what is observed for stress drop, recurrence interval and total slip across dimensions. These findings on similarities and differences of different dimensional models and a corresponding quantification of computational efficiency can guide model design and facilitate result interpretation in future studies.</p>


Sign in / Sign up

Export Citation Format

Share Document