The Himalayan syntaxes, characterized by extreme rates of rock exhumation co-located with major trans-orogenic rivers, figure prominently in the debate on tectonic versus erosional forcing of exhumation. Both the mechanism and timing of rapid exhumation of the Namche Barwa massif in the eastern syntaxis remain controversial. It has been argued that coupling between crustal rock advection and surface erosion initiated in the late Miocene (8–10 Ma). Recent studies, in contrast, suggest a Quaternary onset of rapid exhumation linked to a purely tectonic mechanism. We report new multisystem detrital thermochronology data from the most proximal Neogene clastic sediments downstream of Namche Barwa and use a thermo-kinematic model constrained by new and published data to explore its exhumation history. Modeling results show that exhumation accelerated to ∼4 km/m.y. at ca. 8 Ma and to ∼9 km/m.y. after ca. 2 Ma. This three-stage history reconciles apparently contradictory evidence for early and late onset of rapid exhumation and suggests efficient coupling between tectonics and erosion since the late Miocene. Quaternary acceleration of exhumation is consistent with river-profile evolution and may be linked to a Quaternary river-capture event.
The Nanga Parbat and Namche Barwa massifs, at the respective western and eastern syntaxial terminations of the Himalaya (Fig. 1), share characteristics that have focused research into the coupling between tectonics and surface processes (Zeitler et al., 2001b; Finnegan et al., 2008; Korup et al., 2010; Koons et al., 2013; Wang et al., 2014). Both massifs show young (<10 Ma) high-grade metamorphism and partial melting (Burg et al., 1998; Zeitler et al., 2001a, 2014; Booth et al., 2009), extreme relief, and rapid erosion (Burbank et al., 1996; Finnegan et al., 2008), expressed by exceptionally young thermochronologic ages (Stewart et al., 2008; Enkelmann et al., 2011; Bracciali et al., 2016; King et al., 2016). The two largest Himalayan rivers, the Indus and the Yarlung-Tsangpo-Siang-Brahmaputra, show hairpin bends and kilometer-scale steepened knickzones as they cross these massifs (Fig. 1), sparking a debate on potential erosional controls on tectonics (Zeitler et al., 2001b; Finnegan et al., 2008; Seward and Burg, 2008; Wang et al., 2014; King et al., 2016).
Several models seek to explain these remarkable features. Purely tectonic mechanisms include range-parallel buckling in the indenter-plate corner (Burg et al., 1998), uplift driven by a geometrically stiffened bend in the subducting plate (Bendick and Ehlers, 2014), and orogen-parallel crustal transport arising from velocity and/or strain partitioning (Whipp et al., 2014). In contrast, the tectonic-aneurysm model (Zeitler et al., 2001a, 2001b; Koons et al., 2013) calls for coupling between river incision and rapid exhumation, leading to local crustal weakening and focusing rock pathways into the weakened, rapidly eroding zone. The inflowing material promotes topographic relief growth, localized exhumation, and crustal weakening, creating a positive feedback loop between tectonics and surface processes.
Besides the mechanism, the timing of rapid exhumation is also controversial in the Namche Barwa massif. Early bedrock geochronology and thermochronology studies estimated the onset of rapid exhumation at ca. 4 Ma (Burg et al., 1998; Seward and Burg, 2008), whereas more recent data (Booth et al., 2009; Zeitler et al., 2014) have suggested 8–10 Ma. Detrital thermochronology studies from the Brahmaputra Valley, the Surma Basin (Bangladesh), and the offshore Bengal Fan have proposed rapid syntaxial exhumation starting at either 4–6 Ma (Najman et al., 2019) or <3 Ma (Chirouze et al., 2013; Bracciali et al., 2016). This inconsistency may arise from downstream modification and dilution of characteristic syntaxial exhumation signals (Bracciali et al., 2016; Gemignani et al., 2018); the most robust signal is therefore expected in proximal sedimentary records. Lang et al. (2016) modeled detrital thermochronology data from the proximal Siji section (Fig. 1) to infer an onset of rapid exhumation in Namche Barwa at 5–7 Ma.
To explore the exhumation history of the Namche Barwa massif in more detail, we present new multisystem detrital thermochronology data from Neogene foreland-basin samples collected directly downstream of the syntaxis (Fig. 1) and interpret these using a thermo-kinematic inverse model.
NEW DETRITAL THERMOCHRONOLOGY DATA
We collected ten sandstone samples from three sedimentary sections close to the Siang-Brahmaputra confluence (Fig. 1). These sections are described by Govin et al. (2018), who also determined depositional ages ranging from 0.5 ± 0.3 Ma to 10.0 ± 2.0 Ma (see Tabile S1 in the Supplemental Material1). Provenance data indicate that the source region for these deposits included the Namche Barwa massif (Govin et al., 2018). Here we present new zircon fission-track (ZFT), muscovite 40Ar/39Ar (MAr), and rutile U-Pb (RUPb) data. Closure temperatures of these thermochronometers range from ∼300 °C (ZFT) to >500 °C (RUPb), depending on grain size, composition, and cooling rate (Reiners et al., 2018; Fig. S5 in the Supplemental Material). Because we target the signal from Namche Barwa, inferred to be the most rapidly exhuming part of the sediments’ source area, we employ the minimum-age approach (Galbraith, 2005) to determine the youngest detrital age populations (see the Supplemental Material for details). Sample preparation and analytical methods are reported in the Supplemental Material; single-grain ages are in Tables S2–S4 and Figures S1–S3. All ages are interpreted as cooling ages, as justified in the Supplemental Material.
A plot of the minimum ages of our samples together with literature data as a function of depositional age (Fig. 2) shows two distinct groups: for all three thermochronometers, samples with depositional ages >7.5 Ma have lag times (Bernet et al., 2001) that are >5 m.y., whereas samples with depositional ages ≤7 Ma show short lag times (∼2–3 m.y.). The latter group also shows several age inversions, where the system with lower closure temperature (ZFT) has minimum ages older than those for higher-closure-temperature systems (MAr, RUPb). Such inversions are expected at high exhumation rates in some circumstances (Reiners et al., 2018); alternatively, some of these minimum ages may be unreliable for analytical reasons (i.e., poor counting statistics for grains with low daughter-product abundance; see discussion in the Supplemental Material). We discriminate between (1) “internally consistent” samples, yielding ages ordered with respect to system closure temperatures within a sample and increasing monotonically with depositional age for the same system between samples, and (2) inconsistent samples, which do not meet these criteria.
QUANTIFYING NAMCHE BARWA EXHUMATION
The slope of the lag-time trend indicates whether exhumation rates were steady, increasing, or decreasing through time (Bernet et al., 2001). We used a Bayesian approach (Glotzbach et al., 2011) to fit single- and multi-tier linear regressions to the lag times of internally consistent minimum ages (see the Supplemental Material). The results (Fig. S4) indicate increasing exhumation before ca. 7 Ma, followed by rapid steady exhumation between ca. 7 Ma and 0.5–2.0 Ma, and probable further acceleration since 0.5–2.0 Ma indicated by the youngest two to three samples. However, the onset of rapid exhumation would precede the arrival of young grains in the sedimentary record because of (1) the time required to exhume rocks from the thermochronologic closure depth to the surface, and (2) the time required to re-equilibrate the crustal thermal structure.
To better constrain the exhumation history, we used a one-dimensional version of the thermo-kinematic code Pecube (Braun et al., 2012) to predict a time series of cooling ages resulting from step changes in exhumation rates, accounting for the effect of heat advection during exhumation. Comparison with the detrital thermochronology data is achieved through neighborhood-algorithm inversion; the model inverts for the exhumation rates and the timing of rate changes (Braun et al., 2012). Inversions use either the full data set or only the internally consistent ages, and incorporate uncertainties in both minimum-peak ages and depositional ages. Two-stage and three-stage exhumation scenarios were tested. A full description of the procedure and the different inversions is provided in the Supplemental Material.
Our best-fit inversion uses the internally consistent data set and implies a three-stage exhumation history for the Namche Barwa massif, with an early (8.2 ± 1.8 Ma) onset of rapid exhumation and a late (1.3 ± 0.8 Ma) acceleration (Fig. 3). Initial, intermediate, and final exhumation rates are 0.9 ± 0.4, 4.0 ± 2.0, and 8.6 ± 1.0 km/m.y., respectively. The onset of rapid exhumation at ca. 8 Ma is consistent with metamorphic pressure-temperature-time (P-T-t) paths from the Namche Barwa massif (Palin et al., 2015). Predicted exhumation rates agree with estimates from bedrock thermochronology (Seward and Burg, 2008; Zeitler et al., 2014; Bracciali et al., 2016), including those indicating a recent (<1 Ma) acceleration (King et al., 2016). The total amount of exhumation since ca. 8 Ma predicted by our model is 42 ± 26 km, 1–4 times the ∼15–20 km of exhumation since ca. 8 Ma inferred from P-T-t data (Fig. S8).
Exhumation rates in Namche Barwa prior to ca. 8 Ma are comparable with those elsewhere in the Greater Himalaya during the Neogene (e.g., Thiede and Ehlers, 2013), suggesting similar tectonic processes. P-T-t data (Fig. S8) suggest significant exhumation prior to ca. 8 Ma, captured by the initial phase of our model. In contrast, the clear evidence for accelerating exhumation at ca. 8 Ma and <2–3 m.y. lag times for all systems since that time, which we link to focused rapid exhumation in Namche Barwa, distinguish this easternmost detrital thermochronology record from those elsewhere in the Himalaya (e.g., Szulc et al., 2006; Chirouze et al., 2013). Our finding of sustained rapid exhumation since the late Miocene is consistent with previous work (Lang et al., 2016). However, inclusion of a high-temperature thermochronometer (RUPb) coupled with Pecube inversions allows us to reconstruct a more detailed three-stage exhumation history, reconciling previous apparently contrasting interpretations that emphasized either the earlier (ca. 8 Ma; Booth et al., 2009; Zeitler et al., 2014) or later (<2 Ma; Wang et al., 2014; King et al., 2016) time of exhumation-rate change.
The onset of rapid exhumation at ca. 8 Ma is consistent with the scenario envisaged by Zeitler et al. (2014). The discrepancy between the amount of post–ca. 8 Ma exhumation predicted by our data and that inferred from P-T-t data, for all but our lowest predicted exhumation rates, implies lateral inflow of mid-crustal material, consistent with the tectonic-aneurysm model (Zeitler et al., 2001a, 2001b; Koons et al., 2013). Thus, efficient coupling between crustal rock advection and surface erosion may have initiated at ca. 8 Ma, requiring the existence of a large throughgoing river system at that time. Whereas sedimentary provenance data record a drainage connection between the Yarlung-Tsangpo and the Brahmaputra since the early Miocene (ca. 18 Ma; Lang and Huntington, 2014; Bracciali et al., 2015; Blum et al., 2018), it is unclear when the drainage pathway through the Namche Barwa massif via the Siang was established (Govin et al., 2018).
The trigger for rapid exhumation in Namche Barwa remains debated. It could have been initiated by indenter-corner dynamics (Burg et al., 1998; Bendick and Ehlers, 2014), with coupling between river incision and rapid exhumation developing subsequently, or it could have resulted from capture of the Yarlung-Tsangpo by the Siang shortly before ca. 8 Ma. The latter scenario is consistent with river-incision patterns upstream of Namche Barwa, which have been interpreted to record a wave of incision migrating upstream since ca. 10 Ma (Schmidt et al., 2015). However, that scenario requires prior Yarlung-Tsangpo drainage to the foreland via another, as yet unconstrained, pathway, a prediction that may be tested by provenance analysis of proximal foreland sediment records from candidate fossil trans-orogenic river systems.
Quaternary uplift of the Namche Barwa massif has been inferred from a thick wedge of post–2.6 Ma alluvium preserved immediately upstream (Wang et al., 2014; Fig. 4). This ponded sediment implies that rock uplift temporarily outpaced river incision, steepening the Siang River profile downstream. Quaternary capture of the Parlung River by the Yarlung-Tsangpo-Siang, as suggested by thermochronology (Seward and Burg, 2008; Zeitler et al., 2014) and provenance data (Lang and Huntington, 2014; Govin et al., 2018), would have increased erosional power in the gorge downstream of the capture point. In turn, this may have strengthened the feedback loop and triggered enhanced uplift and exhumation of Namche Barwa. River profiles provide insight into this possibility. The modern Yarlung-Tsangpo-Siang and Indus River profiles differ (Korup et al., 2010; Fig. 4), even though both flow through rapidly exhuming syntaxial massifs (Burbank et al., 1996; Zeitler et al., 2001b; Finnegan et al., 2008; Korup et al., 2010). The inferred pre-Quaternary profile of the Yarlung-Tsangpo-Siang resembles the modern Indus profile, with a more subdued knickzone across Namche Barwa and a morphologic plateau edge located farther upstream. Modeling of river incision (Koons et al., 2013) shows that the differences in modern river profiles can be induced by differing rock-uplift rates in the syntaxial massifs of ∼5 and ∼10 mm/yr (Fig. 4), consistent with the recent acceleration our data imply.
Our new data and modeling reveal a three-stage exhumation history for Namche Barwa, reconciling previous studies focusing on either an early or a late onset of rapid exhumation. Our results suggest that coupling between crustal rock advection and surface erosion initiated in the late Miocene and strengthened during the Quaternary. They suggest a potential role for river capture events in initiating and strengthening tectonic-erosion couplings in a tectonic aneurysm.
This work is part of the Marie Curie Initial Training Network “Investigating Tectonic-Erosion-Climate Couplings” (iTECC) funded by the European Union Research Executive Agency (grant agreement 316966). It was supported by the UK Natural Environment Research Council (NERC; IP-1500–1114). C. Mark is funded by Science Foundation Ireland (18/SIRG/5559). E. Enkelmann, K. Lang, and P. Zeitler provided constructive reviews. Data for this study are available in the Supplemental Material and in the geochron.org database.