Basal control of supraglacial meltwater catchments on the Greenland Ice Sheet
Abstract. Ice surface topography controls the primary routing of surface meltwater on ablation zones of glaciers and ice sheets. Meltwater routing is important for understanding and predicting ice sheet evolution because surface melt can be both a direct source of ice mass loss and an influence on basal sliding and ice advection. Although controls on ice sheet topography at long wavelengths are well known, smaller scale features relevant for meltwater routing are not well understood. Here we examine the effects of two processes that can influence ice sheet surface topography: bed topography transfer and thermal-fluvial incision by supraglacial streams. We implement 2D bed topography and basal sliding transfer functions in seven study regions of the western Greenland Ice Sheet (GIS) ablation zone to study the influence of basal conditions on ice surface topography. Although bed elevation data quality is spatially variable, we find that ∼ 1–10 km scale ice surface features under variable ice thickness, velocity, and surface slope are well predicted by these transfer functions. We then use flow-routing algorithms to extract supraglacial stream networks from 2–5 m resolution digital elevation models, and compare these with synthetic flow networks calculated on ice surfaces predicted by bed topography transfer. Quantitative comparison of these networks reveals that bed topography can explain ∼ 1–10 km surface meltwater routing patterns without significant contributions from thermal-fluvial erosion by streams. We predict how supraglacial internally drained catchment (IDC) patterns on the GIS would change under time-varying ice flow and/or basal sliding regimes. Basal sliding variations exert a significant influence on IDC spatial distribution, and suggest a potential positive feedback between subglacial hydrologic regime to surface IDC patterning. Increased basal sliding will increase IDC spatial density (by decreasing IDC sizes) and cause more disperse meltwater input to the englacial and subglacial environment. This could result in less efficient subglacial channelization and increased basal sliding that would then further increase IDC density.