The Dynamics of the India‐Eurasia Collision: Faulted Viscous Continuum Models Constrained by High‐Resolution Sentinel‐1 InSAR and GNSS Velocities

The distribution and magnitude of forces driving lithospheric deformation in the India-Eurasia collision zone have been debated over many decades. Here we test a two-dimensional (2-D) Thin Viscous Shell approach that has been adapted to explicitly account for displacement on major faults and investigate the impact of lateral variations in depth-averaged lithospheric strength. We present a suite of dynamic models to explain the key features from new high-resolution Sentinel-1 Interferometric Synthetic Aperture Radar as well as Global Navigation Satellite System velocities. Comparisons between calculated and geodetically observed velocity and strain rate fields indicate: (a) internal buoyancy forces from Gravitational Potential Energy acting on a relatively weak region of highest topography (>2,000 m) contribute to dilatation of the high plateau and contraction on the margins; (b) a weak central Tibetan Plateau (∼10²¹ Pa s compared to far-field depth-averaged effective viscosity of at least 10²²–10²³ Pa s) is required to explain the observed long-wavelength eastward velocity variation; (c) localized displacement on fault systems enables strain concentration and clockwise rotation around the Eastern Himalayan Syntaxis. We discuss the tectonic implications for rheology of the lithosphere, distribution of geodetic strain, and partitioning of active faulting and seismicity.