“…The smooth down-glacier increase in debris thickness, and the corresponding decline of the surface ablation rate as discussed above, provide only a first- order description of the debris effects (Benn and Lehmkuhl, 2000; Scherler and others, 2011b; Banerjee and Shankar, 2013). The role of several other complicating factors, e.g., the presence of numerous thermokarst ephemeral ponds and cliffs that increase local melt rate (Reynolds, 2000; Sakai and others, 2000; Miles and others, 2017), vertical and horizontal variations of the thermal properties of debris (Nicholson and Benn, 2013; Rowan and others, 2017), the random short-scale spatial variation of debris thickness (Mihalcea and others, 2006; Zhang and others, 2011; Nicholson and Mertes, 2017; Rounce and others, 2018) and the accumulation contribution from avalanches (Laha and others, 2017) need to be quantified for accurate surface mass-balance estimates on any typical debris-covered Himalayan glacier. The standard glaciological mass-balance measurement protocol (Kaser and others, 2003) may not be designed to handle the above issues.…”