Abstract
East Asia is bounded by the Indian plate to the southwest and the Pacific and Philippine plates to the east, and has undergone complex tectonic evolution since ~55 Ma. In this study, we collect and process three sources of GPS datasets, including GPS observations, GPS positioning time series, and published GPS velocities, to derive unified velocity and strain rate fields for East Asia. We observed southward movement and arc-parallel extension along the Ryukyu Arc and propose that the maximum principal stress axis (striking NEE) in North China could be mainly induced by westward subduction of the Pacific plate and the southward movement of the Ryukyu Arc. The large EW-trending sinistral shear zone that bounds North China has been created by eastward movement of South China to the south and westward subduction of the Pacific plate to the north. GPS velocity profiles and strain rates also demonstrate that crustal deformation in mainland China is controlled by northeastward collision of the Indian plate into Eurasia and westward subduction of the Pacific and Philippine Sea plates beneath Eurasia. In particular, the India-Eurasia continental collision has the most extensive impact, which can reach as far as the southern Lake Baikal. The viscous behavior of the subducting Pacific slab also drives interseismic deformation of North China. The crustal deformation caused by Philippine oceanic subduction is small and is limited to the region between the southeast coast of mainland China and Taiwan island. However, the principal compressional strain around eastern Taiwan is the largest in the region.
Similar content being viewed by others
Introduction
Cenozoic deformation of East Asia is complex due to the collision of India into Eurasia to the southwest and the Western Pacific and Philippine oceanic plates subducting under Eurasia to the east (Fig. 1). Since ~55 million years ago, the continued northward movement of India has uplifted the Himalayan orogen and Tibetan Plateau, and dominates the tectonic evolution of western mainland China1,2,3. In contrast, subduction of the Pacific and Philippine plates has created island arcs, marginal seas, continental rifting, and plays a crucial role in the tectonic stress in eastern China3,4,5.
The convergence between India and Asia causes intercontinental deformation over a large distances of up to ~3000 km, reaching as far as Lake Baikal6. Tomographic images also indicate that the dynamic processes of the Western Pacific and East Asia are responsible for the formation of continental rifts and active intraplate volcanoes in Northeast China7,8. However, the dynamic sources of tectonic motion in East Asia remain elusive. Some have proposed that the tectonics of East Asia can be entirely related to the India-Eurasia collision1,9, while others emphasize that the stress field of East Asia is controlled by the combined convergence of the Indian, Pacific, and Philippine plates10,11. The advent of GPS has provided an accurate and efficient approach for measuring regional and large scale crustal deformation, which are crucial in the attempt to constrain the dynamic processes affecting East Asia.
Despite extensive previous studies that have provided GPS-derived regional crustal motions for regions such as the Tibetan Plateau, Tian Shan, western Pacific subduction zone, and even a sparse velocity field for all of Asia2,12,13,14,15,16,17, a lack of a high resolution and accurate GPS velocity field has prevented the elucidation of crustal motion across the whole of East Asia. Our primary purpose in this study is to generate present-day unified GPS velocity and strain rate fields for East Asia, in order to analyze the patterns of crustal movement across key regions around East Asia. We also discuss reasons for the tectonic stress field direction in North China and compare the influences of continental collision and oceanic subduction on the overall stress field. The GPS velocity field generated and discussed herein is helpful for understanding various regional kinematics and geodynamics that drive continental deformation, and is useful for maintaining modernized geodetic datum.
GPS data
We assembled three sources of GPS datasets to derive a geodetically consistent crustal velocity field for East Asia. The first dataset is raw GPS observations measured in mainland China. The second is continuous GPS time series data available from public online websites (http://www.gsi.go.jp/, Geospatial Information Authority of Japan, GSI, and http://tec.earth.sinica.edu.tw/, Taiwan Earthquake Research Center, TEC), which provide GPS data for the Western Pacific and Philippine subduction zones. The third one is GPS velocities from existing literature, which mainly focus on the India-Asia collision regime.
The GPS data collected and measured in mainland China are from the Crustal Movement Observation Network of China (CMONOC). These GPS stations were surveyed in 1998, 1999, 2001, 2004, and 2007. We did not use any post-2007 data in order to avoid effects of the 2008 Mw 7.9 Wenchuan and 2011 Mw 9.0 Tohoku earthquakes. The GAMIT software package18 was employed to get daily loosely constrained solutions for satellite orbits and coordinates of regional stations. Approximately 100 evenly distributed ITRF core GPS sites were also processed using the GAMIT software. The VMF1 tropospheric mapping function was used to calculate the wet and hydrostatic zenith troposhperic delays. The GLOBK software19 was used to estimate daily positions and uncertainties by combining daily regional and global solutions, and then the daily free network solutions was transformed into ITRF200820. The GPS sites used to define the global ITRF2008 are showed in Fig. S1 in the Supplementary Material. Due to the coseismic deformation of the 2004 Mw 9.1 Sumatra earthquake, we also estimated offsets by using the coseismic slip model to correct the time series. Section A2 and Fig. S2 in the Supplementary Material show the coseismic displacements of the campaign mode GPS stations. We considered the postseismic displacements of this earthquake to be negligible, for the purposes of this study. Then we used a linear function to fit the GPS time series, in order to estimate secular velocity rates with respect to ITRF2008.
The continuous GPS positioning time series used span 1996~2015 and were provided by GSI and TEC. In order to mitigate the seasonal variations in these GPS time series, we only estimated velocities for time series that were at least 2.5 years long. We preferred ruling out parts of GPS time series that also demonstrated prominent transient deformation rather than employing proper models to correct them, because modeling transient deformation, including postseismic deformation, slow slip, or other non-linear processes, can introduce a large number of model dependencies15. Finally, we applied a linear trend and a seasonal variation to fit the GPS time series and infer interseismic velocity. It is well known that errors in daily GPS position estimates are time correlated21. Using a purely white noise model would cause low velocity uncertainties. Therefore, we employed the Maximum Likelihood Estimation technique using the CATS software21 to estimate velocity uncertainties, with the assumption that the error can be categorized by a power law noise model. More details regarding our estimation procedure are included in section A3 of the Supplementary Material.
The GPS velocities derived from previously published studies are mostly distributed in the India-Asia collision regions. Table 1 lists previous studies used in GPS compilation and the regions our study concerned. These secular velocities were originally determined with respect to different global (e.g., ITRF2000, ITRF2005 and ITRF2008) or regional reference frames (e.g., Eurasian and Indian plates).
Because we would like to obtain GPS movement across mainland China and its adjacent regions, a regional stable reference frame within mainland China is suitable. Velocities for 62 GPS stations located on the stable South China block were selected to estimate the Euler pole by applying constraints that minimized motions within the interior of the South China block. We removed potential outliers from the velocities through an iterative procedure: we removed stations with the largest post-fit residual and re-estimated the angular velocity until all the post-fit residuals were less than 2 mm/yr. Finally velocities for 50 stations (blue triangles, Fig. 1) were used to define the South China block. The estimated pole position of the South China block is located at 53.420 ± 1.275°N, −101.102 ± 1.139°W, with a rotation rate of 0.311 ± 0.006°/Ma. Section A4 in Supplementary shows the detailed procedure. Then we transformed the velocity field of mainland China into a regional reference frame with respect to the South China block by using a rigid block rotation.
We chose GPS velocities within mainland China as the reference data set, and velocities obtained from other published papers and estimated from GPS time series were first transformed to this regional reference frame. Here, we employed the two-dimensional Helmert transformation (east and north velocities). The four parameters were determined (two translations, one scaling and one rotation about the Z-axis) for our data and other velocities to be transformed with at least two identical points, using the least squares method. More than two identical points, along with those with sums of post-fit residual squares, were used to estimate the transformation parameters. The post-fit residuals following the Helmert transformation were generally less than 1.0 and 1.0 mm/yr in the eastward and northward velocities of identical points, respectively. We present the GPS velocity field for East Asia with respect to a stable South China block (Fig. 2). Section A6 in the Supplementary Material provides the estimated GPS rates and errors.
We also calculated a strain rate field from GPS velocities using the method proposed by ref.22. This technique employs spherical wavelet-based multi-scale decomposition to calculate a continuous velocity field from discrete GPS observations. The coefficients of the spherical wavelets are then used to directly compute strain rates23. The advantages of this technique are that (1) the velocity field is decomposed into multiple scales at all stations explicitly and consistently, and (2) the minimum scale of the estimated velocity field for a particular location is controlled by the local spatial density of observations.
Results
Velocity field
The horizontal velocity field (Fig. 2) shows first-order deformation in East Asia. To the southwest, the Indian plate moves northward at rates of ~37 mm/year. The Tibetan Plateau experiences NS shortening and EW extension. The Tian Shan orogen, especially the western section, absorbs a large amount of crustal shortening, with rates of 15~18 mm/yr. With respect to South China, the east side of the eastern Himalayan syntaxis (EHS) shows clockwise rotation, while its west side demonstrates movement in the NNE-NE direction.
To the east, due to the subduction of the Pacific plate, northern Japan moves westward. Velocity rates are ~40 mm/yr along the east coast of Japan and decrease to ~20 mm/yr along the coast of the Japan Sea (Fig. 3a). Crossing over the Japan Sea, velocities in southeast Russia and Northeast China trend to the NW, with rates of 6~8 mm/yr. Further south, there is a velocity gradient trending NWW, which is located in the Zhangjiakou-Bohai Seismic Zone (ZBSZ) and forms the northern boundary of North China24. Velocity rates decrease from 4~5 mm/yr in the north to 1~2 mm/yr in the south, showing sinistral slip motion.
To the southeast, the Philippine plate subducts NW under the Ryukyu Arc with a rate of ~80 mm/yr (Fig. 3b). The GPS velocities in the Ryukyu Arc increase from 14 mm/yr in the north to 40 mm/yr in the south, while the velocity directions gradually shift from southwest to south. Eastern Taiwan undergoes significant contraction, at rates of approximately 20 mm/yr. Further south, the northern Luzon arc moves NW at a rate of ~80 mm/yr.
Strain rate fields
The key feature of the principal strain rate field (Fig. 4) is that the contractional strain rates around the three plate boundary zones are the most prominent. Especially the Philippine subduction zone, which generates the highest compressive strain rates of greater than 100 × 10−8/yr in a NWW direction located along the east coast of Taiwan. The principal compressive strain rate along the east coast of Japan caused by Pacific plate subduction is also higher than along the Himalayan Main Frontal Thrust, which is induced by the convergence of the Indian plate. The western Tian Shan also has high compressive strain rates, as does the Himalayan Main Frontal Thrust. The strain rate across the Tibtean Plateau is consistent with ref.25, excluding the magnitudes of compressive strain rate in Nepal. The compressive strain rate25 in this region is larger than ours. By comparing the two GPS velocity fields, we found our GPS velocities are more than ref.25 in the Himalayan orogen and the Indian plate. However, our derived GPS strain rate is similar to ref.26 in the Himalayan orogen.
Figure 5 shows the maximum shear strain rates in East Asia. We observed that maximum shear strain rates of 10 × 10−7/yr, 3 × 10−7/yr, and 1 × 10−7/yr around Taiwan, the Ryukyu Arc, and Japan, respectively. The Indo-Burmese range, eastern and western Himalayan syntaxises, and Sagaing fault also have relatively high shear strain rates.
Discussion
Southward movement along the ryukyu arc
As shown in Fig. 3, crustal deformation in Japan shows westward movement along the Japan Trench and NWW movement along the Nankai Trough, which are both attributed to the subduction of the Pacific plate. However, further south the velocity vectors are in a S-SSW direction along the Ryukyu Arc. The geophysical, geological, and GPS data show that there is substantial contraction along the Nankai Trough, but no contraction along the Ryukyu Trench27. Moment tensor solutions also show that arc-parallel extension exists throughout the Ryukyu Arc region, excluding its northeastern edge28.
GPS-derived strain rates for the Ryukyu Arc show arc-perpendicular contraction and arc-parallel extension (Fig. 4). ref.28 proposed that the arc-parallel extension might not be due to the effects of oblique subduction alone; the back-arc opening process could play a crucial role in this arc-parallel extension28,29,30. Tectonic morphology and geophysical soundings reveal that en echelon central depressions are distributed along the axis of the Okinawa Trough. Both sides of the depressions contain normal faults extending tens of kilometers, indicating that these depressions are grabens. This implies that the trough has been in a state of extension29,31. Independently calculated opening rates for the Okinawa Trough are 10~50 mm/yr29, consistent with our GPS velocity rates of 14~40 mm/yr for the region.
Geodynamics of stress orientation and strain distribution in north china
North China is bounded by the Pacific plate to the east, the Ordos block to the west, stable South China to the south, and Northeast China (belongs to the Amurian block) to the north. ref.32 used focal mechanisms, fault striations inferred from Quaternary fault slip measurements, stress relief, and hydraulic fracturing data to estimate the tectonic stress field for China. Their results indicated that the maximum principal stress axis strikes NEE in North China. However, a compressive stress orientation of NEE cannot be attributed only to the northward collision of the Indian plate and the Pacific plate westward subduction10, since the Pacific plate subduction in a NWW direction seems unlikely to generate such a large deflected angle.
Our GPS velocity field shows that Japan moves westward due to Pacific plate subduction, while the Ryukyu Arc moves southward as a result of Okinawa Trough extension (Fig. 3). Therefore, the westward subduction of the Pacific plate and southward movement along the Ryukyu Arc can generate a principal compressive stress in the NE-SW direction, resulting in the maximum principal stress axis striking NEE in North China.
Components of GPS velocity profiles (Fig. 6a,b) parallel to the ZBSZ across North China show that the ZBSZ has a sinistral slip rate of 3 mm/yr. The GPS velocity field, with respect to South China (Fig. 2), demonstrates that the directions of GPS velocities in Northeast China are primarily to the NWW and the magnitudes decrease gradually from east to west. Consistently, the magnitudes of principal compressive strain rates in Northeast China also decrease from east to west and their azimuths change from EW to NE-SW. This suggests that the westward motion of the Pacific plate plays a crucial role in driving NW-ward movement in Northeast China.
The kinematics of the ZBSZ may also be caused by westward motion of the Pacific plate. Through analyzing GPS velocity profiles, ref.24 found that interseismic deformation in North China absorbs left-lateral shear of 6.0 mm/yr over an ~1100 km distance. They proposed a model that an EW-trending sinistral shear zone (the north and south boundaries being the ZBSZ and Dabie Shan, respectively) bounds a group of NEE-trending dextral strike-slip faults within North China33. ref.24 suggests that the sinistral shear strain across North China related to the stable Amurian block is created by eastward movement of South China, which is driven by eastward expansion of the Tibetan Plateau to the south. Additionally, our results indicate that the sinistral shear across North China could be caused by NW-ward movement of Northeast China caused by westward subduction of the Pacific plate to the north.
Kinematics influence of the Indian plate colliding with eurasia
As mentioned previously, the dynamic sources of tectonic deformation in East Asia are still in question. Some believe that the tectonic features of East Asia are mainly attributable to the collision of the Indian plate with Eurasia1,9, while others suggest that the stress field of East Asia is governed by the combined influences of the Indian, Pacific, and Philippine plates converging with the Eurasian plate10,11.
Our unified GPS velocity and strain rate fields for East Asia can shed new light on the range of crustal deformation caused by the plates colliding with Asia. In order to reveal crustal deformation along the direction of subduction for each plate explicitly, we constructed three velocity profiles across the Himalayan thrust zone (Fig. 6c), Japan subduction zone (Fig. 6d), and Taiwan plate boundary zone (Fig. 6e). The approximate profile locations are shown in Fig. 6f.
Figure 6c shows the velocity components parallel to N30°E, from the Indian plate in the south to Lake Baikal in the north. We observe significant shortening rates of ~30 mm/yr along N30°E across the Tibetan Plateau between the HMFT and the Altyn Tagh fault, over a distance of 1700 km. Further north, the N30°E components of velocity gradually decrease from 5 mm/yr on the Altyn Tagh to −2 mm/yr on southern Lake Baikal, over a distance of 1500 km. The azimuths of the principal compressive strain rates change from N in the HMFT, to NNE in the Tibetan Plateau, and to NE in Mongolia and southern Lake Baikal. The magnitudes of principal compressive strain rates decrease from 2 to 0.1 × 10−8/yr along this profile (Fig. 4).
Geophysical studies reveal that the Indian plate is underthrusting central Tibet at a low angle and reaches north to the Bangong–Nujiang Suture (~500 km north of the HMFT)5,34,35. However, the gradual decrease in velocities and compressive strain rates in the NEE-NE direction suggests that the lithospheric deformation generated by the convergence of the Indian plate has extended ~3000 km away from the HMFT on the Earth’s surface1.
We employed spherical linear block model36,37 to investigate the crustal movement due to elastic strain accumulation from locking on the HMFT. The block geometry, dip angle and locking depth are determined from ref.36. The block model result suggests that crustal deformation result from elastic strain accumulation on the locking HMFT is confined within south of the Jali fault, and directs southward (Fig. 7(f)). The estimated slip rates along faults are provided in the Supplementary Material. ref.36 also used block model constrained by GPS velocities to infer the extent to which relative motion of the Indian and Asian plates is partitioned between localized slip on major faults and distributed deformation in the Tibetan Plateau. Their results suggest that microplate rotation, interseismic elastic strain accumulation, and fault slip on the boundaries of microplates account for surface movement in the Tibetan Plateau. Internal block deformation occurs within some micro-plate, such as the Himalayan Range, Jiali, and west-central plateau blocks. Therefore, the observed northward movement within Tibet and further north is mainly from block rotation and internal strain, which is driven by boundary forces possibly38.
Kinematics influence of the pacific plate subducting under eurasia
Figure 6d shows the velocity components parallel to N65°W, from southern Japan through Korea and the ZBSZ, to the north side of the Ordos block. In general, velocity components that are parallel to N65°W decay with distance from the plate boundary, but the velocity gradient across southern Japan is the largest. Although there are no GPS sites in the Korea Strait or the Yellow Sea, we obtain NW-SE contraction rates of about 4 and 2 mm/yr, respectively, from this velocity profile. Along the ZBSZ in North China, the change in velocity components towards the west is insignificant. The principal compressive strain rates decrease from 15 × 10−8/yr near central southern Japan, to 0.01 × 10−8/yr near the northern side of the Ordos block, and the corresponding azimuths rotate from E-W to NEE-SWW (Fig. 4).
Components of velocity that are parallel to N110°E across eastern North China (Fig. 6a) suggest that eastern North China, from the Dabie Shan to the ZBSZ, experiences a pronounced sinistral shear, while velocity components parallel to N110°E across western North China (Fig. 6b) suggest that sinistral shear strain is primarily concentrated in the vicinity of the ZBSZ, not the entire North China region. Therefore, the influence of Pacific plate subduction on North China weakens from east to west, and is mainly distributed in the eastern region.
Block model constrained by GPS velocity was used to estimate plate movement, fault slip rates and interplate coupling37. However, movement due to elastic strain accumulation from locking on the subduction zone interfaces is NEE, which is opposite to the subduction direction of the Pacific plate. The westward motion of Japan that gradually decreases from east to west reflects the E-W contractional strain built up by inter-plate coupling39. Interseismic deformation induced by subduction is controlled by viscous behavior of the mantle40. Tomographic images show that the subducting Pacific slab is represented by a prominent high-velocity zone, and the flat Pacific slab extends from Japan to Beijing, China, over 2300 km, at a depth of 660 km5,7,41,42. This suggests that Pacific subduction beneath Asia is the major cause for active tectonics and mantle dynamics in a broad region between the Japan Trench and East Asia, rather than the India-Asia collision7. Driving forces from Pacific plate subduction can play a key role in crustal movement decaying from east to west in eastern mainland China.
Kinematics influence of the philippine plate subducting under eurasia
Figure 6e shows the velocity components that are parallel to N60°W from the east coast of Taiwan, through the Taiwan Strait, to South China. From east to west, the velocity components parallel to N60°W demonstrate large convergence rates of ~50 mm/yr across Taiwan, but only a small velocity gradient across South China. The strain rates along the profile show that eastern Taiwan is experiencing NW-SE intensive contraction, with maximum rates of 110 × 10−8/yr, and the southeast coast of South China, across the Taiwan Strait, undergoes E-W contraction with maximum contraction rate of 3.5 × 10−8/yr. The principal strain rates are insignificant in South China.
Taiwan is situated on a plate boundary zone between the Eurasian and Philippine plates. The Philippine plate subducts beneath the Eurasian plate along the Ryukyu Arc in the east, while the Eurasian plate subducts eastward beneath the Philippine plate. P- and S-wave velocity structures and Pn tomography show high velocity anomalies beneath the Ryukyu and Luzon arcs at depths ranging from 100~200 km, and the anomalies related to the Eurasian and Philippine slabs appear to be connected34,43. ref.43 proposed that the origin of the island of Taiwan was induced by collision at the edge point of the two subduction zones, resulting in prominent shortening and strong stresses confined to the eastern Taiwan region. Comparing the geophysics, geology, and GPS data with the theoretical results calculated from the Euler parameter of the Philippine Sea plate, ref.27 proposed that the collision of the Eurasian and Philippine plates in the area of Taiwan between 21.5°N and 24.2°N causes strong compression in the southeast part of the Chinese continent, but with a limited scope.
Range of kinematics influence of three plates colliding with eurasia
The Indian-Eurasian continental collision and continued northward movement of the Indian plate are responsible for crustal deformation that reaches as far as southern Lake Baikal. The Pacific plate subducting under Eurasia may affect tectonic movement in an area from the Japan subduction zone to North China, while the northwestward subduction of the Philippine plate influences the surface deformation of the southeastern coast of the Chinese continent. The northeastward collision of the Indian plate with Asia and the westward subduction of the Pacific and Philippine Sea plates may have played a significant role in the formation and tectonic evolution of mainland China5.
In this paper, we emphasize that the first-order drive forces resulting from India-Eurasia collision and subduction of the Pacific and Philippine plates beneath Eurasia control the deformation of mainland China. The regional tectonic deformation of China’s continent is also characterized by the coupling of rigid block movement and continuous deformation. Regions of high rigidity behave as rigid block-like movement, while those of low rigidity are dominated by continuous deformation. Rheological flow in the lower crust and upper mantle plays an important role in controlling deformation of the upper crust44. Using thin sheet model, ref.38 believed tectonic deformation across Asia is driven by the balance of boundary forces and buoyancy stresses generated by gravitational potential energy gradients. Their result shows that deformation in contractional regions (such as Himalayas, Tibet and Tian Shan) is well result from boundary forces with strong coupling at the India-Eurasia collision zone. But southeastward movement of North and South China requires buoyancy stresses across the Pacific and Philippine subductions.
The GPS velocity filed used in ref.38 is respect to Eurasia, so the GPS velocities observed in South and North China move to east and southeastward. Therefore, buoyancy forces are needed to resist to boundary forces arising from subduction of oceanic plates beneath Asia. However, our GPS-derived velocities are respect to South China, a stable regional block within mainland China, and North China moves northwestward which is consistent with the direction of oceanic plates’ subduction. Whether buoyancy forces still required needs further simulation constrained with our GPS velocity filed across East Asia.
Conclusions
We used GPS measurements to determine the present day unified GPS velocity field for East Asia. The velocities and calculated strain rates show first-order crustal movement induced by the Indian plate penetrating into Eurasia, and the Western Pacific and Philippine Sea plates subducting beneath Eurasia. Our GPS velocity results demonstrate southward movement and arc-parallel extension along the Ryukyu Arc, which may be attributed to back-arc extension. We propose that westward movement of the Pacific plate and southward movement along the Ryukyu Arc could produce a principal compressive stress in the NE-SW direction, resulting in a maximum principal stress axis striking NEE in North China. The large EW-trending sinistral shear zone bounding North China could have possibly been created by eastward movement of South China to the south, and westward subduction of the Pacific plate to the north.
Using GPS velocity profiles along the subduction direction of the Indian, Pacific and Philippine plates separately, we discussed the range of kinematics influence of these subduction boundaries. Our results indicate that tectonic movement in mainland China is controlled by the northeastward collision of the Indian plate and westward subduction of the Pacific and Philippine Sea plates beneath Asia. The Indian-Eurasian continental collision has much more extensive impact than the Pacific and Philippine oceanic plates subducting beneath Eurasia, but the contractional strain caused by the Philippine subduction zone is the largest.
References
Molnar, P. & Tapponnier, P. Cenozoic tectonics of Asia: effects of a continental collision: features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision. Science 189, 419–426 (1975).
Shen, Z. K. et al. Contemporary crustal deformation in east Asia constrained by global positioning system measurements. J. Geophys. Res. 105, 5721–5734 (2000).
Yin, A. Cenozoic tectonic evolution of Asia: A preliminary synthesis. Tectonophysics 488, 293–325 (2010).
Hall, R. Cenozoic geological and plate tectonic evolution of SE Asia and the SW pacific: computer-based reconstructions, model and animations. J. Asian Earth Sci. 20, 353–431 (2002).
Huang, J. & Zhao, D. High-resolution mantle tomography of China and surrounding regions. J. Geophys. Res. 111, B09305 (2006).
Molnar, P. Continental tectonics in the aftermath of plate tectonics. Nature 335, 131–137 (1988).
Zhao, D., Maruyama, S. & Omori, S. Mantle dynamics of Western Pacific and East Asia: insight from seismic tomography and mineral physics. Gondwana Res. 11, 120–131 (2007).
Yin, A. Mode of Cenozoic east-west extension in Tibet suggesting a common origin of rifts in Asia during the Indo-Asian collision. J. Geophys. Res. 105, 21745–21759 (2000).
Jolivet, L., Davy, P. & Cobbold, P. Right-lateral shear along the northwest Pacific margin and the India-Eurasia collision. Tectonics 9, 1409–1419 (1990).
Deng, Q.D., Zhang, Y.M., Xu, G.L. & Fan, F.T. On the tectonic stress field in China and its relation to plate movement. Seismology and Geology, 1, 11–22 (in Chinese with an English abstract) (1979).
Northrup, C. J., Royden, L. H. & Burchfiel, B. C. Motion of the Pacific plate relative to Eurasia and its potential relation to Cenozoic extension along the eastern margin of Eurasia. Geology 23, 719–722 (1995).
Michel, G. W., Becker, M., Angermann, D., Reigber, C. & Reinhart, E. Crustal motion in E-and SE-Asia from GPS measurements. Earth Planets Space 52, 713–720 (2000).
Zubovich, A. V. et al. GPS velocity field for the Tien Shan and surrounding regions. Tectonics 29, TC6014 (2010).
Liang, S. et al. Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements. J. Geophys. Res. Solid Earth 118, 1–11 (2013).
Kreemer, C., Blewitt, G. & Klein, E. C. A geodetic plate motion and Global Strain Rate Model. Geochem. Geophys. Geosyst. 15, 3849–3889 (2015).
Zheng, G. et al. Crustal deformation in the India-Eurasia collision zone from 25 years of GPS measurements. J. Geophys. Res. Solid Earth, 122 (2017).
Calais, E., Dong, L., Wang, M., Shen, Z. & Vergnolle, M. Continental deformation in Asia from a combined GPS solution. Geophys. Res. Lett. 33, L24319 (2006).
Herring, T. A., King, R. W. & McClusky, S. C. GAMIT Reference Manual (Massachusetts Institute of Technology 2010).
Herring, T. A., King, R. W. & McClusky, S. C. GLOBK Reference Manual (Massachusetts Institute of Technology 2010).
Altamimi, X., Collilieux, X. & Metivier, L. ITRF2008: an improved solution of the International Terrestrial Reference Frame. J. Geod. 85, 457–473 (2011).
Williams, S. D. P. CATS: GPS coordinate time series analysis software. GPS Solutions 12, 147–153 (2008).
Vernant, P. et al. Clockwise rotation of the Brahmaputra Valley relative to India: Tectonic convergence in the eastern Himalaya, Naga Hills, and Shillong Plateau. J. Geophys. Res. Solid Earth 119, 6558–6571 (2014).
Tape, C., Muse, P., Simons, M., Dong, D. & Webb, F. Multiscale estimation of GPS velocity fields. Geophys. J. Int. 179, 945–971 (2009).
Zhang, Y. G. et al. Contemporary deformation of the North China plain from Global Positioning System data. Geophys. Res. Lett. 45, 1851–1859 (2018).
Gan, W. et al. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements. J. Geophys. Res. 112, B08416 (2007).
Ge, W.-P., Molnar, P., Shen, Z.-K. & Li, Q. Present-day crustal thinning in the southern and northern Tibetan Plateau revealed by GPS measurements. Geophys. Res. Lett. 42, 5227–5235 (2015).
Zhang, S.X. & Ning, J.Y. Interaction between Philippine Sea plate (PH) and Eurasia (EU) plate and its influence on the movement eastern Asia. Chinese Journal of Geophysics, 45, 188–197 (in Chinese with an English abstract) (2002).
Kubo, A. & Fukuyama, E. Stress field along the Ryukyu Arc and the Okinawa Trough inferred from moment tensors of shallow earthquakes. Earth Planet. Sci. Lett. 210, 305–316 (2003).
Sibuet, J. C. et al. Okinawa trough backarc basin: early tectonic and magmatic evolution. J. Geophys. Res. 103, 30245–30267 (1998).
Lin, S. C. & Kuo, B. Y. Trench-parallel flow in the southern Ryukyu subduction system: effects of progressive rifting of the overriding plate. J. Geophys. Res. Solid Earth 118, 302–315 (2013).
Lee, C.-S., Shor, G. G. Jr., Bibee, L. D., Lee, R. S. & Hilde, T. W. C. Okinawa trough: origin of a back-arc basin. Marine Geology 35, 219–241 (1980).
Xie, F.R., Cui, X.F., Zhao, J.T., Chen, Q.C. & Li, H. Regional division of the recent tectonic stress field in China and adjacent areas. Chinese Journal of Geophysics, 47, 654–662 (in Chinese with an English abstract) (2004).
Yin, A., Yu, X., Shen, Z.-K. & Liu Z.J. A possible seismic gap and high earthquake hazard in the North China Basin. Geology, 4 (2014).
Li, C., Hilst, R. D., Meltzer, A. S. & Engdahl, E. R. Subduction of the Indian lithosphere beneath the Tibetan plateau and Burma. Earth Planet. Sci. Lett. 274, 157–168 (2008).
Wei, W., Xu, J., Zhao, D. & Shi, Y. East Asia mantle tomography: new insight into plate subduction and intraplate volcanism. J. Asian Earth Sci. 60 (2012).
Loveless, J. P. & Meade, B. J. Partitioning of localized and diffuse deformation in the tibetan plateau from joint inversions of geologic and geodetic observations. Earth Planet. Sci. Lett. 303, 1–24 (2011).
Loveless, J. P. & Meade, B. J. Geodetic imaging of plate motions, slip rates, and partitioning of deformation in Japan. J. Geophys. Res. 115, B02410 (2010).
Vergnolle, M., Calais, E. & Dong, L. Dynamics of continental deformation in Asia. J. Geophys. Res. 112, B11403 (2007).
Heki, K. & Mitsui, Y. Accelerated pacific plate subduction following interplate thrust earthquakes at the Japan trench. Earth Planet. Sci. Lett. 363, 44–49 (2013).
Wang, K., Hu, Y. & He, J. Deformation cycles of subduction earthquakes in a viscoelastic earth. Nature 484, 327–332 (2012).
Liu, X., Zhao, D. P., Li, S. Z. & Wei, W. Age of the subducting Pacific slab beneath East Asia and its geodynamic implications. Earth Planet. Sci. Lett. 464, 166–174 (2017).
Ma, J. C. et al. P-wave tomography of Northeast Asia: constraints on the western Pacific plate subduction and mantle dynamics. Physics of the Earth and Planetary Interiors 274, 105–126 (2018).
Koulakov, I. et al. Slab interactions in the Taiwan region based on the P- and S-velocity distributions in the upper mantle. J. Asian Earth Sci. 79, 53–64 (2014).
Thatcher, W. How the continents deform: the evidence from tectonic geodesy. Annu. Rev. Earth Planet Sci. 37, 237–262 (2009).
Bird, P. An updated digital model of plate boundaries. Geochem. Geophys. Geosyst. 4, 1027 (2003).
Bettinelli, P. et al. Plate motion of India and interseismic strain in the Nepal Himalaya from GPS and DORIS measurements. J. Geod. 80, 567–589 (2006).
Simons, W. J. F. et al. A decade of GPS in southeast Asia: resolving Sundaland motion and boundaries. J. Geophys. Res. 112, 623–626 (2007).
Banerjee, P., Bürgmann, R., Nagarajan, B. & Apel, E. Intraplate deformation of the Indian subcontinent. Geophys. Res. Lett. 35, 7–12 (2008).
Maurin, T., Masson, F., Rangin, C., Min, U. T. & Collard, P. First GPS results in northern Myanmar: constant and localized slip rate along the Sagaing fault. Geology 38, 591–594 (2010).
Mohadjer, S. et al. Partitioning of India-Eurasia convergence in the Pamir-Hindu Kush from GPS measurements. Geophys. Res. Lett. 37, 90–98 (2010).
Ashurkov, S. V. et al. GPS geodetic constraints on the kinematics of the Amurian plate. Russian Geology and Geophysics 52, 239–249 (2011).
Shestakov, N. V. et al. Present tectonics of the southeast of Russia as seen from GPS observations. Geophys. J. Int. 184, 529–540 (2011).
Lukhnev, A. V., San’Kov, V. A., Miroshnichenko, A. I., Ashurkov, S. V. & Calais, E. GPS rotation and strain rates in the Baikal–Mongolia region. Russian Geology and Geophysics 51, 785–793 (2010).
Ader, T. et al. Convergence rate across the Nepal Himalaya and interseismic coupling on the main Himalayan thrust: implications for seismic hazard. J. Geophys. Res. Solid Earth 117, 398–399 (2012).
Duong, N. A. et al. Contemporary horizontal crustal movement estimation for northwestern Vietnam inferred from repeated GPS measurements. Earth Planets Space 65, 1399–1410 (2013).
Ischuk, A. et al. Kinematics of the Pamir and Hindu Kush regions from GPS geodesy. J. Geophys. Res. 118, 2408–2416 (2013).
Schiffman, C., Bali, B. S., Szeliga, W. & Bilham, R. Seismic slip deficit in the Kashmir Himalaya from GPS observations. Geophys. Res. Lett. 40, 5642–5645 (2013).
Gahalaut, V. K. et al. Aseismic plate boundary in the Indo-Burmese wedge, northwest Sunda arc. Geology 41, 235–238 (2013).
Devachandra, M., Kundu, B., Catherine, J., Kumar, A. & Gahalaut, V. K. Global Positioning System (GPS) measurements of crustal deformation across the frontal eastern Himalayan syntaxis and seismic-hazard assessment. Bulletin of the Seismological Society of America 104, 1518–1524 (2014).
Jouanne, F. et al. Present-day deformation of northern Pakistan from Salt Ranges to Karakorum Ranges. J. Geophys. Res. Solid Earth 119, 2487–2503 (2014).
Acknowledgements
The GPS data used in this paper were partly obtained from the National Key Scientific Projects “Crustal Movement Observation Network of China”. We thank all participants in constructing the network and conducting the surveys. We thank Jeffrey T. Freymueller for his valuable comments and suggestions. We would like to extend our gratitude to the Geospatial Information Authority of Japan (GSI) and the Taiwan Earthquake Research Center (TEC) for providing GPS time series data. We gratefully acknowledge the editor and anonymous reviewers for their constructive suggestions that significantly improve the original manuscript. This work was financially supported by the National Key Research and Development Program of China (2017YFC1500102), the National Science Foundation of China (41874017, 41590860) and the Science for Earthquake Resilience (XH18068).
Author information
Authors and Affiliations
Contributions
M. Hao was responsible for GPS data processing, interpretation and paper writing. Y. Li conduct strain rates calculation. W. Zhuang collected and analyzed the GPS time series.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hao, M., Li, Y. & Zhuang, W. Crustal movement and strain distribution in East Asia revealed by GPS observations. Sci Rep 9, 16797 (2019). https://doi.org/10.1038/s41598-019-53306-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-53306-y
This article is cited by
-
Non-linear ground deformation detection and monitoring using time series InSAR along the coastal urban areas of Pakistan
Environmental Science and Pollution Research (2024)
-
Thermal ages of the Huatung Basin determined from seismic waveform modeling: insights into Southeast Asia’s evolution
Scientific Reports (2023)
-
Relationship between channel flow initiation and crustal viscosity in convergent settings: an analog modeling approach
International Journal of Earth Sciences (2021)