Hostname: page-component-78c5997874-8bhkd Total loading time: 0 Render date: 2024-11-05T04:19:08.593Z Has data issue: false hasContentIssue false

Small-scale topographically-controlled glacier flow switching in an expanding proglacial lake at Breiðamerkurjökull, SE Iceland

Published online by Cambridge University Press:  31 May 2017

ROBERT D. STORRAR*
Affiliation:
Department of the Natural and Built Environment, Sheffield Hallam University, Sheffield, UK
ANDREW H. JONES
Affiliation:
Department of the Natural and Built Environment, Sheffield Hallam University, Sheffield, UK
DAVID J. A. EVANS
Affiliation:
Department of Geography, Durham University, Durham, UK
*
E-mail: Robert Storrar <[email protected]>
Rights & Permissions [Opens in a new window]

Abstract

Type
Letter
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2017

1. INTRODUCTION

Understanding the internal dynamics of glaciers and ice sheets is important for predicting how they will respond to external forcing factors, such as rising sea level and warming ocean and air temperatures (e.g. Weertman, Reference Weertman1974; Thomas, Reference Thomas1979; Conway and others, Reference Conway2002; Pfeffer, Reference Pfeffer2007; Favier and others, Reference Favier2014; Aschwanden and others, Reference Aschwanden, Fahnestock and Truffer2016). Flow switching, where the trajectory of ice flow significantly alters in time, has been inferred in modern ice sheets (e.g. Retzlaff and Bentley, Reference Retzlaff and Bentley1993; Jacobel and others, Reference Jacobel, Scambos, Raymond and Gades1996; Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997; Conway and others, Reference Conway2002; Hulbe and Fahnestock, Reference Hulbe and Fahnestock2004) and reconstructed in palaeo-ice sheets (e.g. Dowdeswell and others, Reference Dowdeswell, Ottesen and Rise2006; Stokes and others, Reference Stokes, Clark and Storrar2009; Ó Cofaigh and others, Reference Ó Cofaigh, Evans and Smith2010; Stokes and Tarasov, Reference Stokes and Tarasov2010; Winsborrow and others, Reference Winsborrow, Stokes and Andreassen2012; Stokes and others, Reference Stokes, Corner, Winsborrow, Husum and Andreassen2014). Flow switching is associated with changes in ice dynamics (e.g. Payne and Dongelmans, Reference Payne and Dongelmans1997; Conway and others, Reference Conway2002) and could lead to either increases or decreases in velocity. It is therefore of importance to future mass-balance predictions, although few examples of flow switches taking place in modern glaciers have been documented. This letter documents the first observation of a minor flow switch as it happens and provides strong support for a topographic (bathymetric) control.

Winsborrow and others (Reference Winsborrow, Stokes and Andreassen2012) reviewed the literature and proposed five mechanisms that may control flow switching: (1) accumulation of sediment causing accommodation space to be filled (Dowdeswell and others, Reference Dowdeswell, Ottesen and Rise2006); (2) variations in bathymetry causing asynchronous response of ice flow (Stokes and others, Reference Stokes, Clark and Storrar2009); (3) variations in basal thermal regime (Ó Cofaigh and others, Reference Ó Cofaigh, Evans and Smith2010); (4) variations in subglacial meltwater routing and sticky spot location (Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997); and (5) large-scale changes in mass balance (Payne and Dongelmans, Reference Payne and Dongelmans1997). All research into flow switching to date has been carried out on flow switches that have occurred in the past, from scales ranging from the last 150 years (Retzlaff and Bentley, Reference Retzlaff and Bentley1993) to several glacial cycles (Dowdeswell and others, Reference Dowdeswell, Ottesen and Rise2006). Here we document observations of an active small-scale flow switch that is under way at Breiðamerkurjökull, SE Iceland.

2. BREIÐAMERKURJÖKULL

Breiðamerkurjökull is a large (~16 km wide) outlet glacier draining part of the Vatnajökull ice cap in south-east Iceland (Fig. 1). The glacier is composite, comprising four flow units (1–4 from west to east: see Fig. 1a), which are separated by large medial moraines emanating from prominent nunataks. In this paper we focus on the interaction between flow units 3 and 4. The western portion of flow unit 4 terminates in a large proglacial lake named Jökulsárlón and is bounded to the west by a medial moraine named Esjufjallarond. The foreland in this area rises to no more than 50 m above sea level, the most prominent topographic feature of the area being a large trench up to 300 m deep beneath Jökulsárlón, which extends up-glacier for 20 km (Björnsson, Reference Björnsson1996; Björnsson and others, Reference Björnsson, Pálsson and Guðmundsson2001; Fig. 1b). Breiðamerkurjökull has retreated ~6 km since its Little Ice Age maximum position in 1890 and the ice margin is currently retreating over the deepest part of the Jökulsárlón trench (Fig. 1b). Mapping of the glacier and its foreland in 1998 by Evans and Twigg (Reference Evans and Twigg2002) revealed a depression in the ice surface in the vicinity of the deepest part of the Jökulsárlón trench, which is clearly visible in Figure 2a as an area of low ice-surface topography and dense arcuate crevassing.

Fig. 1. (a) 2016 Landsat 8 ETM+ image (7, 3, 2 RGB) showing the location of Breiðamerkurjökull and places referred to in the text. Numbers refer to the flow units. The white box indicates the extent of b. (b) Retreat of Breiðamerkurjökull from 1985 to 2016 superimposed on the elevation of the Jökulsárlón trench, which extends down to deeper shades of blue at ~300 m below sea level (interpolated from contours provided in Björnsson, Reference Björnsson1996).

Fig. 2. Mean surface velocity, elevation and subglacial hydraulic potential changes between 1998 and 2015. Bed DEM below sea-level interpolated from digitised contours from Björnsson and others (Reference Björnsson, Pálsson and Guðmundsson2001). The background image is a 2015 Landsat 8 band 8 image. (a) Landsat-derived mean surface velocity for 1999–2000. Black lines are velocity-derived flowlines. (b) Landsat-derived mean surface velocity for 2014–2015. Black lines are velocity-derived flowlines. (c) Velocity change between a and b. (d) Mean annual surface elevation change between 1998 and 2007. (e) Mean annual surface elevation change between 2007 and 2015. (f) Mean annual surface elevation change between 1998 and 2015. (g) Standard deviation stretch (n = 2.5) of d. (h) Standard deviation stretch (n = 2.5) of e. (i) Standard deviation stretch (n = 2.5) of f. (j) Subglacial flow accumulation in 1998. (k) Subglacial flow accumulation in 2007. (l) Subglacial flow accumulation in 2015.

3. METHODS

Mean annual surface velocity and flow vectors were calculated using an automated normalised cross-correlation feature tracking algorithm (Kääb and Vollmer, Reference Kääb and Vollmer2000). Values were derived from the panchromatic band (15 m resolution) of a pair of Landsat 7 ETM+ scenes (7 October 1999 and 23 September 2000) and a pair of Landsat 8 OLI scenes (8 October 2014 and 25 September 2015).

Individual vertical aerial photographs were used to derive orthophotos and DEMs for 1998 and 2007 using Structure from Motion photogrammetry. The 1998 DEM is 5 m horizontal resolution and is based on the aerial photographs and control points from the (22 August) 1998 glacial geomorphological map of Breiðamerkurjökull (Evans and Twigg, Reference Evans and Twigg2000). The 2007 DEM is 1 m resolution and is based on high resolution (30 cm) NERC ARSF aerial photographs from 3 July 2007 and control points measured in the field (see Storrar and others, Reference Storrar, Evans, Stokes and Ewertowski2015 for details). A 2 m DEM from 10 November 2015, produced by photogrammetry from high-resolution (~0.5 m) optical imagery, was acquired from the ArcticDEM (http://pgc.umn.edu/arcticdem).

Analysis of the terrestrial component of the 2015 DEM reveals that it is horizontally offset from the 1998 and 2007 DEMs by 7 m. It is systematically vertically offset by ~68 m, varying by a range of ~2 m either side of this value. Thus, to measure surface changes, all DEMs were resampled to 50 m horizontal resolution, to conservatively remove any artefact of error in georeferencing. The elevation values of the 2015 DEM were reduced by 68 m to correct for the systematic error. We estimate that the vertical accuracy of the DEMs is ±5 m. Given that the temporal difference between DEMs is 8–17 years, the vertical uncertainty becomes less significant (<1 m a−1) when calculating surface elevation change between years, which was carried out for 1998–2007, 2007–15 and 1998–2015.

A combination of Landsat 4TM, 5TM, 7ETM+ and 8OLI satellite images of the study area from 1985 to 2016 were downloaded from the USGS (http://earthexplorer.usgs.gov/). The temporal resolution is approximately annual from 1985 to 2011, and then monthly (excluding winter) from 2011 to 2016 (details of individual images are provided in the supplementary information). The glacier margin and lateral boundaries of the Esjufjallarond medial moraine were digitised from false colour composite images in the band combination 7, 3, 2 (RGB).

Maps of subglacial hydraulic potential were derived using a GIS approach similar to that of Livingstone and others (Reference Livingstone, Clark, Woodward and Kingslake2013), which is based on Shreve (Reference Shreve1972). Subglacial hydraulic potential (Φ) is derived from elevation and water pressure:

(1) $$\Phi = \rho _{\rm w} gh + P_{\rm w} $$

where ρ w is the density of water (1000 kg m−3), g is the acceleration due to gravity (9.81 m s−2), h is the elevation of the bed and P w is the water pressure. Effective pressure N is the difference between ice pressure P i (917 kg m−3) and subglacial water pressure P w.

(2) $$N = P_{\rm i} - P_{\rm w} $$

P i is controlled by the weight of overlying ice:

(3) $$P_{\rm i} = \rho _{\rm i} gH$$

where H is the ice thickness (m). For simplicity N is assumed to be 0, although this is often not the case, particularly close to the ice margin (cf. Hooke, Reference Hooke1984). However, this simplification is frequently adopted to approximate the generalised pattern of subglacial water flow (e.g. Shreve, Reference Shreve1972; Syverson and others, Reference Syverson, Gaffield and Mickelson1994; Livingstone and others, Reference Livingstone, Clark, Woodward and Kingslake2013), which is the aim in this study. Assuming N to be 0, Eqn (1) can be re-written as:

(4) $$\Phi = \rho _{\rm w} gh + \rho _{\rm i} gH$$

We use the 1998, 2007 and 2015 DEMs for the ice surface, and a DEM interpolated from the subglacial contours published in Björnsson and others (Reference Björnsson, Pálsson and Guðmundsson2001) to derive H and, in turn, Φ. The DEMs were resampled to 200 m resolution to remove the (short-term) effect of ice thickness fluctuations due to crevasses.

4. RESULTS

Mean ice-surface velocity and flow-line patterns (Figs 2a, b) show that areas of higher velocity are generally associated with the deeper parts of the bed. Between 1999–2000 and 2014/15 the fastest flowing parts of the glacier migrated to the west, and a small net increase in velocity was observed in flow unit 4 (Fig. 2c). The direction of flow changed from broadly NW to SE to NNW to SSE, the changes being most notable near to the terminus.

Ice-surface lowering is concentrated around the deepest parts of the Jökulsárlón trench that are coincident with the ice margin (Figs 2d–i). Between 1998 and 2007 this was further to the east than between 2007 and 2015, where the locus of maximum surface lowering was further west. Rates of ice-surface lowering increased after 2007.

Since 1985, Breiðamerkurjökull has undergone asynchronous retreat, with the fastest retreat occurring over the deepest parts of the Jökulsárlón trench (Fig. 1b; Supplementary animation), coincident with the loci of ice-surface lowering. In 1985 the medial moraine Esjufjallarond followed a straight trajectory and terminated on land to the west of Jökulsárlón (Fig. 3b; Supplementary animation). In 2006, a portion of Esjufjallarond close to the terminus began to deflect towards the north-east, indicating flow capture (Fig. 3b). The 2007 DEM shows the early stages of the flow capture: the ice surface in the NE, in the location of the deepest section of the Jökulsárlón trench (Fig. 1b), was ~100 m lower than surrounding areas (Fig. 3a). Esjufjallarond and longitudinal debris folia to the south-west of it, indicate the previous ice flow direction towards the SSE (Fig. 3c). In 2007, Esjufjallarond had not undergone significant displacement, but transverse crevasses oriented SSW to NNE had opened up, indicating a ~45° change in flow towards the east-south-east (Fig. 3c), which is confirmed in the 2015 DEM, which shows only a remnant of Esjufjallarond and extensive crevassing indicating SE to ESE flow. The flow eventually resulted in a large recumbent fold in Esjufjallarond, which was displaced ~300 m by 2016 (Fig. 3). The hinge line of the fold in 2016 was located where ice terminates in Jökulsárlón, while the former terminal section of Esjufjallarond is still on land to the SW (Fig. 3b).

Fig. 3. (a) 2007 SfM DEM of the glacier surface. Black lines demarcate the 2016 margin and medial moraine outlines. The inset graph shows the cross profile indicated by line A-A′ (in metres). The red box indicates the location of c and d. (b) Band 8 Landsat 8 ETM+ image from 13 May 2016 overlain with mapped medial moraine edge positions from 1985 to 2016. Ice margin positions in 1985, 2000, 2007, 2010 and 2016 are shown on the same colour scale. (c) 315° hillshaded 2007 DEM showing the medial moraine (outline in white), longitudinal debris folia (yellow) indicating that the initial flow direction was to the SSE, and transverse crevasses (green) which have opened up indicating the switching of flow to the ESE. (d) 315° hillshaded DEM showing the same area as c in 2015. Note the significant changes compared with the 2007 DEM. The remnant medial moraine is indicated with a red dashed line.

5. DISCUSSION AND IMPLICATIONS

As explained above, it has been suggested that flow switches can be caused by a variety of mechanisms (Winsborrow and others, Reference Winsborrow, Stokes and Andreassen2012). We now assess whether these mechanisms may be responsible for the flow switch observed at Breiðamerkurjökull.

5.1. Accommodation space

Dowdeswell and others (Reference Dowdeswell, Ottesen and Rise2006) suggested that flow switches may occur if accommodation space is filled with sediment. Outwash from the western part of flow unit 4 is deposited in Jökulsárlón (Björnsson and others, Reference Björnsson, Pálsson and Guðmundsson2001), however Jökulsárlón is too large for sediment deposition to be close to filling the available space.

5.2. Changes in basal thermal regime

Breiðamerkurjökull has not undergone significant changes in basal thermal regime over the study period, so we rule out this mechanism.

5.3. Large-scale changes in mass balance

Breiðamerkurjökull has been in a state of negative mass balance since the mid-1990s (Björnsson and Pálsson, Reference Björnsson and Pálsson2008; Pálsson and others, Reference Pálsson, Gunnarsson, Jónsson, Steinþórsson and Pálsson2014). The glacier front has been continuously retreating and the ice surface has been lowering (Figs 1, 2). We suggest that it is likely that this trend of negative mass balance is partly responsible for the observed flow switch, in combination with the subglacial topography and subglacial hydrology.

5.4. Variations in subglacial topography

The capture of part of flow unit 3 by flow unit 4 is likely controlled by bed topography. A large trench up to 300 m below sea level exists along the central axis of Jökulsárlón (Björnsson, Reference Björnsson1996; Björnsson and others, Reference Björnsson, Pálsson and Guðmundsson2001; Fig. 2), the glaciological impacts of which were identified by Evans and Twigg's (Reference Evans and Twigg2002) map of the glacier surface depression (see also Fig. 3a). Retreat and surface lowering have been most marked along this trench (Figs 1b, 2d–i), consistent with the water depth model of iceberg calving, which was modelled at Breiðamerkurjökull by Nick and others (Reference Nick, Van der Kwast and Oerlemans2007). As ice retreats over the Jökulsárlón trench it drives faster ice flow and draws down the ice surface. This is supported by observations of higher ice velocity and surface lowering at the calving front, where the trench is deepest (Fig. 2; Voytenko and others, Reference Voytenko2015). However, it should be noted that calving flux need not necessarily be related to velocity if the glacier is retreating (Benn and others, Reference Benn, Hulton and Mottram2007), which is the case at Breiðamerkurjökull.

It has been demonstrated empirically that water depth exerts a strong control on ice flux, and thereby calving rates (Schoof, Reference Schoof2007) and water depth has been invoked as the mechanism for flow switches in fjord systems (Briner and others, Reference Briner, Bini and Anderson2009), marine-terminating outlet glaciers (Stokes and others, Reference Stokes, Corner, Winsborrow, Husum and Andreassen2014) and palaeo-ice streams (Winsborrow and others, Reference Winsborrow, Stokes and Andreassen2012). We suggest that as Breiðamerkurjökull retreats over the deepest part of the Jökulsárlón trench, flow velocity increases locally where the trench is deepest. This is likely related to calving, which in the water depth model occurs at a greater rate in deeper water (e.g. Nick and others, Reference Nick, Van der Kwast and Oerlemans2007). It is also possible that melting from incursion of warm marine water enhances calving, particularly where cold subglacial meltwater emerges from beneath the glacier (also in the deepest part of the trough: see Figs 2j–l), acting to ‘draw in’ warmer water (Voytenko and others, 2016). Increased flow velocity in turn lowers the ice surface at this location (Fig. 2). The influence of basal topography is highest at the thinnest part of the glacier, where the ratio between water depth and ice thickness is highest, which is why the most significant changes are seen as the terminus retreats over the overdeepening. Ice-surface lowering promotes drawdown of the surrounding ice, which is sufficient to deflect flow from flow unit 3 to the east, into flow unit 4, propagating the switch up-ice. This can be seen by the deviation of flowlines in Figures 2a, b, and in more detail in the deflection of Esjufjallarond (Fig. 3b) and the cross-cutting surface features in the DEM from 2007 (Fig. 3c), which shows the configuration of the glacier during the flow switch. The 2015 DEM of the same area (Fig. 3d) shows the subsequent configuration of the ice-surface post-switch, with extensive crevassing clearly portraying the new flow direction.

5.5. Subglacial hydrology

In addition to calving, subglacial meltwater flow is likely to be heavily influenced by the Jökulsárlón trench. Subglacial meltwater flow is controlled primarily by the bed slope and, to a lesser extent, the ice-surface slope (Shreve, Reference Shreve1972). It has been demonstrated using data from Greenland that subglacial meltwater routing can exhibit temporal variations between catchments and stimulate changes in ice dynamics (Chu and others, 2016). Subglacial hydraulic potential modelling suggests that meltwater beneath Breiðamerkurjökull will drain preferentially into the Jökulsárlón trench (Figs 2j–l). Indeed, strong meltwater currents and suspended sediment plumes indicative of subglacial runoff have been observed near the calving front (Hodgkins and others, Reference Hodgkins, Bryant, Darlington and Brandon2016). Subglacial water pressure in the trench is likely to be high, since the deep water of Jökulsárlón provides a high base-level (Benn and Evans, Reference Benn and Evans2010). High water pressures and preferential routing of water in the Jökulsárlón trench promote increased sliding and may provide an additional explanation for the high velocity of this part of the glacier. Modelled flow accumulation follows the thalweg of the Jökulsárlón trench in all time steps (Figs 2j–l) and so is unlikely to be responsible for the flow switch by itself. Rather, it would have acted to maintain high velocity, the focus of which is ‘steered’ by the basal topography (water depth).The data presented in this paper provide observations from a modern glacier to support the notion that basal topography exerts an important control over ice dynamics (cf. Aschwanden and others, Reference Aschwanden, Fahnestock and Truffer2016). This has implications for understanding the response of ice masses to climatic warming. Basal overdeepenings, which are present beneath many parts of the Greenland (Bamber and others, Reference Bamber2013) and Antarctic (Fretwell and others, Reference Fretwell2013) ice sheets, may alter ice flow vectors and promote flow switching, which may either increase or decrease retreat, depending on the intricacies of the topographic setting. As parts of these ice sheets retreat they will encounter different subglacial topography, making flow switching and concomitant changes in velocity and retreat rates likely. Thus, detailed knowledge of bed topography, and its influence on ice dynamics, is an essential source of information for predicting the fate of water-terminating ice in any setting.

SUPPLEMENTARY MATERIAL

The supplementary material for this article can be found at https://doi.org/10.1017/jog.2017.22.

ACKNOWLEDGEMENTS

We are grateful to Chief Editor Graham Cogley and Scientific Editor Ted Scambos, as well as two anonymous reviewers, whose comments increased the clarity of the work.

References

REFERENCES

Anandakrishnan, S and Alley, RB (1997) Stagnation of ice stream C, West Antarctica by water piracy. Geophys. Res. Lett., 24, 265268 Google Scholar
Aschwanden, A, Fahnestock, MA and Truffer, M (2016) Complex Greenland outlet glacier flow captured. Nat. Commun., 7, 10524 CrossRefGoogle ScholarPubMed
Bamber, J and 9 others (2013) A new bed elevation dataset for Greenland. Cryosphere, 7, 499510 CrossRefGoogle Scholar
Benn, DI and Evans, DJA (2010) Glaciers & glaciation. Hodder Arnold, London Google Scholar
Benn, DI, Hulton, NRJ and Mottram, RH (2007) ‘Calving laws’, ‘sliding laws’ and the stability of tidewater glaciers. Ann. Glaciol., 46, 123130 CrossRefGoogle Scholar
Björnsson, H (1996) Scales and rates of glacial sediment removal: a 20 km long and 300 m deep trench created beneath Breiðamerkurjökull during the Little Ice Age. Ann. Glaciol., 22, 141146 Google Scholar
Björnsson, H and Pálsson, F (2008) Icelandic glaciers. Jökull, 58, 365386 Google Scholar
Björnsson, H, Pálsson, F and Guðmundsson, S (2001) Jökulsárlón at Breiðamerkursandur, Vatnajökull, Iceland: 20th century changes and future outlook. Jökull, 50, 118 Google Scholar
Briner, JP, Bini, AC and Anderson, RS (2009) Rapid early Holocene retreat of a Laurentide outlet glacier through an Arctic fjord. Nat. Geosci., 2, 496499 Google Scholar
Conway, H and 5 others (2002) Switch of flow direction in an Antarctic ice stream. Nature, 419, 465467 CrossRefGoogle Scholar
Dowdeswell, JA, Ottesen, D and Rise, L (2006) Flow switching and large-scale deposition by ice streams draining former ice sheets. Geology, 34, 313316 Google Scholar
Evans, DJA and Twigg, DR (2000) Breiðamerkurjökull 1998. 1: 30,000 Scale Map. University of Glasgow and Loughborough University.Google Scholar
Evans, DJA and Twigg, DR (2002) The active temperate glacial landsystem: a model based on Breiðamerkurjökull and Fjallsjökull, Iceland. Quat. Sci. Rev., 21, 21432177 Google Scholar
Favier, L and 8 others (2014) Retreat of Pine Island Glacier controlled by marine ice-sheet instability. Nat. Clim. Change, 4, 117121 Google Scholar
Fretwell, P and 9 others (2013) Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. Cryosphere, 7, 375393 CrossRefGoogle Scholar
Hodgkins, R, Bryant, R, Darlington, E and Brandon, M (2016) Pre-melt-season sediment plume variability at Jökulsárlón, Iceland, a preliminary evaluation using in-situ spectroradiometry and satellite imagery. Ann. Glaciol., 57, 18 Google Scholar
Hooke, RL (1984) On the role of mechanical energy in maintaining subglacial water conduits at atmospheric pressure. J. Glaciol., 30, 180187 Google Scholar
Hulbe, CL and Fahnestock, MA (2004) West Antarctic ice-stream discharge variability: mechanism, controls and pattern of grounding-line retreat. J. Glaciol., 50, 471484 CrossRefGoogle Scholar
Jacobel, RW, Scambos, TA, Raymond, CF and Gades, AM (1996) Changes in the configuration of ice stream flow from the West Antarctic Ice Sheet. J. Geophys. Res.: Solid Earth, 101, 54995504 Google Scholar
Kääb, A and Vollmer, M (2000) Surface geometry, thickness changes and flow fields on creeping mountain Permafrost: automatic extraction by digital image analysis. Permafrost Periglacial Process., 11, 315326 Google Scholar
Livingstone, SJ, Clark, CD, Woodward, J and Kingslake, J (2013) Potential subglacial lakes and meltwater drainage pathways beneath the Antarctic and Greenland ice sheets. Cryosphere, 7, 17211740 CrossRefGoogle Scholar
Nick, FM, Van der Kwast, J and Oerlemans, J (2007) Simulation of the evolution of Breidamerkurjökull in the late Holocene. J. Geophys. Res.-Solid Earth, 112, B01103/01101B01103/01112 Google Scholar
Ó Cofaigh, C, Evans, DJA and Smith, IR (2010) Large-scale reorganization and sedimentation of terrestrial ice streams during late Wisconsinan Laurentide Ice Sheet deglaciation. Geol. Soc. Am. Bull., 122, 743756 Google Scholar
Pálsson, F, Gunnarsson, A, Jónsson, Þ, Steinþórsson, S and Pálsson, HS (2014) Vatnajökull: mass balance, meltwater drainage and surface velocity of the glacial year 2012–13. Landsvirkjun, LV-2014-068Google Scholar
Payne, A and Dongelmans, P (1997) Self-organization in the thermomechanical flow of ice sheets. J. Geophys. Res.: Solid Earth, 102, 1221912233 Google Scholar
Pfeffer, WT (2007) A simple mechanism for irreversible tidewater retreat. J. Geophys. Res., 112, F03S25 Google Scholar
Retzlaff, R and Bentley, CR (1993) Timing of stagnation of Ice Stream C, West Antarctica, from short-pulse radar studies of buried surface crevasses. J. Glaciol., 39, 553561 Google Scholar
Schoof, C (2007) Ice sheet grounding line dynamics: steady states, stability, and hysteresis. J. Geophys. Res.: Earth Surf. (2003–2012), 112 Google Scholar
Shreve, RL (1972) Movement of water in glaciers. J. Glaciol., 11, 205214 Google Scholar
Stokes, CR and Tarasov, L (2010) Ice streaming in the Laurentide Ice Sheet: a first comparison between data-calibrated numerical model output and geological evidence. Geophys. Res. Lett., 37, L01501 Google Scholar
Stokes, CR, Clark, CD and Storrar, RD (2009) Major changes in ice stream dynamics during deglaciation of the north-western margin of the Laurentide Ice Sheet. Quat. Sci. Rev., 28, 721738 Google Scholar
Stokes, CR, Corner, GD, Winsborrow, MC, Husum, K and Andreassen, K (2014) Asynchronous response of marine-terminating outlet glaciers during deglaciation of the Fennoscandian Ice Sheet. Geology, 42, 455458 Google Scholar
Storrar, RD, Evans, DJA, Stokes, CR and Ewertowski, M (2015) Controls on the location, morphology and evolution of complex esker systems at decadal timescales, Breiðamerkurjökull, SE Iceland. Earth Surf. Process. Landforms, 40, 14211438 Google Scholar
Syverson, KM, Gaffield, SJ and Mickelson, DM (1994) Comparison of esker morphology and sedimentology with former ice-surface topography, Burroughs Glacier, Alaska. Geol. Soc. Am. Bull., 106, 11301142 Google Scholar
Thomas, RH (1979) The dynamics of marine ice sheets. J. Glaciol., 24, 167177 Google Scholar
Voytenko, D and 7 others (2015) Multi-year observations of Breiðamerkurjökull, a marine-terminating glacier in southeastern Iceland, using terrestrial radar interferometry. J. Glaciol., 61, 4254 CrossRefGoogle Scholar
Weertman, J (1974) Stability of the junction of an ice sheet and an ice shelf. J. Glaciol., 13, 311 Google Scholar
Winsborrow, MC, Stokes, CR and Andreassen, K (2012) Ice-stream flow switching during deglaciation of the southwestern Barents Sea. Geol. Soc. Am. Bull., 124, 275290 Google Scholar
Figure 0

Fig. 1. (a) 2016 Landsat 8 ETM+ image (7, 3, 2 RGB) showing the location of Breiðamerkurjökull and places referred to in the text. Numbers refer to the flow units. The white box indicates the extent of b. (b) Retreat of Breiðamerkurjökull from 1985 to 2016 superimposed on the elevation of the Jökulsárlón trench, which extends down to deeper shades of blue at ~300 m below sea level (interpolated from contours provided in Björnsson, 1996).

Figure 1

Fig. 2. Mean surface velocity, elevation and subglacial hydraulic potential changes between 1998 and 2015. Bed DEM below sea-level interpolated from digitised contours from Björnsson and others (2001). The background image is a 2015 Landsat 8 band 8 image. (a) Landsat-derived mean surface velocity for 1999–2000. Black lines are velocity-derived flowlines. (b) Landsat-derived mean surface velocity for 2014–2015. Black lines are velocity-derived flowlines. (c) Velocity change between a and b. (d) Mean annual surface elevation change between 1998 and 2007. (e) Mean annual surface elevation change between 2007 and 2015. (f) Mean annual surface elevation change between 1998 and 2015. (g) Standard deviation stretch (n = 2.5) of d. (h) Standard deviation stretch (n = 2.5) of e. (i) Standard deviation stretch (n = 2.5) of f. (j) Subglacial flow accumulation in 1998. (k) Subglacial flow accumulation in 2007. (l) Subglacial flow accumulation in 2015.

Figure 2

Fig. 3. (a) 2007 SfM DEM of the glacier surface. Black lines demarcate the 2016 margin and medial moraine outlines. The inset graph shows the cross profile indicated by line A-A′ (in metres). The red box indicates the location of c and d. (b) Band 8 Landsat 8 ETM+ image from 13 May 2016 overlain with mapped medial moraine edge positions from 1985 to 2016. Ice margin positions in 1985, 2000, 2007, 2010 and 2016 are shown on the same colour scale. (c) 315° hillshaded 2007 DEM showing the medial moraine (outline in white), longitudinal debris folia (yellow) indicating that the initial flow direction was to the SSE, and transverse crevasses (green) which have opened up indicating the switching of flow to the ESE. (d) 315° hillshaded DEM showing the same area as c in 2015. Note the significant changes compared with the 2007 DEM. The remnant medial moraine is indicated with a red dashed line.

Storrar supplementary material

Storrar supplementary material 1

Download Storrar supplementary material(Video)
Video 25.6 MB