Next Article in Journal
Study on the Smoke Control Performance of Air Curtains in Emergency Stairwells of Multi-Unit Residential Buildings
Next Article in Special Issue
Dynamic Stability Finite Difference Time Domain Analysis of Landfill Based on Hypergravity Test
Previous Article in Journal
Vehicle Multi-Object Detection and Tracking Algorithm Based on Improved You Only Look Once 5s Version and DeepSORT
Previous Article in Special Issue
Urban Seismic Networks: A Worldwide Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Strong Earthquake Risk in Maqin–Maqu Segment of the Eastern Kunlun Fault, Northeast Tibet Plateau

Key Laboratory of Active Tectonics and Volcano, Institute of Geology, China Earthquake Administration, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(7), 2691; https://doi.org/10.3390/app14072691
Submission received: 11 February 2024 / Revised: 11 March 2024 / Accepted: 19 March 2024 / Published: 22 March 2024
(This article belongs to the Special Issue Advanced Research in Seismic Monitoring and Activity Analysis)

Abstract

:
The East Kunlun Fault Zone, as a highly seismically active fault, has witnessed five earthquakes with magnitudes exceeding M7.0 to the west of Animaqing Mountain since 1900. Conversely, the historical records for the Maqin–Maqu segment in the east of Animaqing Mountain show no M7.0 or above earthquakes, designating it as a distinctive seismic gap within this fault zone. We analyzed the tectonic background and structural features of the Maqin–Maqu segment within the East Kunlun Fault Zone to evaluate its potential seismic capacity. Utilizing a new established probability recurrence model, we calculated the seismic hazard for both segments over the next 100 years. The results indicate that the probability of M7.0 or above earthquake occurring in the Maqu segment in the next 100 years is 11.47%, classified as a moderate probability event. The joint probability of at least one M7.0 or above strong earthquake occurring in the entire Maqin–Maqu segment in the next 100 years is 16.14%, also classified as a moderate probability event, while the probability for the Maqin segment alone is 5.36%, classified as a low probability event. Considering the uncertainty of the probability model, a qualitative hazard classification for each segment was further conducted. The comprehensive evaluation suggests a low risk of a major earthquake occurring in the Maqin segment in the next 100 years, while the Maqu segment is assessed to have a higher risk.

1. Introduction

The Eastern Kunlun Fault constitutes a significant fault zone characterized by primarily sinistral strike–slip movement. It originated along the Eastern Kunlun Variscan paleo-tectonic suture within the Tibetan Plateau due to the collision between the Indian and Asian plates. This fault spans approximately 2000 km, following a general west-northwest trend, extending from the west of Bukedaban Peak, situated north of Hoh Xil Lake in Qinghai province, to the Maqu region in Sichuan province [1] (Figure 1). Situated at the northern margin of an upwelling area in the middle of the Tibetan Plateau, the primary portion of this fault zone lies approximately along the 64-km depth contour of the Moho surface, extending in an east–west direction. An isobathic line from the southeast of Maduo-Maqu to the Mohorovicic discontinuity shifts southward, coinciding with a gradual reduction in crustal thickness along this trajectory [2]. Recent seismic refraction data reveal that the crustal thickness of the BaiYan Har Mountain Group, positioned in the southern sector of the fault zone, measures 64 km. This thickness notably decreases to 60 km as one progresses northward, reaching the southern rim of the Qaidam Basin and the Touraine area. This creates a region characterized by variable crustal thickness in the west–northwest direction. The Moho topography and crustal flow mechanism within this fault zone have experienced discernible shifts, delineated by the Aemye Ma-chhen Mountains in an east–west orientation [3,4]. As a result, this paper employs the Aemye Ma-chhen Mountains as a reference to partition the Eastern Kunlun Fault into eastern and western segments, building upon the geometric structural attributes of the fault zone in conjunction with geophysical data.
Since 1900, five earthquakes with magnitudes exceeding 7.0 have occurred in the western part of the Aemye Ma-chhen Mountains in the Eastern Kunlun fault, encompassing the Ms 7.9 Mani earthquake in 1997, Ms 8.1 Kunlun Mountains earthquake in 2001, Ms 7.0 Xiugou earthquake in 1902, Ms 7.1 Akake Lake earthquake in 1963, and Ms 7.5 Tuosuohu earthquake in 1937 (Figure 1). As depicted in Figure 1, there is a notable absence of documented large earthquakes (Ms ≥ 7) in the historical records of the eastern segment of the Eastern Kunlun Fault, marking it as a distinctive seismic gap within this fault zone. Moreover, the Wenchuan earthquake of 2008 heightened the accumulation of Coulomb stress within the eastern segment of the Eastern Kunlun Fault zone [5], consequently elevating the latent risk of earthquakes in this section. Therefore, the Maqin–Maqu segment, constituting the eastern part of the Eastern Kunlun Fault zone, emerges as the seismic gap with the highest population density along the entire fault zone. As such, there is an exigent need to scrutinize the earthquake risk and forecast potential large-scale seismic events.
To predict earthquakes in the medium and long term, a quantitative approach involves expressing the likelihood of earthquake occurrence as a probability [6,7]. Assuming that the time for earthquake occurrence approximately follows a specific recurrence model [8,9], we can employ this model to estimate the probability of recurrent major earthquakes within a defined time frame [10,11,12]. Many research findings indicate that the model may exhibit various distribution types, including normal, logarithmic normal, and Weibull distributions. In the northeastern Tibetan Plateau, occurrences of severe earthquakes are relatively frequent, and research on ancient earthquakes along the fault zone is notably extensive. In this study, through the collection and analysis of data regarding ancient earthquakes, we ascertain sample parameters and derive the probability density function governing earthquake recurrence in this region. This model is subsequently applied in assessing the seismic probability for the Maqin–Maqu segment of the Eastern Kunlun fault zone.

2. Regional Tectonic Setting

The Maqin–Maqu segment, situated at the confluence of Gansu, Sichuan, and Qinghai provinces, constitutes the eastern boundary fault for the Qaidam and Bayan Har blocks. It spans from the northern base of the Aemye Ma-chhen Mountains in the west, passing the north sides of Dongqinggou, Dawutan, Kendingna, Xigongzhou, and Tangdi. Progressing through Keshengtuoluo where it traverses the Yellow River, it extends into the Ruoergai Basin via Maqu. Eventually, it emerges in the swampland north of Luocha, intersecting with the Tacang fault. The total length approximates 310 km. Its overall trend aligns with 295° northwest, with a prevailing southwest inclination [2].
By delineating tectonic indicators along the segmented boundary of the strike–slip fault [6], it can be divided, from west to east, into the Maqin segment and Maqu segment (Figure 2A). The calculated second strain rate invariant indicates that the strain rate is mainly distributed along the Maqin–Maqu segment of the East Kunlun Fault Zone. In the vicinity of the intersection of Animaqing Mountain and the West Gong Co-seismic Fault, there is a high strain rate zone with maximum values of approximately 65 nstrain·a−1 and 70 nstrain·a−1, respectively. East of Maqu, the strain rate is diffusely distributed between the Gongga Shan–Dieshan Fault and the East Kunlun Fault Zone, ranging from 30 to 40 nstrain·a−1 (Figure 2B) [13]. Boundary markers between these segments and the western segment include the uplift of Animaqing Mountain (exhibiting double bending as a marker for segmentation of seismic ruptures) [14,15] (Figure 3A–C), the region of fault intersection and decomposition (triple junction where the Maqin fault intersects with the Awancang fault near the Xike River as a marker for segmentation of seismic ruptures) [16] (Figure 2A–C), and the pull apart basin of the left rank (namely, the Maqu fault aligns with the left rank of the Tacang fault through Ruoergai Basin (Figure 3D,E). Considering the Holocene slip rates for these three segments (Figure 2C), it is evident that each has maintained its distinctive level of activity throughout this period, with the slip rate gradually decreasing from west to east. This trend aligns consistently with the broader Eastern Kunlun fault zone, where activity diminishes in the same west-to-east fashion.
Figure 2. (A) Fault trace distribution of the Maqin–Maqu segments. The fault has been divided into different segments based on geometrical complexity, such as restraining bends, triple junctions, and stepover (marked by pink arrows and dashed circles).(B) The calculated second strain rate invariant indicates that the strain rate is mainly distributed along Maqin–Maqu of the East Kunlun Fault Zone. Data were compiled from Wen Y M et al. (2023) [13]. (C) Compilation of slip rates showing the partitioning of slip rates at the Maqin–Maqu–Awancang Triple Junction [17,18,19].
Figure 2. (A) Fault trace distribution of the Maqin–Maqu segments. The fault has been divided into different segments based on geometrical complexity, such as restraining bends, triple junctions, and stepover (marked by pink arrows and dashed circles).(B) The calculated second strain rate invariant indicates that the strain rate is mainly distributed along Maqin–Maqu of the East Kunlun Fault Zone. Data were compiled from Wen Y M et al. (2023) [13]. (C) Compilation of slip rates showing the partitioning of slip rates at the Maqin–Maqu–Awancang Triple Junction [17,18,19].
Applsci 14 02691 g002
Figure 3. (A) Hillshaded DEM and (B) Oblique Google Earth image showing the geometry and shape of the Anemaqen Transpressive Range and restraining double bend. (C) Topographic profile extracted from SRTM 30 m DEM showing the relief of the range. The profile location is shown in Figure 3A. (D) A simplified geological map illustrating the fault traces of the Maqu Segment and Tazang Fault. The buried Maqu fault segment was determined based on unpublished geophysical profiles. (E) A Google Earth image showing the Tazang Fault obvious sinistral strike–slip displaced mountain ridges, gullies, and Holocene alluvial fans [19].
Figure 3. (A) Hillshaded DEM and (B) Oblique Google Earth image showing the geometry and shape of the Anemaqen Transpressive Range and restraining double bend. (C) Topographic profile extracted from SRTM 30 m DEM showing the relief of the range. The profile location is shown in Figure 3A. (D) A simplified geological map illustrating the fault traces of the Maqu Segment and Tazang Fault. The buried Maqu fault segment was determined based on unpublished geophysical profiles. (E) A Google Earth image showing the Tazang Fault obvious sinistral strike–slip displaced mountain ridges, gullies, and Holocene alluvial fans [19].
Applsci 14 02691 g003

3. Database and Methodology

3.1. Database

This study entails the collection of data pertaining to ancient earthquakes across different segments of 22 active fault zones in the northeastern Tibetan Plateau (refer to Appendix A Research data Table A1). Data selection adheres to three primary criteria: (1) Collection, sorting, and comparative analysis of data to scrutinize ancient earthquake events exposed by various trench profiles within the same fault segment, while excluding repetitive events. (2) Exclusion of ancient earthquake events that transpired an excessive length of time ago. (3) Computation of the recurrence interval (T) for each fault segment to establish the average recurrence interval of ancient earthquakes (Tave). Subsequently, the ancient earthquake data (T/Tave) for each segment are normalized to mitigate data dispersion. Finally, a total of 119 data points (T/Tave) are available for statistical analysis based on the aforementioned principles.

3.2. Methodology and Model

Employing SPSS 26.0 statistical software, we construct a frequency distribution histogram illustrating the recurrence probability (T/Tave) of ancient earthquakes in the northeastern Tibetan Plateau (Figure 4). Figure 4 reveals that while variations in recurrence probabilities among different ancient earthquakes are substantial, they notably diminish after normalization processing, resulting in a distribution skewed towards values near 1.0. The concentration of T/Tave is confined within a range proximate to 1.0. Furthermore, we posit that this frequency distribution may conform to potential fits with a normal distribution, logarithmic normal distribution, and Weibull distribution. Subsequent sections utilize the Q-Q plots to evaluate potential distribution forms Figure 5, Figure 6 and Figure 7.
In statistics, a Q-Q plot (quantile–quantile plot) serves as a probability plot, offering a graphical method for comparing two probability distributions by plotting their quantiles against each other. Each point (x, y) on the plot corresponds to a quantile of the second distribution (y-coordinate) plotted against the same quantile of the first distribution (x-coordinate). This results in a parametric curve where the parameter is the index of the quantile interval [20].
The primary purpose of a Q-Q plot is to compare the shapes of distributions, providing a visual representation of how properties such as location, scale, and skewness are either similar or different in the two distributions. Q-Q plots can be applied to compare collections of data or theoretical distributions. Utilizing Q-Q plots for comparing two data samples can be seen as a non-parametric approach to assessing their underlying distributions. While Q-Q plots are generally more diagnostic than comparing sample histograms, they are less widely known. Commonly, Q-Q plots are employed to compare a dataset to a theoretical model, offering a graphical assessment of goodness of fit rather than reducing to a numerical summary statistic. Additionally, Q-Q plots are used to compare two theoretical distributions to each other. As Q-Q plots compare distributions, there is no requirement for the values to be observed as pairs, as in a scatter plot, or for the numbers of values in the two groups being compared to be equal [21].
Based on the aforementioned analysis, in conjunction with Figure 5, Figure 6 and Figure 7, it is evident that the normal distribution yields the most favorable fit, followed by the Weibull distribution, with the logarithmic normal distribution performing the least satisfactorily. Consequently, this study opts for the normal distribution as the probability density function. T/Tave is expected to adhere to the normal distribution: T/Tave~N(μ, σ2), μ represents the expected value of the normal distribution; it can be interpreted as a position parameter. σ represents the standard deviation; it can be interpreted as the amplitude of the distribution.
Based on the statistical analysis results based on Figure 4, μ = 0.9538, σ = 0.3879, T/Tave is expected to adhere to the normal distribution N(0.9538, 0.38792), and the fitting result successfully passes the χ2 test within the 95% confidence interval. This implies that recurrent events in the same location possess an approximate 95% confidence interval falling within [μ − 2σ, μ + 2σ]. Subsequently, the probability density function for T/Tave conforming to the normal distribution can be derived as follows:
f ( T / T a v e ) = 1 σ 2 π e ( T / T a v e μ ) 2 2 σ 2
The probability distribution function for the northeastern margin of the Tibetan Plateau can be derived as follows:
F ( t ) = Φ t / R 0.9538 0.3879
In Equation (2), Φ represents the standard normal distribution function, R denotes the average recurrence interval of earthquakes at a given location, and t signifies the elapsed time (the duration between the last earthquake and the present).
Employing the aforementioned probability distribution function, the conditional probability (P) of earthquake recurrence can be expressed by the following equation:
P = F ( t + Δ t ) F ( t ) 1 F ( t )
In Equation (3), t represents the elapsed time (the duration between the last earthquake and the present), R denotes the average recurrence interval of earthquakes, while Δt signifies the time increment used to project the probability of future earthquake occurrences. To compute the conditional probability of earthquake recurrence, only the three fundamental parameters—t, R, and Δt—are required.
Regarding the conditional probability model for earthquake recurrence, Zhang (1996) [12] conducted a comprehensive analysis by collating recurrence interval data from 46 ancient earthquakes within the continental interior. Employing the logarithmic normal distribution, Zhang fitted the cumulative frequency and established a universal distribution relationship applicable to all active faults within dynamic continents [12]:
f ( T / T a v e ) = T a v e T σ 2 π exp 1 2 σ 2 ln T T a v e μ 2
Applying the maximum likelihood method, we obtain estimates for the mean (μ = 0.1206) and standard deviation (σ = 0.5054). The resulting probability distribution function is as follows:
F ( t ) = Φ ln ( t / R ) + 0.1206 0.5054
In the subsequent sections, we employ the previously mentioned two distribution models to examine the earthquake risk within the Maqin–Maqu segment of the Eastern Kunlun fault zone. Subsequently, we conduct a comparative analysis and discussion of the results, thereby mutually validating their feasibility.

4. Results and Discussion

4.1. Estimation of Largest Earthquake Magnitudes

In 1994, American scholars D. L. Wells et al. [22] established an empirical relationship between the surface rupture length and the moment magnitude of strike–slip earthquakes based on a large number of samples, which is known briefly as the WC empirical relationship: Mw = 5.16 + 1.12LgL, where Mw represents the moment magnitude, L represents fault length, and the residual standard deviation is 0.28 [23]. When applying the WC empirical relationship, we encounter the problem of conversion between moment magnitude Mw and surface wave magnitude Ms. D. L. Wells suggested in the 1994 work that when Ms is in the range of 5.7–8.0, there is no systematic difference between them, namely, MwMs. However, there is a difference between the magnitude determination method and station data adopted in China and the United States; namely, Ms and Mw are not the same in the continent of China. In 2009, Ran [24] derived a transformation relationship for these two earthquake magnitudes according to seismic data in mainland China during 1973–2008: Ms = 1.412 + 0.845Mw, where Ms represents surface wave magnitude, Mw represents moment magnitude, and the residual standard deviation is 0.11 [23]. In this paper, we use the aforementioned equation of statistical relationship and combine it with the length of the fault to calculate the maximum earthquake magnitude that could occur in each segment (Table 1). Based on the above analysis, regardless of Mw or Ms, the earthquake may be induced by the fault belonging to the strong earthquake (Ms ≥ 7.0).

4.2. Estimation of Elapsed Time since the Last Earthquake (t) and Average Recurrence Interval (R)

For the evaluation of seismic risk, we must comprehensively understand the average recurrence interval of ancient earthquakes in this area and the elapsed time since the last earthquake. Quantitative studies on ancient earthquakes in the east segment of the Eastern Kunlun fault zone began in the 1990s and achieved a succession of research findings. However, the difference between results is relatively large. In this paper, we compare the research results of previous studies and obtain relatively reasonable recurrence intervals for ancient earthquakes and elapsed times since the earthquakes, which provide comparatively accurate source materials for calculating the recurrence probability for each segment.

4.2.1. Maqin Segment

The tectonic geomorphology preserved in the segment at the surface, such as earthquake swells and fault scarps, reflects that the segment had multiple active episodes during the Holocene Epoch. Previous research utilized methods such as excavating trenches and using natural outcrops to identify ancient earthquake events along this segment during different stages (Figure 8). Various researchers [17,24,25,26] determined different periods and times of ancient earthquakes, as shown in Figure 8. Li (2009) established that the most recent age of the buried ancient surface is 358–430 Cal a BP from the footwall (Tc3) of a secondary bench fault profile east of Gequ River, the sampling location of which is an ancient plant layer [18]. However, no records of a relevant earthquake were found after referring to historical documents covering nearly 400 years. Therefore, this age still requires further investigation. It is found from Figure 8 that when the age of an ancient earthquake is 977–1090 Cal a BP, the conclusions of researchers are essentially consistent [17,24] and agree with the time of an earthquake recorded in one passage of the ‘Epic of King Gesar’ historical literature. It is speculated that this age is relatively reliable, and it can be used as the latest elapsed time for the occurrence of an earthquake in this segment. Finally, by comparing and analyzing the main research results, ancient earthquake events of the Maqin segment can be mainly divided into six events: the first is 977–1090 Cal a BP; the second is 2000 ± 300 Cal a BP; the third is 3342–3454 Cal a BP; the fourth is 6600 ± 700 Cal a BP; the fifth is 7971–8050 Cal a BP; and the sixth is 10,000 ± 200 Cal aBP. The calculated average recurrence period of ancient earthquakes is approximately 1765 a.

4.2.2. Maqu Segment

The seismic activity along the Maqu fault is clearly weaker than that along the Maqin fault. Ancient earthquake data of this period were also analyzed intensively by predecessors (Figure 9). There are two trenches near the Xike River and four near Maqu County [2,25,26]. In this paper, we compare the research results of previous studies and derive the elapsed time of earthquakes and the average recurrence interval of earthquakes. By comparing the three profiles of ancient earthquakes near Maqu County, He (2006) used the successive restriction method proposed by Zhang (1996) for ancient earthquake events to obtain comparatively accurate events. Four events occurred during the Holocene Epoch, considering that the first and second times of ancient earthquakes are relatively reliable [12,25]. However, because the number of prospecting trenches is small, the record of ancient earthquakes is still incomplete. Lin (2008) recognized five ancient earthquake events by lamination and sampling of three natural outcrops and two trench sections in the Maqu segment [26]. The location of trench sections among them matches the He (2006) research result, while the obtained age result is also similar. In this paper, by comparing the research results of previous studies, we further combine the sampling positions of outcrops and trench profiles and dating methods to describe seven ancient earthquake events that are relatively reliable: the first is 1210 ± 40 Cal a BP; the second is 1730 ± 50 to 2530 ± 40 Cal a BP; the third is 3736 ± 57 Cal a BP; the fourth is 4850 ± 40 Cal a BP; the fifth is 6100–6700 Cal a BP; the sixth is 8590 ± 70 Cal a BP; and the seventh is 9000–10,000 Cal a BP. We can determine that the past event of the most recent ancient earthquake is 1210 ± 40 Cal a BP, and the average recurrence time period of ancient earthquakes is approximately 1465 a.

4.3. Qualitative Classification of Risk Degrees in Maqin–Maqu Segments

Based on the fact that the average recurrence intervals of earthquakes calculated in this paper are long, we conduct a qualitative classification on the degree of earthquake risk in the Maqin–Maqu segment of the Eastern Kunlun fault zone. We usually use E = t/R to describe the time urgency of earthquake occurrence on the faults, E represents the quantitative parameter of describing the time urgency of earthquake occurrence, t represents the elapsed time since the last earthquake, R represents the average recurrence interval, while, to some extent, we ignore the contribution of R to the time urgency of earthquake occurrence. Even the value of E is the same. In two segments with different R values, the degrees of earthquake risk are different. Generally, the smaller the value of R, the higher is the degree of risk. Therefore, we can use two parameters to divide the degree of risk: (1) it is class C if E < 0.5 and R > 2500 years, which indicates that there is no risk in this section within one thousand years; (2) it is class B if 0.5 < E < 0.8 and R < 2500 years, which indicates that there is a relatively small risk in this segment during the coming several hundred years; and (3) it is class A for 0.8 < E < 1 and R < 2500 years, which indicates that the risk is relatively large in this segment during the coming several hundred years. E takes 0.5, 0.8, and 1.0 as the critical values, which are determined according to the distribution characteristics of earthquake recurrence T/Tave, and R is 2500 years because the recurrence intervals along the two segments are smaller than 2500 years. As the prediction time increases, the increase of risk is also relatively rapid [28]. We divide the segments according to this principle as follows: the Maqin segments belong to class B, and the Maqu segment belongs to class A.

4.4. Probability of Strong Earthquake (Ms ≥ 7) Occurrence in the Next 100 Years

According to the probability distribution function (2) fitted in this paper and the probability distribution function (5) suitable for all active faults in mainland China, we combine the elapsed time and average recurrence interval of the last earthquake for Maqin–Maqu segments in the Eastern Kunlun active fault zone and calculate the conditional probability of earthquake recurrence in the future 20 years, 50 years, and 100 years after 2023 for different segments of this fault zone (Table 2).
In 1984, Wallace et al. [29] analyzed the probability of seismic risk and proposed that recurrence probabilities in the range of 0–10% are considered low; recurrence probabilities in the range of 11–49% are considered intermediate probability events; and recurrence probabilities in the range of 50–100% are considered high probability events. Given the recurrence probabilities calculated in Table 2, the recurrence probability of the Maqu segment in the next 100 years indicates an intermediate probability event, and the Maqin indicates low probability events.
Finally, we find that a comparison result between the probability density function fitted in this study and the conditional probability values calculated using commonly used generic functions reveals that the generic recurrence probability function becomes less sensitive to the dependent variable P (conditional probability) with an increase in the independent variable t/R compared to the fitted function proposed in this paper.

4.5. Joint Probability of the Occurrence of Strong Earthquake with Magnitude Ms ≥ 7 in the Next 100 Years

According to the aforementioned recurrence probability of different fault segments in Table 2 (i.e., the fitting function used in this paper), we calculate the joint probability for the occurrence of at least one large earthquake in the Maqin–Maqu segment of the Eastern Kunlun active fault zone in the future:
P = 1 − (1 − P1) × (1 − P2)…
P1, P2 are the probabilities for the occurrence of an earthquake in the Maqin segment and Maqu segment of the Eastern Kunlun fault.
Table 3 shows that the probability for the occurrence of a strong earthquake in the next 20 years in the east segment of the Eastern Kunlun fault zone is low, and the probabilities for the occurrence of an earthquake in the next 50 and 100 years are intermediate. The probabilities in Table 3 likely represent only the minimum probability values of this active fault zone. Because the data for ancient earthquakes are likely somewhat inadequate, such as the lack of certain active events, they overestimate the recurrence interval and the value of elapsed time.

5. Conclusions

We conducted an analysis of the tectonic background and structural characteristics pertaining to the Maqin–Maqu segment within the East Kunlun Fault Zone. The segment was divided into two seismogenic structures, and an evaluation was performed to determine the potential maximum magnitude of earthquakes that these structures could trigger. Additionally, a recurrence probability model for the region was established by collecting and organizing paleoseismic data from the northeastern part of the northeastern Tibetan Plateau. Subsequently, the model was employed to assess the seismic hazard associated with the Maqin–Maqu segment of the East Kunlun Fault Zone over the next century. The acquired insights and prevailing issues can be summarized as follows:
(1) The calculated recurrence probability values for strong earthquakes in the Maqin–Maqu segment of the Eastern Kunlun fault zone reveal distinct patterns. Over the next 100 years, the Maqu segment presents an intermediate probability of 11.47%, whereas the Maqin segments exhibit values below 10%, signifying events of low probability. Although current research results indicate that the Maqin segment exhibits stronger tectonic activity compared to the Maqu segment, such as a decreasing trend in the sliding rate of the East Kunlun Fault from west to east, the seismic hazard in the next century is closely associated with the time since the last earthquake. The Maqu segment has an elapsed time of around 1200 years, signifying a seismic gap on the East Kunlun Fault throughout the last millennium. Furthermore, it has experienced Coulomb stress loading since the occurrence of the 2008 Wenchuan earthquake. Consequently, there is a critical need to concentrate on monitoring and mitigating the seismic hazard in the Maqu segment.
(2) The sensitivity of conditional probability P, calculated from the strong earthquake recurrence probability model to the prediction time period Δt, is intricately linked with the length of the recurrence cycle R. In general, as R increases, the probability value P experiences a gradual rise with the expanding prediction time interval Δt. However, when R is exceedingly large and Δt is very small, the prediction outcome may be relatively imprecise, underscoring a limitation of the recurrence model in earthquake prediction. The recurrence intervals for different segments in the east segment of the Eastern Kunlun fault zone all hover around one thousand years. While the calculated conditional probabilities are generally modest, they still offer valuable insights into the seismic activity levels of these three segments for several centuries into the future.
(3) We utilized quantitative data on ancient earthquakes pertaining to 22 fault zones in the northeast Tibetan Plateau and subjected it to a normalization process to construct a recurrence probability model for strong earthquakes. A comparative analysis of probability values calculated using the model proposed in our paper and the prevailing universal seismic model reveals that, when the independent variable of the universal model (t/R) approximates 1, the amplitude for the growth of the calculated recurrence probability value (P) is less sensitive compared to the fitting model in this paper. Consequently, for fault zones with a substantial body of studies providing data on ancient earthquakes, it is arguably more accurate to employ the model developed in this paper for assessing the future risk of a large earthquake, particularly in segments characterized by small average recurrence intervals and extended elapsed times. Furthermore, the current universal recurrence model is likely more applicable to faults with a relatively limited number of studies available on ancient earthquakes and segments with relatively long recurrence intervals.

Author Contributions

Z.L.: Formal analysis; Data curation; Drawing diagrams of figures; Writing Original Draft. B.Z.: Writing—Review and Editing; Science and English checking all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key R&D Program project of China (grant number 2022YFC3003501).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available (Table A1). Data analysis was performed using SPSS software (https://www.ibm.com/spss, https://www.ibm.com/docs/en/spss-statistics/beta?topic=features-q-q-plots&mhsrc=ibmsearch_a&mhq=Q-Qplot, accessed on 2 March 2023).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Research Data

Table A1. Data of ancient earthquakes for 22 active fault zones in the northeastern Tibetan Plateau.
Table A1. Data of ancient earthquakes for 22 active fault zones in the northeastern Tibetan Plateau.
No.Name of FaultSituation of SegmentAge of Ancient Earthquake Events (Age)Recurrence Interval (T)Average Recurrence Interval (Tave)Normalization (T/Tave)Data Source
F1Fault on the northern margin of West Qinling MountainsLuanfeng segmentE1: 12,500 ± 500 [30,31]
E2: 7500 ± 5005000 1.2749
E3: 50002500 0.6374
734426639221.0877
Zhangxian segmentE1: 9910–1310
E2: 4250–6930566056601
Guomatan segmentE1: 12,450
E2: 5480 ± 606970 1.1264
1936.2 (74)540661880.8736
F2Haiyuan fault zoneXiaokou-Caixiang fault segment (east segment)E1: 8330–9370 [27,32,33]
E2: 6800–7530800 0.4975
E3: 5690–6100700 0.4353
E4: 590–12604430 2.7550
E5: 192050016080.3109
South and west Huashan fault segmentE1: 8330–9370
E2: 6810–7530800 0.7547
E3: 5640–6100710 0.6698
E4: 4440–5030610 0.5755
E5: 2630–67101810 1.7075
E6: 590–1260137010601.2925
Hasishan-Machangshan fault segment (west segment)E1: 8530–9370
E2: 6150–63502180 1.1315
E3: 3590–41402010 1.0433
E4: 2000159019270.8253
F3Fodongmiao–Hongyazi faultMid-west segment (Q2)E1: >4320 ± 300 [34]
E2: 2100 ± 1001820 1.0344
401169917600.9656
F4Eastern Yumushan FaultWest segmentE1: >10,500 ± 600 [35]
E2: 7200–85001400 0.6432
E3: 3700–39003300 1.5161
180 Fall183021770.8407
F5Eastern Yumushan FaultShanglongwangE1: 10,500 ± 600 [36]
E2: 8600 ± 6001900 0.5672
E3: 3800 ± 600480033501.4328
F6Huangcheng–Shuangta faultShisi segmentE1: 7700 ± 600 [37]
E2: 3400 ± 3004300 1.1292
83331738080.8711
F7Changma fault E1: 9520 ± 90 [38,39]
E2: 9380 ± 150140 0.0593
E3: 5470 ± 603910 1.6564
E4: 3230 ± 552240 0.9490
78315223611.3353
F8Sunan fault segmentMiddle segmentE1: 2200 [40]
E2: 1680520 0.6933
E3: 7009807501.3067
F9Lenglongling faultWest segmentE1: 5926 [41]
E2: 4050 ± 1601876 1.3754
E3: 2900 ± 2701150 0.8431
E4: 1560 ± 3601340 0.9824
470109013640.7991
F10Laohushan-Maomaoshan faultLaohushan segmentE1: 7700 ± 50 [42,43]
E2: 6180 ± 1501520 1.5833
E3: 5200 ± 100980 1.0208
E4: 4250 ± 150950 0.9896
E5: 3050 ± 1501200 1.2500
E6: 2000 ± 3001050 1.0938
E7: 800 ± 1001200 1.2500
122678 0.7063
201029600.1063
Maomaoshan segmentE1: 9000 ± 300
E2: 6600 ± 5002400 1.3333
E3: 5000 ± 3001600 0.8889
E4: 3700 ± 3001300 0.7222
E5: 1800 ± 300190018001.0556
F11Tianqiaogou-Huangyangchuan faultTianqiaogouE1: 27,700 ± 2200 [44]
E2: 21,300 ± 24006400 1.3904
E3: 16,800 ± 14004500 0.9776
E4: 13,700 ± 20003100 0.6735
E5: 10,300 ± 17003400 0.7386
E6: 7590 ± 1002710 0.5887
83750746031.6309
F12Zhongwei–Tongxin faultDaquanshui-Hushanzi segment (west segment)E1: 13,150 ± 800 [45]
E3: 6535 ± 2004584 1.1293
E5: 3000353540590.8709
Xiliangtou-Shuangjingzi segment (Middle segment)E1: 13,150 ± 800
E2: 8566 ± 5004584 1.0703
E4: 5450 ± 2003116 0.7275
301514942831.2022
F13Bayinguole River North Margin fault E1: 32,700 ± 1450 [46]
E2: 15,540 ± 132017,160 1.1652
E3: 3245 ± 33012,29514,7270.8349
F14Elashan fault E1: 12,500 ± 100 [47]
E2: 10,000 ± 1502500 1.0101
E3: 6000 ± 1004000 1.6162
E4: 4100 ± 3001900 0.7677
E5: 2600 ± 400150024750.6061
F15Eastern Kunlun fault zoneKusaihu segment16,865 ± 1018 [18,20,27,48,49,50]
12,935 ± 7743930 1.1158
9730 ± 5923205 0.9100
6955 ± 4252775 0.7879
3100 ± 2013855 1.0945
2001.11.14 (9)384635221.0920
East and west Datan segmentE1: 11,000
E2: 70004000 1.2001
E3: 27004300 1.2901
E4: 1000170033330.5101
Xiugou-Tuosuohu segmentE1: 8000
E2: 60002000 1.0096
E3: 10005000 2.5240
E4: 637363 0.1832
7356419810.2847
Tuosuohu-Xiadawu segmentE1: 10,000
E2: 6100–67003900 1.0400
E3: 2500360037500.9600
Maqin segment9900 ± 200
7971–80501850 1.2107
7200 ± 800771 0.5046
3342–34543746 2.4516
20001342 0.8783
977–109091015280.5955
Maqu segment15,800 ± 2500–24,100 ± 2900
9000–10,0005800 2.9532
7460 ± 60–8690 ± 401540 0.7841
4586 ± 124–7460 ± 602874 1.4633
3736 ± 57–4586 ± 24850 0.4328
1730 ± 50–2530 ± 402006 1.0214
1210 ± 50–1730 ± 50520 0.2648
1055–152415519640.0789
F16Altun fault zoneAndier River-Cheerchen River segmentE1: 3800 [51,52]
E2: 2760–29001040 0.9020
E3: 10771683 1.4597
E4: 34273511530.6375
Huangtuquan-Wuzunxiaoer segmentE1: 6366
E2: 2270–3500286628661.0000
Suoerkuli-Akesai segmentE1: 3900–4130
E2: 1800–21702100 1.4094
E3: 92088014900.5906
E1: 16,400 ± 100
E2: 13,500 ± 180290029001.0000
Subei-Kuantanshan segmentE1: 18,620–18,780
E2: 11,330–12,5906030 1.2048
E3: 7080–7350398050050.7952
F17Maduo-Gande faultGande Segment (Q3–Q4)E1: 4000–4700 [53]
E2: 3000–32301000 0.8333
E3: 1600140012001.1667
F18Tuolaishan faultLiuhuanggou-Youhuli segment (east segment)E1: 13,240–13,960 [54]
E2: >6080 ± 4506710 2.6733
E3: 3570–4390169025100.6733
F19Changma-Ebo faultChangma fault segment (Q4)E1: 12,000–13,000 [55]
E2: 9520 ± 902480 1.3333
E3: 8020–8700820 0.4409
E4: >6670 ± 801270 0.6828
E5: >5470 ± 601200 0.6452
E6: 32302240 1.2043
78315218601.6946
F20Yellow River-Linwu faultLinwu segmentE1: 27,000 ± 7800 [56]
E2: 20,0007000 1.3333
E3: 13,000 ± 6007000 1.3333
E4: 10,600 ± 5002400 0.4571
E5: 6000460052500.8762
F21Luoshan east faultE1: 8200 ± 600 [57]
E2: 5000 ± 703200 1.2384
E3: 3300 ± 12 01700 0.6579
449285125841.1033
F22Helanshan east faultE1: 8400 [58]
E2: 57002700 0.9963
E3: 26003100 1.1439
271232927100.8594

References

  1. Deng, Q.; Ran, Y.; Yang, X. Active Tectonics Map of China (1:4,000,000); Seismological Press: Beijing, China, 2007. (In Chinese) [Google Scholar]
  2. Seismological Bureau of Qinghai Province. East Kunlun Active Fault Zone; Institute of Crustal Deformarion, China Seismological Bureau, Seismological Press: Beijing, China, 1999; pp. 1–186. (In Chinese) [Google Scholar]
  3. Zhang, Z.; Klemperer, S.; Bai, Z.; Chen, Y.; Teng, J. Crustal structure of the Paleozoic Kunlun orogeny from an active-source seismic profile between Moba and Guide in East Tibet, China. Gondwana Res. 2011, 19, 994–1007. [Google Scholar] [CrossRef]
  4. Karplus, M.S.; Zhao, W.; Klemperer, S.L. Injection of Tibetan crust beneath the south Qaidam Basin: Evidence from INDEPTH IV wide-angle seismic data. J. Geophys. Res.-Solid Earth 2011, 116, 1–23. [Google Scholar] [CrossRef]
  5. Parsons, T.; Chen, J.; Kirby, E. Sress changes from the 2008 wenchuan earthquake and increased hazard in the Sichuan basin. Nature 2008, 454, 509–510. [Google Scholar] [CrossRef] [PubMed]
  6. Sykes, L.R.; Nishenko, S.P. Probabilities of occurrence of large plate rupturing earthquakes for San Andreas, San Jacinto and Imperial fault, California, 1983~2003. J. Geophys. Res. 1984, 89, 5905–5927. [Google Scholar] [CrossRef]
  7. Nishenko, S.P. Seismic potential for large and great interplate earthquakes along the Chilean and southern Peruvian margins of South America: A quantitative reappraisal. J. Geophys. Res. 1985, 90, 3589–3615. [Google Scholar] [CrossRef]
  8. Nishenko, S.P.; Buland, R. A generic recurrence interval distribution for earthquake forecasting. Bull. Seismol. Soc. Am. 1987, 77, 1382–1399. [Google Scholar]
  9. Nishenko, S.P.; Singh, S.K. Conditional probabilities for the recurrence of large and great iner-plate earthquakes along the Mexican Subduction zone. Bull. Seismol. Soc. Am. 1987, 77, 2094–2114. [Google Scholar]
  10. Wen, X.Z. Conditional Probabilities for the Recurrence of Earthquakes on the Xianshuihe Fault Zone within the Coming Three Decades. Earthq. Res. China 1990, 6, 8–14. (In Chinese) [Google Scholar]
  11. Liu, J.; Wang, L.M. Evaluate the middle-and longterm seismic risk in the Fen-Wei seismic belt by using the active fault data. Acta Seismol. Sin. 1999, 18, 427–435. (In Chinese) [Google Scholar] [CrossRef]
  12. Zhang, P.Z.; Mao, F.Y. Active Faulting and Fault Specific Probabilistic Seismic Hazard Assessment. In The Active Fault in China; Seismological Press: Beijing, China, 1996; pp. 12–31. (In Chinese) [Google Scholar]
  13. Wen, Y.; Fang, Z.; He, K.; Yang, J.; Xiong, L.; Xu, C. Present-day crustal deformation of the Maqin-Maqu segment in the East Kunlun fault zone revealed by Sentinel-1 images. Chin. J. Geophys. 2023, 66, 4517–4532. (In Chinese) [Google Scholar] [CrossRef]
  14. Elliott, A.J.; Oskin, M.E.; Liu, Z.J.; Shao, Y. Rupture termination at restraining bends: The last great earthquake on the Altyn Tagh fault. Geophys. Res. Lett. 2015, 42, 2164–2170. [Google Scholar] [CrossRef]
  15. Duan, B. Multicycle dynamics of the Aksay Bend along the Altyn Tagh fault in Northwest China: 1. A simplified double bend. Tectonics 2019, 38, 1101–1119. [Google Scholar] [CrossRef]
  16. Xu, X.; Wang, F.; Zheng, R.; Chen, W.; Ma, W.; Yu, G.; Chen, G.; Tapponnier, P.; Van Der Woerd, J.; Meriaux, A.S.; et al. Late Quaternary sinistral slip rate along the Altyn Tagh fault and its structural transformation model. Sci. China Ser. D Earth Sci. 2005, 48, 384–397. [Google Scholar] [CrossRef]
  17. Li, C.X. The Long-Term Faulting Behavior of the Eastern Segment (Maqin-Maqu) of the East Kunlun Fault since the Late Quaternary. Ph.D. Thesis, Institute of Geology, China Earthquake Administration, Beijing, China, 2009. [Google Scholar]
  18. Harkins, N.; Kirby, E.; Shi, X.; Wang, E.; Burbank, D.; Chun, F. Millennial slip rates along the eastern Kunlun fault: Implications for the dynamics of intracontinental deformation in Asia. Lithosphere 2010, 2, 247–266. [Google Scholar] [CrossRef]
  19. Ren, J.; Xu, X.; Yeats, R.S.; Zhang, S. Millennial slip rates of the Tazang fault, the eastern termination of Kunlun fault: Implications for strain partitioning in eastern Tibet. Tectonophysics 2013, 608, 1180–1200. [Google Scholar] [CrossRef]
  20. Makkonen, L. Bringing closure to the plotting position controversy. Commun. Stat. Theory Methods 2008, 37, 460–467. [Google Scholar] [CrossRef]
  21. Wilk, M.B.; Gnanadesikan, R. Probability plotting methods for the analysis of data. Biometrika 1968, 55, 1–17. [Google Scholar] [CrossRef] [PubMed]
  22. Wells, D.L.; Coppersmith, K.J. Updated empirical relationships among magnitude, rupture length, rupture area, and surface displacement. Bull. Seism. Soc. Am. 1994, 84, 947–1002. [Google Scholar] [CrossRef]
  23. Ran, H.L. Research on uncertainty of upper limit earthquake magnitude in potential seismic source zone. Acta Seimol. Sin. 2009, 31, 396–402. (In Chinese) [Google Scholar]
  24. Li, C.F.; He, Q.L.; Zhao, G.G. Paleo-earthquake studies on the eastern section of the kunlun fault. Acta Seismol. Sin. 2005, 27, 60–67. (In Chinese) [Google Scholar] [CrossRef]
  25. He, W.G.; Xiong, Z.; Yuan, D.Y. Palaeo-earthquake Study on the Maqu Fault of East Kunlun Active Fault. Earthq. Res. China 2006, 22, 126–134. (In Chinese) [Google Scholar]
  26. Lin, A.M.; Guo, J.M. Non-uniform Slip Rate and Millennial Recurrence Interval of Large Earthquakes along the Eastern Segment of the Kunlun Fault, Northern Tibet. Bull. Seismol. Soc. Am. 2008, 98, 2866–2878. [Google Scholar] [CrossRef]
  27. Institute of Geology, State Seismological Bureau, Seismological Bureau of Ningxia Hui Autonomous Region. Research on the Active Haiyuan Fault Zone; Seismological Press: Beijing, China, 1990. [Google Scholar]
  28. Matsuda, T. Estimation of future destructive earthquakes from active faults on land in Japan. J. Phys. Earth Suppl. 1977, 25, 795–855. [Google Scholar] [CrossRef] [PubMed]
  29. Wallace, R.E. Terms for expressing earthquake potential prediction, and probability. Bull. Seism. Soc. Am. 1984, 74, 1819–1925. [Google Scholar]
  30. Teng, R.Z.; Jin, Y.Q.; Li, X.H. Characteristics of recent activity along the fault zone at the northern edge of western Qinling Mt. Northwestern Seismol. J. 1994, 1, 85–90. [Google Scholar]
  31. Cao, J.J.; Liu, B.Z.; Wen, X.Z. Determination of the Average Recurrence Intervals of Characteristic Earthquake and Estimate of Earthquake Rick on Northern Xiqinling Faults. Earthq. Res. 2003, 26, 372–381. [Google Scholar]
  32. Zhang, P.Z.; Min, W.; Deng, Q.D. Paleoearthquake and strong earthquake recurrence interval in Haiyuan Active Fault. Sci. China Ser. D 2003, 33, 705–713. [Google Scholar]
  33. Ou, Q.; Kulikova, G.; Yu, J.; Elliott, A.; Parsons, B.; Walker, R. Magnitude of the 1920 Haiyuan Earthquake Reestimated Using Seismological and Geomorphological Methods. J. Geophys. Res. Solid Earth 2020, 125, 8. [Google Scholar] [CrossRef]
  34. Zhen, W.J. Geometric Pattern and Active Tectonics of the Hexi Corridor and Its Adjacent Regions. Ph.D. Thesis, Institute of Geology, China Earthquake Administration, Beijing, China, 2009. [Google Scholar]
  35. Chen, B.L.; Wang, C.Y.; Gong, H.L. A new understanding of the characteristics of Late Quaternary activity of the northern Yumushan marginal fault in the Hexi corridor, northwestern China. Geol. Bull. China 2007, 26, 976–983. [Google Scholar]
  36. Zhou, J.S.; Mao, S.M.; Chen, Z.T. Preliminary study on the Late Quaternary Active feature of The Yumushan Eastern Marginal Fault Zone. Earthq. Res. Plateau 1993, 5, 31–35. [Google Scholar]
  37. Wang, Y.C.; Liu, B.C. Analysis on Seismic Risk for Faults in The Mid-Esstern Qilianshan Area. Northwestern Seismol. J. 2001, 23, 330–338. [Google Scholar]
  38. Kang, L.X. Analysis of paleoearthquake on Changma Fault. J. Seismol. 1986, 4, 16–21. [Google Scholar]
  39. Guo, J.X.; Hou, Z.Q.; Hou, K.M. Active Character and Palaeo-earthquake on median Section of Changma-Qilian Fault Zone in Holocene. Northwestern Seismol. J. 1990, 12, 38–43. [Google Scholar]
  40. Liu, B.C.; Cao, J.J.; Yuan, D.Y. Quantitative data of Acitive Faults within the Active Tectonic Block in North Qinghai-xizang Plateau. Seismol. Geol. 2008, 30, 161–175. [Google Scholar]
  41. Hen, W.G.; Liu, B.J.; Yuan, D.Y. Preliminary research on the paleoearthquakes along the Leng Long Ling active fault zone. In Research on Active Fault (8); Institute of Geology, China Seismological Bureau, Ed.; Seismological Press: Beijing, China, 2006; pp. 64–74. (In Chinese) [Google Scholar]
  42. Institute of Geology, State Seismological Burean, Lanzhou Institute of Seismology, State Seismological Bureau. The Qilian Moutain-Hexi Corridou Active Fault System; Seismological Press: Beijing, China, 1993. (In Chinese) [Google Scholar]
  43. Liu, X.F.; Liu, B.Z.; Lv, T.Y. The research on the Laohushan active fault South China. J. Seismol. 1994, 14, 9–16. [Google Scholar]
  44. Zen, W.J.; Yuan, D.Y. Characteristics of Palaeoearthquake Activity along The Active Tianqiaogou Huangyangchuan Fault on The Eastern Section of The Qilanshan mountains. Seismol. Geol. 2004, 26, 645–657. [Google Scholar]
  45. Min, W.; Zhang, P.Z.; Deng, Q.D. Study of Holocene paleoearthquake on Zhongwei-Tongxin Fault zone. Seismol. Geol. 2001, 23, 357–366. [Google Scholar]
  46. Liu, X.L.; Yuan, D.Y. Study on The Active features of Bsy Guole River Active Fault, Delingha, Qinghai. Northwestern Seismol. J. 2004, 26, 303–308. [Google Scholar]
  47. Yuan, D.Y.; Zhang, P.Z.; Liu, X.L. The tectonic activity and deformation features during the Late Quaternary of ElashanMountain active fault zone in Qinghai Province and its implication for the deformationof the northeastern margins of the Qinghai-Tibet Platea. Earth Sci. Front. 2004, 11, 393–402. [Google Scholar]
  48. Hu, D.; Wu, Z.; Wu, Z.; Zhao, X.; Ye, P. Late quaternary paleoseismic history on the kusaihu segment of East Kunlun fault zone in northern Tibet. Quat. Sci. 2007, 27, 27–34. (In Chinese) [Google Scholar]
  49. Hu, D.; Ye, P.; Wu, Z.; Wu, Z.; Zhao, X.; Liu, Q. Research onHolocene paleoearthquakes on the Xidatan segment of the EastKunlun fault zone in Northern Tibet. Quat. Sci. 2006, 26, 1012–1020. (In Chinese) [Google Scholar]
  50. Ren, J.; Xu, X.; Yeats, R.S.; Zhang, S.; Ding, R.; Gong, Z. Holocene pa-leoearthquakes of the Maoergai fault, eastern Tibet. Tectonophysics 2013, 590, 121–135. [Google Scholar] [CrossRef]
  51. Task Group of Altun Active Fault Zone; State Seismological Bureau. Altun Active Fault Zone; Seismological Press: Beijing, China, 1992. (In Chinese) [Google Scholar]
  52. Shao, Y.; Zou, X.; Yuan, D.; Yao, Y.; Liu, X. Late Quaternary Slip along Yangguan Fault at Northeastern Section of Altyn Tagh Fault and Implications for Seismic Risk. Earth Sci. 2021, 46, 683–696. [Google Scholar] [CrossRef]
  53. Yu, W.; Chen, Z.; Chen, B.; Wang, Y.; Meng, L.; Sun, Y.; Zhang, W. Early Paleozoic Tectonic Deformation in Qiashenkansayigou Area, North Altun, and Implication for Tectonic Evolution. J. Geomech. 2019, 25, 301–312, (In Chinese with English abstract). [Google Scholar]
  54. Xiong, R.W.; Ren, J.W.; Zhang, J.L. Late Quaternary Active Characteristics of the Gande Segment in the Maduo-Gande Fault Zone. Earthquake 2010, 30, 65–73. [Google Scholar]
  55. The Research Group on Active Fault System around Ordos Massif, State Seismological Bureau. Active Fault System around Ordos Massif; Seismological Press: Beijing, China, 1988. (In Chinese) [Google Scholar]
  56. Cai, Z.Z.; Liao, Y.H.; Zhang, W.J. Late Quaternary Paleoearthquakes and Their Rupture Features along the Lingwu Fault. Seismol. Geol. 2001, 23, 15–23. [Google Scholar]
  57. Min, W.; Chai, Z.Z.; Wang, P. The study on the paleoearthquakes on the eastern pienmont fault of the luoshan mountains in Holocene. Earthq. Res. Plateau 1993, 5, 97–102. [Google Scholar]
  58. Du, P.; Chai, Z.Z.; Liao, Y.H. Study on Holocene Acticity of The South Segment of The Eastern Piedmont fault of the South Segment of The Eastern Piedmont Fault of Helan mountains between Taomengou and Yushugou. Seismol. Geol. 2009, 31, 256–264. [Google Scholar]
Figure 1. Topographic map of the northern Tibetan Plateau region. Major active faults are adapted from Deng et al., 2007 [1]. The upper-left figure shows the location of Northern Tibet. The Blue dashed box indicates the location of Figure 2. The Eastern Kunlun Fault zone is represented by bold colorful lines, with historic earthquake surface ruptures marked by red or yellow lines. Purple lines represent the faults for which paleoseismology records are used in this study. Detailed earthquake timing for these faults is listed in Table A1. Fault name abbreviations: BYF, North Bayinguole River fault; EHF, Eastern Helanshan piedmont fault; ELF, Eastern Luoshan piedmont fault; EYF, Eastern Yumushan piedmont fault; FHF, Fodongmiao–Hongyazi Fault; HSF, Huangcheng–Shuangta fault; THF, Tianqiaogou–Huangyangchuan fault; TLSF, Tuolaishan fault; YF, Yumushan fault; ZTF, Zhongwei–Tongxin fault.
Figure 1. Topographic map of the northern Tibetan Plateau region. Major active faults are adapted from Deng et al., 2007 [1]. The upper-left figure shows the location of Northern Tibet. The Blue dashed box indicates the location of Figure 2. The Eastern Kunlun Fault zone is represented by bold colorful lines, with historic earthquake surface ruptures marked by red or yellow lines. Purple lines represent the faults for which paleoseismology records are used in this study. Detailed earthquake timing for these faults is listed in Table A1. Fault name abbreviations: BYF, North Bayinguole River fault; EHF, Eastern Helanshan piedmont fault; ELF, Eastern Luoshan piedmont fault; EYF, Eastern Yumushan piedmont fault; FHF, Fodongmiao–Hongyazi Fault; HSF, Huangcheng–Shuangta fault; THF, Tianqiaogou–Huangyangchuan fault; TLSF, Tuolaishan fault; YF, Yumushan fault; ZTF, Zhongwei–Tongxin fault.
Applsci 14 02691 g001
Figure 4. Histogram of normalized data of ancient earthquake recurrence intervals in the northeastern Tibetan Plateau T/Tave and the normal distribution fitting curve. T represents recurrence interval of ancient earthquake, Tave represents average recurrence interval, T/Tave represents normalized data of ancient earthquake recurrence intervals.
Figure 4. Histogram of normalized data of ancient earthquake recurrence intervals in the northeastern Tibetan Plateau T/Tave and the normal distribution fitting curve. T represents recurrence interval of ancient earthquake, Tave represents average recurrence interval, T/Tave represents normalized data of ancient earthquake recurrence intervals.
Applsci 14 02691 g004
Figure 5. Diagram of normal Q-Q plot and reverse trend normal Q-Q plot. A normal Q-Q plot comparing randomly generated, independent standard normal data on the vertical axis to a standard normal population on the horizontal axis. The linearity of the points suggests that the data are normally distributed.
Figure 5. Diagram of normal Q-Q plot and reverse trend normal Q-Q plot. A normal Q-Q plot comparing randomly generated, independent standard normal data on the vertical axis to a standard normal population on the horizontal axis. The linearity of the points suggests that the data are normally distributed.
Applsci 14 02691 g005
Figure 6. Diagram of logarithmic normal Q-Q plot and reverse trend logarithmic normal Q-Q plot. A Q-Q plot of a sample of data versus a logarithmic distribution. The deciles of the distributions are shown in the plot. Many outliers are evident at the high end of the range.
Figure 6. Diagram of logarithmic normal Q-Q plot and reverse trend logarithmic normal Q-Q plot. A Q-Q plot of a sample of data versus a logarithmic distribution. The deciles of the distributions are shown in the plot. Many outliers are evident at the high end of the range.
Applsci 14 02691 g006
Figure 7. Diagram of Weibull distribution Q-Q plot and reverse trend Weibull distribution Q-Q plot. A Q-Q plot of a sample of data versus a Weibull distribution. The deciles of the distributions are shown in the plot. many outliers are evident at the high end of the range.
Figure 7. Diagram of Weibull distribution Q-Q plot and reverse trend Weibull distribution Q-Q plot. A Q-Q plot of a sample of data versus a Weibull distribution. The deciles of the distributions are shown in the plot. many outliers are evident at the high end of the range.
Applsci 14 02691 g007
Figure 8. (A) Distribution map of the Maqin segment and adjacent segments of the Eastern Kunlun fault. (B) The occurrence ages of paleoearthquakes since 12 ka are compiled [25,26,27].
Figure 8. (A) Distribution map of the Maqin segment and adjacent segments of the Eastern Kunlun fault. (B) The occurrence ages of paleoearthquakes since 12 ka are compiled [25,26,27].
Applsci 14 02691 g008
Figure 9. (A) Distribution map of the Maqu segment and adjacent segments of the Eastern Kunlun fault. (B) The occurrence ages of paleoearthquakes since 12 ka are compiled [25,26,27].
Figure 9. (A) Distribution map of the Maqu segment and adjacent segments of the Eastern Kunlun fault. (B) The occurrence ages of paleoearthquakes since 12 ka are compiled [25,26,27].
Applsci 14 02691 g009
Table 1. Estimate of the magnitude of Maqin and Maqu Segment for single or joint rupture.
Table 1. Estimate of the magnitude of Maqin and Maqu Segment for single or joint rupture.
Name of SegmentLength of Segment (km)Boundary Mark of Segment Division (with West Segment)Type of RuptureEquation of Statistical RelationshipEstimate of the Magnitude
(Mw)
Estimate of the Magnitude
(Ms)
Maqin segment132Animaqing Mountain upliftSingle-segment ruptureMw = 5.16 + 1.12LgL
[22]
Ms = 1.412 + 0.845Mw
[23]
7.57.7
Maqu segment180Triple junction (fault intersection decomposition zone)Single-segment rupture7.67.8
Maqin–Maqu segment 312Animaqing Mountain upliftTwo-segment joint rupture7.98.1
Table 2. Conditional probability of strong earthquake recurrence in the future 20 years, 50 years, and 100 years for Maqin–Maqu segment of the Eastern Kunlun active fault zone (calculated from 2023).
Table 2. Conditional probability of strong earthquake recurrence in the future 20 years, 50 years, and 100 years for Maqin–Maqu segment of the Eastern Kunlun active fault zone (calculated from 2023).
Name of SegmentElapsed Time since Last Earthquake (t)Average Recurrence Interval of Ancient Earthquakes (R)Qualitative Classification of RiskRelationship of Probability Distribution Function(Ms ≥ 7) Conditional Probability P (%)
2023–2043
(Δt = 20)
2023–2073
(Δt = 50)
2023–2123
(Δt = 100)
Maqin977–1090a1765 aB(2)0.82–1.00%2.09–2.56%4.33–5.27%
(5)1.11–1.19%2.65–2.72%5.19–5.34%
Maqu1210 ± 40a1465 aA(2)2.26%5.70%11.47%
(5)1.192.955.77
Table 3. Joint probability for the occurrence of at least one large earthquake in Maqin–Tazang segment of the Eastern Kunlun fault zone during the future 20 years, 50 years, and 100 years (calculated from 2023).
Table 3. Joint probability for the occurrence of at least one large earthquake in Maqin–Tazang segment of the Eastern Kunlun fault zone during the future 20 years, 50 years, and 100 years (calculated from 2023).
Name of SegmentMaqin–Maqu Segment of the Eastern Kunlun Fault
Future time interval/a2023–2043
(Δt = 20)
2023–2073
(Δt = 50)
2023–2123
(Δt = 100)
Joint probability3.24%8.11%16.14%
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

Li, Z.; Zhou, B. Assessment of Strong Earthquake Risk in Maqin–Maqu Segment of the Eastern Kunlun Fault, Northeast Tibet Plateau. Appl. Sci. 2024, 14, 2691. https://doi.org/10.3390/app14072691

AMA Style

Li Z, Zhou B. Assessment of Strong Earthquake Risk in Maqin–Maqu Segment of the Eastern Kunlun Fault, Northeast Tibet Plateau. Applied Sciences. 2024; 14(7):2691. https://doi.org/10.3390/app14072691

Chicago/Turabian Style

Li, Zhengfang, and Bengang Zhou. 2024. "Assessment of Strong Earthquake Risk in Maqin–Maqu Segment of the Eastern Kunlun Fault, Northeast Tibet Plateau" Applied Sciences 14, no. 7: 2691. https://doi.org/10.3390/app14072691

APA Style

Li, Z., & Zhou, B. (2024). Assessment of Strong Earthquake Risk in Maqin–Maqu Segment of the Eastern Kunlun Fault, Northeast Tibet Plateau. Applied Sciences, 14(7), 2691. https://doi.org/10.3390/app14072691

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