Home Characteristics of the present crustal deformation in the Tibetan Plateau and its relationship with strong earthquakes
Article Open Access

Characteristics of the present crustal deformation in the Tibetan Plateau and its relationship with strong earthquakes

  • Haiping Ma EMAIL logo , Hui Zhang , Minjuan Li , Shanyi Wu , Pengtao Wang , Qian Wang EMAIL logo , Jing Zhao EMAIL logo and Zhiqiang Ma
Published/Copyright: July 15, 2023
Become an author with De Gruyter Brill

Abstract

To study the characteristics of the present crustal movement in the Tibetan Plateau and explore its relationship between strong earthquakes with magnitudes of 8 and above, the velocity field size was analyzed based on the global position system (GPS) campaign observations and the time series of site north-ward displacement and long baseline were discussed using the GPS continuous observations. The results show that the velocity field size in the Tibetan Plateau decreases from southwest to north, northeast, and southeast, and the value of the velocity in the west is significantly greater than that in the east in the same dimension. The maximum value is located in the southwest and the minimum value is located in the east. The Wenchuan earthquake is located in the mutation region, where the rate and the direction of the crustal movement are quite different. The crustal deformation extent is large in the region close to the seismic source before the earthquake, reflecting that the regional stress accumulation is fast and its required time is relatively short. However, the crustal deformation extent is relatively small in the region away from the seismic source before the earthquake, reflecting that the regional stress accumulation is slow and its required time is relatively long. The N-ward movement became significantly strong after Nepal M S 8.1 earthquake; the occurrence of this earthquake may have caused the unlocking of large-scale faults near the seismic source, which further intensified the NE-ward subduction and collision of the Indian Plate. The compression of the Indian Plate to the Tibetan Plateau slowed down after the 2008 Wenchuan M S 8.0 earthquake, and increased significantly after 2015, which boosted strain accumulation in the Tibetan Plateau, and attention needs to be paid continuously to strong earthquake risk in this region.

1 Introduction

Located in the southern part of the Asian continent and formed by the extrusion of the Indian Ocean plate and the Eurasian Plate, the Tibetan Plateau is the highest, largest, and youngest plateau, whose crustal thickness is about twice the normal crustal thickness and which is one of the regions with the strongest lithospheric deformation and the most concentrated earthquake activity in Asia and even in the world [1,2,3]. The Tibetan Plateau has been uplifting for the past two to three million years, whose strong uplift is the result of the direct collision and extrusion of the Indian Plate to the north and the Asian Plate [4]. Due to its vast territory, a large number of blocks, rapid uplift, thick crust, strong deformation, and complex internal tectonic activities, the Tibetan Plateau has always been a research hotspot for geographers and geophysicists all over the world to study its tectonic background, uplift causes, deformation characteristics, and dynamic mechanism. Abundant achievements have been made in studying crustal deformation characteristics in the Tibetan Plateau at home and abroad. Shen et al. studied the movement and deformation characteristics of the main fault zones on the southeastern boundary of the Tibetan Plateau [5]. Zuza et al. estimated the amount and rate of Cenozoic shortening at the northern margin of the Tibetan Plateau [6]. Copley inferred two obvious normal fault seismic areas in the crust of the southeast edge of the Qinghai Tibet Plateau by using seismic activity, global position system (GPS), and gravity data [7]. The Tibetan Plateau is divided into five blocks structurally from south to north: the Lhasa block, the Qiangtang block, the Bayan Har block, the Qaidam block, and the Qilian block, which were separated by a huge suture zone in the near east–west direction [8,9,10,11,12].

GPS observation can provide the relative movement information of the crust in various spatial scales. Thus, it is widely used in seismic geology and earthquake prediction. Moreover, the density and accuracy of the velocities make it possible to better quantify the crustal deformation characteristic [13]. Previous research on present crustal movements in Tibetan Plateau using GPS campaign data has mainly focused on the change state of velocity field and strain rate field, which has supported a series of studies on the relationship between active blocks, deformation mechanisms of the Tibetan Plateau, strain fields, and strong earthquakes processes [5,11,14,15,16,17,18]. Other researchers conducted GPS time series analysis of the dynamic change of the time process of crustal movement based on the GPS continuous observation results, including GPS station time series analysis and baseline time series between GPS stations, analyzed the state of crustal movement before the occurrence of moderately strong earthquakes, and discussed the relationship between the characteristics [12,19,20,21,22]. However, the relationship between the current crustal movement and the large earthquake with a magnitude of 8.0 or greater is still lacking.

In this paper, the velocity field in the study area was analyzed based on the GPS campaign observation data of the northeastern margin of the Tibetan Plateau from 1999 to 2016. The displacement variation of GPS stations and the baseline variation of extension between stations were analyzed using the GPS continuous observation data. Besides, the characteristics and dynamic changes of crustal movement in the Tibetan Plateau were studied, and the relationship between the temporal and spatial variation of crustal movement and strong earthquakes was discussed. The innovational work include (1) the most densely distributed velocity field of the crust movement in mainland China and its surrounding areas and the latest GPS continuous observation data used to illustrate the motion characteristics of the study area, (2) the method for evaluating the N-direction displacement using the GPS data from large-scale continuous stations, and (3) the state of crustal movement before the occurrence of moderately strong earthquakes and its relations of the strong earthquake with a magnitude of 8.0 and above. The results can provide a reference for determining the dangerous locations of strong earthquakes and importance for mitigating seismic geo-hazards.

2 Study area, GPS data, and methodology

2.1 Study area

The Tibetan Plateau and its marginal areas are the high occurrence areas of large earthquakes with magnitudes greater than 7. The seismogenic structures of these are the large strike slip faults at the boundaries of active blocks, the reverse faults of the crustal shortening at the front of the active blocks, and the normal faults with locally spreading in the trailing edge of active blocks (Figure 1), indicating that the seismicity of the Tibetan Plateau is closely related to the block movement. The blocks are tectonic units with different sizes, motion states and strain states inside the Plate, which are formed from the division of plates by tectonic faults or fault zones. The Tibetan Plateau is divided into five active blocks including the Qilian block, Qaidam block, Bayan Har block, Qiangtang block, and Lhasa block [9].

Figure 1 
                  The spatial distribution of active-tectonic blocks of Tibetan Plateau and the historical strong earthquakes (The red circles are the earthquakes above magnitude 7 since 1900.).
Figure 1

The spatial distribution of active-tectonic blocks of Tibetan Plateau and the historical strong earthquakes (The red circles are the earthquakes above magnitude 7 since 1900.).

2.2 GPS data

According to the observation mode, GPS measurements can be divided into continuous observation and campaign observation modes. The continuous observation mode refers to the installation of GPS equipment on the permanent ground stations, and the data are continuously stored in the receiver storage device or transmitted back to the data center through a network. The main output is GPS time series results. However, the campaign observation adopted the mode of multi-stage retest at certain time intervals (generally longer than a year) on the permanent ground station, and each observation phase is generally approximately 96 h. The main output is the GPS velocity field results. The models of GPS observation equipment are shown in Table 1. Before 2009, the receiver and antenna models of campaign observation and continuous observation modes are the same, the receiver model is ASHTECH Z-XII3, and the antenna model is ASH701945B_M. However, the receiver model of the campaign observation mode after 2009 is mainly TPS NET-G3A, and the antenna model is mainly TPSCR.G3. The receiver model for the continuous observation mode after 2009 is mainly TRIMBLE NETR9 and the antenna model is mainly TRM59800.00. The campaign observation mode ensures the spatial resolution of the crustal deformation research. Although the spatial resolution of the continuous observation mode is low, it makes up for the insufficient temporal resolution of campaign observation and reflects the temporal dynamic process of crustal deformation objectively [12,23].

Table 1

Models of GPS observation equipment

Observation method Receiver Antenna Observing session
Campaign ASHTECH Z-XII3 ASH701945B_M Before 2009
TPS NET-G3A TPSCR.G3 After 2009
Continuous ASHTECH Z-XII3 ASH701945B_M Before 2009
TRIMBLE NETR9 TRM59800.00 After 2009

The GPS velocity field data used in this paper come from the literature [15], which are gathered GPS data from Crustal Movement Observation Network of China (CMONOC) I and II, and more data are added from densified regional campaign GPS networks and regional continuous GPS sites located in seismically active areas. To better analyze the dynamic changes of the crustal deformation in the Tibetan Plateau before earthquakes of magnitude above 8 since GPS observation, 26 GPS continuous sites from CMONOC I established by the Chinese Government were selected to study the site displacement and baseline. Velocity field data processing adopted a unified strategy to process or reprocess all raw GPS data using GAMIT [24]. Besides, approximately 100 globally distributed International GPS Service (IGS) sites are incorporated in data processing. Finally, the results of the same day are combined using GLOBK [25].

The GPS continuous data used in this paper were adopted using GAMIT/GLOBK and QOCA software [24,25,26,27], as well as the site coordinate time series results under the ITRF2014 reference framework processed by combining with the IGS sites around the Chinese Mainland [28]. For the parameter setting of the reference stations, loosely constrained algorithms are used for the station positions in the processing of GAMIT software; the value of the constraint is 50 mm. In the adjustment with GLOBK software, 75 reference stations from two global networks are selected to estimate the GPS coordinates and covariance with respect to the IGS14 frame [28,29]. The station motion model is established. To reduce the interference of periodic signals in the identification of time series anomaly information, the single-site time series are all the results after removing the period, correcting the earthquake events, replacing the instruments, and the step caused by unknown reasons. The methods are in accordance with the literature [28,30].

2.3 Methodology

GPS velocity field can reflect the crustal relative movement and deformation directly and clearly [31,32]. When the crustal movement and deformation are analyzed in a region, the rigid motion without deformation will be removed in the study area generally. The South China block is the largest relatively stable block on the Chinese Mainland. When the movement and deformation of the Tibetan Plateau are used by GPS data, taking the stable South China block as the reference datum is the best choice [33,34,35,36].

Time series analysis of GPS sites can study the position time series of continuous GPS sites. Due to this, the microdynamic change information of crustal movement can be extracted and the accurate movement trend of sites can be obtained [21,37,38,39]. The baseline time series between GPS stations refers to using the data of two GPS continuous stations to calculate the time-varying curve of the distance between them. The time series of the baselines can directly reflect the shortening or elongation of the crust between two GPS stations and can be used to describe the extension, compression, and strike-slip characteristics of the blocks spanned [22,40,41,42,43,44,45,46].

3 The characteristics of regional horizontal movement in the Tibetan Plateau

This article discussed the results of the GPS velocity field relative to the South China block based on the velocity field results of the ITRF2008 reference framework [15]. Figure 2 shows the GPS velocity field from 1999 to 2016 and the block boundary of the Tibetan Plateau and its surrounding areas. From the velocity field results relative to the South China block, it can be seen that because the south of the Tibetan Plateau is affected by NE-ward compression of the Indian Plate and resisted by three rigid blocks of Alashan, Ordos, and South China in the northeast and east, and also of the Tarim Basin in the northwest, Figure 2 shows that the size and direction of the velocity field of Tibetan Plateau are significantly different from the southwest to the north, northeast, east, and southeast.

Figure 2 
               GPS velocity field for the investigated area spanning the 1999–2016 interval.
Figure 2

GPS velocity field for the investigated area spanning the 1999–2016 interval.

Figure 2 shows that the velocity values of the Tibetan Plateau ranging from 36–40 mm/a in the southwest Tibetan Plateau to 6–8 mm/a in the southeast border of the Tarim basin and 11–25 mm/a in the southwest border of the Tarim basin, to 3–7 mm/a in the southeast border of the Alashan block, to 0.2–2 mm/a in the west boundary of the Ordos block, to 0.3–1.9 mm/a in the west boundary of South China. The characteristic of velocity field size is that it decreases gradually from southwest to north, northeast, and southeast in the study area, and the velocity value in the west is significantly greater than that in the east in the same dimension. Compared with the northwest boundary of the Tibetan Plateau, the velocity values are quite different on both sides of the northeast boundary, the east boundary, and the southeast boundary of the Tibetan Plateau. Because the region with large far-field velocity and small near-field velocity is easy to accumulate strain energy, this kind of region is easy to generate an earthquake.

On the other hand, the motion vector direction of the Tibetan Plateau changes from the near NS-ward in the southwest to NW-ward in the northwest relative to the South China block. East of the E90° takes NE-ward, and the movement tends to be a bifurcate horizontal expansion in the east margin of the Tibetan Plateau, which shows as the anticlockwise rotating Gan Qing vortex and clockwise rotating Yunnan–Tibet vortex [47,48]. This kind of rotation is more obvious and intense in the south than that in the north. The direction of the motion vector is from west to east and from south to north in the northern region (E90–E104°, N34–40°), which shows NEE-NE-NNE-NS-NNW-NW-ward in general, showing an anticlockwise rotating tendency. The direction of the motion vector is from west to east and from north to south in the southern region (E90–E104°, N21–34°), which shows NNE-NE-EW-SEE-SE-S-SW-SWW-ward in general, showing a clockwise rotating tendency. According to previous studies, the mutation region, where the velocity and direction of the crustal movement are quite different in the eastern and northeastern margins of the Tibetan Plateau, was an abnormal region of crustal deformation, and it was also a region with stress accumulation and the earthquake might happen here [47].

4 The North-ward displacement time series of GPS continuous sites

As shown in Figure 3, the GPS continuous sites in the interior and boundary of the Tibetan block are LHAS, KMIN, XIAG, LUZH, DLHA, DXIN, XNIN, XIAA, and YANC. Due to the serious interference of XNIN and XIAA, they were not analyzed in this paper. Figure 4 shows the results of the north-ward displacement time series of the remaining seven sites. When tracking and analyzing the time series of GPS coordinates of a single site, the time series graphs of Figure 4 all removed the linearity to highlight the dynamic information. Table 2 shows the movement rates of these seven sites in different periods.

Figure 3 
               GPS continuous station and baselines from the Indian Plate to Tibetan Plateau.
Figure 3

GPS continuous station and baselines from the Indian Plate to Tibetan Plateau.

Figure 4 
               The north-ward displacement detrend time series of GPS continuous site: (a) LHAS, (b) KMIN, (c) XIAG, (d) LUZH, (e) DLHA, (f) DXIN, and (g) YANC.
Figure 4

The north-ward displacement detrend time series of GPS continuous site: (a) LHAS, (b) KMIN, (c) XIAG, (d) LUZH, (e) DLHA, (f) DXIN, and (g) YANC.

Table 2

The north-ward displacement rates of GPS continuous sites in Figure 4

Station Rate (mm/a)
Jan. 1, 2002–Dec. 30, 2005 Dec. 31, 2005–May 12, 2008 May 13, 2008–Jan. 1, 2013 Jan. 2, 2013–Apr. 25, 2015 Jan. 1, 2017–May. 19, 2021
LHAS −0.10 ± 0.04 1.54 ± 0.09 0.67 ± 0.05 −3.30 ± 0.01 −1.17 ± 0.05
KMIN −3.99 ± 0.07 −0.80 ± 0.07 4.26 ± 0.05 −2.56 ± 0.09 1.58 ± 0.04
XIAG −4.19 ± 0.07 3.43 ± 0.04 0.40 ± 0.05 −0.43 ± 0.02 0.63 ± 0.06
LUZH −1.74 ± 0.05 −0.81 ± 0.08 0.79 ± 0.03 0.17 ± 0.07 0.60 ± 0.04
DLHA −0.87 ± 0.04 −0.13 ± 0.09 0.55 ± 0.03 1.45 ± 0.08 1.63 ± 0.04
DXIN −0.53 ± 0.03 0.19 ± 0.07 −0.73 ± 0.04 0.07 ± 0.01 1.11 ± 0.03
YANC −1.35 ± 0.03 −0.84 ± 0.05 0.03 ± 0.03 −0.36 ± 0.06 0.78 ± 0.03

Among the seven GPS sites mentioned above, the site LHAS is located at the south Tibetan block, KMIN, XIAG, and LUZH are located at the southeastern boundary of the Tibetan block, and DLHA, YANC, and DXIN are located at the northeast Tibetan Plateau. The movement trend of the three regional sites is similar: there were turning changes before both the 2008 Wenchuan M S 8.0 earthquake and the 2015 Nepal M S 8.1 earthquake. Since the observation start time of CMONOC phase Ⅰ is close to the 2001 West Kunlun Mountains M S 8.1 earthquake, the trend before the earthquake could not be judged well; it was not analyzed in this article. In Figure 4(a)–(d), the N-ward time series curves of LHAS, KMIN, XIAG, and LUZH tend to go down from 2002 to 2005, and the downward trend changes from 2005 to May 12, 2008 (occurred the Wenchuan earthquake). Before the Wenchuan earthquake, in which the trend of LHAS and XIAG changes from downward to upward, the movement rates before the turning are −0.10 ± 0.04 and −4.19 ± 0.07 mm/a and after the turning are 1.54 ± 0.09 and 3.43 ± 0.04 mm/a, respectively. Also, the downward trend of KMIN and LUZH obviously slows down; the movement rates before the slowing down are −3.99 ± 0.07 and −1.74 ± 0.05 mm/a and after the slowing down are −0.80 ± 0.07 and −0.81 ± 0.08 mm/a, respectively. The changing range of LUZH rates near the Wenchuan epicenter is relatively small among the four sites, and the changing range of KMIN, LHAS, and XIAG, which are far from the Wenchuan epicenter, are all relatively large. In Figure 4(e)–(g), the north-ward time series curves of DLHA, DXIN, and YANC tend to go down from 2002 to 2005, and the downward trend changes from 2005 to May 12, 2008, before the Wenchuan earthquake, showing that the downward trend tends to be slow down, the movement rates before the slowing down are −0.87 ± 0.04, −0.53 ± 0.03, and −1.35 ± 0.03 mm/a, and after the slowing down are −0.13 ± 0.09, 0.19 ± 0.07, and −0.84 ± 0.05 mm/a, respectively. It can be seen from the velocity field diagram in Figure 2 that the epicenter of the Wenchuan earthquake is located in the northeast corner of the outermost circle of the clockwise rotating Yunnan Tibet vortex from west to east in the south of the Tibetan Plateau, and LHAS, KMIN, XIAG, and LUZH are all located in the clockwise rotating Yunnan Tibet vortex. From the perspective of power sources, they were more related to the seismogenic of the Wenchuan earthquake. Therefore, the abnormal changes are obviously larger than those sites of DLHA, YANC, and DXIN in the northeast.

As can be seen in Figure 4(b)–(d), the north-ward time series curves of KMIN, XIAG, and LUZH tended to rise after the 2008 Wenchuan earthquake and fall down about two years before the Nepal earthquake on April 24, 2015, and the movement rates changed from 4.26 ± 0.05, 0.40 ± 0.05, and 0.79 ± 0.03 to −2.56 ± 0.09, −0.43 ± 0.02, and 0.17 ± 0.07 mm/a, respectively, indicating that the north-direction motion slowed down before the Nepal earthquake. The changes in site LHAS were unable to be analyzed due to the data missing in Figure 4(a) before the Nepal earthquake. As can be seen in Figure 4(e)–(g), the north-ward time series curves of DLHA, DXIN, and YANC showed no trend turning change before the Nepal earthquake. Therefore, the abnormal range of the Nepal earthquake before the earthquake was still in the southwest margin of the Tibetan Plateau, and there were not enough sites in the south margin to determine the abnormal condition before the earthquake. It can be seen from Figure 4 that the GPS continuous sites KMIN, XIAG, LUZH, DLHA, DXIN, and YANC in the interior and boundary of the Tibetan block are turning upward in the time curves in the last five years at least, and their movement rates are 1.58 ± 0.04, 0.63 ± 0.06, 0.60 ± 0.04, 1.63 ± 0.04, 1.11 ± 0.03, and 0.78 ± 0.03 mm/a, showing that the north-ward movement is significantly enhanced, which needs to be focused on.

5 The GPS baseline time series

The calculation result of a single station has a systemic influence such as a reference frame [49,50,51]; thus, to reduce this influence, the characteristics of crustal deformation of the Tibetan Plateau for nearly 20 years are analyzed in this article through the baseline time series method by selecting GPS continuous observation stations on the Indian Plateau and different blocks of the Tibetan Plateau. The research results have certain reference significance for the studies and evaluations of the strong earthquake activities in the Tibetan Plateau and surrounding areas.

Figure 3 shows the seven baselines from the GPS station named IISC in the Indian Plate to seven GPS stations as TASH, WUSH, DXIN, DLHA, LHAS, YANC, and XIAG in or around the Tibetan Plateau (where the stations of WUSH, DXIN, and YANC are located in the stable blocks; however, others are located in the active blocks). The crustal deformation characteristics of the Tibetan Plateau under the pushing effect of the Indian Plate since 1999 are analyzed via the detrended baselines time series, which are shown in Figures 5 and 6. Table 2 shows the statistics of the baseline rate, the annual variation ratio (AVR), and the changes of turning in Figures 5 and 6. The baseline is in a shortened state while the rate value in the table is negative.

Figure 5 
               The GPS detrend baseline time series: (a) IISC-WUSH, (b) IISC-DXIN, and (c) IISC-YANC.
Figure 5

The GPS detrend baseline time series: (a) IISC-WUSH, (b) IISC-DXIN, and (c) IISC-YANC.

Figure 6 
               The GPS detrend baseline time series: (a) IISC-DLHA, (b) IISC-LHAS, (c) IISC-TASH, and (d) IICS-XIAG.
Figure 6

The GPS detrend baseline time series: (a) IISC-DLHA, (b) IISC-LHAS, (c) IISC-TASH, and (d) IICS-XIAG.

In Figure 5, baselines IISC-WUSH, IISC-DXIN, and IISC-YANC are the three baselines from Bangalore, India, to the western boundary of the Tarim Basin, the Alexa block, and the Ordos block, respectively. Table 2 shows that their shortening rates are −20.24 ± 0.01, −34.70 ± 0.01, and −33.21 ± 0.01 mm/a, and the AVR is −18.05 × 10−7/a, −28.85 × 10−7/a, and −27.99 × 10−7/a, respectively; the rate values are all negative, indicating that the baselines are all shortened. Figure 5(a)–(c) show that the baseline shortening rate has changed a little in recent 20 years, indicating that the pushing rate of the Indian Plate to the north of the Tibetan Plateau is relatively stable, which may be related to the location of the GPS station in the three baselines, the movement of the station WUSH is blocked by the Junggar Basin, and DXIN and YANC are located in rigid blocks.

Table 2 shows that the shortening rates of baseline IISC-DLHA, IISC-LHAS, IISC-TASH, and IISC-XIAG are −27.29 ± 0.01, −11.22 ± 0.01, −28.14 ± 0.01, and −35.30 ± 0.01 mm/a, respectively, and the AVR are −13.83 × 10−7/a and −6.60 × 10−7/a, −14.20 × 10−7/a and −13.87 × 10−7/a respectively, and the rate values are all negative, indicating that the baselines are all shortened.

As seen from from Figure 6(a) and (b), the time series curves of the baseline IISC-DLHA and IISC-LHAS fall rapidly from 1999 to 2005. From 2005 to the Wenchuan earthquake, the downward trend of baseline IISC-DLHA turns to be slow and the trend of baseline IISC-LHAS turns to go up; they all reflect that the crustal shortening rate slowed down about three years before the Wenchuan earthquake. The slowdown rate of the baseline IISC-LHAS is more obvious, which may be related to the strain accumulation to a certain extent before the Wenchuan earthquake.

It can be seen from Figure 6(a)–(d) and Table 2 the four baselines IISC-DLHA, IISC-LHAS, IISC-TASH, and IISC-XIAG changed after the 2015 Nepal earthquake. The upward trend of baseline IISC-DLHA slows down, and the trend of baseline IISC-LHAS, IISC-TASH, and IISC-XIAG turns to go down; they all indicate that the crustal shortening rate accelerated after the 2015 Nepal earthquake, which may reflect the adjustment of the Indian Plate compression movement after the Nepal earthquake.

The last column of data in Table 2 is the absolute value of the shortening rate in the two periods of the seven baselines in Figure 3: the absolute value of the rate after the 2015 Nepal M S 8.1 earthquake minus the absolute value of the rate between the 2008 Wenchuan M S 8.0 earthquake and the 2015 Nepal M S 8.1 earthquake. Figure 5 and Table 2 both show that the first three baselines IISC-WUSH, IISC-DXIN, and IISC-YANC have relatively small changes, with absolute differences of 0.94, 0.12, and 0.18, respectively. Figure 6 and Table 2 show that baselines IISC-DLHA, IISC-TASH, and IISC-XIAG have significant turning changes before and after the Nepal earthquake, and the absolute difference has also significantly increases, with values of 1.02, 2.54, and 1.64. Baseline IISC-LHAS has been in the state of shortening acceleration before and after the Nepal earthquake, so its absolute difference is small, at 0.24. The site LHAS is closer to the epicenter of Nepal than other sites, and the rate change of the baseline IISC-LHAS is the smallest after the earthquake, which may indicate that the area near the site LHAS is in a significant invariant zone in the dynamic adjustment process caused by the large earthquake, and where a certain degree of strain has accumulated [22].

It can be seen from Figure 3 that baselines IISC-DLHA and IISC-LHAS are in the same northeast-ward. DLHA and LHAS are located in the northeast and south margins of the Tibetan Plateau, respectively. The annual motion rate of baseline IISC-DLHA is significantly larger than that of baseline IISC-LHAS, indicating that there is significant deformation in the interior of the Tibetan Plateau, and it promotes the strong earthquakes to be bred and occurred in Bayan Har and its surrounding blocks.

6 Discussion

The energy of earthquake preparation comes from the accumulation of strain energy, which is generated in the process of crustal movement under the dynamic action of the tectonic boundary [52]. The generation and preparation of earthquakes lead to crustal deformation directly. Therefore, earthquakes are the result of long-term accumulation and sudden release of energy in the process of crustal movement [53]. GPS observation can provide relative motion information of various spatial scales, so these observation results are widely used to study the relationship between the temporal and spatial distribution dynamics of crustal movement and strong earthquakes, the deformation process and mechanism of earthquakes, etc., which can provide a reference for the judgment and exploration of danger sites of a strong earthquake. It may become an important aspect of the expansion of earthquake prediction research to physical prediction.

6.1 Significant difference in the area of crustal movement

Strong earthquakes usually occur in regions with significant differences in crustal movement, that is, the regions with significant differences in vector size and direction of motion velocity. Only the difference in crustal movement can produce strain accumulation. When the regional crustal movement is consistent, it will not produce strain accumulation, and its strong earthquake risk is not high [54]. According to the results of recent research, almost 100% of strong earthquakes of magnitude 8 and 80% of strong earthquakes of magnitude above 7 are located in the boundary of the block in the Chinese Mainland (Figure 1) [55]. The active block can produce movement and deformation relatively under the dynamic action caused by compression between plates and intraplate mantle convection. The velocity and direction of crustal movement in the eastern and northeastern margins of the Tibetan Plateau are quite different, which is an abnormal region of crustal deformation and a frequent occurrence of strong earthquakes. It can be seen from Figure 2 that the seismic source of the Wenchuan earthquake is located in the eastern boundary of the Tibetan Plateau, which is the lowest velocity region of the whole Tibetan Plateau.

6.2 Seismogenic process of near-field and far-field and post-earthquake influence

Four to five years before the 2008 Wenchuan earthquake, the seven GPS continuous stations in the interior and boundary of the Tibetan Plateau block all have abnormal changes. Compared with the northeastern margin of the Tibetan Plateau [Figure 4(e)–(g)], the abnormal variation range of the sites in the southern and southwestern margins was obviously larger [Figure 4(a)–(d)]. The sites in the southern Tibetan Plateau experienced an abnormal change and tended to go down during the two years before the 2015 Nepal earthquake [Figure 4(b)–(d)]. It reflects that the stress accumulation was fast and the required time of accumulation was relatively short in the region of the seismic source area close to the seismic source that the power source needed for earthquake preparation. However, when the power source was far away from the seismic source area, the stress accumulation was slow and took relatively a long time.

The north-ward movement of the Tibetan Plateau block became significantly stronger after the 2015 Nepal M S 8.1 earthquake (Figure 4 and Table 2). The earthquake may have resulted in the unlocking of large-scale faults near the seismic source area, and the plate subduction and collision would be further intensified. Therefore, over the long term, the northeast-ward subduction and collision of the Indian Plate show an accelerating trend, which promotes the strain accumulation inside of the Tibetan Plateau, and attention needs to be paid continuously to the strong earthquake risk in this region.

6.3 Influence of deformation difference

The crustal deformation state and strong earthquake activity of the Chinese Mainland are closely related to the dynamic process of plate boundary. Because the GPS observation has a unified reference framework, it, which crosses the plate boundary, can effectively identify the dynamic process of the above boundary. From 1999 to 2019, Ma thought that the AVR of baseline DXIN-LHAS (located between site DXIN in the Alashan block and site LHAS in the Lhasa block) was 5.6 times that of the baseline DXIN-DLHA (located between site DXIN in the Alashan block and site DLHA in the Qaidam block) and the AVR of the baseline YANC-LHAS (located between site YANC in Ordos block and site LHAS in Lhasa block) was 7.4 times that of the baseline YANC-DLHA (located between site YANC in Ordos block and site DLHA in the Qaidam block) [56]. Due to the above evidence, they showed that the AVR between the Lhasa block and two rigid blocks of Alashan and Ordos was significantly faster than that between the Qaidam block and two rigid blocks of Alashan and Ordos, indicating that the internal deformation of the Qiangtang block and the Bayan Har block between LHAS-DLHA was strong relatively. Combined with the time series results of site and baseline in this paper (Figures 4 and 5 and 6, and Table 3), it can be found that the characteristics of shortening enhanced features between the Indian Plate and the Tibetan Plateau block have been obvious since 2015, which may reflect that the northeast movement of the Tibetan Plateau continues to strengthen under the compression of the Indian Plate.

Table 3

The rates and AVR of GPS baselines in Figure 2

Baseline Rate (mm/a) AVR (10−7/a) Rate (mm/a) Rate (mm/a) Absolute difference
1999–2020 1999–2020 ①–②
IISC-WUSH −20.24 ± 0.01 −18.05 −20.833 ± 0.03 −19.896 ± 0.03 0.94
IISC-DXIN −34.70 ± 0.01 −28.85 −34.397 ± 0.02 −34.514 ± 0.03 0.12
IISC-YANC −33.21 ± 0.01 −27.99 −32.897 ± 0.02 −33.076 ± 0.02 0.18
IISC-DLHA −27.29 ± 0.01 −13.83 −26.213 ± 0.02 −27.235 ± 0.02 1.02
IISC-LHAS −11.22 ± 0.01 −6.60 −11.930 ± 0.03 −12.173 ± 0.03 0.24
IISC-TASH −28.14 ± 0.01 −14.20 −27.532 ± 0.02 −30.069 ± 0.07 2.54
IISC-XIAG −35.30 ± 0.01 −13.87 −34.300 ± 0.05 −35.943 ± 0.07 1.64

AVR, the annual variation ratio. ① From 2008 Wenchuan M S 8.0 Earthquake to 2015 Nepal M S 8.1 earthquake, and ② after 2015 Nepal M S 8.1 earthquake.

6.4 Identification of the uncoordinated regions in the process of large earthquake adjustment

The 2008 Wenchuan M S 8.0 earthquake and the 2013 Lushan M S 7.0 earthquake occurred in the significant invariant zone (strain energy accumulation area) in the dynamic adjustment process caused by the 2001 Kunlun Mountain M S 8.1 earthquake and the 2004 Sumatra M S 9.0 earthquake [22]. Although there may be a phenomenon that the Indian Plate’s pushing effect on the Tibetan plateau after the 2015 Nepal M S 8.1 earthquake is enhanced from a large-scale perspective (Figure 4), which also reflects the internal deformation enhancement of the Tibetan Plateau. We need to use other data and methods to find the invariant zone in adjustment or the area with high locking characteristics from a small-scale perspective, to determine the specific risk area.

7 Conclusions

The results of the velocity field show that the velocity field of the Tibetan Plateau decreases from southwest to north, northeast, and southeast, and the velocity value in the west is significantly greater than that in the east in the same dimension. The maximum value is about 40 mm/a in the southwest margin, and the minimum value is about 0.2 mm/a in the east margin. The direction of the motion vector changes from NS-ward in the southeast to NW-ward in the northwest, the east of E90° turns into the NE direction, and the movement tends to be a bifurcate horizontal expansion in the east margin of the Tibetan Plateau, which shows anticlockwise rotating of Gan Qing vortex and clockwise rotating of Yunnan–Tibet vortex (Figure 2). The eastern and northeastern margins of the Tibetan Plateau are the mutation regions with great differences in the velocity and direction of crustal movement, and they are the areas of stress accumulation and the areas where earthquakes may occur.

The results of the north-ward displacement detrend time series (Figure 4 and Table 2) of GPS continuous sites show that the 7 GPS continuous sites of LHAS, KMIN, XIAG, LUZH, DLHA, DXIN, and YANC in the interior and boundary of the Tibetan Plateau block all had abnormal changes before the 2008 Wenchuan earthquake, showing a downward trend from 2002 to 2005 and the downward tendency became slow or changed upward from 2005 to the occurrence of the earthquake. In addition, the variation range of the sites in the southern part of the Tibetan Plateau was obviously larger than that in the other areas, and the southern sites of the Tibetan Plateau tendency changed abnormally and decreased in the two years before the Nepal earthquake. The continuous GPS sites of KMIN, XIAG, DLHA, LUZH, YANC, and DXIN in the interior and boundary of the Tibetan block have shown a significant increase in the north-ward movement in at least five years.

The results of the site’s baseline detrend time series (Figures 5 and 6, Table 3) show that the movement rate of the Indian Plate relative to the rigid block on the edge of the Tibetan Plateau was relatively stable. Under the influence of the 2008 Wenchuan M S 8.0 earthquake, the movement of the Indian Plate relative to the Tibetan Plateau had been changed from 2008 to 2015, showing as the compression slowed down. The compression of the India plate against the NE-ward of the Tibetan Plateau was significantly stronger than it was 7 years ago after 2015, and the enhancement promoted the preparation and occurrence of strong earthquakes in west China.

Acknowledgments

The authors are grateful to the First Monitoring and Application Center, China Earthquake Administration for providing the time series data of GPS continuous stations, to Dr. Yanqiang Wu, Dr. Hongbao Liang, and Dr. Wei Zhan for the support in the calculating process of the time series, and to Mr. Xiaowei Xu for his contribution to language editing. The authors also appreciate the editors and reviewers for their detailed and valuable comments and suggestions.

  1. Funding information: This work was funded by Spark Program of Earthquake Sciences (XH21035Y), the Special Fund for Innovation Team, Gansu Earthquake Agency (Grant No. 2020TD-01-01), the Science and Technology Project of Gansu Province (20JR10RA502), Spark Program of Earthquake Sciences (XH22039D, XH20057), Grant of NSFC (51778590, 41904092), and the Earthquake Tracking Task of CEA (2021010213, 2020010211, 2018010204).

  2. Author Contributions: Conceptualization: Haiping Ma, and Qian Wang. Methodology” Haiping Ma and Jing Zhao. Analysis” Haiping Ma, Jing Zhao, Hui Zhang and Minjuan Li. Investigation: Haiping Ma, Hui Zhang, Qian Wang, Jing Zhao and Shanyi Wu. Data processing: Haiping Ma, Hui Zhang, Pengtao Wang and Jing Zhao. Writing - original draft preparation: Haiping Ma, Hui Zhang, Jing Zhao, Qian Wang, Minjuan Li, Pengtao Wang and Shanyi Wu. Writing - revised draft preparation: Haiping Ma, Qian Wang, Jing Zhao and Zhigiang Ma. Project administration: Haiping Ma and Qian Wang. Funding acquisition: Haiping Ma and Qian Wang.

  3. Conflicts of interest: The authors declare that they have no competing interest.

  4. Data availability statement: All data included in this study are available upon request by contact with the corresponding author.

References

[1] Royden LH, Burchfiel BC, Hilst RD. The geological evolution of the tibetan plateau. Science. 2008;321:1054–9.10.1126/science.1155371Search in Google Scholar PubMed

[2] England P, Molnar P. Active deformation of asia: from kinematics to dynamics. Science. 1997;278(5338):647–50.10.1126/science.278.5338.647Search in Google Scholar

[3] Chen ZH, Burchfiel BC, Liu Y, King RW, Royden LH, Tang W, et al. Global positioning system measurements from eastern tibet and their implications for India/Eurasia intercontinental deformation. J Geophys Res Solid Earth. 2000;105(7):16215–27.10.1029/2000JB900092Search in Google Scholar

[4] Tapponnier P, Peltzer G, Armijo R. On the mechanics of the collision between India and Asia. Geological Society. Vol. 19. London, UK: Special Publications; 1986. p. 113–57.10.1144/GSL.SP.1986.019.01.07Search in Google Scholar

[5] Shen Z, Lü J, Wang M, Bürgmann R. Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau. J Geophys Res. 2005;110(B11409):1–17. 10.1029/2004JB003421.Search in Google Scholar

[6] Zuza AV, Cheng X, Yin A. Testing models of Tibetan Plateau formation with Cenozoic shortening estimates across the Qilian Shan-Nan Shan thrust belt. Geosphere. 2016;12(2):501–32.10.1130/GES01254.1Search in Google Scholar

[7] Copley A. Kinematics and dynamics of the southeastern margin of the Tibetan Plateau. Geophys J Int. 2008;174(3):1081–100.10.1111/j.1365-246X.2008.03853.xSearch in Google Scholar

[8] Deng Q, Zhang P, Ran Y, Yang X, Min W, Chu Q. Basic characteristics of active tectonics of China. Sci China Ser D Earth Sci. 2003;46:356–72.10.1360/03yd9032Search in Google Scholar

[9] Zhang P, Deng Q, Zhang G. Active tectonic blocks and strong earthquakes in the continent of China. Sci China Earth Sci. 2003;z2:13–24.10.1360/03dz0002Search in Google Scholar

[10] Tapponnier P, Zhiqin X, Roger F, Meyer B, Arnaud N, Wittlinger G, et al. Oblique stepwise rise and growth of the Tibet Plateau. Science. 2001;294:1671–7.10.1126/science.105978Search in Google Scholar PubMed

[11] Zhang P, Shen Z, Wang M, Gan W, Bürgmann R, Molnar P, et al. Continuous deformation of the Tibetan Plateau from global positioning system data. Geology. 2004;32(9):809–12.10.1130/G20554.1Search in Google Scholar

[12] Zhang F, Wu Y, Sun D. Analysis of the baseline time series from continuous gps stations in the Mainland of China. Earthquake. 2011;31(4):86–97.Search in Google Scholar

[13] Briole P, Ganas A, Elias P, Dimitrov D. The GPS velocity field of the Aegean. new observations, contribution of the earthquakes, crustal blocks model. Geophys J Int. 2021;226(1):468–92.10.1093/gji/ggab089Search in Google Scholar

[14] Gan W, Zhang P, Shen Z, Niu Z, Wang M, Wan Y, et al. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements. J Geophys Res Solid Earth. 2007;112:B08416.10.1029/2005JB004120Search in Google Scholar

[15] Wang M, Shen Z. Present‐day crustal deformation of continental china derived from GPS and its tectonic implications. J Geophys Res Solid Earth. 2020;125(2).10.1029/2019JB018774Search in Google Scholar

[16] Wang W, Qiao X, Yang S, Wang D. Present-day velocity field and block kinematics of Tibetan Plateau from GPS measurements. Geophys J Int. 2017;208(2):1088–102.10.1093/gji/ggw445Search in Google Scholar

[17] Liang S, Gan W, Shen C, Xiao G, Liu J, Chen W, et al. Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements. J Geophys Res: Solid Earth. 2013;118(10):5722–32.10.1002/2013JB010503Search in Google Scholar

[18] Wu Y, Jiang Z, Zhao J, Liu X, Wei W, Liu Q, et al. Crustal deformation before the 2008 Wenchuan MS 8.0 earthquake studied using GPS data. J Geodynamics. 2015;85:11–23.10.1016/j.jog.2014.12.002Search in Google Scholar

[19] Liu X, Wu Y, Jiang Z, Zhan W, Li Q, Wei W, et al. Preseismic deformation in the seismogenic zone of the Lushan MS7.0 earthquake detected by GPS observations, Fixation of sulfur dioxide into small molecules. Sci China Earth Sci. 2015;58:1592–601.Search in Google Scholar

[20] Zhan W, Wu Y, Liang H, Zhu S, Zhang F, Liu J. Characteristics of the seismogenic model for the 2015 Nepal MW 7.8 earthquake derived from GPS data. Chin J Geophysics. 2015;58(5):1818–26.Search in Google Scholar

[21] Department of monitoring and prediction, China seismological bureau. Research on short-term prediction method of GPS, satellite remote sensing and earth changing magnetic field earthquake. Beijing: Seismological Press. 2006.Search in Google Scholar

[22] Wu Y, Jiang Z, Liang H, Pang Y, Zhu S, Chang L, et al. 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 Sens. 2018;10(6):894.10.3390/rs10060894Search in Google Scholar

[23] Ma H, Feng J, Wu S, Xu R, Huang L, Li N, et al. The analysis of the abnormal characteristics of gps data before the 2013 minxian-zhangxian MS 6.6 earthquake. Earth Environ Sci. 2008;153:062044.10.1088/1755-1315/153/6/062044Search in Google Scholar

[24] Herring TA, King RW, McClusky SC. GAMIT reference manual, GPS analysis at MIT. Release 10.4. Cambridge, MA, USA: Massachussetts Institute Technology; 2010a. Available online: edu/∼simon/gtgk/index.htm (accessed on October 5 2011).Search in Google Scholar

[25] Herring TA, King RW, McClusky SC. GLOBK reference manual,” global kalman filter VLBI and GPS analysis program. Release 10.4. Cambridge, MA, USA: Massachussetts Institute Technology; 2010b. Available online: http://www-gpsg.mit.edu/∼simon/gtgk/index.htm (accessed on October 5 2018).Search in Google Scholar

[26] Dong D, Hering TA, King RW. Estimating regional deformation from a combination of space and terrestrial geodetic data. J Geophys Res. 1998;72:200–14.10.1007/s001900050161Search in Google Scholar

[27] Zhao J, Zhan W, Ren J, Jiang Z, Gu T, Liu J, Niu A, Yuan Z, et al. GPS time series inversion of the healing process of the middle segment of the longmenshan fault after the 2008 wenchuan earthquake. Acta Geodaetica et Cartographica Sin. 2021;50(1):37–51.Search in Google Scholar

[28] Liang H, Zhan W, Li J. Vertical surface displacement of mainland China from GPS using the multisurface function method. Adv Space Res. 2021;68(12):4898–915.10.1016/j.asr.2021.02.024Search in Google Scholar

[29] Altamimi AB, Ramadan RA. Towards internet of things modeling: a gateway approach. Complex Adapt Syst Modeling. 2016;4(1):25.10.1186/s40294-016-0038-3Search in Google Scholar

[30] Nikolaidis R. Observation of geodetic and seismic deformation with the global positioning system. Doctoral dissertation. University of California, San Diego; 2002.Search in Google Scholar

[31] Holt WE, Chamot‐Rooke N, Le Pichon X, Haines AJ, Shen‐Tu B, Ren J. Velocity field in Asia inferred from quaternary fault slip rates and global positioning system observations. J Geophys Res Solid Earth. 2000;105(B8):19185–209.10.1029/2000JB900045Search in Google Scholar

[32] Özarpacı S, Doğan U, Ergintav S, Çakır Z, Özdemir A, Floyd M, et al. Present GPS velocity field along 1999 Izmit rupture zone: evidence for continuing afterslip 20 years after the earthquake. Geophys J Int. 2021;224(3):2016–27.10.1093/gji/ggaa560Search in Google Scholar

[33] Wu Y, Jiang Z, Liu X, Wei W, Zhu S, Zhang L, et al. A comprehensive study of gridding methods for GPS horizontal velocity fields. Pure Appl Geophys. 2017;174(3):1201–17.10.1007/s00024-016-1456-zSearch in Google Scholar

[34] Jiang Z, Ma Z, Zhang X, Chen B, Hu Y, Wang S, et al. Analysis of recent horizontal crustal strain field and tectonic deformation in the northeast margin of Qingahai Tibet block. Seismology Geol. 2001;23(003):337–45.Search in Google Scholar

[35] Ding X, Zhan W, Wang Z. Influence on results calculated for regional GPS continuous data with different reference frames. J Geodesy Geodyn. 2014;34(2):161–8.Search in Google Scholar

[36] Zhao J, Ren J, Liu J, Jiang Z, Liu X, Liang H, et al. Coupling fraction and relocking process of the Longmenshan Fault Zone following the 2008 Mw7.9 Wenchuan earthquake – ScienceDirect. J Geodyn. 2020;137:101730.10.1016/j.jog.2020.101730Search in Google Scholar

[37] Guns K, Xu X, Bock Y, Sandwell D. GNSS-corrected InSAR displacement time series spanning the 2019 Ridgecrest, CA earthquakes. Geophys J Int. 2022;230(2):1358–73.10.1093/gji/ggac121Search in Google Scholar

[38] Kudłacik I, Kapłon J, Lizurek G, Crespi M, Kurpiński G. High-rate GPS positioning for tracing anthropogenic seismic activity: the 29 January 2019 mining tremor in Legnica- Gogów Copper District, Poland. Measurement. 2020;168(15):108396.10.1016/j.measurement.2020.108396Search in Google Scholar

[39] Hall K, Schmidt D, Houston H. Peak tremor rates lead peak slip rates during propagation of two large slow earthquakes in Cascadia. Geochem Geophys Geosyst. 2019;20(2):4665–75.10.1029/2019GC008510Search in Google Scholar

[40] Jiang Z, Fang Y, Wu Y, Wang M, Du F, Ping J. The dynamic process of regional crustal movement and deformation before Wenchuan MS8.0 earthquake. Chin J Geophys. 2009;52(2):505–18.Search in Google Scholar

[41] Amiri-Simkooei R, Tiberius C. Assessing receiver noise using GPS short baseline time series. GPS Solut. 2007;11(1):21–35.10.1007/s10291-006-0026-8Search in Google Scholar

[42] Iz HB. Differencing reveals hidden changes in baseline length time-series. J Geodesy. 2006;80(5):259–69.10.1007/s00190-006-0066-4Search in Google Scholar

[43] Bergstrand S, Scherneck HG, Milne GA, Johansson JM. Upper mantle viscosity from continuous GPS baselines in Fennoscandia. J Geodynamics. 2005;39(2):91–109.10.1016/j.jog.2004.08.004Search in Google Scholar

[44] Guo L, Hu X, Zhang F, Chen C, Sun D. Baseline variation of GPS continuous sites and current activity of Qinghai-Tibet block. J Geodesy Geodynamics. 2009;29(4):10–4.Search in Google Scholar

[45] Ma H, Zhang H, Feng J, Wu S, Wang P, Dou X, et al. Characteristics of the stress-strain variation before the 2019 Xiahe MS 5.7 earthquake. Prog Geophysics(Chin). 2021;36(3):933–9.Search in Google Scholar

[46] Zhu S, Zhou W. Regional Crustal Deformation Analysis before and after the Minxian-Zhangxian Earthquake. China Earthq Eng J. 2015;37(3):731–8.Search in Google Scholar

[47] Chen Z, Liu Y, Zhang X. GPS monitoring, crustal deformation and geodynamics in the East Qinghai-Tibet Plateau. Guangzhou: Guangdong Sci Technol Press. 2015.Search in Google Scholar

[48] Ma H, Wu Y, Xu R, Wu S, Wang Q. Research of recent GPS crustal deformation characteristics in the Northeastern Edge of Qinghai-Tibet Plateau. J Phys Conf Ser. 2017;910:012028.10.1088/1742-6596/910/1/012028Search in Google Scholar

[49] Dong D, Fang P, Bock Y, Cheng M, Miyazaki SI. Anatomy of apparent seasonal variations from GPS-derived site position time series. J Geophys Res. 2002;107:2075.10.1029/2001JB000573Search in Google Scholar

[50] Liu X, Wu Y, Jiang Z, Zhan W, Li Q, Wei W, et al. Preseismic deformation in the seismogenic zone of the Lushan MS7.0 earthquake detected by GPS observations. Sci China: Earth Sci. 2015;58:1592–601.10.1007/s11430-015-5128-0Search in Google Scholar

[51] Zhan W, Yang B, Zhang F, Liu ZH, Meng X. Effect of Selecting GNSS Frame sits on The Time Series of Displacement. J Geodesy Geodynamics. 2011;31(1):91–4.Search in Google Scholar

[52] Wang S. Consideration on both the geodetic deformation monitoring of the north-south seismic belts and the big earthquake predictions. Earthq Res Sichan. 2011;36(4):13–5.Search in Google Scholar

[53] Wang S, Zhang X, Zhu Y, Zhang X, Xue F. Study on the characteristics of strong earthquake deformation anomalies and for ecasting index in the north segment of the south—north seismic belt of China. Northwest Seismol J. 2012;34(2):105–13.Search in Google Scholar

[54] Jiang Z, Yang G, Fang Y, Zhang X, Wu Y, Ping J. Study on the relations between distribution of crustal movement and strong earthquakes using GPS survey. Recent Dev World Seismology. 2007;7:32–42.Search in Google Scholar

[55] Shao Z, Feng W, Wang P, Yin X. A study review on characteristic of seismic activity of active-tectonic block boundaries in mainland China. Seismol Geol. 2020;2:271–82.Search in Google Scholar

[56] Ma H, Wu Y, Feng J, Li X, Xu H, Wang Q. Crustal movement characteristics in the northeastern region of the qinghai-tibet plateau based on GPS reference station. Earth Environ Sci. 2019;384:012209.10.1088/1755-1315/384/1/012209Search in Google Scholar

Received: 2021-07-20
Revised: 2022-05-05
Accepted: 2022-06-22
Published Online: 2023-07-15

© 2022 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Diagenesis and evolution of deep tight reservoirs: A case study of the fourth member of Shahejie Formation (cg: 50.4-42 Ma) in Bozhong Sag
  3. Petrography and mineralogy of the Oligocene flysch in Ionian Zone, Albania: Implications for the evolution of sediment provenance and paleoenvironment
  4. Biostratigraphy of the Late Campanian–Maastrichtian of the Duwi Basin, Red Sea, Egypt
  5. Structural deformation and its implication for hydrocarbon accumulation in the Wuxia fault belt, northwestern Junggar basin, China
  6. Carbonate texture identification using multi-layer perceptron neural network
  7. Metallogenic model of the Hongqiling Cu–Ni sulfide intrusions, Central Asian Orogenic Belt: Insight from long-period magnetotellurics
  8. Assessments of recent Global Geopotential Models based on GPS/levelling and gravity data along coastal zones of Egypt
  9. Accuracy assessment and improvement of SRTM, ASTER, FABDEM, and MERIT DEMs by polynomial and optimization algorithm: A case study (Khuzestan Province, Iran)
  10. Uncertainty assessment of 3D geological models based on spatial diffusion and merging model
  11. Evaluation of dynamic behavior of varved clays from the Warsaw ice-dammed lake, Poland
  12. Impact of AMSU-A and MHS radiances assimilation on Typhoon Megi (2016) forecasting
  13. Contribution to the building of a weather information service for solar panel cleaning operations at Diass plant (Senegal, Western Sahel)
  14. Measuring spatiotemporal accessibility to healthcare with multimodal transport modes in the dynamic traffic environment
  15. Mathematical model for conversion of groundwater flow from confined to unconfined aquifers with power law processes
  16. NSP variation on SWAT with high-resolution data: A case study
  17. Reconstruction of paleoglacial equilibrium-line altitudes during the Last Glacial Maximum in the Diancang Massif, Northwest Yunnan Province, China
  18. A prediction model for Xiangyang Neolithic sites based on a random forest algorithm
  19. Determining the long-term impact area of coastal thermal discharge based on a harmonic model of sea surface temperature
  20. Origin of block accumulations based on the near-surface geophysics
  21. Investigating the limestone quarries as geoheritage sites: Case of Mardin ancient quarry
  22. Population genetics and pedigree geography of Trionychia japonica in the four mountains of Henan Province and the Taihang Mountains
  23. Performance audit evaluation of marine development projects based on SPA and BP neural network model
  24. Study on the Early Cretaceous fluvial-desert sedimentary paleogeography in the Northwest of Ordos Basin
  25. Detecting window line using an improved stacked hourglass network based on new real-world building façade dataset
  26. Automated identification and mapping of geological folds in cross sections
  27. Silicate and carbonate mixed shelf formation and its controlling factors, a case study from the Cambrian Canglangpu formation in Sichuan basin, China
  28. Ground penetrating radar and magnetic gradient distribution approach for subsurface investigation of solution pipes in post-glacial settings
  29. Research on pore structures of fine-grained carbonate reservoirs and their influence on waterflood development
  30. Risk assessment of rain-induced debris flow in the lower reaches of Yajiang River based on GIS and CF coupling models
  31. Multifractal analysis of temporal and spatial characteristics of earthquakes in Eurasian seismic belt
  32. Surface deformation and damage of 2022 (M 6.8) Luding earthquake in China and its tectonic implications
  33. Differential analysis of landscape patterns of land cover products in tropical marine climate zones – A case study in Malaysia
  34. DEM-based analysis of tectonic geomorphologic characteristics and tectonic activity intensity of the Dabanghe River Basin in South China Karst
  35. Distribution, pollution levels, and health risk assessment of heavy metals in groundwater in the main pepper production area of China
  36. Study on soil quality effect of reconstructing by Pisha sandstone and sand soil
  37. Understanding the characteristics of loess strata and quaternary climate changes in Luochuan, Shaanxi Province, China, through core analysis
  38. Dynamic variation of groundwater level and its influencing factors in typical oasis irrigated areas in Northwest China
  39. Creating digital maps for geotechnical characteristics of soil based on GIS technology and remote sensing
  40. Changes in the course of constant loading consolidation in soil with modeled granulometric composition contaminated with petroleum substances
  41. Correlation between the deformation of mineral crystal structures and fault activity: A case study of the Yingxiu-Beichuan fault and the Milin fault
  42. Cognitive characteristics of the Qiang religious culture and its influencing factors in Southwest China
  43. Spatiotemporal variation characteristics analysis of infrastructure iron stock in China based on nighttime light data
  44. Interpretation of aeromagnetic and remote sensing data of Auchi and Idah sheets of the Benin-arm Anambra basin: Implication of mineral resources
  45. Building element recognition with MTL-AINet considering view perspectives
  46. Characteristics of the present crustal deformation in the Tibetan Plateau and its relationship with strong earthquakes
  47. Influence of fractures in tight sandstone oil reservoir on hydrocarbon accumulation: A case study of Yanchang Formation in southeastern Ordos Basin
  48. Nutrient assessment and land reclamation in the Loess hills and Gulch region in the context of gully control
  49. Handling imbalanced data in supervised machine learning for lithological mapping using remote sensing and airborne geophysical data
  50. Spatial variation of soil nutrients and evaluation of cultivated land quality based on field scale
  51. Lignin analysis of sediments from around 2,000 to 1,000 years ago (Jiulong River estuary, southeast China)
  52. Assessing OpenStreetMap roads fitness-for-use for disaster risk assessment in developing countries: The case of Burundi
  53. Transforming text into knowledge graph: Extracting and structuring information from spatial development plans
  54. A symmetrical exponential model of soil temperature in temperate steppe regions of China
  55. A landslide susceptibility assessment method based on auto-encoder improved deep belief network
  56. Numerical simulation analysis of ecological monitoring of small reservoir dam based on maximum entropy algorithm
  57. Morphometry of the cold-climate Bory Stobrawskie Dune Field (SW Poland): Evidence for multi-phase Lateglacial aeolian activity within the European Sand Belt
  58. Adopting a new approach for finding missing people using GIS techniques: A case study in Saudi Arabia’s desert area
  59. Geological earthquake simulations generated by kinematic heterogeneous energy-based method: Self-arrested ruptures and asperity criterion
  60. Semi-automated classification of layered rock slopes using digital elevation model and geological map
  61. Geochemical characteristics of arc fractionated I-type granitoids of eastern Tak Batholith, Thailand
  62. Lithology classification of igneous rocks using C-band and L-band dual-polarization SAR data
  63. Analysis of artificial intelligence approaches to predict the wall deflection induced by deep excavation
  64. Evaluation of the current in situ stress in the middle Permian Maokou Formation in the Longnüsi area of the central Sichuan Basin, China
  65. Utilizing microresistivity image logs to recognize conglomeratic channel architectural elements of Baikouquan Formation in slope of Mahu Sag
  66. Resistivity cutoff of low-resistivity and low-contrast pays in sandstone reservoirs from conventional well logs: A case of Paleogene Enping Formation in A-Oilfield, Pearl River Mouth Basin, South China Sea
  67. Examining the evacuation routes of the sister village program by using the ant colony optimization algorithm
  68. Spatial objects classification using machine learning and spatial walk algorithm
  69. Study on the stabilization mechanism of aeolian sandy soil formation by adding a natural soft rock
  70. Bump feature detection of the road surface based on the Bi-LSTM
  71. The origin and evolution of the ore-forming fluids at the Manondo-Choma gold prospect, Kirk range, southern Malawi
  72. A retrieval model of surface geochemistry composition based on remotely sensed data
  73. Exploring the spatial dynamics of cultural facilities based on multi-source data: A case study of Nanjing’s art institutions
  74. Study of pore-throat structure characteristics and fluid mobility of Chang 7 tight sandstone reservoir in Jiyuan area, Ordos Basin
  75. Study of fracturing fluid re-discharge based on percolation experiments and sampling tests – An example of Fuling shale gas Jiangdong block, China
  76. Impacts of marine cloud brightening scheme on climatic extremes in the Tibetan Plateau
  77. Ecological protection on the West Coast of Taiwan Strait under economic zone construction: A case study of land use in Yueqing
  78. The time-dependent deformation and damage constitutive model of rock based on dynamic disturbance tests
  79. Evaluation of spatial form of rural ecological landscape and vulnerability of water ecological environment based on analytic hierarchy process
  80. Fingerprint of magma mixture in the leucogranites: Spectroscopic and petrochemical approach, Kalebalta-Central Anatolia, Türkiye
  81. Principles of self-calibration and visual effects for digital camera distortion
  82. UAV-based doline mapping in Brazilian karst: A cave heritage protection reconnaissance
  83. Evaluation and low carbon ecological urban–rural planning and construction based on energy planning mechanism
  84. Modified non-local means: A novel denoising approach to process gravity field data
  85. A novel travel route planning method based on an ant colony optimization algorithm
  86. Effect of time-variant NDVI on landside susceptibility: A case study in Quang Ngai province, Vietnam
  87. Regional tectonic uplift indicated by geomorphological parameters in the Bahe River Basin, central China
  88. Computer information technology-based green excavation of tunnels in complex strata and technical decision of deformation control
  89. Spatial evolution of coastal environmental enterprises: An exploration of driving factors in Jiangsu Province
  90. A comparative assessment and geospatial simulation of three hydrological models in urban basins
  91. Aquaculture industry under the blue transformation in Jiangsu, China: Structure evolution and spatial agglomeration
  92. Quantitative and qualitative interpretation of community partitions by map overlaying and calculating the distribution of related geographical features
  93. Numerical investigation of gravity-grouted soil-nail pullout capacity in sand
  94. Analysis of heavy pollution weather in Shenyang City and numerical simulation of main pollutants
  95. Road cut slope stability analysis for static and dynamic (pseudo-static analysis) loading conditions
  96. Forest biomass assessment combining field inventorying and remote sensing data
  97. Late Jurassic Haobugao granites from the southern Great Xing’an Range, NE China: Implications for postcollision extension of the Mongol–Okhotsk Ocean
  98. Petrogenesis of the Sukadana Basalt based on petrology and whole rock geochemistry, Lampung, Indonesia: Geodynamic significances
  99. Numerical study on the group wall effect of nodular diaphragm wall foundation in high-rise buildings
  100. Water resources utilization and tourism environment assessment based on water footprint
  101. Geochemical evaluation of the carbonaceous shale associated with the Permian Mikambeni Formation of the Tuli Basin for potential gas generation, South Africa
  102. Detection and characterization of lineaments using gravity data in the south-west Cameroon zone: Hydrogeological implications
  103. Study on spatial pattern of tourism landscape resources in county cities of Yangtze River Economic Belt
  104. The effect of weathering on drillability of dolomites
  105. Noise masking of near-surface scattering (heterogeneities) on subsurface seismic reflectivity
  106. Query optimization-oriented lateral expansion method of distributed geological borehole database
  107. Petrogenesis of the Morobe Granodiorite and their shoshonitic mafic microgranular enclaves in Maramuni arc, Papua New Guinea
  108. Environmental health risk assessment of urban water sources based on fuzzy set theory
  109. Spatial distribution of urban basic education resources in Shanghai: Accessibility and supply-demand matching evaluation
  110. Spatiotemporal changes in land use and residential satisfaction in the Huai River-Gaoyou Lake Rim area
  111. Walkaway vertical seismic profiling first-arrival traveltime tomography with velocity structure constraints
  112. Study on the evaluation system and risk factor traceability of receiving water body
  113. Predicting copper-polymetallic deposits in Kalatag using the weight of evidence model and novel data sources
  114. Temporal dynamics of green urban areas in Romania. A comparison between spatial and statistical data
  115. Passenger flow forecast of tourist attraction based on MACBL in LBS big data environment
  116. Varying particle size selectivity of soil erosion along a cultivated catena
  117. Relationship between annual soil erosion and surface runoff in Wadi Hanifa sub-basins
  118. Influence of nappe structure on the Carboniferous volcanic reservoir in the middle of the Hongche Fault Zone, Junggar Basin, China
  119. Dynamic analysis of MSE wall subjected to surface vibration loading
  120. Pre-collisional architecture of the European distal margin: Inferences from the high-pressure continental units of central Corsica (France)
  121. The interrelation of natural diversity with tourism in Kosovo
  122. Assessment of geosites as a basis for geotourism development: A case study of the Toplica District, Serbia
  123. IG-YOLOv5-based underwater biological recognition and detection for marine protection
  124. Monitoring drought dynamics using remote sensing-based combined drought index in Ergene Basin, Türkiye
  125. Review Articles
  126. The actual state of the geodetic and cartographic resources and legislation in Poland
  127. Evaluation studies of the new mining projects
  128. Comparison and significance of grain size parameters of the Menyuan loess calculated using different methods
  129. Scientometric analysis of flood forecasting for Asia region and discussion on machine learning methods
  130. Rainfall-induced transportation embankment failure: A review
  131. Rapid Communication
  132. Branch fault discovered in Tangshan fault zone on the Kaiping-Guye boundary, North China
  133. Technical Note
  134. Introducing an intelligent multi-level retrieval method for mineral resource potential evaluation result data
  135. Erratum
  136. Erratum to “Forest cover assessment using remote-sensing techniques in Crete Island, Greece”
  137. Addendum
  138. The relationship between heat flow and seismicity in global tectonically active zones
  139. Commentary
  140. Improved entropy weight methods and their comparisons in evaluating the high-quality development of Qinghai, China
  141. Special Issue: Geoethics 2022 - Part II
  142. Loess and geotourism potential of the Braničevo District (NE Serbia): From overexploitation to paleoclimate interpretation
Downloaded on 11.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/geo-2022-0387/html
Scroll to top button