Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation

Constraints on crustal compositional architecture across the North China-Altaids transition and implications for craton margin reworking

The phase of secular evolution of continents is manifested as the degree of compositional differentiation, modification and maturation of continental crusts, which is vital in understanding the mechanism of continental evolution but is difficult to quantify. Here we use integrated passive- and active-source seismic profiling to conduct joint analysis and inversion and derive Vs and Vp/Vs section models across the North China Craton (NCC) to southeastern Altaids boundary zone that bears a tectonic transition from a reworked ancient craton margin to a Phanerozoic accretionary orogen. We systematically exploited the imaged multiple physical properties as precise and delicate proxies to constrain the compositional architecture in the crust across this important tectonic transition subject to various crustal evolutional phases. Our Vs and Vp/Vs imaging, together with the existing isotopic data, characterizes the Yin Shan-Yan Shan belt as the northern NCC margin with layered homogeneous compositions that point to an evolved crust. However, the lower-crustal low-Vs/high-Vp/Vs signature that overlaps the shallowly dipping to horizontal reflective fabrics suggests that the crust of the northern NCC margin has undergone considerable reworking through lower-crustal-stretching-assisted melt migration and mixing since the late Paleozoic to Mesozoic eras. The process probably involved crust-mantle interaction and thus resulted in a compositionally modified ancient crustal basement. On the contrary, the southeastern Altaids domain manifests crustal complexity in compositions and structures inferred to be indicative of a juvenile crust of the Phanerozoic accretionary orogen. Our results provide deep physical-property constraints that shed new light on the crustal evolution of a complex craton margin.

New constraints on Cenozoic subduction between India and Tibet

The type of lithosphere subducted between India and Tibet since the Paleocene remains controversial; it has been suggested to be either entirely continental, oceanic, or a mixture of the two. As the subduction history of this lost lithosphere strongly shaped Tibetan intraplate tectonism, we attempt to further constrain its nature and density structure with numerical models that aim to reproduce the observed history of magmatism and crustal thickening in addition to present-day plateau properties between 83°E and 88°E. By matching time-evolving geological patterns, here we show that Tibetan tectonism away from the Himalayan syntaxis is consistent with the initial indentation of a craton-like terrane at 55 ± 5 Ma, followed by a buoyant tectonic plate with a thin crust, e.g., a broad continental margin (Himalandia). This new geodynamic scenario can explain the seemingly contradictory observations that had led to competing hypotheses like the subduction of Greater India versus largely oceanic subduction prior to Indian indentation.

Subduction Evolution Controlled Himalayan Orogenesis: Implications from 3-D Subduction Modeling

Himalayan orogenesis remains enigmatic in terms of Tibetan Plateau geodynamics originating from the Cenozoic India–Eurasian continental collision. India underthrusts below Tibet to the Yarlung–Tsangpo suture, which has been identified as the northernmost boundary for underplating. However, the way in which the historical evolution of continental subduction induces plateau uplift and the way it controls the variation in uplift between outboard and inboard areas is still unclear. To interpret the evolutionary mechanisms involved in the Himalayan growth history, we constructed different 3-D dynamic models at important stages to address these questions related to the formation of the Himalayas on the basis of paleoenthalpy evidence encoded in fossil leaves from recently documented assemblages in southern Tibet. The results show that (1) the effect of crustal thickening was the predominant factor in the early evolution from the Paleocene to the early Eocene, which resulted in a moderate growth rate. (2) The consecutive slab break-off eastward from the western syntaxis and the associated slab rebound significantly accelerated orogenesis from the late Eocene to the Oligocene. The upwelling asthenospheric flow was a key control of increasing crustal buoyancy, which resulted in the fastest growth of the Himalayas during the early Miocene. (3) Thereafter, the gradually enhanced monsoon and surface erosion during accompanying the increasing mountain height resulted in a slowdown of the orogenic rate, which counterbalanced the buoyant force produced by asthenospheric flow driving continuous Himalayan growth.

Topography of mountain belts controlled by rheology and surface processes

It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics^1-4. However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process^5,6 and mantle-scale tectonic model⁷. End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state⁸, strength controlled (Bm ≈ 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.

Constraining the Timing of Evolution of Shear Zones in Two Collisional Orogens: Fusing Structural Geology and Geochronology

In recent decades, constraining the timing of shear activity has been one of the main topics of research about the tectono-metamorphic evolution of orogenic belts. We present a review of a combined structural and geochronological approach to two major ductile regional shear zones, in two collisional orogens: the first one affecting the Variscan basement in northern Sardinia (Italy) and the External Crystalline Massifs of the Alps (East Variscan Shear Zone; EVSZ), and the second one deforming the medium- to high-grade rocks of the metamorphic core of the Himalaya (High Himalayan Discontinuity). High-resolution, texturally and chemically controlled monazite geochronology applied in separated shear zones of the Variscan belt allowed recognizing a similar timing of activity ranging between c. 340–330 and 300 Ma. This approach led to a better understanding of the evolution of the EVSZ, supporting a model where several branches were active according to a growth by linkage model. Following a similar approach, in situ U-Th-Pb analysis of monazite constrained the timing of top-to-the-S/SW shearing of a regional-scale High Himalayan Discontinuity in the Himalayan belt to between c. 28 Ma and 17 Ma. Earlier exhumation of the hanging wall was triggered by shear zone activity, whereas at the same time, the footwall was still experiencing burial with increasing P-T conditions. The timing of shearing of this shear zone fits with an in-sequence shearing tectonic model for the exhumation of the Himalayan mid-crust.

Indian plate paleogeography, subduction and horizontal underthrusting below Tibet: paradoxes, controversies and opportunities

The India–Asia collision zone is the archetype to calibrate geological responses to continent–continent collision, but hosts a paradox: there is no orogen-wide geological record of oceanic subduction after initial collision around 60–55 Ma, yet thousands of kilometers of post-collisional subduction occurred before the arrival of unsubductable continental lithosphere that currently horizontally underlies Tibet. Kinematically restoring incipient horizontal underthrusting accurately predicts geologically estimated diachronous slab break-off, unlocking the Miocene of Himalaya–Tibet as a natural laboratory for unsubductable lithosphere convergence. Additionally, three endmember paleogeographic scenarios exist with different predictions for the nature of post-collisional subducted lithosphere but each is defended and challenged based on similar data types. This paper attempts at breaking through this impasse by identifying how the three paleogeographic scenarios each challenge paradigms in geodynamics, orogenesis, magmatism or paleogeographic reconstruction and identify opportunities for methodological advances in paleomagnetism, sediment provenance analysis, and seismology to conclusively constrain Greater Indian paleogeography.

Ongoing formation of felsic lower crustal channel by relamination in Zagros collision zone revealed from regional tomography

Complex interaction of rheologically contrasting layers within the lithosphere during the collision of continental plates leads to active faulting, which represents a serious hazard to the population and infrastructure. One of the collision scenarios presumes the existence of a middle-lower crustal channel composed of subducted silicic upper crustal rocks, which is thought to exist in the Tibetan-Himalayan system. Based on the results of seismic tomography, we argue that a similar mechanism of crustal channeling takes place beneath the Zagros mountain system in southwestern Iran. The 3D seismic velocity model reveals an inverted crustal architecture of the collision zone, in which the low-velocity felsic (granitic and sedimentary) upper crustal rocks of the Arabian plate form a seismically inactive lower crustal channel below the higher-velocity mafic (basaltic) middle-upper crustal layer of the Iranian crust. Based on existing numerical models, we suggest that the formation of the felsic channel is likely governed by separation (delamination) of the weak felsic upper crust of the subducting Arabian lithosphere and its ductile underplating under rheologically stronger upper-middle crust of the Iranian plate.

Influence of the Thickness of the Overriding Plate on Convergence Zone Dynamics

The important role played by the upper plate in convergence zones dynamics has long been underestimated but is now more and more emphasized. However, the influence of its thickness and/or strength on orogenic systems evolution remains largely unknown. Here we present results from 3D thermo-mechanical numerical simulations of convergence zones (including oceanic subduction followed by continental subduction/collision), in which we vary the rheological profile of the overriding plate (OP). For this, we systematically modify the crustal thickness of the overriding lithosphere and the temperature at the Moho to obtain a thermal thickness of the overriding lithosphere ranging from 80 to 180 km. While all models share a common global evolution (i.e., slab sinking, interaction between slab and the 660 km discontinuity, continental subduction/collision, and slab breakoff), they also highlight first-order differences arising from the variations in the OP strength (thermal thickness). With a thin/weak OP, slab rollback is favored, the slab dip is low, the mantle flow above the slab is vigorous, and the trench migrates at a high rate compared to a thick/strong OP. In addition, slab breakoff and back-arc basin formation events occur significantly earlier than in models involving a thick OP. Our models therefore highlight the major role played by the thickness/strength of the OP on convergence zone dynamics and illustrate its influence in a quantitative way.

Multi-terrane structure controls the contrasting lithospheric evolution beneath the western and central-eastern Tibetan plateau

The Tibetan plateau is manifested by contrasting along-strike lithospheric structures, but its formation mechanism and the relationship with the heterogeneous multi-terrane configuration is a challenging problem. Here we conduct systematic numerical modeling to explore the roles of width, density, and rheological properties of the multiple terranes in the lithospheric evolution of the Tibetan plateau, which reveals two distinct collision modes. In Mode-I, the lithospheric mantles of both the strong and weak terranes in the Tibetan plate are completely detached, followed by the underthrusting of Indian lithosphere beneath the whole plateau. Alternatively, Mode-II is characterized by full detachment of the weak terranes, but (partial) residue of the strong terranes during collision. These two contrasting modes, broadly consistent with the lithospheric structures of western and central–eastern Tibetan plateau, respectively, are strongly dependent on the along-strike variation of the width of the strong Lhasa–Qiangtang terranes.

The Tibetan plateau is manifested by contrasting along-strike lithospheric structures, but the correlation with multi-terrane configuration remains challenging. Here, the authors show the crucial roles of the original geometric shape of accreted terranes in regulating the lithospheric evolution of Tibetan plateau.

Melting conditions in the modern Tibetan crust since the Miocene

Abundant granitic rocks exposed in ancient mountain belts suggest that crustal melting plays a major role in orogenic processes. However, complex field relations and superposition of multiple tectonic events make it difficult to determine the role of melting in orogenesis. In contrast, geophysical measurements image present-day crustal conditions but cannot discriminate between partial melt and aqueous fluids. Here we connect pressure–temperature paths of Himalayan Miocene crustal rocks to the present-day conditions beneath the Tibetan plateau imaged with geophysical data. We use measurements of electrical conductivity to show that 4–16% water-rich melt is required to explain the crustal conductivity in the north-western Himalaya. In southern Tibet, higher melt fractions >30% reflect a crust that is either fluid-enriched (+1% H₂O) or hotter (+100 °C) compared to the Miocene crust. These melt fractions are high enough for the partially molten rocks to be significantly weaker than the solid crust.

Crustal melting may play a fundamental role in orogenic processes, but quantifying crustal melt remains difficult. Here, the authors combine pressure-temperature paths, electrical conductivity and geophysical data to elucidate the melting conditions in Tibet since the Miocene.

Oblique convergence and strain partitioning in the outer deformation front of NE Himalaya

Himalayan-Tibetan orogeny has considered as a natural black box in the context of geodynamic evolution and tectonic complexity. The eastward extrusion model of Tibetan crust contradicts with the oblique convergence model in the NE-Himalaya (Bhutan/Arunachal region), where the overall convergence rate accommodated in the Himalaya is about 20–25% less than that in the neighbouring central Himalaya and Eastern Himalayan syntaxis (EHS). We propose that instead of partitioning in the backarc, the NE-Himalaya has developed an active sliver along the Assam-Brahmaputra valley in the outer deformation front, in order to accommodate the deficiency in long-term plate convergence between Himalaya and southern Tibet. We argue that the strong eastward extrusion of Tibetan crust along NE-Himalaya is the main driving force for the unusual development of the Assam-Brahmaputra sliver. This new hypothesis can explain active convergence along EHS, low convergence and subdued topography in Bhutan and Arunachal Himalaya, kinematic and space-problem of Indo-Burmese wedge, and finally solves the contradiction between Tibetan extrusion and oblique convergence model of the HimalayanTibetan orogeny.

Slow Long-Term Exhumation of the West Central Andean Plate Boundary, Chile

We present a regional analysis of new low-temperature thermochronometer ages from the Central Andean fore arc to provide insights into the exhumation history of the western Andean margin. To derive exhumation rates over 10 million-year timescales, 38 new apatite and zircon (U-Th)/He ages were analyzed along six ~500-km long near-equal-elevation, coast parallel, transects in the Coastal Cordillera (CC) and higher-elevation Precordillera (PC) of the northern Chilean Andes between latitudes 18.5°S and 22.5°S. These transects were augmented with age-elevation profiles where possible. Results are synthesized with previously published thermochronometric data, corroborating a previously observed trenchward increase in cooling ages in Peru and northern Chile. One-dimensional thermal-kinematic modeling of all available multichronometer equal-elevation samples reveals mean exhumation rates of <0.2 km/Myr since ~50 Ma in the PC and ~100 Ma in the CC. Regression of pseudovertical age-elevation transects in the CC yields comparable rates of ~0.05 to ~0.12 km/Myr between ~40 and 80 Ma. Differences between the long-term mean 1-D rates and shorter-term age-elevation-derived rates indicate low variability in the exhumation history. Modeling results suggest similar background exhumation rates in the CC and PC; younger ages in the PC are largely a function of increased heat flow and consequently an elevated geothermal gradient near the arc. Slow exhumation rates are suggestive of semiarid conditions across the region since at least the Eocene and deformation and development of the Andean fore arc around this time.

Crustal rheology controls on the Tibetan plateau formation during India-Asia convergence

The formation of the Tibetan plateau during the India-Asia collision remains an outstanding issue. Proposed models mostly focus on the different styles of Tibetan crustal deformation, yet these do not readily explain the observed variation of deformation and deep structures along the collisional zone. Here we use three-dimensional numerical models to evaluate the effects of crustal rheology on the formation of the Himalayan-Tibetan orogenic system. During convergence, a weaker Asian crust allows strain far north within the upper plate, where a wide continental plateau forms behind the orogeny. In contrast, a stronger Asian crust suppresses the plateau formation, while the orogeny accommodates most of the shortening. The stronger Asian lithosphere is also forced beneath the Indian lithosphere, forming a reversed-polarity underthrusting. Our results demonstrate that the observed variations in lithosphere deformation and structures along the India-Asia collision zone are primarily controlled by the strength heterogeneity of the Asian continental crust.

Vertical crustal motions across Eastern Tibet revealed by topography-dependent seismic tomography

Using a topography-dependent tomographic scheme, the seismic velocity structure of the Eastern Tibetan Plateau, including the uplifted Longmenshan (LMS) orogenic belt, is accurately imaged in spite of the extreme topographic relief in the LMS region and thick sedimentary covers in the neighbouring Sichuan Basin. The obtained image shows a high-resolution upper crustal structure on a 500 km-long profile that is perpendicular to the LMS. The image clearly shows that the crystalline basement was uplifted within the LMS orogenic belt, and that the neighbouring Songpan-Ganzi Terrane was covered by a thick flysch belt, with evidence of near-surface thrust faults caused by convergence between Eastern Tibet and the Sichuan Basin. The indication that the lower crust beneath the LMS was folded and pushed upwards and the upper crust was removed by exhumation, supports the concept of a lower crustal channel flow beneath Eastern Tibet. The image also reveals that the destructive Wenchuan earthquake of year 2008 occurred in the upper crust, directly at the structural discontinuity between Eastern Tibet Plateau and the Sichuan Basin.

Pliocene-Quaternary crustal melting in central and northern Tibet and insights into crustal flow

There is considerable controversy over the nature of geophysically recognized low-velocity–high-conductivity zones (LV–HCZs) within the Tibetan crust, and their role in models for the development of the Tibetan Plateau. Here we report petrological and geochemical data on magmas erupted 4.7–0.3 Myr ago in central and northern Tibet, demonstrating that they were generated by partial melting of crustal rocks at temperatures of 700–1,050 °C and pressures of 0.5–1.5 GPa. Thus Pliocene-Quaternary melting of crustal rocks occurred at depths of 15–50 km in areas where the LV–HCZs have been recognized. This provides new petrological evidence that the LV–HCZs are sources of partial melt. It is inferred that crustal melting played a key role in triggering crustal weakening and outward crustal flow in the expansion of the Tibetan Plateau.

The role of the low velocity-high conductivity zones (LV–HCZs) in developing the Tibetan Plateau has remained controversial. Here, Wang et al. present new geochemical and petrological data that show the LV–HCZs are sources of partial melt thus giving insight into the development of the Tibetan Plateau.

Mid-Pleistocene climate transition drives net mass loss from rapidly uplifting St. Elias Mountains, Alaska

Erosion, sediment production, and routing on a tectonically active continental margin reflect both tectonic and climatic processes; partitioning the relative importance of these processes remains controversial. Gulf of Alaska contains a preserved sedimentary record of the Yakutat Terrane collision with North America. Because tectonic convergence in the coastal St. Elias orogen has been roughly constant for 6 My, variations in its eroded sediments preserved in the offshore Surveyor Fan constrain a budget of tectonic material influx, erosion, and sediment output. Seismically imaged sediment volumes calibrated with chronologies derived from Integrated Ocean Drilling Program boreholes show that erosion accelerated in response to Northern Hemisphere glacial intensification (∼ 2.7 Ma) and that the 900-km-long Surveyor Channel inception appears to correlate with this event. However, tectonic influx exceeded integrated sediment efflux over the interval 2.8-1.2 Ma. Volumetric erosion accelerated following the onset of quasi-periodic (∼ 100-ky) glacial cycles in the mid-Pleistocene climate transition (1.2-0.7 Ma). Since then, erosion and transport of material out of the orogen has outpaced tectonic influx by 50-80%. Such a rapid net mass loss explains apparent increases in exhumation rates inferred onshore from exposure dates and mapped out-of-sequence fault patterns. The 1.2-My mass budget imbalance must relax back toward equilibrium in balance with tectonic influx over the timescale of orogenic wedge response (millions of years). The St. Elias Range provides a key example of how active orogenic systems respond to transient mass fluxes, and of the possible influence of climate-driven erosive processes that diverge from equilibrium on the million-year scale.

Erosion in southern Tibet shut down at ∼10 Ma due to enhanced rock uplift within the Himalaya

Exhumation of the southern Tibetan plateau margin reflects interplay between surface and lithospheric dynamics within the Himalaya-Tibet orogen. We report thermochronometric data from a 1.2-km elevation transect within granitoids of the eastern Lhasa terrane, southern Tibet, which indicate rapid exhumation exceeding 1 km/Ma from 17-16 to 12-11 Ma followed by very slow exhumation to the present. We hypothesize that these changes in exhumation occurred in response to changes in the loci and rate of rock uplift and the resulting southward shift of the main topographic and drainage divides from within the Lhasa terrane to their current positions within the Himalaya. At ∼17 Ma, steep erosive drainage networks would have flowed across the Himalaya and greater amounts of moisture would have advected into the Lhasa terrane to drive large-scale erosional exhumation. As convergence thickened and widened the Himalaya, the orographic barrier to precipitation in southern Tibet terrane would have strengthened. Previously documented midcrustal duplexing around 10 Ma generated a zone of high rock uplift within the Himalaya. We use numerical simulations as a conceptual tool to highlight how a zone of high rock uplift could have defeated transverse drainage networks, resulting in substantial drainage reorganization. When combined with a strengthening orographic barrier to precipitation, this drainage reorganization would have driven the sharp reduction in exhumation rate we observe in southern Tibet.

Topographic Evolution and Climate Aridification during Continental Collision: Insights from Computer Simulations

How do the feedbacks between tectonics, sediment transport and climate work to shape the topographic evolution of the Earth? This question has been widely addressed via numerical models constrained with thermochronological and geomorphological data at scales ranging from local to orogenic. Here we present a novel numerical model that aims at reproducing the interaction between these processes at the continental scale. For this purpose, we combine in a single computer program: 1) a thin-sheet viscous model of continental deformation; 2) a stream-power surface-transport approach; 3) flexural isostasy allowing for the formation of large sedimentary foreland basins; and 4) an orographic precipitation model that reproduces basic climatic effects such as continentality and rain shadow. We quantify the feedbacks between these processes in a synthetic scenario inspired by the India-Asia collision and the growth of the Tibetan Plateau. We identify a feedback between erosion and crustal thickening leading locally to a <50% increase in deformation rates in places where orographic precipitation is concentrated. This climatically-enhanced deformation takes place preferentially at the upwind flank of the growing plateau, specially at the corners of the indenter (syntaxes). We hypothesize that this may provide clues for better understanding the mechanisms underlying the intriguing tectonic aneurisms documented in the Himalayas. At the continental scale, however, the overall distribution of topographic basins and ranges seems insensitive to climatic factors, despite these do have important, sometimes counterintuitive effects on the amount of sediments trapped within the continent. The dry climatic conditions that naturally develop in the interior of the continent, for example, trigger large intra-continental sediment trapping at basins similar to the Tarim Basin because they determine its endorheic/exorheic drainage. These complex climatic-drainage-tectonic interactions make the development of steady-state topography at the continental scale unlikely.

Erosion influences the seismicity of active thrust faults

Assessing seismic hazards remains one of the most challenging scientific issues in Earth sciences. Deep tectonic processes are classically considered as the only persistent mechanism driving the stress loading of active faults over a seismic cycle. Here we show via a mechanical model that erosion also significantly influences the stress loading of thrust faults at the timescale of a seismic cycle. Indeed, erosion rates of about ~0.1-20 mm yr(-1), as documented in Taiwan and in other active compressional orogens, can raise the Coulomb stress by ~0.1-10 bar on the nearby thrust faults over the inter-seismic phase. Mass transfers induced by surface processes in general, during continuous or short-lived and intense events, represent a prominent mechanism for inter-seismic stress loading of faults near the surface. Such stresses are probably sufficient to trigger shallow seismicity or promote the rupture of deep continental earthquakes up to the surface.

New geochronological constraints on the thermal and exhumation history of the Lesser and Higher Himalayan Crystalline Units in the Kullu-Kinnaur area of Himachal Pradesh (India)

New geochronological, petrological and structural data from the Beas-Sutlej area of Himachal Pradesh (India) are used to reconstruct the tectonothermal and exhumation history of this part of the Himalayan orogen. Sm-Nd garnet ages at 40.5 ± 1.3 Ma obtained on a pegmatoid from the inverse metamorphic High Himalayan Crystalline (HHC) in the Malana-Parbati area probably mark local melting during initial decompression. Ongoing exhumation in ductilely deformed leuco-gneiss is constrained by Sm-Nd garnet ages at 29 ± 1 Ma and white mica Rb-Sr ages around 24-20 Ma, while Bt Rb-Sr ages indicate a drop of regional metamorphic temperatures below 300 °C between 15 and 12 Ma. The deep Sutlej gorge exposes medium-grade paragneisses and Proterozoic orthogneisses of the Lesser Himalayan Crystalline (LHC), overthrust by the HHC along the Main Central Thrust (MCT). Mica cooling ages in the HHC are in the range of 14-11 Ma. Above the extruded wedge of the HHC, the Leo Pargil leucogranite and associated dykes intrude the Haimanta Unit (HU) below the weakly metamorphic Palaeo-Mesozoic sediments of the Tethyan Himalayas (TH). The Leo Pargil leucogranite yielded a mean Sm-Nd garnet age of 19 ± 1 Ma and Rb-Sr muscovite and biotite cooling ages between 16.4 and 11.6 Ma. Marked young extrusion of LHC units resulted in differentiated exhumation/cooling of more frontal parts of the orogen. Very young ductile deformation of LHC gneisses near Wangtu is constrained by late-kinematic pegmatite intrusions crosscutting the main mylonitic foliation. Sm-Nd garnet and Rb-Sr muscovite ages of these pegmatites range between 7.9 ± 0.9 and 5.5 ± 0.1 Ma. Published apatite FT ages down to 0.6 Ma also document accelerated diachronous sub-recent exhumation of different parts of the orogen. Together with geochronological data from the literature, the new results demonstrate that the HHC and the HU were deformed by shortening and crustal thickening during the Eohimalayan phase (Late Eocene-Oligocene), followed by a strong thermal overprint and intrusions of granitoids during the Neohimalayan Phase (Early to Middle Miocene). The LHC experienced amphibolite facies metamorphic conditions in the Late Miocene prior to extrusion between the HHC and the very low-grade Lesser Himalayan sediments. In conjunction with climate changes, young tectonic activity in this central part of the Himalayan orogen may have strongly influenced fluvial incision and erosion, and therefore, contributed to the accelerated uplift, as indicated by extensive accumulation of Late Miocene to Early Pleistocene fluviatile-lacustrine sediments in the Zanda basin, the Transhimalayan headwaters of the Sutlej, in Western Tibet.

Evidence for mechanical coupling and strong Indian lower crust beneath southern Tibet

How surface deformation within mountain ranges relates to tectonic processes at depth is not well understood. The upper crust of the Tibetan Plateau is generally thought to be poorly coupled to the underthrusting Indian crust because of an intervening low-viscosity channel. Here, however, we show that the contrast in tectonic regime between primarily strike-slip faulting in northern Tibet and dominantly normal faulting in southern Tibet requires mechanical coupling between the upper crust of southern Tibet and the underthrusting Indian crust. Such coupling is inconsistent with the presence of active 'channel flow' beneath southern Tibet, and suggests that the Indian crust retains its strength as it underthrusts the plateau. These results shed new light on the debates regarding the mechanical properties of the continental lithosphere, and the deformation of Tibet.

Indication for clockwise rotation in the Siang window south of the eastern Himalayan syntaxis and new geochronological constraints for the area

Palaeomagnetic, rock magnetic and geochronological investigations were carried out on the Abor volcanics of Arunachal Pradesh, NE India. A Late Palaeozoic formation age for part of the Abor volcanics cannot be excluded based on K–Ar whole rock dating. Low-temperature thermochronometers – zircon (U–Th)/He and fission track analyses – yield a maximum burial temperature of c. 150–170 °C during Late Miocene. ZFT thermochronology of the Yinkiong and Miri Fms. indicates a post-Paleocene and post-Jurassic deposition age, respectively. This infers that the volcanic rocks intercalating or intruding them are not part of the Late Palaeozoic sequence but represent one or more, latest Cretaceous to Tertiary event(s). Therefore the Abor volcanics are connected to at least two separate events of volcanism. From palaeomagnetic sites, two characteristic magnetic remanence components were separated: a low-coercivity-component demagnetized below 20 mT and a high-coercivity-component demagnetized between 15 and 100 mT. Fold tests support a secondary origin of both components. Thermochronological and rock magnetic analyses indicate a low-grade overprint event between India–Asia collision and Miocene, which probably represents the time of remanence acquisition. The high-coercivity-component shows a trend of clockwise declinations, which is likely related to vertical-axis rotations of the eastern Himalayas due to eastward extrusion of the Tibetan Plateau.

Tectonic geomorphology along the eastern margin of Tibet: insights into the pattern and processes of active deformation adjacent to the Sichuan Basin

We present a review and synthesis of the tectonic geomorphology along the eastern margin of the Tibetan Plateau adjacent to and north of the Sichuan Basin. Re-evaluation of spatial variations in the form of fluvial longitudinal profiles provides a refined image of the distribution of anomalously steep channels. Three new analyses demonstrate that these variations in channel steepness reflect variations in the locus and rate of differential rock uplift. First, measurements of channel width along trunk streams reveal systematic co-variations in channel hydraulic geometry and slope that suggests channels are dynamically adjusted to spatial variations in erosion rate. Second, recent determinations of the functional relationship between channel steepness indices and erosion rate allow a quantitative estimation of erosion rate from channel profile form. Third, comparison of rock uplift patterns to variations in the distribution of slip associated with the 2008 Wenchuan earthquake confirms that channel gradients reflect differential rock uplift. Our analysis suggests that reactivated fault systems adjacent to the Sichuan Basin are primarily responsible for accommodating differential rock uplift, but that rock uplift northward along the margin is not associated with active faults and is likely sustained by flow and thickening in the deep crust.

Landscape development of the Himalayan–Tibetan orogen: a review

The Himalayan–Tibetan orogen provides one of the best natural laboratories in which to examine the nature and dynamics of landscape development within continent–continent collision zones. Many new tectonic–climatic–geomorphological theories and models have emerged and/or have been greatly influenced as a consequence of the study of the region and the quest to understand its geomorphological development. These include models of the interactions between tectonics, climate and surfice processes, notably, the influence of climate on surface uplift by denudational unloading; the limiting of topography by glaciation (the glacial buzz-saw model); localized uplift at syntaxes by enhanced fluvial and glacial erosion that, in turn, weaken the lithosphere, enhancing surface uplift and exhumation (the tectonic aneurysm model); climate-driven out-of-sequence thrusting and crustal channel flow; glacial damming leading to differential erosion and uplift; paraglaciation; and the influence of extreme events such as earthquakes, landslides, and floods as major formative processes. The development of new technologies, including satellite remote sensing and global positioning systems, and analytical methods such as numerical dating is now allowing these theories and models to be tested and will inevitably lead to new paradigms.

Monsoon control over erosion patterns in the Western Himalaya: possible feed-back into the tectonic evolution

The Indus Delta is constructed of sediment eroded from the western Himalaya and since 20 ka has been subjected to strong variations in monsoon intensity. Provenance changes rapidly at 12–8 ka, although bulk and heavy mineral content remains relatively unchanged. Bulk sediment analyses shows more negative εNd and higher 87Sr/86Sr values, peaking around 8–9 ka. Apatite fission track ages and biotite Ar–Ar ages show younger grains ages at 8–9 ka compared to at the Last Glacial Maximum (LGM). At the same time δ13C climbs from –23 to –20‰, suggestive of a shift from terrestrial to more marine organic carbon as Early Holocene sea level rose. U–Pb zircon ages suggest enhanced erosion of the Lesser Himalaya and a relative reduction in erosion from the Transhimalaya and Karakoram since the LGM. The shift in erosion to the south correlates with those regions now affected by the heaviest summer monsoon rains. The focused erosion along the southern edge of Tibet required by current tectonic models for the Greater Himalaya would be impossible to achieve without a strong summer monsoon. Our work supports the idea that although long-term monsoon strengthening is caused by uplift of the Tibetan Plateau, monsoon-driven erosion controls Himalayan tectonic evolution.

Supplementary material:

A table of the population breakdown for zircons in sands and the predicted Nd isotope composition of sediments based on the zircons compared to the measured whole rock value is available at http://www.geolsoc.org.uk/SUP18412

The mountain-lowland debate: deforestation and sediment transport in the upper Ganga catchment

The Himalaya-Gangetic Plain region is the iconic example of the debate about the impact on lowlands of upland land-use change. Some of the scientific aspects of this debate are revisited by using new techniques to examine the role of deforestation in erosion and river sediment transport. The approach is whole-of-catchment, combining a history of deforestation with a history of sediment sources from well before deforestation. It is shown that deforestation had some effect on one very large erosional event in 1970, in the Alaknanda subcatchment of the Upper Ganga catchment, but that both deforestation and its effects on erosion and sediment transport are far from uniform in the Himalaya. Large magnitude erosional events occur for purely natural reasons. The impact on the Gangetic Plain of erosion caused by natural events and land cover change remains uncertain.

Near-isothermal conditions in the middle and lower crust induced by melt migration

The thermal structure of the crust strongly influences deformation, metamorphism and plutonism. Models for the geothermal gradient in stable crust predict a steady increase of temperature with depth. This thermal structure, however, is incompatible with observations from high-temperature metamorphic terranes exhumed in orogens. Global compilations of peak conditions in high-temperature metamorphic terranes define relatively narrow ranges of peak temperatures over a wide range in pressure, for both isothermal decompression and isobaric cooling paths. Here we develop simple one-dimensional thermal models that include the effects of melt migration. These models show that long-lived plutonism results in a quasi-steady-state geotherm with a rapid temperature increase in the upper crust and nearly isothermal conditions in the middle and lower crust. The models also predict that the upward advection of heat by melt generates granulite facies metamorphism, and widespread andalusite-sillimanite metamorphism in the upper crust. Once the quasi-steady-state thermal profile is reached, the middle and lower crust are greatly weakened due to high temperatures and anatectic conditions, thus setting the stage for gravitational collapse, exhumation and isothermal decompression after the onset of plutonism. Near-isothermal conditions in the middle and lower crust result from the thermal buffering effect of dehydration melting reactions that, in part, control the shape of the geotherm.

Hydrogen isotopes in Eocene river gravels and paleoelevation of the Sierra Nevada

We determine paleoelevation of the Sierra Nevada, California, by tracking the effect of topography on precipitation, as recorded in hydrogen isotopes of kaolinite exposed in gold-bearing river deposits from the Eocene Yuba River. The data, compared with the modern isotopic composition of precipitation, show that about 40 to 50 million years ago the Sierra Nevada stood tall (>/=2200 meters), a result in conflict with proposed young surface uplift by tectonic and climatic forcing but consistent with the Sierra Nevada representing the edge of a pre-Eocene continental plateau.

Provenance of the Greater Himalayan Sequence and associated rocks: predictions of channel flow models

Numerical models for channel flow in the Himalayan—Tibetan system are compatible with many tectonic and metamorphic features of the orogen. Here we compare the provenance of crustal material in two channel flow models (HT1 and HT111) with observations from the Himalaya and southern Tibet. Thirty million years after the onset of channel flow, the entire model crust south of the India—Asia suture still consists only of ‘Indian’ material. The model Greater Himalayan Sequence (‘GHS’) is derived from Indian middle crust originating ≤1000 km south of the initial position of the suture, whereas the Lesser Himalayan Sequence (‘LHS’) is derived mainly from crust originating ≥1400 km south of the suture. Material tracking indicates little or no mixing of diverse crustal elements in the exhumed region of the model ‘GHS’, which is derived from originally contiguous materials that are transported together in the top of the channel flow zone. These results are compatible with provenance data indicating a clear distinction between GHS and LHS protoliths, with the GHS originating from a more distal position (relative to cratonic India) than the LHS. In model HT111, domes formed between the suture and the orogenic front are cored by ‘Indian’ middle crust similar to the ‘GHS’, consistent with data from the north Himalayan gneiss domes. Material tracking shows that plutons generated south of the suture should have ‘Indian’ crustal signatures, also compatible with observations. Model ‘GHS’ pressure—temperature—time (P-T-t) paths pass through the dehydration melting field between 30 and 15 Ma, consistent with observed leucogranite ages. Finally, exposure of midcrustal ‘GHS’ and ‘LHS’ material at the model erosion front is consistent with the observed appearance of sedimentary detritus in the Lesser Himalaya. We conclude that channel flow model results are compatible with provenance data from the Himalaya and southern Tibet.

Growth of continental plateaus by channel injection: models designed to address constraints and thermomechanical consistency

Weak, possibly partially molten, middle crust may exist and deform by channel flow beneath continental plateaus, thereby significantly influencing their dynamics. The role of channel flows in the transition zone between the plateau and the foreland is, however, unclear. We develop successively more complete approximate models for the channel injection (CI) mode in which differential pressure pumps channel material from beneath the plateau into the transitional crust, which thickens it and widens the plateau. The motivation is to improve our understanding of the controls on the growth of continental plateaus and the interactions in the transition zone, and to gain more insight into the results of more complex numerical models. In model CI-1, a channel with constant viscosity and thickness exists in the transitional crust and the pumped material accretes/freezes above and below the channel. Although results compare favourably with the geometry of some natural examples, this model is incomplete because the connection between the transition zone and the plateau is not considered. Model CI-2 includes a decrease in channel viscosity when the channel depth exceeds a critical value, D*, a proxy for onset of melt weakening or low viscosities at high temperatures. The model completes the connection to the plateau, but relies on the arbitrary choice of D*. Model CI-3 is more physically based, and considers the channel viscosity and thickness to depend on temperature, calculated by an associated thermal model that includes radioactive self-heating, and advection of heat by channel material. This model demonstrates self-consistent plateau widening if the channel viscosity decreases at the critical temperature, T*. Acceptable comparisons with the topography of Tibet are achieved with transition zone viscosities that decrease from 1019−1022 Pa s to subplateau values of 1018–1019 Pa s, with T* of 700−750°C. Additional analyses and tests are used to determine the range of parameter values for which CI models are both internally consistent and compatible with observations. Additional modes of deformation in the transition zone, viscous thickening (VT) and plastic translation (PT), may also be important.