Abstract. Mountain glaciers provide us a window into past climate change and landscape evolution, but the pattern of glacier evolution at centennial or suborbital timescale remains elusive, especially in monsoonal Himalayas. We simulated the glacier evolution in Bhutanese Himalaya, a typical monsoon influenced region, during the Little Ice Age (LIA), using the Open Global Glacier Model and six paleo-climate datasets. Compared with the mapped glacial landforms, the model can well capture the glacier length changes, especially for the experiment driving by the GISS climate dataset, but overestimates the changes in glacier area. Simulation results reveal four glacial substages at 1270s–1400s, 1470s–1520s, 1700s–1710s, and 1820s–1900s in the study area. From further analysis, a negative correlation between the number of the substages and glacier length was found, which suggests that the number and occurrence of glacial substages are regulated by the heterogeneous responses of glaciers to climate change. In addition, the changes in summer temperature dominated the glacier evolution in this region during the LIA.