Simulation of Pulse-Like Ground Motions during the 2015 Mw 8.3 Illapel Earthquake with a New Source Model Using Corrected Empirical Green’s Functions

Pulse-like near-source ground motions were observed by the local network during the 2015 Mw 8.3 Illapel, Chile earthquake. Such ground motions can be quite damaging to a wide range of infrastructures. The primary objective of this study is to provide a source model that can explain such ground motions. The isolated nature of the pulses indicated that the rupture of some small isolated region on the fault contributed to the generation of the pulse. Therefore, we considered such regions and termed them as Strong Motion Pulse Generation Areas (SPGAs). We used the corrected empirical Green's function (EGF) method because this method has been successfully applied to near-source pulse-like ground motions in previous studies. We simulated synthetic waveforms using the frequency dependent quality factor Q=239f0.71 and empirical site amplification factors, which we obtained by applying a generalized inversion technique to local weak-motion data. The result indicated that the observed ground motions from the Mw 8.3 Illapel earthquake can readily be explained with a source model that involves two SPGAs with dimensions of several kilometers in spite of the huge rupture zone of the earthquake. The source model can reproduce velocity waveforms, acceleration Fourier amplitude spectra (FAS) and pseudoacceleration response spectra. It also reproduces the duration of strong ground motions quite accurately. No significant bias was found with respect to distance and frequency. In conclusion, the corrected EGF method is a very efficient tool to simulate near-source ground motions of a subduction earthquake when it is combined with higher stress-drop subevents whose sizes are adjusted to the observed pulse widths.

Nonlinear Site Effects from the 30 November 2018 Anchorage, Alaska, Earthquake

ABSTRACT Anchorage, Alaska, is a natural laboratory for recording strong ground motions from a variety of earthquake sources. The city is situated in a tectonic region that includes the interface and intraslab earthquakes related to the subducting Pacific plate and crustal earthquakes from the upper North American plate. The generalized inversion technique was used with a local rock reference station to develop site response at >20 strong-motion stations in Anchorage. A database of 94 events recorded at these sites from 2005 to 2019 was also compiled and processed to compare their site response with those in the 2018 Mw 7.1 event (main event). The database is divided into three datasets, including 75 events prior to the main event, the main event, and 19 aftershocks. The stations were subdivided into the site classes defined in the National Earthquake Hazards Reduction Program based on estimated average shear-wave velocity in of the upper 30 m (VS30), and site-response results from the datasets were compared. Nonlinear site response was observed at class D and DE sites (VS30 of 215–300 and 150–215 m/s, respectively) but not at class CD and C sites (VS30 of 300–440 and 440–640 m/s, respectively). The relationship of peak ground acceleration versus peak ground velocity divided by VS30 (shear-strain proxy) was shown to further support the observation that sites with lower VS30 experienced nonlinear site response.

Stress-Drop Estimates for Induced Seismic Events in the Fort Worth Basin, Texas

ABSTRACT Earthquakes in the Fort Worth basin (FWB) have been induced by the disposal of recovered wastewater associated with extraction of unconventional gas since 2008. Four of the larger felt earthquakes, each on different faults, prompted deployment of local distance seismic stations and recordings from these four sequences are used to estimate the kinematic source characteristics. Source spectra and the associated source parameters, including corner frequency, seismic moment, and stress drop, are estimated using a modified generalized inversion technique (GIT). As an assessment of the validity of the modified GIT approach, corner frequencies and stress drops from the GIT are compared to estimates using the traditional empirical Green’s function (EGF) method for 14 target events. For these events, corner-frequency residuals (GIT−EGF) have a mean of −0.31 Hz, with a standard deviation of 1.30 Hz. We find consistent mean stress drops using the GIT and EGF methods, 9.56 and 11.50 MPa, respectively, for the common set of target events. The GIT mean stress drop for all 79 earthquakes is 5.33 MPa, similar to estimates for global intraplate earthquakes (1–10 MPa) as well as other estimates for induced earthquakes near the study area (1.7–9.5 MPa). Stress drops exhibit no spatial or temporal correlations or depth dependency. In addition, there are no time or space correlations between estimated FWB stress drops and modeled pore-pressure perturbations. We conclude that induced earthquakes in the FWB occurring on normal faults in the crystalline basement release pre-existing tectonic stresses and that stress drops on the four sequences targeted in this study do not directly reflect perturbations in pore-fluid pressure on the fault.

Determination of path attenuation and site characteristics of the North-west Himalayan region and adjoining regions within the Indian Territory using Generalized inversion method

North-west Himalayas and its adjoining regions have been experiencing deadly earthqaukes from time to time and are home for a large portion of population of Indian subcontinent. Knowledge of regional path attenuation and site parameters are prerequisite while attempting seismic hazard studies towards minimizing damages during future earthqaukes for a region. Present work focuses on the determination of path attenuation and site characteristics of earthqaukes recording stations, located in the north-west Himalayas and its adjoining regions, within India. It is done using two- step generalized inversion technique. In the first step of inversion, non-parametric attenuation curves are developed by constraining attenuation to be a smooth decaying function with hypocentral distance. Qs = (105 ± 11)f (0.94 ± 0.08) as S wave quality factor is obtained indicating that the region is seismically active having high degree of heterogeneities in the crustal medium. In the second step of generalized inversion, site amplification curve, at each recording station, is computed as the ratio of site spectral amplitude of horizontal and vertical components. In addition, based on Horizontal to vertical spectral ratio (HVSR) method, predominant frequency of each recording station is calculated. Values of predominant frequencies based on HVSR and generalized inversion are found matching for each of the recording station. Based on obtained predominant frequency, site class of 101 recording stations, which at present are absent, are determined in this work. Determined path attenuation as well as site parameters can be collectively used for developing regional ground motion models and subsequently for seismic hazard studies for the selected region.

Spectral Characteristics of Ground Motion from Induced Earthquakes in the Fort Worth Basin, Texas, Using the Generalized Inversion Technique

ABSTRACT A generalized inversion technique (GIT) is applied to local seismic data from 90 induced earthquakes (ML 2.0–3.9) in the Fort Worth Basin (FWB) of north Texas to separate path, site, and source characteristics and to improve local seismic hazard assessment. Seismograms from three earthquake sequences on spatially separated basement faults are recorded on 66 temporary stations. Because of the lack of hard-rock recording sites within the sedimentary basin, we developed a site correction method for the appropriate GIT process. At about 30 km distance from the hypocenters, we observed a change in spectral attenuation and thus focus data analysis within this distance range. The estimated quality factors for S and P waves result in a QS that is larger than QP which we interpret as a result of concentrations of crustal pore fluids or partial fluid-saturated material along the path; an interpretation consistent with fluid-rich sedimentary rocks in the FWB. Strong site amplifications as much as five times on horizontal components reflect the thick sediments in the basin. A limited number of sites exhibit amplification or deamplification on the vertical component that limits the use of horizontal-to-vertical spectral ratio methods for characterizing the site effect relative to the site effects estimated by GIT. Stress drops for all earthquakes range from 1.18 and 21.73 MPa with a mean of 4.46 MPa, similar to values reported for tectonic intraplate events. The stress-drop values suggest that strong motion and seismic hazard from the injection-induced earthquake in the FWB are comparable to those for tectonic earthquakes. The strong site amplification and fluid effects on propagation attenuation may be crucial factors to take into account for estimating seismic hazards of induced earthquakes in sedimentary basins.

Estimation of Source Spectra, Attenuation, and Site Responses from Strong-Motion Data Recorded in the 2019 Changning Earthquake Sequence

ABSTRACT This study used 306 accelerograms recorded at 22 strong-motion stations to investigate the source parameters, quality factor (Q), and site effects of S-wave Fourier acceleration amplitude spectrum (FAS) of the 2019 MS 6.0 Changning earthquake sequence in China with surface-wave magnitudes (MS) of 4.1–6.0. The generalized inversion technique (GIT) was adopted. The inverted stress drop of the mainshock was 1.15 MPa, and those of the aftershocks varied from 0.11 to 1.04 MPa with an average value of 0.43 MPa. The MS of these earthquakes were larger than Mw with an average magnitude difference of 0.22. The inverted Q values increase rapidly with frequencies at 0.5–4.0 Hz from 62 to 2920 and become less dependent at 4.0–25.0 Hz. Such a phenomenon indicates that the propagation path attenuation mechanism transited to intrinsic at high frequencies. A bilinear Q(f) model for which Q(f)=237.6f1.27 (Q<1280) and Q=1280 at higher frequencies was obtained. The high-frequency attenuation model of the study area was κ=0.0420+0.0001262R. The inverted site responses of the 22 stations were compared with those calculated using the horizontal-to-vertical spectral ratio (HVSR) method. In general, the amplification curves of most stations obtained with the GIT were similar to those of HVSR, and the amplification levels were relatively higher. Contrarily, obvious discrepancies were observed between the results estimated from the two methods at several stations. Such effects were attributed to the limitation that the majority of the stations were distributed along the boundary of the basin and mountainous areas, and the inverted Q values were representative of the specific area rather than the pure basin and mountainous areas. Finally, a nonlinear soil site effect was observed at 51GXT in earthquakes with peak ground acceleration greater than 300 cm/s2. The nonlinearity obviously aggravated the site amplification at 1.0–5.0 Hz.

Non-parametric spectral modelling of source parameters, path attenuation and site effects from broad-band waveforms of the Alborz earthquakes (2005–2017)

SUMMARY S-wave spectral amplitudes from 312 crustal earthquakes recorded at the Iranian National Broadband Seismic Network in the Alborz region between 2005 and 2017 are analysed in order to evaluate earthquake source parameters, path attenuation and site amplification functions using the non-parametric generalized inversion technique (GIT). We exploit a total number of 1117 seismograms with ML 3–5.6 in the frequency range 0.3–20 Hz. The evaluated non-parametric attenuation functions decay uniformly with distance for the entire frequency range and the estimated S-wave quality factor shows low Q values with relatively strong frequency dependence. We assume the omega-square source model to retrieve earthquake source parameters from the inverted source spectra. The obtained stress drops range from 0.02 to 16 MPa with a mean value of 1.1 MPa. Stress drop and radiated energy show fairly self-similar scaling with seismic moment over the available magnitude range; however, the magnitude range of this study is too narrow to draw a definite conclusion on source scaling characteristics. The obtained moment magnitude Mw and the local magnitude ML are linearly correlated and approximately equivalent in the range of Mw 3–4. For larger events, Mw generally underestimates ML by about 0.1–0.5 magnitude units. The estimated site amplification functions for horizontal component (GIT H) are nearly flat with no obvious pre-dominant frequency peaks for most stations, as expected for the sites of permanent broad-band seismic stations located on rock, though a few stations show amplification peaks from 1 to 8 Hz, with a maximum amplification of about a factor of 7 with respect to the reference site. The evaluated site responses for the vertical components present remarkable amplification or deamplification, leading to differences of the H/V amplitude levels in comparison with the GIT H amplification curves. The results of this study provide a valuable basis for predicting appropriate ground motions in a context of seismic hazard assessment.