Partially Molten Middle Crust Beneath Southern Tibet: Synthesis of Project INDEPTH Results

Crustal melting and continent uplift by mafic underplating at convergent boundaries

The thick crust of the southern Tibetan and central Andean plateaus includes high-conductivity, low-velocity zones ascribed to partial melt. The melt origin and effect on plateau uplift remain speculative, in particular if plateau uplift happens before continental collision. The East Anatolian Plateau (EAP) has experienced similar, more recent uplift but its structure is largely unknown. Here we present an 80 km deep geophysical model across EAP, constrained by seismic receiver functions integrated with interpretation of gravity data and seismic tomographic, magnetotelluric, geothermal, and geochemical models. The results indicate a 20 km thick lower crustal layer and a 10 km thick mid-crustal layer, which both contain pockets of partial melt. We explain plateau uplift by isostatic equilibration following magmatism associated with roll-back and break-off of the Neo-Tethys slab. Our results suggest that crustal thickening by felsic melt and mafic underplate are important for plateau uplift in the EAP, Andes and Tibet.

The crustal juvenile and reworking in eastern Tibet due to tectonism and magmatism

Crustal thickening in Tibet results from the collision of India and Eurasia. The fine structure of the crust is vital information to detect the tectonic genesis of crustal growth and reworking. Here, we constructed the crustal structure of eastern Tibet and its neighboring areas using depth-domain receiver function imaging. The imaging results reveal widespread low-velocity layers in the mid-crust, alternating high- and low-velocity structures in the lower crust, and regional variations in the velocity at the bottom of the lower crust and the thickness of the entire crust. We propose that the lower crust is juvenile crust, formed by multiphase magmatism and tectonic activity during continental collision. The regional variations in crustal structure can be attributed to differences in the tectonic position of subduction-collision system and the distinct crustal-derived materials that contribute to the formation of the juvenile crust. The low-velocity layers in the mid-crust are more likely to be the product of arc magmatism during the assembly of micro-continents. Based on the constraints of the unique crustal structure, we suggest that tectonic compression and magmatism caused the juvenile and reworking of the Tibetan crust.

Upward and outward growth of north-central Tibet: Mechanisms that build high-elevation, low-relief plateaus

Large orogenic plateaus, such as the Tibetan Plateau, are characterized by high-elevation, low-relief topography, in contrast to the rugged terrains of narrower mountain belts. A key question is how low-elevation hinterland basins, characteristic of broad regions of shortening, were raised while regional relief was flattened. This study uses the Hoh Xil Basin in north-central Tibet as an analogue for late-stage orogenic plateau formation. The precipitation temperatures of lacustrine carbonates deposited between ~19 and ~12 million years ago record an early to middle Miocene phase of surface uplift of 1.0 ± 0.7 km. The results of this study demonstrate the contribution of sub-surface geodynamic processes in driving regional surface uplift and redistribution of crustal material to flatten plateau surfaces during the late stage of orogenic plateau formation.

Characteristics, sources, and risk assessment of thallium and associated with metal(loid)s in the Yarlung Tsangpo River Basin, southern Tibetan Plateau

The Tibetan Plateau (TP) is known as the water tower of Asia, and the water quality has long been a focus of public concern, especially in the Yarlung Tsangpo River Basin (YTRB), a unique area that is climate-sensitive, geologically complex, eco-fragile, and densely populated. Thallium (Tl) is a typical metal that is more toxic than Pb, Cd, and As and often occurs in sulfide minerals. Although large-scale polymetallic sulfide mineralization developed in the YTRB, the geochemical dispersion and potential risk of Tl in aquatic environments of the YTRB remain poorly understood. In this study, the concentration, distribution, source, and health risk of Tl and associated metal(loid)s in the hot springs and surface water in the YTRB were systematically analyzed. The results showed that the trace elements (Cd, Cr, Zn, Cu, Al, Sr, Ni, Co, Mn, Pb) in water environments are within the recommended limits, except for Tl and As. Principal component analysis (PCA) and correlation analysis (CA) showed that the elements of Tl and As were positively related to each other in either both hot spring water and surface water, indicating their common origin. Spatial variations suggested that high levels of Tl and As observed in the north YTRB, which may be relevant to the reduction-dissolution of Tl (As)-bearing minerals and the magmatic hydrothermal system formed in the shallow part of the northern YTRB. Furthermore, source apportionment identified natural sources of Cu, Ni, Cr, Co, Mn, Zn, and Cd and anthropogenic inputs of Al and Pb. Exposure assessment studies have found that ingestion is the primary route of As and Tl exposure to local population, and balneological and bathing purposes do not constitute a human health concern. This study offers valuable insights into the risk of naturally occurring Tl enrichment being hidden in As-rich hydrosphere in the YTRB and other regions with similar geoenvironmental contexts.

Crustal melting in orogenic belts revealed by eclogite thermal properties

Partial melting in the continental crust may play a critical role on the behavior of continents during collision. However, the occurrence of partial melt in orogenic continental crust is not well understood. Since the temperature of the orogen is controlled by the thermal properties of constituent rocks, we measured the thermal conductivity and diffusivity of eclogite, the most important ultrahigh pressure metamorphic rocks, as a function of pressure, temperature, composition, and water content, and simulated the thermal structure of the Sulu and Himalaya-Tibet orogens in eastern and southwestern China, respectively. Our results show that the temperature at ~30-km depth beneath the orogens reaches the solidus of wet granite and phengite (~940 K), therefore, the partial melting in the orogenic continental crust is well explained. The melt may facilitate the exhumation of subducted crust, produce the low seismic-velocity zone, and cause the high-conductivity anomaly in the shallow depth of orogenic belts.

By measuring the thermal properties of eclogite at high pressures, the authors found that temperature of orogenic continental crust is sufficient to melt granite and phengite, leading to low-velocity and high-conductivity anomalies in orogenic belts.

The Hindu Kush slab break-off as revealed by deep structure and crustal deformation

Break-off of part of the down-going plate during continental collision occurs due to tensile stresses built-up between the deep and shallow slab, for which buoyancy is increased because of continental-crust subduction. Break-off governs the subsequent orogenic evolution but real-time observations are rare as it happens over geologically short times. Here we present a finite-frequency tomography, based on jointly inverted local and remote earthquakes, for the Hindu Kush in Afghanistan, where slab break-off is ongoing. We interpret our results as crustal subduction on top of a northwards-subducting Indian lithospheric slab, whose penetration depth increases along-strike while thinning and steepening. This implies that break-off is propagating laterally and that the highest lithospheric stretching rates occur during the final pinching-off. In the Hindu Kush crust, earthquakes and geodetic data show a transition from focused to distributed deformation, which we relate to a variable degree of crust-mantle coupling presumably associated with break-off at depth.

Here, the authors document active slab break-off and the crustal response during continental collision under the Hindu Kush, a rarely observed process since it happens over geologically short time spans.

Revised chronology of central Tibet uplift (Lunpola Basin)

Knowledge of the topographic evolution of the Tibetan Plateau is essential for understanding its construction and its influences on climate, environment, and biodiversity. Previous elevations estimated from stable isotope records from the Lunpola Basin in central Tibet, which indicate a high plateau since at least 35 Ma, are challenged by recent discoveries of low-elevation tropical fossils apparently deposited at 25.5 Ma. Here, we use magnetostratigraphic and radiochronologic dating to revise the chronology of elevation estimates from the Lunpola Basin. The updated ages reconcile previous results and indicate that the elevations of central Tibet were generally low (<2.3 km) at 39.5 Ma and high (3.5 to 4.5 km) at ~26 Ma. This supports the existence in the Eocene of low-elevation longitudinally oriented narrow regions until their uplift in the early Miocene, with potential implications for the growth mechanisms of the Tibetan Plateau, Asian atmospheric circulation, surface processes, and biotic evolution.

Seismic observations, numerical modeling, and geomorphic analysis of a glacier lake outburst flood in the Himalayas

Glacial lake outburst floods (GLOFs) are a substantial hazard for downstream communities in vulnerable regions, yet unpredictable triggers and remote source locations make GLOF dynamics difficult to measure and quantify. Here, we revisit a destructive GLOF that occurred in Bhutan in 1994 and apply cross-correlation-based seismic analyses to track the evolution of the GLOF remotely (~100 kilometers from the source region). We use the seismic observations along with eyewitness reports and a downstream gauge station to constrain a numerical flood model and then assess geomorphic change and current state of the unstable lakes via satellite imagery. Coherent seismic energy is evident from 1 to 5 hertz beginning approximately 5 hours before the flood impacted Punakha village, which originated at the source lake and advanced down the valley during the GLOF duration. Our analysis highlights potential benefits of using real-time seismic monitoring to improve early warning systems.

Active strike-slip faults and an outer frontal thrust in the Himalayan foreland basin

The Himalayan foreland basin formed by flexure of the Indian Plate below the advancing orogen. Motion on major thrusts within the orogen has resulted in damaging historical seismicity, whereas south of the Main Frontal Thrust (MFT), the foreland basin is typically portrayed as undeformed. Using two-dimensional seismic reflection data from eastern Nepal, we present evidence of recent deformation propagating >37 km south of the MFT. A system of tear faults at a high angle to the orogen is spatially localized above the Munger-Saharsa basement ridge. A blind thrust fault is interpreted in the subsurface, above the sub-Cenozoic unconformity, bounded by two tear faults. Deformation zones beneath the Bhadrapur topographic high record an incipient tectonic wedge or triangle zone. The faults record the subsurface propagation of the Main Himalayan Thrust (MHT) into the foreland basin as an outer frontal thrust, and provide a modern snapshot of the development of tectonic wedges and lateral discontinuities preserved in higher thrust sheets of the Himalaya, and in ancient orogens elsewhere. We estimate a cumulative slip of ∼100 m, accumulated in <0.5 Ma, over a minimum slipped area of ∼780 km² These observations demonstrate that Himalayan ruptures may pass under the present-day trace of the MFT as blind faults inaccessible to trenching, and that paleoseismic studies may underestimate Holocene convergence.

Direct structural evidence of Indian continental subduction beneath Myanmar

Indian continental subduction can explain Cenozoic crustal deformation, magmatic activity and uplift of the Tibetan Plateau following the India-Asia collision. In the western Himalayan syntaxis and central Himalaya, subduction or underthrusting of the Indian Plate beneath the Eurasian Plate is well known from seismological studies. However, because information on the deep structure of the eastern Himalayan syntaxis is lacking, the nature of the Indian subduction slab beneath Myanmar and the related tectonic regime remain unclear. Here, we use receiver function common conversion point imaging from a densely spaced seismic array to detect direct structural evidence of present-day Indian continental subduction beneath Asia. The entire subducting Indian crust has an average crustal thickness of ~30 km, dips at an angle of ~19°, and extends to a depth of 100 km under central Myanmar. These results reveal a unique continental subduction regime as a result of Indian-Eurasian continental collision and lateral extrusion.

Indian continental subduction can explain crustal deformation, magmatic activity and uplift of the Tibetan Plateau following collision, however, the nature of the Indian subducting slab beneath Myanmar and the related tectonic regime remain unclear. Here, the authors present direct structural evidence of present-day Indian continental subduction beneath Asia.

Lithospheric delamination and upwelling asthenosphere in the Longmenshan area: insight from teleseismic P-wave tomography

We apply teleseismic P-wave tomography to reconstruct the velocity structure of the Longmenshan area. Our results show possible large-scale delamination beneath the Songpan-Ganzi and Qiangtang terranes, which induced upwelling asthenosphere. Upwelling asthenosphere might have led to lower crust heating, facilitating eastward extrusion of the Songpan Ganzi terrane resulting in localized deformation and uplift along the Longmenshan orogenic belt. We suggest that the eastward extrusion of the Songpan-Ganzi terrane against the rigid lithospheric root of the Sichuan Basin results in stress accumulation and release, leading to large earthquakes in the Longmenshan area.

Bimodal seismicity in the Himalaya controlled by fault friction and geometry

There is increasing evidence that the Himalayan seismicity can be bimodal: blind earthquakes (up to Mw ~ 7.8) tend to cluster in the downdip part of the seismogenic zone, whereas infrequent great earthquakes (Mw 8+) propagate up to the Himalayan frontal thrust. To explore the causes of this bimodal seismicity, we developed a two-dimensional, seismic cycle model of the Nepal Himalaya. Our visco-elasto-plastic simulations reproduce important features of the earthquake cycle, including interseismic strain and a bimodal seismicity pattern. Bimodal seismicity emerges as a result of relatively higher friction and a non-planar geometry of the Main Himalayan Thrust fault. This introduces a region of large strength excess that can only be activated once enough stress is transferred upwards by blind earthquakes. This supports the view that most segments of the Himalaya might produce complete ruptures significantly larger than the 2015 Mw 7.8 Gorkha earthquake, which should be accounted for in future seismic hazard assessments.

There is increasing evidence that the seismicity of large Himalayan earthquakes can be bimodal, but the underlying mechanisms are unclear. Here, the authors present a model and show that the bimodal seismicity results from a relatively higher friction and a non-planar geometry of the Himalayan megathrust.

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.

Outlining tectonic inheritance and construction of the Min Shan region, eastern Tibet, using crustal geometry

The ongoing collision between India and Eurasia has created the Tibetan Plateau, which features high elevations and large crustal thicknesses. The easternmost portion of the plateau has long been a key region for studying the uplift mechanism of the Tibetan Plateau, especially after the 2008 Ms. 7.9 Wenchuan earthquake. However, previous studies have assumed that easternmost Tibet is tectonically homogeneous, and the tectonic significance of the Min Shan has been overshadowed by that of its more conspicuous neighbour, the Longmen Shan region. Here, we describe the crustal geometry of the Min Shan region using two newly obtained deep seismic reflection profiles. In this study, we identify an upper-lower crust mechanical decoupling within the Min Shan region; the Min Shan region is tectonically delineated by an inherited boundary fault zone, the Huya fault zone, which was responsible for triggering the 2017 Jiuzhaigou M 7.0 earthquake. Together with the gravity dataset and previous studies in this area, the outlined crustal geometry indicated that crustal-scale shortening at the eastern plateau margin is a primary mechanism driving uplift, although extensive uplift might have occurred due to the decoupled shortening between the upper and lower crust.

Nonuniform subduction of the Indian crust beneath the Himalayas

Himalayan tectonic activity is triggered by downward penetration of the Indian plate beneath the Asian plate. The subsurface geometry of this interaction has not been fully investigated. This study presents novel constraints on this geometry provided by two newly obtained, deep seismic reflection profiles. The profiles cover 100- and 60-km transects across the Yarlung-Zangbo suture of the Himalaya-Tibet orogen at c. 88°E. Both profiles show a crustal-scale outline of the subducting Indian crust. This outline clearly shows Indian understhrusting southern Tibet, but only to a limited degree. When combined with a third seismic reflection profile of the western Himalayas, the new profiles reveal progressive, eastward steepening and shortening in the horizontal advance of the subducting Indian crust.

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.

Telomere Dysfunction Triggers Palindrome Formation Independently of Double-Strand Break Repair Mechanisms

Inverted chromosome duplications or palindromes are linked with genetic disorders and malignant transformation. They are considered by-products of DNA double-strand break (DSB) repair: the homologous recombination (HR) and the nonhomologous end joining (NHEJ). Palindromes near chromosome ends are often triggered by telomere losses. An important question is to what extent their formation depends upon DSB repair mechanisms. Here we addressed this question using yeast genetics and comparative genomic hybridization. We induced palindrome formation by passaging cells lacking any form of telomere maintenance (telomerase and telomere recombination). Surprisingly, we found that DNA ligase 4, essential for NHEJ, did not make a significant contribution to palindrome formation induced by telomere losses. Moreover RAD51, important for certain HR-derived mechanisms, had little effect. Furthermore RAD52, which is essential for HR in yeast, appeared to decrease the number of palindromes in cells proliferating without telomeres. This study also uncovered an important role for Rev3 and Rev7 (but not for Pol32) subunits of polymerase ζ in the survival of cells undergoing telomere losses and forming palindromes. We propose a model called short-inverted repeat-induced synthesis in which DNA synthesis, rather than DSB repair, drives the inverted duplication triggered by telomere dysfunction.

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.

Moho topography, ranges and folds of Tibet by analysis of global gravity models and GOCE data

The determination of the crustal structure is essential in geophysics, as it gives insight into the geohistory, tectonic environment, geohazard mitigation, etc. Here we present the latest advance on three-dimensional modeling representing the Tibetan Mohorovičić discontinuity (topography and ranges) and its deformation (fold), revealed by analyzing gravity data from GOCE mission. Our study shows noticeable advances in estimated Tibetan Moho model which is superior to the results using the earlier gravity models prior to GOCE. The higher quality gravity field of GOCE is reflected in the Moho solution: we find that the Moho is deeper than 65 km, which is twice the normal continental crust beneath most of the Qinghai-Tibetan plateau, while the deepest Moho, up to 82 km, is located in western Tibet. The amplitude of the Moho fold is estimated to be ranging from −9 km to 9 km with a standard deviation of ~2 km. The improved GOCE gravity derived Moho signals reveal a clear directionality of the Moho ranges and Moho fold structure, orthogonal to deformation rates observed by GPS. This geophysical feature, clearly more evident than the ones estimated using earlier gravity models, reveals that it is the result of the large compressional tectonic process.

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.

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

Trace element hydrochemistry indicating water contamination in and around the Yangbajing geothermal field, Tibet, China

Thirty-eight water samples were collected at Yangbajing to investigate the water contamination resulting from natural geothermal water discharge and anthropogenic geothermal wastewater drainage. The results indicate that snow or snow melting waters, Yangbajing River waters and cold groundwaters are free from geothermal water-related contamination, whereas Zangbo river waters are contaminated by geothermal wastewaters. Moreover, there may exist geothermal springs under the riverbed of a tributary stream of Zangbo River as shown by its Cd, Li, Mo and Pb concentrations. The efforts made in this study show trace element hydrochemistry can well indicate water quality degradation related to geothermal water exploitation.

Underplating in the Himalaya-Tibet collision zone revealed by the Hi-CLIMB experiment

We studied the formation of the Himalayan mountain range and the Tibetan Plateau by investigating their lithospheric structure. Using an 800-kilometer-long, densely spaced seismic array, we have constructed an image of the crust and upper mantle beneath the Himalayas and the southern Tibetan Plateau. The image reveals in a continuous fashion the Main Himalayan thrust fault as it extends from a shallow depth under Nepal to the mid-crust under southern Tibet. Indian crust can be traced to 31 degrees N. The crust/mantle interface beneath Tibet is anisotropic, indicating shearing during its formation. The dipping mantle fabric suggests that the Indian mantle is subducting in a diffuse fashion along several evolving subparallel structures.

The Appalachian Inner Piedmont: an exhumed strike-parallel, tectonically forced orogenic channel

The Appalachian Inner Piedmont (IP) extends along orogenic strike some 700 km from North Carolina to Alabama. Its physical attributes contrast with those of other Appalachian tectonic elements: gentle dip of dominant foliation; imbricate stack of fold nappes; dominant sillimanite-grade metamorphism and near ubiquitous migmatization; heterogeneous, non-plane deformation; and earlier S-foliations transposed to C-foliations southeast of the mid-Palaeozoic Brevard fault zone forming a 10–20 km wide amphibolite-facies shear zone along the western flank of the IP. The IP contains west- and SW-directed thrust sheets and mineral stretching lineation, sheath folds on all scales, and other indicators that define a curved crustal flow pattern throughout the belt. Field and modern geochronologic data confirm that the IP is not exotic. It contains a Laurentian component (eastern Tugaloo terrane) and an internal terrane (Cat Square) that contains both Laurentian and Gondwanan detrital zircons, separated by the Brindle Creek fault. Cat Square terrane rocks likely accumulated in a Devonian remnant ocean that closed beginning c. 400 Ma.

The complex but consistently asymmetric, NW- to west- to SW-directed flow pattern throughout the IP reflects confinement beneath a > 15 km thick overburden produced during subduction of Cat Square and Laurentian components beneath the approaching Carolina superterrane along the Central Piedmont suture. Oblique NE-to-SW transpressive subduction to > 15 km depth initiated partial melting, forcing escape from the collision zone in an along-strike orogenic channel. The IP detached from rocks to the west of the mid-Palaeozoic Brevard fault zone as the collision zone tightened and the IP mass flowed c. 200 km southwestward in the channel. The top of the channel is preserved at the NE end of the IP, and the base (Brevard fault zone) is preserved to the west and SW. As an exhumed orogenic channel, the curved IP flow paths may provide insight for middle to lower crustal deformation and flow in modern orogens.

Crustal flow modes in large hot orogens

Crustal-scale channel flow numerical models support recent interpretations of Himalayan—Tibetan tectonics proposing that gravitationally driven channel flows of low-viscosity, melt-weakened, middle crust can explain both outward growth of the Tibetan Plateau and ductile extrusion of the Greater Himalayan Sequence. We broaden the numerical model investigation to explore three flow modes: homogeneous channel flow (involving laterally homogeneous crust); heterogeneous channel flow (involving laterally heterogeneous lower crust that is expelled and incorporated into the mid-crustal channel flow); and the hot fold nappes style of flow (in which mid-/lower crust is forcibly expelled outward over a lower crustal indentor to create fold nappes that are inserted into the mid-crust). The three flow modes are members of a continuum in which the homogeneous mode is driven by gravitational forces but requires very weak channel material. The hot fold nappe mode is driven tectonically by, for example, collision with a strong crustal indentor and can occur in crust that is subcritical for homogeneous flows. The heterogeneous mode combines tectonic and gravitationally driven flows. Preliminary results also demonstrate the existence and behaviour of mid-crustal channels during advancing and retreating dynamical mantle lithosphere subduction. An orogen temperature—magnitude (T-M) diagram is proposed and the positions of orogens in T-M space that may exhibit the flow modes are described, together with the characteristic positions of a range of other orogen types.

Channel flow, ductile extrusion and exhumation in continental collision zones: an introduction

The channel flow model aims to explain features common to metamorphic hinterlands of some collisional orogens, notably along the Himalaya-Tibet system. Channel flow describes a protracted flow of a weak, viscous crustal layer between relatively rigid yet deformable bounding crustal slabs. Once a critical low viscosity is attained (due to partial melting), the weak layer flows laterally due to a horizontal gradient in lithostatic pressure. In the Himalaya-Tibet system, this lithostatic pressure gradient is created by the high crustal thicknesses beneath the Tibetan Plateau and ‘normal’ crustal thickness in the foreland. Focused denudation can result in exhumation of the channel material within a narrow, nearly symmetric zone. If channel flow is operating at the same time as focused denudation, this can result in extrusion of the mid-crust between an upper normal-sense boundary and a lower thrust-sense boundary. The bounding shear zones of the extruding channel may have opposite shear sense; the sole shear zone is always a thrust, while the roof shear zone may display normal or thrust sense, depending on the relative velocity between the upper crust and the underlying extruding material. This introductory chapter addresses the historical, theoretical, geological and modelling aspects of channel flow, emphasizing its applicability to the Himalaya-Tibet orogen. Critical tests for channel flow in the Himalaya, and possible applications to other orogenic belts, are also presented.

Crustal rheology of the Himalaya and Southern Tibet inferred from magnetotelluric data

The Cenozoic collision between the Indian and Asian continents formed the Tibetan plateau, beginning about 70 million years ago. Since this time, at least 1,400 km of convergence has been accommodated by a combination of underthrusting of Indian and Asian lithosphere, crustal shortening, horizontal extrusion and lithospheric delamination. Rocks exposed in the Himalaya show evidence of crustal melting and are thought to have been exhumed by rapid erosion and climatically forced crustal flow. Magnetotelluric data can be used to image subsurface electrical resistivity, a parameter sensitive to the presence of interconnected fluids in the host rock matrix, even at low volume fractions. Here we present magnetotelluric data from the Tibetan-Himalayan orogen from 77 degrees E to 92 degrees E, which show that low resistivity, interpreted as a partially molten layer, is present along at least 1,000 km of the southern margin of the Tibetan plateau. The inferred low viscosity of this layer is consistent with the development of climatically forced crustal flow in Southern Tibet.