In this study, we apply a glacier mass balance and ice redistribution model to simulate historical and future glacier change in the Everest region of Nepal. High-resolution temperature and precipitation fields derived from gridded APHRODITE data, and validated against independent station observations from the EVK2CNR network, are used to drive the historical model from 1961 to 2007. The model is calibrated against geodetically derived estimates of net glacier mass change from 1992 to 2008, termini position of four large glaciers at the end of the calibration period, average velocities observed on selected debris-covered glaciers, and total glacierized area. We integrate field-based observations of glacier mass balance and ice thickness with remotely-sensed observations of decadal glacier change to validate the model. Between 1961 and 2007, the mean modelled volume change over the Dudh Kosi basin is −6.4 ± 1.5 km3, a decrease of 15.6% from the original estimated ice volume in 1961. Modelled glacier area change between 1961 and 2007 is −101.0 ± 11.4 km2, a decrease of approximately 20% from the initial extent. Scenarios of future climate change, based on CMIP5 RCP4.5 and RCP8.5 end members, suggest that glaciers in the Everest region will continue to lose mass through the 21st century. Glaciers in the basin are concentrated between 5000 and 6000 m of elevation, and are thus expected to be sensitive to changes in temperature and equilibrium line altitude (ELA). Glacier volume reductions between −35 to −62% are possible by 2050, and sustained temperature increases to 2100 may result in total glacier volume losses of between −73 and −96%.