The Haiyuan-Liupanshan fault, an active tectonic feature at the Tibetan Plateau's northeastern boundary, was ruptured by two M8 earthquakes (1920 and 1927) bracketing an unbroken section (the Tianzhu seismic gap). A high seismic hazard is expected along the gap.To monitor deformation characteristics and do a seismic risk assessment, we made measurements at two newly built campaign-mode Global Positioning System (GPS) stations and 13 pre-existing stations in 2013 and 2014.Adding existing data from 1999 to 2014, we derived a new velocity field.Based on the horizontal velocity, we used three block models to invert the deformation of four crustal blocks. The results suggest non-uniform deformation in the interior of the Lanzhou block, the Ordos block and the Alaxan block, but uniform deformation in the Qilian block.Fault slip rates derived from block models show a decreasing trend from west to east, (2.0-3.2 mm/a on the Haiyuan fault to 0.9-1.5 mm/a on the Liupanshan fault).The Haiyuan fault evidences sinistral striking-slip movement, while the Liupanshan fault is primarily thrusting due to transformation of the displacement between the strike-slip and crustal shortening.The locking depth of each segment along the Haiyuan fault obtained by fitting the fault parallel velocities varies drastically from west to east (21.8-7.1 km).The moment accumulation rate, calculated using the slip rate and locking depth, is positively correlated with the locking depth.Given the paucity of large seismic events during the previous millennium, the Tuolaishan segment and the Maomaoshan segment have higher likelihood of nucleation for a future event.
The Haiyuan fault system is part of a major left-lateral fault system in the northeastern margin of the Tibetan Plateau that connects the seismically active Qilian Shan in the west to the tectonically active Liupan Shan in the east, which abuts against the relatively stable Ordos block (Molnar and Tapponnier 1975; Tapponnier and Molnar 1977; Li et al. 2009) (Fig. 1). The fault accommodates part of the deformation induced by the Indian-Asian collision, but no consensus has been reached about the best description of the widespread deformation that occurs there (Cavalié et al. 2008; Thatcher 2007). Neotectonic (Burchfiel et al. 1991;
Zhang et al. 1991; Gaudemer et al. 1995; Yuan et al. 1998;
Lasserre et al. 1999; He et al. 2000; Tian et al. 2001; Li et al. 2009) and geodetic studies (Meade 2007; Thatcher 2007; Cavalié et al. 2008; Wang et al. 2009; Duvall and Clark 2010; Loveless and Meade 2011; Jolivet et al. 2013;
Li et al. 2013) of the Haiyuan fault estimated the Holocene and present-day slip rate, which contributes to a better understanding of the mechanisms of continental deformation and the seismic hazard of the fault. Although the crustal deformation in the northeastern Tibetan Plateau has been studied and characterized, issues still exist, e.g., the fault slip rates obtained by different methods are different from each other and the largest differences are of an order of magnitude, the near-field deformation is relatively poorly constrained by previous Global Positioning System (GPS) data and block models seldom consider non-uniform internal strain. In addition, four historical earthquakes that occurred 3500–3900 years ago suggest an average recurrence interval of approximately 1000 years, corresponding to the time elapsed since the last major event (Liu-Zeng et al. 2007). Thus the Tianzhu gap is the most probable site for the next large rupture on the fault and an assessment of its seismic hazard is an important goal (Cavalié et al. 2008).
Figure
1.
Tectonic map of the northeastern Tibetan Plateau. The inset map on the upper right shows the study area and division of the blocks. Black circles represent earthquakes with magnitude equal to or less than 5.0 and white circles correspond to earthquakes with a magnitude greater than 5.0 during the period 1920–2015 (USGS). Orange lines represent surface ruptures of the 1920 M8 and 1927 M8–8.3 earthquakes, respectively.The red line follows the Tianzhu seismic gap (Gaudemer et al. 1995). Rectangles represent segments with high seismic hazard. TLS Tuolanshan, LLL Lenglongling, JQH Jinqianghe, MMS Maomaoshan, LHS Laohushan, HYW Western Haiyan, HYM Middle Haiyuan, HYE Eastern Haiyuan, LPN Northern Liupanshan, LPS Southern Liupanshan; QL Qilian block, LZ Lanzhou block
For this paper, we collected GPS data from Crustal Movement Observation Network of China (CMONOC Ⅰ) and Continental Tectonic Environmental Monitoring Network of China (CMONOC Ⅱ) from 1999–2014. Sixteen campaign-mode GPS stations in the near-field of the Haiyuan fault were established in 2013 and two of these were surveyed during 2013–2014. Moreover, 13 campaign-mode GPS stations from CMONOC were occupied two epochs in the period of 2013–2014 (Fig. 2). A denser velocity field in the northeastern Tibetan Plateau was derived. A rigid model, a uniform strain model and a linear strain model were applied to invert the deformation of the blocks. The fault slip rate, the locking depth, and the moment accumulation rate of each segment were sequentially inverted. Block deformation and its implication for tectonics in the northeastern Tibetan Plateau are discussed.
Figure
2.
A map showing the locations of the GPS stations around the northeastern margin of the Tibetan Plateau. The black squares are CMONOC campaign stations, and the red stars are the CMONOC continuous stations. The green stars are stations we established in 2013 and monitored in 2013 and 2014, and the pink stars are stations we monitored in 2014
The present-day tectonic features of the ~1000 km Haiyuan-Liupanshan fault are the result of geodynamic processes related to the India-Eurasian convergence some 50 Myr ago (Molnar and Tapponnier 1975; Tapponnier and Molnar 1977). The deformation of the Haiyuan fault is characterized by a left-lateral strike-slip movement (Gaudemer et al. 1995), while the Liupanshan fault is mainly thrusting (Li et al. 2013). The left-lateral strike-slip on the Haiyuan fault has been transformed into crustal shortening and contraction in the eastern end along the Liupanshan (Zheng et al. 2013).
Two large shocks have ruptured the Haiyuan fault during the twentieth century: the M8 Haiyuan earthquake in 1920 and the M8–8.3 Gulang thrust earthquake in 1927 (Zhang et al. 1987; Cavalié et al. 2008) (Fig. 1). Gaudemer et al. (1995) identified a 260 km-long stretch of the fault extending between those two rupture zones as a seismic gap, designated as the Tianzhu gap, which can be divided into four segments (Liu-Zeng et al. 2007) from west to east; (1) the Lenglongling segment, (2) the Jinqianghe segment, (3) the Maomaoshan segment, and (4) the Laohushan segment. The western and eastern ends of the Haiyuan-Liupanshan fault are divided into the Tuolaishan segment, the Haiyuan segment (composed of the western, middle and eastern segments), and the Liupanshan segment (consisting of northern and southern segments) (Li et al. 2009). The Haiyuan segment is the segment ruptured during the 1920 event (Fig. 1).
The Holocene slip rate of the Haiyuan fault was estimated from offset measurements and the dating of morphological markers. The slip rate decreases from 19 ± 5 mm/a along the Lenglongling segment to 12 ± 4 mm/a along the Maomaoshan-Laohushan segments (Lasserre et al. 1999), to 3.5–6.5 mm/a along the Haiyuan segment (4.5 ± 1.1 mm/a (Li et al. 2009), 5–6.5 mm/a (Tian et al. 2001), 4–6 mm/a (Zhang et al. 1991)). However, He et al. (2000) and Yuan et al. (1998) estimated the slip rate along the Lenglongling and Maomaoshan segment is only 4–5 mm/a. Slip rate estimates using block models based on GPS data suggest a 3.2–9 mm/a strike-slip rate along the Haiyuan fault (3.6–5.5 mm/a (Li et al. 2013), 4.6–4.7 mm/a (Loveless and Meade 2011), 7–9 mm/a (Meade 2007), 5–6 mm/a (Thatcher 2007), 3.4 ± 0.2 mm/ a (Wang et al. 2009)), while the slip rates derived from dislocation models tend to be smaller between 1.2 and 5.5 mm/a [2-4.5 mm/a (Duvall and Clark 2010), 1.2-2.2 mm/a (Zhang et al. 2005), 1.4-4.8 mm/a (Ge et al. 2013), 2-5.5 mm/a (Zhang et al. 2011)]. Cavalié et al.
(2008) obtained a 4.2–8 mm/a slip rate by modeling the fault parallel velocity profile along the Maomaoshan-Laohushan segment using ERS InSAR data and inferred the presence of creep at a shallow depth. Additionally, the slip rate of the Liupanshan fault may be less than 1.7 mm/a (Li et al. 2013; Zhang et al. 2011). The cause for the discrepancy between the short term (i.e., geodetic) and the average long-term (i.e., Quaternary) slip rates remains controversial (Jolivet et al. 2012).
3.
Methodology
3.1
GPS measurement and data processing
The principal data used for this study are from the Crustal Movement Observation Network of China (CMONOC Ⅰ) and the Tectonic and Environmental Observation Network of Mainland China (CMONOC Ⅱ) (Gan et al. 2012) and were collected during 1999 and 2014, including 24 continuous stations, and 161 campaign-mode stations with an occupation of seven epochs over the time spanning from 1999 to 2013. Moreover, we conducted observations at 11 stations (including two newly built (2013) stations and nine pre-existing CMONOC stations) in 2013 and 29 stations (including two newly built (2013) stations, 13 pre-existing CMONOC stations and 14 newly built stations) in 2014 across the Haiyuan fault (Fig. 2).
The GPS data were processed together with 20 IGS sites around China using GAMIT/GLOBK10.4 software (Herring et al. 2010; King and Bock 2010). We processed the data in four steps (Shen et al. 2001): (1) All the observational data for a given day were processed together with the data from 20 IGS stations to solve for the daily loosely constrained station coordinates and the satellite orbits, (2) the daily solution for the local stations was combined with loosely constrained global solutions produced at the Scripps Orbital and Position Analysis Center (SOPAC,
http://sopac.ucsd.edu/) using GLOBK software. Time series for the site coordinates were derived. Daily solutions with errors more than 20 mm or three sigma were deleted after removing the linearization tendency, (3) the station positions and velocities were estimated in the ITRF2008 using the GLOBK software. We imposed a frame of reference by minimizing the position and velocity deviations of 34 globally distributed IGS tracking stations with respect to the ITRF2008 while estimating the orientation, translation, and scale transformation parameters, and (4) the velocity solution was transformed into a Eurasia-fixed frame of reference using the Euler vector for Eurasia with respect to the ITRF2008 (Altamimi et al. 2012) (Fig. 3).
Figure
3.
GPS velocities with respect to the stable Eurasian plate. Error ellipses indicate the 95 % confidence levels. Red arrows are sites used to define the self-fixed reference frame. Black rectangles indicate the velocity profiles for Fig. 5. A Tuolaishan segment, B Lenglongling segment, C Jinqianghe segment, D Maomaoshan segment, E Laohushan segment, F Haiyuan western and middle segment, G Haiyuan eastern segment, H Liupanshan northern segment, I Liupanshan southern segment
The northeastern Tibetan Plateau is divided into four blocks based on previous studies (Zhang et al. 2003; Wang et al. 2011), the Ordos block, the Alaxan block, the Qilian block and the Lanzhou block (Figs. 1, 2). The block boundaries coincide to the Haiyuan-Liupanshan fault.
The traditional theory of plate tectonics assumes that a block is a rigid body. The block kinematic can be described by the rigid model (RRM) (Li et al. 2004) using Euler's theorem:
(1)
where () are the longitude and latitude of a point within the block, (Ve, Vn) are the eastern and northern components of the velocity vector, r is the radius of the Earth, and (wx, wy, wz) is the Euler vector.
However, other studies (Li et al. 2001, 2004; Qu et al. 2014) have shown that blocks more closely resemble elastic or visco-elastic entities rather than purely rigid bodies. Deformation exists not only over the boundary zones of the blocks but also within the blocks. GPS velocities are a result of the combined effect of block rotation and internal deformation. If the deformation within a block is homogeneous, the velocity vector of a point within the block could be described using a uniform strain model (HSM) as follows:
(2)
where εe, εen, εne, εn are the parameters that describe the uniform strain within a block and x = cosφ(λ–λ0) and y = r(φ–φ0), (λ0, φ0) are the longitude and latitude of the block's geometric center.
If strain is uneven within a block, the motion of a point within the block could be described by rigid rotation plus the non-uniform deformation of the block. The observed GPS velocity can be defined by a linear strain model (LSM) (Li et al. 2001) as shown below:
(3)
where A0, A1, A2, B0, B1, B2, C0, C1, C2 are the strain parameters.
According to statistical theory (Tao 2007), the unbiasedness and effectiveness of a model are parameters that can be used to evaluate which model is more suitable to delineate block deformation. We used ΔV and SΔv as the criteria for unbiasedness and optimization, respectively:
(4)
(5)
where Δvei and Δvni are the residual velocity components of a GPS point in the eastern and northern directions, n is the number of GPS stations and R is the number of undetermined parameters (Qu et al. 2014).
4.
Results and analyses
4.1
Characteristics of GPS velocity field
The velocity field of the Tibetan Plateau relative to a stable Eurasia includes a rigid rotational component and an interior deformation component (Gan et al. 2007) (Fig. 3). In order to highlight the interior deformation of the northeastern Tibetan, we solved for the Euler vector of 16 GPS stations (Fig. 2) which are stable and evenly distributed within the region. The selection of well-behaved, evenly distributed, and stable stations follows criteria proposed by Rebischung et al. (2012) (Table 1). We also rejected sites with residuals beyond the 1.5 sigma [Eq. (1)]. The resulting Euler pole has a location and magnitude of (61.46° ± 3.12°, 99.19° ± 1.13°, 0.156° ± 0.018°/Ma).
Figure 4 shows the new velocity field after taking out the rigid rotation of the northeastern Tibetan Plateau. Velocities to the south of the Haiyuan-Liupanshan fault show a clockwise rotation, while velocities to the north of the fault zone exhibit a counter-clockwise rotation. The results demonstrate that the Haiyuan-Liupanshan fault is experiencing a left-lateral striking-slip movement with a compression component, which is in agreement with previous studies (Gan et al. 2007; Ge et al. 2013; Wang et al. 2011).
Figure 4 also shows that the fault movement from the Tuolaishan segment to the Laohushan segment is strikeslip with a thrusting component, the Haiyuan segment has mainly left-lateral movement and the Liupanshan segment is dominated by thrusting.
Table
1.
Criteria used to select the northeastern Tibetan-fixed reference stations
Figure
4.
GPS velocity field in a northeastern Tibetan Plateau fixed reference frame. The thick gray solid lines indicate active faults. The gray ellipse at the tip of each velocity vector is the 95 % confidence level
Horizontal velocities show an obvious reduction from the south to the north (about 7.0 mm/a in the south, decreasing to less than 2.0 mm/a near the Haiyuan fault), demonstrating that the crustal strain is accumulating along this deformation zone. Velocities decrease significantly from the outside to the center of the Lanzhou block (Fig. 4), suggesting that an internal strain exists within the block. Velocities also change slightly within the Alaxan and the Qilian blocks, indicating that these blocks are deformable. However, whether an internal strain exists in the Ordos block could not be ascertained from the data in Fig. 4.
4.2
Block deformation and fault slip rate
We calculated SΔv and ΔV for the RRM, HSM, and LSM, respectively (Table 2). The results suggest that the deformation of the Alaxan and the Lanzhou blocks could be described using the HSM and that of the Qilian block using the LSM. However, it is difficult to choose a deformation model to invert the Ordos block between the HSM and the LSM. The applicability of the two models for a specific block can be distinguished with an F test (Qu et al. 2014). We carried out the F test and concluded that the LSM is more suitable for the Ordos block at a significance level of 0.05.
The directions of the residuals shown in Fig. 5 are stochastic in every block for these three block models. Meanwhile, the residuals derived from the LSM are smaller than those from the HSM in the Ordos, the Alaxan and the Lanzhou blocks but larger in the Qilian block, and residuals derived from the RRM are the largest. As a result, it is reasonable to model the Ordos, the Alaxan, and the Lanzhou blocks with the LSM and the Qilian block with the HSM. The residuals histogram conforms to the Gaussian normal distribution, indicating that the models fit the GPS velocities well without systematic errors. It should be noted that the residuals around the block boundaries have a relatively larger magnitude due to the concentration of the crustal deformation (Thatcher 2007; Wang et al. 2011). In other words, strain from interseismic coupling on blockbounding faults has significant impact on sites in the nearfield, especially near block boundaries. Complicated crustal deformation around the blocks boundaries may be beyond the modeling of the current three block models.
Figure
5.
Residuals of the GPS velocities and its histogram. The histogram consists of residuals derived from LSM within the Ordos, the Alaxan, and the Lanzhou block, and from HSM within the Qilian block
The slip rate of each segment was estimated based on block models in terms of the relative movements of faults (Table 4). The RRM was defined as model A in order to obtain slip rates that matched long-term geological results (Zhang et al. 2007). The LSM was used to model the Ordos, the Alaxan, and the Lanzhou blocks, and the HSM was used to delineate the Qilian block. We named the mixed latter as model B (Table 3). Fault parallel and normal slip rate profiles for every segment were also mapped (Fig. 6).
Figure
6.
Profiles along the fault (a) and perpendicular to the fault (b). Blue arc tangent curves are fitted velocity profiles. Red and green rectangles are used to estimate relative movement across faults and black rectangles are used to describe velocities tendency. a Tuolaishan segment, b Lenglongling segment, c Jinqianghe segment, d Maomaoshan segment, e Laohushan segment, f Haiyuan western and middle segment, g Haiyuan eastern segment, h Liupanshan northern segment, i Liupanshan southern segment
The slip rates in Table 4 show that the overall characteristic of the Haiyuan-Liupanshan fault is dominated by a left-lateral striking-slip movement with a compression component. The slip rates derived from model A suggest that the fault activity decreases from the Tuolaishan segment to the southern Liupanshan segment with a reduction in the striking-slip rate from ~4.8 mm/a to 0.3 mm/a and in the shorting rate from ~3 mm/a to 0.7 mm/a. There is an increase in the compression component along the Liupanshan segment due to the transformation of the displacement between Haiyuan and Liupanshan (Zheng et al. 2013). Estimated slip rates from profiles (Table 4; Fig. 6) across the fault are consistent with model A. Results from model B exhibit similar trends to those from model A, with a smaller magnitude. Studies (Zhang et al. 2007; Wei et al. 2012) suggest that the slip rate obtained from the RRM and the velocity profiles probably reflects the far-field deformation of the fault, and yet the slip rate derived from the HSM and the LSM reflects the near-field crustal deformation of the fault.
Velocity gradients across the Haiyuan segment (Fig. 6a, f,g) indicate that the fault has a shallow locking depth and may still suffer from the impact of the 1920 M8.0 earthquake (Li et al. 2013; Meade and Hager 2005; Savage and Prescott 1978). Velocity gradients are also found on the western side of the Liupanshan segment (Fig. 6b, h, i). The compression component decreases from far-field (100–200 km) to near-field of the fault (50–10 km) with a slip rate from ~8 mm/a to ~6 mm/a, demonstrating that the left-lateral strike-slip on the Haiyuan fault is transforming into a thrusting component of the Liupanshan fault and crustal shortening in the region. The relative movement (Fig. 6b, h, i) is very small between the two sides (<50 km) of the Liupanshan segment, suggesting that strain might have accumulated in this area.
Our preferred slip rate for the Haiyuan fault derived from model B (from the Tuolaishan segment to the Haiyuan eastern segment) is about 2.0–3.2 mm/a, which is consistent with the estimate of Ge et al. (2013) of 1.4–3.5 mm/a. Using a finite method, Zhang et al. (2011) modeled GPS velocities in the northeastern Tibetan Plateau and obtained a value for the strike-slip of the Haiyuan fault of 2–5.5 mm/a. A similar slip rate of 2–4.5 mm/a was obtained by analyzing GPS velocities and the Quaternary fault slip of the main faults in northeastern Tibet (Duvall and Clark 2010). Moreover, the inversion of the 1683 GPS velocities within 31 blocks based on a block model by Wang et al. (2009) revealed that the strike-slip rate of the Haiyuan fault is 3.0–3.4 mm/a. Those slip rates correspond well with our results. Moreover, the slip rate of the Liupanshan fault is similar to that found in previous studies (~1.0 mm/a) (Zhang et al. 2011; Li et al. 2013), but with a smaller fault normal component. Discrepancies may arise from different models and datasets.
4.3
Locking depth and moment accumulation rate of faults
A rough estimate of relative motion of the plate across a section where the fault is locked can be obtained from a dislocation model (Savage and Burford 1973):
(6)
where D is the locking depth, b is the long-term slip rate of the fault, v is fault parallel velocity, and x is the distance to the fault.
The fault locking depth was estimated by a least-squares fit of the fault parallel velocities (Table 2). The locking depth in the Liupanshan segment (northern and southern Liupanshan) is none for the model and is not suitable for low-angle faults (Fig. 6a, h, i). The overall characteristic of locking depth from the Tuolaishan segment to the Haiyuan eastern segment is a decrease from 21.8 to 7.1 km, which is consistent with the decreasing trend for crustal thickness along the fault from west to east (Cui et al. 2009a, b). Cui et al. (2009a, b) gave a locking depth of 22.2 and 11.4 km for the Maomaoshan and the Laohushan segments, respectively, and 3.6–8.5 km for the Haiyuan segment based on the Smith-3D dislocation model. A locking depth of 15 km for the Jinqianghe segment is preferred by Li et al. (2013) based on GPS inversion. Our result for the locking depth is consistent with those studies. Moreover, the focal depth of the 1990 Tianzhu-Jingtai M 6.2 earthquake (the largest in decades) (Fig. 1) that occurred on the Maomaoshan is 12 km, which matches our locking depth well. The focal depth of the 1920 Haiyuan M 8 earthquake is generally believed to be between 17 and 20 km, which is different from our fitting locking depth for the Haiyuan segment. We speculate that the lower crust of the Haiyuan segment was not fully stuck until now, which may be the reason that no earthquakes with magnitude greater than 6.0 have occurred here during the past 95 years.
The seismic moment accumulation rate per unit length of the fault was also calculated. The moment accumulation rate M depends on the fault length l, the locking depth dj, the slip rate vj, and the rock shear modulus μ:
(7)
An estimated locking depth was used to calculate the seismic moment accumulation rate per unit length for each segment (Table 4). As expected, the rate of moment accumulation is positively correlated with the locking depth. Those rates can be compared with the stress accumulation rate and the recurrence interval to establish the seismic hazard (Smith and Sandwell 2003). Fault segments that have a high seismic moment accumulation rate are associated with deep locking depths, while faults that have a shallow locking depth have a lower seismic moment accumulation rate and a correspondingly reduced hazard potential (Burgmann 2000). Regarding our understanding of earthquakes, all of the strain that builds up on the fault in the period between earthquakes must be released solely during large earthquakes. Tight fault coupling, high seismic moment accumulation rate, and the paucity of large seismic events over previous decades along the Tuolaishan segment and the Maomaoshan segment may imply more strain accumulation. Considering that the Tianzhu seismic gap corresponding to the time elapsed since the last major event (Liu-Zeng et al. 2007), we infer that a greater possibility of nucleation for a future event on the Tuolaishan segment and the Maomaoshan segment (lying along the Tianzhu seismic gap) exists. It should be emphasized that a complete scientific knowledge of the structure, palaeoseismic history, and dynamics of one fault are needed to our sense of security in terms of its seismic hazard (Zhang 2013).
5.
Discussion
5.1
Characteristics of block deformation
GPS measurements that uniquely quantify the present-day deformation of the Tibetan Plateau can be simply and usefully described by the relative motions of the active blocks and the fault slip across the block boundaries (Meade and Hager 2005; Thatcher 2007). Our results show that the block models fit the horizontal GPS velocities well for the northeastern Tibetan Plateau. The Lanzhou block and the Qilian block can be best described by the LSM and the HSM, respectively. The Qilian block is smaller in size and has fewer secondary faults within the block than the Lanzhou block, and it deforms more uniformly. Qu et al.
(2014) used regional plate kinematic models to fit GPS data in the Fenwei basin and suggested that the Ordos block deforms rigidly. On the basis of 1350 GPS velocities,
Zhang et al. (2005) inferred that the Ordos block behaved as a coherent block similar to the oceanic rigid block without internal deformation. Wang et al. (2011) came to a similar conclusion for the Ordos block. In addition, Wang et al. (2009) asserted that the Ordos block and the Alaxan block are obviously rigid and the strength of these blocks is comparable to that of the lithospheric plate. However, our preferred models for the Ordos block (LSM) and the Alaxan block (LSM) are different from those in these studies because our GPS sites within the Ordos and the Alaxan blocks are primarily distributed around the block boundaries where complicated crustal deformation is concentrated (Thatcher 2007; Wang et al. 2011).
Although the relative motion located at or near the block boundaries agrees to the expectations based on the structural mapping of the active faults, the predicted slip rates are generally considerably less than those obtained by geological methods. Moreover, we should keep a rigorous attitude that these three models maybe improved more in mathematical sense than crustal block deformation.
5.2
Implication for tectonic deformation
Two alternative end-number models were proposed to describe the tectonic deformation occurring in the Tibetan and elsewhere on the continent: continuum versus block models (England and Mckenzie 1982; Avouac and Tapponnier 1993). The key divergence of the two models is whether the deformation is distributed broadly or localized on a few major faults (England and Molnar 2005). Kinematic analysis of the Haiyuan fault provides some clues. Although velocities fit block models well, attention should be paid that the low slip rates estimated by the space geodetic methods along the Haiyuan-Liupanshan fault imply that continuum models should be preferred for the northeastern Tibetan Plateau, but this is still controversial (Thatcher 2007; Cavalié et al. 2008). Zhang et at. (2005) argued that the present-day tectonic deformation of Chinese mainland can be described in terms of a coupling model of rigid block movement and continuous deformation. Meanwhile, same set of GPS data could be fitted well by microplate model (Thatcher 2007) and visco-elastic model (Cao et al. 2009). A compromise is that both the GPS data and the low fault slip rates implied are quite compatible with the block motions in the northeastern Tibetan Plateau (Zhang et al. 2005; Thatcher 2007).
6.
Conclusions
On the basis of the measurements made from 185 GPS stations around the northeastern margin of the Tibetan Plateau since 1999, we derived a new velocity field, which plays a more important role in constraining the near-field crustal deformation of the Haiyuan-Liupanshan fault. After removing the rigid rotation of the northeastern Tibetan Plateau, interior deformation demonstrates quite different features bounded by the Haiyuan-Liupanshan fault. Clockwise rotation at the southern side of the fault and counter-clockwise at the northern side as indicated by the GPS velocities suggest that the Haiyuan-Liupanshan fault is undergoing a left-lateral striking-slip movement with a compression component. Block models used to invert the present-day deformation of the northeastern Tibetan Plateau show that active deformation in this region can be explained by crustal block rotation with an internal strain. The slip rates inferred from the block models vary slightly along the Haiyuan fault (2.0–3.2 mm/a) and the Liupanshan fault (0.9–1.5 mm/a), while the locking depth derived from the fault parallel velocities exhibit a large discrepancy along the Haiyuan fault (7.1–21.8 km). Compared to the adjacent area in the seismic gap, the Tuolaishan segment and the Maomaoshan segment with its high rate of moment accumulation coupled with the paucity of large earthquakes during the past~1000 years constitute a high seismic risk in the future.
Acknowledgments
We express our gratitude and thanks to all participants who helped to construct the network and take the GPS measurements. We are grateful to Guo Liming, Xu Xiaobo, Wang Jiaqing, Zuo Ronghu, Zhang Yingfeng and Jiang Yu for taking the GPS measurements in the summer of 2013 and 2014. We express our appreciation to three anonymous reviewers for their pertinent comments. Partial figures were generated with the GMT software (Wessel and Smith 1995). This study is supported by the National Natural Science Foundation of China (No. 41461164002, No. 41374015), the State Key Laboratory of Earthquake Dynamics (LED2013A02), and the Fundamental Research Funds for the Central Universities (R1401038A, 14CX02110A).
Open Access
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Avouac J, Tapponnier P (1993) Kinematic model of active deformation in central Asia. Geophys Res Lett 20(10):895-898 doi: 10.1029/93GL00128
Burchfiel BC, Zhang PZ, Wang YP (1991) Geology of the Haiyuan fault zone, Ningxia-Hui autonomous Region, China, and its relation to the evolution of the northeastern margin of the Tibetan Plateau. Tectonics 10(6):1091-1110. doi: 10.1029/92TC00001
Burgmann R (2000) Earthquake potential along the northern Hayward fault, California. Science 289(5482):1178-1182 doi: 10.1126/science.289.5482.1178
Cao JL, Shi YL, Zhang H, Wang H (2009) Numerical simulation of GPS observed clockwise rotation around the eastern Himalayan syntax in the Tibetan Plateau. Chin Sci Bull 54(8):1398-1410 http://www.cnki.com.cn/Article/CJFDTotal-JXTW200908020.htm
Cavalié O, Lasserre C, Doin MP, Peltzer G, Sun J, Xu X, Shen ZK (2008) Measurement of interseismic strain across the Haiyuan fault (Gansu, China), by InSAR. Earth Planet Sci Lett 275(3-4):246-257 doi: 10.1016/j.epsl.2008.07.057
Duvall AR, Clark MK (2010) Dissipation of fast strike-slip faulting within and beyond northeastern Tibet. Geology 38(3):223-226 doi: 10.1130/G30711.1
England P, Mckenzie D (1982) A thin viscous sheet model for continental deformation. Geophys J R Astron Soc 70(2):295-321 doi: 10.1111/j.1365-246X.1982.tb04969.x
England P, Molnar P (2005) Late Quaternary to decadal velocity fields in Asia. J Geophys Res 110(B12):501-509 doi: 10.1029-2004JB003541/
Gan W, Zhang P, Shen Z, 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(B8):582-596 doi: 10.1029-2005JB004120/
Gan W, Li Q, Zhang R, Shi HB (2012) Construction and application of tectonics and environmental observation network of mainland China. J Eng Stud 4:324-331. (in Chinese with English abstract)
Gaudemer Y, Tapponnier P, Meyer B, Peltzer G, Guo S, Chen Z (1995) Partitioning of crustal slip between linked, active faults in the eastern Qilianshan, and evidence for a major seismic gap, the 'Tianzhu gap', on the western Haiyuan fault, Gansu (China). Geophys J Int 120(3):599-645 doi: 10.1111/j.1365-246X.1995.tb01842.x
Ge WP, Wang M, Shen ZK, Yuan DY, Zheng WJ (2013) Interseismic kinematics and deformation patterns on the upper crust of Qaidam-Qilianshan block. Chin J Geophys 56:2994-3010
Jolivet R, Lasserre C, Doin MP, Peltzer G, Avouac JP, Sun J, Dailu R (2013) Spatio-temporal evolution of aseismic slip along the Haiyuan fault, China: implications for fault frictional properties. Earth Planet Sci Lett 377-378(5):23-33 http://www.sciencedirect.com/science/article/pii/S0012821X13003956
King RW, Bock Y (2010) Documentation for the GAMIT GPS analysis software, release 10.0. Massachusetts Institute of Technology, Cambridge
Lasserre C, Morel PH, Gaudemer Y (1999) Post glacial left slip-rate and past occurrence of M > 8 earthquakes on western Haiyuan fault, Gansu, China. J Geophys Res 104(B8):17633-17651. doi: 10.1029/1998JB900082. (in Chinese with English abstract)
Li YX, Li Z, Zhang JH, Huang C, Zhu WY, Wang M, Guo LQ, Zhang ZF, Yang GH (2004) Horizontal strain field in the Chinese mainland and its surrounding areas. Chin J Geophys 47(2):222-231. (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTOTAL-DQWX200402007.htm
Liu-Zeng J, Klinger Y, Xu X, Lasserre C, Chen G, Chen W, Chen WB, Tapponnier P, Zhang B (2007) Millennial recurrence of large earthquakes on the Haiyuan fault near Songshan, Gansu Province, China. Bull Seismol Soc Am 97(1B):14-34 doi: 10.1785/0120050118
Loveless JP, Meade BJ (2011) Partitioning of localized and diffuse deformation in the Tibetan Plateau from joint inversions of geologic and geodetic observations. Earth Planet Sci Lett 303(1-2):11-24. doi: 10.1016/j.epsl.2010.12.014
Meade BJ (2007) Present-day kinematics at the India-Asia collision zone. Geology 35(1):81 doi: 10.1130/G22924A.1
Molnar P, Tapponnier P (1975) Cenozoic tectonics of Asia: effects of a continental collision. Science 189(4201):419-426 doi: 10.1126/science.189.4201.419
Rebischung P, Griffiths J, Ray J, Schmid R, Collilieux X, Garayt B (2012) IGS08: the IGS realization of ITRF2008. GPS Solut 16(4):483-494 doi: 10.1007/s10291-011-0248-2
Savage JC, Burford RO (1973) Geodetic determination of relative plate motion in central California. J Geophys Res 78(5):832-845 doi: 10.1029/JB078i005p00832
Savage JC, Prescott WH (1978) Asthenosphere readjustment and the earthquake cycle. J Geophys Res 83(B7):3369-3376 doi: 10.1029/JB083iB07p03369
Shen ZK, Wang M, Li Y, Jackson DD, Yin A, Dong DN, Fang P (2001) Crustal deformation along the Altyn Tagh fault system, western China, from GPS. J Geophys Res 106(B12):30607-30621 doi: 10.1029/2001JB000349
Smith B, Sandwell D (2003) Coulomb stress accumulation along the San Andreas fault system. J Geophys Res 108(B6):2296 doi: 10.1029/2002JB002136/pdf
Tao BZ (2007) The statistical theory and methods of the survey measurement data processing. Surveying and Mapping Press, Beijing (in Chinese)
Tapponnier P, Molnar P (1977) Active faulting and tectonics in China. J Geophys Res 82(20):2905-2930 doi: 10.1029/JB082i020p02905
Tian QJ, Ding GY, Shen XH (2001) Pull-apart basins and the total lateral displacement along the Haiyuan fault zone in Cenozoic. Earthq Res China 17(2):167-175 http://www.cqvip.com/Main/Detail.aspx?id=1000531763
Wang W, Yang S, Wang Q (2009) Crustal block rotations in Chinese mainland revealed by GPS measurements. Earthq Sci 22(6):639-649 doi: 10.1007/s11589-009-0639-4
Wang H, Liu M, Cao J, Shen X, Zhang G (2011) Slip rates and seismic moment deficits on major active faults in mainland China. J Geophys Res 116(B2):1161-1172 doi: 10.1029-2010JB007821/
Zhang PZ, Burchfiel BC, Molnar P, Zhang WQ, Jiao DC, Deng QD, Wang YP, Leigh Royden, Song FM (1991) Amount and style of late Cenozoic deformation in the Liupan Shan area, Ningxia autonomous region, China. Tectonics 10(6):1111-1129 doi: 10.1029/90TC02686
Zhang PZ, Molnar P, Xu X (2007) Late quaternary and present-day rates of slip along the Altyn Tagh fault, northern margin of the Tibetan Plateau. Tectonics 26(5):1-8 doi: 10.1029-2006TC002014/
Zheng WJ, Zhang PZ, He WG, Yuan DY, Shao YX, Zheng DW, Ge WP, Min W (2013) Transformation of displacement between strike-slip and crustal shortening in the northern margin of the Tibetan Plateau: evidence from decadal GPS measurements and late Quaternary slip rates on faults. Tectonophysics 584:267-280 doi: 10.1016/j.tecto.2012.01.006
Peng, S., Ding, K., Li, T. et al. The refined slip rates of the Haiyuan fault system from GNSS and levelling data. Geophysical Journal International, 2025, 240(3): 1690-1704.
DOI:10.1093/gji/ggaf013
2.
Li, Y., Shao, Y., Wang, T. et al. Assessing the data quality and seismic monitoring capabilities of the Belt and Road GNSS network. Earthquake Science, 2025, 38(1): 56-66.
DOI:10.1016/j.eqs.2024.09.007
3.
Wang, S., Xu, K., Zhang, M. Past century stress evolution and seismic hazards in the Haiyuan Fault in northeastern Tibetan Plateau: A case of the 2022 Menyuan MS6.9 earthquake. Advances in Space Research, 2025.
DOI:10.1016/j.asr.2025.03.065
4.
Su, X., Bao, Q., Gao, Z. et al. Three-dimensional interseismic crustal deformation in the northeastern margin of the Tibetan Plateau using GNSS and InSAR. Journal of Asian Earth Sciences, 2024.
DOI:10.1016/j.jseaes.2024.106328
5.
Xu, Y., Zeng, X., Luo, G. Relocation and seismogenic structure of the 2021 Guyuan earthquake swarm | [2021 年固原震群重定位及发震构造分析]. Acta Seismologica Sinica, 2024, 46(5): 767-786.
DOI:10.11939/jass.20230006
6.
He, L., Feng, G., Li, Z. et al. Interseismic kinematics along the Tuolaishan-Lenglongling fault determined by Sentinel-1 InSAR observations. Tectonophysics, 2024.
DOI:10.1016/j.tecto.2023.230152
7.
Liu, H., Xie, L., Zhao, G. et al. A joint InSAR-GNSS workflow for correction and selection of interferograms to estimate high-resolution interseismic deformations. Satellite Navigation, 2023, 4(1): 14.
DOI:10.1186/s43020-023-00105-6
8.
Dong, X., Li, Z., Jing, X. et al. Stratigraphic sequence characteristics and geochronology research progress of the Cenozoic in the arcuate tectonic belt on the northeastern margin of the Tibet Plateau | [青藏高原东北缘弧形构造带新生代地层沉积序列及年代学研究进展]. Journal of Geomechanics, 2023, 29(4): 465-484.
DOI:10.12090/j.issn.1006-6616.2023048
9.
Li, X., Hammond, W.C., Pierce, I.K.D. et al. Present-Day Strike-Slip Faulting and Intracontinental Deformation of North China: Constraints From Improved GPS Observations. Geochemistry, Geophysics, Geosystems, 2023, 24(7): e2022GC010781.
DOI:10.1029/2022GC010781
10.
Yu, S., Su, X. A Crustal Deformation Pattern on the Northeastern Margin of the Tibetan Plateau Derived from GPS Observations. Remote Sensing, 2023, 15(11): 2905.
DOI:10.3390/rs15112905
11.
Wu, D., Qu, C., Zhao, D. et al. Slip Models of the 2016 and 2022 Menyuan, China, Earthquakes, Illustrating Regional Tectonic Structures. Remote Sensing, 2022, 14(24): 6317.
DOI:10.3390/rs14246317
12.
Kou, L., Dong, X., Zhenhong, L.I. et al. Initiation and Development of the Late Cenozoic Uplift of Daluo Mountains, Northeastern Margin of the Tibetan Plateau. Acta Geologica Sinica (English Edition), 2022, 96(6): 1917-1931.
DOI:10.1111/1755-6724.14891
13.
Wang, S., Shi, Y., Feng, X. et al. Late Quaternary sinistral strike-slipping of the Liupanshan-Baoji fault zone: Implications for the growth of the northeastern Tibetan Plateau. Geomorphology, 2021.
DOI:10.1016/j.geomorph.2021.107628
14.
Zhu, L., Ji, L., Jiang, F. Variations in Locking Along the East Kunlun Fault, Tibetan Plateau, China, Using GPS and Leveling Data. Pure and Applied Geophysics, 2020, 177(1): 215-231.
DOI:10.1007/s00024-019-02231-2
15.
Duan, H., Zhou, S., Li, R. Estimation of dip angle of Haiyuan faults based on seismic data | [基于地震活动性资料估计海原断裂倾角]. Acta Geophysica Sinica, 2018, 61(9): 3713-3721.
DOI:10.6038/cjg2018L0601
16.
Shi, F., Shao, Z., Zhan, W. et al. Numerical modeling of the shear modulus and stress state of active faults in the northeastern margin of the Tibetan plateau | [青藏高原东北缘活动断裂剪切模量及应力状态数值模拟]. Acta Geophysica Sinica, 2018, 61(9): 3651-3663.
DOI:10.6038/cjg2018L0631
17.
Li, Y., Shan, X., Qu, C. et al. Elastic block and strain modeling of GPS data around the Haiyuan-Liupanshan fault, northeastern Tibetan Plateau. Journal of Asian Earth Sciences, 2017.
DOI:10.1016/j.jseaes.2017.10.010
18.
Hao, M., Li, Y.-H., Qin, S.-L. Spatial and Temporal Distribution of Slip Rate Deficit Across Haiyuan-Liupan Shan Fault Zone Constrained by GPS Data. Dizhen Dizhi, 2017, 39(3): 471-484.
DOI:10.3969/j.issn.0253-4967.2017.03.003
19.
Li, Y., Shan, X., Qu, C. et al. Fault locking and slip rate deficit of the Haiyuan-Liupanshan fault zone in the northeastern margin of the Tibetan Plateau. Journal of Geodynamics, 2016.
DOI:10.1016/j.jog.2016.07.005
DownLoad:
