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

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.

100 Years of Paper Seismograms from Denmark and Greenland, 1907–2008

The paper seismograms from 100 years of observations in Denmark and Greenland has since October 2021 been made available through the Danish National Archives. Five case stories illustrate the quality and variation of the seismograms, and the historical context of operation of the stations. (1) The earliest recorded earthquake in the archive is recorded at GDH station in Greenland, where the 1907 Mw 7.2 earthquake in Tajikistan is recorded on smoked paper. (2) The first Danish earthquake is a local event close to Copenhagen in 1930. (3) We have illustrated the 50 megaton nuclear explosion in Novaya Zemlya in 1961—the largest nuclear test explosion ever. (4) The M 9.2 earthquake in Alaska in 1964 recorded on several instruments at COP. (5) A local earthquake in northeast Greenland recorded both on paper on World-Wide Standard Seismographic Network instruments and digitally on a modern broadband instrument.

Cartesian Meshing Spherical Earth (CMSE): A Code Package to Incorporate the Spherical Earth in SPECFEM3D Cartesian Simulations

The SPECFEM3D_Cartesian code package is widely used in simulating seismic wave propagation on local and regional scales due to its computational efficiency compared with the one-chunk version of the SPECFEM3D_Globe code. In SPECFEM3D_Cartesian, the built-in meshing tool maps a spherically curved cube to a rectangular cube using the Universal Transverse Mercator projection (UTM). Meanwhile, the geodetic east, north, and up directions are assigned as the local x–y–z directions. This causes coordinate orientation issues in simulating waveform propagation in regions larger than 6° × 6° or near the Earth’s polar regions. In this study, we introduce a new code package, named Cartesian Meshing Spherical Earth (CMSE), that can accurately mesh the 3D geometry of the Earth’s surface under the Cartesian coordinate frame, while retaining the geodetic directions. To benchmark our new package, we calculate the residual amplitude of the CMSE synthetics with respect to the reference synthetics calculated by SPECFEM3D_Globe. In the regional scale simulations with an area of 1300 km × 1300 km, we find a maximum of 5% amplitude residual for the SPECFEM3D_Cartesian synthetics using the mesh generated by the CMSE, much smaller than the maximum amplitude residual of 100% for the synthetics based on its built-in meshing tool. Therefore, our new meshing tool CMSE overcomes the limitations of the internal mesher used by SPECFEM3D_Cartesian and can be used for more accurate waveform simulations in larger regions beyond one UTM zone. Furthermore, CMSE can deal with regions at the south and north poles that cannot be handled by the UTM projection. Although other external code packages can be used to mesh the curvature of the Earth, the advantage of the CMSE code is that it is open-source, easy to use, and fully integrated with SPECFEM3D_Cartesian.

Tectonic Context and Possible Triggering of the 2019–2020 Puerto Rico Earthquake Sequence

An historically unprecedented seismic moment was released by crustal events of the 2019–2020 earthquake sequence near southwest Puerto Rico. The sequence involved at least two, and perhaps three interacting fault systems. The largest Mw 6.4 event was likely triggered by left lateral strike-slip events along the eastern extension of the North Boquerón Bay-Punta Montalva fault zone. The mainshock occurred in a normal fault zone that extends into a region where previous studies documented extensional deformation, beyond the Ponce fault and the Bajo Tasmanian fault. Coulomb stress changes by the mainshock may have triggered further normal-faulting aftershocks, left lateral strike-slip events in the region where these two fault zones interacted, and possibly right lateral strike-slip aftershocks along a third structure extending southward, the Guayanilla fault zone. Extension directions of the seismic sequence are consistently north-northwest–south-southeast-oriented, in agreement with the Global Navigation Satellite Systems-inferred motion direction of eastern Hispaniola relative to western Puerto Rico, and with crustal stress estimates for the overriding plate boundary zone.

Active Faulting of Landforms along the Tuosuhu-Maoniushan Fault and Its Seismotectonic Implications in Eastern Qaidam Basin, China

The fold-and-thrust belt along the northern margin of the Qaidam basin is a typical active tectonic belt located in the northeast Tibetan Plateau. This belt is at a high risk of strong earthquakes with magnitudes larger than 6, as shown by multiple recorded events during 1962–2009. The lack of detailed late Quaternary surficial faulting data and systematic seismotectonic studies has posed difficulties in properly assessing the seismic risks and understanding the ongoing geodynamics in this region. In this study, we mapped the geomorphic features and fault traces from high-resolution satellite images and field investigations of the Tuosuhu-Maoniushan fault (TMF). Field photogrammetry was conducted to obtain deformation measurements using a DJI M300 real-time kinematic (RTK) drone. The TMF displaces the Holocene and late Pleistocene alluvial terraces in the eastern Qaidam basin. This fault dips to the south in the west and central segments (as a boundary of the Denan depression) and to the north in the eastern segment along the piedmont of the Maoniushan Mountains. The vertical slip rate is estimated to be 0.37 ± 0.08 mm/yr, which is similar to that of the active southern Zongwulongshan fault. By integrating our investigations with the previously published studies on deep structures and Cenozoic geology of the region, we propose a deep-seated thrust model for the seismotectonics of the northern margin of the Qaidam basin. The Aimunike, Tuosuhu-Maoniushan, southern Zongwulongshan, and Zongwulong faults, along with many folds, form an active compressional zone. The complex across-strike structures and along-strike segmentation could facilitate the release of strain through earthquakes of magnitude 6–7 in this broad seismotectonics belt, rather than through strong surface-rupturing events resulting from a single mature large fault.

Fluids Triggered the 2021 Mw 6.1 Yangbi Earthquake at an Unmapped Fault: Implications for the Tectonics at the Northern End of the Red River Fault

On 21 May 2021 a magnitude Mw 6.1 earthquake occurred in Yangbi region, Yunan, China, which was widely felt and caused heavy casualties. Imaging of the source region was conducted using our improved double-difference tomography method on the huge data set recorded by 107 temporary stations of ChinArray-I and 62 permanent stations. Pronounced structural heterogeneities across the rupture source region are discovered and locations of the hypocenters of the Yangbi earthquake sequence are significantly improved as the output of the inversion. The relocated Yangbi earthquake sequence is distributed at an unmapped fault that is almost parallel and adjacent (∼15 km distance) to the Tongdian–Weishan fault (TWF) at the northern end of the Red River fault zone. Our high-resolution 3D velocity models show significant high-velocity and low-VP/VS ratios in the upper crust of the rupture zone, suggesting the existence of an asperity for the event. More importantly, low-VS and high-VP/VS anomalies below 10 km depth are imaged underlying the source region, indicating the existence of fluids and potential melts at those depths. Upward migration of the fluids and potential melts into the rupture zone could have weakened the locked asperity and triggered the occurrence of the Yangbi earthquake. The triggering effect by upflow fluids could explain why the Yangbi earthquake did not occur at the adjacent TWF where a high-stress accumulation was expected. We speculate that the fluids and potential melts in the mid-to-lower crust might have originated either from crustal channel flow from the southeast Tibet or from local upwelling related to subduction of the Indian slab to the west.

Time-Lapse Monitoring of Daily Velocity Changes in Binchuan, Southwestern China, Using Large-Volume Air-Gun Source Array Data

Temporal changes of seismic velocities in the Earth’s crust can be induced by stress perturbations or material damage from reasons such as strong ground motion, volcanic activities, and atmospheric effects. However, monitoring the temporal changes remains challenging, because most of them generally exist in small travel-time differences of seismic data. Here, we present an excellent case of daily variations of the subsurface structure detected using a large-volume air-gun source array of one-month experiment in Binchuan, Yunnan, southwestern China. The seismic data were recorded by 12 stations within ∼10 km away from the source and used to detect velocity change in the crust using the deconvolution method and sliding window cross-correlation method, which can eliminate the “intercept” error when cutting the air-gun signals and get the real subsurface variations. Furthermore, the multichannel singular spectral analysis method is used to separate the daily change (∼1 cycle per day) from the “long-period” change (<1 cycle per day) or noise. The result suggests that the daily velocity changes at the two nearest stations, 53277 (offset ∼700 m) and 53278 (offset ∼2.3 km), are well correlated with air temperature variation with a time lag of 5.0 ± 1.5 hr, which reflects that the velocity variations at the subsurface are likely attributed to thermoelastic strain. In contrast, both daily and long-period velocity changes at distant stations correlate better with the varying air pressure than the temperature, indicating that the velocity variations at deeper depth are dominated by the elastic loading of air pressure. Our results demonstrate that the air-gun source is a powerful tool to detect the velocity variation of the shallow crust media.

Reply to “Comment on ‘High-Definition Mapping of the Gutenberg–Richter b-Value and Its Relevance: A Case Study in Italy’ by M. Taroni, J. Zhuang, and W. Marzocchi” by Laura Gulia, Paolo Gasperini, and Stefan Wiemer

Taroni et al. (2021; hereafter TZM21) proposed a method to perform a spatial b-value mapping based on the weighted-likelihood estimation and applied this method to the Italian region as a tutorial example. In the accompanying comment, Gulia et al. (2021; hereafter GGW21) did not challenge the TZM21’s method, but they argued that the catalog used by TZM21 is contaminated by quarry blasts, introducing a bias that may impact any seismotectonic or hazard interpretations. Although in TZM21 the application to the Italian territory was only a tutorial example and we purposely did not make any thorough discussion on the meaning of the results in terms of seismotectonic or seismic hazards (that would have required many more analyses), we acknowledge the potential role of the quarry blasts, and we add some further analysis here. We thank GGW21 for giving us this opportunity. Here, removing the part of the catalog contaminated by quarry blasts and applying the same analysis as in TZM21, we obtain results that are very similar to the ones reported in TZM21; specifically, only one region that is characterized by low natural seismicity rate shows a marked effect of the quarry blasts on the b-value.

FrosPy: A Modular Python Toolbox for Normal Mode Seismology

We present free oscillations Python (FrosPy), a modular Python toolbox for normal mode seismology, incorporating several Python core classes that can easily be used and be included in larger Python programs. FrosPy is freely available and open source online. It provides tools to facilitate pre- and postprocessing of seismic normal mode spectra, including editing large time series and plotting spectra in the frequency domain. It also contains a comprehensive database of center frequencies and quality factor (Q) values based on 1D reference model preliminary reference Earth model for all normal modes up to 10 mHz and a collection of published measurements of center frequencies, Q values, and splitting function (or structure) coefficients. FrosPy provides the tools to visualize and convert different formats of splitting function coefficients and plot these as maps. By giving the means of using and comparing normal mode spectra and splitting function measurements, FrosPy also aims to encourage seismologists and geophysicists to learn about normal mode seismology and the study of the Earth’s free oscillation spectra and to incorporate them into their own research or use them for educational purposes.

Using LSTM Neural Networks for Onsite Earthquake Early Warning

Onsite earthquake early warning (EEW) systems determine possible destructive S waves solely from initial P waves and issue alarms before heavy shaking begins. Onsite EEW plays a crucial role in filling in the blank of the blind zone near the epicenter, which often suffers the most from disastrous ground shaking. Previous studies suggest that the peak P-wave displacement amplitude (Pd) may serve as a possible indicator of destructive earthquakes. However, the attempt to use a single indicator with fixed thresholds suffers from inevitable errors because the diversity in travel paths and site effects for different stations introduces complex nonlinearities. In addition, the short warning time poses a threat to the validity of EEW. To conquer the aforementioned problems, this study presents a deep learning approach employing long short-term memory (LSTM) neural networks, which can produce a highly nonlinear neural network and derive an alert probability at every time step. The proposed LSTM neural network is then tested with two major earthquake events and one moderate earthquake event that occurred recently in Taiwan, yielding the results of a missed alarm rate of 0% and a false alarm rate of 2.01%. This study demonstrates promising outcomes in both missed alarms and false alarms reduction. Moreover, the proposed model provides an adequate warning time for emergency response.