As debris cover affects the meltwater production and behaviour of glaciers it is important to understand how, and over what timescales, it forms. Here we develop an advanced 3-D numerical model that describes transport of sediment through a glacier to the point where it emerges at the surface. The numerical performance of the model is satisfactory and it reproduces debris structures observed within real-world glaciers, thereby offering a useful tool for future studies of debris-covered glaciers.