Published September 2025 | Version Supplemental material
Journal Article Open

Ongoing India–Eurasia collision predominantly driven by Sumatra–Java slab pull

  • 1. ROR icon Southern University of Science and Technology
  • 2. ROR icon California Institute of Technology
  • 3. ROR icon Institute of Geology and Geophysics
  • 4. ROR icon University of Chinese Academy of Sciences

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

Supplementary Information

Source Data Fig. 1

Source Data Fig. 2

Source Data Fig. 3

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. 5 Predicted plate motion and the second invariant of strain rate tensor from Model 1

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

Extended Data Fig. 9 Predictions of basal drag beneath the Indian plate and line forces at plate boundaries and inside the plate for Model 1 and Model 5

Files

41561_2025_1771_MOESM1_ESM.pdf

Files (28.2 MB)

Name Size Download all
md5:cf3288d2b2e749b8b9f63bce7149e9b9
107.2 kB Download
md5:fc76504f57d5437700961bbf4084dbd8
865.3 kB Download
md5:dad854bfefb8347ffb655c55f165bab6
3.1 MB Download
md5:2af247620da2b0b701f84380982a636a
378.9 kB Preview Download
md5:467233efd735abd024da3f4598a9679c
5.6 MB Download
md5:ef79513c20eca01ee08ea34235f54d1e
7.9 MB Download
md5:66c9f630dd742649b6e2258de73b65ba
89.4 kB Download
md5:af1eebb1821f6bec5adb121e2ee85574
101.1 kB Download
md5:0df9c76b24c064116f255b19ecc7a267
2.1 MB Download
md5:b76209da1aebf021a56b7f23054371a1
50.9 kB Download
md5:00e933e5090081ac8488e3b53b1bc0f5
6.8 MB Download
md5:5fd12af215031a062e74c0bff4968eb4
1.1 MB Download

Additional details

Related works

Describes
Journal Article: https://rdcu.be/eGx1V (ReadCube)
Is supplemented by
Dataset: 10.5061/dryad.d7wm37qd4 (DOI)
Software: https://github.com/johannrudi/rhea (URL)

Funding

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

Dates

Accepted
2025-07-16
Available
2025-08-19
Published online

Caltech Custom Metadata

Caltech groups
Seismological Laboratory, Division of Geological and Planetary Sciences (GPS)
Publication Status
Published