Ongoing India–Eurasia collision predominantly driven by Sumatra–Java slab pull
Abstract
Continued India–Eurasia convergence since the early Palaeogene has led to the formation of the Tibetan Plateau. Yet the primary driving mechanisms of this protracted convergence remain debated, limiting our understanding of continental collision dynamics. Here we provide a holistic quantification of various driving forces to this convergence by integrating high-resolution, plate-boundary-resolving global convection models with observational constraints. Whereas different forces can produce the observed plate motion, we show that the primary driving force can be definitively constrained when Indo-Australian intraplate stress and strain rates are used as constraints in addition to plate motions. Specifically, we identify that the position of the transition in stress orientation within the Indo-Australian plate is highly sensitive to the relative strength of plate-boundary forces. When the plate motion and this stress-orientation transition are fit simultaneously, we find slab pull from Sumatra–Java subduction is the predominant driving force of India–Eurasia convergence with continental collision exerting an overall resisting force and rule out mantle basal drag playing more than a secondary role. We suggest slab pull from adjacent subduction zones has been the primary driver of the uplift of the Tibetan Plateau since its onset and so this may be an exceptional event in Earth's history.
Copyright and License
© 2025 Springer Nature Limited. Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
Acknowledgement
This study is supported by the National Key R&D Program of China 2023YFF0803200 (J.H.) and the National Natural Science Foundation of China (NSFC) 42174106, 92155307 and 92355302 (J.H.). Q.Z. is supported by the National Key R&D Program of China 2023YFF0803404.
Data Availability
The World Stress Map data can be accessed at www.world-stress-map.org/. The focal mechanism solution data from International Seismological Centre can be accessed at www.isc.ac.uk. Data generated for this study are available on Dryad at https://doi.org/10.5061/dryad.d7wm37qd4 (ref. 62). Source data are provided with this paper.
Code Availability
The adaptive nonlinear Stokes solver (Rhea) and scripts related to force calculations are available on Github at https://github.com/johannrudi/rhea. The code and software used to make the figures can be downloaded at www.soest.hawaii.edu/gmt/ and www.paraview.org/. The Matlab script (stress2grid) used to estimate mean SHmax orientations on a regular grid can be downloaded at www.world-stress-map.org/software.
Supplemental Material
Source Data Extended Data Fig. 1
Source Data Extended Data Fig. 3
Source Data Extended Data Fig. 4
Source Data Extended Data Fig. 5
Source Data Extended Data Fig. 6
Source Data Extended Data Fig. 7
Source Data Extended Data Fig. 8
Source Data Extended Data Fig. 9
Extended Data Fig. 1 SHmax records and mean orientations for the Indo-Australian plate
Extended Data Fig. 2 Viscosity fields in geodynamic models
Extended Data Fig. 3 Diagram of the calculation method for plate-boundary forces
Extended Data Fig. 4 Effective ridge push for the Central and Southeast Indian Ridge
Extended Data Fig. 6 Predicted velocity and viscosity fields at 300 km depth for Model 1 and Model 6
Extended Data Fig. 7 Predicted plate motions for models
Extended Data Fig. 8 Predicted second invariant of the strain rate tensor
Files
Name | Size | Download all |
---|---|---|
md5:66c9f630dd742649b6e2258de73b65ba
|
89.4 kB | Download |
md5:0df9c76b24c064116f255b19ecc7a267
|
2.1 MB | Download |
md5:cf3288d2b2e749b8b9f63bce7149e9b9
|
107.2 kB | Download |
md5:2af247620da2b0b701f84380982a636a
|
378.9 kB | Preview Download |
md5:fc76504f57d5437700961bbf4084dbd8
|
865.3 kB | Download |
md5:467233efd735abd024da3f4598a9679c
|
5.6 MB | Download |
md5:5fd12af215031a062e74c0bff4968eb4
|
1.1 MB | Download |
md5:ef79513c20eca01ee08ea34235f54d1e
|
7.9 MB | Download |
md5:dad854bfefb8347ffb655c55f165bab6
|
3.1 MB | Download |
md5:af1eebb1821f6bec5adb121e2ee85574
|
101.1 kB | Download |
md5:b76209da1aebf021a56b7f23054371a1
|
50.9 kB | Download |
md5:00e933e5090081ac8488e3b53b1bc0f5
|
6.8 MB | Download |
Additional details
- Ministry of Science and Technology of the People's Republic of China
- National Key R&D Program of China 2023YFF0803200
- National Natural Science Foundation of China
- 42174106
- National Natural Science Foundation of China
- 92155307
- National Natural Science Foundation of China
- 92355302
- Ministry of Science and Technology of the People's Republic of China
- National Key R&D Program of China 2023YFF0803404
- Accepted
-
2025-07-16
- Available
-
2025-08-19Published online
- Caltech groups
- Seismological Laboratory, Division of Geological and Planetary Sciences (GPS)
- Publication Status
- Published