India-Asia Collision Research Articles

Genesis of a sediment-hosted vein-type Cu deposit in the Lanping Basin, SW China: New insights from geology and LA–ICP–MS analysis of fluid inclusions

Sediment-hosted Cu deposits, including stratiform and epigenetic vein-type Cu deposits, are the second largest source of Cu worldwide. However, many facets of these vein-type Cu deposits, such as the nature and origin of ore fluids and geodynamic mechanisms that produced these fluids, remain controversial, limiting our interpretations of the ore genesis and hampering the selection of exploration strategies. Numerous sediment-hosted vein-type Cu deposits are present in the Lanping Basin, SW China, of which the Jinman deposit is the largest. Copper orebodies at Jinman consist of Cu-sulfide–quartz–carbonate veins and altered wall-rocks hosted mainly in faults and fracture zones in metamorphosed and deformed clastic rocks. The deposit was formed in three stages, pre-, syn- and post-ore, with the first two stages producing hydrothermal veins. The veins in pre-ore stage consist of quartz + ankerite ± sericite ± pyrite, whereas those in syn-ore stage consist of Cu-sulfides + quartz + carbonate ± sericite. During the syn-ore stage, tennantite, quartz, and ankerite crystallized early, then bornite was precipitated, followed by chalcopyrite and chalcocite and finally calcite. Only minor sericite was formed during the post-ore stage. Two major types of fluid inclusions—CO2-rich and aqueous inclusions—are present in both pre- and syn-ore quartz, whereas only aqueous inclusions are present in syn-ore calcite. In pre- and syn-ore quartz and calcite, CO2-rich inclusions have higher homogenization temperatures (275–323 °C) and pressures (307–3495 bar) and lower salinities (0.22–2.7 wt% NaCl eqv.) than aqueous inclusions (164–299 °C, 6.5–80.4 bar and 3.7–20.9 wt% NaCl eqv.). However, the small variations in Cs/Na and Cs/(Na + K) ratios [log Cs/Na = –3.3 ∼ –1.9, log Cs/(Na + K) = –3.3 ∼ –2.1] of the fluid inclusions suggest a same source for the ore fluids. We propose that during both stages, CO2-rich inclusions resulted from unmixing of initial deep-seated, single-phase CO2-rich fluids, whereas aqueous inclusions likely represent the fluids that evolved from these unmixed fluids by extensive CO2 degassing. This interpretation, combined with microthermometric results of CO2-rich inclusions, suggests a mesothermal, low-salinity H2O–NaCl–CO2 ore fluid system. The fluid inclusions from Jinman have similar compositions to those of metamorphic fluids, suggesting a metamorphic origin for the ore-forming fluids. The evolution of fluid compositions suggests that fluid degassing likely acted as a major mechanism for Cu precipitation. We propose that the ore fluids were generated through metamorphic devolatilization of basement of the Lanping Basin in the Biluoxueshan–Chongshan metamorphic zone under a compressional or transpressional regime caused by Cenozoic India–Asia collision. The Jinman deposit has many similarities to orogenic gold deposits worldwide and thus can be thought of as an orogenic deposit. Based on a compilation of the geological characteristics of global sediment-hosted epigenetic vein-type Cu deposits, we suggest that other sediment-hosted epigenetic vein-type Cu deposits with similar genesis to the Jinman deposit are likely present in other regions worldwide. This highlights the exploration potential of this type of deposit in these regions and the need for further research.

Tectonic geomorphology and Quaternary slip history of the Fuyun fault, southwestern Altai Mountains, central Asia

The northwest-trending Altai Mountains of central Asia expose a complex network of thrust and strike-slip faults that are key features accommodating intracontinental crustal shortening related to the Cenozoic India-Asia collision. In this study, we investigated the Quaternary slip history of the Fuyun fault, a right-lateral strike-slip fault bounding the southwestern margin of the Altai Mountains, through geologic mapping, geomorphic surveying, and optically stimulated luminescence (OSL) geochronology. At the Kuoyibagaer site, the Fuyun fault displaces three generations of Pleistocene–Holocene fill-cut river terraces (i.e., T3, T2, and T1) containing landslide and debris-flow deposits. The right-lateral offsets are magnified by erosion of terrace risers, suggesting that river course migration has been faster than slip along the Fuyun fault. The highest Tp2 terrace was abandoned in the middle Pleistocene (150.4 ± 8.1 ka uppermost OSL age) and was displaced 145.5 +45.6/–12.1 m along the Fuyun fault, yielding a slip rate of 1.0 +0.4/–0.1 mm/yr since the middle Pleistocene. The lower Tp1 terrace was abandoned in the late Pleistocene and aggraded by landslides and debris flows in the latest Pleistocene–Holocene (36.7 ± 1.6 ka uppermost OSL age). Tp1 was displaced 67.5 +14.2/–6.1 m along the Fuyun fault, yielding a slip rate of 1.8 +0.5/–0.2 mm/yr since the late Pleistocene. Our preferred minimum slip rate of ~1 mm/yr suggests the Fuyun fault accommodates ~16% of the average geodetic velocity of ~6 mm/yr across the Altai Mountains. Integration of our new Fuyun slip rate with other published fault slip rates accounts for ~4.2 mm/yr of convergence across the Chinese Altai, or ~70% of the geodetic velocity field.

Largest known madtsoiid snake from warm Eocene period of India suggests intercontinental Gondwana dispersal

Here we report the discovery of fossils representing partial vertebral column of a giant madtsoiid snake from an early Middle Eocene (Lutetian, ~ 47 Ma) lignite-bearing succession in Kutch, western India. The estimated body length of ~ 11–15 m makes this new taxon (Vasuki indicus gen et sp. nov.) the largest known madtsoiid snake, which thrived during a warm geological interval with average temperatures estimated at ~ 28 °C. Phylogenetically, Vasuki forms a distinct clade with the Indian Late Cretaceous taxon Madtsoia pisdurensis and the North African Late Eocene Gigantophis garstini. Biogeographic considerations, seen in conjunction with its inter-relationship with other Indian and North African madtsoiids, suggest that Vasuki represents a relic lineage that originated in India. Subsequent India-Asia collision at ~ 50 Ma led to intercontinental dispersal of this lineage from the subcontinent into North Africa through southern Eurasia.

An updated chronostratigraphic framework for the Cenozoic sediments of southeast margin of the Tibetan Plateau: Implications for regional tectonics

As one of the most important accommodation zones during the India-Asia collision, the southeast margin of the Tibetan Plateau (SEMTP) is characterized by many large-scale strike-slip faults, a regional low-relief relict that has been deeply incised by several continental scale rivers that originate from central Tibet, and numerous plant and vertebrate fossil biotas residing in Cenozoic sedimentary basins. Therefore, the SEMTP not only provides a meaningful constraint on the geodynamic evolution of the Tibetan Plateau, but also serves as an exemplar to understand the various interactions among surface uplift, drainage network reorganization, climate change, and biodiversity. Precise dating of Cenozoic sediments in the basins of the SEMTP is crucial for understanding these processes. Previously, the geochronologic framework of Cenozoic sediments in the SEMTP largely relied on fossils and regional stratigraphic correlations, which often involved circular reasoning and large uncertainties. In the past two decades, high-resolution magnetostratigraphy, detrital zircon geochronology, and isotopic dating of interbedded tephra layers have been carried out on SEMTP Cenozoic sedimentary basins in Yunnan and eastern Xizang provinces, SW China. Using these new data, we build a revised chronostratigraphic framework for Cenozoic sediments in SEMTP, and use this updated temporal framework to contextualize regional tectonics and climate change. The new chronostratigraphic framework shows that (1) many of the “Neogene” sedimentary basins based on plant fossils and regional lithostratigraphic correlation in the SEMTP were actually formed in the late Eocene and early Oligocene; (2) most of the Paleogene sedimentary successions in the SEMTP ended at the late Eocene-early Oligocene and were unconformably overlaid by Middle-Late Miocene sediments. Oligocene-Early Miocene sediments are often regionally absent. Sedimentary structures within the basin successions show that Paleogene basins were generally formed in compressional settings while the Late Miocene basins were mostly developed in extensional settings. The new chronostratigraphic framework, together with analyses of sedimentary basin structures, suggests that the SEMTP began to experience crustal shortening early in the India-Asia collision (∼50 Ma), producing Paleogene sedimentary basins and giving rise to the subsequent widespread surface uplift and regional erosion in the Oligocene-Middle Miocene. Tectonic inversion from compression to extension since the Middle Miocene precludes significant uplift of the SEMTP after this time, which may be related to geodynamic changes in the evolution of the Tibetan Plateau.

Crustal Deformation in Eastern Himalayan Syntaxis Constrained by Ambient Noise Tomography

ABSTRACT As the leading edge of the Indian–Asian collision, the eastern Himalayan syntaxis region has experienced extensive tectonic activities, resulting in complex crustal uplift and deformation in the corner area of the southeastern pathway for the extrusion of Tibetan plateau materials. Despite considerable efforts, the corresponding deformation mechanisms remain uncertain. This study presents a new 3D high-resolution azimuthal anisotropic shear-wave velocity model in the crust and uppermost mantle derived from ambient noise dispersion data. Results show that the upper crustal anisotropy aligns with the geological boundaries and major faults nearby, suggesting shape-preferred orientations. The upper crustal low velocity and weak anisotropy beneath the core of the eastern Himalayan syntaxis (EHS) are closely associated with the high fragmentation of shallow rocks and the upwelling of hot materials during the ongoing subduction of the Indian plate. Our model also reveals relatively complex anisotropic patterns in the midlower crust. The eastern Lhasa terrane, in particular, exhibits low velocity and strong anisotropy with a northwest–southeast-oriented fast axis, supporting the local scale midlower crustal “channel flow” model. In addition, a conspicuous, elongated low-velocity zone along the northwest–southeast direction is observed in the midlower crust and uppermost mantle beneath the Bangong–Nujiang suture. The anisotropy in this region increases with depth, and the fast directions are consistently parallel to the northeast subduction of the Indian plate. We infer that this low-velocity zone may result from partial melting under local compression driven by the Indian–Asian collision. On the basis of newly revealed anisotropic model and previous studies, we construct a new dynamic model, which reveals that the migration of mechanically weak material in the midlower crust and the significant contribution of the northeast subduction of the Indian plate jointly control the crustal deformation of the EHS region.

Segmentation of the Tashkurgan normal fault in the eastern Pamir: Insights from geomorphology and thermochronology and implications for fault-slip transfer

At the northwestern end of the India–Asia collision zone, in the eastern Pamir interior, the Kongur Shan extensional system extends for ∼250 km as a composite system of normal faults. As the southernmost segment of the extensional system, the Tashkurgan fault can be divided into a northern and southern segment by the intersection of the Tahman and Tashkurgan faults. To evaluate the tectonic activity along the Tashkurgan fault and explore the extensional characteristics of the southern Kongur Shan extensional system, we integrate geomorphologic analysis with thermochronologic data and modeling. Geomorphic indices, including local relief, normalized channel steepness index (ksn), and river χ-plots, show that the southern Tashkurgan fault is more active than the northern Tashkurgan fault. On the footwall of the northern Tashkurgan fault, most thermochronologic data, which include biotite 40Ar/39Ar (Arbt), zircon (U-Th)/He (ZHe), apatite fission track (AFT), and apatite (U-Th-Sm)/He (AHe) dates, record emplacement-related cooling, and only a few AHe dates possibly record exhumation of <2 km since ∼7 Ma ago. On the footwall of the southern Tashkurgan fault, almost all of the ZHe, AFT, and AHe dates record significant late Miocene to present exhumation; thus, we perform three-dimensional thermokinematic modeling for the normal slip along this fault segment. Modeling results indicate a probably constant slip rate on the southern Tashkurgan fault at a value of 1.4–1.5 mm/a, which corresponds to 6–7 km of both horizontal extension and footwall uplift since 6.5 Ma ago. All together, the geomorphic and thermochronologic evidences imply that the southern Tashkurgan fault is responsible for the ∼E–W extension along the southernmost portion of the Kongur Shan extensional system, while the northern fault segment plays a minor role in the extension. This result, combined with the previously published magnitude of extension along the southern Kongur Shan fault, indicates that the extension at any site along the southern Kongur Shan extensional system can be primarily attributed to normal faulting along a single border fault, which is a typical characteristic of an interbasin transfer zone with a transfer fault.

Petrogenesis and Tectonic Implications of the Oligocene Dalongtan Shoshonitic Syenite Porphyry in Central Yunnan, Southeastern Tibetan Plateau: Constraints from Geochronology, Geochemistry and Sr-Nd-Hf Isotopes

Shoshonitic rocks are widely distributed in post-collisional settings and provide key information on deep geodynamic mechanisms and magmatic evolution. In this paper, we present petrographic, zircon U-Pb age-related, trace elemental, Hf isotopic, bulk-rock elemental, and Sr-Nd isotopic data of the Dalongtan shoshonitic syenite porphyries (DSSPs) in central Yunnan, southeastern Tibet. The DSSPs formed at 33.2 ± 0.3 Ma in a post-collisional setting. They define linear trends on Harker diagrams, and they display similar trace element patterns and enriched bulk-rock Sr-Nd isotopes [(87Sr/86Sr)i = 0.70964–0.70968, εNd(t) = −12.9 to −12.7] and zircon Hf isotopes (εHf(t) = −15.7 to −13.1) to the coeval mantle-derived potassic mafic rocks. This suggests that the DSSPs were fractionated from the lithospheric mantle-derived mafic magmas. The DSSPs, along with the coeval felsic and mafic magmatic rocks (37.2–32.3 Ma), exhibit a planar distribution on the SE Tibet and predate the left-lateral shearing of the Ailaoshan–Red River shear zone (ARSZ) (32–22 Ma), suggesting that there are no genetic relationships between them. The DSSPs have geochemical characteristics similar to those of A-type granites, with high total alkalinity (10.39–11.17 wt.%), HFSE concentrations (Zr + Nb + Ce + Y = 890.2–1054.3 ppm), Ga/Al ratios (10,000 × Ga/Al = 2.95–3.46), whole-rock zircon saturation temperatures (906–947 °C), and oxygen fugacity (ΔFMQ = +3.30–+4.65), indicating that they are products of the high-temperature melting of the lithosphere as a result of asthenosphere upwelling in extensional settings. Based on our data and regional observations, it is proposed that the generation of the DSSPs may be linked to the convective thinning of the thickened lithospheric mantle following the India–Asia collision.

Structure and provenance of the Cretaceous Pingshanhu Basin in the Hexi Corridor: Implications for Mesozoic tectonics in the northern Tibetan Plateau

Abstract The construction of Earth’s largest highland, the Tibetan Plateau, is generally considered to have been generated by the Cenozoic India-Asia collision. However, the extent to which high topography existed prior to the Cenozoic remains unclear. The Hexi Corridor foreland basin of the northern Tibetan Plateau is an ideal region in which to investigate this history, given its widespread exposure of Early Cretaceous sedimentary sequences. In this study, we examined the Early Cretaceous strata in the northern Hexi Corridor to understand the relationships between pre-Cenozoic sedimentation and tectonic deformation and constrain the late Mesozoic tectonic setting of the adjacent Qilian Shan and Alxa blocks bordering the northern Tibetan Plateau. Results of sandstone petrology analyses, paleocurrent observations, and U-Pb geochronology suggest that the oldest Early Cretaceous sediments deposited in the northern Hexi Corridor were sourced from the southern Alxa block during the earliest Cretaceous. By the late Early Cretaceous, Hexi Corridor sediments were sourced from both the southern Alxa block to the north and the Qilian Shan to the south. Sandstone petrologic results indicate that the northern Hexi Corridor experienced a tectonic transition from contraction to extension during the Early Cretaceous. These findings suggest that the northern Tibetan Plateau region was partially uplifted to a high elevation during the late Mesozoic before the India-Asia collision.

Early Cretaceous magmatism archives crustal architecture of northern Lhasa Terrane, central Tibetan Plateau

Identification of continental crustal architecture in orogenic belts is essential for understanding the geodynamic evolution of global continental collision systems. However, ongoing debates surrounding the crustal nature of the Lhasa Terrane in the central Tibetan Plateau impede our comprehension of its tectonic evolution. Mesozoic intermediate-felsic magmatism in the Lhasa Terrane offers key insights into the crustal architecture before the India-Asia collision. Here we report zircon U-Pb-Hf isotope, whole-rock element and SrNd isotope data from Early Cretaceous granitoids in the Baingoin batholith, northern Lhasa Terrane. These granitoids (127–115 Ma) exhibit variable SiO2 (62.99–73.86 wt%), K2O (1.15–4.74 wt%), and MgO (0.49–2.48 wt%) contents. They are classified as low- to high-K calc-alkaline and metaluminous to peraluminous series, enriched in large ion lithophile elements (e.g., Rb) and light rare earth elements ((La/Yb)N = 7.08–18.3), and depleted in high field strength elements (e.g., Nb, Ta, and Ti). These granitoids display geochemical affinities with I-type granitoid rocks and can be further subdivided into three groups. Group 1 granitoids exhibit the highest zircon εHf(t) (−1.57 to 5.85) and whole-rock εNd(t) (−3.37 to −2.45), and lowest (87Sr/86Sr)i (0.706018 to 0.706922) values, whereas Group 2 granitoids have wide ranges of zircon Hf (εHf(t) = −7.78 to 2.26) and enriched whole-rock SrNd isotope compositions ((87Sr/86Sr)i = 0.707214 to 0.707918; εNd(t) = −7.38 to −6.14). These suggest that Group 1 and Group 2 granitoids may be dominantly derived from the juvenile and ancient crustal sources, respectively. Group 3 granitoids are characterized by zircon Hf (εHf(t) = −2.85 to 3.52) and whole-rock SrNd isotope compositions ((87Sr/86Sr)i = 0.706473 to 0.707692; εNd(t) = −4.50 to −3.01) falling between the ranges of the previous two groups, likely represent a mixture of juvenile crustal-derived melt with 10–15% of reworked ancient crustal materials. Integrating previous research, we propose that the persistent ancient crustal basement beneath the northern Lhasa Terrane was concurrently undergoing significant modification through juvenile crustal materials during the Early Cretaceous. These Early Cretaceous Baingoin batholith preserves crucial records of crustal evolution and highlights the influence of upwelling asthenosphere during the transition from Bangong-Nujiang oceanic subduction to Lhasa-Qiangtang collision.

Protracted and Progressive Crustal Melting during Continental Collision in the Pamir and Plateau Growth

Abstract Determining crustal melting in parallel with geodynamic evolution provides critical information on plateau crustal thickening and uplift. Here we investigate the timing and duration of crustal melting through in situ analysis of zircon U–Pb ages, trace elements and Hf–O isotopes, and whole-rock elements and Sr–Nd–Hf isotopes for the granites and high-grade metamorphic rocks from the Pamir Plateau. Zircon dates record protracted crustal melting for both Central Pamir (43–33 and 22–12 Ma) and South Pamir (28–10 Ma). The Pamir Cenozoic granites are characterized by significant elemental and isotopic heterogeneity. The elemental variability within the Pamir Cenozoic granites is dominantly attributed to fractional crystallization of K-feldspar and plagioclase with subordinate biotite from a variably fractionated melt, and the accumulation of early crystallized feldspar during magma ascent. Peritectic mineral entrainment and accessory mineral crystallization had some influence on the geochemical characteristics of the garnet-bearing leucogranite dikes. Zircon Hf isotopes and whole-rock Sr–Nd–Hf isotopes show secular variations for both Central and South Pamir granites. The Central Pamir granites show a mild decrease in whole-rock ƐNd(t) values from Eocene (−4.3 to −4.9) to Miocene (−6.2 to −7.7), and the zircon ƐHf(t) values decrease from c. 40 Ma (+2 to −5) to c. 10 Ma (−4 to −8). In contrast, the South Pamir granites have highly variable whole-rock Sr–Nd–Hf (87Sr/86Sr(i) = 0.7053 to 0.7830; ƐNd(t) = −31.5 to +0.2; ƐHf(t) = −40.0 to +8.2) and zircon Hf isotopes (ƐHf(t) = −31.7 to +7.5) and display a strong decrease in ƐNd(t) and ƐHf(t) values from c. 25 Ma to c. 13 Ma. Geochemical and isotopic data indicate that both the Central and South Pamir experienced crustal melting from juvenile lower crust to ancient lower-middle crustal materials, and Indian crustal materials were incorporated into the melt region of the South Pamir leucogranites from c. 20 Ma. Our study highlights a causal link between a chain of events that includes magma underplating induced by lithosphere thinning and slab breakoff, lithosphere delamination and underthrusting of Indian lithosphere, and formation of the Cenozoic granites in Pamir. This series of processes are incorporated here into a comprehensive model for the geodynamic evolution of the Pamir during the India–Asia collision.

Linking Porphyry Cu Formation to Tectonic Change in Postsubduction Settings: A Case Study from the Giant Yulong Belt, Eastern Tibet

Abstract Porphyry deposits in magmatic arcs form coincident with changes to steady-state oceanic subduction conditions, such as changes in plate convergence rate and vector or angle of subduction. However, it remains unclear whether such processes also operated during formation of postsubduction porphyry deposits. The Yulong magmatic belt in the eastern Tibetan Plateau consists of middle to late Eocene igneous rocks (~51–35 Ma) that formed during the India-Asia collision, whereas all known porphyry deposits are associated with late Eocene rocks (43–35 Ma). A synthesis of new and published geochemical data shows marked variations from the middle to late Eocene, including increasing whole-rock La/Yb, Sr/Y, and EuN/EuN* values, as well as zircon EuN/EuN* values. These geochemical variations, together with petrographic observations, indicate a transition from plagioclase-dominated to amphibole-dominated fractionation from the middle to late Eocene. Coupled changes of magma compositions and porphyry Cu metallogeny from the middle to late Eocene coincided with, or were slightly preceded by, the onset of regional uplift and crustal thickening, triggered by the India-Asia hard collision and rapid deceleration of the India-Asia convergence rate at ca. 50 to 44 Ma. Crustal thickening may have caused prolonged magma differentiation at greater depths and accumulation of dissolved H2O, both of which contributed to amphibole-dominated fractionation and generation of hydrous melt that are prospective for porphyry Cu mineralization. Our study highlights the importance of tectonic changes in the formation of the Yulong and other postsubduction porphyry Cu belts—a scenario similar to that operated in subduction-related settings such as the Andes.

Eocene crustal thickening in the Tethyan Himalaya: Insights from Barrovian metamorphism and granite geochemistry from the Ramba area

Magmatism, structures, and metamorphism in the Ramba dome of the Tethyan Himalaya were investigated to shed light on orogenic processes during the early stages of the India-Asia collision. Deformed granite dikes in the dome envelope yield zircon U-Pb ages of ca. 45 Ma. These Eocene granites have adakitic, Na-rich compositions (K2O/Na2O = 0.20−0.61), weak to no Eu anomaly, enrichment in Sr, depletion in heavy rare earth elements and Y, and low MgO and Mg# contents. These characteristics contrast with the Miocene potassic granites in the core of the dome and suggest that the Eocene adakites were derived from the high-pressure melting of crustal amphibolites in a thick crust. The mica schists of the dome envelope have an early foliation (S1) that is overprinted by upright folds (F2). Phase-equilibria modeling of garnet and staurolite mica schists suggests a Barrovian-type, prograde P-T evolution in association with S1, with peak conditions of 6.7−7.2 kbar/590−605 °C and 7.3−7.8 kbar/650−670 °C, respectively, which are typical of crustal thickening metamorphism. Monazites from S1-dominated staurolite mica schists yield metamorphic ages of ca. 51−49 Ma, while those from the late foliation (S2) that transposed S1 give younger ages of ca. 10 Ma. The integration of geochemical, structural, metamorphic, and geochronological data suggests that peak Barrovian D1 metamorphism and adakitic magmatism occurred in the Eocene in response to crustal thickening. The results provide critical constraints for addressing the crustal shortening deficit of the region.

Polyphase Exhumation of the East Kunlun Orogenic Belt: Evidence from Modern River Detrital Zircon and Apatite Fission Track Dating

Abstract The East Kunlun Orogenic Belt (E-KOB) stands out as one of the most prominent basin-mountain geomorphic features in the northern interior of the Tibetan Plateau. It records a series of accretion-collision events from the Mesozoic to the Cenozoic. In particular, with the uplifting of the Tibetan Plateau, the E-KOB experienced intracontinental deformation and exhumation in the Cenozoic. Clarifying the exhumation history of the E-KOB is crucial to define the growth time and mechanism of the Tibetan Plateau. In this study, we apply detrital zircon fission-track (ZFT) and apatite fission-track (AFT) analyses on modern river sands in order to constrain the regional exhumation history of the eastern E-KOB. Four peak ages have been identified and interpreted as results of rapid exhumation correlated with intracontinental deformation. Two older peak ages at 144.7–141.0 and 114.6–82.1 Ma are in good accordance with the collision time of the north-south Lhasa-Qiangtang Block along the Bangong-Nujiang suture zone and the subsequent progressive deformation stage toward the north. Peak age at 60.9–45.3 Ma is coeval with the initial timing of the India-Asia collision. The youngest peak age at 25.1–18.3 Ma matches well with the extensive outward and upward growth of the Tibetan Plateau during the Oligocene to Miocene time. The Cretaceous and early Cenozoic rapid exhumations suggest that the E-KOB has been involved in the intracontinental deformation induced by collisions of the Lhasa-Qiangtang and India-Asia from the south. It implies that the northern Tibetan Plateau likely has been elevated or was a structural high before the Eocene. In addition, some of our detrital samples show a younger ZFT peak age than the AFT peak age. We attributed this data bias to the contribution of hydrodynamic sorting and/or lithological difference. The combination of ZFT and AFT dating has advantages in eliminating interfering age signals in detrital thermochronology.

Timing of rifting of the Dongkaco microcontinent (Central Tibet) and implications for Neo-Tethyan evolution

Rifting of continental blocks away from northern Gondwana resulted in the development of several Tethyan oceanic branches, the number of which and the timing of tectonic processes are actively debated. Within the Mesozoic Bangong-Nujiang suture, locally exposed Paleozoic strata are potential witnesses of the early development of Tethyan basins, yet their origin and tectonic significance remain poorly understood. In this study, Paleozoic limestone and sandstone units exposed in the Dongkaco area are identified as part of an independent Dongkaco microcontinent, differentiated from the South Qiangtang terrane by the Dongqiao suture on the north and from the Lhasa Block by the Beila suture on the south. The Upper Paleozoic deltaic to shallow-marine feldspatho-quartzose and quartzose sandstones within the microcontinent are interbedded with glaciogenic sequences and exhibit provenance characteristics that align with the Lhasa Block. These stratigraphic and provenance features suggest that during the Carboniferous to Early Permian the Dongkaco microcontinent was attached to the Lhasa Block. The discovery of Middle Triassic (Anisian) radiolarian chert associated with the Beila ophiolite indicates that sea-floor spreading was active within the Beila Ocean at that time. Data presented here show that the detachment of the Dongkaco microcontinent from the Lhasa Block resulted in the Beila Ocean, a branch of Bangong-Nujiang Ocean that expanded through Middle Permian to Early Triassic times. These findings shed new light on the paleogeographic evolution that preceded India-Asia collision and final amalgamation of the Himalayan-Tibetan orogen.

Timing of India-Asia suturing: Evidence from a remnant peripheral foreland basin in Xigaze, South Tibet

The absence of a fully preserved foreland basin system in the central-eastern Yarlung-Zangpo Suture Zone (YZSZ) of Tibet has hindered the understanding of the initial timing of the India-Asia collision along strike. To address this issue, remnant foredeep and wedge-top deposits of the foreland basin system are described from the Tethyan Himalaya sequence and the YZSZ in the Xigaze area. The foredeep is defined by the olistostrome of the Zongzhuo Formation, which consists of matrix and blocks, including chert, limestone, sandstone and volcanic rocks. Based on the matrix properties, the Zongzhuo Formation is divided into two units: a lower member (Zhilong unit) with a matrix composed of black shale and thin layers of siltstone and an upper member (Yongla unit) with a matrix composed of distinctly red shale or yellow-green siliceous mudstone with interbedded siltstone and limestone. The Zongzhuo Formation is dated to 65–57 Ma by detrital zircon U - Pb chronology and radiolarian assemblages. In the YZSZ, the Puxia and Jiabi units are interpreted as wedge-top deposits consisting of mudstone matrix and exotic blocks, such as sandstone, chert and limestone. The maximum depositional ages (MDAs) of the sandstone blocks from both units are ca. 65–59 Ma. Detrital zircon provenance analysis and trace element data exclude an intraoceanic arc origin and point to a source from the Gangdese arc and recycled accretionary wedge in the suture zone. Therefore, we propose that the initial collision between India and Asia occurred at 65 Ma in the central-eastern YZSZ.

Late Oligocene to Early Pliocene Exhumation and Structural Development in the Western Himalaya, Northern Pakistan: Implications for the Cenozoic Metamorphic Overprint

Abstract New middle Miocene to Pliocene (~14–3 Ma) apatite fission track (AFT) cooling ages combined with published K–Ar/Ar–Ar and zircon fission track (ZFT) ages from the Hazara and Swat regions of Pakistan are used to explain the Oligocene to Pliocene structural evolution in the Western Himalaya. The structural model explains the distribution of K–Ar/Ar–Ar ages in three distinct age groups (Proterozoic, Paleozoic-Mesozoic, and Eocene to Oligocene). The Proterozoic to Mesozoic sequence of northern Hazara and Swat experienced elevated temperature and pressure conditions, evident by reset Eocene to Oligocene K–Ar/Ar–Ar hornblende and Eocene to Miocene muscovite ages, caused by Kohistan overthrusting the Indian margin during and after the India–Asia collision. Samples from the Indus syntaxis with Paleo to Mesoproterozoic K–Ar/Ar–Ar hornblende ages and Eocene to Oligocene Ar–Ar muscovite ages show no signs of Cenozoic metamorphism; these samples were thermally imprinted up to the Ar–Ar muscovite closure temperature. Neoproterozoic to Lower Paleozoic rocks from the southern parts of Hazara and Swat show Mesozoic to Oligocene partially reset Ar–Ar muscovite ages and preservation of Ordovician metamorphism. The combined analysis of published K–Ar/Ar–Ar (muscovite), ZFT, and new AFT ages (~14–12 Ma) suggests that the Main Central thrust/Panjal thrust was active from Oligocene to early Miocene (~30–18 Ma), and the Nathia-Gali and Main Boundary thrusts were active from the middle to late Miocene (~14–9 Ma) in the Hazara area. New and published AFT ages (~6–3 Ma) from the Indus syntaxis suggest that early Pliocene tectonic thickening in the hinterland formed the N–S trending Indus anticline, creating an erosional half window in the Main Mantle thrust, forming the Indus syntaxis, and dividing the Main Central thrust sheet into the Hazara and Swat segments.

Successive accretions of future allochthonous terranes and multiple subduction zone jumps: Implications for Tethyan evolution

Abstract The accretion of future allochthonous terranes (e.g., microcontinents or oceanic plateaus) onto the southern margin of Asia occurred repeatedly during the evolution and closure of the Tethyan oceanic realm, but the specific geodynamic processes of this protracted convergence, successive accretion, and subduction zone initiation remain largely unknown. Here, we use numerical models to better understand the dynamics that govern multiple terrane accretions and the polarity of new subduction zone initiation. Our results show that the sediments surrounding the future terranes and the structural complexity of the overriding plate are important factors that affect accretion of multiple plates and guide subduction polarity. Wide (≥400 km) and buoyant terranes with sediments behind them and fast continental plate motions are favorable for multiple unidirectional subduction zone jumps, which are also referred to as subduction zone transference, and successive terrane-accretion events. The jumping times (~3–20 + m.y.) are mainly determined by the convergence rates and rheology of the overriding complex plate with preceding terrane collisions, which increase with slower convergence rates and/or a greater number of preceding terrane collisions. Our work provides new insights into the key geodynamic conditions governing multiple subduction zone jumps induced by successive accretion and discusses Tethyan evolution at a macro level. More than 50 m.y. after India-Asia collision, subduction has yet to initiate along the southern Indian plate, which may be the joint result of slower plate convergence and partitioned deformation across southern Asia.

Multi‐channel seismic reflection study of tectonic–sedimentary features and subduction initiation in the middle Kyushu–Palau Ridge and adjacent basins

AbstractThe Kyushu–Palau Ridge and adjacent basins are ideal locations for investigating the formation and evolution of marginal seas and initiation of plate subduction. In this study, the tectonic–sedimentary features and crustal structure of the Kyushu–Palau Ridge and adjacent basins were investigated using newly obtained deep seismic reflection and borehole data. The initial mechanism of subduction in the West Philippine Sea and its tectonic evolution are discussed. Two sets of sedimentary strata with different provenances occur in the eastern Philippine Basin. The thickness of the lower strata is variable, and most of the sediment was sourced from island arc volcanism on the Kyushu–Palau Ridge. These strata thicken towards the Kyushu–Palau Ridge, and volcaniclastic rock aprons are developed near the foot of the Kyushu–Palau Ridge. The upper strata have a relatively uniform thickness and comprise fine‐grained, deep‐water marine sediments. The crustal thickness of the West Philippine and Parece Vela Basins is 6–7 km, which is similar to the global average thickness of oceanic crust. The Moho beneath the West Philippine Basin has a broad and undulating form that mimics the thickness of the oceanic crust beneath the sediments. The depth to the Moho decreases towards the Central Basin Spreading Center. The Moho beneath the Parece Vela Basin is planar, which is in contrast to the undulating nature of the oceanic crust base in this area. The West Philippine Basin may have been located at the northern margin of Australia in the Southern Hemisphere during the Mesozoic. It developed gradually in a continental margin arc as a result of Palaeogene inter‐arc and oceanic extension and contains fragments of continental crust. Seismic profiles and drillhole data in the West Philippine Basin reveal that compression occurred during the Eocene. Subduction along the paleo‐Izu–Bonin–Mariana arc may have been caused by the far‐field effects of the India–Asia collision. Subduction was accompanied by lateral propagation and a compressive stress field until the break‐up of the island arc at ca. 30 Ma.

Comparison of Sn-related granitoids in subduction and collision settings by accessory mineral geochemistry: A case study in the Tengchong-Lianghe tin belt, SW China

Magmatic-hydrothermal Sn deposits can develop in both subduction and collision settings. However, a systematic comparison of the geochemical fingerprints and petrogenetic evolution of Sn-related granitoids from both settings is still lacking. The Tengchong-Lianghe tin belt hosts numerous Sn deposits but of different metallogenic ages, such as the Xiaolonghe and Lailishan Sn deposits, which provide an ideal study area for addressing this issue. Zircon U–Pb dating suggests that the Xiaolonghe and Lailishan Sn-related granitoids were emplaced at 74.2–75.5 and 50.0–51.1 Ma, respectively. Given the initial India-Asia collision commenced at ca. 65–60 Ma, these ages, consistent with available Sn metallogenic ages, suggest the formation of the Xiaolonghe and Lailishan Sn deposits in subduction and collisional settings, respectively. New zircon and apatite compositions, and apatite Sr isotopic data, combined with previously reported whole-rock geochemical results, all indicate that the Xiaolonghe and Lailishan Sn-related granitoids share similar magma sources, evolution processes, redox states, and volatile contents. The presence of hornblende, high apatite NdN/NdN* (mostly > 1), and whole-rock evolutionary trends, together indicate that they represent fractionated I-type granites. Distinctly elevated apatite 87Sr/86Sr ratios (0.70988–0.71430) suggest that they originated from an ancient lower crust. Subsequently, they underwent an extensive fractional crystallization dominated by feldspar, as revealed by whole-rock, apatite, and zircon elemental variation trends. The low zircon CeN/CeN* (mainly < 200) and ΔFMQ (–1.9 to 1.7), together with elevated apatite F content (2.30–3.69 wt%), reveal that they crystallized from reduced and F-rich magmas. The relatively high whole-rock zircon saturation temperatures (mainly > 800 ℃) and Ba/Pb ratios imply that they were produced at a higher temperature by biotite-dehydration melting, which requires additional heat from the mantle. The specific mechanism that triggers mantle upwelling (oceanic slab subduction or break-off) could be the most significant difference between Sn-related granitoids formed during subduction and collision processes.

Active CO2 emissions from thermal springs in the Karakoram fault system and adjacent regions, western Tibetan Plateau

Large-scale intra-continental strike-slip faults provide the pathways for groundwater circulation and the release of CO2-rich fluids through thermal springs, which is a significant type of deep CO2 emissions into the atmosphere and hydrosphere. India-Asia collision created numerous strike-slip faults and thermal springs in the Tibetan Plateau and adjacent regions, but their CO2 origins and fluxes are still not well understood. Here, we present a quantitative study on CO2 emissions from thermal springs in the Karakoram fault system (KKFS) and adjacent regions, western Tibetan Plateau, with the aim of revealing the origins of CO2-rich fluids and estimating the dissolved and gaseous CO2 fluxes. Hydrogen and oxygen isotopic data show that the geothermal waters were recharged by local precipitation (e.g., rainfall and snow-melt water) and experienced variable degrees of water-rock interaction, giving rise to the dominant Na–HCO3 type and minor Na–Cl and Ca–HCO3 type waters discharging from thermal springs. Spatially, the thermal springs within or close to the KKFS (Group 1 and Group 2) have larger circulation depths of 1.3–5.1 km, while those away from the KKFS (Group 3) exhibit shallow (0.6−0.8 km) groundwater circulation. Mass balance calculations reveal high deeply-sourced carbon inputs (60.4−94.3%) to dissolved inorganic carbon (DIC) of the Group 1 and Group 2 samples, while the Group 3 samples have higher shallow carbon contributions (50.4−60.9%) to DIC via dissolution of carbonate and organic carbon. This is consistent with higher Δ14CDIC values (−788.6‰ to −692.1‰) in Group 3 samples and the lack of radioactive 14C in Group 1 and Group 2 samples (Δ14CDIC = −999.9‰ to −985.9‰). Based on chemical compositions and fluxes of spring water and bubbling gases, the total CO2 output from thermal springs (including dissolved and gaseous CO2) in the KKFS is estimated to be ∼2332 t a−1, in which the endogenic carbon such as metamorphic CO2 accounts for ∼1863 t a−1. Our results provide the first estimate of deep CO2 flux for the KKFS and would contribute to a better understanding of active CO2 emissions facilitated by strike-slip faulting processes in continental collision zone.