Search

Subglacial meltwater drainage can enhance localized melting along grounding zones and beneath the ice shelves of marine-terminating glaciers. Efforts to constrain the evolution of subglacial hydrology and the resulting influence on ice stability in space and on decadal to millennial timescales are lacking. Here, we apply sedimentological, geochemical, and statistical methods to analyze sediment cores recovered offshore Thwaites Glacier, West Antarctica to reconstruct meltwater drainage activity through the pre-satellite era. We find evidence for a long-lived subglacial hydrologic system beneath Thwaites Glacier and indications that meltwater plumes are the primary mechanism of sedimentation seaward of the glacier today. Detailed core stratigraphy revealed through computed tomography scanning captures variability in drainage styles and suggests greater magnitudes of sediment-laden meltwater have been delivered to the ocean in recent centuries compared to the past several thousand years. Fundamental similarities between meltwater plume deposits offshore Thwaites Glacier and those described in association with other Antarctic glacial systems imply widespread and similar subglacial hydrologic processes that occur independently of subglacial geology. In the context of Holocene changes to the Thwaites Glacier margin, it is likely that subglacial drainage enhanced submarine melt along the grounding zone and amplified ice-shelf melt driven by oceanic processes, consistent with observations of other West Antarctic glaciers today. This study highlights the necessity of accounting for the influence of subglacial hydrology on grounding-zone and ice-shelf melt in projections of future behavior of the Thwaites Glacier ice margin and marine-based glaciers around the Antarctic continent., Frontiers in Earth Science, 10, ISSN:2296-6463

The mid-Cretaceous was one of the warmest intervals of the past 140 million years (Myr) driven by atmospheric CO2 levels around 1000 ppmv. In the near absence of proximal geological records from south of the Antarctic Circle, it remains disputed whether polar ice could exist under such environmental conditions. Here we present results from a unique sedimentary sequence recovered from the West Antarctic shelf. This by far southernmost Cretaceous record contains an intact ~3 m-long network of in-situ fossil roots. The roots are embedded in a mudstone matrix bearing diverse pollen and spores, indicative of a temperate lowland rainforest environment at a palaeolatitude of ~82°S during the Turonian–Santonian (93–83 Myr). A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric CO2 contents of 1120–1680 ppmv and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo in high-CO2 climate worlds.

Presence/absence data for vent taxa compiled from published literature for Kemp Caldera and 15 well-studied vent fields in neighbouring oceanic regions.

Summary of ROV Isis deployments in the Kemp Caldera.

We will present new multibeam bathymetry data that make the Anvers-Hugo Trough west of the Antarctic Peninsula one of the most completely surveyed palaeo-ice stream pathways in Antarctica. We interpret landforms revealed by these data as indicating that subglacial water availability played an important role in facilitating ice stream flow in the trough during late Quaternary glacial periods. Specifically, we observe a set of northward-shoaling valleys that are eroded into the upstream edge of a sedimentary basin, extend northwards from a zone containing landforms typical of erosion by subglacial water flow, and coincide spatially with the onset of mega-scale glacial lineations. Water was likely supplied to the ice stream bed episodically as a result of outbursts from a subglacial lake previously hypothesized to have been located in the Palmer Deep basin on the inner continental shelf. In a palaeo-ice stream confluence area, close juxtaposition of mega-scale glacial lineations with landforms that are characteristic of slow, dry-based ice flow, suggests that water availability was also an important control on the lateral extent of these palaeo-ice streams. These interpretations are consistent with the hypothesis that subglacial lakes or areas of elevated geothermal heat flux play a critical role in the onset of many large ice streams. The interpretations also have implications for the dynamic behaviour of the Anvers-Hugo Trough palaeo-ice stream and, potentially, of several other Antarctic palaeo-ice streams. Keywords: multibeam bathymetry, ice stream, subglacial water, landform

The majority of glaciers draining the Antarctic Peninsula Ice Sheet are thinning and retreating rapidly1. It is widely understood that these changes are driven by both a warming ocean and atmosphere. However, there are other mechanisms, including pinning points created by bathymetric highs and a reverse bed gradient, that are thought to have an important control on ice stream behaviour (Weertman, 1974; Jamieson et al., 2012). Our understanding of the interplay between these mechanisms and time-scales over which they are important is currently limited in time to the advent of satellite monitoring. By reconstructing the cause and style of ice stream retreat following the Last Glacial Maximum (LGM; 25-19 ka BP), it is possible to gain a greater insight into the mechanisms which drive glacier retreat (Ó Cofaigh et al., 2014). Sedimentary sequences deposited during the LGM and the subsequent deglaciation on polar continental shelves, provide an important archive of past changes (Ó Cofaigh et al., 2014). Previous studies have typically identified three sediment facies assemblages; sub-glacial, transitional and open marine (Ó Cofaigh et al., 2014; Domack et al., 1988; Smith et al., 2011). Transitional sediment facies are deposited at the grounding line and are often targeted for radiocarbon dating, as they represent the onset of glaciomarine sedimentation following the retreat of grounded ice (Domack et al., 1988; Smith et al., 2014; Heroy et al., 1996). Despite the development of depositional models to help explain the processes occurring at grounding lines (Powell et al., 1995 and 1996), there is still significant uncertainty about the temporal and spatial variations in grounding line sedimentation along and across a palaeo-ice stream trough. Here we use a multi-proxy approach (water content, shear strength, magnetic susceptibility, density, contents of biogenic opal, Total Organic Carbon and CaCO3, grain size distribution and X-radiographs) on marine sediment cores recovered from the Anvers-Hugo Palaeo-Ice Stream Trough (AHT), western Antarctic Peninsula shelf, to identify variability in transitional sediment facies deposited along and across the trough. We discuss possible controls on the variability in transitional sediment facies and how this is related to the rate and style of ice stream retreat. Our data reveal systematic variability in the types and volume of transitional sediments deposited during the last deglaciation of AHT. A detailed analysis of the transitional sediment facies shows that this variability reflects different phases of ice stream behaviour. Large volumes of ice proximal sediment facies recovered seawards of grounding zone wedges are indicative of episodes of grounding line still-stands. Re-advances of the grounding line, concurrent with a shallowing of the reverse bed gradient and a narrowing of the trough, appear to have occurred during the final stages of deglaciation. This is indicated by interlaminated ice-proximal and ice-distal sediment facies within inner shelf cores. Transitional sediment variability additionally captures the evolution of the ice stream during deglaciation, including the formation of a small ice shelf on the inner shelf. Keywords: Antarctic Peninsula, Last Glacial Maximum, ice stream, sediment cores References Cook, A. J., Holland, P. R., Meredith, M. P., Murray, T., Luckman, A. & Vaughan, D. G, 2016. Ocean forcing of glacier retreat in the western Antarctic Peninsula. Science, 353, 283-286. Weertman, J, 1974. Stability of the Junction of an Ice Sheet and an Ice Shelf. Journal of Glaciology, 13, 3-11. Jamieson, S. S. R., Vieli, A., Livingstone, S. J., Cofaigh, C. O., Stokes, C., Hillenbrand, C.-D. & Dowdeswell, J. A, 2012. Ice-stream stability on a reverse bed slope. Nature Geoscience, 5, 799-802. Ó Cofaigh, C., Davies, B. J., Livingstone, S. J., Smith, J. A., Johnson, J. S., Hocking, E. P., Hodgson, D. A., Anderson, J. B., Bentley, M. J., Canals, M., Domack, E., Dowdeswell, J. A., Evans, J., Glasser, N. F., Hillenbrand, C.-D., Larter, R. D., Roberts, S. J. & Simms, A. R, 2014. Reconstruction of ice-sheet changes in the Antarctic Peninsula since the Last Glacial Maximum. Quaternary Science Reviews, 100, 87-110. Domack, E. W. & Harris, P. T, 1998. A new depositional model for ice shelves, based upon sediment cores from the Ross Sea and the Mac. Robertson shelf, Antarctica. Annals of Glaciology, 27, 281-284. Smith, J. A., Hillenbrand, C.-D., Kuhn, G., Larter, R. D., Graham, A. G. C., Ehrmann, W., Moreton, S. G. & Forwick, M, 2011. Deglacial history of the West Antarctic Ice Sheet in the western Amundsen Sea Embayment. Quaternary Science Reviews, 30, 488-505. Smith, J. A., Hillenbrand, C.-D., Kuhn, G., Klages, J. P., Graham, A. G. C., Larter, R. D., Ehrmann, W., Moreton, S. G., Wiers, S. & Frederichs, T, 2014. New constraints on the timing of West Antarctic Ice Sheet retreat in the eastern Amundsen Sea since the Last Glacial Maximum. Global and Planetary Change, 122, 224-237. Heroy, D. C. & Anderson, J. B, 1996. Radiocarbon constraints on Antarctic Peninsula Ice Sheet retreat following the Last Glacial Maximum (LGM). Quaternary Science Reviews, 26, 3286-3297. Powell, R. D., Dawber, M., McInnes, J. N. & Pyne, A. R, 1996. Observations of the Grounding-line Area at a Floating Glacier Terminus. Annals of Glaciology, 22, 217-223. 1Powell, R. D. & Domack, E, 1995. Modern Glacimarine Environments. In: Glacial Environments, Volume 1 (ed. J Menzies). Butterworth-Heinemann, 445-486.

During expedition PS104 with RV Polarstern in February and March 2017 the MARUM MeBo 70 seabed drilling system was deployed at nine sites on the continental shelf of the Amundsen Sea Embayment, West Antarctica. A total of 57 meters of sediment core were recovered from 11 boreholes located in Pine Island Bay, Pine Island Trough, Bear Ridge and Cosgrove-Abbot Trough with recovery rates ranging from 7 to 76%. The main scientific objective of the drilling was to reconstruct the Late Mesozoic to Quaternary environmental history in this part of the Antarctic continental margin, with a special focus on the past dynamics of the marine based West Antarctic Ice Sheet (WAIS) from its inception to the last glacial cycle. Another main goal of the expedition was to test the suitability of the MeBo drill system for operating on the Antarctic continental shelf and recovering pre-glacial and glacially influenced sedimentary sequences. Here we will present the first results of sedimentological investigations carried out on the drill cores. These comprise (i) visual lithological descriptions, (ii) CT-scanning records of core stratigraphy, sedimentary structures, and possible artefacts induced by the drilling process, (iii) measurements of physical properties performed with a multi-sensor core logger, and (iv) characterisation of the geochemical composition of the drilled sedimentary strata using X-ray fluorescence (XRF) scanner data. Preliminary biostratigraphic investigations conducted on board ship indicated that the recovered sedimentary strata were deposited during various time slices spanning from the Late Cretaceous–Palaeocene to the Late Quaternary. We will provide an update of these initial chronological findings. Keywords: Drill cores, shelf sediments, West Antarctic Ice Sheet.

Grounding-zone wedges (GZW) have been mapped on the sea floor in various sectors of the formerly glaciated continental shelf around Antarctica. In most cases, these wedges record periods of grounding-line stillstands during ice-sheet retreat following the Last Glacial Maximum (~26-19 ka BP). The presence of GZWs along the axis of a palaeo-ice stream trough therefore indicates episodic retreat of the grounding line from its LGM to modern position. However, information about their internal structure is sparse, and precise chronological constraints for both the onset and the duration of the stillstands they represent are still lacking. Consequently, the role of GZW formation in modulating post-LGM ice-sheet retreat cannot be reliably quantified. This information is vital, however, for calculating reliable retreat rates during the past, which are essential for evaluating and understanding the significance of modern retreat rates, particularly for the rapidly changing Amundsen Sea sector. Here we present a novel combination of swath bathymetric, reflection seismic, and sub-bottom sediment profiler data from a newly discovered stacked GZW in the Cosgrove-Abbot palaeo-ice stream trough in the eastern Amundsen Sea Embayment. In total, six generations of overlapping GZWs were mapped over a distance of ~40 km. We will present first estimates of GZW volumes through integration of the different geophysical datasets. Additionally, we recovered eight sediment cores, sampling most of the individual GZWs within the stack, which may allow us to establish age constraints for each grounding-line retreat episode. Together with the estimated GZW volumes, the ages from sediment cores may also enable the calculation of sediment flux rates at grounding lines, which remain elusive for Antarctic grounding lines. This knowledge will help refine available post-LGM retreat chronologies for the Amundsen Sea Embayment, which, in turn, serve as a basis for validating and improving ice-sheet models in an area where precise simulations of future retreat are urgently needed.

This paper compiles and reviews marine and terrestrial data constraining the dimensions and configuration of the Antarctic Peninsula Ice Sheet (APIS) from the Last Glacial Maximum (LGM) through deglaciation to the present day. These data are used to reconstruct grounding-line retreat in 5 ka time-steps from 25 ka BP to present. Glacial landforms and subglacial tills on the eastern and western Antarctic Peninsula (AP) shelf indicate that the APIS was grounded to the outer shelf/shelf edge at the LGM and contained a series of fast-flowing ice streams that drained along cross-shelf bathymetric troughs. The ice sheet was grounded at the shelf edge until ∼20 cal ka BP. Chronological control on retreat is provided by radiocarbon dates on glacimarine sediments from the shelf troughs and on lacustrine and terrestrial organic remains, as well as cosmogenic nuclide dates on erratics and ice moulded bedrock. Retreat in the east was underway by about 18 cal ka BP. The earliest dates on recession in the west are from Bransfield Basin where recession was underway by 17.5 cal ka BP. Ice streams were active during deglaciation at least until the ice sheet had pulled back to the mid-shelf. The timing of initial retreat decreased progressively southwards along the western AP shelf; the large ice stream in Marguerite Trough may have remained grounded at the shelf edge until about 14 cal ka BP, although terrestrial cosmogenic nuclide ages indicate that thinning had commenced by 18 ka BP. Between 15 and 10 cal ka BP the APIS underwent significant recession along the western AP margin, although retreat between individual troughs was asynchronous. Ice in Marguerite Trough may have still been grounded on the mid-shelf at 10 cal ka BP. In the Larsen-A region the transition from grounded to floating ice was established by 10.7–10.6 cal ka BP. The APIS had retreated towards its present configuration in the western AP by the mid-Holocene but on the eastern peninsula may have approached its present configuration several thousand years earlier, by the start of the Holocene. Mid to late-Holocene retreat was diachronous with stillstands, re-advances and changes in ice-shelf configuration being recorded in most places. Subglacial topography exerted a major control on grounding-line retreat with grounding-zone wedges, and thus by inference slow-downs or stillstands in the retreat of the grounding line, occurring in some cases on reverse bed slopes.

Multibeam swath bathymetry datasets collected over the past two decades have been compiled to identify palaeoice stream pathways in the easternmost Amundsen Sea Embayment. We mapped 3010 glacial landforms to reconstruct palaeo-ice flow in the �250 km-long Abbot Glacial Trough that was occupied by a large palaeo-ice stream, fed by two tributaries (Cosgrove and Abbot) that reached the continental shelf edge during the last maximum ice-sheet advance. The mapping has enabled a clear differentiation between glacial landforms interpreted as indicative of wet- (e.g. mega-scale glacial lineations) and cold-based ice (e.g. hill-hole pairs) during the last glaciation of the continental shelf. Both the regions of fast palaeo-ice flow within the palaeo-ice stream troughs, and the regions of slow palaeo-ice flow on adjacent seafloor highs (referred to as inter-ice stream ridges) additionally record glacial landforms such as grounding-zone wedges and recessional moraines that indicate grounding line stillstands of the ice sheet during the last deglaciation from the shelf. As the palaeo-ice stream flowed along a trough with variable geometry and variable subglacial substrate, it appears that trough sections characterized by constrictions and outcropping hard substrate that changes the bed gradient, led the pace of grounding-line retreat to slow and subsequently pause, resulting in the deposition of grounding-zone wedges. The stepped retreat recorded within the Abbot Glacial Trough corresponds well to post-glacial stepped retreat interpreted for the neighbouring Pine Island-Thwaites Palaeo-Ice Stream trough, thus suggesting a uniform pattern of episodic retreat across the eastern Amundsen Sea Embayment. The correlation of episodic retreat features with geological boundaries further emphasises the significance of subglacial geology in steering ice stream flow. Our new geomorphological map of the easternmost Amundsen Sea Embayment resolves the pathways of palaeo-ice streams that were probably all active during the last maximum extent of the ice sheet on this part of the shelf, and reveals the style of postglacial grounding-line retreat. Both are important input variables in ice sheet models and therefore can be used for validating the reliability of these models.

The Weddell Sea sector is one of the main formation sites for Antarctic Bottom Water and an outlet for about one fifth of Antarctica’s continental ice volume. Over the last few decades, studies on glacialegeological records in this sector have provided conflicting reconstructions of changes in ice-sheet extent and ice-sheet thickness since the Last Glacial Maximum (LGM at ca 23e19 calibrated kiloyears before present, cal ka BP). Terrestrial geomorphological records and exposure ages obtained from rocks in the hinterland of the Weddell Sea, ice-sheet thickness constraints from ice cores and some radiocarbon dates on offshore sediments were interpreted to indicate no significant ice thickening and locally restricted grounding-line advance at the LGM. Other marine geological and geophysical studies concluded that subglacial bedforms mapped on theWeddell Sea continental shelf, subglacial deposits and sediments over-compacted by overriding ice recovered in cores, and the few available radiocarbon ages from marine sediments are consistent with major ice-sheet advance at the LGM. Reflecting the geological interpretations, different icesheet models have reconstructed conflicting LGM ice-sheet configurations for the Weddell Sea sector. Consequently, the estimated contributions of ice-sheet build-up in the Weddell Sea sector to the LGM sealevel low-stand of w130 m vary considerably. In this paper, we summarise and review the geological records of past ice-sheet margins and past icesheet elevations in the Weddell Sea sector. We compile marine and terrestrial chronological data constraining former ice-sheet size, thereby highlighting different levels of certainty, and present two alternative scenarios of the LGM ice-sheet configuration, including time-slice reconstructions for post- LGM grounding-line retreat. Moreover, we discuss consistencies and possible reasons for inconsistencies between the various reconstructions and propose objectives for future research. The aim of our study is to provide two alternative interpretations of glacialegeological datasets on Antarctic Ice- Sheet History for the Weddell Sea sector, which can be utilised to test and improve numerical icesheet models

High-resolution swath bathymetry data collected during several research cruises over the past two decades reveal a palaeo-ice stream trough (Abbot Glacial Trough) crossing the middle and outer shelf of the easternmost Amundsen Sea Embayment, east of the main Pine Island Trough. Regions of both fast palaeo-ice flow (within the central trough) and slow palaeo-ice flow (on adjacent seafloor highs referred to as inter-ice stream ridges) bear glacial landforms indicative of phases of grounding-line stabilization of the ice sheet. We associate a grounding-zone wedge situated within the outer Abbot Glacial Trough with a grounding-zone wedge in outer Pine Island Trough and suggest a synchronous grounding-line halt in both troughs. New sediment echosounder and sediment core data collected from outer Abbot Glacial Trough, between the seaward limit of the grounding-zone wedge and the shelf edge, reveal an up to 6 m-thick well stratified drape that is composed of unconsolidated glaciomarine sediments occasionally bearing calcareous microfossils. In order to decipher whether this unusually thick sediment drape might indicate sub-ice shelf and/or seasonal-open marine deposition throughout or since the Last Glacial Maximum, we used a multi-proxy approach to characterize its lithofacies and applied radiocarbon dating of calcareous microfossils. Here we present our initial results and discuss since when the outer shelf in the eastern Amundsen Sea has been free of grounded-ice. Such information will 1) improve ice sheet models that aim to reconstruct the flow and extent of the West Antarctic Ice Sheet during the Last Glacial Maximum, 2) help to quantify the ice volume of the West Antarctic Ice Sheet during this time, and 3) prove or reject the possibility that Antarctic benthic biota endured glacial periods in outer shelf refugia.

Ice loss from the marine-based, inherently unstable West Antarctic Ice Sheet (WAIS) contributes to the currently observed rise in sea-level and may raise it by up to 3.3-5 metres in the future. Over the last few decades, glaciers draining the WAIS into the Amundsen Sea Embayment (ASE), in particular into Pine Island Bay, have shown thinning, grounding-line retreat and ice-flow acceleration at dramatic rates. These changes are mainly attributed to significant ice-shelf melting by upwelling warm deep water. A critical unknown, limiting our ability to accurately predict future WAIS behaviour, is the poorly constrained long-term context of ice-sheet retreat in the ASE. Here we present a new pre-Holocene to present chronology for WAIS retreat in Pine Island Bay (PIB) based on radiocarbon dating of marine sediment cores. The dates give evidence that grounded ice had retreated close to its modern-day position by ~10 ka BP. Maximum average retreat rates calculated from the deglaciation ages suggest, that the current rapid WAIS retreat in Pine Island Bay is unprecedented over the last ~10 ka and originates in recent changes in regional climate, ocean circulation or ice-sheet dynamics. However, our data and previously published ages for grounding-line retreat from the wider ASE further demonstrate, that, other than in the Ross Sea, the WAIS did not retreat continuously since the LGM. A unique assemblage of glacial morphological features mapped on the eastern ASE shelf suggest a more complex deglacial history, with ice masses slowly flowing and/or stagnating on topographic highs (’Inter-ice stream ridges’) adjacent to main palaeo-ice stream troughs. The incorporation of our results into ice-sheet models will improve predictions of future sea-level rise.

During recent years interest in the West Antarctic Ice Sheet (WAIS) increased among the geoscientific community because: (1) its bed is located mainly below sea level and drops towards the interior of the Antarctic continent, making the WAIS inherently unstable in a future warming world, (2) extensive thinning and associated grounding-line retreat of the Pine Island and Thwaites glaciers draining into the Amundsen Sea Embayment (ASE) have suggested major changes are underway in the sector, increasing contributions to current sea level rise, (3) a complete collapse of the WAIS in the future would raise global eustatic sea level by 3.4-5 m, while melting of the ASE drainiage basin alone would raise sea level by 1.2-1.5 m. A sea level rise of that magnitude would cause major global socio-economical and ecological problems. Detailed knowledge of the long-term behaviour of the WAIS in the ASE during the recent geological past (Last Glacial Maximum [LGM] to present) will contribute to a better understanding of current glacier dynamics and help to improve numerical ice-sheet models, which aim to predict future sea-level rise. Previous marine geoscientific studies in the ASE focused on the main palaeo-ice stream troughs to reconstruct the LGM extent of the WAIS on the continental shelf and its subsequent retreat history. However, little is known about elevated marinal areas of the palaeo-ice streams in the ASE, where ice retreat most likely lagged behind that in the troughs. Here we present results from multibeam swath bathymetry surveys, high-resolution seismic and sedimentological investigations from the former bed in an inter-ice stream area between the Pine Island-Thwaites palaeo-ice stream and a fast-flow tributary emanating from the area now occupied by the Cosgrove Ice Shelf. The data show an unusual assemblage of glacial morphological features including crevasse-squeeze ridges, large-scale hummocks perpendicular to the palaeo-ice flow, associated recessional moraines, and hill-hole pairs. This combination of bedforms has not been described before from the Antarctic shelf and indicates a more complex ice flow behaviour for the eastern ASE than suggested by the pattern of bedforms in the palaeo-ice stream troughs alone. Our data indicate that slow flowing ice masses covered the topographical highs adjacent to the Pine Island-Thwaites palaeo-ice stream (PITPIS) during the LGM. These ice masses most likely stagnated during a phase of general stillstand of the PITPIS. Here we introduce a six-phase formation model based on these interpretations. New radiocarbon ages indicate a pre-Holocene deglaciation of the inter-ice stream ridge between ~11.6 and 16 ka BP. This new information can be used as a reference dataset for interpreting more inter-ice stream areas in future studies, since they are key areas for stabilising ice streams, and form a large part of the ice sheet in general. New insights into ice dynamics here may help improve ice flow models.

The geological record of landforms and sediments produced beneath deglaciating ice sheets offers insights into inaccessible glacial processes. Large subglacial valleys formed by meltwater erosion of sediments (tunnel valleys) are widespread in formerly glaciated regions such as the North Sea. Obtaining a better understanding of these features may help with the parameterisation of basal melt rates and the interplay between basal hydrology and ice dynamics in numerical models of past, present, and future ice-sheet configurations. However, the mechanisms and timescales over which tunnel valleys form remain poorly constrained. Here, we present a series of numerical modelling experiments, informed by new observations from high-resolution 3D seismic data (6.25 m bin size, ∼4 m vertical resolution), which test different hypotheses of tunnel valley formation and calculate subglacial water routing, seasonal water discharges, and the rates at which tunnel valleys are eroded beneath deglaciating ice sheets. Networks of smaller or abandoned channels, pervasive slump deposits, and subglacial landforms are imaged inside and at the base of larger tunnel valleys, indicating that these tunnel valleys were carved through the action of migrating smaller channels within tens of kilometres of the ice margin and were later widened by ice-contact erosion. Our model results imply that the drainage of extensive surface meltwater to the ice-sheet bed is the dominant mechanism responsible for tunnel valley formation; this process can drive rapid incision of networks of regularly spaced subglacial tunnel valleys beneath the fringes of retreating ice sheets within hundreds to thousands of years during deglaciation. Combined, our observations and modelling results identify how tunnel valleys form beneath deglaciating mid-latitude ice sheets and have implications for how the subglacial hydrological systems of contemporary ice sheets may respond to sustained climate warming.