Mega‐scale glacial lineations (MSGLs) are a characteristic landform on ice stream beds. Solving the puzzle of their formation is key to understanding how ice interacts with its bed and how this, in turn, influences the dynamics of ice streams. However, a comprehensive and detailed characterization of this landform's size, shape and spatial arrangement, which might serve to test and refine formational theories, is largely lacking. This paper presents a detailed morphometric analysis and comparison of 4043 MSGLs from eight palaeo‐ice stream settings: three offshore (Norway and Antarctica), four onshore (Canada), and one from under a modern ice stream in West Antarctica. The length of MSGLs is lower than previously suggested (mode 1000–2000 m; median 2892 m), and they initiate and terminate at various locations on an ice stream bed. Their spatial arrangement reveals a pattern that is characterized by an exceptional parallel conformity (80% of all mapped MSGLs have an azimuth within 5° from the mean values), and a fairly constant lateral spacing (mode 200–300 m; median 330 m), which we interpret as an indication that MSGLs are a spatially self‐organized phenomenon. Results show that size, shape and spatial arrangement of MSGLs are consistent both within and also generally between different ice stream beds. We suggest this results from a common mechanism of formation, which is largely insensitive to local factors. Although the elongation of MSGLs (mode 6–8; median 12.2) is typically higher than features described as drumlins, these values and those of their width (mode 100–200 m; median 268 m) overlap, which suggests the two landforms are part of a morphological continuum and may share a similar origin. We compare their morphometry with explicit predictions made by the groove‐ploughing and rilling instability theories of MSGL formation. Although the latter was most compatible, neither is fully supported by observations. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.

Atmospheric warming is increasing surface melting across the Antarctic Peninsula, with unknown impacts upon glacier dynamics at the ice-bed interface. Using high-resolution satellite-derived ice velocity data, optical satellite imagery and regional climate modelling, we show that drainage of surface meltwater to the bed of outlet glaciers on the Antarctic Peninsula occurs and triggers rapid ice flow accelerations (up to 100% greater than the annual mean). This provides a mechanism for this sector of the Antarctic Ice Sheet to respond rapidly to atmospheric warming. We infer that delivery of water to the bed transiently increases basal water pressure, enhancing basal motion, but efficient evacuation subsequently reduces water pressure causing ice deceleration. Currently, melt events are sporadic, so efficient subglacial drainage cannot be maintained, resulting in multiple short-lived (<6 day) ice flow perturbations. Future increases in meltwater could induce a shift to a glacier dynamic regime characterised by seasonal-scale hydrologically-driven ice flow variations.

Understanding the evolution of the ice-bed interface is fundamentally important for gaining insight into the dynamics of ice masses and how subglacial landforms are created. However, the formation of the suite of landforms generated at this boundary — subglacial bedforms — is a contentious issue that is yet to be fully resolved. Bedforms formed in aeolian, fluvial, and marine environments either belong to separate morphological populations or are thought to represent a continuum of forms generated by the same governing processes. For subglacial bedforms, a size and shape continuum has been hypothesised, yet it has not been fully tested. Here we analyse the largest data set of subglacial bedform size and shape measurements ever collated (96,900 bedforms). Our results show that flutes form a distinct population of narrow bedforms. However, no clear distinction was found between drumlins and megascale glacial lineations (MSGLs), which form a continuum of subglacial lineations. A continuum of subglacial ribs also exists, with no clear size or shape distinctions indicating separate populations. Furthermore, an underreported class of bedform with no clear orientation to ice flow (quasi-circular bedforms) overlaps with the ribbed and lineation continua and typically occurs in spatial transition zones between the two, potentially merging these three bedform types into a larger continuum

Ice streams drain large portions of ice sheets and play a fundamental role in governing their response to atmospheric and oceanic forcing, with implications for sea-level change. The mechanisms that generate ice stream flow remain elusive. Basal sliding and/or bed deformation have been hypothesized, but ice stream beds are largely inaccessible. Here we present a comprehensive, multi-scale study of the internal structure of mega-scale glacial lineations (MSGLs) formed at the bed of a palaeo ice stream. Analyses were undertaken at macro- and microscales, using multiple techniques including X-ray tomography, thin sections and ground penetrating radar (GPR) acquisitions. Results reveal homogeneity in stratigraphy, kinematics, granulometry and petrography. The consistency of the physical and geological properties demonstrates a continuously accreting, shallow-deforming, bed and invariant basal conditions. This implies that ice stream basal motion on soft sediment beds during MSGL formation is accommodated by plastic deformation, facilitated by continuous sediment supply and an inefficient drainage system.

The Greenland Ice Sheet is a large contributor to global sea level rise, and current mass losses are projected to accelerate. However, model projections of future ice sheet evolution are limited by the fact that the ice sheet is not in equilibrium with present-day climate but is still adjusting to past changes that occurred over thousands of years. While the influence of such committed adjustments on future ice sheet evolution remains unquantified, it could be addressed by calibrating numerical ice sheet models over larger timescales and, importantly, against empirical data on ice margin positions. To enable such paleo data–model interactions, we need Greenland-wide empirical reconstructions of past ice sheet extent that combine geomorphological and geochronological evidence. Despite an increasing number of field studies producing new chronologies, such a reconstruction is currently lacking in Greenland. Furthermore, a time slice reconstruction can help to (i) answer open questions regarding the rate and pattern of ice margin evolution in Greenland since the glacial maximum, (ii) develop a standardised record of empirical data, and (iii) identify new sites for future field campaigns. Based on these motivations, we here present PaleoGrIS 1.0, a new Greenland-wide isochrone reconstruction of ice sheet extent evolution through the Late Glacial and early- to mid-Holocene informed by both geomorphological and geochronological markers. Our isochrones have a temporal resolution of 500 years and span ∼ 7.5 kyr from approximately 14 to 6.5 kyr BP. We describe the resulting reconstruction of the shrinking ice sheet and conduct a series of ice-sheet-wide and regional analyses to quantify retreat rates, areal extent change, and their variability across space and time. During the Late Glacial and early- to mid-Holocene, we find the Greenland Ice Sheet has lost about one-third of its areal extent (0.89 million km²). Between ∼ 14 and ∼ 8.5 kyr BP, it experienced a near-constant rate of areal extent loss of 170 ± 27 km² yr⁻¹. We find that the ice-sheet-scale pattern of margin retreat is well correlated to atmospheric and oceanic temperature variations, which implies a high sensitivity of the ice sheet to deglacial warming. However, during the Holocene, we observe inertia in the ice sheet system that likely caused a centennial- to millennial-scale time lag in ice extent response. At the regional scale, we observe highly heterogeneous deglacial responses in ice extent evident in both the magnitude and rate of retreat. We hypothesise that non-climatic factors, such as the asymmetrical nature of continental shelves and onshore bed topographies, play important roles in determining the regional- to valley-scale dynamics. PaleoGrIS 1.0 is an open-access database designed to be used by both the empirical and numerical modelling communities. It should prove a useful basis for improved future versions of the reconstruction when new geomorphological and geochronological data become available.