Rupture analysis of the 2020 Petrinja earthquake based on seismological observations

2021 ◽  
Author(s):  
Mathieu Causse
Keyword(s):  
Spectral Analysis ◽  
Stress Drop ◽  
Kinematic Model ◽  
Source Model ◽  
Corner Frequency ◽  
Rupture Propagation ◽  
Rupture Velocity ◽  
Geophysical Research ◽  
Apparent Source ◽  
Rupture Duration

<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>

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

2022 ◽  
Author(s):  
Chen Ji ◽  
Ralph J. Archuleta
Keyword(s):  
Stress Drop ◽  
Wave Speed ◽  
Dynamic Stress ◽  
Corner Frequency ◽  
Rupture Velocity ◽  
Scaling Relation ◽  
Shear Wave Speed ◽  
Dynamic Stress Drop ◽  
Static Stress Drop ◽  
Frequency Source

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.

Source parameter analysis from strong motion records of the Friuli, Italy, earthquake sequence (1976-1977)

1987 ◽  
Vol 77 (4) ◽  
pp. 1127-1146
Author(s):  
Giuseppe De Natale ◽  
Raul Madariaga ◽  
Roberto Scarpa ◽  
Aldo Zollo
Keyword(s):  
Frequency Domain ◽  
Stress Drop ◽  
High Frequency ◽  
Epicentral Distance ◽  
Strong Motion ◽  
Source Model ◽  
Parameter Analysis ◽  
Corner Frequency ◽  
Spectral Parameters ◽  
Single Station

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.

Near-field spectral analysis of data-integrative dynamic rupture earthquake simulations of the 1992 Landers earthquake

2020 ◽  
Author(s):  
Nico Schliwa ◽  
Alice-Agnes Gabriel
Keyword(s):  
Spectral Analysis ◽  
Ground Motion ◽  
Near Field ◽  
Slip Distribution ◽  
Corner Frequency ◽  
Dynamic Rupture ◽  
Geophysical Research ◽  
Acoustic Sensing ◽  
Distributed Acoustic Sensing ◽  
Large Earthquakes

<p>The rise of observations from Distributed Acoustic Sensing (e.g., Zhan 2020) and high-rate GNSS networks (e.g., Madariaga et al., 2019) highlight the potential of dense ground motion observations in the near-field of large earthquakes. Here, spectral analysis of >100,000 synthetic near-field strong motion waveforms (up to 2 Hz) is presented in terms of directivity, corner frequency, fall-off rate, moment estimates and static displacements.</p><p>The waveforms are generated in 3‐D large-scale dynamic rupture simulations which incorporate the interplay of complex fault geometry, topography, 3‐D rheology and viscoelastic attenuation (Wollherr et al., 2019). A preferred scenario accounts for off-fault deformation and reproduces a broad range of observations, including final slip distribution, shallow slip deficits, and spontaneous rupture termination and transfers between fault segments. We examine the effects of variations in modeling parameterization within a suite of scenarios including purely elastic setups and models neglecting viscoelastic attenuation. </p><p>First, near-field corner frequency mapping implementing a novel spectral seismological misfit criterion reveals rays of elevated corner frequencies radiating from each slipping fault at 45 degree to rupture forward direction. The azimuthal spectral variations are specifically dominant in the vertical components indicating we map rays of direct P-waves prevailing (Hanks, 1980). The spatial variation in corner frequencies carries information on co-seismic fault segmentation, slip distribution, focal mechanisms and stress drop. Second, spectral fall-off rates are variably inferred during picking the associated corner frequencies to identify the crossover from near-field to far-field spectral behaviour in dependence on distance and azimuth. Third, we determine static displacements with the help of near-field seismic spectra.</p><p>Our findings highlight the future potential of spectral analysis of spatially dense (low frequency) ground motion observations for inferring earthquake kinematics and understanding earthquake physics directly from near-field data; while synthetic studies are crucial to identify "what to look for" in the vast amount of data generated.</p><p><em>References:</em></p><p>Hanks, T.C., 1980. The corner frequency shift, earthquake source models and Q.</p><p>Madariaga, R., Ruiz, S., Rivera, E., Leyton, F. and Baez, J.C., 2019. Near-field spectra of large earthquakes. Pure and Applied Geophysics, 176(3), pp.983-1001.</p><p>Wollherr, S., Gabriel, A.-A. and Mai, P.M., 2019.  Landers 1992 “reloaded”: Integrative dynamic earthquake rupture modeling. Journal of Geophysical Research: Solid Earth, 124(7), pp.6666-6702.</p><p>Zhan, Z., 2020. Distributed Acoustic Sensing Turns Fiber‐Optic Cables into Sensitive Seismic Antennas. Seismological Research Letters, 91(1), pp.1-15.</p>

Source physics interpretation of non-self-similar double-corner frequency source spectral model JA19_2S

2021 ◽  
Author(s):  
Chen Ji ◽  
Ralph Archuleta
Keyword(s):  
Aspect Ratio ◽  
Stress Drop ◽  
Dynamic Stress ◽  
Static Stress ◽  
Corner Frequency ◽  
Rupture Velocity ◽  
Dynamic Stress Drop ◽  
Self Similar ◽  
Static Stress Drop ◽  
Frequency Source

<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>

Far-field pulse shapes from circular sources with variable rupture velocities

1997 ◽  
Vol 87 (5) ◽  
pp. 1288-1296
Author(s):  
Nicholas Deichmann
Keyword(s):  
Stress Drop ◽  
Seismic Moment ◽  
Total Duration ◽  
Source Model ◽  
P Wave ◽  
Far Field ◽  
Rupture Velocity ◽  
P Waves ◽  
S Waves ◽  
Rate Functions

Abstract Recently, Sato (1994) developed a simple earthquake source model of a circular rupture expanding outward from the center of a fault with constant stress drop. In contrast to previous models, the rupture velocity is allowed to vary over the duration of faulting. This model is used to synthesize apparent moment-rate functions for a three-stage source process: first, the rupture starts out with a gradually increasing velocity, then, it continues to expand uniformly until, finally, it slows to a gradual stop. Synthetic velocity seismograms are obtained from a convolution of the apparent moment-rate functions with a causal Q-operator and an appropriate instrument response. Comparisons with an example of an earthquake signal show that, in the context of the proposed model, the observed emergent P-wave onset, which is not compatible with a constant rupture velocity, can be explained by a gradually accelerating rupture front. Systematic departures from the generally expected scaling relationship between seismic moment and rupture duration are often interpreted as evidence for a dependence of stress drop on seismic moment. However, the trade-off between stress drop and rupture velocity inherent in all kinematic source models implies that such deviations can just as well be attributed to systematic variations of rupture velocity. Whereas, in general, the total duration of the far-field displacement pulse is shorter for P waves than for S waves, the model predicts that the rise time, τ1/2, of the displacement pulse should be longer for P waves than for S waves. This feature could constitute a critical test of the model and also provide a constraint on the rupture velocity.

Variation of seismic scalar moment–corner frequency relationship during development of a hydraulic fracture system

2019 ◽  
Vol 38 (2) ◽  
pp. 123-129
Author(s):  
Takashi Mizuno ◽  
Joel Le Calvez ◽  
Jim Rutledge
Keyword(s):  
Hydraulic Fracture ◽  
Stress Drop ◽  
Static Stress ◽  
Corner Frequency ◽  
Fracture System ◽  
Rupture Velocity ◽  
Static Stress Drop ◽  
Two Parameter ◽  
Parameter Relation ◽  
Scalar Moment

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.

A comparison of two methods for earthquake source inversion using strong motion seismograms

1994 ◽  
Vol 37 (6) ◽  
Author(s):  
B. P. Cohee ◽  
G. C. Beroza
Keyword(s):  
Rise Time ◽  
Seismic Moment ◽  
Strong Motion ◽  
Rupture Time ◽  
Rupture Propagation ◽  
Rupture Velocity ◽  
Strong Motion Data ◽  
Motion Data ◽  
Inversion Methods ◽  
Window Method

In this paper we compare two time-domain inversion methods that have been widely applied to the problem of modeling earthquake rupture using strong-motion seismograms. In the multi-window method, each point on the fault is allowed to rupture multiple times. This allows flexibility in the rupture time and hence the rupture velocity. Variations in the slip-velocity function are accommodated by variations in the slip amplitude in each time-window. The single-window method assumes that each point on the fault ruptures only once, when the rupture front passes. Variations in slip amplitude are allowed and variations in rupture velocity are accommodated by allowing the rupture time to vary. Because the multi-window method allows greater flexibility, it has the potential to describe a wider range of faulting behavior; however, with this increased flexibility comes an increase in the degrees of freedom and the solutions are comparatively less stable. We demonstrate this effect using synthetic data for a test model of the Mw 7.3 1992 Landers, California earthquake, and then apply both inversion methods to the actual recordings. The two approaches yield similar fits to the strong-motion data with different seismic moments indicating that the moment is not well constrained by strong-motion data alone. The slip amplitude distribution is similar using either approach, but important differences exist in the rupture propagation models. The single-window method does a better job of recovering the true seismic moment and the average rupture velocity. The multi-window method is preferable when rise time is strongly variable, but tends to overestimate the seismic moment. Both methods work well when the rise time is constant or short compared to the periods modeled. Neither approach can recover the temporal details of rupture propagation unless the distribution of slip amplitude is constrained by independent data.

A theoretical study of the radiation from small strikeslip earthquakes at close distances

1983 ◽  
Vol 73 (1) ◽  
pp. 83-96 ◽  
Author(s):  
Michel Campillo ◽  
Michel Bouchon
Keyword(s):  
Ground Motion ◽  
Epicentral Distance ◽  
Homogeneous Medium ◽  
Source Model ◽  
Strike Slip ◽  
Peak Ground Velocity ◽  
Corner Frequency ◽  
Crack Model ◽  
Sharp Change ◽  
S Wave

abstract We present a study of the seismic radiation of a physically realistic source model—the circular crack model of Madariaga—at close distance range and for vertically heterogeneous crustal structures. We use this model to represent the source of small strike-slip earthquakes. We show that the characteristics of the radiated seismic spectra, like the corner frequency, are strongly affected by the presence of the free surface and by crustal layering, and that they can be considerably different from the ones of the homogeneous-medium far-field solution. The vertical and radial displacement spectra are the most strongly affected. We use this source model to calculate the decay of peak ground velocity with epicentral distance and source depth for small strike-slip earthquakes in California. For distances between 10 and 80 km, the peak horizontal velocity decay is of the form r−1.25 for a 4-km hypocentral depth and r−1.65 for deeper sources. The predominance of supercritically reflected arrivals beyond epicentral distances of 70 to 80 km produces a sharp change in the rate of decay of the ground motion. For most of the cases considered, the peak ground velocity increases between 80 and 100 km. We also show that the S-wave velocity in the source layer is the lower limit of phase velocities associated with significant ground motion.

A dynamic model for far-field acceleration

1982 ◽  
Vol 72 (4) ◽  
pp. 1049-1068
Author(s):  
John Boatwright
Keyword(s):  
Stress Drop ◽  
High Frequency ◽  
Dynamic Stress ◽  
The Body ◽  
Body Waves ◽  
Far Field ◽  
Rupture Velocity ◽  
Frequency Level ◽  
Average Change ◽  
Dynamic Stress Drop

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.

The relations between the corner frequency, seismic moment and source dynamic parameters derived from the spontaneous rupture of a circular fault

2021 ◽  
Vol 228 (1) ◽  
pp. 134-146
Author(s):  
Jian Wen ◽  
Jiankuan Xu ◽  
Xiaofei Chen
Keyword(s):  
Stress Drop ◽  
Seismic Moment ◽  
Dynamic Models ◽  
Scaling Law ◽  
Boundary Integral ◽  
Corner Frequency ◽  
Dynamic Parameters ◽  
Zone Size ◽  
Dynamic Rupture ◽  
Scaling Relationships

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.

Sign in / Sign up

Export Citation Format

Share Document