Since the Cenozoic, a series of extensional south-north normal faults and gneiss-granite domes evolved in the southern Tibetan Plateau, the formation mechanism of which is of scientific interest and which has implications for the tectonic dynamics of the plateau. Typical of such features are the Xainza-Dinggye rift and the Mabja gneiss dome, which are located in the Xainza-Xietongmen-Dinggye region in the central Tibetan Plateau. In this study, Magnetotelluric measurements across this region are used to generate a high-resolution 3-D electrical resistivity model of the subsurface and to analyze the cause of the conductive zones. The large-scale conductive zones identified in the middle-lower crust may result from aqueous melt partial melting, whereas the smaller-scale conductive zones in the upper-middle crust may result from saline fluids, possibly with varying minor volumes of melts. Subsequently, based on the electrical resistivity model and combined with the spatiotemporal coupling of the geological, geochemical, and geophysical data, the state and migration features of crustal materials are discussed. The results show that the upwelling of mantle materials along subduction channels and slab-windows related to the tearing of the Indian lithospheric plate contributed to the partial melting of the middle-lower crust in the Lhasa terrane. Furthermore, partial melting of the upper-middle crust in the Tethys-Himalaya terrane resulted from southern extrusion of crustal materials in the Lhasa terrane. These two mechanisms can significantly reduce the effective viscosity. We speculate that the deformation of the brittle upper crust that is controlled by large-scale ductile layers characterized by weak rheology is the main dynamic mechanism of rift evolution. Meanwhile, the metamorphism and anatexis in the upper-middle crust of the Tethys-Himalaya terrane related to the southern extrusion of materials contributed to the evolution of the Mabja gneiss dome. During the middle Miocene, the southern extrusion of crustal materials may have been influenced by the cooling events beneath the Mabja gneiss dome, which can explain why the deep areas beneath the Mabja gneiss dome have middle-high resistivity. In addition, our study region is located in the Mediterranean-Himalayan seismic belt, and mainly includes shallow-focus earthquakes and intermediate-depth earthquakes. In the north, shallow-focus earthquakes are mainly controlled by the accumulation of stress in the brittle layer of the overlying crust related to the ductile layer of the middle and lower crust. In the south, shallow-focus earthquakes (e.g., Dingri MS6.8 earthquake) mainly occur in the rigid, resistive block, which is surrounded by conductive zones, possibly because fluid migration may be hindered by these resistive blocks. The intermediate-depth earthquakes are mainly controlled by the response in the subsurface area, which is related to the detachment of the Indian lithospheric mantle from the Indian crust.

Both low resistivity zones and low velocity zones are distributed in the middle-lower crust of the western Lhasa terrane, Tibetan Plateau, China. Some estimates from electrical resistivity data suggest large volume fractions of silicate melts that are difficult to reconcile with seismic velocity data that prefer lower volumes. A second conductive phase, such as saline fluids, that drastically reduces the conductivity but does not significantly affect the seismic velocity because of its low volume may be able to explain these differences. In this study, a 3-D model of the electrical resistivity structure is generated on a profile along longitude 85°E from a latitude of 29°N to 32.5°N. Based on experimental measurement of melts and alkali-rich fluids (e.g., H ₂O-NaCl), we estimate the volume fraction of each phase that is required to explain the conductive anomalies observed in the geophysical model. The model reveals that the maximum bulk conductivity of the mid-lower crust in the south (1.52 S/m) is much higher than the conductivity of the mid-lower crust in the north (0.18 S/m) when taking 31°N as a rough boundary, near Coqen region. We hypothesize that the conductive zones in the south of the Coqen region may result from a silicate melt and alkali-rich fluid (multicomponent) system. In contrast, partial melting alone can explain the conductive zones in the north. The hypothesis can reconcile the predictions from electrical resistivity data and seismic data, and it corresponds well with zircon Hf isotope data. For example, a combination such as the presence of <1% NaCl-bearing aqueous fluids in addition to 5-10% partial melt can reconcile electrical conductivity data and seismic data. We propose that the contributions from partial melt or saline fluids are controlled by the distinct tectonic dynamics in each region. Furthermore, the model compatible with the idea that the Indian lower crust subducted northwards beneath the Lhasa terrane and may not extend far beyond the Indus-Yarlung Zangbo suture (approximately 30-31°N). The widespread distribution and interconnection of crustal conductors at different depths is consistent with the lateral migration of materials. However, both geophysical data sets agree that some anomalies are discontinuous along the profile. Furthermore, the low-angle subducted Indian Plate with no obvious tearing feature and a low volume of melts may have contributed to the absence of long, continuous, N-S-trending normal faults in this region.

The tectonic dynamics of the Lhasa-Gangdese terrane, southern Tibet, are still unclear and open questions persist regarding the structure, physical properties, and rheology of the lithosphere. A three-dimensional electrical resistivity model was generated beneath the Lhasa terrane, 79°−94°E and 28°–33°N, an area of ∼1,500 by ∼600 km, using data from 537 magnetotelluric measurements. To overcome computational challenges, we present an approach in which multiple high-resolution models from overlapping subregions are combined by means of an error-weighted averaging scheme to construct a continuous electrical resistivity model. Based on the electrical structure, the rheology of the mid-lower crust was investigated by constraining the melt fraction, in order to explore controls on dynamic mechanisms from a whole-system perspective. The models shed light on the vertical and lateral migration of crustal materials, magma transportation, and continental deformation in the Lhasa terrane. The model shows resistive features in the south that may represent the Indian plate, and conductive features above this that may represent the migration of materials beneath the Lhasa terrane, both of which vary substantially along the Indus-Yarlung Zangbo suture. In the mid-lower crust, conductive features (with a likely corresponding decrease of the effective viscosity) may be related to underplating, magma reservoirs, mineralization sources, and vertical motion of buoyant materials. A weakened mid-lower crust has consequences for surface deformation and may have contributed to the formation of (north-south) rift zones. Furthermore, the inhomogeneous distribution of weakened areas in the mid-lower crust has implications for the local lateral migration of materials.