Fast and Robust Earthquake Source Spectra and Moment Magnitudes from Envelope Inversion

ABSTRACT With the present study, we introduce a fast and robust method to calculate the source displacement spectra of small earthquakes on a local to regional scale. The work is based on the publicly available Qopen method of full envelope inversion, which is further tuned for the given purpose. Important source parameters—seismic moment, moment magnitude, corner frequency, and high-frequency fall off—are determined from the source spectra by fitting a simple earthquake source model. The method is demonstrated by means of a data set comprising the 2018 West Bohemia earthquake swarm. We report moment magnitudes, corner frequencies, and centroid moment tensors inverted from short-period body waves with the Grond package for all earthquakes with a local magnitude larger than 1.8. Moment magnitudes calculated by envelope inversion show a very good agreement to moment magnitudes resulting from the probabilisitc moment tensor inversion. Furthermore, source displacement spectra from envelope inversion show a good agreement with spectra obtained by multiple taper analysis of the direct onsets of body waves but are not affected by the large scatter of the second. The seismic moments obtained with the envelope inversion scale with corner frequencies according to M0∝fc−4.7. Earthquakes of the present data set result in a smaller stress drop for smaller magnitudes. Self-similarity of earthquake rupture is not observed. In addition, we report frequency-dependent site amplification at the used stations.

Earthquake Source Characteristics and S-Wave Propagation Attenuation in the Junction of the Northwest Tarim Basin and Kepingtage Fold-and-Thrust Zone

We separated the propagation path attenuation and source spectra from the S-wave Fourier amplitude spectra of the observed ground motions recorded during 46 small-to-moderate earthquakes in the junction of the northwest Tarim Basin and Kepingtage fold-and-thrust zone, mainly composed of two Jiashi seismic sequences in 2020 and 2018. Slow seismic wave decay was observed as the distance increased, while the quality factor regressed as 60.066 f0.988 for frequency f = 0.254–30 Hz reflects the strong anelastic attenuation in the study region. We estimated the stress drops for the 46 earthquakes under investigation from the preferred corner frequencies and seismic moments by fitting the inverted source spectra and the theoretical ω-square model. The relationship between seismic moment and corner frequency and the dependence of the stress drop on the moment magnitude reveal the breakdown of earthquake self-similar scaling for the events in this study. The temporal variation in stress drops indicates that the mainshock plays a short-term role in the source characteristics of the surrounding earthquakes. Aftershocks immediately following the mainshock show a low stress release and then gradually recover in a short time. The healing process for the fractured fault in the mainshock may be one reason for the stress drop recovery of the aftershock. The foreshock with the low stress release occurring in the high-heterogeneity fault zone may motivate the following occurrence of the largest magnitude mainshock with a high stress drop. We inferred that the foreshock-mainshock behavior, including several moderate events, may be predisposed to occur in our study region characterized by an inhomogeneous crust.

Reliability of Source Parameters for Small Events in Central Italy: Insights from Spectral Decomposition Analysis Applied to Both Synthetic and Real Data

ABSTRACT We apply a spectral decomposition approach to isolate the source spectra from propagation and site effects and, in turn, to estimate the source parameters of small-to-moderate earthquakes that occurred in central Italy. The data set is composed of about 400,000 waveforms relevant to 4111 earthquakes in the moment magnitude range 1.5–6.5, recorded by a high-density network of stations installed in the study area. We first investigate the reliability of the source parameters for small magnitudes through numerical simulations. We generate synthetic spectra for different source scaling models and near-surface attenuation effects, considering the source–station geometry and the data availability of the central Italy data set. Our analysis with synthetics shows that the spectral decomposition is effective in isolating the source contributions from other factors. Moreover, the analysis of the residual distributions suggests that moment magnitude 1.8 is the lower bound for the retrieval of reliable Brune’s source parameters, although we observe an increase of residual’s variability below magnitude 3, and the estimated source parameters could be biased below magnitude 2.3. Remarkably, the assessment of the stress drop Δσ for small events is strongly hampered by site-specific attenuation near the surface. In view of the results with synthetics, we analyze the source parameters of earthquakes recorded in central Italy. The corner frequency versus seismic moment relationship describes a source scaling in which Δσ increases with increasing moment magnitude Mw, the mean Δσ varying from 0.1 MPa for Mw<2 to 7.9 MPa for Mw>5. In particular, Δσ increases mainly for Mw in the ranges 2.5–3 and 4.5–5.2. The corner frequencies estimated from the apparent source spectra do not show any dependence on hypocentral distance and magnitude, confirming that uncorrected anelastic attenuation effects do not significantly bias the results.

Rupture Characteristics and Bedrock Structural Control of the 2016 Mw 6.0 Intraplate Earthquake in the Petermann Ranges, Australia

ABSTRACT The 20 May 2016 surface-rupturing intraplate earthquake in the Petermann Ranges is the largest onshore earthquake to occur in the Australian continent in 19 yr. We use in situ and Interferometric Synthetic Aperture Radar surface observations, aftershock distribution, and the fitting of P-wave source spectra to determine source properties of the Petermann earthquake. Surface observations reveal a 21-km-long surface rupture trace (strike=294°±29°) with heterogeneous vertical displacements (<0.1–0.96 m). Aftershock arrays suggest a triangular-shaped rupture plane (dip ≈ 30°) that intersects the subsurface projection of the major geophysical structure (Woodroffe thrust [WT]) proximal to the preferred location of the mainshock hypocenter, suggesting the mainshock nucleated at a fault junction. Footwall seismicity includes apparent southwest-dipping Riedel-type alignments, including possible activation of the deep segment of the WT. We estimate a moment magnitude (Mw) of 6.0 and a corner frequency (fc) of 0.2 Hz, respectively, from spectral fitting of source spectra in the 0.02–2 Hz frequency band. These translate into a fault area of 124 km2 and an average slip of 0.36 m. The estimated stress drop of 2.2 MPa is low for an intraplate earthquake; we attribute this to low-frictional slip (effective coefficient of friction >0.015) along rupture-parallel phyllosilicate-rich surfaces within the host rock fabric with possible additional contributions from elevated pore-fluid pressures.

The characteristic of source spectra and stress drop of earthquakes in the Bucaramanga nest

<p>Spectra analysis is helpful to understand earthquake rupture processes and estimate source parameters like stress drop. Obtaining real source spectra and source time function isn’t easy, because the station recordings contain path effect and we usually can’t get precise path information. Empirical Green’s function (EGF) method is a popular way to cancel out the path effect, main two of which are the stacking spectra method (Prieto et al, 2006) and the spectral ratio method (Viegas et al, 2010; Imanishi et al, 2006). In our study, we apply the latter with multitaper spectral analysis method (Prieto et al, 2009) to calculate relative source spectra and relative source time function. Target event and EGFs must have similar focal mechanism and be collocated, so we combine correlation coefficient of wave at all stations and focal mechanism similarity to select proper EGFs.</p><p>The Bucaramanga nest has very high seismicity, so it’s suitable to calculate source spectra by using EGF method. We calculate the source spectra and source time function of about 1540 earthquakes (3-5.7ml, 135-160km depth) at Bucaramanga nest in Colombia. Simultaneously we also estimate corner frequency by fitting spectral source model (Brune, 1970; Boatwright, 1980) and stress drop using simple model (Eshelby, 1957) of earthquakes with multiple station recordings or EGFs. We obtain about 30000 events data with 12 stations from National Seismological Network of Colombia (RSNC).</p><p>The result show that the source spectra of most earthquakes fitted well by omega-square model are smooth, and the source spectra of some have obvious ‘holes’ near corner frequency, and the source time function of a few earthquakes appear two separate peeks. The first kind of earthquakes are style of self-arresting ruptures (Xu et al. 2015), which can be autonomously arrested by itself without any outside interference. Abercrombie (2014) and Wen et al. (2018) both researched the second kind of earthquakes and Wen think that this kind of earthquakes are style of the runaway ruptures including subshear and supershear ruptures. The last kind of earthquakes maybe be caused by simultaneous slip on two close rupture zone. Stress drop appear to slightly increase with depth and are very high (assuming rupture velocity/s wave velocity is 0.9). We also investigate the high-frequency falloff n, usually 2, of Brune model and Boatwright model by fitting all spectra, and find that the best value of n for Boatwright model is 2 and for Brune model is 3.5.</p>