Published September 2025 | 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: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

Created:
September 16, 2025
Modified:
September 16, 2025