New zircon LA-ICPMS U–Pb age data, whole-rock major and trace elements, and zircon Hf isotopic data of intrusive rocks from Larelaxin and Caima plutons in the southern margin of the Western Qiangtang subterrane provide important evidence of northward subduction of the Bangong–Nujiang Tethyan Ocean lithosphere. Host granitoids including quartz diorites, granodiorites, and syenogranites and associated mafic varieties including gabbroic enclaves, dioritic dikes and enclaves, and quartz dioritic enclaves are investigated in this study. Five host granitoid samples are dated at 163–160 Ma, and one dioritic dike and three dioritic enclave samples are dated at 162–158 Ma, which indicates that these rocks were contemporaneous. The quartz diorites and granodiorites are normal calc-alkaline I-type granitoids. The syenogranites are characterized by high SiO2 (74–77 wt.%) and differentiation index (DI = 92–97) and marked depletion in Ba, Nb, Ta, Sr, P, Ti, and Eu, thus indicating that they are highly fractionated I-type granitoids. The host granitoids exhibit uniform zircon εHf(t) values (− 1.4 to + 1.9). Considering the current data and their negative whole-rock εNd(t) values (− 4.7 to − 3.5) reported in recent studies, these rocks can be interpreted to have resulted from the partial melting of ancient mafic lower crust with varying contributions from mantle-derived or juvenile crust-derived components. One gabbroic enclave sample exhibits high Al2O3 (18.3 wt.%), low MgO (3.3 wt.%), high TiO2 (1.4 wt.%), and high Nb (24 ppm) and is geochemically similar to high-alumina basalt (HAB) and high-Nb basalt (HNB). This indicates that the gabbroic enclave originated from the partial melting of mantle wedge peridotite that was metasomatized by slab melting and subsequently experienced significant fractional crystallization of olivine and clinopyroxene. The dioritic dike and enclave samples yield εHf(t) values (− 1.6 to + 3.1) similar to those of the host granitoids. They most likely originated from magma mixing between mantle-derived and crust-derived melts at or close to the Moho. The intrusive rocks with geochemical diversity of the Larelaxin and Caima plutons from Bangong Tso in western Tibet can be interpreted to have resulted from melting, assimilation, storage, and homogenization (a MASH process) above subduction zones, thus indicating that the Bangong–Nujiang Tethyan Ocean lithosphere was subducted northward beneath the Western Qiangtang subterrane during the Late Jurassic.