Next Article in Journal
Study of Ionospheric Bending Angle and Scintillation Profiles Derived by GNSS Radio-Occultation with MetOp-A Satellite
Previous Article in Journal
Coupled Coordination Analysis between Urbanization and Eco-Environment in Ecologically Fragile Areas: A Case Study of Northwestern Sichuan, Southwest China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

An Aftershock Deletion Method Based on Fault Buffer Zone

1
Key Laboratory of Earthquake Geodesy, Institute of Seismology, China Earthquake Administration, Wuhan 430071, China
2
Institute of Disaster Prevention, Langfang 065201, China
3
Zachry Department of Civil and Environmental Engineering, Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(6), 1662; https://doi.org/10.3390/rs15061662
Submission received: 7 January 2023 / Revised: 16 March 2023 / Accepted: 16 March 2023 / Published: 19 March 2023
(This article belongs to the Section Earth Observation Data)

Abstract

:
The existence of aftershocks in an earthquake sequence can impact the analysis of the mainshock. In this study, we present a method for deleting an aftershock sequence based on the spatial relationship between earthquakes and faults. This method improves the performance of space window selection in the classical K-K method by eliminating aftershocks with an ideal fault buffer zone. The determination of fault buffer zones is based on a trial-and-error analysis of 69,714 earthquake records from the China Seismic Network Center (CENC) collected between 1980 and 2020. We selected 20 typical big earthquakes ( M L 7.0–8.0 or ~ M s 6.6–8.0; for earthquakes above magnitude M s 7 or M L 7.2 , M L is approximately equal to M s ) as the mainshocks to establish the fault buffer zones. We also propose an empirical formula to determine the distance of the fault buffer zone by counting the aftershock deletion effect at different buffer distances. Compared with the classical K-K method, our method considers the correlation between the spatial distribution of aftershocks and faults, eliminates earthquake groups that are not related to the mainshock, greatly reduces the spatial range of aftershocks, improves the performance of deleting aftershocks of different magnitudes, and provides a new rule and reference for aftershock deletion.

1. Introduction

Earthquakes are one of the natural disasters that have a significant impact on life and property [1,2] and are often viewed as independent events. Earthquakes’ statistical characteristics have always been an important tool for researchers to predict earthquakes and evaluate earthquake risk. Aftershocks are responses to stress changes produced by large earthquakes and represent the most common phenomenon triggered by earthquakes [3]. Aftershock analysis is important for the probabilistic seismic hazard analysis (PSHA) [4,5], as it involves a series of earthquakes that occur after the main earthquake. After the mainshock, the stress distribution of the fault changes significantly, resulting in a series of aftershocks with the process of stress readjustment. According to Shilen and Toksoz [6], aftershocks are entirely inconsistent with the Poisson process, and they can obscure the statistical characteristics of an earthquake as an independent event to some extent [7], such as the unknown recurrence period [8]. Although the impact of aftershocks on the recurrence cycle is small, it cannot be ignored when evaluating local earthquake risk [6]. Therefore, deleting aftershocks is crucial for earthquake statistical analysis and earthquake prediction.
Numerous studies have been conducted to develop effective and efficient methods for removing aftershocks. Knopoff and Gardner [9] deleted aftershocks by eliminating all earthquake sequences within a particular distance and time period of a large earthquake event within a square area. Shlien and Toksdz [10] considered the spatial distance, time interval, and regional seismic activity to determine the mainshock-aftershock relationship between two earthquake events. Console et al. [11] used an empirical formula that factors in the spatial distance and magnitudes of two earthquake events to remove aftershocks. The K-K method, proposed by Keilis-Borok and Knopoff L et al. [12], adjusts the size of the spatial and time window used to remove aftershocks based on the magnitude of the mainshock. Davis and Cliff [13] presented a unified relationship between earthquakes and the parameters of the space-time window, while Chen et al. [8] determined aftershocks by considering the length of the fault and the time window specified by Console et al. [11]. In recent years, some other methods have also been proposed, such as stochastic decluttering based on point processes [14], nearest-neighbor distances of events in the space-time-energy domain [15,16,17], narrow spatial aftershock zones [18], and nonparametric network-tree aftershock identification [19]. Most of these methods are based on temporal and spatial relations and can quickly remove aftershocks from the entire earthquake catalog.
The methods mentioned above typically use circular or square areas centered on the main earthquake’s epicenter to determine the spatial range of aftershocks, without considering the fault characteristics closely related to earthquakes. However, aftershocks usually occur on the same fracture plane as the mainshock, and stress changes resulting from Coulomb stress transfer can also trigger aftershocks in other fault zones. Stress drops and normalized rupture width most strongly correlate with aftershock productivity [20,21,22,23]. Nevertheless, most earthquakes are concentrated around the fault zone near the mainshock’s epicenter. It is important to note that earthquakes are the result of tectonic movements, which release the accumulated high stress in the crust through fault movements. Faults are a macroscopic feature of tectonic movements, and the occurrence of earthquakes is highly correlated with the fault scale, including secondary faults, faults, plane width, and fault types, especially its length. The fault length of an active fault is positively correlated with the stress accumulation in the crust.
The spatial location of earthquakes is closely related to fault zones. To reduce the spatial range of aftershocks selected by the K-K method, a new quick method called the “fault buffer zone” is proposed. This method is mainly based on the positive correlation between the crustal stress and the fault length. As a quick method, it temporarily ignores other attributes of the fault, such as fault dislocation mode and fault plane, which require careful geological verification. The basic idea is to select the closest related fault of a specific mainshock as the main fault and use the line elements determined by the main fault to extend a certain distance outward to determine the buffer zone of the main fault, which is the aftershock space. Fault length is a comprehensive response to geological structure and stress accumulation. Therefore, this method considers the correlation between aftershocks, main faults, and differences in geological tectonic environments in different regions.
The introductory section of this paper provides an overview of the fault buffer zone method for deleting aftershocks. The subsequent section describes the methodology for determining the main parameters of the fault buffer zone method. The third section presents the empirical formula for calculating the fault buffer zone, which is based on a statistical analysis of earthquakes in the Chinese mainland since 1980 and tests on major earthquakes. The final section summarizes the main conclusions of the study.

2. Aftershock Deletion Algorithm Based on “Fault Buffer Zone”

The K-K algorithm [12] employs an empirical method to remove aftershocks by defining a spatial and temporal based on the magnitude of the mainshock. The mainshock magnitude is denoted as M , the aftershock magnitude as m , the distance between two earthquakes as r , the time interval as t , the spatial window distance as R 0 , and the time window interval as T 0 . If the following condition is met, then m is classified as the aftershock of M :
r R 0 M ; t T 0 M ; m < M ,
The empirical value of the aftershock space and time window of the K-K algorithm is presented in Table 1.
The K-K method’s spatial window is typically a circle centered on the epicenter. In this study, we replace the circle with a “fault buffer zone” that considers the correlation between earthquakes and their corresponding seismogenic faults. The fault buffer zone is a closed polygon that extends a certain distance outward from the seismogenic fault. Figure 1 shows the flowchart of the aftershock deletion algorithm based on a fault buffer zone.
Figure 1 illustrates the aftershock deletion algorithm based on the fault buffer zone. The seismogenic fault determined by geological survey or focal mechanism solution and closest to the epicenter of the mainshock is considered the target active fault and is simplified into a linear unit. Active faults are those that have been active since the late Quaternary and may still be active today or be reactivated in the future. The proposed method for aftershock deletion only requires the identification of faults and the determination of their existence through the analysis and interpretation of aerial and satellite photos, geophysical maps, and geological maps, without considering fault type. This is because the aftershocks are distributed in the buffer zone, which is established with the seismogenic fault as the center, regardless of the fault type, including normal fault, reverse fault, and strike-slip fault.
The fault buffer zone replaces the circular spatial window of the K-K method by extending a certain distance outward from the seismogenic fault, which is simplified into a linear unit, to form a closed polygon.
This aftershock deletion algorithm based on the fault buffer zone requires determining the following parameters for the earthquake (mainshock) being studied:
(1)
The seismogenic fault closest to the mainshock was identified through a geological survey or focal mechanism solution, simplified into a linear unit.
(2)
The length of the seismogenic fault.
(3)
The buffer distance of the fault buffer zone.

2.1. Determining the Seismogenic Fault and Its Length Related to the Mainshock

Earthquakes may be associated with a single fault, multiple faults, multisegmented faults, or concealed faults [24]. To identify the seismogenic fault related to the mainshock, postearthquake investigation and focal mechanism solutions can be used. If the earthquake is related to multiple faults, the faults can be simplified into a linear unit. The next step is to determine the length of the seismogenic fault. The empirical Formula (2) by Guo et al. [25] shows the relationship between earthquake magnitude (surface wave magnitude) and fault length, while Wells and Coppersmith [26] proposed a formula relating earthquake moment magnitude ( M w ) and fault rupture length. However, China’s earthquake catalog uses near earthquake magnitude to define the magnitude ( M L ) in the early stage. In this paper, we still use the calculation formula of M L and fracture length, because China earthquakes catalog with M magnitude by default and will introduce deviation if historical earthquakes are converted to Mw magnitude and it is not possible to use Mw for small earthquakes. For earthquakes above magnitude M s 7.0 or M L 7.2 , M L is approximately equal to M s . In this formula, M represents the earthquake magnitude, and L represents the length of the fault in kilometers. There is a significant positive correlation between the logarithm of earthquake magnitude and fault length.
M = 3.3 + 2.1 log L ,

2.2. Determination of Fault Buffer Distance

Based on the aforementioned methodology, the fault buffer zone is determined by taking the seismogenic fault as the center and extending a certain distance, which is determined by the Joyner-Boore distance. The Joyner-Boore distance refers to the shortest distance from an observation point on the ground to the projection of the fault on the ground. The size of the buffer zone is determined by the magnitude of the mainshock. A larger magnitude corresponds to a larger Joyner-Boore distance and a larger buffer zone.
There are 69,714 earthquake records of magnitude greater than M L 3.0 collected from the China Seismic Network Center (CENC) between 1980 and 2020. While strong earthquakes generally receive great attention and provide rich observation data, such as seismogenic faults, focal mechanism solutions, and geological observations, it is important to find an appropriate fault buffer distance for earthquakes of different magnitudes. To accomplish this, we selected almost all 20 typical strong earthquakes with magnitudes ranging from M L 7.0– M L 8.0( ~ M s 6.6–8.0) as the main shocks and carried out trial-and-error analysis. The time and space windows provided by the K-K method were used as a reference to determine the initial value of the spatial location and length of the main faults.
To determine the appropriate buffer distance for earthquakes of different magnitudes, we adjusted the buffer distance, i.e., the Joyner-Boore distance, according to the magnitude of the earthquake (Figure 2), with the initial value set to the space window distance R 0 of the K-K method. Proper adjustment of the buffer distance is necessary to effectively handle the aftershocks of strong earthquakes.

3. Statistical Results and Discussion

The trial-and-error analysis of aftershock deletion is carried out using earthquake catalog data (69,714 events greater than M L 3.0) collected by CENC from 1 January 1980 to 23 December 2020. Among them, there are 2772 earthquakes with magnitudes greater than or equal to 5.0 and 461 earthquakes with magnitudes of 6.0 and above. The fault data used in this analysis is sourced from the active tectonic map of China (1:4 million) developed by Academician Deng Qidong [27] and provided by the National Seismological Science Data Center. The fault data includes the geographical location and fault attributes of major faults in China, such as fault name, length, strike, dip, dip angle, fault nature, and the latest active age.
In Section 2, we used the fault buffer zone method to analyze the typical aftershocks of M s 7.0~8.0 earthquakes. We obtained appropriate buffer distances for 20 great earthquakes with magnitudes between 7.0 and 8.0. These earthquakes were divided into six groups according to their magnitudes: M L 7.0, M L 7.1, M L 7.2, M L 7.4, M L 7.6, and M L 8.0 (for earthquakes above magnitude M s 7.0 or M L 7.2 , M L is approximately equal to M s ). The relationship between an earthquake’s magnitude and fault buffer distance was obtained for each group using least square error.
To examine the relationship between the magnitude of the mainshock and buffer distance, and to test its effectiveness in deleting aftershocks for earthquakes of other magnitudes, we fitted a function between the magnitude and the buffer distance (Figure 3) and analyzed the results in Formula (3). The fitting effect was excellent, with a low SSE (sum of squares of residuals), a low RMSE (root mean square error), and a coefficient of determination (R-square) of 0.99. In Section 3.1, we provide a detailed statistical analysis of aftershock deletion for mainshocks of different magnitudes, including M L 7.0, M L 7.1, M L 7.2, M L 7.4, M L 7.6, and M L 8.0.
l o g R = 0.2621 M 0.2682 ,

3.1. Verification of M L 7.0~ M L 7.1 Earthquake Fault Buffer Zone

Table 2 shows the M L 7.0~7.1 earthquakes that occurred in mainland China since 1980. There was a total of six earthquakes, and the fault buffer zone calculated according to Formula (3) for each earthquake is shown in the table. These buffer zones can be used for aftershock deletion in future earthquakes with similar magnitudes.
Figure 4 illustrates the distribution of aftershocks and target active faults of the M L 7.0 and M L 7.1( ~ M s 6.6– M s 6.9) earthquakes. The red dashed area represents the fault buffer zone calculated using Formula (3), while the black circular area denotes the spatial range of the deleted aftershocks selected by the K-K method corresponding to the mainshock. In Figure 4a, the fault buffer of a shallow strike-slip M L 7.0 earthquake near the western end of the Muji graben basin in the northern part of the Pamir syntax [28] is shown. The buffer distance, which is equal to 36 km, is obtained by applying Formula (3). The number of aftershocks falls within the K-K method prediction, with aftershocks mainly between M L 3.0 and M L 4.8. The number of aftershocks decreases year by year, and the focal depth ranges from 4 to 16 km.
On the other hand, Figure 4b–f show the fault buffer zone of several M L 7.1 earthquakes, with a calculated distance of 39 km using Formula (3). The dotted red line represents the buffer zone, which basically captures the aftershock distribution range. Aftershock magnitudes are generally small, concentrated in M L 3.0–4.0.
The figure demonstrates that the distribution of aftershocks determined by the fault buffer zone generally conforms to the structural geological distribution. In Figure 4b, a strike-slip fault with a slight overthrust component [29] occurred between the southwest of Qinghai Province and the edge of the Tuotuo River. The aftershocks are mainly concentrated in the area near the fault. The mainshock (shown in Figure 4c) occurred near the Tanggemu Farm on the southwestern edge of the Gonghe Graben Basin, Qinghai, and the magnitude of aftershocks ranges from M L 3.0 to M L 4.0. The aftershocks in Figure 4d mainly occur at the junction of China, Tajikistan, and Kyrgyzstan, with a broad range of magnitudes from M L 3.0 to M L 6.4 [30]. The mainshock in Figure 4e occurred in Jinggu, Yunnan, with its epicenter located at the junction of various subduction, collision, suture, overlap, and strong crustal accretion [31]. The magnitude of aftershocks is mostly M L 3.0–4.0. The aftershocks in Figure 4f are located around the mainshock and extended along the northwest and southeast sides of the mainshock [32], with the largest aftershock magnitude being M L 5.6.
In summary, the fault buffer distance of 39 km is appropriate for M L 7.0– M L 7.1 earthquakes. As shown in Figure 4, the aftershocks within the fault buffer zone are more concentrated and better reflect the correlation between the aftershock sequence and fault strike compared to the aftershocks filtered out by the K-K method.

3.2. Verification of M L 7.2 Earthquake Fault Buffer Zone

Table 3 displays the M L 7.2 earthquakes that have occurred in mainland China since 1980, with a total of six earthquakes selected for the statistics of aftershock deletion, using them as the mainshock.
According to Formula (3), the fault buffer distance of 42 km is suitable for removing aftershocks from all the M L 7.2 earthquakes. The aftershocks of the Daofu earthquake, in Sichuan Province, shown in Figure 5a, are distributed on the northwest strip and have generally small magnitudes and shallow focal depths, except for a few earthquakes. The mainshock in Figure 5b occurred at the Lijiang-Daju fault, which is a normal fault with a left-handed movement component, and the aftershocks are distributed in the north-south direction, with the maximum aftershock magnitude of 6.3. The mainshock in Figure 5c occurred in the Karakoram Pass southwest of Hotan, Xinjiang. Most aftershocks are distributed in the southwest of the mainshock with a low magnitude, and the earthquake distribution shows no obvious characteristics [33]. The lack of earthquake data may be the reason for this, as it is in the border area between two countries, some earthquakes may not have been recorded.
The mainshock in Figure 5d occurred in Zhongba, Tibet, with a generally small aftershock magnitude. The mainshock in Figure 5e occurred in Lushan County, Yaan City, Sichuan Province, which is similar in focal nature to the Wenchuan earthquake (in 2008). In Figure 5f, the mainshock occurred in Jiuzhaigou County, Aba Tibetan Autonomous Prefecture, and the northern Sichuan Province. The aftershocks of this earthquake are characterized by a small number of moderate and strong aftershocks and a large difference in magnitude between the largest aftershock and the mainshock. The sequence of aftershocks showed a NW-SE trend [34].

3.3. Verification of M L 7.4( M s 7.4) Earthquake Fault Buffer Zone

Table 4 shows four mainshocks with a magnitude of M L 7.4, all of which occurred in Western China.
The fault buffer distance of 48 km is suitable for all mainshocks. The mainshock in Figure 6a happened in Wuqia County, Xinjiang Uygur Autonomous Region. The mainshock in Figure 6b is close to the Kunlun Mountain Watershed at the junction of Yutian, Cele, and Tibet. This earthquake is believed to have been caused by a tensile rupture with a strike-slip component of the Guozhacuo fault of the Altyn Tagh fault zone under the action of the NS direction force [30], and the aftershock sequence attenuates quickly, mainly in Tibet, Ritu, and parts of Xinjiang. The mainshock in Figure 6c occurred in Yushu, Qinghai Province, and the largest aftershock magnitude is M L 6.6. The region with the most concentrated aftershocks is located west of the mainshock, and the spatial distribution of aftershocks gradually shrank westward over time [35]. The mainshock in Figure 6d occurred in Yutian, Xinjiang, and its aftershocks are distributed in a southwest direction and spread along one side of the fault.

3.4. Verification of M L 7.6( M s 7.6) Earthquake Fault Buffer Zone

There are two mainshocks with M L 7.6, and the fault buffer distance is 54 km (Table 5).
The mainshock in Figure 7a, which occurred in Mani, Tibet, has a buffer distance of 54 km. The aftershocks have small magnitudes, and the sequence decays faster with time. The mainshock in Figure 7b occurred in Lancang County, at the junction of Gengma and Cangyuan Counties in southwest Yunnan. The aftershocks are densely distributed along the Longling-Lancang fault, which includes the surrounding dense earthquake clusters. It was found that the maximum aftershock magnitude was 7.2. This earthquake can also be considered a double mainshock aftershock type.

3.5. Verification of M L 8.0( M s 8.0) Earthquake Fault Buffer Zone

Since 1980, there have been two earthquakes in mainland China with a magnitude greater than M L 8.0( M s 8.0). The first was the M L 8.1( M s 8.1) earthquake that occurred west of Kunlun Mountain Pass in 2001. The second was the Wenchuan earthquake on May 12, 2008, which caused many casualties. The Wenchuan earthquake occurred in Sichuan Province; its epicenter was at 103.4°E, 31°N, and had a focal depth of 14 km [36]. The main fault responsible for this earthquake is the Yingxiu-Beichuan fault [37].
Figure 8 shows the buffer distance of the M L 8.0 Wenchuan earthquake, which is 66 km. Aftershocks mainly occurred in areas near the mainshock, such as Wenchuan, Beichuan, Mianzhu and the number of aftershocks was large. However, the magnitude of aftershocks was generally small, and it decreased slowly over time. The aftershocks were densely distributed along the decreased direction of the mainshock.

4. Test of Fault Buffer Zone of Mainshock and Discussion

Formula (3) was obtained using the trial-and-error method based on an earthquake with a magnitude range of M L 7.0– M L 8.0( ~ M s 5.6– M s 8.0). Tests have shown that the improved algorithm, based on the fault buffer zone, is also suitable for different magnitudes, such as M L 8.1 ( M s 8.1), M L 6.8( ( M s 6.6)), and M L 5.9 ( M s 5.6) (Table 6). The results of aftershocks selected by both methods are presented in Figure 9 and Figure 10. When compared with the original K-K algorithm, the fault buffer zone method selects aftershocks that are clustered and centered on the target fault, which is consistent with the correlation between earthquake occurrence and geological structure. The earthquake shown in Figure 9a occurred in the west of the Kunlun Mountain Pass, and it is the largest earthquake in China in 70 years. Formula (3) was used to calculate the buffer distance, which was found to be 71.58 km. The earthquake shown in Figure 9b occurred in Xiaojin County, Aba Tibetan, and Qiang Autonomous Prefecture, Sichuan Province, with a buffer distance of 32.67 km. The earthquake in Figure 9c occurred in Hejing, Xinjiang, with a buffer distance of 18.98 km. Figure 9d is the aftershock deletion of Changning earthquake, Sichuan, China. The aftershock sequence has a low magnitude and a concentrated distribution, indicating that this method can accurately and quickly delete aftershocks, even for earthquakes with lower magnitudes. The three applications demonstrate that our method takes the fault strike into account. Figure 10 shows the precise location of the recent Changning earthquake in Sichuan and the range of deleted aftershocks using this algorithm. The figure shows that this algorithm can effectively distinguish between aftershocks of different mainshocks.
The improved K-K aftershock deletion algorithm based on the fault buffer zone is more effective and accurate than the original K-K method. It selects aftershock sequences that are more concentrated and distributed along the effective faults, leading to a reduction in the number of aftershocks with time while the magnitude of the aftershocks is generally small. This approach is applicable to earthquakes above M L 6.3 and provides a more accurate aftershock sequence.

5. Conclusions

Removing aftershocks, that is, separating clustered events (i.e., aftershocks) from background earthquakes in the epicenter area. Accurately obtaining the aftershock of a mainshock can reveal changes in seismic activity in time and space [38]. Based on this, the accurate Coulomb stress and b value can be calculated, and the change characteristics of the fault stress and subsequent earthquake occurrence trend can be obtained.
In this paper, we propose a novel method for removing aftershocks based on the fault buffer zone. Our method builds on the basic assumption that there is a high correlation between earthquake magnitude and a seismogenic fault. Instead of using a circular area centered on the earthquake epicenter, our approach employs the fault buffer zone to accurately identify and remove aftershocks. The conclusions of this paper are as follows:
  • We conducted a statistical analysis of earthquakes with magnitudes of M L 7.0~8.0 in mainland China since 1980 and divided the 20 earthquakes into six groups for analysis, ranging from M L 7.0, M L 7.1, M L 7.2, M L 7.4, M L 7.6, and M L 8.0 ( o r   ~ M s 6.6–8.0 and for earthquakes above magnitude M s 7 or M L 7.2 , M L is approximately equal to M s ). By using a trial-and-error method, we established an empirical formula for the fault buffer distance. Our trials on the aftershock deletion of other magnitude earthquakes above M L 5.9 show that this empirical formula is also effective, and the deleted aftershock sequences are accurately clustered using the fault buffer zone-based spatial window method. We compared our method with the K-K method and found that our method is applicable for deleting aftershocks at different magnitudes.
  • The test results show that the proposed fault buffer zone-based method improves on the traditional K-K method by taking into account relationships between earthquake magnitude and fault rupture length, increasing the accuracy of aftershock selection. The results of aftershock deletion using this method show good agreement with the actual earthquake occurrence distribution.
However, this method still has room for improvement, as the presence of a concealed fault may affect the accuracy of aftershock deletion. Furthermore, it is necessary to verify the applicability of the empirical formula in other regions of the world.

Author Contributions

Conceptualization, G.Y. and B.Z.; methodology, G.Y.; statistical analysis, Y.L.; validation, G.Y., Y.L. and B.Z.; formal analysis, G.Y.; writing—original draft preparation, G.Y. and Y.L.; writing—review and editing, G.Y. and B.Z.; visualization, Y.L.; supervision, G.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant No. 42174104, U1939204, and 41204014) and Hubei provincial Natural Science Foundation of China (2022CFB350).

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the editors and the anonymous reviewers for their comments and constructive suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Chaudhary, M.T.; Piracha, A. Natural Disasters—Origins, Impacts, Management. Encyclopedia 2021, 1, 1101–1131. [Google Scholar] [CrossRef]
  2. Kamranzad, F.; Memarian, H.; Zare, M. Earthquake Risk Assessment for Tehran, Iran. ISPRS Int. J. Geoinf. 2020, 9, 430. [Google Scholar] [CrossRef]
  3. Devries, P.M.R.; Viégas, F.; Wattenberg, M.; Meade, B.J. Deep learning of aftershock patterns following large earthquakes. Nature 2018, 560, 632–634. [Google Scholar] [CrossRef] [PubMed]
  4. Gardner, J.K.; Knopoff, L. Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian? Bull. Seismol. Soc. Am. 1974, 64, 1363–1367. [Google Scholar] [CrossRef]
  5. Taroni, M.; Akinci, A. Good practices in PSHA: Declustering, b-value estimation, foreshocks and aftershocks inclusion; A case study in Italy. Geophys. J. Int. 2021, 224, 1174–1187. [Google Scholar] [CrossRef]
  6. Shilen, S.; Toksoz, M.N. A statistical method of identifying dependent events and earthquake aftershocks. Seismol. Res. Lett. 1974, 45, 3–16. [Google Scholar] [CrossRef]
  7. Sheng, J.Q.; Ma, Q.H.; Ren, X.M.; Zhao, W.M. Aftershock deletion method based on comparison of regional seismic activity frequency. Earthquake 2007, 27, 85–94. [Google Scholar]
  8. Chen, L.; Liu, J.; Chen, Y.; Chen, L.S. Deletion of aftershocks in seismic activity analysis. Chin. J. Geophys. 1998, 41 (Suppl. 1), 244–252. (In Chinese) [Google Scholar]
  9. Knopoff, L.; Gardner, J.K. Higher seismic activity during local night on the raw worldwide earthquake catalogue. Geophys. J. Int. 1972, 3, 3. [Google Scholar] [CrossRef] [Green Version]
  10. Shlien, S. A clustering model for earthquake occurrences. Bull. Seismol. Soc. Am. 1970, 60, 1765–1787. [Google Scholar]
  11. Console, R.; Gasparini, C.; Simoni, B.D.; Marcelli, L.; Vecchi, M. Preambolo al catalogo sismico nazionale (csn). i criteri di informazione del catalogo sismico nazionale (csn). Ann. Geophys. 1979, 32, 37–77. [Google Scholar] [CrossRef]
  12. Keilis-Borok, V.I.; Knopoff, L.; Rotvain, I.M. Bursts of aftershocks, long-term precursors of strong earthquakes. Nature 1980, 283, 259–263. [Google Scholar] [CrossRef]
  13. Davis, S.D.; Cliff, F. Single-link cluster analysis, synthetic earthquake catalogues, and aftershock identification. Geophys. J. Int. 1991, 2, 289–306. [Google Scholar] [CrossRef] [Green Version]
  14. Zhuang, J.; Ogata, Y.; Vere-Jones, D. Stochastic declustering of space-time earthquake occurrences. J. Am. Stat. Assoc. 2002, 97, 369–380. [Google Scholar] [CrossRef]
  15. Zaliapin, I.; Ben-Zion, Y. Earthquake clusters in southern California I: Identification and stability. J. Geophys. Res. Solid Earth. 2013, 118, 2847–2864. [Google Scholar] [CrossRef]
  16. Zaliapin, I.; Ben-Zion, Y. Earthquake declustering using the nearest-neighbor approach in space-time-magnitude domain. J. Geophys. Res. Solid Earth. 2020, 125, e2018JB017120. [Google Scholar] [CrossRef]
  17. Zaliapin, I.; Ben-Zion, Y. Perspectives on clustering and declustering of earthquakes. Seismol. Res. Lett. 2022, 93, 386–401. [Google Scholar] [CrossRef]
  18. Rosson, Z.; Walter, J.I.; Goebel, T.; Chen, X. Narrow spatial aftershock zones for induced earthquake sequences in Oklahoma. Geophys. Res. Lett. 2019, 46, 10358–10366. [Google Scholar] [CrossRef] [Green Version]
  19. Baiesi, M.; Paczuski, M. Scale-free networks of earthquakes and aftershocks. Phys. Rev. E—Stat. Phys. Plasmas Fluids Relat. Interdiscip. Top. 2004, 69, 8. [Google Scholar] [CrossRef] [Green Version]
  20. Dascher-Cousineau, K.; Brodsky, E.E.; Lay, T.; Goebel, T. What controls variations in aftershock productivity? J. Geophys. Res. Solid Earth. 2020, 125, e2019JB018111. [Google Scholar] [CrossRef]
  21. Zhan, Z.; Jin, B.; Wei, S.; Graves, R.W. Coulomb stress change sensitivity due to variability in mainshock source models and receiving fault parameters: A case study of the 2010–2011 Christchurch, New Zealand, earthquakes. Seismol. Res. Lett. 2011, 82, 800–814. [Google Scholar] [CrossRef] [Green Version]
  22. Sboras, S.; Lazos, I.; Mouzakiotis, E.; Karastathis, V.; Pavlides, S.; Chatzipetros, A. Fault modelling, seismic sequence evolution and stress transfer scenarios for the July 20, 2017 (MW 6.6) Kos–Gökova Gulf earthquake, SE Aegean. Acta Geophys. 2020, 68, 1245–1261. [Google Scholar] [CrossRef]
  23. Laurenti, L.; Tinti, E.; Galasso, F.; Franco, L.; Marone, C. Deep learning for laboratory earthquake prediction and autoregressive forecasting of fault zone stress. Earth Planet Sci. Lett. 2022, 598, 117825. [Google Scholar] [CrossRef]
  24. Lazos, I.; Pikridas, C.; Chatzipetros, A.; Pavlides, S. Determination of local active tectonics regime in central and northern Greece, using primary geodetic data. Appl. Geomat. 2021, 13 (Suppl. 1), 3–17. [Google Scholar] [CrossRef]
  25. Guo, Z.J.; Qin, B.Y.; Xu, W.Y.; Tang, Q. Preliminary study on a model for the development of the focus of an earthquake. Chin. J. Geophys. 1973, 16, 43–48. [Google Scholar]
  26. Wells, D.; Coppersmith, K. New Empirical Relationships among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement. Bull. Seismol. Soc. Am. 1994, 84, 974–1002. [Google Scholar]
  27. Deng, Q.D.; Zhang, P.Z.; Ran, Y.K.; Yang, X.P.; Min, W.; Chen, L.C. Active tectonics and earthquake activities in China. Earth Sci. Front. 2003, 10 (Suppl. 1), 66–73. [Google Scholar]
  28. Chen, J.; Li, T.; Sun, J.B.; Fang, J.H.; Yao, Y.; Li, Y.H.; Wang, H.R.; Fu, B. Coseismic surface ruptures and seismogenic Muji fault of the 25 November 2016 Arketao Mw6. 6 earthquake in northern Pamir. Seismol. Geol. 2016, 38, 1160–1174. [Google Scholar]
  29. Zhou, H.L. Source mechanism of earthquake (Nov.5, 1988) located between Xiji and Totohe Bark in Qinghai Province. Earthq. Res. China. 1995, 11, 84–91. (In Chinese) [Google Scholar]
  30. Li, Z.H.; Ma, H.S.; Qu, Y.J. Study on seismogenic structure and seismic activity characteristics before the Yutian M7.3 earthquake on March 21, 2008, Xinjiang. Earthq. Res. China. 2009, 25, 199–205. (In Chinese) [Google Scholar]
  31. Xu, P.S.; Liu, Z.F.; Zhang, Z.Q.; Li, J.; Liu, L.F.; Su, Y.J. Double difference relocation and focal mechanisms of the Jinggu Ms6.6 earthquakes sequences in Yunnan Province in 2014. Earth Sci.—J. China Univ. Geosci. 2015, 40, 1741–1754. (In Chinese) [Google Scholar]
  32. Wei, W.; Xie, C.; Zhou, B.; Zhi, G.; Wang, Y. Location of the mainshock and aftershock sequences of the M6.9 Mainling earthquake, Tibet. Chin. Sci. Bull. 2018, 63, 1493–1501. (In Chinese) [Google Scholar] [CrossRef] [Green Version]
  33. Luo, F.Z.; Qu, Y.J.; Wang, J. Study on the M7.1 earthquake at the mouth of Kalakunlun MT. in Hetian, Xinjiang on Nov.19, 1996. Inland Earthq. 2003, 17, 33–38. (In Chinese) [Google Scholar]
  34. Wu, W.W.; Wei, Y.L.; Long, F.; Liang, M.J.; Chen, X.F.; Sun, W.; Zhao, J. Study on source parameters of the 8 August 2017 M7.0 Jiuzhaigou earthquake and its aftershocks, Northern Sichuan. Seismol. Geol. 2020, 42, 492–512. [Google Scholar]
  35. Shao, Z.G.; Ma, H.S.; Zhang, L.P.; Huang, J.P.; Wen, X.Z. The characteristics of co-seismic slip and aftershocks distribution of the Ms7.1 earthquake at Qinghai Yushu in 2010 and its relationship with tectonics. Chin. J. Geophys. 2013, 56, 3800–3810. (In Chinese) [Google Scholar]
  36. Lei, J.; Zhao, D. Structural heterogeneity of the Longmenshan fault zone and the mechanism of the 2008 Wenchuan earthquake (Ms8.0). Geochem. Geophys. Geosyst. 2009, 10. [Google Scholar] [CrossRef]
  37. Klimenko, M.V.; Klimenko, V.V.; Zakharenkova, I.E.; Pulinets, S.A.; Zhao, B.; Tsidilina, M.N. Formation mechanism of great positive TEC disturbances prior to Wenchuan earthquake on May 12, 2008. Adv. Space Res. 2011, 48, 488–499. [Google Scholar] [CrossRef]
  38. Chen, M. Aftershocks and Background Seismicity in Tangshan and the Rest of North China. J. Geophys. Res. Solid Earth 2021, 126, e2020JB021395. [Google Scholar] [CrossRef]
Figure 1. Flow chart of aftershock deletion algorithm based on the fault buffer zone.
Figure 1. Flow chart of aftershock deletion algorithm based on the fault buffer zone.
Remotesensing 15 01662 g001
Figure 2. Method for establishing and adjusting the fault buffer zone (The red five-pointed star represents the main earthquake, while the black line represents its corresponding seismogenic fault. The black line on the left of the figure above represents the seismogenic fault. The circles represent the buffer range of each point on the fault, and the red envelope represents the buffer area of the fault. The multiple envelope areas on the right of the figure above represent different buffer areas corresponding to different earthquakes). (a) The buffer distance for a earthquake; (b) adjustment of the buffer distance.
Figure 2. Method for establishing and adjusting the fault buffer zone (The red five-pointed star represents the main earthquake, while the black line represents its corresponding seismogenic fault. The black line on the left of the figure above represents the seismogenic fault. The circles represent the buffer range of each point on the fault, and the red envelope represents the buffer area of the fault. The multiple envelope areas on the right of the figure above represent different buffer areas corresponding to different earthquakes). (a) The buffer distance for a earthquake; (b) adjustment of the buffer distance.
Remotesensing 15 01662 g002
Figure 3. Illustrates the logarithmic relationship between the magnitude and the fault buffer distance (The red triangle indicates the statistical optimal fracture buffer distance, while the blue dotted line represents the fitted trend of fault buffer distance. The fitting correlation coefficient is shown in the figure, with the SSE and RMSE being very low, and the R-square value reaching 0.99; SSE: sum of squares of residuals; RMSE: root mean square error; R-square: coefficient of determination).
Figure 3. Illustrates the logarithmic relationship between the magnitude and the fault buffer distance (The red triangle indicates the statistical optimal fracture buffer distance, while the blue dotted line represents the fitted trend of fault buffer distance. The fitting correlation coefficient is shown in the figure, with the SSE and RMSE being very low, and the R-square value reaching 0.99; SSE: sum of squares of residuals; RMSE: root mean square error; R-square: coefficient of determination).
Remotesensing 15 01662 g003
Figure 4. The fault buffer zone of M L 7.0 and M L 7.1 earthquakes, as well as the distribution of aftershocks (the red dotted line indicates the fault buffer zone, while the black circle indicates the range of K-K aftershocks; the red five-pointed star represents the mainshock; and the time of the mainshock is displayed in the upper left corner of the figure; the legend in the lower right corner of the figure explains the symbols used in the figure). (a) The Aktau M L 7.0 ( M s 6.7) earthquake; (b) the Xiji and toto river M L 7.1 ( M s 6.8) earthquake; (c) the Tanggemu M L 7.1 ( M s 6.9) earthquake; (d) the Wuqia M L 7.1 ( M s 6.8) earthquake; (e) the Jinggu M L 7.1 ( M s 6.6) earthquake; (f) the Linzhi M L 7.1 ( M s 6.9) earthquake.
Figure 4. The fault buffer zone of M L 7.0 and M L 7.1 earthquakes, as well as the distribution of aftershocks (the red dotted line indicates the fault buffer zone, while the black circle indicates the range of K-K aftershocks; the red five-pointed star represents the mainshock; and the time of the mainshock is displayed in the upper left corner of the figure; the legend in the lower right corner of the figure explains the symbols used in the figure). (a) The Aktau M L 7.0 ( M s 6.7) earthquake; (b) the Xiji and toto river M L 7.1 ( M s 6.8) earthquake; (c) the Tanggemu M L 7.1 ( M s 6.9) earthquake; (d) the Wuqia M L 7.1 ( M s 6.8) earthquake; (e) the Jinggu M L 7.1 ( M s 6.6) earthquake; (f) the Linzhi M L 7.1 ( M s 6.9) earthquake.
Remotesensing 15 01662 g004
Figure 5. Establishment of M L 7.2 earthquake fault buffer zone. (a) The Daofu M L 7.2 ( M s 6.8) earthquake; (b) the Lijiang M L 7.2 ( M s 7.0) earthquake; (c) the Hetian M L 7.2 ( M s 7.1) earthquake; (d) the Zhongba M L 7.2 ( M s 6.8) earthquake; (e) the Lushan M L 7.2 ( M s 7.0) earthquake; (f) the Jiuzhaigou M L 7.2 ( M s 7.0).
Figure 5. Establishment of M L 7.2 earthquake fault buffer zone. (a) The Daofu M L 7.2 ( M s 6.8) earthquake; (b) the Lijiang M L 7.2 ( M s 7.0) earthquake; (c) the Hetian M L 7.2 ( M s 7.1) earthquake; (d) the Zhongba M L 7.2 ( M s 6.8) earthquake; (e) the Lushan M L 7.2 ( M s 7.0) earthquake; (f) the Jiuzhaigou M L 7.2 ( M s 7.0).
Remotesensing 15 01662 g005
Figure 6. Establishment of M L 7.4 earthquake fault buffer zone. (a) The Wuqia M L 7.4 ( M s 7.4) earthquake; (b) the Yutian M L 7.4 ( M s 7.4) earthquake; (c) the Yushu M L 7.4 ( M s 7.1) earthquake; (d) the Yutian M L 7.3 ( M s 7.3) earthquake.
Figure 6. Establishment of M L 7.4 earthquake fault buffer zone. (a) The Wuqia M L 7.4 ( M s 7.4) earthquake; (b) the Yutian M L 7.4 ( M s 7.4) earthquake; (c) the Yushu M L 7.4 ( M s 7.1) earthquake; (d) the Yutian M L 7.3 ( M s 7.3) earthquake.
Remotesensing 15 01662 g006
Figure 7. Establishment of M L 7.6) earthquake fault buffer zone. (a) the Mani M L 7.6 ( M s 7.9) earthquake. (b) the Lancang-gengma M L 7.6 ( M s 7.6) earthquake.
Figure 7. Establishment of M L 7.6) earthquake fault buffer zone. (a) the Mani M L 7.6 ( M s 7.9) earthquake. (b) the Lancang-gengma M L 7.6 ( M s 7.6) earthquake.
Remotesensing 15 01662 g007
Figure 8. Establishment of M L 8.0 ( M s 8.0) earthquake fault buffer zone.
Figure 8. Establishment of M L 8.0 ( M s 8.0) earthquake fault buffer zone.
Remotesensing 15 01662 g008
Figure 9. Establishment of three earthquakes’ fault buffer zones in Table 6. (a) the Kunlun Mountain Pass M L 8.1 ( M s 8.1) earthquake. (b) the Xiaojin M L 6.8 ( M s 6.6) earthquake; (c) the Hejing M L 5.9 ( M s 5.6) earthquake; (d) the Changning M L 6.3 ( M s 6.0) earthquake.
Figure 9. Establishment of three earthquakes’ fault buffer zones in Table 6. (a) the Kunlun Mountain Pass M L 8.1 ( M s 8.1) earthquake. (b) the Xiaojin M L 6.8 ( M s 6.6) earthquake; (c) the Hejing M L 5.9 ( M s 5.6) earthquake; (d) the Changning M L 6.3 ( M s 6.0) earthquake.
Remotesensing 15 01662 g009
Figure 10. The aftershock deletion of Changning earthquake, Sichuan, China using the fault buffer zone algorithm.
Figure 10. The aftershock deletion of Changning earthquake, Sichuan, China using the fault buffer zone algorithm.
Remotesensing 15 01662 g010
Table 1. Aftershock space and time window of the K-K method.
Table 1. Aftershock space and time window of the K-K method.
M R 0 ( M ) / k m T 0 ( M ) / d M R 0 ( M ) / k m T 0 ( M ) / d
2.0–2.53065.0–5.550183
2.5–3.530125.5–6.550365
3.5–4.040236.5–7.0100548
4.0–4.540467.0–7.5100730
4.5–5.040927.5–8.0150913
Table 2. Fault buffer distance of the M L 7.0 and M L 7.1 earthquakes.
Table 2. Fault buffer distance of the M L 7.0 and M L 7.1 earthquakes.
ID M L DateLocationTarget FaultBuffer Distance (km)
LongitudeLatitudeDepth
a7.025 November 201674.1°39.2°10Muji fault36
b7.15 November 198890.7°32.9°/Chibu zhangcuo
-zhasa fault
39
c7.126 April 1990100.13°36.12°29Northern margin fault of Gonghe basin39
d7.15 October 200874.03°39.45°10Kazikearte fault39
e7.17 October 2014100.55°23.4°10Sanlinyang-siyongjie fault39
f7.118 November 201795°29.75°10Ani bridge fault39
Table 3. Fault buffer distance of M L 7.2 earthquake.
Table 3. Fault buffer distance of M L 7.2 earthquake.
ID M L DateLocationTarget FaultBuffer Distance (km)
LonLatDepth
a7.224 January 1981101.17°31°12Xianshuihe fault42
b7.23 February 1996100.22°27.3°10Lijiang-Daju fault42
c7.219 November 199678.35°35.43°/Eslak karaur42
d7.225 August 200883.65°30.93°10Yare fault42
e7.220 April 2013102.98°30.3°17Guanxian-Jiangyou fault42
f7.28 August 201795°29.75°10Tazang fault42
Table 4. Fault buffer distance of M L 7.4 earthquake.
Table 4. Fault buffer distance of M L 7.4 earthquake.
ID M L DateLocationMain FaultBuffer Distance (km)
LonLatDepth
a7.423 August 198575.6°39.58°/Kazikealte fault48
b7.421 March 200881.43°35.8°/Altun Mountain’s south margin fault48
c7.414 April 201096.58°33.22°14Ganzi-Yushu- Fenghuoshan48
d7.412 February 201482.52°36.13°10Yare fault48
Table 5. Fault buffer distance of M L 7.6 earthquake.
Table 5. Fault buffer distance of M L 7.6 earthquake.
ID M L DateLocationTarget FaultBuffer Distance (km)
LonLatDepth
a7.618 November 199787.3°35.2°/Altun Mountain’s south margin fault54
b7.66 November 198899.72°22.83°13longling-lancang fault54
Table 6. Fault buffer distance of other earthquakes.
Table 6. Fault buffer distance of other earthquakes.
ID M L DateLocationTarget FaultInterception Method
LonLatDepth
a8.114 November 200190.53°35.93°10Muziluk-Whale LakeFar-end
b6.822 September 1989102.38°31.55°10Malkang faultFar-end
c5.930 August 200883.85°42.72°6Southern margin fault of Greater Yuludus BasinBoth sides
d6.317 June
2019
104.97°104.97°16104.97°Far-end
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, G.; Zhao, B.; Liu, Y. An Aftershock Deletion Method Based on Fault Buffer Zone. Remote Sens. 2023, 15, 1662. https://doi.org/10.3390/rs15061662

AMA Style

Yang G, Zhao B, Liu Y. An Aftershock Deletion Method Based on Fault Buffer Zone. Remote Sensing. 2023; 15(6):1662. https://doi.org/10.3390/rs15061662

Chicago/Turabian Style

Yang, Guangliang, Bingjie Zhao, and Yijun Liu. 2023. "An Aftershock Deletion Method Based on Fault Buffer Zone" Remote Sensing 15, no. 6: 1662. https://doi.org/10.3390/rs15061662

APA Style

Yang, G., Zhao, B., & Liu, Y. (2023). An Aftershock Deletion Method Based on Fault Buffer Zone. Remote Sensing, 15(6), 1662. https://doi.org/10.3390/rs15061662

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop