Abstract. Debris-covered glaciers are common in rapidly eroding alpine landscapes. When thicker than a few centimeters, surface debris suppresses melt rates. If continuous debris cover is present, ablation rates can be significantly reduced leading to increases in glacier length. In order to quantify feedbacks in the debris-glacier-climate system, we developed a 2-D long-valley numerical glacier model that includes englacial and supraglacial debris advection. We ran 120 simulations on a linear bed profile in which a hypothetical steady state debris-free glacier responds to a step increase of surface debris deposition. Simulated glaciers advance to steady states in which ice accumulation equals ice ablation, and debris input equals debris loss from the glacier terminus. Our model and parameter selections can produce 2-fold increases in glacier length. Debris flux onto the glacier and the relationship between debris thickness and melt rate strongly control glacier length. Debris deposited near the equilibriumline altitude, where ice discharge is high, results in the greatest glacier extension when other debris-related variables are held constant. Debris deposited near the equilibrium-line altitude re-emerges high in the ablation zone and therefore impacts melt rate over a greater fraction of the glacier surface. Continuous debris cover reduces ice discharge gradients, ice thickness gradients, and velocity gradients relative to initial debris-free glaciers. Debris-forced glacier extension decreases the ratio of accumulation zone to total glacier area (AAR). Our simulations reproduce the "general trends" between debris cover, AARs, and glacier surface velocity patterns from modern debris-covered glaciers. We provide a quantitative, theoretical foundation to interpret the effect of debris cover on the moraine record, and to assess the effects of climate change on debris-covered glaciers.