Supraglacial meltwater reaching a glacier bed can increase ice surface velocities via hydraulic jacking and enhanced basal sliding. However, the relationships between the structure of supraglacial drainage systems, sink-point distributions, glacier flow processes and the magnitude of interannual velocity variability are poorly understood. To explore the hypothesis that spatial variations in the rate and mechanisms of glacier flow are linked to variations in supraglacial drainage system structure and sink-point distribution across an ice cap, we mapped supraglacial drainage systems on Devon Ice Cap from Landsat-7 ETM+ imagery. Spatial patterns of surface velocity and interannual velocity variability were determined using gradient correlation applied to Landsat-7 ETM+ images. Velocity variability is greater in areas close to sink-point locations, presumably because hydrologically forced basal sliding and/or bed deformation are enhanced in such areas. The distribution and characteristics of supraglacial drainage systems may play an important role in determining the distribution of regions of basal sliding, highlighting the need for knowledge of the supraglacial drainage system structure and sink-point distribution to inform efforts to model the dynamic response of Arctic ice caps to future climate warming.