Joint inversion of ground gravity data and satellite gravity gradients between Nepal and Bhutan: New insights on structural and seismic segmentation of the Himalayan arc

Along-strike variation in the geometry of lithospheric structures is a key control parameter for the occurrence and propagation of major interplate earthquakes in subduction and collision zones. The lateral segmentation of the Himalayan arc is now well-established from various observations, including topography, gravity anomalies, exhumation rates, and present-day seismic activity. Good knowledge of the main geometric features of these segments and their boundaries is thus the next step to improve seismic hazard assessment in this area. Following recent studies, we focus our approach on the transition zone between Nepal and Bhutan where both M > 8 earthquakes and changes in the geometry of the Indian plate have been documented. Ground gravity data sets are combined with satellite gravity gradients provided by the GOCE mission (Gravity Field and Steady-State Ocean Circulation Explorer) in a joint inversion to assess the location and the geometry of this transition. We obtain a ca. 10 km wide transition zone located at the western border of Bhutan that is aligned with the Madhupur fault in the foreland and coincides with the Dhubri-Chungthang fault zone and the Yadong-Gulu rift in Himalaya and southern Tibet, respectively. This sharp segment boundary at depth can act as a barrier to earthquake rupture propagation. It can possibly restrict the size of large earthquakes and thus reduce the occurrence probability of M > 9 earthquakes along the Main Himalayan Thrust.