Evolution of the melt source during protracted crustal anatexis: An example from the Bhutan Himalaya

The chemical compositions of magmatic zircon growth zones provide powerful insight into evolving magma compositions due to their ability to record both time and the local chemical environment. In situ U-Pb and Hf isotope analyses of zircon rims from Oligocene-Miocene leucogranites of the Bhutan Himalaya reveal, for the first time, an evolution in melt composition between 32 and 12 Ma. The data indicate a uniform melt source from 32 Ma to 17 Ma, and the progressive addition of an older source component to the melt from at least ca. 17 Ma. Age-corrected epsilon(Hf) ratios decrease from between -10 and -15 down to values as low as -23 by 12 Ma. Complementary whole-rock Nd isotope data corroborate the Hf data, with a progressive decrease in epsilon(Nd(t)) from ca. 18 to 12 Ma. Published zircon and whole-rock Nd data from different lithotectonic units in the Himalaya suggest a chemical distinction between the younger Greater Himalayan Series (GHS) and the older Lesser Himalayan Series (LHS). The time-dependent isotopic evolution shown in the leucogranites demonstrates a progressive increase in melt contribution from older lithologies, suggestive of increasing LHS involvement in Himalayan melting over time. The time-resolved data are consistent with LHS material being progressively accreted to the base of the GHS from ca. 17 Ma, facilitated by deformation along the Main Central thrust. From 17 Ma, decompression, which had triggered anatexis in the GHS since the Paleogene, enabled melting in older sources from the accreted LHS, now forming the lowermost hanging wall of the thrust.