Earthquake Science最新文献

High-quality datasets are critical for the development of advanced machine-learning algorithms in seismology. Here, we present an earthquake dataset based on the ChinArray Phase I records (X1). ChinArray Phase I was deployed in the southern north-south seismic zone (20° N–32° N, 95° E–110° E) in 2011–2013 using 355 portable broadband seismic stations. CREDIT-X1local, the first release of the ChinArray Reference Earthquake Dataset for Innovative Techniques (CREDIT), includes comprehensive information for the 105,455 local events that occurred in the southern north-south seismic zone during array observation, incorporating them into a single HDF5 file. Original 100-Hz sampled three-component waveforms are organized by event for stations within epicenter distances of 1,000 km, and records of ≥ 200 s are included for each waveform. Two types of phase labels are provided. The first includes manually picked labels for 5,999 events with magnitudes ≥ 2.0, providing 66,507 Pg, 42,310 Sg, 12,823 Pn, and 546 Sn phases. The second contains automatically labeled phases for 105,442 events with magnitudes of −1.6 to 7.6. These phases were picked using a recurrent neural network phase picker and screened using the corresponding travel time curves, resulting in 1,179,808 Pg, 884,281 Sg, 176,089 Pn, and 22,986 Sn phases. Additionally, first-motion polarities are included for 31,273 Pg phases. The event and station locations are provided, so that deep learning networks for both conventional phase picking and phase association can be trained and validated. The CREDIT-X1local dataset is the first million-scale dataset constructed from a dense seismic array, which is designed to support various multi-station deep-learning methods, high-precision focal mechanism inversion, and seismic tomography studies. Additionally, owing to the high seismicity in the southern north-south seismic zone in China, this dataset has great potential for future scientific discoveries.

Large earthquakes frequently occur along complex fault systems. Understanding seismic rupture and long-term fault evolution requires constraining the geometric and material properties of fault zone structures. We provide a comprehensive overview of recent advancements in seismological methods used to study fault zone structures, including seismic tomography, fault zone seismic wave analysis, and seismicity analysis. Observational conditions limit our current ability to fully characterize fault zones, for example, insufficient imaging resolution to discern small-scale anomalies, incomplete capture of crucial fault zone seismic waves, and limited precision in event location accuracy. Dense seismic arrays can overcome these limitations and enable more detailed investigations of fault zone structures. Moreover, we present new insights into the structure of the Anninghe-Xiaojiang fault zone in the southeastern margin of the Qinghai-Xizang Plateau based on data collected from a dense seismic array. We found that utilizing a dense seismic array can identify small-scale features within fault zones, aiding in the interpretation of fault zone geometry and material properties.

The Indo-Gangetic Plain (IGP) is one of the most seismically vulnerable areas due to its proximity to the Himalayas. Geographic information system (GIS)-based seismic characterization of the IGP was performed based on the degree of deformation and fractal dimension. The zone between the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT) in the Himalayan Mountain Range (HMR) experienced large variations in earthquake magnitude, which were identified by Number-Size (NS) fractal modeling. The central IGP zone experienced only moderate to low mainshock levels. Fractal analysis of earthquake epicenters reveals a large scattering of earthquake epicenters in the HMR and central IGP zones. Similarly, the fault fractal analysis identifies the HMR, central IGP, and south-western IGP zones as having more faults. Overall, the seismicity of the study region is strong in the central IGP, south-western IGP, and HMR zones, moderate in the western and southern IGP, and low in the northern, eastern, and south-eastern IGP zones.

We developed a modified stochastic finite-fault method for estimating strong ground motions. An adjustment to the dynamic corner frequency was introduced, which accounted for the effect of the location of the subfault relative to the hypocenter and rupture propagation direction, to account for the influence of the rupture propagation direction on the subfault dynamic corner frequency. By comparing the peak ground acceleration (PGA), pseudo-absolute response spectra acceleration (PSA, damping ratio of 5%), and duration, the results of the modified and existing methods were compared, demonstrating that our proposed adjustment to the dynamic corner frequency can accurately reflect the rupture directivity effect. We applied our modified method to simulate near-field strong motions within 150 km of the 2008 M_W7.9 Wenchuan earthquake rupture. Our modified method performed well over a broad period range, particularly at 0.04–4 s. The total deviations of the stochastic finite-fault method (EXSIM) and the modified EXSIM were 0.1676 and 0.1494, respectively. The modified method can effectively account for the influence of the rupture propagation direction and provide more realistic ground motion estimations for earthquake disaster mitigation.

The scientific goal of the Anninghe seismic array is to investigate the detailed geometry of the Anninghe fault and the velocity structure of the fault zone. This 2D seismic array is composed of 161 stations forming sub-rectangular geometry along the Anninghe fault, which covers 50 km and 150 km in the fault normal and strike directions, respectively, with ∼ 5 km intervals. The data were collected between June 2020 and June 2021, with some level of temporal gaps. Two types of instruments, i.e. QS-05A and SmartSolo, are used in this array. Data quality and examples of seismograms are provided in this paper. After the data protection period ends (expected in June 2024), researchers can request a dataset from the National Earthquake Science Data Center.

On August 6, 2023, a magnitude M_W5.5 earthquake struck Pingyuan County, Dezhou City, Shandong Province, China. This event was significant as no large earthquakes had been recorded in the region for over a century, and no active fault had been previously identified. This study collects 1309 P-wave arrival times and 866 S-wave arrival times from 74 seismic stations less than 200 km to the epicenter to constrain the spatial distribution of the mainshock and its 125 early aftershocks by the double difference earthquake relocation method, and selects 864 P-waveforms from 288 stations located within 800 km of the epicenter to constrain the focal mechanism solution of the mainshock through centroid moment tensor inversion. The relocation and the inversion indicate, the Pingyuan M_W5.5 earthquake was caused by a rupture on a buried fault, likely an extensive segment of the Gaotang fault. This buried fault exhibited a dip of approximately 75° to the northwest, with a strike of 222°, similar to the Gaotang fault. The rupture initiated at the depth of 18.6 km and propagated upward and northeastward. However, the ground surface was not broken. The total duration of the rupture was ∼6.0 s, releasing the scalar moment of 2.5895 × 10¹⁷ N·m, equivalent to M_W5.54. The moment rate reached the maximum only 1.4 seconds after the rupture initiation, and the 90% scalar moment was released in the first 4.6 s. In the first 1.4 seconds of the rupture process, the rupture velocity was estimated to be 2.6 km/s, slower than the local S-wave velocity. As the rupture neared its end, the rupture velocity decreased significantly. This study provides valuable insights into the seismic characteristics of the Pingyuan M_W5.5 earthquake, shedding light on the previously unidentified buried fault responsible for the seismic activity in the region. Understanding the behavior of such faults is crucial for assessing seismic hazards and enhancing earthquake preparedness in the future.

Wave propagation in horizontally layered media is a classical problem in seismic-wave theory. In semi-infinite space, a nondispersive Rayleigh wave mode exists, and the eigendisplacement decays exponentially with depth. In a layered model with increasing layer velocity, the phase velocity of the Rayleigh wave varies between the S-wave velocity of the bottom half-space and that of the classical Rayleigh wave propagated in a supposed half-space formed by the parameters of the top layer. If the phase velocity is the same as the P- or S-wave velocity of the layer, which is called the critical mode or critical phase velocity of surface waves, the general solution of the wave equation is not a homogeneous (expressed by trigonometric functions) or inhomogeneous (expressed by exponential functions) plane wave, but one whose amplitude changes linearly with depth (expressed by a linear function). Theories based on a general solution containing only trigonometric or exponential functions do not apply to the critical mode, owing to the singularity at the critical phase velocity. In this study, based on the classical framework of generalized reflection and transmission coefficients, the propagation of surface waves in horizontally layered media was studied by introducing a solution for the linear function at the critical phase velocity. Therefore, the eigenvalues and eigenfunctions of the critical mode can be calculated by solving a singular problem. The eigendisplacement characteristics associated with the critical phase velocity were investigated for different layered models. In contrast to the normal mode, the eigendisplacement associated with the critical phase velocity exhibits different characteristics. If the phase velocity is equal to the S-wave velocity in the bottom half-space, the eigendisplacement remains constant with increasing depth.

Seismic attenuation is a fundamental property of the Earth's media. Attenuation structure for the complicated geological structures with strong seismicity in the Sichuan-Yunnan region is poorly studied. In this study, we collected 108,399 waveforms of 11,517 local small earthquakes with magnitudes between 1.5 and 3.5 from January 2014 to September 2021 in the Sichuan-Yunnan region and its adjacent areas. We employed an envelope inversion technique for separating the intrinsic and scattering attenuations of the S coda wave, and obtained the intrinsic and scattering attenuation structures for frequencies between 0.25 and 8.00 Hz. The attenuation structures correlate well with the geological units, and some major faults mark the attenuation variations where historic large earthquakes have occurred. The regional average attenuation shows a negative frequency dependence. The average scattering attenuation has a faster descending rate than the average intrinsic attenuation, and is dominant at low frequencies, while at high frequencies the average intrinsic attenuation is stronger. The lateral variation in the intrinsic attenuation is consistent with the variation in heat flow, the scattering attenuation may be related to the scatter distribution and size. The total attenuation is consistent with the previous studies in this region, and the separate intrinsic and scattering attenuation may be useful in understanding regional tectonics and important in earthquake prevention and disaster reduction.

In 2023, two consecutive earthquakes exceeding a magnitude of 7 occurred in Türkiye, causing severe casualties and economic losses. The damage to critical urban infrastructure and building structures, including highways, railroads, and water supply pipelines, was particularly severe in areas where these structures intersected the seismogenic fault. Critical infrastructure projects that traverse active faults are susceptible to the influence of fault movement, pulse velocity, and ground motions. In this study, we used a unique approach to analyze the acceleration records obtained from the seismic station array (9 strong ground motion stations) located along the East Anatolian Fault (the seismogenic fault of the M_W7.8 mainshock of the 2023 Türkiye earthquake doublet). The acceleration records were filtered and integrated to obtain the velocity and displacement time histories. We used the results of an on-site investigation, jointly conducted by China Earthquake Administration and Türkiye’s AFAD, to analyze the distribution of PGA, PGV, and PGD recorded by the strong motion array of the East Anatolian Fault. We found that the maximum horizontal PGA in this earthquake was 3.0 g, and the maximum co-seismic surface displacement caused by the East Anatolian Fault rupture was 6.50 m. As the fault rupture propagated southwest, the velocity pulse caused by the directional effect of the rupture increased gradually, with the maximum PGA reaching 162.3 cm/s. We also discussed the seismic safety of critical infrastructure projects traversing active faults, using two case studies of water supply pipelines in Türkiye that were damaged by earthquakes. We used a three-dimensional finite element model of the PE (polyethylene) water pipeline at the Islahiye State Hospital and fault displacement observations obtained through on-site investigation to analyze pipeline failure mechanisms. We further investigated the effect of the fault-crossing angle on seismic safety of a pipeline, based on our analysis and the failure performance of the large-diameter Thames Water pipeline during the 1999 Kocaeli earthquake. The seismic method of buried pipelines crossing the fault was summarized.