Seismogenic Faulting of the 2016 Mw 6.0 Hutubi Earthquake in the Northern Tien Shan Region: Constraints From Near-Field Borehole Strain Step Observations and Numerical Simulations

Abstract

Constraining the seismotectonic faulting that occurred as a result of the 2016 Mw 6.0 Hutubi earthquake provides valuable information about the deformation in the northern Tien Shan region. However, due to the lack of surface rupturing and high quality near-field teleseismic data, the exact nature of the faulting remains controversial. In our study, we analyze the coseismic strain time series of the Mw 6.0 Hutubi earthquake using strain data collected from nearby borehole stations. The tensile and compressive coseismic strain steps account for most of the recorded borehole data in this region. Employing a numerical model that is based on elastic dislocation theory, we reproduce the observed tensile and compressive coseismic strain steps using source parameters that were generated through seismic wave inversion, seismic reflection data, and aftershock relocation. By conducting a comparative analysis between the predicted and observed coseismic strain steps, we study the seismogenic faulting of the Mw 6.0 Hutubi earthquake. The results indicate that when the source parameters are 292°/62°/80° (strike/dip/rake), the predicted tensile and compressive characteristics for 13/16 channels are consistent with the observational data. Based on these results, we infer that the seismogenic faulting, which is located near the Horgos-Manas-Tugulu fault, can be characterized as a high-angle blind back-thrust fault with a north-dipping fault plane. Providing constraints on the seismogenic faulting associated with the 2016 Mw 6.0 Hutubi earthquake also yields to understand the mechanism of the overall deformation pattern in the northern Tien Shan region.

Introduction

On December 8, 2016, the Mw 6.0 Hutubi earthquake occurred within the active fold-and-thrust belts of the northern Tien Shan foreland basin, northwestern China (Li et al., 2018; Lu et al., 2018). However, as the absence of coseismic surface rupturing and high-quality far-field seismic observation record, the seismogenic faulting of this earthquake remains controversial.

In previous work, this earthquake fault has been described as either a low-angle south-dipping thrust fault or a high-angle north-dipping back-thrust fault (Liu et al., 2017; Liu et al., 2018; Lu et al., 2018; Gong et al., 2019; Wang et al., 2019; Xu et al., 2019; Yang et al., 2019). From seismic reflection data and seismic wave inversion results, it is postulated that the mainshock is controlled by a large-scale south-dipping thrust fault, while the aftershocks relocation results indicate the presence of a secondary back-thrust fault (Lu et al., 2018; Yang et al., 2018). However, the aftershock relocation suggests that this seismogenic faulting is a north-dipping secondary fault (Liu et al., 2018; Xu et al., 2019).

Earthquakes transfer stress to the surrounding rock, which can cause observable changes in geodetic data (Nykolaishen et al., 2015). Borehole strain observations record robust coseismic deformation information that not only provide constraints for the tensile and compressive coseismic stress and strains (Gladwin et al., 1994; Linde et al., 1996; Qiu and Shi, 2003; Johnston et al., 2006; Qiu et al., 2013; Barbour et al., 2015; Gong et al., 2019), but they can also act as parameters in numerical simulations of fault geometry. Using the observed four-component borehole strain data at the Changping station, Qiu and Shi (2003) analyzed the coseismic strain steps and estimated the corresponding stress changes induced by two Zhangbei earthquakes. By quantifying spatial differences in the borehole strain data that captured changes created by the 2004 Kii peninsular earthquakes, Asai et al. (2005) determined the coseismic strain steps in the Nankai Trough. Barbour et al. (2015) employed the data recorded on nine borehole strainmeters to characterize the coseismic strain caused by earthquakes with a range of magnitudes and distances. Based on the static-dynamic strain response of the area affected by the Mw 6.0 Hutubi earthquake, Gong et al. (2019) identified four stations where the borehole data capture resolvable static strain changes.

In order to analyze and characterize the coseismic strain step that was created by the Mw 6.0 Hutubi earthquake, we collected data from four near-field four-component borehole strainmeters that reside in observational stations located throughout the northern Tien Shan region. Based on elastic dislocation theory (Okada, 1985), we calculated the coseismic stress field using the different sets of source parameters that were inferred via seismic wave inversion, seismic reflection data, and aftershock relocation. After using these source parameters to generate synthetic strain step data, we then constrained the seismogenic faulting of the Mw 6.0 Hutubi earthquake by identifying the source parameters that minimized the difference between the predicted and observed tensile and compressive coseismic strain steps.

Tectonic Setting

The Tien Shan orogenic belt is bounded by the Tarim Basin to the south and the Junggar Basin to the north (Figure 1) (Zhang et al., 2003). As one of the largest orogenic belts within the Asian continent, this region was formed by multiple instances of plate collision and subduction (Xiao et al., 2004; Charvet et al., 2011; Chen et al., 2014; Qiu et al., 2019). Since the Cenozoic, the Tien Shan tectonic zone has experienced strong shortening and uplift due to continued convergence between the Indian and Eurasian plates (Molnar and Tapponnier, 1975; Allen et al., 1999; Qiu et al., 2019; Zhou et al., 2020); this ongoing collision has both created new faults and reactivated old faults (Yin et al., 1998). GPS data indicates that the crustal shortening rate in the Tien Shan tectonic zone increases from ∼5–10 mm/yr in the east to ∼20 mm/yr in the west (Yang et al., 2008; Liu et al., 2016). Deep seismic reflection profiles have revealed the thin-skinned tectonics of the crust beneath the northern Tien Shan piedmont, which consists of (from south to north) the piedmont Qigu fold-and-thrust zone, the Horgos-Manas-Tugulu active fold-and-thrust zone, and the Dushanzi-Anjihai active fold-and-thrust zone (Zhang et al., 1994; Yang et al., 2013; Gong et al., 2019). These active reverse fault-fold zones are the primary causes of the frequent seismicity in this region (Avouac et al., 1993; Zhang et al., 1996; Deng et al., 2000; Wang et al., 2004).

FIGURE 1

In 2016, the Mw 6.0 Hutubi earthquake occurred in the active fold-and-thrust belt of the northern Tien Shan foreland basin. This region has accommodated ∼50% of the predicted crustal shortening (Yin et al., 1998; Zubovich et al., 2010; Yang et al., 2019). Seismic reflection profiles indicated that the seismogenic faulting created by the Mw 6.0 Hutubi earthquake is closely related to the south-dipping Horgos-Manas-Hutubi fault (Lu et al., 2018). However, an aftershock relocation analysis provides support for the hypothesis that this earthquake occurred on a north-dipping back-thrust fault that exhibits the ‘y-type’ distribution with the main rupture plane of the 1906 Manas earthquake (Lu et al., 2018). If this is true, then this fold-type earthquake was caused by the continuous stress accumulation that arose during the ongoing collision between the Tien Shan region and the Junggar Basin (Liu et al., 2018).

Data

It has been repeatedly demonstrated that four-component borehole strainmeters yield high-quality and high-precision (10^–10) recordings of the static and dynamic strains caused by earthquakes (Qiu and Shi, 2003; Ouyang et al., 2004; Qiu et al., 2013; Gong et al., 2019). To better monitor the ongoing crustal deformation in this area and to explore the relationship between borehole coseismic strain steps and strong earthquakes, we deployed eleven RZB-type four-component borehole strainmeters in the northern Tien Shan region since 2015.

In this paper, we analyzed the tensile and compressive coseismic strain steps recorded at stations QEG, SC, BLT, and YSG. All four strainmeters were installed directly into the bedrock of the southern Junggar thrust fault, the Borokonu strike-slip fault, and thrust fault located along the northern margin of Bogda Mountain at depths of 35–157 m (Table 1).

The coseismic strain steps and subsequent changes recorded at stations QEG, SC, BLT, and YSG are shown in Figure 2. According to the definition of a tensional or compressional coseismic strain step (Qiu and Shi, 2003), the coseismic strain step is tensional when the step change △ε ＞ 0, as shown in the S1 component of the QEG observation station (Figure 2). The maximum tensile and compressive changes of the coseismic strain steps observed at stations QEG, SC, BLT, and YSG reach values of 1.6 × 10^–8, −3.4 × 10^–8, 2.5 × 10^–8, and −3.6 × 10^–9, respectively (Figure 2; Table 2).

FIGURE 2

Models and Methods

The source parameters of the Mw 6.0 Hutubi earthquake, according to the published seismic wave, artificial seismic reflection profiles, and aftershocks relocation studies are summarized in Table 3. Where there is a consensus on the strike and rake of the seismogenic fault (all are oriented nearly E-W), there are contradictory characterizations of the fault’s dip angle as either a south-dipping large-scale thrust fault or a north-dipping secondary back-thrust fault (Liu et al., 2018; Lu et al., 2018; Xu et al., 2019).

For the slip model, InSAR observations reveal that the slip distributions are dominated by a nearly pure-thrust fault with no apparent surface rupture, a rupture length of 20 km, and peak slip values ranging from ∼0.1 m to 0.56 m at depths of 12–18 km (Liu et al., 2017; Wang et al., 2019; Yang et al., 2019). According to these published inversion results, we determined that the rupture length, rupture width, and average slip magnitude are 20 km, 15 km, and 0.1 m, respectively. Furthermore, we characterized the seismogenic fault as a blind thrust fault (Lu et al., 2018; Wang et al., 2019; Yang et al., 2019) with a fault tipline at a depth of 7.5 km.

Based on the depth profiles of converted waves from earthquakes in the northern Tien Shan region (Shao et al., 1996), we determined the parameters necessary for our stratified model (Figure 3). By employing PSGRN/PSCMP (Wang et al., 2006), a program that consists of a dislocation source model embedded in a mixed elastic/inelastic layered half-space, we generated six independent stress tensors based on the three different sets of source parameter values for the fault strike, dip, and rake (Table 3). Then we calculated the principal stresses and the direction cosines from the eigenvalues and the eigenvectors of the stress tensors, respectively. From this data, we quantified the stress fields at the four borehole stations.

FIGURE 3

Results

Characteristics of the Coseismic Principal Stress Field

Figure 4 shows the characteristic coseismic principal stress field for the three different sets of source parameters. In model 1, the orientations of the maximum horizontal principal compressive stress at the QEG, SC, YSG, and BLT observation stations are NW-SE, NW-SE, NEE-SWW, and NEE-SWW, respectively. In model 2, the orientations of the maximum horizontal principal compressive stress at the QEG, SC, YSG, and BLT observation stations are NWW-SEE, nearly E-W, NW-SE, and NEE-SWW, respectively. In model 3, the orientations of the maximum horizontal principal compressive stress at the QEG, SC, YSG, and BLT observation stations are NE-SW, nearly E-W, NNE-SSW, and NEE-SWW, respectively.

FIGURE 4

Characteristics of the Coseismic Strain Field

By combining the maximum coseismic principal stress field and the known azimuths of the four components of the borehole strainmeters (Table 2), we quantified the tensile and compressive characteristics of the coseismic strain steps (Table 4) and compared those results to the coseismic steps recorded at the borehole stations (Table 2). In model 1, the predicted tensile and compressive step characteristics match the observed step characteristics on 11/16 channels. In model 2, the predicted tensile and compressive step characteristics match the observed step characteristics on 10/16 channels. In model 3, the predicted tensile and compressive characteristics match the observed step characteristics on 13/16 channels (Table 4). Based on our analysis, we conclude source parameters of 292°/62°/80° (strike/dip/rake) best reproduce the observed strain steps, suggesting that the Mw 6.0 Hutubi earthquake seismogenic fault is a north-dipping high-angle blind thrust fault.

Discussion

Near-field geodetic observations provide the basic data necessary to characterize coseismic deformation (Johnston et al., 2006; Nykolaishen et al., 2015) and numerically simulate both the fault rupture plane and earthquake-induced deformation (Shen et al., 2009; Jiang et al., 2020).

Reliability Analysis of the Results

Our observed tensile and compressive coseismic strain steps are consistent with those in published studies (Gong et al., 2019). To validate the reliability of our simulation results, we calculated the coseismic vertical displacement field, the maximum principal stress field at the surface, and relations between the Coulomb stress change caused by the mainshock for the best-fit source parameters and aftershocks distribution. The coseismic vertical displacement field shows that when the seismogenic fault experienced bilateral rupturing, the hanging wall exhibited coseismic vertical uplift while the footwall exhibited coseismic subsidence (Figure 5). The projected hanging wall displacement is greater than that of the footwall (Figure 5A); these displacement patterns are consistent with those inferred from cross-fault profiles based on InSAR observations (Wang et al., 2019; Yang et al., 2019). When the fault plane dips to the south, the modeled vertical coseismic displacement is a poor fit for the InSAR observations (Figure 5B). From these results, we conclude that the seismogenic fault is a north-dipping thrust fault.

FIGURE 5

The maximum principal stress field shows that the coseismic stress changes are relatively large at stations QEG, SC, and BLT, while the changes recorded at station YSG are relatively small (Figure 6). These coseismic stress changes are consistent with the observed coseismic strain steps (Table 2), and we attribute the spatial variation in the observed borehole strains to the heterogeneity of the coseismic stress changes.

FIGURE 6

Using the aftershocks relocation depths (Xu et al., 2019) in conjunction with the best-fit source parameters, we determined that the static coseismic Coulomb stress change at the two ends of the rupture plane at a depth of 9 km exceeds 0.01 MPa, suggesting that the mainshock hastens or triggers the aftershocks (Figure 7). However, for some aftershocks that occurred on or near the rupture plane, this seemingly contradictory relation may be driven by afterslip (Ross et al., 2017; Perfettini et al., 2019). This observation is consistent with the findings of Wang et al. (2019) and provides further support for the theory that the seismogenic fault generated by the Mw 6.0 Hutubi earthquake is a north-dipping blind thrust fault.

FIGURE 7

Tectonic Implications

Due to compression caused by the collision of the NS-trending Indian plate and the rigid Junggar Basin, the zone between the northern Tien Shan and the Junggar Basin contains many fold-and-thrust structures (Burchfiel et al., 1999; Deng et al., 1999), including listric south-dipping thrust faults and a set of north-dipping high-angle back-thrust faults at the Frontier of the nappe (Liu et al., 2007; Guan et al., 2016). Similar to the categorization of the seismogenic Horgos-Manas-Tugulu fault of the 1906 Mw 7.7 Manas earthquake as a large-scale listric thrust fault (Zhang et al., 1994; Wang et al., 2004), structural studies indicate that the back-thrust faults related to the anticline located in the fold belts near the Hutubi earthquake epicenter are all blind thrust faults (Chen et al., 2007; Lu et al., 2018). Modeling results indicate that the continued India-Eurasia collision resulted in the reactivation of the Tien Shan mountain range and its related fault system (Neil and Houseman, 1997; Li et al., 2018); the Mw 6.0 Hutubi earthquake occurred on a reactivated rift-related structure that is connected to the Horgos-Manas-Tugulu thrust fault.

Double-difference relocation of the aftershocks demonstrate that the Hutubi earthquake seismogenic fault, as a north-dipping blind back-thrust fault that has a thrust direction that opposes the general sense of thrust in the piedmont thrust nappe of the northern Tien Shan, is located near the Horgos-Manas-Tugulu fault (Liu et al., 2018; Xu et al., 2019). This analysis is essentially consistent with our findings where the north-dipping fault plane generates a coseismic principal stress field that generates the observed changes in the near-field coseismic tensile and compressive strain steps. Thus, it is likely that the Hutubi earthquake seismogenic fault is a high-angle blind back-thrust fault with a north-dipping fault plane.

A back-thrust fault, which can rupture independently (Chuang et al., 2014; He et al., 2018), is regarded as a secondary structure in the thrust nappe tectonic belt (Stagpoole and Nicol, 2008; Xu et al., 2015; Zelilidis et al., 2016; Martinod et al., 2020). Based on the occurrence of the 1906 Mw 7.7 Manas earthquake and the 2016 Mw 6.0 Hutubi earthquake, we conclude that the Late Cenozoic crustal deformation caused by the continuous collision between Indian and Eurasia plates was accommodated by slip on the Horgos-Manas-Tugulu thrust fault and also on the secondary back-thrust faults.

In order to analyze the influence of the secondary back-thrust fault on crustal deformation, we employ the ANSYS commercial software to construct a two-dimensional elastic finite element model (Figure 8A). The model includes the lateral heterogeneity, the large-scale thrust fault, and the secondary back-thrust fault. The fault geometry and regional properties are referred to the previous studies (Wang et al., 2004; Liu et al., 2018). First, we model the vertical deformation with only the large-scale fault included (Figure 8B); then we analyze the influence of the back-thrust fault on vertical deformation. As shown by the modeled vertical surface displacement in Figure 8C, this back-thrust fault contributes to the significant uplift observed in this region.

FIGURE 8

Furthermore, with the dip angle of the large-scale thrust fault increasing from the deep to the shallow, this geometry change inhibits the slip on the thrust fault (Bilham and England, 2001) and benefits the strain energy accumulation on the back-thrust faults (Figure 9) (Xu et al., 2015; Zhang et al., 2019), although this process is influenced by other factors such as confining pressure, fault frictional coefficient and rheological property (Neil and Houseman, 1997; Collettini et al., 2019). The two fault systems jointly absorbed or facilitated the crustal shortening (Zubovich et al., 2010; Poblet and Lisle, 2011) and are responsible for the present-day uplift of the northern Tien Shan region that expressed as a ‘pop-up’ structure.

FIGURE 9

Conclusion

In this paper, we used the near-field borehole data to analyze the tensile and compressive characteristics of the coseismic strain steps for the Mw 6.0 Hutubi earthquake. After constructing a model based on elastic dislocation theory, we simulated the coseismic strain step characteristics for three different sets of source parameters. Our model results indicate that when the seismogenic fault source parameters are 292°/62°/80° (strike/dip/slip), the predicted and observed tensile and compressive characteristics of the coseismic strain steps are fairly consistent with one another. With these source parameters, the Hutubi earthquake is likely a high-angle blind back-thrust fault with a north-dipping fault plane. Based on additional data, we inferred that this seismogenic fault is located near the Horgos-Manas-Tugulu thrust fault, which is the seismogenic fault of the 1906 Mw 7.7 Manas earthquake.

Our results, which are informed by both near-field borehole strain observational data and numerical simulations, provide new insight into the analysis of the geometry of seismogenic fault. With more high-quality geodetic and seismic data, it may be possible to elucidate the seismotectonic and tectonic deformation of the northern Tien Shan region using high-resolution three-dimensional models.

Funding

This work was jointly supported by the National Natural Science Foundation of China (Grant Nos. 41874116, 41574081), the Research Grant from the Institute of Crustal Dynamics, China Earthquake Administration (No. ZDJ2019-16).

Acknowledgments

We appreciate Renqi Lu for providing the seismotectonic model of the Hutubi earthquake and Shanshan Liang for the relocated aftershocks. We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript. Figures were plotted using the Generic Mapping Tools (Wessel et al., 2013).

References

• 1

  AllenM. B.VincentS. J.WheelerP. J. (1999). Late Cenozoic tectonics of the Kepingtage thrust zone: interactions of the Tien Shan and Tarim Basin, northwest China. Tectonics18, 639–654. 10.1029/1999tc900019

  • CrossRef
  • Google Scholar

• 2

  AsaiY.OkuboM.IshiiH.AokiH.YamauchiT.KitagawaY.et al (2005). Co-seismic strain-steps associated with the 2004 off the Kii peninsula earthquakes-Observed with Ishii-type borehole strainmeters and quartz-tube extensometers. Earth Planets Space.57, 309–314. 10.1186/bf03352568

  • CrossRef
  • Google Scholar

• 3

  AvouacJ. P.TapponnierP.BaiM.YouH.WangG. (1993). Active thrusting and folding along the northern Tien Shan and late Cenozoic rotation of the Tarim relative to Dzungaria and Kazakhstan. J. Geophys. Res.98, 6755–6804. 10.1029/92jb01963

  • CrossRef
  • Google Scholar

• 4

  BarbourA. J.AgnewD. C.WyattF. K. (2015). Coseismic strains on plate boundary observatory borehole strainmeters in southern California. Bull. Seismol. Soc. Am.105, 431–444. 10.1785/0120140199

  • CrossRef
  • Google Scholar

• 5

  BilhamR.EnglandP. (2001). Plateau ‘pop-up’ in the great 1897 Assam earthquake. Nature410, 806–809. 10.1038/35071057https://pubmed.ncbi.nlm.nih.gov/11298446/

  • CrossRef
  • Google Scholar

• 6

  BurchfielB. C.BrownE. T.QidongD.XianyueF.JunL.MolnarP.et al (1999). Crustal shortening on the margins of the Tien Shan, Xinjiang, China. Int. Geol. Rev.41, 665–700. 10.1080/00206819909465164

  • CrossRef
  • Google Scholar

• 7

  CharvetJ.ShuL.Laurent-CharvetS.WangB.FaureM.CluzelD.et al (2011). Palaeozoic tectonic evolution of the Tianshan belt, NW China. Sci. China Earth Sci.54, 166–184. 10.1007/s11430-010-4138-1

  • CrossRef
  • Google Scholar

• 8

  ChenK.GumiauxC.AugierR.ChenY.MeiY. H.LinW.et al (2014). A multidiscipline method of geological survey, seismic profile and gravity measurement applied to fold-and-thrust belt：a case study along the Hutubi River in the northern piedmont of Tianshan. Chin. J. Geophys.57, 75–87. 10.6038/cjg20140108

  • CrossRef
  • Google Scholar

• 9

  ChenS. P.QiJ. F.YuF. S.YangQ. (2007). Deformation charateristics in the southern margin og the Junggar Basin and their controlling factors. Acta Geol. Sin.81, 151–157.

  • Google Scholar

• 10

  ChuangR. Y.JohnsonK. M.KuoY.-T.WuY.-M.ChangC.-H.KuoL.-C. (2014). Active back thrust in the eastern Taiwan suture revealed by the 2013 Rueisuei earthquake: evidence for a doubly vergent orogenic wedge?Geophys. Res. Lett.41, 3464–3470. 10.1002/2014gl060097

  • CrossRef
  • Google Scholar

• 11

  CollettiniC.TeseiT.ScuderiM. M.CarpenterB. M.VitiC. (2019). Beyond Byerlee friction, weak faults and implications for slip behavior. Earth Planet Sci. Lett.519, 245–263. 10.1016/j.epsl.2019.05.011

  • CrossRef
  • Google Scholar

• 12

  DengQ.FengX.ZhangP.XuX.YangX.PengS.et al (2000). Active tectonics of Tian shan. Beijing: Seismological Press.

  • Google Scholar

• 13

  DengQ.FengX.ZhangP.YangX.XuX.PengS.et al (1999). Reverse fault and fold zone in the Urumqi range-front depression of the northern Tianshan and its genetic mechanism. Earth Sci. Front.6, 191–201.

  • Google Scholar

• 14

  GladwinM. T.GwytherR. L.HartR. H. G.BreckenridgeK. S. (1994). Measurements of the strain field associated with episodic creep events on the san andreas fault at San Juan Bautista, California. J. Geophys. Res.99, 4559–4565. 10.1029/93jb02877

  • CrossRef
  • Google Scholar

• 15

  GongZ.JingY.LiH.LiL.FanX.LiuZ. (2019). Static-dynamic strain response to the 2016 M6.2 Hutubi earthquake (eastern Tien Shan, NW China) recorded in a borehole strainmeter network. J. Asian Earth Sci.183, 103958. 10.1016/j.jseaes.2019.103958

  • CrossRef
  • Google Scholar

• 16

  GuanS.StockmeyerJ. M.ShawJ. H.PleschA.ZhangJ. (2016). Structural inversion, imbricate wedging, and out-of-sequence thrusting in the southern Junggar fold-and-thrust belt, northern Tian Shan, China. Bulletin100, 1443–1468. 10.1306/04041615023

  • CrossRef
  • Google Scholar

• 17

  HeP.HetlandE. A.NiemiN. A.WangQ.WenY.DingK. (2018). The 2016 Mw 6.5 Nura earthquake in the Trans Alai range, northern Pamir: possible rupture on a back-thrust fault constrained by Sentinel-1A radar interferometry. Tectonophysics749, 62–71. 10.1016/j.tecto.2018.10.025

  • CrossRef
  • Google Scholar

• 18

  JiangG.QiaoX.WangX.LuR.LiuL.YangH.et al (2020). GPS observed horizontal ground extension at the Hutubi (China) underground gas storage facility and its application to geomechanical modeling for induced seismicity. Earth Planet Sci. Lett.530, 115943. 10.1016/j.epsl.2019.115943

  • CrossRef
  • Google Scholar

• 19

  JohnstonM. J. S.BorcherdtR. D.LindeA. T.GladwinM. T. (2006). Seismomagnetic effects from the long-awaited 28 september 2004 M 6.0 parkfield earthquake. Bull. Seismol. Soc. Am.96, S206–S220. 10.1785/0120050822

  • CrossRef
  • Google Scholar

• 20

  KingG. C.SteinR. S.LinJ. (1994). Static stress changes and the triggering of earthquakes. Bull. Seismol. Soc. Am.84, 935–953.

  • Google Scholar

• 21

  LiY.WeiD.TianH.JiaD.FangL.YanB.et al (2018). The 3-D structural model of an out-of-sequence earthquake in China: Implication for the reactivation of positive inversion structures along the Northern Tianshan fold-and-Thrust belt. Tectonics37, 4359–4376. 10.1029/2018tc005075

  • CrossRef
  • Google Scholar

• 22

  LindeA. T.GladwinM. T.JohnstonM. J. S.GwytherR. L.BilhamR. G. (1996). A slow earthquake sequence on the San Andreas fault. Nature383, 65–68. 10.1038/383065a0

  • CrossRef
  • Google Scholar

• 23

  LiuB. J.ShenJ.ZhangX. K.ChenY.FangS. M.SongH. P.et al (2007). The crust structures and tectonics of ürümqi depression revealed by deep seismic reflection profile in the northern margin of tianshan mountains. Chin. J. Geophys.50, 1464–1472. 10.1002/cjg2.1146

  • CrossRef
  • Google Scholar

• 24

  LiuD. Q.LiuM.WangH. T.LiJ.ChengJ.WangX. Q. (2016). Slip rates and seismic moment deficits on major faults in the Tianshan region. Chin. J. Geophys.59, 1647–1660. 10.6038/cjg20160510

  • CrossRef
  • Google Scholar

• 25

  LiuG.QiaoX.XiongW.ZhouY.NieZ.XiaC. (2017). Source models for the 2016 M w6.0 Hutubi earthquake, Xinjiang, China: a possible reverse event. Geodesy Geodyn.8, 311–318. 10.1016/j.geog.2017.05.005

  • CrossRef
  • Google Scholar

• 26

  LiuJ. M.WangQ.LiJ.WuC. Y.ZhaoB. B.KongX. Y. (2018). Relocation of the Hutubi Ms6.2 earthquake sequence on 8 December 2016 and analysis of the seismogenic structure. Seismol. Geol.40, 566–578. 10.3969/j.issn.0253-4967.2018.03.005

  • CrossRef
  • Google Scholar

• 27

  LuR.HeD.XuX.WangX.TanX.WuX. (2018). Seismotectonics of the 2016 M 6.2 Hutubi earthquake: implications for the 1906 M 7.7 Manas earthquake in the Northern Tian Shan Belt, China. Seismol Res. Lett.89, 13–21. 10.1785/0220170123

  • CrossRef
  • Google Scholar

• 28

  MartinodJ.GéraultM.HussonL.RegardV. (2020). Widening of the Andes: an interplay between subduction dynamics and crustal wedge tectonics. Earth Sci. Rev.204, 103170. 10.1016/j.earscirev.2020.103170

  • CrossRef
  • Google Scholar

• 29

  MolnarP.TapponnierP. (1975). Cenozoic tectonics of Asia: effects of a continental collision: features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision. Science189, 419–426. 10.1126/science.189.4201.419https://pubmed.ncbi.nlm.nih.gov/17781869/

  • CrossRef
  • Google Scholar

• 30

  NeilE. A.HousemanG. A. (1997). Geodynamics of the Tarim Basin and the Tian Shan in central Asia. Tectonics16, 571–584. 10.1029/97tc01413

  • CrossRef
  • Google Scholar

• 31

  NykolaishenL.DragertH.WangK.JamesT. S.SchmidtM. (2015). GPS observations of crustal deformation associated with the 2012 Mw 7.8 Haida Gwaii earthquake. Bull. Seismol. Soc. Am.105, 1241–1252. 10.1785/0120140177

  • CrossRef
  • Google Scholar

• 32

  OkadaY. (1985). Surface deformation due to shear and tensile faults in a half-space. Bull. Seismol. Soc. Am.75, 1135–1154.

  • Google Scholar

• 33

  OuyangZ.ZhangZ.ShuG. L. (2004). Review and outlook for west China borehole strainmeter netorks. Chin. J. Rock Mech. Eng.23, 4058–4062.

  • Google Scholar

• 34

  PerfettiniH.FrankW. B.MarsanD.BouchonM. (2019). Updip and along‐strike aftershock migration model driven by afterslip: application to the 2011 Tohoku‐Oki aftershock sequence. J. Geophys. Res. Solid Earth.124, 2653–2669. 10.1029/2018jb016490

  • CrossRef
  • Google Scholar

• 35

  PobletJ.LisleR. J. (2011). Kinematic evolution and structural styles of fold-and-thrust belts, Geol. Soc., London, Spec. Publ.349, 1–24. 10.1144/sp349.1

  • CrossRef
  • Google Scholar

• 36

  QiuJ.RaoG.WangX.YangD.XiaoL. (2019). Effects of fault slip distribution on the geometry and kinematics of the southern Junggar fold-and-thrust belt, northern Tian Shan. Tectonophysics772, 228209. 10.1016/j.tecto.2019.228209

  • CrossRef
  • Google Scholar

• 37

  QiuZ.ShiY. (2003). Observations of remote coseismic stress step-changes. Sci. China Earth Sci.46, 75–81. 10.1360/03dz0007

  • CrossRef
  • Google Scholar

• 38

  QiuZ.TangL.ZhangB.GuoY. (2013). In situ calibration of and algorithm for strain monitoring using four-gauge borehole strainmeters (FGBS).J. Geophys. Res. Solid Earth.118, 1609–1618. 10.1002/jgrb.50112

  • CrossRef
  • Google Scholar

• 39

  RossZ. E.RollinsC.CochranE. S.HaukssonE.AvouacJ.-P.Ben-ZionY. (2017). Aftershocks driven by afterslip and fluid pressure sweeping through a fault-fracture mesh. Geophys. Res. Lett.44, 8260–8267. 10.1002/2017gl074634

  • CrossRef
  • Google Scholar

• 40

  ShaoX.ZhangJ.FanH.ZhengJ. D.XuY.ZhangH. Q.et al (1996). The crust structures of Tianshan orogenic belt: a deep sounding work by converted waves of earthquakes along Urumqi-Korla profile. Chin. J. Geophys.39, 336–346.

  • Google Scholar

• 41

  ShenZ.-K.SunJ.ZhangP.WanY.WangM.BürgmannR.et al (2009). Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake. Nat. Geosci.2, 718–724. 10.1038/ngeo636

  • CrossRef
  • Google Scholar

• 42

  StagpooleV.NicolA. (2008). Regional structure and kinematic history of a large subduction back thrust: Taranaki Fault, New Zealand. J. Geophys. Res.: Solid Earth113, B01403. 10.1029/2007jb005170

  • CrossRef
  • Google Scholar

• 43

  WangC.YangZ. E.LuoH.MooneyW. D. (2004). Crustal structure of the northern margin of the eastern Tien Shan, China, and its tectonic implications for the 1906 M?7.7 Manas earthquake. Earth Planet Sci. Lett.223, 187–202. 10.1016/j.epsl.2004.04.015

  • CrossRef
  • Google Scholar

• 44

  WangM.ShenZ.-K. (2020). Present-day crustal deformation of continental China derived from GPS and its tectonic implications. J. Geophys. Res.: Solid Earth.125, e2019JB018774. 10.1029/2019jb018774

  • CrossRef
  • Google Scholar

• 45

  WangR.Lorenzo-MartínF.RothF. (2006). PSGRN/PSCMP-a new code for calculating co- and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory. Comput. Geosci.32, 527–541. 10.1016/j.cageo.2005.08.006

  • CrossRef
  • Google Scholar

• 46

  WangX.ZhangJ.JiangW.WangD.XuG.XiaoZ.et al (2019). Gravity field and deep seismogenic environment in the Longmen Shan and adjacent regions, Eastern Tibetan Plateau. J. Asian Earth Sci.176, 79–87. 10.1016/j.jseaes.2019.01.011

  • CrossRef
  • Google Scholar

• 47

  WesselP.SmithW. H. F.ScharrooR.LuisJ.WobbeF. (2013). Generic mapping tools: improved version released. Eos Trans. AGU.94, 409–410. 10.1002/2013eo450001

  • CrossRef
  • Google Scholar

• 48

  XiaoW.-J.ZhangL. C.QinK. Z.SunS.LiJ. L. (2004). Paleozoic accretionary and collisional tectonics of the Eastern Tianshan (China): implications for the continental growth of central Asia. Am. J. Sci.304, 370–395. 10.2475/ajs.304.4.370

  • CrossRef
  • Google Scholar

• 49

  XuZ.-Q.DilekY.YangJ.-S.LiangF.-H.LiuF.BaD.-Z.et al (2015). Crustal structure of the Indus-Tsangpo suture zone and its ophiolites in southern Tibet. Gondwana Res.27, 507–524. 10.1016/j.gr.2014.08.001

  • CrossRef
  • Google Scholar

• 50

  XuZ. G.LiangS. S.ZouL. Y.LiuJ. G. (2019). Mechanism of the 2016 Hutubi, Xinjiang, Ms6.2 mainshock and relocation of its aftershock sequences. Seismol. Geol.41, 44–57. 10.3969/j.issn.0253-4967.2019.01.003

  • CrossRef
  • Google Scholar

• 51

  YangS.LiJ.WangQ. (2008). The deformation pattern and fault rate in the tianshan mountains inferred from GPS observations. Sci. China Ser. D-Earth Sci.51, 1064–1080. 10.1007/s11430-008-0090-8

  • CrossRef
  • Google Scholar

• 52

  YangW.ChengJ.YaoQ.CuiR. S.LongH. Y.HanY. Y. (2018). The preliminary study on the seismogenic structure of the Hutubi Ms6.2 earthquake. Seismol. Geol.40, 1100–1114.

  • Google Scholar

• 53

  YangX.LiA.HunagW. (2013). Uplift differential of active fold zones during the late quaternary, northern piedmonts of the tianshan mountains, China. Sci. China Earth Sci.56, 794–805. 10.1007/s11430-012-4531-z

  • CrossRef
  • Google Scholar

• 54

  YangY. H.HuJ. C.ChenQ.WangZ. G.TsaiM. C. (2019). A blind thrust and overlying folding earthquake of the 2016 Mw 6.0 Hutubi earthquake in the Northern Tien Shan fold‐and‐thrust belts, China. Bull. Seismol. Soc. Am.109, 770–779. 10.1785/0120180150

  • CrossRef
  • Google Scholar

• 55

  YinA.NieS.CraigP.HarrisonT. M.RyersonF. J.XianglinQ.et al (1998). Late Cenozoic tectonic evolution of the southern Chinese Tian Shan. Tectonics17, 1–27. 10.1029/97tc03140

  • CrossRef
  • Google Scholar

• 56

  ZelilidisA.PapatheodorouG.MaravelisA. G.ChristodoulouD.TserolasP.FakirisE.et al (2016). Interplay of thrust, back-thrust, strike-slip and salt tectonics in a fold and thrust belt system: an example from Zakynthos Island, Greece. Int. J. Earth Sci.105, 2111–2132. 10.1007/s00531-016-1299-y

  • CrossRef
  • Google Scholar

• 57

  ZhangP.DengQ.YangX.PengS.XuX.FengX. (1996). Late cenozoic tectonic deformation and mechanism along the tianshan mountain, northwestern China. Earthq. Res. China.12, 127–140.

  • Google Scholar

• 58

  ZhangP.DengQ.ZhangG.MaJ.GanW.MinW.et al (2003). Active tectonic blocks and strong earthquakes in the continent of China. Sci. China Earth Sci.46, 13–24. 10.1360/03dz0002

  • CrossRef
  • Google Scholar

• 59

  ZhangP. Z.DengQ. D.XuX. W.FengX. Y.PengS. Z.YangX. P.et al (1994). Blind thrust, folding earthquake, and the 1906 Manas earthquake, Xinjing. Seismol. Geol.16, 193–204.

  • Google Scholar

• 60

  ZhangZ.ZhangH.WangL.ChengH.ShiY. (2019). Late cenozoic structural deformation and evolution of the central-southern Longmen Shan fold-and-thrust belt, China: insights from numerical simulations. J. Asian Earth Sci.176, 88–104. 10.1016/j.jseaes.2019.01.026

  • CrossRef
  • Google Scholar

• 61

  ZhouY.WuC.YuanB.WangJ.ZhouT.WangY.et al (2020). Cenozoic tectonic patterns and their controls on growth strata in the northern Tianshan fold and thrust belt, northwest China. J. Asian Earth Sci.198, 104237. 10.1016/j.jseaes.2020.104237

  • CrossRef
  • Google Scholar

• 62

  ZubovichA. V.WangX. Q.ScherbaY. G.SchelochkovG. G.ReilingerR.ReigberC.et al (2010). GPS velocity field for the Tien Shan and surrounding regions. Tectonics29, TC6014. 10.1029/2010tc002772

  • CrossRef
  • Google Scholar