Mantle xenoliths carried by the Cenozoic basalts in Huinan are composed of harzburgite, dunite, lherzolite, wehrlite, and pyroxenite. These rocks provide an opportunity to trace the composition and evolution of the lithospheric mantle beneath the northeastern North China Craton (NCC). Major and trace element analyses reveal that the primary compositions of these xenoliths are widely modified by partial melting and multistage metasomatism. Three types of metasomatic agents are identified based on the geochemical features of clinopyroxene. The high Sr-87/Sr-86 (>0.7065), high Ti/Eu but low Ca/Al ratios of clinopyroxene core in one refractory harzburgite indicate an early metasomatic event caused by subducted Paleo-Asian Oceanic slab-derived high-Sr-87/Sr-86 silicate melts. The variation trend of clinopyroxenes in harzburgite, dunite, and wehrlite xenoliths in the plot of Ti/Eu-(La/Yb)(N) indicate the overprinting of subsequent carbonatite metasomatism. The low Sr-87/Sr-86 ratio (<0.7035) of carbonatite melts suggest an origin of the asthenosphere mantle. In addition, Sr-87/Sr-86 ratios of clinopyroxenes in some pyroxenite and dunite xenoliths increase with increasing U/Zr and Ti/Eu ratios, while in some harzburgites, they decrease with increasing U/Zr and Ti/Eu ratios. These features together suggest a later metasomatic event caused by fluid-bearing silicate melts with moderately high Sr-87/Sr-86 ratio (similar to 0.7048), which were suggested to be related to the subduction of the Paleo-Pacific Oceanic slab. Considering the multistage mantle metasomatism widely recorded in the NCC, we speculate that the early mantle metasomatism could have played an important role in the thinning of the NCC. (C) 2020 Elsevier B.V. All rights reserved.