Oblique stepwise rise and growth of the Tibet plateau

Hydrogeochemical mechanism of Li-Cs-rich thermal springs in the Karakoram Fault, Qinghai-Xizang Plateau: New insights from multiple stable isotopes

High lithium (Li) and cesium (Cs) concentrations in the Qinghai-Xizang Plateau thermal springs pose environmental and health challenges, but their origins and enrichment mechanisms remain unclear. This study focuses on the Sogdoi geothermal field, located along the southern Karakoram Fault, to investigate these processes. Multi-isotope analyses (H, O, Li, Sr) reveal that Li and Cs predominantly originate from the host rocks, especially granitoids and meta-sedimentary rocks, rather than from magmatic fluids. We propose a two-stage enrichment process mechanism that includes water-rock interactions and steam loss (~22 %) as water ascends. The geothermal field contains three distinct spring groups, all sourced from a common parent geothermal fluid at ~7 km depth with temperatures around 200 °C. Geochemical variations across groups reflect different degrees of conductive cooling, steam loss, and cold water dilution. Group I (Li: 45.21-51.93 mg/L, Cs: 27.37-29.39 mg/L, 87Sr/86Sr: 0.717-0.725, δ7Li: 1.45-1.74 ‰) represents deep fluids with minor dilution. Group II (Li: 32.16-39.98 mg/L, Cs: 28.72-33.91 mg/L, 87Sr/86Sr: 0.710-0.712, δ7Li: 1.64-1.79 ‰) and Group III (Li: 21.49-22.37 mg/L, Cs: 6.56-6.79 mg/L, 87Sr/86Sr: ~0.709, δ7Li: 4.22 ‰) show increasing dilution effects. Our findings indicate that Sr isotopes are highly sensitive to cold water mixing, while Li isotopes demonstrate greater resistance but are more susceptible to mineral adsorption. The binary Li-Sr isotope mixing model is more effective than single-isotope models in discerning mixing effects. It is crucial to select samples that are minimally affected by cold water intrusion to accurately trace fluid sources and their evolution. The anomalous Cs concentrations across groups are likely due to adsorption processes, as evidenced by decreasing Cs/Na ratios. We propose that deep faults are crucial in forming Li-Cs-rich thermal springs by enabling the deep circulation of water. This study provides insights into the hydrogeochemical mechanisms of Li and Cs in geothermal fluids, contributing to resource exploitation and environmental management.

B-Li differential enrichment of geothermal systems in the Da Qaidam and Gonghe-Guide Basin, northeastern Tibetan Plateau: Evidence from water chemistry and H-O-B-Li isotopes

Typical geothermal systems in the Da Qaidam (DQ) and Gonghe-Guide Basin (GGB) on the northeastern Tibetan Plateau (NETP) discharged different BLi contents. A widely accepted metallogenic model is that the salt-lake type BLi deposits in the TP are recharged by geothermal fluids with B-Li-rich, carried by rivers and enriched in the terminal salt lakes. The B-Li-rich geothermal water is the key source of mineralization in salt lakes, however, enrichment mechanism governing differential BLi contents in DQ and GGB geothermal systems remains ambiguous. This study systematically deciphers water chemistry and isotope characteristics (δD, δ18O, δ11B, δ7Li) of river waters, geothermal waters, sinters and surrounding rocks to discuss the enrichment process of BLi elements in the DQ and GGB geothermal systems. The δD and δ18O values of geothermal systems in the DQ (δD = -79.60 ∼ -82.40 ‰, δ18O = -10.97 ∼ -11.38 ‰) and GGB (δD = -64.00 ∼ -97.10 ‰, δ18O = -8.70 ∼ -13.00 ‰) are close to the GMWL and LMWL, indicating meteoric origin. The δD and high Cl- values (300-900 mg/L) of geothermal waters in the DQ and Qiabuqia, Qunaihai, Zhacangsi along Riyueshan of the GGB imply that these geothermal waters mixed by magmatic fluid and meteoric water. The hot springs in the DQ (B = 38.35-46.29 mg/L, Li = 3.11-3.72 mg/L) and GGB (B = 0.17-8.16 mg/L, Li = 0.08-10.49 mg/L) exhibit differential BLi geochemical characteristics. B and Li contents are higher in DQ hot springs and in hot springs along Riyueshan of the GGB, respectively. Comparison of the BLi contents and δ11B-δ7Li values of geothermal waters, sinters and host rocks reveals that BLi contents of geothermal waters are controlled by B-rich HP-UHP metamorphic rocks formed by metamorphism and Li-rich granites or pegmatites formed by magmatism in the Qilian Shan, respectively. Moreover, metamorphic and magmatic processes, combining with deep circulation reactivated by the thrust or strike-slip faults, create differential enrichment of BLi elements in the geothermal systems, such as DQ and GGB on the NETP. This study highlights into understanding the differential enrichment of BLi in the geothermal system on the TP. Furthermore, the resource elemental abundance of geothermal waters can be applied as an important prospecting indicator of endogenous BLi deposits.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-Transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

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.

Paleogeography of the southwestern Ordos Basin and exhumation history of the Liupan Shan area

The Liupan Shan, situated on the southwestern Ordos Basin, underwent Cenozoic uplift caused by the India-Asia collision and subsequent northeastward tectonic movements. The strata in this region record both the paleogeography of the southwestern Ordos Basin and the uplift history of the Tibetan Plateau. However, past studies have rarely analyzed the strata comprehensively, resulting in overlooked information within them. We present the first detrital zircon (U-Th)/He (ZHe) data from the Lower Cretaceous deposits of the Liupan Shan and the Upper Miocene red clay of Chaona, constraining a maximum burial depth of less than 6-7 km at 60 Ma for the Liupan area. By integrating zircon U-Pb ages, paleocurrent data, and sedimentary facies, we found a primary ZHe age peak (∼210 Ma) indicating that the Lower Cretaceous deposits were sourced from the Qinling Orogenic Belt, while the red clay was likely originated from the Liupan Shan. Although the Chaona strata are not preserved in the Liupan Shan, our results suggest that the Miocene Chaona section is crucial for understanding the Late Cenozoic uplift of the Liupan Shan. The distinct stratigraphic ages and source areas reflect the complex depositional and tectonic history of the region. Thermal modeling results revealed three stages of tectonic events (Late Cretaceous, Late Paleocene-Early Eocene, and Late Cenozoic), which are linked to plate interaction. Our findings offer new insights into the long-term tectonic evolution of the Tibetan Plateau.

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.

Variability in interseismic strain accumulation rate and style along the Altyn Tagh Fault

Major strike-slip faults that develop between strong and weaker regions are thought to focus along narrow shear zones at the rheological boundary. Here we present the InSAR-derived velocity field spanning almost the entire length of one such fault, the 1600 km-long Altyn Tagh Fault (ATF), and analyse the strain distribution. We find that localisation of strain is actually variable, in contrast to other major strike-slip faults that show little variation, with strain concentrated at the fault for some sections and distributed over broad (>100 km) shear zones for others. Slip rate along the ATF is also variable, decreasing along the fault from 11.6 ± 1.6 mm/yr in the west to 7.2 ± 1.4 mm/yr in the central portion, before increasing again to 11.7  ± 0.9 mm/yr over the eastern portion. We show that the variable shear zone width may be linked to geological variability and the influence of heat flow, and the results imply that sub-parallel faults play an important role in the overall deformation field. This demonstrates the significance of accurately characterising strain rates over a broad region when assessing seismic hazard.

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae)

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Impact of the 2008 M W 7.9 Great Wenchuan earthquake on South China microplate motion

Tectonic plate motions drive the earthquake cycle, as they result in the slow accrual and sudden release of energy along plate boundaries. Steadiness of plate motions over the earthquake cycle is a central tenet of the plate tectonics theory and has long been a main pillar in models of earthquake genesis, or of plate-margins seismic potential inferred from slip-deficit estimates. The advent of geodesy in the geosciences and the availability of multi-year-long series of position measurements permit tracking the motions of tectonic plates from before to after the time of significant seismic events that occur along their margins. Here, we present evidence that large earthquakes are capable of modifying the motions of entire microplates. We use high precision Global Navigation Satellite System (GNSS) position time-series covering the periods 2001-2004 and 2014-2017 to demonstrate that, contrary to the tenet above, the South China microplate motion changed after the 2008 M W 7.9 Great Wenchuan earthquake. The GNSS data and associated uncertainties indicate a plate motion slowdown of up to 20% that is beyond the possible impact of data noise and is thus tectonically meaningful. We use quantitative models of torque balance to show that generating this kinematic change requires a force upon the South China microplate compatible with that imparted by the Great Wenchuan earthquake of 2008. The existence of a kinematic signal linked to the earthquake cycle that impacts an entire microplate might offer an additional, novel perspective to assessing the hazards of earthquake-prone tectonic regions.

Uplift of the Tibetan Plateau driven by mantle delamination from the overriding plate

The geodynamic evolution of the Tibetan Plateau remains highly debated. Any model of its evolution must explain the plateau's growth as constrained by palaeo-altitude studies, the spatio-temporal distribution of magmatic activity, and the lithospheric mantle removal inferred from seismic velocity anomalies in the underlying mantle. Several conflicting models have been proposed, but none of these explains the first-order topographic, magmatic and seismic features self-consistently. Here we propose and test numerically an evolutionary model of the plateau that involves gradual peeling of the lithospheric mantle from the overriding plate and consequent mantle and crustal melting and uplift. We show that this model successfully reproduces the successive surface uplift of the plateau to more than 4 km above sea level and is consistent with the observed migration of magmatism and geometry of the lithosphere-asthenosphere boundary resulting from subduction of the Indian plate and delamination of the mantle lithosphere of the Eurasian plate. These comparisons indicate that mantle delamination from the overriding plate is the driving force behind the uplift of the Tibetan Plateau and, potentially, orogenic plateaus more generally.

REVIEW: Evidence supporting the 'preparation for oxidative stress' (POS) strategy in animals in their natural environment

Hypometabolism is a common strategy employed by resilient species to withstand environmental stressors that would be life-threatening for other organisms. Under conditions such as hypoxia/anoxia, temperature and salinity stress, or seasonal changes (e.g. hibernation, estivation), stress-tolerant species down-regulate pathways to decrease energy expenditures until the return of less challenging conditions. However, it is with the return of these more favorable conditions and the reactivation of basal metabolic rates that a strong increase of reactive oxygen and nitrogen species (RONS) occurs, leading to oxidative stress. Over the last few decades, cases of species capable of enhancing antioxidant defenses during hypometabolic states have been reported across taxa and in response to a variety of stressors. Interpreted as an adaptive mechanism to counteract RONS formation during tissue hypometabolism and reactivation, this strategy was coined "Preparation for Oxidative Stress" (POS). Laboratory experiments have confirmed that over 100 species, spanning 9 animal phyla, apply this strategy to endure harsh environments. However, the challenge remains to confirm its occurrence in the natural environment and its wide applicability as a key survival element, through controlled experimentation in field and in natural conditions. Under such conditions, numerous confounding factors may complicate data interpretation, but this remains the only approach to provide an integrative look at the evolutionary aspects of ecophysiological adaptations. In this review, we provide an overview of representative cases where the POS strategy has been demonstrated among diverse species in natural environmental conditions, discussing the strengths and weaknesses of these results and conclusions.

Source Models of the 2016 and 2022 Menyuan Earthquakes and Their Tectonic Implications Revealed by InSAR

The Haiyuan fault system plays a crucial role in accommodating the eastward expansion of the Tibetan Plateau (TP) and is currently slipping at a rate of several centimeters per year. However, limited seismic activities have been observed using geodetic techniques in this area, impeding the comprehensive investigation into regional tectonics. In this study, the geometric structure and source models of the 2022 Mw 6.7 and the 2016 Mw 5.9 Menyuan earthquakes were investigated using Sentinel-1A SAR images. By implementing an atmospheric error correction method, the signal-to-noise ratio of the 2016 interferometric synthetic aperture radar (InSAR) coseismic deformation field was significantly improved, enabling InSAR observations with higher accuracy. The results showed that the reliability of the source models for those events was improved following the reduction in observation errors. The Coulomb stress resulting from the 2016 event may have promoted the strike-slip movement of the western segment of the Lenglongling fault zone, potentially expediting the occurrence of the 2022 earthquake. The coseismic slip distribution and the spatial distribution of aftershocks of the 2022 event suggested that the seismogenic fault may connect the western segment of the Lenglongling fault (LLLF) and the eastern segment of the Tuolaishan fault (TLSF). Additionally, the western segment of the surface rupture zone of the northern branch may terminate in the secondary branch close to the Sunan-Qilian fault (SN-QL) strike direction, and the earthquake may have triggered deep aftershocks and accelerated stress release within the deep seismogenic fault.

Leaf ecological stoichiometry and anatomical structural adaptation mechanisms of Quercus sect. Heterobalanus in southeastern Qinghai-Tibet Plateau

BACKGROUND

With the dramatic uplift of the Qinghai-Tibet Plateau (QTP) and the increase in altitude in the Pliocene, the environment became dry and cold, thermophilous plants that originally inhabited ancient subtropical forest essentially disappeared. However, Quercus sect. Heterobalanus (QSH) have gradually become dominant or constructive species distributed on harsh sites in the Hengduan Mountains range in southeastern QTP, Southwest China. Ecological stoichiometry reveals the survival strategies plants adopt to adapt to changing environment by quantifying the proportions and relationships of elements in plants. Simultaneously, as the most sensitive organs of plants to their environment, the structure of leaves reflects of the long-term adaptability of plants to their surrounding environments. Therefore, ecological adaptation mechanisms related to ecological stoichiometry and leaf anatomical structure of QSH were explored. In this study, stoichiometric characteristics were determined by measuring leaf carbon (C), nitrogen (N), and phosphorus (P) contents, and morphological adaptations were determined by examining leaf anatomical traits with microscopy.

RESULTS

Different QSH life forms and species had different nutrient allocation strategies. Leaves of QSH plants had higher C and P and lower N contents and higher N and lower P utilization efficiencies. According to an N: P ratio threshold, the growth of QSH species was limited by N, except that of Q. aquifolioides and Q. longispica, which was limited by both N and P. Although stoichiometric homeostasis of C, N, and P and C: N, C: P, and N: P ratios differed slightly across life forms and species, the overall degree of homeostasis was strong, with strictly homeostatic, homeostatic, and weakly homeostatic regulation. In addition, QSH leaves had compound epidermis, thick cuticle, developed palisade tissue and spongy tissue. However, leaves were relatively thin overall, possibly due to leaf leathering and lignification, which is strategy to resist stress from UV radiation, drought, and frost. Furthermore, contents of C, N, and P and stoichiometric ratios were significantly correlated with leaf anatomical traits.

CONCLUSIONS

QSH adapt to the plateau environment by adjusting the content and utilization efficiencies of C, N, and P elements. Strong stoichiometric homeostasis of QSH was likely a strategy to mitigate nutrient limitation. The unique leaf structure of the compound epidermis, thick cuticle, well-developed palisade tissue and spongy tissue is another adaptive mechanism for QSH to survive in the plateau environment. The anatomical adaptations and nutrient utilization strategies of QSH may have coevolved during long-term succession over millions of years.

Characteristics and Geological Significance of Lacustrine Hydrothermal Sedimentary Rocks in the Yingejing Sag in Bayingebi Basin, Inner Mongolia, Northwestern China

In order to promote site screening for high-level radioactive waste (HLW) disposal purposes, the characteristics of argillaceous rock (potential host rock) from the Yingejing Sag of the Bayingebi Basin, Northwest China have been well discussed. Results show that (1) Unlike argillaceous host rocks in foreign countries, the argillaceous rock mainly consists of analcite, dolomite, and albite; the contents of clay minerals are only about 10%. Five typical structures could be categorized, dominating by the massive structure. (2) Geochemical characteristics have the characteristics of abundance in deep source gas and liquid trace elements, a gentle right dip in the distribution pattern of rare-earth elements. (3) Petrological and geochemical characteristics determine the argillaceous rock as the genesis of lacustrine hydrothermal sedimentary rock. The hydrothermal sedimentary model also has been constructed, mainly controlled by tectonic activity of the Altyn Tagh fault from 100 to 120 Ma. The massive argillites with analcite and dolomite would lay the foundation for confirming the site for HLW disposal purposes in China.

Identifying the source and fate of boron in geothermal water: Evidence from B/Na and B isotopes

High level dissolved B, which poses risks to human health, has been widely observed in geothermal water. In the Guide Basin, NW China, a series of geothermal water samples along a fault show a wide range of B contents ranging from 3.14 to 8.33 mg/L, which are higher than the WHO Guideline value equaling 2.4 mg/L in drinking water. To identify the sources and fate of B, we conduct a comprehensive analysis of hydrochemistry and stable isotopes (D, 18O and 11B) of three thermal fields representing three stages of hydrogeochemical evolution (stages I, II and III). From stage I to III, there are trends of increasing mineral dissolution, which is supported by increasing mean reservoir temperature and concentrations of conservative elements (Cl, Na, K, Li and Si). Geothermal water in stage I with meteoric origin and the lowest reservoir temperature has the highest B/Na resulting from silicate dissolution and falls on the mixing line between granitoids and cold water on the plot of δ11B versus 1/B, showing the control of silicate dissolution. However, geothermal water in stage III has lower Ca, B Sr and B/Na than that in stage II. Because of the occurrence of other processes, geothermal water in stages II and III deviates from the LMWL. Compared with geothermal water in stage I, the increased Sr/Ca and decreased B/Ca show that B are removed by both coprecipitation and vapor separation. With the aid of B isotopes, we find vapor separation dominates in stage II, whereas carbonate precipitation dominates in stage III. Overall, a combined use of three isotopes (H, O and B) and three element ratios (B/Na, B/Ca and Sr/Ca) leads to a complete understanding of B cycle and hydrogeochemical evolution in hydrothermal systems.

Geodynamic models of Indian continental flat slab subduction with implications for the topography of the Himalaya-Tibet region

The slab structure and high elevation of the Himalaya-Tibet region and their underlying mechanisms have been widely discussed. Many studies interpret a flat slab segment of Indian continental lithosphere located below the overriding plate, but interpretations of the northward extent of the flat slab differ substantially, with minimum estimates placing the boundary at the northern margin of the Himalaya (Indus-Yarlung Tsangpo suture), and maximum estimates placing it at the northern boundary of Tibet. In this study, we investigate for the first time if a flat slab segment of subducted buoyant Indian continental lithosphere below the Himalaya-Tibet region is geodynamically feasible and we quantify its northward extent, as well as its contribution to the high topography of the region. We conduct three large-scale fully-dynamic (buoyancy-driven) analogue experiments to simulate the subduction of the Indian continent. Our preferred, and geodynamically most feasible, model shows a continental flat slab extending northward up to ~ 320 km from the Himalayan thrust front, in agreement with recent estimates. Furthermore, it suggests that the positively buoyant flat slab segment of the Indian continent contributes some ~ 1.5-2 km to the high topography of the Himalaya-Southern Tibet region by providing an upward force to elevate the overriding Eurasian plate.

Genetic Diversity and Phylogenetic Analysis of Zygophyllum loczyi in Northwest China's Deserts Based on the Resequencing of the Genome

In order to study the genetics of local adaptation in all main deserts of northwest China, whole genomes of 169 individuals were resequenced, which covers 20 populations of Zygophyllum loczyi (Zygophyllales: Zygophylaceae). We describe more than 15 million single nucleotide polymorphisms and numerous InDels. The expected heterozygosity and PIC values associated with local adaptation varied significantly across biogeographic regions. Variation in environmental factors contributes largely to the population genetic structure of Z. loczyi. Bayesian analysis performed with STRUCTURE defined four genetic clusters, while the results of principle component analysis were similar. Our results shows that the Qaidam Desert group appears to be diverging into two branches characterized by significant geographic separation and gene flow with two neighboring deserts. Geological data assume that it is possible that the Taklamakan Desert was the original distribution site, and Z. loczyi could have migrated later on and expanded within other desert areas. The above findings provide insights into the processes involved in biogeography, phylogeny, and differentiation within the northwest deserts of China.

Analyzing the horizontal orientation of the crustal stress adjacent to plate boundaries

The spatial analysis of horizontal stress orientation is important to study stress sources and understand tectonics and the deformation of the lithosphere. Additional to the stress sources, the geometry of stress fields depends on the underlying coordinate reference system, which causes spatial distortions that bias the analysis and interpretation of stresses. The bias can be avoided when the stress field is decomposed and transformed into the reference frame of its first-order stress source. We present a modified and extended theory based on the empirical link between the orientation of first-order stresses and the trajectories of lateral plate boundary forces. This link is applied to analyze the orientation of horizontal stresses, their patterns, and tectonic structures from the perspective of their first-order source or cause. By using only parameters for the relative motion between two neighboring plates, we model the first-order orientation of the maximum horizontal stress that statistically fits the orientation of [Formula: see text]80% of the global stress data adjacent to plate boundaries. Considerable deviations of the observed stress from the predicted first-order stress direction can reveal the geometry of second-order stresses and confine areas where other stress sources dominate. The model's simple assumptions, independence from the sample size, potential application to regional to global scale analysis, and compatibility with other spatial interpolation algorithms make it a powerful method for analyzing stress fields. For immediate use, the presented method is implemented in the free and open-source software package tectonicr, which is written in the computer language R.

Molecular phylogeny of mega-diverse Carabus attests late Miocene evolution of alpine environments in the Himalayan-Tibetan Orogen

The timing, sequence, and scale of uplift of the Himalayan-Tibetan Orogen (HTO) are controversially debated. Many geoscientific studies assume paleoelevations close to present-day elevations and the existence of alpine environments across the HTO already in the late Paleogene, contradicting fossil data. Using molecular genetic data of ground beetles, we aim to reconstruct the paleoenvironmental history of the HTO, focusing on its southern margin (Himalayas, South Tibet). Based on a comprehensive sampling of extratropical Carabus, and ~ 10,000 bp of mitochondrial and nuclear DNA we applied Bayesian and Maximum likelihood methods to infer the phylogenetic relationships. We show that Carabus arrived in the HTO at the Oligocene-Miocene boundary. During the early Miocene, five lineages diversified in different parts of the HTO, initially in its southern center and on its eastern margin. Evolution of alpine taxa occurred during the late Miocene. There were apparently no habitats for Carabus before the late Oligocene. Until the Late Oligocene elevations must have been low throughout the HTO. Temperate forests emerged in South Tibet in the late Oligocene at the earliest. Alpine environments developed in the HTO from the late Miocene and, in large scale, during the Pliocene-Quaternary. Findings are consistent with fossil records but contrast with uplift models recovered from stable isotope paleoaltimetry.

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.

Climate heterogeneity shapes phylogeographic pattern of Hippophae gyantsensis (Elaeagnaceae) in the east Himalaya-Hengduan Mountains

The interaction of recent orographic uplift and climate heterogeneity acted as a key role in the East Himalaya-Hengduan Mountains (EHHM) has been reported in many studies. However, how exactly the interaction promotes clade diversification remains poorly understood. In this study, we both used the chloroplast trnT-trnF region and 11 nuclear microsatellite loci to investigate the phylogeographic structure and population dynamics of Hippophae gyantsensis and estimate what role geological barriers or ecological factors play in the spatial genetic structure. The results showed that this species had a strong east-west phylogeographic structure, with several mixed populations identified from microsatellite data in central location. The intraspecies divergence time was estimated to be about 3.59 Ma, corresponding well with the recent uplift of the Tibetan Plateau. Between the two lineages, there was significant climatic differentiation without geographic barriers. High consistency between lineage divergence, climatic heterogeneity, and Qingzang Movement demonstrated that climatic heterogeneity but not geographic isolation drives the divergence of H. gyantsensis, and the recent regional uplift of the QTP, as the Himalayas, creates heterogeneous climates by affecting the flow of the Indian monsoon. The east group of H. gyantsensis experienced population expansion c. 0.12 Ma, closely associated with the last interglacial interval. Subsequently, a genetic admixture event between east and west groups happened at 26.90 ka, a period corresponding to the warm inter-glaciation again. These findings highlight the importance of the Quaternary climatic fluctuations in the recent evolutionary history of H. gyantsensis. Our study will improve the understanding of the history and mechanisms of biodiversity accumulation in the EHHM region.

Evidence for structural control of mare volcanism in lunar compressional tectonic settings

One of the long-standing enigmas for lunar tectonic-thermal evolution is the spatiotemporal association of contractional wrinkle ridges and basaltic volcanism in a compressional regime. Here, we show that most of the 30 investigated volcanic (eruptive) centers are linked to contractional wrinkle ridges developed above preexisting basin basement-involved ring/rim normal faults. Based on the tectonic patterns associated with the basin formation and mass loading and considering that during the subsequent compression the stress was not purely isotropic, we hypothesize that tectonic inversion produced not only thrust faults but also reactivated structures with strike-slip and even extensional components, thus providing a valid mechanism for magma transport through fault planes during ridge faulting and folding of basaltic layers. Our findings suggest that lunar syn-tectonic mare emplacement along reactivated inherited faults provides important records of basin-scale structure-involved volcanism, which is more complex than previously considered.

The constancy of chemical weathering intensity on hillslopes in the arid to semiarid Qilian Mountains, NE Tibetan Plateau

Quantifying the relationship between the chemical weathering and denudation rates of active orogenic belts over a range of climates is the key to addressing the controversy over the uplift-weathering hypothesis. However, studies have focused on warm and humid environments and have not examined cold and arid environments. Here, we present a new dataset of the chemical depletion fraction (CDF: ratio of the chemical weathering rate to the total denudation rate) across the arid to semiarid Qilian Mountains on the northeastern Tibetan Plateau, where the uplift-weathering hypothesis has been proposed. We selected 60 points from 12 catchments in the middle Qilian Mountains. At each point, we collected three samples (soil, saprolite, and bedrock samples) and calculated the CDF values based on their Zr concentrations. We found no clear correlation between the CDF and climatic factors (temperature, precipitation, and normalized difference vegetation index (NDVI)), topographic factors (slope and local relief), and denudation rate. The calculated chemical weathering rates, nevertheless, are positively correlated with precipitation, NDVI, and denudation rate, and negatively correlated with temperature. This result indicates that the Qilian Mountains are under supply-limited conditions, even at high denudation rates (>800 t km^-2 y^-1). We speculate that low temperatures (<0 °C) could intensify near-surface chemical weathering by promoting the physical breakdown of the bedrock and increasing soil water availability. This mechanism causes a compensation effect maintaining the supply-limited conditions in landscapes with high denudation rates. Combing a worldwide dataset regarding the correlations between CDF and climatic factors and denudation rates, we argue that intensified denudation since the late Cenozoic contributed to global cooling.

Distinct lake sedimentary imprints of earthquakes, floods and human activities in the Xiaojiang Fault zone: Towards a quantitative paleoseismograph in the southeastern Tibetan Plateau

Lake sediments that widely distributed in the active and complicated fault zones have been recently showing great potential for paleoseismic reconstruction. However, flood events and human activities may make the seismic signal unrecognizable. In this study, high-resolution analyses of sedimentary structure, physical and chemical proxies, as well as absolutely radioactive dating were conducted on seven representative sediment cores from the depocenter, nearshore and inlet areas of Yangzong Lake, a typical fault lake in the Xiaojiang Fault zone, southeastern Tibetan Plateau (TP). These new data were calibrated by historical documents, suggesting that seismically induced mass-transport deposits (MTDs, i.e., turbidites) were massive and/or amalgamated (earthquake doublet), became fining and thickening towards the lake center (without changing lake morphology), and occasionally exhibited soft sediment deformation structures (SSDs, i.e., microfaults). These sediments were relatively poorly sorted and instantaneously deposited from slope failures within the lake. An extremely strong earthquake could cause coseismic subsidence of the lake basin and destruct the local hydrological system, resulting in exceptionally high Mn and total inorganic carbon (TIC) contents in the lake center. In contrast, flood deposits were thinner with horizontal beddings, had higher terrestrial organic matter (higher C/N ratios), and distributed locally in the lake inlet area. Human activities-induced sediments were inversely graded, poorly sorted and gradually deposited, had horizontal beddings and no erosive base, and exhibited high carbon, Pb and Zn contents and low C/N ratios. In addition, macroseismic investigations and statistical results from intensity prediction equations (IPEs) provided a conservative threshold of ∼8 Modified Mercalli Intensities (MMI) for triggering turbidites, and a ∼ 10 MMI for inducing coseismic subsidence and hydrological destruction. This study was among the first attempts to establish a quantitative lacustrine paleoseismograph in the southeastern TP, and the new results would greatly improve the valid assessment of geohazard risks.

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.

Revisiting the mechanisms of mid-Tertiary uplift of the NE Tibetan Plateau

Contrasting views exist on timing and mechanisms of Tertiary crustal uplift in the NE Tibetan Plateau based on different approaches, with many models attributing surface uplift to crustal shortening. We carry out a comprehensive investigation of mid-Tertiary stratigraphy, sedimentology, and volcanism in the West Qinling, Hoh Xil and Qaidam basin, and the results challenge previous views. It was held that the discordance between Oligocene and Miocene strata is an angular unconformity in the West Qinling, but our field observations show that it is actually a disconformity, indicative of vertical crustal uplifting rather than crustal shortening at the Oligocene to Miocene transition. Widespread occurrence of synsedimentary normal faults in mid-Tertiary successions implicates supracrustal stretching. Miocene potassic-ultrapassic and mafic-ultramafic volcanics in the Hoh Xil and West Qinling suggest a crucial role of deep thermomechanical processes in generating crust- and mantle-sourced magmatism. Also noticeable are the continuity of mid-Tertiary successions and absence of volcanics in the Qaidam basin. Based on a holistic assessment of stratigraphic-sedimentary processes, volcanic petrogenesis, and spatial variations of lithospheric thicknesses, we speculate that small-sale mantle convection might have been operating beneath northeast Tibet in the mid-Tertiary. It is assumed that northward asthenospheric flow was impeded by thicker cratonic lithosphere of the Qaidam and Alxa blocks, thereby leading to edge convection. The edge-driven convection could bring about surface uplift, induce supracrustal stretching, and trigger vigorous volcanism in the Hoh Xil and West Qinling in the mid-Tertiary period. This mechanism satisfactorily explains many key geologic phenomena that are hardly reconciled by previous models.

Leaf morphology, functional trait and altitude response in perennial vetch (Vicia unijuga A. Braun), alfalfa (Medicago sativa L.) and sainfoin (Onobrychis viciifolia Scop.)

Species have plasticity across altitude gradients in leaf morphology and function, and their response to high altitude conditions was mainly reflected in leaf cell metabolism and gas exchange. Leaf morphological and functional adaptation to altitude has received research attention in recent years, but there are no studies for forage legumes. Here we report differences in 39 leaf morphology and functional traits of three leguminous forages (alfalfa, sainfoin and perennial vetch) at three sites in Gansu Province, China, ranging from 1768 to 3074 m altitude to provide information for potential use in breeding programmes. With increasing altitude, plant water status increased, reflecting increase in soil water content and decreased average temperature, which lead to leaf intercellular CO₂ concentration. Stomatal conductance and evapotranspiration increased significantly but water-use efficiency decreased. At high altitude, Φ_PSII decreased but non-photochemical quenching and chlorophyll a:b ratio increased while spongy mesophyll tissue and leaf thickness increased. These changes may be due to UV or low-temperature damage of leaf protein and metabolic cost of plant protection or defence responses. Contrary to many other studies, leaf mass per area decreased significantly at higher altitude. This was consistent with predictions under the worldwide leaf economic spectrum on the basis that soil nutrients increased with increasing altitude. The key species differences were more irregularly shaped epidermal cells and larger stomatal size in perennial vetch compared to alfalfa or sainfoin that enhanced gas exchange and photosynthesis by generating mechanical force, increasing guard cell turgor, and promoting stomatal operation. The lower adaxial stomatal density also enhanced water-use efficiency. These adaptations might confer perennial vetch an advantage in environments with extreme diurnal temperature fluctuation or in frigid conditions.

A new biologic paleoaltimetry indicating Late Miocene rapid uplift of northern Tibet Plateau

The uplift of the Tibet Plateau (TP) during the Miocene is crucial to understanding the evolution of Asian monsoon regimes and alpine biodiversity. However, the northern Tibet Plateau (NTP) remains poorly investigated. We use pollen records of montane conifers ( Tsuga , Podocarpus , Abies , and Picea ) as a new paleoaltimetry to construct two parallel midrange paleoelevation sequences in the NTP at 1332 ± 189 m and 433 ± 189 m, respectively, during the Middle Miocene [~15 million years ago (Ma)]. Both midranges increased rapidly to 3685 ± 87 m in the Late Miocene (~11 Ma) in the east, and to 3589 ± 62 m at ~7 Ma in the west. Our estimated rises in the east and west parts of the NTP during 15 to 7 Ma, together with data from other TP regions, indicate that during the Late Miocene the entire plateau may have reached a high elevation close to that of today, with consequent impacts on atmospheric precipitation and alpine biodiversity.

The stepwise Indian-Eurasian collision and uplift of the Himalayan-Tibetan plateau drove the diversification of high-elevation Scytodes spiders

The Cenozoic Indian-Eurasian collision and uplift of the Himalayan-Tibetan Plateau (HTP) are among the most important geological events in the world. They have affected the diversification of regional biota of many taxonomic groups on and around the HTP. However, the exact timing and model of the collision and uplift events and speciation on and around the HTP are still in debate. The Himalayas group of Scytodes spitting spiders (Araneae: Scytodidae) are distributed at high elevations of the HTP and northern Indochina. Here, we reconstruct a dated molecular phylogeny for pan-Himalayan Scytodes spiders, including the Himalayas group, with full geographical sampling of the species from the HTP and Indochina. We test a hypothesis to explain that the rich montane biodiversity of the region is uplift-driven diversification-that orogeny drives rapid in situ speciation of the resident Scytodes lineages. Our findings revealed that the separation of the Himalayas clade from the Myanmar clade took place during the middle Oligocene, reflecting the final collision of India with Eurasia. The deep divergences among three clades (the Himalayas, the Myanmar and the Indochina clades) occurred from the middle Eocene to the middle Oligocene, corresponding to two early uplift events of the HTP. The evolutionary split between the Himalayas + Myanmar and Indochina clades were simultaneous with the rapid lateral extrusion of Indochina by the initial Himalayan uplift around the Eocene. This study highlights the importance of the diversification of dispersal-limited, high-elevation invertebrates as independent lines of evidence to reflect key tectonic events in the Himalayan-Indochina region, supporting the stepwise model for the Indian-Eurasian collision and uplift of HTP.

An overlooked dispersal route of Cardueae (Asteraceae) from the Mediterranean to East Asia revealed by phylogenomic and biogeographical analyses of Atractylodes

BACKGROUND AND AIMS

The East Asian-Tethyan disjunction pattern and its mechanisms of formation have long been of interest to researchers. Here, we studied the biogeographical history of Asteraceae tribe Cardueae, with a particular focus on the temperate East Asian genus Atractylodes DC., to understand the role of tectonic and climatic events in driving the diversification and disjunctions of the genus.

METHODS

A total of 76 samples of Atractylodes from 36 locations were collected for RAD-sequencing. Three single nucleotide polymorphism (SNP) datasets based on different filtering strategies were used for phylogenetic analyses. Molecular dating and ancestral distribution reconstruction were performed using both chloroplast DNA sequences (127 Cardueae samples) and SNP (36 Atractylodes samples) datasets.

KEY RESULTS

Six species of Atractylodes were well resolved as individually monophyletic, although some introgression was identified among accessions of A. chinensis, A. lancea and A. koreana. Dispersal of the subtribe Carlininae from the Mediterranean to East Asia occurred after divergence between Atractylodes and Carlina L. + Atractylis L. + Thevenotia DC. at ~31.57 Ma, resulting in an East Asian-Tethyan disjunction. Diversification of Atractylodes in East Asia mainly occurred from the Late Miocene to the Early Pleistocene.

CONCLUSIONS

Aridification of Asia and the closure of the Turgai Strait in the Late Oligocene promoted the dispersal of Cardueae from the Mediterranean to East China. Subsequent uplift of the Qinghai-Tibet Plateau as well as changes in Asian monsoon systems resulted in an East Asian-Tethyan disjunction between Atractylodes and Carlina + Atractylis + Thevenotia. In addition, Late Miocene to Quaternary climates and sea level fluctuations played major roles in the diversification of Atractylodes. Through this study of different taxonomic levels using genomic data, we have revealed an overlooked dispersal route between the Mediterranean and far East Asia (Japan/Korea) via Central Asia and East China.

Does Large-Scale Crustal Flow Shape the Eastern Margin of the Tibetan Plateau? Insights From Episodic Magmatism of Gongga-Zheduo Granitic Massif

The mechanisms driving crustal deformation and uplift of orogenic plateaus are fundamental to continental tectonics. Large-scale crustal flow has been hypothesized to occur in eastern Tibet, but it remains controversial due to a lack of geologic evidence. Geochemical and isotopic data from Cenozoic igneous rocks in the eastern Tibet-Gongga-Zheduo intrusive massif, provide a way to test this model. Modeling results suggest that Cenozoic magmas originated at depths of ∼30-40 km, the depth that crustal flow has been postulated to occur at. Detailed isotopic analyses indicate that the igneous rocks are derived from partial melting of the local Songpan-Ganzi crust, arguing against a long-distance crustal flow. Episodic magmatism during the Cenozoic showing a repeated shifting of magmatic sources can be correlated with crustal uplift. The continued indentation of the Indian Block and upwelling of the asthenosphere contribute to the crustal deformation, magmatism, and uplift.

Crustal Electrical Structure of the Ganzi Fault on the Eastern Tibetan Plateau: Implications for the Role of Fluids in Earthquakes

The initiation and evolution of seismic activity in intraplate regions are controlled by heterogeneous stress and highly fractured rocks within the rock mass triggered by fluid migration. In this study, we imaged the electrical structure of the crust beneath the Ganzi fault using a three-dimensional magnetotelluric inversion technique, which is host to an assemblage of resistive and conductive features extending into the lower crust. It presents a near-vertical low-resistance zone that cuts through the brittle ductile transition zone, extends to the lower crust, and acts as a pathway for fluid migration from the crustal flow to the upper crustal depths. Conductors in the upper and lower crust are associated with saline fluids and 7% to 16% partial melting, respectively. The relationship between the earthquake epicenter and the surrounding electrical structure suggests that the intraplate seismicity is triggered by overpressure fluids, which are dependent on fluid volume changes generated by the decompression dehydration of partially molten material during upwelling and native fluid within the crustal flow.

The Southwestern Boundary of Cenozoic Qaidam Basin: Constraints from Heavy Mineral Analysis

The formation of the Qaidam Basin plays an important role in unraveling the growth history of the Tibetan plateau. An extraordinary thick Cenozoic sedimentary succession of the Qaidam Basin is a great contributor to the study of the basin’s evolution history. To date, there has been disagreement on the southwestern boundary of the Paleogene Qaidam Basin. In this study, the method of heavy mineral analysis was adopted to reconstruct the southwestern boundary of the Qaidam Basin. The stable heavy minerals which represent the maturity of detrital sediments can roughly reflect the distance between the source and the deposit area. Therefore, the isogram of the stable heavy mineral index (ZTR = 20) was compiled to infer the location of the source area of the southwestern Qaidam Basin. The isogram shows that the boundary of the southwestern Qaidam Basin stretched southwesterly to the present-day Qiman Tagh Eastern Kunlun Mountains during the Paleogene. Additionally, the isolines present a remarkable northward migration since the late Eocene, which indicates the boundary of the Cenozoic Qaidam Basin that withdrew northward since the late Eocene. The specific location of the southern source area of the Qaidam Basin can be deduced at the Adatan fault, the middle of the present-day Eastern Kunlun Mountains. This result also supports the idea that the Qaidam Basin was an independent basin during the early Cenozoic era, and the Eastern Kunlun Mountains have already been exhumed during that time, serving as a prominent source of clastic sediments in the southwestern Qaidam Basin.

Rupture process of the 2021 M7.4 Maduo earthquake and implication for deformation mode of the Songpan-Ganzi terrane in Tibetan Plateau

The deformation mode of the Tibetan Plateau is of crucial importance for understanding its construction and extrusion processes, as well as for the assessment of regional earthquake potential. Block motion and viscous flow models have been proposed to describe the deformation field but are not fully supported by modern geophysical observations. The 2021 Mw 7.4 Maduo earthquake, which occurred inside the Songpan-Ganzi terrane (SGT) in central-east Tibet, provides a chance to evaluate the associated deformation mode of the region. We conduct a joint inversion for this earthquake and resolve a bilateral rupture process, which is characterized by super- and subshear rupture velocities, respectively. We interpret this distinct rupture behavior to be the result of the respective slip concentration depths of the two ruptured segments. We analyze geological, seismic, and geodetic evidence and find that the SGT upper crust shows distributed shear deformation and distinct transverse anisotropy, which are associated with folded structures originating from compression of the paleo-Tethys ocean accretional prism realigned by following shear deformation. The SGT receives lateral shear loading from its NS boundary and accommodates a right-step sinistral motion across the terrane boundary faults. The unique tectonic setting of the SGT defines locations and behaviors of internal faulting and strong earthquakes such as the 2021 Maduo earthquake, with the latter occurring on slow-moving faults at intervals of several thousands of years.

Strain Field Features and Three-Dimensional Crustal Deformations Constrained by Dense GRACE and GPS Measurements in NE Tibet

The continuing impact between the Eurasia Plate and India results in the thickening and shortening of the N-S Tibetan Plateau. There has been strong tectonic movement along the boundary of the zones of deformation of the NE corner of the Tibetan plateau (NET) since the new tectonic period, with its dynamic mechanisms remaining controversial. Here, we use observations of 39 Continuous Global Positioning System (CGPS) gauges and 451 Crustal Movement Observation Network of China (CMONOC) campaign-mode stations to detect the three-dimensional deformation of the crust in the NET. Improved processing procedures are implemented to strengthen the patterns of strain throughout the NET. The principal component analysis (PCA) technique is introduced to decompose the time series into spatial eigenvectors and principal components (PCs), and the first three PCs are used to estimate and rectify common mode errors (CMEs). In addition, GRACE observations are used to detect deformation changes that account for non-tidal oceanic mass loading, hydrological loading, and surface pressure. The rectified deformation of the crust indicates the anisotropic nature of both the subsidence and uplift, and that the highest uplift rate of the Longmen Shan fault uplift reaches 7.13 ± 0.53 mm/yr. Finally, the horizontal velocity is further used to enumerate the strain rates throughout the NET. The results show that the shear band retained property in line with the strike-slip fault along the Altyn Tagh fault, the Qilian Shan faults, the Haiyuan fault, the West Qinling fault, the East Kunlun fault, and the Longmen Shan fault. In addition, the results further indicate that the whole NET shows a strong relationship with the mean principal rates of horizontal shortening strain. Extension and compression of the crust reasonably describe its sinking and uplifting.

Geochemistry of Mudstones/Silty Mudstones from the Qigequan Formation and Shizigou Formation in Yuejin-II Area, Southwestern Area of the Qaidam Basin: Implications for Sedimentary Environment and Sandstone-Type Uranium Mineralization

The Qaidam Basin has been the focus of sandstone-type uranium prospecting since the 1950s. In recent years, relying on the uranium geological survey project supported by the China Geological Survey, and the cooperation with the Qinhai oil company, drilling work at the Yuejin-II area in this basin has achieved breakthroughs on industrial-level sandstone-type uranium exploration. In this study, we present major, trace and REE geochemical analysis of the Qigequan Formation and the Shizigou Formation mudstones/silty mudstones collected from an industrial uranium ore drillholes in the Yuejin-II area. The Shizigou and Qigequan Formations exhibit signatures of non-intense alteration, low rock maturity and proximal provenance. The overall arid paleoclimatic conditions controlled the sedimentation of large volume of uranium rich materials. The period of relative humidity prompted the sedimentation of reduced agents. Aided by the neotectonic-induced slopes and tectonic windows, oxygenated fluids migrated along permeable layers and extracted the hexavalent uranium, transported in the form of uranyl ion (UO22+). When the oxygen–uranium-rich fluids finally infiltrated into the reductive sand body, the hexavalent uranium was reduced to tetravalent uranium and deposited in the form of uranium compounds. Large-scale and centrally distributed reductive sand bodies provided favorable ore storage space for the sandstone-type uranium mineralization in the Yuejin-II area.

Geological Significance of Late Permian Magmatic Rocks in the Middle Section of the Ailaoshan Orogenic Belt, SW China: Constraints from Petrology, Geochemistry and Geochronology

The Ailaoshan orogenic belt, located in the SE margin of the Qinghai–Tibet Plateau, is an important Paleo-Tethys suture zone in the eastern margin of the Sanjiang Tethys tectonic domain. The areas of Mojiang and Zhenyuan, located in the middle part of the Ailaoshan orogenic belt, are the key parts of the Ailaoshan Paleo-Tethys Ocean closure and collision orogeny. The rhyolites outcropped in the Mojiang area, and the granite porphyries outcropped in Zhenyuan area, are systematically studied for petrology, isotope geochemistry and geochronology. The Zircon U-Pb geochronology of rhyolites and granite porphyries give weighted average ages of 253.4 ± 4.2 Ma and 253.3 ± 2.0 Ma, respectively, both of which were formed in the late Permian period. The rhyolites belong to potassic calc-alkaline to subalkaline series. The patterns of the rare earth elements (REE) show a right-inclined seagull-type distribution, and the trace elements plot is right-inclined. The granite porphyries are high potassic calc-alkaline to subalkaline. The REE patterns show a right-inclined distribution, and the trace elements plot is right-inclined, which is consistent with the typical patterns observed in the crust. The peraluminous, highly differentiated and high ASI values suggest that rhyolites and granite porphyries are S-type granites. The zircon εHf(t) of the rhyolites range from −7.22 to −0.72, and two-stage Hf zircon model ages are (TDMC) 1771–2352 Ma, indicating that the magma source area is mainly crust-derived. The zircon εHf(t) of the granite porphyries range from −0.97 to 4.08, and two-stage Hf zircon model ages are (TDMC) 1336–1795 Ma, indicating that the magma is derived from a depleted mantle source and the partial melting of ancient crustal materials. The rhyolites and granite porphyries were possibly formed in the syn-collisional tectonic setting during the late Permian, and their ages limited the time of the final closure of the Ailaoshan Ocean and the initiation of collisional orogeny.

Distribution and Characteristics of Damming Landslides Triggered by 1920 M~8 Haiyuan Earthquake, NW China

Earthquake-triggered damming landslides threaten downstream residents and affect the regional landscape by disrupting water and sediment fluxes. Therefore, it is essential to study the distribution characteristics and distinctive controlling factors of earthquake-triggered damming landslides to provide a reference for treating landslide dams caused by damming landslides. This study uses the 1920 M~8 Haiyuan earthquake-triggered landslides as an example to study the characteristics and topographic effects of damming landslides in the Loess Plateau in Northwestern China. A detailed Haiyuan-earthquake-triggered damming landslide inventory was established. The statistics of terrain, geology, seismic factors, and information gain rankings were used to quantify the significance of the controlling factors. The aspect ratio, equivalent coefficient of friction, area, and slope position was calculated. Damming landslides’ distinctive geomorphic and morphological characteristics were summarized through comparisons with non-damming landslides. The results showed that damming landslides were concentrated in areas with thick loess sediment, low relief, and close proximity to a river. Loess thickness was the most critical control factor among them. Damming landslides have the geomorphological characteristics of a large ratio of length to width (L/W), a low ratio of height to length (H/L), large scales, and entire-slope failure. Moreover, damming landslides can transform the topography of the Loess Plateau through their long-term effects. These findings highlight the characteristics of damming landslides in the Loess Plateau and supplement the global landslide dam inventory. They provide a reference for assisting in earthquake-triggered damming landslides treatments in the Loess Plateau.

Asymmetric Interseismic Strain across the Western Altyn Tagh Fault from InSAR

As the northern boundary of the Tibetan Plateau, the long Altyn Tagh fault (ATF) controls the regional tectonic environment, and the study of its long-term fault slip rate is key to understanding the tectonic evolution and deformation of the northern Tibetan Plateau. In this paper, we measure the fault slip rate of the western segment of the ATF using InSAR observations between 2015 to 2020. The Multi-Temporal Interferometric InSAR analysis is applied to obtain the two-dimensional fault-parallel and vertical displacement fields. The spatially dense InSAR observations clearly illustrate the asymmetrical pattern of displacement fields across the fault. Constrained by our InSAR observations, the fault slip rate and locking depth of the western segment of the ATF are inverted using four different models in a Bayesian framework. The two-layer viscoelastic model incorporating lateral heterogeneity of rheology in the lower crust indicates that the fault slip rate of the western ATF is estimated to be 9.8 ± 1.1 mm/yr (at 83.8°E across the ATF) and 8.6 ± 1.1 mm/yr (at 85.1°E), respectively, and the locking depth is 15.8 ± 4.3 km and 14.8 ± 4.9 km. Our new estimates generally agree with the previous estimates of fault slip rate constrained by GPS observations. We conclude that the contrast between the thickness of the elastic layer and the shear modulus of the Tibetan plateau and the Tarim basin jointly contribute to the asymmetric interseismic strain accumulation on the ATF.

Phylogeny and Biogeographic History of Parnassius Butterflies (Papilionidae: Parnassiinae) Reveal Their Origin and Deep Diversification in West China

Simple Summary

Butterflies of the genus Parnassius are distributed in the mountains across the Northern Hemisphere. Studies have shown that this genus originated in the regions of West China–Central Asia. To further explore the spatiotemporal pattern and driving mechanism of Parnassius diversification, we reconstructed the phylogeny and biogeographic history of Parnassius based on 45 species. Ancestral area reconstruction obtained by using the statistical dispersal–extinction cladogenesis model indicated that Parnassius originated in West China (Qinghai–Tibet Plateau and Xinjiang) during the Middle Miocene. Paleoenvironment changes and host plants were probably influenced by the dispersal of Parnassius butterflies from West China to East Asia, Europe, and North America. Furthermore, ancient gene introgression might have contributed to the spread of Parnassius butterflies from the high mountains of the Qinghai–Tibet Plateau to the low-altitude areas of Central East China. This study will provide an understanding of the phylogeny and biogeographic history of the genus Parnassius.

Abstract

We studied 239 imagoes of 12 Parnassius species collected from the mountains of the Qinghai–Tibet Plateau (QTP) and its neighbouring areas in China. We selected three mitochondrial gene (COI, ND1, and ND5) sequences, along with the homologous gene sequences of other Parnassius species from GenBank, to reconstruct the phylogenetic tree and biogeographic history of this genus. Our results show that Parnassius comprises eight monophyletic subgenera, with subgenus Parnassius at the basal position; the genus crown group originated during the Middle Miocene (ca. 16.99 Ma), and species diversification continued during sustained cooling phases after the Middle Miocene Climate Optimum (MMCO) when the QTP and its neighbouring regions experienced rapid uplift and extensive orogeny. A phylogenetic network analysis based on transcriptomes from GenBank suggests that ancient gene introgression might have contributed to the spread of the Parnassius genus to different altitudes. Ancestral area reconstruction indicates that Parnassius most likely originated in West China (QTP and Xinjiang) and then spread to America in two dispersal events as subgenera Driopa and Parnassius, along with their host plants Papaveraceae and Crassulaceae, respectively. Our study suggests that extensive mountain-building processes led to habitat fragmentation in the QTP, leading to the early diversification of Parnassius, and climate cooling after MMCO was the driving mechanism for the dispersal of Parnassius butterflies from West China to East Asia, Europe, and North America.

The Potential of Sentinel-1A Data for Identification of Debris-Covered Alpine Glacier Based on Machine Learning Approach

Microwave remote sensing is one of the main approaches to glacier monitoring. This paper provides a comparative analysis of how different types of radar information differ in identifying debris-covered alpine glaciers using machine learning algorithms. Based on Sentinel-1A data, three data suites were designed: A backscattering coefficient (BC)-based data suite, a polarization decomposition parameter (PDP)-based data suite, and an interference coherence coefficient (ICC)-based data suite. Four glaciers with very different orientations in different climatic zones of the Tibetan Plateau were selected and classified using an integrated machine learning classification approach. The results showed that: (1) The boosted trees and subspace k-nearest neighbor algorithms were optimal and robust; and (2) the PDP suite (63.41–99.57%) and BC suite (55.85–99.94%) both had good recognition accuracy for all glaciers; notably, the PDP suite exhibited better rock and debris recognition accuracy. We also analyzed the influence of the distribution of glacier surface aspect on the classification accuracy and found that the more asymmetric it was about the sensor orbital plane, the more difficult it was for the BC and PDP suites to recognize the glacier, and a large slope could further reduce the accuracy. Our results suggested that during the inventory or classification of large-scale debris-covered alpine glaciers, priority should be given to polarization decomposition features and elevation information, and it is best to divide the glaciers into multiple subregions based on the spatial relationship between glacier surface aspect and radar beams.

Phylogeny and biogeography of the northern temperate genus Dracocephalum s.l. (Lamiaceae)

The northern temperate genus Dracocephalum consists of approximately 70 species mainly distributed in the steppe-desert biomes of Central and West Asia and the alpine region of the Qinghai-Tibetan Plateau (QTP). Previous work has shown that Dracocephalum is not monophyletic and might include Hyssopus and Lallemantia. This study attempts to clarify the phylogenetic relationships, diversification patterns, and the biogeographical history of the three genera (defined as Dracocephalum s.l.). Based on a sampling of 66 taxa comprising more than 80% from extant species of Dracocephalum s.l., morphological, phylogenetic (maximum parsimony, likelihood, and Bayesian inference based on nuclear ITS and ETS, plastid rpl32-trnL, trnL-trnF, ycf1, and ycf1-rps15, and two low-copy nuclear markers AT3G09060 and AT1G09680), molecular dating, diversification, and ancestral range estimation analyses were carried out. Our results demonstrate that both Hyssopus and Lallemantia are embedded within Dracocephalum and nine well-supported clades can be recognized within Dracocephalum s.l. Analyses of divergence times suggest that the genus experienced an early rapid radiation during the middle to late Miocene with major lineages diversifying within a relatively narrow timescale. Ancestral area reconstruction analyses indicate that Dracocephalum s.l. originated in Central and West Asia and southern Siberia, and dispersed from Central and West Asia into the QTP and adjacent areas twice independently during the Pliocene. The aridification of the Asian interior possibly promoted the rapid radiation of Dracocephalum within this region, and the uplift of the QTP appears to have triggered the dispersal and recent rapid diversification of the genus in the QTP and adjacent regions. Combining molecular phylogenetic and morphological evidence, a revised infrageneric classification of Dracocephalum s.l. is proposed, which recognizes nine sections within the genus.

Pulsed rise and growth of the Tibetan Plateau to its northern margin since ca. 30 Ma

The onset of mountain building along margins of the Tibetan Plateau provides a key constraint on the processes by which the high topography in Eurasia formed. Although progressive expansion of thickened crust underpins most models, several studies suggest that the northern extent of the plateau was established early, soon after the collision between India and Eurasia at ca. 50 Ma. This inference relies heavily on the age and provenance of Cenozoic sediments preserved in the Qaidam basin. Here, we present evidence in the northern plateau for a considerably younger inception and evolution of the Qaidam basin, based on magnetostratigraphies combined with detrital apatite fission-track ages that date the basin fills to be from ca. 30 to 4.8 Ma. Detrital zircon-provenance analyses coupled with paleocurrents reveal that two-stage growth of the Qilian Shan in the northeastern margin of the Tibetan Plateau began at ca. 30 and at 10 Ma, respectively. Evidence for ca. 30 and 10 to 15 Ma widespread synchronous deformation throughout the Tibetan Plateau and its margins suggests that these two stages of outward growth may have resulted from the removal of mantle lithosphere beneath different portions of the Tibetan Plateau.

Coseismic Deformation Field Extraction and Fault Slip Inversion of the 2021 Yangbi MW 6.1 Earthquake, Yunnan Province, Based on Time-Series InSAR

An earthquake of moderate magnitude (MW 6.1) occurred in Yangbi County, Dali, Yunnan Province, China, on 21 May 2021. Compared to strong earthquakes, the measurement of the deformation fields of moderate earthquakes is more susceptible to errors associated with atmospheric, orbital, and topographic features. We adopted a new time-series InSAR method to process preseismic and postseismic Sentinel-1A SAR time-series images and separated the coseismic deformation signals from various error signals. This method uses preseismic time-series interferograms to estimate the spatially correlated look angle error induced by the digital elevation model and the atmospheric and orbital errors in the master image. The preseismic and postseismic time-series interferograms were then segmented for spatio-temporal filtering to provide a precise estimate of the atmospheric and orbital errors in slave images. Such time-series processing accurately separates various errors from the coseismic deformation signal and prevents the coseismic deformation signal from being included as noise in the error estimation during filtering. Based on this approach, we effectively eliminated the masking of the deformation signal by the errors and extracted coseismic deformation field of the Yangbi MW 6.1 earthquake with high precision. The maximum LOS displacement in the ascending and descending tracks were determined to be −74 and −62 mm, respectively. Subsequently, we used the Geodetic Bayesian Inversion Software to invert the fault geometric parameters of this earthquake, and based on this inverted the rupture slip distribution using the least-squares method. The results showed that the fault orientation is 133.43°, dip angle is 76.98°, source depth is 5.5 km, fault sliding mode is right-lateral strike-slip. The moment magnitude (MW) was calculated to be 6.07.

The rise and demise of the Paleogene Central Tibetan Valley

Reconstructing the Paleogene topography and climate of central Tibet informs understanding of collisional tectonic mechanisms and their links to climate and biodiversity. Radiometric dates of volcanic/sedimentary rocks and paleotemperatures based on clumped isotopes within ancient soil carbonate nodules from the Lunpola Basin, part of an east-west trending band of basins in central Tibet and now at 4.7 km, suggest that the basin rose from <2.0 km at 50 to 38 million years (Ma) to >4.0 km by 29 Ma. The height change is quantified using the rates at which wet-bulb temperatures (T_w) decline at land surfaces as those surface rise. In this case, T_w fell from ~8°C at ~38 Ma to ~1°C at 29 Ma, suggesting at least ~2.0 km of surface uplift in ~10 Ma under warm Eocene to Oligocene conditions. These results confirm that a Paleogene Central Tibetan Valley transformed to a plateau before the Neogene.