Emergence of Low-Frequency Aftershocks of the 2019 Ridgecrest Earthquake Sequence

ABSTRACT Low-frequency earthquakes (LFEs) generally have relatively stronger spectral components in the lower frequency range compared with what is expected for regular earthquakes based on their magnitude. LFEs generally occur in volcanic systems or deep (>∼15 km) in plate boundary fault zones; however, LFEs have also been observed in nonvolcanic, upper crustal settings. Because there are few studies that explore the spatiotemporal behaviors of LFEs in the shallow crust, it remains unclear whether the shallow-crustal LFEs reflect local attenuation in their immediate vicinity or differences in their source mechanism. Therefore, it is important to identify shallow-crustal LFEs and to characterize their spatiotemporal activity, which may also improve our understanding of LFEs. In this study, we focus on detecting shallow-crustal LFEs and explore the possible generation mechanisms. We analyze 29,646 aftershocks in the 2019 Ridgecrest, California, earthquake sequence, by measuring the frequency index (FI) to identify candidate low-frequency aftershocks (LFAs), while accounting for the magnitude dependency of the FI. Using small earthquakes (ML 1–3) recorded in the borehole stations to minimize the attenuation effects in near-surface layers, we identify 68 clear LFAs in total. Based on their distribution and comparisons with other seismic parameters measured by Trugman (2020), the LFAs possess distinct features from regular events in the same depths range, including low corner frequencies and low stress drops. Events in the close vicinity of LFAs exhibit lower average FI values than regular aftershocks, particularly if the hypocentral distance between an LFA and its neighbors is less than 1 km. Our results suggest that LFAs are related to local heterogeneity or a highly fractured fault zone correlated with an abundance of cross faults induced by the aftershock sequence at shallow depths. Zones of high pore-fluid pressure in intensely fractured fault zones could cause the bandlimited nature of LFAs and LFEs in general.

Vibration Effect Induced by Rock Breaking Technology Based on Dry Ice and Energy-gathered Agent in Trench Excavation

Rock breaking technology based on dry ice and energy-gathered agent has been developed and successfully applied in trench excavation for construction of oil pipeline. The vibration velocity waveform induced by this technology was monitored in site test to determine the attenuation law of vibration velocity with hypocentral distance. The results show that this rock breaking technology is effective method of trench excavation. It does not excessively damage the adjacent rock mass, ensuring the integrity of ditch walls. The vibration velocity induced by this technology is decay with the increase of hypocentral distance. At the hypocentral distance of 10m, the vibration velocity reduces to less than 20mm/s, which meets the requirements of the safety standard of blasting vibration in general buildings engineering. The results of this experiment have an important guiding effect on the field engineering practice and application of rock breaking technology based on dry ice and energy-gathered agent.

Recalibration of the Local Magnitude (ML) Scale for Earthquakes in the Yellowstone Volcanic Region

ABSTRACT The Yellowstone volcanic region is one of the most seismically active areas in the western United States. Assigning magnitudes (M) to Yellowstone earthquakes is a critical component of monitoring this geologically dynamic zone. The University of Utah Seismograph Stations (UUSS) has assigned M to 46,767 earthquakes in Yellowstone that occurred between 1 January 1984 and 31 December 2020. Here, we recalibrate the local magnitude (ML) distance and station corrections for the Yellowstone volcanic region. This revision takes advantage of the large catalog of earthquakes and an increase in broadband stations installed by the UUSS since the last ML update in 2007. Using a nonparametric method, we invert 7728 high-quality, analyst-reviewed amplitude measurements from 1383 spatially distributed earthquakes for 39 distance corrections and 20 station corrections. The inversion is constrained with four moment magnitude (Mw) values determined from time-domain inversion of regional-distance broadband waveforms by the UUSS. Overall, the new distance corrections indicate relatively high attenuation of amplitudes with distance. The distance corrections decrease with hypocentral distance from 3 km to a local minimum at 80 km, rise to a broad peak at 110 km, and then decrease again out to 180 km. The broad peak may result from superposition of direct arrivals with near-critical Moho reflections. Our ML inversion doubles the number of stations with ML corrections in and near the Yellowstone volcanic region. We estimate that the additional station corrections will nearly triple the number of Yellowstone earthquakes that can be assigned an ML. The new ML distance and station corrections will also reduce uncertainties in the mean MLs for Yellowstone earthquakes. The new MLs are ∼0.07 (±0.18) magnitude units smaller than the previous MLs and have better agreement with 12 Mws (3.15–4.49) determined by the UUSS and Saint Louis University.

New Magnitude Proxy for Earthquake Early Warning Based on Initial Time Series and Frequency

A new characteristic parameter Sdτ is proposed to improve the performance of magnitude estimation in earthquake early warning (EEW). Sdτ is the product of summation of absolute displacement multiplied by the maximum predominant period (τmaxP) for the first arriving seconds of a seismic wave. About 30,725 underground records at borehole stations for 3645 earthquakes with magnitude between 4.0 and 9.0 from the Japanese KiK-net were used to compare the magnitude proxy performance based on the proposed Sdτ with that based on either τmaxP or peak displacement Pd. The comparison results show that for a magnitude between 4.0 and 7.3, Sdτ has a better correlation with magnitude and higher estimated accuracy than either τmaxP or Pd. Hypocentral distance is not required when using Sdτ, but it can be used to further improve the accuracy of magnitude estimate. These results confirm that Sdτ can significantly improve the accuracy and timeliness of continuous magnitude estimation in an EEW system.

New Ground Motion to Intensity Conversion Equations for China

In this study, we use the strong motion records and seismic intensity data from 11 moderate-to-strong earthquakes in the mainland of China since 2008 to develop new conversion equations between seismic intensity and peak ground motion parameters. Based on the analysis of the distribution of the dataset, the reversible conversion relationships between modified Mercalli intensity (MMI) and peak ground acceleration (PGA), peak ground velocity (PGV), and pseudo-spectral acceleration (PSA) at natural vibration periods of 0.3 s, 1.0 s, 2.0 s, and 3.0 s are obtained by using the orthogonal regression. The influence of moment magnitude, hypocentral distance, and hypocentral depth on the residuals of conversion equations is also explored. To account for and eliminate the trends in the residuals, we introduce a magnitude-distance-depth correction term and obtain the improved relationships. Furthermore, we compare the results of this study with previously published works and analyze the regional dependence of conversion equations. To quantify the regional variations, a regional correction factor for China, suitable for adjustment of global relationships, has also been estimated.

Stress Drop Derived from Spectral Analysis Considering the Hypocentral Depth in the Attenuation Model: Application to the Ridgecrest Region, California

ABSTRACT We investigate the impact of considering a depth-dependent attenuation model on source parameters assessed through a spectral decomposition. In particular, we evaluate the effect of considering the hypocentral depth as an additional variable for the attenuation model, using as the target the tendency of the average stress drop to increase with depth, as observed in recent studies. We analyze the Fourier spectra of S-wave windows for about 1900 earthquakes with a magnitude above 2.5 recorded in the Ridgecrest region, southern California. Two different parameterizations of the attenuation term are implemented in the spectral decomposition, either as a function of the hypocentral distance alone or as a function of both epicentral distance and depth. The comparison of the spectral attenuation curves shows that, although the hypocentral model describes, on average, the range of values spanned by the attenuation curve for different depths, systematic differences with distance, depth, and frequency are observed. These differences are transferred to the source spectra and, in turn, to the source parameters extracted from the best-fitting ω−2 models. In particular, stress drops for events deeper than 7 km are, on average, almost double even when depth is introduced explicitly in the attenuation model. The increase of stress drop with depth is confirmed also after accounting for the increase of the shear velocity with depth, which absorbs about 30%–40% of the total increase. Moreover, a qualitative comparison with a model for the gradient of the effective normal stress confirms the reliability of the observed trend. Finally, the coherent spatial patterns shown by a simplified 2D tomographic representation of the spectral residuals highlights the impact on ground-shaking variability of the lateral variability of the crustal attenuation properties in the region.

Shear-Wave Velocity Site Characterization in Oklahoma from Joint Inversion of Multimethod Surface Seismic Measurements: Implications for Central U.S. Ground-Motion Prediction

ABSTRACT We analyze multimethod shear (SH)-wave velocity (VS) site characterization data acquired at three permanent and 25 temporary seismograph stations in Oklahoma that recorded M 4+ earthquakes within a 50 km hypocentral distance of at least one of the 2016 M 5.1 Fairview, M 5.8 Pawnee, or M 5.0 Cushing earthquakes to better constrain earthquake ground-motion modeling in the region. We acquired active-source seismic data for time-averaged VS to 30 m depth (VS30) at 28 seismograph stations near the Fairview, Pawnee, and Cushing epicentral areas. The SH-wave refraction travel times coupled with Rayleigh- and Love-wave phase velocity dispersion were extracted and modeled in a nonlinear least-squares (L2) joint inversion to obtain a best-fit 1D VS versus depth profile for each site. At a subset of sites where the preferred L2 inverse model did not optimally fit each of the Love, Rayleigh, and SH travel-time datasets, we explore application of simulated annealing in a joint inversion to find a more global solution. VS30 values range from 262 to 807 m/s for the preferred measured (in situ) VS profiles, or National Earthquake Hazards Reduction Program (NEHRP) site class D to B, and are broadly comparable with estimates from previous data reports in the region. Site amplification estimates were calculated next from 1D SH transfer functions of the preferred VS profiles and then compared against observed horizontal-to-vertical spectral ratios (HVSRs) from nearby seismograph stations. We generally see good agreement between the predicted in situ model and the observed HVSR resonant frequencies, with nominal amplifications between 2 and 10 within the 2–15 Hz frequency band. Next, using 40 known in situ VS30 measurements in the region, we demonstrate that the in situ VS30 values improve the fit for selected suites of ground-motion models (GMMs) for M 4+ earthquakes within a 50 km hypocentral distance when compared with proxy methods, arguing for future development of GMMs implementing in situ VS profiles.

Station Correction of P-Alert Network to Improve Magnitude Estimation for Earthquake Early Warning

<p>Magnitude estimation for earthquake early warning has been shown that it can be achieved by utilizing the relationship among the first three seconds P-wave amplitude, hypocentral distance and magnitude. However, the regression models in previous studies about P-Alert didn't include station correction factors, which may cause non-negligible effects. Thus, to improve the precision of magnitude estimation, we take station corrections into consideration when building the regression model. For the reason that station corrections are the unobserved latent variables of the model, we adopt the iteration regression method, which is based on the expectation-maximization algorithm, to determine them. By using this method, we are able to approach the values of both the station corrections and the coefficients of the regression model after several iterations. Our preliminary results show that after utilizing the iteration regression method, the standard deviation reduces from 0.30 to 0.26, and the station corrections we get range from -0.70 to 0.66.</p>

Energetic characteristics of High Frequency (HF) and Very High Frequency (VF) Martian events

<p>The InSight seismometer SEIS recorded tens of high-frequency (1.5-5Hz; HF) and Very-high frequency (1.5-15Hz, VF) Martian events. They are characterized by two temporally separated arrivals with a gradual beginning, a broad maximum and a very long decay. This observation is consistent with a long-range propagation of seismic P and S waves in a heterogeneous crust (Van Driel et al., accepted). To examine this hypothesis, first, we employ basic multiple-scattering concepts on the two groups of events. Then, we propose a full envelope modeling based on elastic radiative transport in a half-space. The model parametrization and the radiative transfer equations are presented in (Lognonné, P., et al. (2020) and Margerin, L., (2017)). We find that both HF and VF signals are depolarized and verify Gaussian statistics, at the exception of the ballistic primary and secondary arrivals. These properties agree with a multiple-scattering origin. For VF events, the energy partitioning ratio V<sup>2</sup>/H<sup>2</sup> between horizontal and vertical components is frequency dependent. We observe that V<sup>2</sup>/H<sup>2 </sup>is maximum at the so-called ‘2.4Hz resonance’ (~2) and decreases rapidly at frequencies higher than 5Hz (~0.1) then i remains relatively low up to frequencies of 15Hz at least. HF events do not exhibit a decrease of V<sup>2</sup>/H<sup>2 </sup>at high frequencies however further analysis reveals a strong correlation between energy partitioning and signal-to-noise (S/N) ratio for HF events. This observation suggests that a part of the difference between the HF and VF events can to some extent be explained by noise contamination. The generally low V<sup>2</sup>/H<sup>2 </sup>ratio of VF events is reminiscent of the response of unconsolidated layers, as observed at Pinyon Flats Observatory on Earth (Margerin, L., et al. (2009)). Unlike earthquakes and moonquakes observed in the same frequency band, the delay time measured from onset to peak of the secondary arrival of HF and VF events is frequency-independent. This suggests that the spectrum of heterogeneity of the Martian crust is smooth. We observe that, for HF and VF events, the delay time is weakly dependent on hypocentral distance. This observation cannot be reconciled with the predictions of multiple-scattering theories in a statistically homogeneous medium however it suggests a stratification of heterogeneity in the Martian lithosphere. The coda quality factor Q<sub>c</sub> of VF events is high and shows a linear increase with frequency. Q<sub>c</sub> of HF events is higher but it may be overestimated due to the noise contamination. The linear frequency dependence of Q<sub>c</sub> is strongly reminiscent of the leakage effect in a crustal scattering waveguide and suggests that part of the observed coda attenuation may be of structural origin. The full envelope modeling of the S0334a VF event results shows that the estimated value of the diffusivity (≃ 619 km<sup>2</sup>/s) is almost 6 times greater than for the S0128a VF event (≃ 90 km<sup>2</sup>/s). This observation again suggests a stratification of heterogeneity. In future works, we will perform the full envelope modeling of all the VF selected events at different frequencies to constrain a 1D attenuation and diffusion model of the Martian crust.</p>

Is an Earthquake Felt Inside a Car?

The analysis of how an earthquake is felt was addressed with the data provided by citizens through a website dedicated to the perception of earthquakes in Italy (Data and Resources). The analysis focused on the perception of earthquakes by observers inside both parked and moving cars. These situations were compared with outdoor ones. The felt percentage of each situation was quantified for epicentral distance ranges and European Macroseismic Scale (EMS) degree. One of the main findings was the greatest sensitivity to shaking for people inside parked cars due to resonance phenomena of the automobile–observer system. The distribution of the intensity of perception in the car was analyzed as a function of the hypocentral distance and the magnitude of the earthquake. It was possible to define the attenuation trends of these intensities. The comparison of these trends with those of the equations for estimation of response spectral ordinates allowed us to have an evaluation of the frequency values of the seismic waves that caused the vibrations felt, which were found to agree with the typical frequencies of the car–observer system, as highlighted by independent studies. The results of this analysis show the possibility to include the perception of the earthquake inside a parked and moving car among the diagnostics used in the definition of macroseismic intensity degree of the EMS.