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

Block versus continuum description of lithospheric deformation in the India-Eurasia collision zone has been hotly debated over many decades. Here we apply the adapted two-dimensional (2-D) Thin Viscous Shell (TVS) approach explicitly accounting for displacement on major faults in Tibet (Altyn Tagh, Haiyuan, Kunlun, Xianshuihe, Sagaing, and Main Pamir Thrust Faults) and investigate the impact of lateral variations in depth-averaged lithospheric strength. We present a suite of dynamic models to explain the key observations from new high-resolution Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) as well as Global Navigation Satellite System (GNSS) velocities. Comparisons between calculated and observed velocity and strain rate fields indicate: (a) internal buoyancy forces from Gravitational Potential Energy (GPE) acting on a relatively weak region of high topography (~2,000 m) contribute to dilatation of high plateau and contraction on the margins; (b) a weak central Tibet (~10²¹ Pa s relative to far-field depth-averaged effective viscosity of 10²² to 10²³ Pa s) yields the observed long-wavelength eastward velocity variation away from major faults; (c) slip resistance on faults produces strain localization and clockwise rotation around the Eastern Himalayan Syntaxis (EHS). We discuss the tectonic implications for rheology of the lithosphere, distribution of geodetic strain, and partitioning of active faulting and seismicity in light of our best-fit geodynamic solutions.