Reconstructing subglacial meltwater dynamics from the spatial and temporal variation in the form and pattern of eskers

Meltwater drainage beneath glaciers and ice sheets is intimately linked to their dynamics. Meltwater may increase ice velocity if it acts to lubricate the bed; conversely, an efficient subglacial meltwater drainage system may preclude meltwater induced ice acceleration by limiting the amount of water available to facilitate sliding. Thus, understanding the nature of meltwater flow beneath ice masses is crucial for predicting how ice sheets and glaciers will react to increased meltwater input. However, direct observation of subglacial meltwater drainage systems is extremely difficult, meaning that indirect methods such as remote sensing, numerical modelling, dye tracing and geophysical survey are the only way to observe this environment. These methods often suffer from excessive uncertainty and poor spatial and, particularly, temporal resolution. This thesis presents the results of an alternative approach, using the geomorphological record of eskers to understand the former behaviour of meltwater beneath the Laurentide Ice Sheet (LIS) in Canada, and at Breiðamerkurjӧkull in Iceland. Eskers are elongate, sinuous ridges of glaciofluvial sand and gravel deposited in glacial drainage channels. Despite a large body of research on eskers, no systematic analysis of the large-scale properties of eskers, or the implications this may have for understanding subglacial meltwater, has yet been undertaken. Eskers are mapped at the ice sheet (continental) scale in Canada from 678 Landsat ETM+ images and at high resolution (~30 cm) from 407 aerial photographs of the Breiðamerkurjӧkull foreland, in order to address three outstanding questions: (i) What controls the pattern and morphology of eskers? (ii) How did subglacial drainage systems evolve during ice sheet deglaciation? (iii) How can eskers be used to further our understanding of subglacial hydrology? Over 20,000 eskers are mapped in Canada, revealing that esker systems are up to 760 km long, and are surprisingly straight. The spacing between eskers is relatively uniform and they exhibit little change in elevation from one end to another. As the LIS deglaciated between 13 cal ka and 7 cal ka, eskers increased in frequency, which is interpreted to represent an increase in meltwater drainage in channelized, rather than distributed, systems. Eskers are abundant over the resistant rocks of the Canadian Shield and also show a strong preference for formation in areas covered with till. Esker length, sinuosity and spacing appear to be unrelated to the underlying geology. Finally, two types of complex esker systems are proposed: esker fan complexes and topographically constrained esker complexes. The formation of esker complexes is dependent on sediment and meltwater supply and the pre-existing topography controls the overall shape of the esker systems

A model for interaction between conduits and surrounding hydraulically connected distributed drainage based on geomorphological evidence from Keewatin, Canada

© 2020 Author(s). We identify and map visible traces of subglacial meltwater drainage around the former Keewatin Ice Divide, Canada, from high-resolution Arctic Digital Elevation Model (ArcticDEM) data. We find similarities in the characteristics and spatial locations of landforms traditionally treated separately (i.e. meltwater channels, meltwater tracks and eskers) and propose that creating an integrated map of meltwater routes captures a more holistic picture of the large-scale drainage in this area. We propose the grouping of meltwater channels and meltwater tracks under the term meltwater corridor and suggest that these features in the order of 10s-100sm wide, commonly surrounding eskers and transitioning along flow between different types, represent the interaction between a central conduit (the esker) and surrounding hydraulically connected distributed drainage system (the meltwater corridor). Our proposed model is based on contemporary observations and modelling which suggest that connections between conduits and the surrounding distributed drainage system within the ablation zone occur as a result of overpressurisation of the conduit. The widespread aerial coverage of meltwater corridors (5%-36% of the bed) provides constraints on the extent of basal uncoupling induced by basal water pressure fluctuations. Geomorphic work resulting from repeated connection to the surrounding hydraulically connected distributed drainage system suggests that basal sediment can be widely accessed and evacuated by meltwater

Glacial geomorphology of the northern Kivalliq region, Nunavut, Canada, with an emphasis on meltwater drainage systems

This paper presents a glacial geomorphological map of glacial lineations, ribbed terrain, moraines, meltwater channels (subglacial and ice-marginal/proglacial), eskers, glaciofluvial deposits, ice-contact outwash fans and deltas and abandoned shorelines on the bed of the former Laurentide Ice Sheet in northern Canada. Mapping was compiled from satellite imagery and digital elevation data and landforms were digitised directly into a Geographical Information System. The map reveals a complex glacial history characterised by multiple ice-flow events, including fast-flowing ice-streams. Moraines record a series of pauses or re-advances during overall SE retreat towards the Keewatin Ice Divide. The distribution of subglacial meltwater landforms indicates that several distinctive scales and modes of drainage system operated beneath the retreating ice sheet. This includes a large (>100 km) integrated network of meltwater channels, eskers, ice-contact outwash fans and deltas and glaciofluvial deposits that originates at the northern edge of Aberdeen Lake. The map comprises zone 66 of the Canadian National Topographic System, which encompasses an area of 160,000 km2. It is presented at a scale of 1:500,000 and is designed to be printed at A0 size

3D morphometry of De Geer Moraines and Crevasse-Squeeze Ridges: Differentiating between pushing and squeezing mechanisms from remotely sensed data

De Geer Moraines (DGM) and Crevasse-Squeeze Ridges (CSR) are important landforms that can provide useful insights regarding palaeo-glacial processes. Specifically, these landforms can provide information concerning ice-marginal dynamics, and/or subglacial processes, depending on the context in which they are formed. The extraction of 3D morphometric data from these ridges can help to elucidate their formational processes, and potentially enable landform differentiation. We develop a new Python-based ArcGIS toolbox that can automatically extract 3D morphometric data from large sample sets of linear features. The morphometry toolbox may be applied to a wide range of research disciplines that are concerned with quantifying the morphometry of any elongated landforms. This is particularly useful for DGM and CSR studies, where visual similarities can result in confusion over landform type and/or formation. Here we present a case study from southwest Finland and the Northwest Territories, Canada, whereby high-resolution 3D morphometric data is used to analyse and classify DGMs and CSRs. The results reveal key differences in morphometric properties between the landforms which enables a quantified foundation by which to differentiate them. The studied CSRs are found to be higher, wider, steeper, more symmetrical, less sinuous and more voluminous than the studied prominent DGM. In contrast, a tendency for cross-sectional asymmetry in DGM supports an origin by ice-marginal pushing, rather than basal squeeze-up into crevasses. This is further supported by CSRs being less sinuous than DGM due to them being constrained to the dimensions and planform of the (relatively straight) host crevasses, whereas DGM follow a more sinuous path related to the ice margin shape. Future work should include sedimentological and geophysical studies to constrain DGM internal architecture and formation processes. The results may then be used to validate the application of DGM for detailed ice marginal reconstructions

A quasi-annual record of time-transgressive esker formation: implications for ice sheet reconstruction and subglacial hydrology

We identify and map chains of esker beads (series of aligned mounds) up to 15 m high and on average ~ 65 m wide across central Nunavut, Canada from the high-resolution (2 m) ArcticDEM. Based on the close one-to-one association with regularly spaced, sharp crested ridges interpreted as De Geer moraines, we interpret the esker beads to be quasi-annual ice-marginal deposits formed time-transgressively at the mouth of subglacial conduits during deglaciation. Esker beads therefore preserve a high-resolution record of ice-margin retreat and subglacial hydrology. The well-organised beaded esker network implies that subglacial channelised drainage was relatively fixed in space and through time. Downstream esker bead spacing constrains the typical pace of deglaciation in central Nunavut between 7.2 and 6 ka 14C BP to 165–370 m yr−1, although with short periods of more rapid retreat (> 400 m yr−1). Under our time-transgressive interpretation, the lateral spacing of the observed eskers provides a true measure of subglacial conduit spacing for testing mathematical models of subglacial hydrology. Esker beads also record the volume of sediment deposited in each melt season, thus providing a minimum bound on annual sediment fluxes, which is in the range of 103–104 m3 yr−1 in each 6–10 km wide subglacial conduit catchment. We suggest the prevalence of esker beads across this predominantly marine terminating sector of the former Laurentide Ice Sheet is a result of sediment fluxes that were unable to backfill conduits at a rate faster than ice-margin retreat. Esker ridges, conversely, are hypothesised to form when sediment backfilling of the subglacial conduit outpaced retreat resulting in headward esker growth close to but behind the margin. The implication, in accordance with recent modelling results, is that eskers in general record a composite signature of ice-marginal drainage rather than a temporal snapshot of ice-sheet wide subglacial drainage

Eskers on Mars: Morphometric comparisons to eskers on Earth and implications for sediment-discharge dynamics of subglacial drainage

Mars’ present climate is extremely cold and arid. Until recently, it was widely thought that debris-covered glaciers in Mars’ mid-latitudes have been pervasively cold-based since their formation 10s–100s Myr ago. However, we recently discovered eskers associated with ~110–150 Myr old glaciers in the Phlegra Montes [1] and NW Tempe Terra [2] regions of Mars’ northern mid-latitudes. Eskers are sinuous ridges comprising sediments deposited in glacial meltwater conduits. Therefore, eskers associated with existing mid-latitude glaciers on Mars indicate that localised wet-based glaciation did occur during Mars’ most recent geological period. Eskers are important tools for reconstructing the nature, extent, and dynamics of wet-based glaciation on Earth, and have similar potential for Mars. We used 1–2 m/pixel resolution digital elevation models derived from 25–50 cm/pixel High Resolution Imaging Science Experiment stereo-pair images to measure the planform and 3D morphometries of the Phlegra Montes and NW Tempe Terra eskers, and compare them with the morphometries of Quaternary-aged eskers in Canada [3] and SW Finland [4]. We found that the Martian eskers have remarkably similar lengths, sinuosities and heights to terrestrial eskers, but that the Martian eskers are typically wider and have lower side slopes. Large width-height ratios of the Martian eskers are consistent with our previous measurements of ancient (~3.5 Ga) eskers close to Mars’ south pole [5], and may arise from differences in either: esker degradation state, or fundamental glacio-hydrological controls on esker formation between Mars and Earth. Portions of the two Martian eskers with comparable crest morphologies (e.g., sharp- or round-crested) have similar width-height relationships, suggesting that glacio-hydrological processes may exert controls upon the observed relationships between esker morphology and morphometry. Our morphometric analyses also reveal that the Martian esker in NW Tempe Terra has a ‘stacked’ morphology: the crest of a wide, round-crested underlying ridge is superposed by a narrow, sharp- to multi-crested ridge. Based on morpho-sedimentary relationships observed along terrestrial eskers [6], we interpret this transition to represent waning sediment supply and meltwater discharge towards the end of the esker-forming drainage episode(s). Direct sedimentary insights into Martian eskers are not yet possible so we emphasise that such inferences should be rigorously grounded in observations of analogous landforms on Earth. This work was funded by STFC grant ST/N50421X/1. References: [1] Gallagher, C., and Balme, M.R., (2015), Earth. Planet. Sci. Lett. 431, 96-109, [2] Butcher, F.E.G., et al. (2017), J. Geophys. Res. Planets. 122(12), 2445-2468, [3] Storrar, R.D., et al. (2014) Quat. Sci. Rev. 105, 1-25, [4] Storrar, R.D., and Jones, A., Unpublished, [5] Butcher, F.E.G., et al. (2016), Icarus 275, 65-84, [6] Burke, M.J., et al. (2010) Geol. Soc. Am. Bull. 122, 1637-1645

Post-little ice age glacial geomorphology of contrasting topographic settings at Skálafellsjökull, southeast Iceland

Glacial geomorphological mapping from the southern margin of Skálafellsjökull, southeast Iceland, depicts a topographically diverse mountainside, influencing glacier dynamics, landform formation and glacier retreat since the Little Ice Age maximum in ∼1890. The glacial landforms present are typical of southeast Icelandic temperate glaciers, comprising recessional push moraines, including sawtooth moraines, and associated fluting. Study area A demonstrates an abandoned lobe confined by steep V-shaped topography, displaying moraines and minimal fluting, suggesting low preservation of landforms, and changes in glacier behaviour. At study area B, the sawtooth moraine morphology demonstrates changes in the glacier margin as the ice interacted with a series of topographic benches during active recession. The steep-sided valley at study area C illustrates densely spaced arcuate moraines, reflecting subtle changes in ice elevation. This mapping provides a framework for further investigations into glacier retreat rates and the influence of local topography and climate