Testing erosional and kinematic drivers of exhumation in the central Himalaya

Abstract

This study investigates the influence of fault geometry, kinematics, and displacement on the exhumation history of the central Himalaya using geologic mapping constraints, a new interpreted cross-section, and a suite of 176 thermochronometer ages through the Marsyangdi valley in central Nepal. Guided by the cross-section, we integrate a forward model of fold-thrust belt evolution with a 2D thermokinematic model. Model-predicted and measured thermochronometer ages were compared to evaluate the sensitivity of thermochronometer ages to the geometry and location of structures, and their rates of deformation. Results indicate 84% of the measured data can be reproduced with a largely in-sequence system of faulting where displacement occurs on the Main Central thrust (MCT) from 23-16 Ma, the Ramgarh-Munsiari thrust (RMT) from 16-7.5 Ma, the Trishuli thrust (TT) from 7.5-6 Ma and the Main Boundary thrust (MBT) from 6-3 Ma. Our cross-section solution shows the development of a duplex that initiates at 4 Ma with the TT and MBT as the roof thrust. The duplex is translated over the Main Himalayan thrust (MHT) ramp, concurrent with forward propagation of faulting in the synorogenic Siwaliks from 3 Ma to present. Notably, the 2-0.5 Ma apatite fission-track ages between the MCT and Tethyan strata 45 km north of the MCT are not reproduced by the wide range of in-sequence deformation scenarios we explored. Instead, these data are consistent with simulations that include significant displacement on two out-of-sequence (OOS) faults in the last ∼1 Ma: 10 km of OOS faulting near the MCT followed by 5 km of displacement on a fault 15 km south of the MCT. Modeled OOS thrusting both builds significant topography in the hinterland and flexurally suppresses topography in the foreland. Consequently, our preferred kinematic solution has ∼10 km of final fault motion in the Siwaliks to rebuild topography at the MBT and the active MFT from ∼0.5 Ma to the present.Modeled exhumation rates are highly variable through time, and are highest during translation of rocks over high (∼10 km) ramps and during significant changes in architecture, such as duplex formation and OOS faulting. Notably high exhumation rates include 10.5 mm/yr during the first 4 Myr of MCT motion, 4-6 mm/yr at ∼7 Ma, due to the translation over an ∼10 km high TT ramp, and 7-12 mm/yr from 4 Ma to present driven by the development of the duplex and its translation over the MHT ramp, followed by OOS thrusting. The high (7-12 mm/yr) exhumation rates documented here from 4 Ma, correlate with a time of distinct climate change, including strengthening of the monsoon, and northern hemisphere glaciation, lending support to potential climate-tectonic feedbacks.