Temperature rise (∆T ) associated with shear-banding of metallic glasses is of great importance for their performance. However, experimental measurement of ∆T is difficult due to temporal and spatial localization of shear bands and, as a result, our understanding of the mechanism of ∆T is limited. Here, based on molecular dynamics simulations we observe a spectrum of ∆T , which depends on both sample size and strain rate, in the shear bands of CuZr metallic glass under tension. More importantly, we find that the maximum sliding velocity of the shear bands correlates linearly with the corresponding ∆T , ranging from ∼25 K up to near the melting point for the samples studied. Taking heat diffusion into account, we expect ∆T to be lower than 25 K for the lower end of sliding velocity. At high temperature, shear band bifurcation and/or multiplication can occur as a negative feedback mechanism that prevents temperature rising well above the melting point.