Debris-covered glaciers are ubiquitous in the Himalaya, and supraglacial debris significantly alters how glaciers respond to climate forcing. Estimating debris thickness at the glacier scale, however, remains a challenge. This study inverts a subdebris melt model to estimate debris thickness for three glaciers in the Everest region from digital elevation model difference-derived elevation change. Flux divergences are estimated from ice thickness and surface velocity data. Monte Carlo simulations are used to incorporate the uncertainties associated with debris properties, flux divergence, and elevation change. On Ngozumpa Glacier, surface lowering data from 2010 to 2012 and 2012 to 2014 are used to calibrate and validate the method, respectively. The debris thickness estimates are consistent with existing in situ measurements. The method performs well over both actively flowing and stagnant parts of the glacier and is able to accurately estimate thicker debris ({ extgreater}0.5 m). Uncertainties associated with the thermal conductivity and elevation change contribute the most to uncertainties of the debris thickness estimates. The surface lowering associated with ice cliffs and supraglacial ponds was found to significantly reduce debris thickness, especially for thicker debris. The method is also applied to Khumbu and Imja-Lhotse Shar Glaciers to highlight its potential for regional application.