Numerical Study on Earthquake Energy Partitioning: Relationships Among Radiated Energy, Seismic Moment, and Stress Drop

2020 ◽  
Vol 125 (1) ◽  
Author(s):  
Z. Khademian ◽  
M. Nakagawa ◽  
U. Ozbay
2020 ◽  
Author(s):  
Marie Violay ◽  
Federica Paglialunga ◽  
François X. Passelègue

<p>Earthquakes correspond to a sudden release of elastic energy stored during inter-seismic period by tectonic loading around fault. The earthquake energy budget consists of four non-independent terms: the energy release rate (by unit crack length), the fracture energy, the heat term and finally the radiated energy. These terms depend on the rupture and sliding velocities, the amount of slip and the stress drop. Because of the impossibility to access to stress and strain conditions at depth, the earthquake energy budget cannot be fully constrained from seismological data, limiting our understanding of its influence on rupture propagation.</p><p>To address this issue, we conducted stick-slip experiments with large samples in a biaxial configuration apparatus. By imposing constant normal load and increasing shear load, seismic events were produced on a 20 cm long fault, for which the energy budget was estimated using different methods.</p><p>Fracture energy was estimated by recording the strain field around the crack tip through high frequency (2 MHz) strain gage rosettes array and comparing it to the theoretical LEFM strain field predictions obtained for same conditions (i.e. rupture velocity, distance from the fault). Fracture energies were then inverted and found to range in between 1 and 10 J/m<sup>2</sup>. At the same time the energy partitioning was estimated through stress-slip evolution during rupture. The fracture energies obtained from this method are almost one order of magnitude larger than the ones inverted from LEFM and range in between 1 and 90 J/m<sup>2</sup>. Moreover, the energy partitioning shows the radiated energy ranging between 80 and 300 J/m<sup>2</sup> and finally the heat/thermal energy as the largest fraction of the energy partitioning with values ranging from 200 to 2500 J/m<sup>2</sup>. Our preliminary results highlight the importance of understanding the contribution of heat energy in frictional processes, since this term cannot be estimated from seismological data.</p>


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.


2009 ◽  
Vol 6 (6) ◽  
pp. 10849-10881
Author(s):  
J. Hong ◽  
J. Kim

Abstract. The Tibetan Plateau is a critical region in the research of biosphere-atmosphere interactions on both regional and global scales due to its relation to Asian summer monsoon and El Niño. The unique environment on the Plateau provides valuable information for the evaluation of the models' surface energy partitioning associated with the summer monsoon. In this study, we investigated the surface energy partitioning on this important area through comparative analysis of two biosphere models constrained by the in-situ observation data. Indeed, the characteristics of the Plateau provide a unique opportunity to clarify the structural deficiencies of biosphere models as well as new insight into the surface energy partitioning on the Plateau. Our analysis showed that the observed inconsistency between the two biosphere models was mainly related to: 1) the parameterization for soil evaporation; 2) the way to deal with roughness lengths of momentum and scalars; and 3) the parameterization of subgrid velocity scale for aerodynamic conductance. Our study demonstrates that one should carefully interpret the modeling results on the Plateau especially during the pre-monsoon period.


2020 ◽  
Vol 110 (1) ◽  
pp. 345-356 ◽  
Author(s):  
Itzhak Lior ◽  
Alon Ziv

ABSTRACT Currently available earthquake early warning systems employ region-specific empirical relations for magnitude determination and ground-motion prediction. Consequently, the setting up of such systems requires lengthy calibration and parameter tuning. This situation is most problematic in low seismicity and/or poorly instrumented regions, where the data available for inferring those empirical relations are scarce. To address this issue, a generic approach for real-time magnitude, stress drop, and ground-motion prediction is introduced that is based on the omega-squared model. This approach leads to the following approximate expressions for seismic moment: M0∝RT0.5Drms1.5/Vrms0.5, and stress drop: Δτ∝RT0.5Arms3/Vrms2, in which R is the hypocentral distance; T is the data interval; and Drms, Vrms, and Arms are the displacement, velocity, and acceleration root mean squares, respectively, which may be calculated in the time domain. The potential of these relations for early warning applications is demonstrated using a large composite data set that includes the two 2019 Ridgecrest earthquakes. A quality parameter is introduced that identifies inconsistent earthquake magnitude and stress-drop estimates. Once initial estimates of the seismic moment and stress drop become available, the peak ground velocity and acceleration may be estimated in real time using the generic ground-motion prediction equation of Lior and Ziv (2018). The use of stress drop for ground-motion prediction is shown to be critical for strong ground accelerations. The main advantages of the generic approach with respect to the empirical approach are that it is readily implementable in any seismic region, allows for the easy update of magnitude, stress drop, and shaking intensity with time, and uses source parameter determination and peak ground motion predictions that are subject to the same model assumptions, thus constituting a self-consistent early warning method.


2020 ◽  
Author(s):  
Davide Scafidi ◽  
Daniele Spallarossa ◽  
Matteo Picozzi ◽  
Dino Bindi

<p>Understanding the dynamics of faulting is a crucial target in earthquake source physics (Yoo et al., 2010). To study earthquake dynamics it is indeed necessary to look at the source complexity from different perspectives; in this regard, useful information is provided by the seismic moment (M0), which is a static measure of the earthquake size, and the seismic radiated energy (ER), which is connected to the rupture kinematics and dynamics (e.g. Bormann & Di Giacomo 2011a). Studying spatial and temporal evolution of scaling relations between scaled energy (i.e., e = ER/M0) versus the static measure of source dimension (M0) can provide valuable indications for understanding the earthquake generation processes, single out precursors of stress concentrations, foreshocks and the nucleation of large earthquakes (Picozzi et al., 2019). In the last ten years, seismology has undergone a terrific development. Evolution in data telemetry opened the new research field of real-time seismology (Kanamori 2005), which targets are the rapid determination of earthquake location and size, the timely implementation of emergency plans and, under favourable conditions, earthquake early warning. On the other hand, the availability of denser and high quality seismic networks deployed near faults made possible to observe very large numbers of micro-to-small earthquakes, which is pushing the seismological community to look for novel big data analysis strategies. Large earthquakes in Italy have the peculiar characteristic of being followed within seconds to months by large aftershocks of magnitude similar to the initial quake or even larger, demonstrating the complexity of the Apennines’ faults system (Gentili and Giovanbattista, 2017). Picozzi et al. (2017) estimated the radiated seismic energy and seismic moment from P-wave signals for almost forty earthquakes with the largest magnitude of the 2016-2017 Central Italy seismic sequence. Focusing on S-wave signals recorded by local networks, Bindi et al. (2018) analysed more than 1400 earthquakes in the magnitude ranges 2.5 ≤ Mw ≤ 6.5 of the same region occurred from 2008 to 2017 and estimated both ER and M0, from which were derived the energy magnitude (Me) and Mw for investigating the impact of different magnitude scales on the aleatory variability associated with ground motion prediction equations. In this work, exploiting first steps made in this direction by Picozzi et al. (2017) and Bindi et al. (2018), we derived a novel approach for the real-time, robust estimation of seismic moment and radiated energy of small to large magnitude earthquakes recorded at local scales. In the first part of the work, we describe the procedure for extracting from the S-wave signals robust estimates of the peak displacement (PDS) and the cumulative squared velocity (IV2S). Then, exploiting a calibration data set of about 6000 earthquakes for which well-constrained M0 and theoretical ER values were available, we describe the calibration of empirical attenuation models. The coefficients and parameters obtained by calibration were then used for determining ER and M0 of a testing dataset</p>


1996 ◽  
Vol 86 (5) ◽  
pp. 1623-1626
Author(s):  
Z. L. Wu ◽  
Y. T. Chen ◽  
S. G. Kim

Abstract The physical significance of “earthquake quanta” is discussed in the perspective of a proposed concept “seismon.” It is pointed out that the physical significance of earthquake quanta may be understood in the way similar to that of phonon in solid-state physics. Using the statistical properties of the seismons, the property such that for large earthquakes the mean stress drop is approximately constant, while for small events the mean stress drop increases with the seismic moment, may be obtained, without detailed considerations of the failure criteria and the interaction between different earthquake quanta.


Sign in / Sign up

Export Citation Format

Share Document