Detection of widespread fluids in the Tibetan crust by magnetotelluric studies

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.

Correlated crustal and mantle melting documents proto-Tibetan Plateau growth

The mechanism that causes the rapid uplift and active magmatism of the Hoh-Xil Basin in the northern Tibetan Plateau and hence the outward growth of the proto-plateau is highly debated, more specifically, over the relationship between deep dynamics and surface uplift. Until recently the Hoh-Xil Basin remained uncovered by seismic networks due to inaccessibility. Here, based on linear seismic arrays across the Hoh-Xil Basin, we present a three-dimensional S-wave velocity (VS) model of the crust and uppermost mantle structure beneath the Tibetan Plateau from ambient noise tomography. This model exhibits a widespread partially molten crust in the northern Tibetan Plateau but only isolated pockets in the south manifested as low-VS anomalies in the middle crust. The spatial correlation of the widespread low-VS anomalies with strong uppermost mantle low-VS anomalies and young exposed magmatic rocks in the Hoh-Xil Basin suggests that the plateau grew through lithospheric mantle removal and its driven magmatism.

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.

Southern Tibetan rifting since late Miocene enabled by basal shear of the underthrusting Indian lithosphere

Syncontractional extension is prominent in present-day Tibet, but its origin remains vigorously debated. Several deep-seated geodynamic processes (e.g., Indian underthrusting, horizontal flow, and mantle upwelling) have been linked to Tibetan rifting. Indian underthrusting is a good candidate because it can well explain why surface rifts are more prominent south of the Bangong-Nujiang suture; however, how Indian underthrusting causes extension is not well understood and lacks observational constraints. Seismic anisotropy, measured by exploiting the birefringence effect of shear waves, can be indicative of the deformation styles within the crust. Here, we unveil the dominant convergence-parallel alignment of anisotropic fabrics in the deep crust of the southern Tibetan rifts using seismic recordings collected from our recently deployed and existing seismic stations. This finding suggests that the strong north-directed shearing exerted by the underthrusting Indian plate is key to enabling present-day extension in southern Tibet.

A Review of Subsurface Electrical Conductivity Anomalies in Magnetotelluric Imaging

After 70 years of development, magnetotelluric (MT), a remote sensing technique for subsurface electrical resistivity imaging, has been widely applied in resource exploration and the deep tectonic evolution of the Earth. The electrical resistivity anomalies and their quantitative interpretation are closely related to or even controlled by the interconnected high-conductivity phases, which are frequently associated with tectonic activity. Based on representative electrical resistivity studies mainly of the deep crust and mantle, we reviewed principal electrical conduction mechanisms, generally used conductivity mixing models, and potential causes of high-conductivity including the saline fluid, partial melting, graphite, sulfide, and hydrogen in nominally anhydrous minerals, and the general methods to infer the water content of the upper mantle through electrical anomaly revealed by MT.

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.

Barium isotope evidence for the role of magmatic fluids in the origin of Himalayan leucogranites

As an important post-collisional magmatic product in the orogenic belt, the Himalayan leucogranites are the critical host rocks for a number of rare-metal mineralization (such as Li, Be, Cs, Rb, Nb, Ta, and Sn). However, there is still a lack of good understanding on the formation and evolution of the leucogranites. Particularly, the role of the magmatic fluids in transporting and enriching the rare elements is not clear. Here we measure Ba isotope compositions for leucogranites from the Kampa Dome of the Himalayan belt to understand the fluid activity and behavior of fluid-mobile elements during leucogranite formation. Our results show that the δ^138/134Ba of leucogranites range from -1.32‰ to +0.12‰, much lower than the literature values for S-type granites and various sedimentary materials, suggesting that the Ba isotope compositions of the leucogranites does not reflect the sedimentary source signatures. Instead, their low δ^138/134Ba is accompanied by non-charge-and-radius-controlled (CHARAC) twin-element (such as Nb/Ta) behaviors, clearly showing the involvement of magmatic fluids during magma evolution. Experimental studies suggest that the low δ^138/134Ba of the magmatic fluids most likely results from exsolution from a large deep magma reservoir. Such fluids not only modified Ba isotope compositions of the leucogranites, but also transported many fluid-mobile metal elements which may help form the rare metal ore deposits. Therefore, Ba isotope data provide new insights into formation and evolution of magmatic fluids and exploration of the rare-metal mineralization.

No mafic layer in 80 km thick Tibetan crust

All models of the magmatic and plate tectonic processes that create continental crust predict the presence of a mafic lower crust. Earlier proposed crustal doubling in Tibet and the Himalayas by underthrusting of the Indian plate requires the presence of a mafic layer with high seismic P-wave velocity (V_p > 7.0 km/s) above the Moho. Our new seismic data demonstrates that some of the thickest crust on Earth in the middle Lhasa Terrane has exceptionally low velocity (V_p < 6.7 km/s) throughout the whole 80 km thick crust. Observed deep crustal earthquakes throughout the crustal column and thick lithosphere from seismic tomography imply low temperature crust. Therefore, the whole crust must consist of felsic rocks as any mafic layer would have high velocity unless the temperature of the crust were high. Our results form basis for alternative models for the formation of extremely thick juvenile crust with predominantly felsic composition in continental collision zones.

Extensive Sills in the Continental Basement from Deep Seismic Reflection Profiling

Crustal seismic reflection profiling has revealed the presence of extensive, coherent reflections with anomalously high amplitudes in the crystalline crust at a number of locations around the world. In areas of active tectonic activity, these seismic “bright spots” have often been interpreted as fluid magma at depth. The focus in this report is high-amplitude reflections that have been identified or inferred to mark interfaces between solid mafic intrusions and felsic to intermediate country rock. These “frozen sills” most commonly appear as thin, subhorizontal sheets at middle to upper crustal depths, several of which can be traced for tens to hundreds of kilometers. Their frequency among seismic profiles suggest that they may be more common than widely realized. These intrusions constrain crustal rheology at the time of their emplacement, represent a significant mode of transfer of mantle material and heat into the crust, and some may constitute fingerprints of distant mantle plumes. These sills may have played important roles in overlying basin evolution and ore deposition.

Multi-Parametric Climatological Analysis Reveals the Involvement of Fluids in the Preparation Phase of the 2008 Ms 8.0 Wenchuan and 2013 Ms 7.0 Lushan Earthquakes

A multi-parametric approach was applied to climatological data before the Ms 8.0 2008 Wenchuan and Ms 7.0 2013 Lushan earthquakes (EQs) in order to detect anomalous changes associated to the preparing phase of those large seismic events. A climatological analysis for seismic Precursor Identification (CAPRI) algorithm was used for the detection of anomalies in the time series of four parameters (aerosol optical depth, AOD; skin temperature, SKT; surface latent heat flux, SLHF and total column water vapour, TCWV). Our results show a chain of processes occurred within two months before the EQs: AOD anomalous response is the earliest, followed by SKT, TCWV and SLHF in the EQs. A close spatial relation between the seismogenic Longmenshan fault (LMSF) zone and the extent of the detected anomalies indicates that some changes occurred within the faults before the EQs. The similarity of time sequence of the anomalies between the four parameters may be related to the same process: we interpret the observed anomalies as the consequence of the upraising of gases from a fluid-rich middle/upper crust along pre-existing seismogenic faults, and of their release into the atmosphere. Our multi-parametric analytical approach is able to capture phenomena related to the preparation phase of strong EQs.

Normal faulting and viscous buckling in the Tibetan Plateau induced by a weak lower crust

Flow of weak lower crust has been invoked to reconcile observed topographic gradients, uniform elevations, slow seismic velocity, and high conductivity measured in the Tibetan Plateau, with viscosity estimates of 10¹⁶-10²¹ Pa·s. Here we investigate the dynamic response resulting from a range of lower crust viscosities in a 3-D lithospheric-scale geodynamic simulation of the India-Eurasia collision zone to determine bounds of physically viable lower crustal strengths. We show that thickening of the plateau is accommodated through viscous buckling of the upper crust in response to lower crustal flow for a lower crustal viscosity on the order of 10²⁰ Pa·s. This generates two east-west trending bands of surface subsidence and dilatation consistent with observed normal faulting and estimates of vertical velocity. These results suggest viscous buckling of the upper crust, induced by lower crustal flow from gravitational pressure gradients due to high topography, is responsible for the observed extension in Tibet.

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.

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.

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.

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.

Seismic images of crust and upper mantle beneath Tibet: evidence for Eurasian plate subduction

Seismic data from central Tibet have been combined to image the subsurface structure and understand the evolution of the collision of India and Eurasia. The 410- and 660-kilometer mantle discontinuities are sharply defined, implying a lack of a subducting slab beneath the plateau. The discontinuities appear slightly deeper beneath northern Tibet, implying that the average temperature of the mantle above the transition zone is about 300 degrees C hotter in the north than in the south. There is a prominent south-dipping converter in the uppermost mantle beneath northern Tibet that might represent the top of the Eurasian mantle lithosphere underthrusting the northern margin of the plateau.

Normal faulting in central Tibet since at least 13.5 Myr ago

Tectonic models for the evolution of the Tibetan plateau interpret observed east-west thinning of the upper crust to be the result of either increased potential energy of elevated crust or geodynamic processes that may be unrelated to plateau formation. A key piece of information needed to evaluate these models is the timing of deformation within the plateau. The onset of normal faulting has been estimated to have commenced in southern Tibet between about 14 Myr ago and about 8 Myr ago and, in central Tibet, about 4 Myr ago. Here, however, we report a minimum age of approximately 13.5 Myr for the onset of graben formation in central Tibet, based on mineralization ages determined with Rb-Sr and 40Ar-39Ar data that post-date a major graben-bounding normal fault. These data, along with evidence for prolonged activity of normal faulting in this and other Tibetan grabens, support models that relate normal faulting to processes occurring beneath the plateau. Thinning of the upper crust is most plausibly the result of potential-energy increases resulting from spatially and temporally heterogeneous changes in thermal structure and density distribution within the crust and upper mantle beneath Tibet. This is supported by recent geophysical and geological data, which indicate that spatial heterogeneity exists in both the Tibetan crust and lithospheric mantle.