The Xianshuihe Fault Zone is one of the most historically seismically active regions in mainland China. However, the seismicity along this fault zone has been quiescent for the past 40 years, since the Daofu M6.9 earthquake in 1981. Understanding its current deformation patterns and fault coupling characteristics is of great significance to estimate the potential risk of strong earthquakes. In this study, we analyzed the dynamic deformation and fault coupling characteristics along the Xianshuihe Fault Zone using Global Positioning System (GPS) data for 1999–2007 and 2016–2020. The results show that the deformation pattern of the Xianshuihe fault zone underwent a dynamic adjustment after the Wenchuan and Lushan earthquakes, i.e., the maximum shear strain accumulation rates of the Luhuo and Daofu sections significantly decreased from 6.0×10-8/a to 3.2×10-8/a, while that of the southeastern segment (i.e., Kangding and Moxi sections) increased from 4.5×10-8/a to 6.2×10-8/a. Additionally, the slip rate and deformation width of the Xianshuihe Fault Zone also changed during these two periods. Combined with the near-field cross-fault observation data, we suggest that the surrounding strong earthquakes 2008 Wenchuan Mw7.9 and 2013 Lushan Mw6.6 had evident differential impacts on the deformation pattern of the Xianshuihe Fault Zone. The fault-coupling inversion results show that the locking degree of the Xianshuihe Fault Zone continued to increase after the Mw7.9 Wenchuan and Mw6.6 Lushan earthquakes, especially the Qianning and Moxi sections increased significantly, with an average coupling coefficient of greater than 0.9 and left-lateral slip-rate deficits of ~5 mm/a and ~8 mm/a, respectively. In contrast, the locking degree of the Kangding section decreased with almost no slip-rate deficit, which may be due to the partial energy release caused by the Mw5.9 and Mw5.6 Kangding earthquakes in 2014. The analysis of the recent rupture history and strain accumulation characteristics of the Xianshuihe Fault Zone indicates that both the Qianning and Moxi sections have a high seismic potential for the next strong earthquake in the Xianshuihe Fault Zone.
Research Article
Dynamic Deformation And Fault Locking of The Xianshuihe Fault Zone, Southeastern Tibetan Plateau: Implications For Seismic Hazards
https://doi.org/10.21203/rs.3.rs-783486/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 22 Feb, 2022
You are reading this latest preprint version
Xianshuihe Fault Zone
Strain rate accumulation
Fault locking degree
Slip-rate deficit
Seismic hazards
The Xianshuihe Fault Zone, located in the southeastern margin of the Tibetan Plateau, is a block boundary fault that plays a significant role in accommodating the Late Cenozoic crustal deformation of the eastern Tibetan Plateau (Wang et al., 1998; Tapponnier et al., 1982; Avouac and Tapponnier, 1993). This fault is characterized by a high left-lateral strike-slip rate and frequent seismic activity. Neotectonic and geodetic studies have estimated that the left-lateral slip on the northwestern and southeastern segments of the Xianshuihe Fault are approximately 15 ± 5 mm/a and 5–9 mm/a (Allen et al., 1991; Gan et al., 2007; Wang et al., 2009; Wang and Shen, 2020; Ji et al., 2020). In the past ~ 300 years, at least eight M ≥ 7.0 earthquakes and fifteen M ≥ 6.5 earthquakes initiated in the Xianshuihe Fault Zone (Wen et al., 2008; Yi et al., 2015). However, the strong seismic activity on the Xianshuihe Fault Zone has been relatively quiescent for ~ 40 years, with only two moderate earthquakes, the Mw5.9 and Mw5.6 Kangding earthquakes in 2014. In contrast, a series of strong earthquakes have successively occurred along the other boundary faults around the Bayan Har Block, including the Mw7.8 Kunlun earthquake in 2001, the Mw7.9 Wenchuan earthquake in 2008, and the Mw6.9 Yushu earthquake in 2010 (Deng et al., 2015). Large earthquakes can modulate the crustal deformation, stress field, and strain accumulation of the surrounding faults (King et al., 1994; Stein, 2003; Nalbant and McCloskey, 2011; Wu et al., 2018). Therefore, the potential seismic hazard posed by the Xianshuihe Fault deserves additional attention.
The strain rate accumulation, fault coupling, and slip-rate deficit provide insights into the seismic hazard assessment on faults (Miura et al., 2002; McCaffrey et al., 2002; Ader et al., 2012; Riguzzi et al., 2012; Jiang et al., 2015; Dogru et al., 2018). Based on joint inversions of Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data, Jiang et al. (2015) inverted the distribution of the interseismic fault locking of the Xianshuihe–Zemuhe Fault system and found that a potential asperity exists in the Qianning–Kangding section. Using GPS data from 1999 to 2017, Guo et al. (2018) inverted the slip rate and the locking degree of the Xianshuihe Fault Zone based on a screw dislocation model (Weertman et al., 1966; Wang et al., 2009) and suggested that the Daofu–Kangding section is at high risk for Mw6.6 or greater earthquakes in the future. Using GPS data from 2004 to 2017, Li et al. (2020) also inverted the slip deficit rate and locking degree of the Xianshuihe Fault Zone based on a two-dimensional elastic dislocation model (Savage et al., 1999; Segall, 2010). Their results revealed that the Daofu–Qianning section is experiencing creep at the surface, and it will not be able to generate earthquakes with magnitudes of larger than 7.0 in the near future. Moreover, they determined that the Selaha Fault is most likely to generate a Mw > 7 earthquake. Based on 2014–2019 InSAR images, Qiao et al. (2021) estimated the deformation field and determined the slip distribution of the Xianshuihe Fault Zone. They concluded that the Qianning section is capable of generating a Mw6.6 earthquake. Although these studies agree that the Xianshuihe Fault Zone presents a high seismic hazard, the potential location of the strong earthquakes in the Xianshuihe Fault Zone remains controversial.
The dynamic changes in the crustal movement obtained via geodetic observations not only respond to the rupture process of great earthquakes but also reflect the seismogenic processes of strong earthquakes and the adjustment process after an earthquake, which is an essential information for determining the subsequent seismic potentials of large earthquakes (Jiang et al., 2009; Zou et al., 2015; Wu et al., 2018; Zhao et al., 2020). In this study, we focused on the dynamic deformation characteristics of the Xianshuihe Fault Zone using two periods of GPS observations and estimated the potential seismic hazards it poses. First, we analyzed the characteristics of the historical seismic activity along the Xianshuihe Fault Zone, and then, we analyzed the dynamic deformation characteristics by calculating the strain rate fields and GPS velocity profiles. Finally, based on the negative dislocation model, we inverted the fault coupling and slip-rate deficit distribution to provide a detailed estimation of the potential seismic hazards posed by the Xianshuihe Fault Zone.
2.1 Tectonic setting
The Xianshuihe Fault Zone is a highly active left-lateral strike-slip boundary fault between the Bayan Har Block and the Sichuan-Yunnan Block within the eastern Tibetan Plateau (Zhang et al., 2003). The Xianshuihe Fault is approximately 400 km long, strikes 40°N–50°W, and has a steep dip angle. It extends from Donggu, Ganzi County in the northwest through Luhuo County, Daofu County, and Kangding County to Shimian County in the south (Qian et al., 1988; Wen et al., 1989; Xu et al., 2003). The Xianshuihe Fault Zone is divided into northwestern and southeastern segments by a geometric discontinuity near the Qianning Basin (Qian et al., 1988). The northwestern segment has a simple fault geometry and consists of the Luhuo, Daofu, and Qianning sections, whereas the structure of the southeastern segment is relatively complicated (Bai et al., 2018). The Kangding section, which extends from the Qianning Basin to Kangding County, consists of three nearly parallel secondary faults, the Yala River, Selaha, and Zheduotang branches. The Moxi section extends from south of Kangding to Shimian, and the surface trace of this fault is distinct (Fig. 1). The detailed geometry of each fault section is presented in Table 1.
|
Number |
Fault section |
Strike |
Dip |
|---|---|---|---|
|
S1 |
Luhuo |
N45°W |
NE80° |
|
S2 |
Daofu |
N45°W |
NE80° |
|
S3 |
Qianning |
N40°W |
NE80° |
|
S4 |
Kangding |
N30°W |
NE70° |
|
S5 |
Moxi |
N20°W |
NE70° |
2.2 Historic seismic activity
The Xianshuihe Fault Zone has experienced numerous historic strong earthquakes. Here, we compiled a catalog of the M ≥ 6.0 earthquakes during 1700–2020 from two sources, including records of historical earthquakes in China and data from the China National Seismic Network. Earthquake rupture data from Wen et al. (2008) were also collected. The magnitude–time frequency diagram and spatiotemporal distribution of the strong earthquake (M ≥ 6.0) ruptures in the Xianshuihe Fault Zone are shown in Fig. 2. We identified two seismically active periods, 1700–1816 and 1893–present (Fig. 2a). During 1700–1816, at least nine M ≥ 6.0 earthquakes occurred in the Xianshuihe Fault Zone, including three M ≥ 7.0 earthquakes, which led to the rupture of most of the sections of the Xianshuihe Fault Zone (Fig. 2b). During the second active period, i.e., since 1893, seven M ≥ 6.5 earthquakes occurred in the Xianshuihe Fault Zone, including five M ≥ 7.0 earthquakes. However, most of these earthquakes occurred in the northwestern segment. Only the M 7.5 Zheduotang earthquake in 1955 occurred in the Kangding section, and the seismicity in the southeastern segment was sparse. Although the Mw5.9 and Mw5.6 Kangding earthquakes occurred in 2014, their rupture scales were small, only about 30 km long (Yi et al., 2015). In addition, there have been no strong earthquakes in the Moxi section since 1786, and 234 years have passed since the last earthquake. Thus, the seismicity of this section is also worthy of attention.
3.1 Data collection and processing
The GPS data used in this study were collected from the Crustal Movement Observation Network of China (CMONOC I) and the Tectonic and Environmental Observation Network of Mainland China (CMONOC II). Most of the stations have operated for over a 10-year span or have had at least five observational campaigns with the campaign sites being observed for 4 days during each campaign, which allowed us to obtain high-accuracy velocities. In general, the crustal deformation recorded by the GPS stations contains information about the steady block motion and the post-seismic effects. It should be noted that three moderate-strong earthquakes, including the Mw7.9 Wenchuan earthquake in 2008, the Mw6.6 Lushan earthquake in 2013, and the Mw5.9 and Mw5.6 Kangding earthquake in 2014, occurred in the study area. We divided the GPS data into two periods, 1999–2007 and 2016–2020, corresponding to the periods before and after these three moderate-strong earthquakes to obtain the regional dynamic deformation characteristics. The first period mainly consists of the data from the campaign stations, which were collected in 1999, 2001, 2004, and 2007. The second period includes continuous and campaign stations and the campaign stations’ data were collected in 2016, 2017, 2018, 2019, and 2020.
We processed the GPS data using GAMIT/GLOBK (version 10.61) (Herring et al., 2015a, 2015b). The process mainly included three steps. (1) We processed the original observation data using GAMIT to obtain the loosely constrained daily solutions and satellite orbits. (2) Then, we performed the joint adjustment of the loosely constrained daily solutions of all the stations to obtain their position time series and estimations of the velocities in ITRF2014 using the software GLOCK10.61. For the details of the method used, refer to Liang et al. (2018, 2021) and Su and Zhan (2021). (3) We transformed the velocity solutions into a fixed-Eurasian reference frame (Fig. 3) using the Euler vector for Eurasia with respect to the International Terrestrial Reference System (Altamimi et al., 2017).
3.2 Methods and model
3.2.1 Strain rate calculations
The strain rate field is an important indicator for describing the regional deformation. It is not affected by the reference frame and can reflect the crustal deformation characteristics on different spatial scales (Ward, 1994; Reddy et al., 2000; Jin and Park, 2006; Zhang and Sagiya, 2017). Considering the observation errors of the actual observation data, we used the least-squares collocation method (Wu et al., 2009, 2011) to derive the strain rate fields, which has a relatively good robustness. In this method, the functional relationship between the GPS velocity field and the station locations is constructed, and then, the partial derivative function of the displacement and strain are used to directly obtain the strain rate tensors ελ, εφ, and ελφ of any grid point. The details of the calculation process have been described by Wu et al. (2011).
Then, using Eq. (2), additional strain rate vectors can be calculated.
where ε1 is the principal tensional strain rate, ε2 is the principal compressive strain rate, γmax is the maximum shear strain rate, and εsurface is the surface strain rate.
3.2.2 GPS profiles and fault slip rates
The GPS profile can be used to reveal the differential movements and the relationship between the displacement distribution and faults, and it is widely used in the analysis of fault deformation characteristics (Jiang et al., 2009; Maurin et al., 2010; Zhang, 2008). First, we selected the GPS stations within a certain range on both sides of the fault zone. Then, we directly projected the GPS velocities onto the directions perpendicular and parallel to the fault, and thus, the velocity component was obtained. In addition, we used an elastic dislocation model (Savage and Burford, 1973) to estimate the long-term slip rate.
3.2.3 DEFNODE negative dislocation inversion
We used the elastic block model (McCaffrey et al., 2002). This modeling was implemented using the open source DEFNODE program, and it has mostly been used to estimate fault slip deficit. The model combines block rotation, uniform block strain, and Okada dislocations (Eq. 1), and it can simultaneously estimate the Block’s Euler poles, the fault slip rate, and the fault-locking coefficient to provide the best fit of the geodetic, geologic, and seismic data. In this approach, the spatial variability of the fault coupling is defined by the parameter Phi (the coupling coefficient). If Phi = 0, the fault is fully creeping or decoupled for a long-term slip rate. If Phi = 1, the fault is completely locked. A value of Phi between 0 and 1 indicates that the fault is partially locked. The slip rate deficit on the fault plane was derived by multiplying the fault slip vector by Phi, which provides the elastic deformation around the fault.
(3)where V ̅sf is the observed surface velocity, V ̅ br is the velocity caused by the rotation of the block, V ̅ is is the velocity caused by the internal strain of the block, and V ̅ fs is the velocity caused by the negative dislocation of the fault locking.
Compared with general inversion programs of the same type, the DEFNODE program increases the inversion of the permanent strain inside the block, and the permanent strain inside the block can replace the partial postseismic viscoelastic relaxation effect. This program can be used for related inversions after the earthquake (Zhao et al., 2018).
3.2.4 Inversion model
The study area is divided into three tectonic blocks: the Bayan Har Block, the Sichuan-Yunnan Block, and the South China Block (Zhang et al., 2003). The Xianshuihe Fault Zone, the Longmenshan Fault Zone, and the Anninghe–Xiaojiang fault system were included in our inversion. The block properties (rigid body or elasticity) need to be set before the inversion. Based on previous studies (Zhao et al., 2015), we defined the Sichuan-Yunnan and Bayan Har Blocks as elastic during the inversion, which means that the locking degree (locking coefficient), the Euler vector of the block motion, and the uniform strain tensor inside the block are inverted simultaneously. The interior of the South China Block is relatively stable with weak deformation, and thus, it was modeled as rigid and the components of the strain rate were not calculated in the inversion.
Eleven nodes were set along the fault’s strike, and the average distance between these nodes was approximately 35–40 km. The depth of the fault model was set to 20 km based on the results of the following previous studies. Zhu et al. (2005) relocated the small earthquakes that occurred in western Sichuan and found that 90% of the small earthquake events in the Xianshuihe Fault Zone occurred in the upper crust at depths of 0–15 km. Liu et al. (2014) concluded that the stress concentration was generally within 20 km of the upper crust based on joint inversion of P-wave receiver functions and background noise. The details of the models structure are shown in Fig. 4.
4.1 Strain rate fields
Since the differences in the station spacing during 1999–2007 and 2016–2020 would have a certain effect on the strain rate fields, we selected the data for the same GPS station in both time periods and used the least-squares collocation method to calculate the dynamic strain rate fields, which covered the Xianshuihe Fault Zone. The calculation results are shown in Fig. 5.
The strain rate fields show that (Fig. 5a) before the Wenchuan earthquake, the maximum shear strain rate was relatively high in the Luhuo and Daofu sections in the northwestern segment of the Xianshuihe Fault Zone, with maximum values of 6.0×10− 8/a and 5.7×10− 8/a, respectively, which was significantly higher than that of the southeastern section (Table 2). However, after the Wenchuan and Lushan earthquakes, the maximum shear strain rates of the Luhuo and Daofu sections decreased significantly (Fig. 5b), while that of the southeastern segment (Kangding and Moxi sections) increased, that is, the maximum shear strain accumulation rate increased from 4.5×10− 8/a to 6.2×10− 8/a (Table 2), reflecting the shear deformation pattern of the Luhuo and Daofu section underwent a dynamic adjustment after the Wenchuan and Lushan earthquakes. In addition, the Qianning section also deserves special attention, because the shear strain accumulation rate of this section remained unchanged during the two periods, indicating that this section tends to be locked.
|
Number |
Fault section |
1999–2007 |
2016–2020 |
||
|---|---|---|---|---|---|
|
value (a− 1) |
error |
value (a− 1) |
error |
||
|
S1 |
Luhuo |
6.0×10− 8 |
0.34 |
5.0×10− 8 |
0.18 |
|
S2 |
Daofu |
5.7×10− 8 |
0.29 |
3.2×10− 8 |
0.20 |
|
S3 |
Qianning |
4.3×10− 8 |
0.34 |
3.2×10− 8 |
0.19 |
|
S4 |
Kangding |
4.5×10− 8 |
0.33 |
6.2×10− 8 |
0.21 |
|
S5 |
Moxi |
4.5×10− 8 |
0.39 |
6.2×10− 8 |
0.20 |
The principal strain rates calculated for 1999–2007 and 2016–2020 (Fig. 5) demonstrate that the roughly E–W directed compression and N–S directed extension were dominant in the northwestern segment of the Xianshuihe Fault, and from NW to SE along the strike of the Xianshuihe Fault, the compression gradually transitioned to NW–SE directed and the extension transitioned to NE–SW directed. Furthermore, we noted that the values of the principal strain rates in the northwestern segment were larger than those in the southeastern segment during 1999–2007. However, after the Wenchuan and Lushan earthquakes, the values of the principal strain rate in the northwestern segment decreased, while the values in the southeastern segment significantly increased, indicating that the tectonic stress on the southeastern section enhanced
4.2 GPS velocity profile
We created two velocity profiles (A and B) across the Xianshuihe Fault Zone. The locations of profiles are shown in Fig. 3. The profiles are 400 km long and 130 km (A) and 110 km (B) wide.
The velocity profile results in Fig. 6 show that an obvious left-lateral strike deformation characteristics dominates in the Xianshuihe Fault Zone. Before the Mw7.9 Wenchuan earthquake, the left-lateral strike-slip rate of the northwestern segment of the Xianshuihe Fault Zone was ~ 11.5 mm/a, but it decreased to ~ 7.8 mm/a after the Wenchuan and Lushan earthquakes. Moreover, the near-field deformation width of the fault was about 125 km before the Wenchuan earthquake, but it decreased to within the vicinity of the fault (~ 50 km) after the earthquake.
The velocity profile results (Figs. 6c, 6d) show that before the Mw7.9 Wenchuan earthquake, the left-lateral strike-slip rate of the southeastern segment of the Xianshuihe Fault Zone was ~ 11 mm/a, but it hardly decreased after the Mw6.6 Lushan and Mw5.9 Kangding earthquakes, kept a left-lateral strike-slip rate of ~ 10.8 mm/a. Moreover, the near-field deformation width of the fault almost unchanged during the two periods.
4.3 Fault locking
Based on the fault model above, we obtained the optimal model parameters through multiple calculations. Figure 7 shows the GPS velocity residuals for the two time periods for the best-fit model. The predicted velocities fit the GPS velocities well without systematic error. The majority of the residuals are within the error range, and the directions of the residuals are stochastic throughout the entire study area. Only the residuals near the Zemuhe Fault Zone are slightly larger. 7
The inversion results show that before the Wenchuan earthquake (Fig. 8a), the Luhuo section of the Xianshuihe Fault Zone was coupled, with a maximum fault-locking depth of ~ 16 km, and the average fault-coupling coefficient was greater than 0.85 at depths of 0–16 km. The Daofu section was shallowly coupled at depths of < 5 km and the Qianning section was coupled only at the surface. Moreover, the southeastern segment of the Xianshuihe Fault Zone exhibited a high degree of locking, especially in the Kangding section, which was locked to a depth of ~ 16 km. The locking depth in the southeastern segment gradually weakened from northwest to southeast. However, after the Wenchuan and Lushan earthquakes (Fig. 8b), the locking degree of the Kangding section greatly decreased, and it became almost completely decoupled and freely creeped. In contrast, the locking degrees of the Daofu, Qianning, and Moxi sections increased, with average coupling coefficients of greater than 0.9 at depths 8 km, 8 km and 16 km, respectively. The details are shown in Fig. 8b.
The fault-parallel slip-rate deficit demonstrates that the Xianshuihe Fault Zone experienced a left-lateral slip-rate deficit (positive value) (Fig. 9). Before the Wenchuan earthquake, the left-lateral slip-rate deficit of the Luhuo and Daofu sections was the highest (8–10 mm/a), followed by the Moxi section (6–8 mm/a). After the Wenchuan and Lushan earthquakes, the left-lateral slip-rate deficit decreased to 2–6 mm/a in the Luhuo, Daofu, and Qianning sections. After the Wenchuan and Lushan earthquakes, the Kangding section exhibited almost no left-lateral slip-rate deficit, and that of the northern tip of the Moxi section remained at approximately 6–8 mm/a.
The spatial pattern and rate of the strain accumulation during the interseismic phase are very important for interpreting the mechanisms of earthquakes and for evaluating the seismic potentials of faults (Navarro et al., 2003; Fialko, 2006; Moreno et al., 2010).
5.1 Different impacts of major earthquakes on the Xianshuihe Fault Zone
The current deformation characteristics of the Xianshuihe Fault Zone can be described by the strain rate field, the fault slip across the block boundaries, and the fault coupling based on geodetic data (Zhao et al., 2018; Li et al., 2019; Wang and Shen, 2020). Our results show that the high shear strain rates were concentrated in the Xianshuihe Fault Zone, and the values in the northwestern section were significantly larger than those in the southeastern section before the Wenchuan earthquake, indicating that the strain accumulation rate in the northern section is larger than that in the southern section, which is consistent with the high slip rate in the northwestern section and the relative low slip rate in the southeastern section (Meade, 2007; Qu et al., 2018; Zhang et al., 2019), exhibiting the basic features of the block-like tectonic extrusion model (Avouac and Tapponnier, 1993; Tapponnier et al., 1982).
After the strong Wenchuan and Lushan earthquakes, the shear strain rate in the Luhuo and Daofu sections significantly decreased. The left-lateral slip-rate and the slip-rate deficit in these sections also decreased (Figs. 6b, 9b). We conclude that these variations were mainly caused by the large-scale decoupling of the central and southern sections of the Longmenshan Fault Zone accompanied by the Wenchuan and Lushan earthquakes. After the Wenchuan earthquake, the extrusion rate of the Bayan Har Block toward the Sichuan Basin accelerated which decreased the far-field left-lateral strike-slip rate and the fault slip-rate deficit and increased the tensile movement of the northwestern section of the Xianshuihe Fault. These major earthquakes may have alleviated the strain accumulation rate in the northwestern segment of the Xianshuihe Fault Zone.
To further verify the impacts of the Wenchuan earthquake and the other regional strong earthquakes on the Xianshuihe Fault Zone; we also collected near-field cross-fault short-baseline observational data for the northwestern segment of the Xianshuihe Fault Zone from 1982 to 2018. The distributions of the fault-crossing sites are shown in Fig. 3. Fault-crossing deformation observations can directly reflect the local tectonic deformation and changes in the stress state of active faults, and they have been used to analyze near-field fault activity (Bo et al., 1998; Li et al., 2016; Li et al., 2017). In this study, the calculation method of Bo et al. (1998) was used to quantitatively calculate the horizontal strike-slip component and the horizontal tension and compression components. The positive slope of the curve indicates right-lateral and tension movements (Fig. 10).
The results of the fault activity analysis (Fig. 10) reveal that the northwestern segment of the Xianshuihe Fault Zone has predominantly experienced left-lateral strike-slip movement (the slope of the curve was negative), with a low tension and compression rate, before the Wenchuan earthquake. However, after the Wenchuan earthquake, and especially after the Lushan earthquake, the tensile activity in the northwestern segment of the Xianshuihe Fault Zone increased significantly. Furthermore, it was found that the near-field left-lateral strike-slip movement at some of the sites, such as Zhuwo, Goupu, and Laoqianning, decreased slightly. Therefore, the results of the near-field cross-fault observations also indicate that the Wenchuan and Lushan earthquakes affected the movement features and strain accumulation rate of the northwestern section of the Xianshuihe Fault Zone.
In addition, we noted that the shear strain rates of the southeastern section of the Xianshuihe Fault increased after the series of earthquakes (Fig. 5b). The Kangding section has almost no left-lateral slip-rate deficit (Fig. 9b) and still has a high left-lateral slip-rate, while the Moxi section remains strongly coupled, with a high left-lateral slip-rate deficit of 8 mm/a in the second period (Figs. 8b, 9b). This is consistent with the latest InSAR results (Qiao et al., 2021), which also indicate a high shallow slip rate (16.3–19.8 mm/a) in the western Kangding section during 2014–2019, and a locking depth of 18.5 km in the Moxi section. The southeastern section of the Xianshuihe Fault Zone is located in the triple junction area where the three major fault zones interact and restrict one another. The Wenchuan and Lushan earthquakes occurred in the central northern segment and southern segment of the Longmenshan Fault Zone, respectively. The post-seismic geological survey and precise aftershock positioning results reveal that the ruptures of the Lushan earthquake and the Wenchuan earthquake did not pass through, and the Lushan earthquake did not cause the entire southern section of the Longmenshan fault to decouple (Fang et al., 2013; Gao et al., 2014; Xu et al., 2013). The triple junction area is still in a strong compressional deformation environment. More than a year after the Lushan earthquake, the Mw5.9 and Mw5.6 Kangding earthquakes on November 22, 2014, occurred on the Selaha fault between the Kangding and Qianning sections of the Xianshuihe Fault Zone, which also confirms that this area had a high level of stress accumulation after the Wenchuan and Lushan earthquakes. However, owing to the small-scale rupture and the limited energy released by the Mw5.9 Kangding earthquake, the event was not sufficient to significantly change the strain throughout the entire triple junction area. Therefore, the southeastern section of the Xianshuihe Fault still had a high strain accumulation rate after the series of earthquakes. All of these results indicate that the surrounding strong earthquakes, i.e., the Mw7.9 Wenchuan earthquake in 2008 and the Mw6.6 Lushan earthquake in 2013, had evident differential impacts on the strain accumulation pattern of the Xianshuihe Fault Zone.
5.2 Seismic hazard analysis
The seismic gaps, locking degree, low b-value section of the background seismic activity, and Coulomb stress changes are effective indicators for determining the locations where strong earthquakes may occur in the future (Hubert-Ferrari et al., 2000; Nalbant et al., 1998; Shan et al., 2015; Stein et al., 1997; Wen et al., 2008; Yi et al., 2008).
The spatial and temporal distributions of the earthquake ruptures caused by the major earthquakes along the Xianshuihe Fault Zone reveal that no strong earthquakes have occurred in the Qianning section since the M7.0 Qianning earthquake in 1893, and the elapsed time exceeds the recurrence interval of the historic strong earthquakes, indicating that this seismic gap is a potential location for strong earthquakes in the future. Following the Wenchuan and Lushan earthquakes, the locking degree and locking depth of the Qianning section increased. In addition, the GPS-derived strain rates show that the Qianning section is an invariant deformation adjustment zone. These results indicate that the Qianning section is in a coupled state with little deformation. Yi et al. (2008) calculated the b-value from seismic data (from the China Seismic Network and Sichuan Regional Seismic Network) for 1976–2006 and found a large-scale asperity with a low b-value from Daofu to Qianning (defined as the Qianning section here) of the Xianshuihe Fault Zone. The weak deformation of this section may be related to the large-scale asperity. In addition, the Coulomb stress change shows that the occurrences of the Wenchuan and Lushan earthquakes caused increased stress in the Qianning section (Li et al., 2012; Luo and Liu, 2018; Parsons et al., 2008; Shan et al., 2013; Toda et al., 2008; Wan and Shen, 2010),which is consistent with our inversion results. Therefore, we suggest that the Qianning section may be a potential seismogenic source for a future strong earthquake.
Another potential seismogenic source of a strong future earthquake is the Moxi section, which is located at the triple junction of the Xianshuihe Fault Zone, the Anninghe Fault Zone, and the Longmenshan Fault Zone. The tectonics of this area are especially complex. It has been over 200 years since the M7.75 Kangding-Moxi earthquake occurred in 1786 in the Moxi section (Wen et al., 2008). Owing to the influence of long-term tectonic loading, this section already has a high strain accumulation (Wang et al., 2016). In addition, our fault-coupling inversion results show that this section has a high degree of fault-locking and a high left-lateral rate deficit. Moreover, the Coulomb stress change in this section has been enhanced by historical strong earthquakes (Li et al., 2020; Shao et al., 2016). Therefore, the risk of future strong earthquakes in this area also requires more attention. In contrast, the locking degree of the Kangding section has almost no slip-rate deficit. This may be related to moderate Mw5.9 and Mw5.6 Kangding earthquakes in 2014, which caused part of the energy in the Kangding section to be released. High-resolution InSAR velocity maps have been used to identify postseismic afterslip near the 2014 Kangding earthquake rupture (Qiao and Zhou, 2021; Li and Bürgmann, 2021). However, owing to a lack of high precision, near-field, intensive geodetic observations of the study area, we cannot clearly determine the deformation mechanism.
In this study, we focus on exploiting the dynamic deformation characteristics of the Xianshuihe Fault Zone and its potential risk of strong earthquakes occurring using GPS data for 1999–2007 and 2016–2020. The main conclusions drawn are as follows.
1) The high shear strain rate was mainly concentrated in the Xianshuihe Fault Zone, and the values of the northwestern section were greater than those of the southeastern section before the Wenchuan earthquake, which is consistent with the spatial differences in the slip rate.
2) The deformation mode of the Xianshuihe Fault Zone underwent a dynamic adjustment after the Wenchuan and Lushan earthquakes. The shear strain accumulation rates of the Luhuo and Daofu sections decreased from 6.0×10− 8/a to 5.0×10− 8/a and from 5.7×10− 8/a to 3.2×10− 8/a, respectively, while that of the southeastern segment increased from 4.5×10− 8/a to 6.2×10− 8/a. The slip rate and deformation width of the Xianshuihe Fault Zone also changed during these two periods, reflecting the differential impacts of the surrounding major earthquakes on the deformation pattern of the Xianshuihe Fault Zone.
3) The locking degree of the Qianning and Moxi sections increased significantly after the Wenchuan earthquakes, while no slip-rate deficit occurred in the Kangding section, which is probably because part of the energy was released by the Mw5.9 and Mw5.6 Kangding earthquakes in 2014. Strain accumulation characteristics and recent rupture history of the Xianshuihe Fault Zone suggest that both the Qianning and Moxi sections have a high potential for future strong earthquakes.
GPS: Global positioning system; InSAR: Interferometric Synthetic Aperture Radar; CMONOC I: Crustal Movement Observation Network of China; CMONOC II: Tectonic and Environmental Observation Network of Mainland China.
Availability of data and materials
The horizontal GPS velocity data used in this study can be accessed at ftp://60.30.77.19.
Competing Interest
This manuscript has not been published or presented elsewhere in part or in entirety and is not under consideration by another journal. We have read and understood your journal’s policies, and we believe that neither the manuscript nor the study violates any of them. There are no conflicts of interest to declare.
Funding
This study was supported by the National Natural Science Foundation of China (Nos. 41974011, 41874116), the Science for Earthquake Resilience Program (No. XH18067Y), and the National Key Research and Development Program of China (No. 2018YFC1503606).
Author Contributions
Layue Li, Yanqiang Wu, and Yujiang Li designed the research, analyzed the results and wrote the manuscript. Wei Zhan and Xinzhong Liu performed the GPS data processing and participated in the creation of the graphics. All authors have read and agreed to the published version of the manuscript.
Acknowledgements
We acknowledge the CMONOC I (Crustal Movement Observation Network of China) and CMONOC II (Tectonic and Environmental Observation Network of Mainland China) for providing GPS observations. We would like to thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript. Most figures were plotted using the Generic Mapping Tools (Wessel et al., 2013).
- Ader T, Avouac JP, Liu-Zeng J, Lyon-Caen H, Bollinger L, Galetzka J, Genrich J, Thomas M, Chanard K, Sapkota SN (2012) Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard. J Geophys Res 117(B4):B04403
- Allen CR, Luo ZL, Qian H, Wen XZ, Zhou HW, Huang WS (1991) Field study of a highly active fault zone: The Xianshuihe fault of southwestern China. Geol Soc Am Bull 103(9):1178–1199
- Altamimi Z, Métivier L, Rebischung P, Rouby H, Collilieux X (2017) ITRF2014 plate motion model. Geophys J Int 209:1906–1912. doi:10.1093/gji/ggx136
- Avouac JP, Tapponnier P (1993) Kinematic model of active deformation in central Asia. Geophys Res Lett 20(10):895–898
- Bai M, Chevalier ML, Pan J, Replumaz A, Leloup PH, Métois M, Li H (2018) Southeastward increase of the late Quaternary slip–rate of the Xianshuihe fault, eastern Tibet. Geodynamic and seismic hazard implications. Earth Planet Sci Lett 485:19–31
- Bo W, Xie J, Guo L (1998) Analysis and exploration of deformation of Babaoshan fault belt. Earthquake 18(1):63–68
- Deng QD, Cheng SP, Ma J, Du P (2015) Seismic activities and earthquake potential in the Tibetan Plateau. Chin J Geophys 57(5):678–697
- Dogru A, Gorgun E, Aktug B, Ozener H (2018) Seismic hazard assessment of the central North Anatolian Fault (Turkey) from GPS–derived strain rates and b–values. Geomatics Natural Hazards Risk 9:356–367
- Fang LH, Wu JP, Wang WL, Lv ZY, Wang CZ, Yang T, Cai Y (2013) Relocation of mainshock and aftershock sequences of Ms7.0 Sichuan Lushan earthquake. Chin Sci Bull 58:1901–1909
- Fialko Y (2006) Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system. Nature 441:968–971
- Gan W, Zhang P, Shen ZK, Niu Z, Wang M, Wan Y, Zhou D, Cheng J (2007) Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements. J. Geophys. Res. 112
- Gao Y, Wang Q, Zhao B, Shi Y (2014) A rupture blank zone in middle south part of Longmenshan Faults: Effect after Lushan Ms 7.0 earthquake of 20 April 2013 in Sichuan, China. SciChina (Earth Sci) 57:2036–2044
- Guo RM, Zheng Y, Tian W, Xu JQ, Zhang WT (2018) Locking status and earthquake potential hazard along the Middle–South Xianshuihe Fault. Remote Sensing 10(12):2048
- Herring TA, King RW, McClusky SC (2010a) GAMIT Reference Manual. GPS Analysis at MIT. Release 10.4. Massachusetts Institute Technology, http://www–gpsg.mit.edu/~simon/gtgk/index.htm
- Herring TA, King RW, McClusky SC (2010b) GLOBK Reference Manual. Global Kalman filter VLBI and GPS analysis program. Release 10.4. Massachusetts Institute Technology, http://www–gpsg.mit.edu/~simon/gtgk/index.htm
- Hubert–Ferrari A, Barka A, Jacques E, Nalbant SS, Meyer B, Armijo R, Tapponnier P, King GC (2000) Seismic hazard in the Marmara Sea region following the 17 August 1999 Izmit earthquake. Nature 404:269–273
- Ji L, Zhang W, Liu C, Zhu L, Xu J, Xu X (2020) Characterizing interseismic deformation of the Xianshuihe fault, eastern Tibetan Plateau, using Sentinel–1 SAR images. Adv Space Res 66:378–394
- Jiang G, Xu X, Chen G, Liu Y, Fukahata Y, Wang H, Yu G, Tan X, Xu C (2015) Geodetic imaging of potential seismogenic asperities on the Xianshuihe–Anninghe–Zemuhe fault system, southwest China, with a new 3–D viscoelastic interseismic coupling model. J Geophys Res 120(3):1855–1873
- Jiang ZS, Fang Y, Wu YQ, Wang M, Du F, Ping JJ (2009) The dynamic process of regional crustal movement and deformation before Wenchuan Ms8.0 earthquake. Chin J Geophys 52(2):505–518
- Jin S, Park PH (2006) Strain accumulation in South Korea inferred from GPS measurements. Earth Planets Space 58(5):529–534
- King GC, Stein RS, Lin J (1994) Static stress changes and the triggering of earthquakes. Bull Seismol Soc Am 84:935–953
- Li J, Zhou B, Li T, Yang Y, Li Z (2020) Locking depth, slip rate, and seismicity distribution along the Daofu–Kangding segment of the Xianshuihe fault system. eastern Tibetan Plateau JAsian Earth Sci 193:104328
- Li LY, Ji LY, Li YJ, Wei Z (2016) Current activity characteristics of Xianshuihe fault zone and its relationship with strong earthquakes based on fault–crossing measurements. Acta Seism Sin 38(5):739–750
- Li Y, Chen L, Liu S, Zhan W (2017) Discussion on mechanism of cross–fault baseline variations along the northwestern segment of the Xianshuihe fault zone. Chin J Geophys 60(2):554–563
- Li Y, Chen L, Lu Y (2012) Numerical simulation on the influences of Wenchuan earthquake on the surrounding faults. Earthquake Sci 25(2):143–150
- Li Y, Liu M, Li Y, Chen L (2019) Active crustal deformation in southeastern Tibetan Plateau: The kinematics and dynamics. Earth Planet Sci Lett 523:115708
- Li YJ, Shi FQ, Zhang H, Wei WX, Xu J, Shao ZG (2020) Coulomb stress change on active faults in Sichuan–Yunnan region and its implications for seismic hazard. Seismol Geol 42(2):526–546
- Li YX, Bürgmann R (2021) Partial coupling and earthquake potential along the Xianshuihe Fault, China. JGeophysRes:Solid earth. doi:10.1029/2020JB021406
- Liang H, Wu Y, Chen C, Xiong W (2018) Coseismic displacement and dislocation inversion of 2014 Yutian M(s) 7.3 earthquake in Xinjiang. Chin J Geophys 61(12):4817–4826
- Liang HB, Zhan W, Li JW (2021) Vertical surface displacement of mainland China from GPS using the multisurface function method. Adv Space Res. https://doi.org/10.1016/j.asr.2021.02.024
- Liu QY, Van Der Hilst RD, Li Y, Yao HJ, Chen JH, Guo B, Qi SH, Wang J, Huang H, Li SC (2014) Eastward expansion of the Tibetan Plateau by crustal flow and strain partitioning across faults. Nat Geosci 7(5):361–365
- Luo G, Liu M (2018) Stressing Rates and Seismicity on the Major Faults in Eastern Tibetan Plateau. J Geophys Res 123:10968–10986
- Maurin T, Masson F, Rangin C, Min UT, Collard P (2010) First global positioning system results in northern Myanmar: Constant and localized slip rate along the Sagaing fault. Geology 38(7):591–594
- McCaffrey R, Stein S, Freymueller J (2002) Crustal block rotations and plate coupling. Plate Boundary Zones Geodyn Ser 30:101–122
- Meade BJ (2007) Present–day kinematics at the India–Asia collision zone. Geology 35:81–84
- Miura S, Sato T, Tachibana K, Satak Y (2002) Strain accumulation in and around ou backbone range, northeastern japan as observed by a dense gps network. Earth Planets Space 54(11):1071–1076
- Moreno M, Rosenau M, Oncken O (2010) 2010 Maule earthquake slip correlates with pre–seismic locking of Andean subduction zone. Nature 467:198–202
- Nalbant SS, Hubert A, King GC (1998) Stress coupling between earthquakes in northwest Turkey and the north Aegean Sea. J Geophys Res 103:24469–24486
- Navarro A, Ca Talo J, Fernandes R, Miranda M, Bastos L (2003) The Terceira island (Azores) crustal deformations from GPS data. Earth Planets Space 55:637–642
- Parsons T, Ji C, Kirby E (2008) Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin. Nature 454:509–510
- Qian H, Allen C, Luo Z, Wen X, Zhou H, Huang W (1988) The active characteristics of Xianshuihe fault in Holocene. Earthq Res China 4(2):9–18
- Qiao X, Zhou Yu (2021) Geodetic imaging of shallow creep along the Xianshuihe fault and its frictional properties. Earth Planet Sci Lett 567:117001
- Qu W, Lu Z, Zhang Q, Hao M, Wang Q, Qu F, Zhu W (2018) Present–day crustal deformation characteristics of the southeastern Tibetan Plateau and surrounding areas by using GPS analysis. J Asian Earth Sci 163:22–31
- Reddy CD, El-Fiky Gamal, Teruyuki K, Seiichi Shimada KV, Kumar (2000) Crustal strain field in the deccan trap region, western india, derived from gps measurements. Earth Planets Space 52:965–969
- Riguzzi F, Crespi M, Devoti R, Doglioni C, Pietrantonio G, Pisani AR (2012) Geodetic strain rate and earthquake size: new clues for seismic hazard studies. Phys Earth Planet Inter 206–207,67–75, http://dx.doi.org/10.1016/ j.pepi.2012.07.005
- Savage JC, Burford RO (1973) Geodetic determination of relative plate motion in central California. J geophys res 78(5):832–845
- Savage JC, Svarc JL, Prescott WH (1999) Geodetic estimates of fault slip rates in the San Francisco Bay area. J Geophys Res: Solid Earth 104:4995–5002
- Segall P (2010) Earthquake and volcano deformation. Princeton University Press, pp 1–432
- Stein RS (2003) Earthquake conversations. Sci Am 288:72–79
- Shan B, Xiong X, Wang R, Zheng Y, Yadav R (2015) Stress evolution and seismic hazard on the Maqin–Maqu segment of East Kunlun Fault zone from co-, post-and interseismic stress changes. Geophys J Int 200:244–253
- Shan B, Xiong X, Wang R, Zheng Y, Yang S (2013) Coulomb stress evolution along Xianshuihe–Xiaojiang Fault System since 1713 and its interaction with Wenchuan earthquake, May 12, 2008. Earth Planet. Sci. Lett. 377–378, 199–210
- Shao Z, Xu J, Ma H, Zhang L (2016) Coulomb stress evolution over the past 200years and seismic hazard along the Xianshuihe fault zone of Sichuan, China. Tectonophysics 670:48–65
- Stein RS, Barka AA, Dieterich JH (1997) Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering. Geophys J Int 128:594–604
- Su G, Zhan W (2021) Seasonal and long-term vertical land motion in Southwest China determined using GPS, GRACE, and surface loading model. Earth Planets Space73,131. https://doi.org/10.1186/s40623-021-01459-4
- Tapponnier P, Peltzer G, Dain A, Armijo R, Cobbold P (1982) Propagating extrusion tectonics in Asia: New insights from simple experiments with plasticine. Geology 10(12):611–616
- Toda S, Lin J, Meghraoui M, Stein RS (2008) 12 May 2008 M = 7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems. Geophys Res Lett 35(17). doi:10.1029/2008GL034903
- Wan Y, Shen ZK (2010) Static Coulomb stress changes on faults caused by the 2008 Mw 7.9 Wenchuan, China earthquake. Tectonophysics 491:105–118
- Wang E, Burchfiel BC, Royden LH, Chen L, Chen J, Li W, Chen Z (1998) Late Cenozoic Xianshuihe–Xiaojiang, Red River, and Dali fault systems of southwestern Sichuan and central Yunnan, China. Geol Soc Am Spec Pap 327:1–108
- Wang H, Wright TJ, Biggs J (2009) Interseismic slip rate of the northwestern Xianshuihe fault from InSAR data. Geophys Res Lett 36(3):139–145
- Wang M, Shen ZK (2020) Present–Day crustal deformation of continental China derived from GPS and its tectonic implications. J. Geophys. Res. 125, e2019JB018774
- Wang X, Xiao J, Xu H, He J (2016) Dynamic responses of the Xianshuihe and Longmenshan fault zones to regional tectonic loading. Chin J Geophys 59:1403–1413
- Ward SN (1994) A multidisciplinary approach to seismic hazard in southern California. Bull Seismol Soc Am 84(5):1293–1309
- Wen X, Allen C, Luo Z, Qian H, Zhou H, Huang W (1989) Segmentation, geometric features, and their seismotectonic implications for the Holocene Xianshuihe fault zone. Acta Seism Sin 11:362–372
- Wen XZ, Ma SL, Xu XW, He YN (2008) Historical pattern and behavior of earthquake ruptures along the eastern boundary of the Sichuan–Yunnan faulted–block, southwestern China. Phys Earth Planet Inter 168(1–2):16–36
- Weertman J, Weertman JR (1966) Elementary Dislocation Theory. Oxford University Press, Oxford
- Wessel P, Smith WH, Scharroo R, Luis J, Wobbe F (2013) Generic mapping tools: improved version released. Eos Trans Am Geophys Union 94:409–410
- Wu YQ, Jiang ZS, Yang GH, Fang Y, Wang WX (2009) The application and method of GPS strain calculation in whole mode using least–squares collocation in sphere surface. Chin J Geophys 52(7):1707–1711
- Wu YQ, Jiang ZS, Yang GH, Wei W, Liu X (2011) Comparison of GPS strain rate computing methods and their reliability. GeophysJInt 185(2):703–717
- Wu YQ, Jiang ZS, Liang HB, Pang YJ, Zhu S, Liu C, Chen CY, Li JW (2018) Deformation response of seismogenic faults to the wenchuan Ms 8.0 earthquake: a case study for the southern segment of the longmenshan fault zone. Remote Sensing 10(6):894
- Xu XW, Wen XZ, Han ZJ, Chen GH, Li CY, Zheng WJ, Zhang SM, Ren ZQ, Xu C, Tan XB (2013) Lushan Ms7.0 earthquake: A blind reserve–fault event. Chin Sci Bull 58:3437–3443
- Xu XW, Wen XZ, Zheng RZ, Ma WT, Song FM, Yu GH (2003) Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan–Yunnan region, China. SciChina Ser D 46:210–226
- Yi G, Long F, Wen X, Liang M, Wang S (2015) Seismogenic structure of the M6.3 Kangding earthquake sequence on 22 Nov. 2014, Southwestern China. Chin J Geophys 58(4):1205–1219
- Yi GX, Wen XZ, Su YJ (2008) Study on the potential strong-earthquake risk for the eastern boundary of the Sichuan–Yunnan active faulted-block, China. Chin J Geophys 51(6):719–1725
- Zhang L, Cao D, Zhang J, Sui L (2019) Interseismic Fault Movement of Xianshuihe Fault Zone Based on Across–Fault Deformation Data and InSAR. Pure Appl Geophys 176:649–667
- Zhang P (2008) Present tectonic deformation, strain distribution and deep dynamic process in west of Sichuan, eastern margin of the Tibetan Plateau. Sci China Ser D 38:1041–1056
- Zhang P, Deng Q, Zhang G, Ma J, Gan W, Min W, Mao F, Wang Q (2003) Active tectonic blocks and strong earthquakes in the continent of China. Sci China Ser D Earth Sci 46:13–24
- Zhang X, Sagiya T (2017) Shear strain concentration mechanism in the lower crust below an intraplate strike-slip fault based on rheological laws of rocks. Earth Planets Space 69(1):82. DOI 10.1186/s40623-017-0668-5
- Zhao J, Jiang ZS, Niu AF, Liu J, Wu YQ, Wei WX, Liu XX, Yan W (2015) Study on dynamic characteristics of fault locking and fault slip deficit in the eastern boundary of the Sichuan–Yunnan rhombic block. Chin J Geophys 58(3):872–885
- Zhao J, Ren J, Jiang Z, Yue C (2018) Three–dimensional deformation characteristics of the Xianshuihe Fault Zone. Seismol Geol 40(4):818–831
- Zhao J, Ren J, Liu J, Liu X, Liang H, Niu A, Yue C, Yuan Z (2020) Coupling fraction and relocking process of the Longmenshan Fault Zone following the 2008 Mw7.9 Wenchuan earthquake–ScienceDirect. J Geodyn, 137. doi:https://doi.org/10.1016/j.jog.2020.101730
- Zou ZY, Jiang ZS, Wu YQ, Wei WX, Liu XX (2015) Dynamic characteristics of crustal movement in north–south seismic belt from GPS velocity field before and after the Wenchuan Earthquake. Chin J Geophys 58(5):1597–1609
- Zhu AL, Xu XW, Zhou YS, Yin JY, Gan WJ, Chen GH (2005) Relocation of small earthquakes in Western Sichuan, China and its implications for active tectonics. Chin J Geophys 48(3):692–700
