To understand the formation conditions of debris-covered glaciers, we examined the dimension and shape of debris-covered areas and potential debris-supply (PDS) slopes of 213 glaciers in the Bhutan Himalaya. This was undertaken using satellite images with 2.5 m spatial resolution for manual delineation of debris-covered areas and PDS slopes. The most significant correlation exists between surface area of southwest-facing PDS slopes and debris-covered area. This result suggests that the southwest-facing PDS slopes supply the largest quantity of debris mantle. The shape of debris-covered areas is also an important variable, quantitatively defined using a geometric index. Elongate or stripe-like debris-covered areas on north-flowing glaciers are common throughout the Bhutan Himalaya. In contrast, south-flowing glaciers have large ablation zones, entirely covered by debris. Our findings suggest that this difference is caused by effective diurnal freeze–thaw cycles rather than seasonal freeze–thaw cycles, permafrost degradation, or snow avalanches. In terms of geographic setting, local topography also contributes to glacier debris supply and the proportion of debris cover on the studied glaciers is suppressed by the arid Tibetan climate, whereas the north-to-south asymmetric topography of the Bhutan Himalaya has less influence on the proportion of debris cover.