Variable exhumation rates and variable displacement rates: Documenting recent slowing of Himalayan shortening in western Bhutan

Nadine McQuarrie; Tobgay Tobgay; Sean P. Long; Peter W. Reiners; Michael A. Cosca

doi:10.1016/j.epsl.2013.10.045

Home/Archive/USGS

USGS · 70188812

Variable exhumation rates and variable displacement rates: Documenting recent slowing of Himalayan shortening in western Bhutan

Nadine McQuarrie·Tobgay Tobgay·Sean P. Long·Peter W. Reiners·Michael A. Cosca

View original source ↗Export metadata ↓

Abstract

We link exhumational variability in space and time to the evolving geometry of the Himalayan fold–thrust belt in western Bhutan. By combining new and published geochronologic and thermochronologic data we document the burial age, peak temperatures and complete cooling history from 20 Ma to the present over an across-strike distance of ∼125 km. These integrated cooling curves highlight windows of fast exhumation that vary spatially and temporally. We propose that pulses of fast exhumation are a result of structures that facilitate the vertical motion of material, illustrated in sequentially-restored cross sections. Due to a range of permissible geometries at depth, we explore and evaluate the impact of geometry on kinematics and rates of deformation. The linked cooling history and cross sections provide estimates of both magnitude and timing of thrust sheet displacement and highlight temporal variability in potential shortening rates. Structural and chronologic data illustrate a general north to south progression of Himalayan deformation, with emplacement of the Main Central thrust (MCT), Paro thrust and Shumar thrust by 12 to no later than 9 Ma. Two different geometries and kinematic scenarios for the Lesser Himalayan duplex are proposed. A north to south propagating duplex system requires that the southern portion of that system, south of the MCT, deformed and cooled by 9 Ma, leaving only the southernmost thrust sheets, including the Main Boundary and Main Frontal thrusts, to deform between 9 and 0 Ma. This limited post 9 Ma shortening would necessitate a marked slowdown in convergence accommodated on the Main Himalayan thrust. A two-tiered duplex system, which allows for the Paro window duplex and the southern Baxa duplex to form simultaneously, permits duplex formation and accompanying exhumation until 6 Ma. Limited cooling from ∼200 °C to the surface post 6 Ma suggests either a decrease in shortening rates from 6 to 0 Ma or that duplex formation and exhumation are temporally decoupled. Our combined cooling curves highlight that the youngest cooling ages may not mark the fastest thrusting rates or the window of fastest exhumation. Instead, temporal variations in exhumation are best viewed through identifying transients in exhumation rate. We suggest that the strongest control on exhumation magnitude and variability is fold–thrust belt geometry, particularly the locations and magnitudes of footwall ramps, which can change over 10ʼs of km distance. Balanced cross sections predict the location and magnitude of these ramps and how they vary in space and time, providing an untapped potential for testing permissible cross-section geometries and kinematics against measured cooling histories.

Explore related subjects

Earth and Planetary Science Letters Article

Source-reported bounding extent: 26.67° to 27.8° latitude; 88.83° to 90.085° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Nadine McQuarrie, Tobgay Tobgay, Sean P. Long, Peter W. Reiners, Michael A. Cosca. 2014. Variable exhumation rates and variable displacement rates: Documenting recent slowing of Himalayan shortening in western Bhutan. https://doi.org/10.1016/j.epsl.2013.10.045

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗

USGS2026

Patterns of rift basin development and the fidelity of the subsidence record: Insights through Bayesian modeling of rapid tectonic subsidence in a Rio Grande rift basin, Socorro, NM, U.S.A

Characterizing tectonic subsidence rates within depositional sequences provides direct insight into the driving mechanism(s) of accommodation in a basin. However, the temporal resolution of this record is often stymied by a lack of high-precision and high-resolution ages, which enable a more complete description of basin subsidence drivers. We explore the effect of high-precision and high-resolution ages in modeling accommodation for the Miocene La Jencia Basin of the central Rio Grande rift (RGR) and interpret driving mechanism behavior from these models (e.g., lithospheric thinning). We present a new geochronologic dataset of both laser ablation-inductively coupled plasma mass spectrometry (LA-ICPMS) and chemical abrasion-isotope dilution-thermal ionization mass spectrometry (CA-ID-TIMS) data and use these ages in Bayesian accommodation modeling. Models constrained by high-precision and high-accuracy TIMS ages yield peak tectonic subsidence rates exceeding 220 m/Myr, and an average Miocene subsidence of ∼120 m/Myr. While timing and magnitude vary, all models suggest two pulses of rapid Miocene tectonic subsidence, which we interpret to reflect basin-bounding fault movement. Prior to peak subsidence, there was an initial period of fault linkage and organization that occurred over <1–3 Myr that produced the basin-bounding La Jencia-Cerro Colorado fault zone. A comparison of published tectonic subsidence rates to those modeled here shows that while tectonic subsidence during continental rifting is highly variable, the La Jencia Basin rates appear relatively high. However, the significant difference between peak and average La Jencia Basin rates modeled here highlights the potential for underestimation of many records of tectonic subsidence due to a lack of high-precision and high-resolution age constraints. Furthermore, age data and modeling results presented here document fault movement and consequent rates of tectonic subsidence that lower-resolution data would not, providing a high-fidelity case study of continental rift basin development.

New Mexico↗

USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗