Geometry and kinematics of the Main Himalayan Thrust and Neogene crustal exhumation in the Bhutanese Himalaya derived from inversion of multithermochronologic data

Both climatic and tectonic processes affect bedrock erosion and exhumation in convergent orogens, but determining their respective influence is difficult. A requisite first step is to quantify long‐term (~10 6  year) erosion rates within an orogen. In the Himalaya, past studies suggest long‐term erosion rates varied in space and time along the range front, resulting in numerous tectonic models to explain the observed erosion rate distribution. Here, we invert a large data set of new and existing thermochronological ages to determine both long‐term exhumation rates and the kinematics of Neogene tectonic activity in the eastern Himalaya in Bhutan. New data include 31 apatite and five zircon (U‐Th)/He ages, and 49 apatite and 16 zircon fission‐track ages along two north‐south oriented transects across the orogen in western and eastern Bhutan. Data inversion was performed using a modified version of the 3‐D thermokinematic model Pecube, with parameter ranges defined by available geochronologic, metamorphic, structural, and geophysical data. Among several important observations, our three main conclusions are as follows: (1) Thermochronologic ages do not spatially correlate with surface traces of major fault zones but appear to reflect the geometry of the underlying Main Himalayan Thrust; (2) our data are compatible with a strong tectonic influence, involving a variably dipping Main Himalayan Thrust geometry and steady state topography; and (3) erosion rates have remained constant in western Bhutan over the last ~10 Ma, while a significant decrease occurred at ~6 Ma in eastern Bhutan, which we partially attribute to convergence partitioning into uplift of the Shillong Plateau.