Temporal and Spatial Variations of Enriched Source Components in Linzizong Volcanic Succession, Tibet, and Implications for the India–Asia Collision

Abstract The temporal and spatial distribution of enriched source components at sites of continent–continent collision provides critical insights into mantle dynamic processes associated with subduction- and collision-related events. However, determining the origin of such enriched components remains a significant challenge. We report a comprehensive dataset of the Linzizong volcanic succession (LVS) from four locations with varying distance across-strike to the Indus–Yarlung suture in southern Tibet, which marks the exposed surface expression of the India–Asia collision zone. The LVS rocks in this study can be divided into two groups: a calc-alkaline Group 1 (69–55 Ma), mainly including basaltic–andesitic varieties, and a shoshonitic Group 2 (52–50 Ma), consisting predominantly of silicic rocks with minor mafic compositions. Group 1 samples are likely derived from the fractional crystallization of primitive basaltic melts as a result of the partial melting of a metasomatized mantle wedge. These samples are decoupled in Nd–Hf isotopic compositions, suggesting an incorporation of subducting sediment-derived melts into the mantle wedge. The influence of sediment-derived melt on the mantle source increases away from the suture zone toward Asia (i.e., from the south to the north) as indicated by the more enriched Sr, Nd, Pb, and Hf isotopic compositions, as well as elevated Th/La and La/Sm ratios. The heavy δ26Mg values, and high Ba/Th and Sr/Th ratios of samples close to the suture coincide with the dehydration of the subducting Neo-Tethyan slab. Group 2 mafic samples have depleted and coupled εNd–εHf isotopic compositions and are characterized by elevated (La/Yb)N and Dy/Yb ratios as well as low Zr/Nb ratios, indicating an origin of enriched garnet-bearing lithospheric mantle with contributions from asthenosphere-derived materials. Group 2 silicic samples are isotopically enriched both near and far away from the suture, which can be attributed to the involvement of ancient lower crust-derived melt from Tethyan Himalaya and central Lhasa subterrane basement, respectively. Our results show that the LVS are the magmatic response to late subduction (Group 1), and to initial India–Asia collision and slab breakoff (Group 2). Negative trends in the whole-rock Nd and zircon Hf isotopic compositions at ~51 Ma should be treated with caution for geodynamic interpretations, depending on the distance between the samples and the India–Asia suture.