Dr Jenna Sutherland

Ice-marginal lakes influence the dynamic behaviour of glaciers and ice sheets, impacting the rate at which they lose mass. In Greenland, accelerated ice loss over recent decades had led to an increase in the number of lakes bordering the ice sheet margin. This landscape evolution has sparked a growing field of research focused on quantitatively understanding the interactions between lakes and glaciers, so that ice-marginal lakes can be accounted for in models of ice sheet change. Ice loss from the Greenland Ice Sheet directly contributes to global sea level rise; understanding the drivers of this mass loss is important for accurately predicting future sea level. This article outlines recent advances in our understanding of lake−glacier interactions across Greenland during the past, present and future, and discusses key priorities for further research. We conclude by suggesting a series of activities that introduce Post-16 students to relevant datasets and techniques.

Coastal environments in the Arctic are increasingly affected by the rapidly changing climate. Significant and complex impacts of atmospheric warming have been intensifying, with changes observed both in terrestrial and marine environments. Here, we describe the overview and highlight the study results of multidisciplinary research activities performed under the ArCS II project (Arctic Challenge for Sustainability II) in the Qaanaaq coastal region of northwestern Greenland. The Japanese Arctic projects GRENE-Arctic and ArCS have conducted research at this study site since 2012. In continuity with these previous efforts, field and satellite measurements were carried out to quantify glacier and ice sheet changes. Fish, marine mammals and seabirds, which are key natural resources to human livelihoods, were studied in collaboration with local fishermen and hunters to examine habitat use and clarify the potential responses of marine ecosystems to the changing environments. Greenlandic villages are also directly affected by the flooding of glacial streams and landslides, which were monitored to better understand the driving mechanisms and risks to Arctic societies in the future. Research was also carried out in Qaanaaq village to investigate waste management and housing conditions. The study results were shared with residents through workshops that took place in Qaanaaq and nearby smaller villages. Our results show that coastal environments in northwestern Greenland are changing with increasingly evident impact on human livelihoods. Further collaboration with the villagers, notably in co-designing research questions and interests, is crucial to anticipate, reduce and mitigate the impacts of environmental changes on Arctic communities.

Nascent peatlands represent an emerging, nature‐based carbon sink in the global climate system. A warming climate and changing precipitation regime could drive peat initiation beyond the current latitudinal and altitudinal boundaries of the peatland bioclimatic envelope, through increases in plant productivity and moisture availability, with potential implications for global radiative forcing. However, contemporaneous observations of new peat formation remain scarce. We investigate peat initiation within the deglaciating Rob Roy valley in the Southern Alps, Aotearoa/New Zealand. We find that montane peats have developed across the head of the valley since ∼1949 C.E., coinciding with regional climate warming and glacial retreat. Further, we identify a common ecological succession, characterized by a rise in brown mosses (mainly Bryum) beginning around ∼1963 C.E. Our findings indicate the potential for wider peat expansion in increasingly warm and wet montane landscapes. However, further bioclimatic modeling is required to elucidate where future peatland developments may occur.

We present the first systematic inventory of surge-type glaciers for the whole of Greenland compiled from published datasets and multitemporal satellite images and digital elevation models. The inventory allows us to define the spatial and climatic distribution of surge-type glaciers and to analyse the timing of surges from 1985 to 2019. We identified 274 surge-type glaciers, an increase of 37% compared to previous work. Mapping surge-type glacier distribution by temperature and precipitation variables derived from ERA5-Land reanalysis data shows that the west and east clusters occur in well-defined climatic envelopes. Analysis of the timing of surge active phases during the periods ~1985 to 2000 (T1) and ~2000 to 2019 (T2) suggests that overall surge activity is similar in T1 and T2, but there appears to be a reduction in surging in the west cluster in T2. Our climate analysis shows a coincident increase in mean annual and mean winter air temperature between T1 and T2. We suggest that as glaciers thin under current warming, some surge-type glaciers in the west cluster may be being prevented from surging due to (1) their inability to build-up sufficient mass and (2) a switch from a polythermal to a largely cold-based thermal regime.

The centennial response of land-terminating glaciers in Greenland to climate change is largely unknown. Yet, such information is important to understand ongoing changes and for projecting the future evolution of Arctic subpolar glaciers, meltwater runoff, and sediment fluxes. This paper analyses the topography, geomorphology, and sedimentology of prominent moraine ridges and the proglacial areas of ice cap outlet glaciers on the Qaanaaq peninsula (Piulip Nunaa). We determine geometric changes of glaciers since the neoglacial maximum; the Little Ice Age (LIA), and we compare glacier behaviour during the LIA with that of the present day. There has been very little change in the rate of volume loss of each outlet glacier since the LIA compared with the rate between 2000 and 2019. However, the percentage of each glacier that is likely composed of cold-based ice has increased since the LIA, typically by 20%. The LIA moraines comprise subrounded, striated, and faceted clasts that evidence subglacial transport, and outwash plains, flutes, kames, and eskers that evidence subglacial motion and meltwater within temperate ice. Contrastingly, contemporary ice margins and their convex ice surfaces comprise pronounced primary foliation, ephemeral supraglacial drainage, sediment drapes from thrust plane fractures, and an absence of open crevasses and moulins. These calculations and observations together lead us to interpret that these outlet glaciers have transitioned towards an increasingly cold-based thermal regime despite a warming regional climate. Thermal regime transitions control glacier dynamics and therefore should be incorporated into glacier evolution models, especially where polythermal glaciers prevail and where climate is changing rapidly.

Glaciers and ice caps (GICs) are important contributors of meltwater runoff and to global sea level rise. However, knowledge of GIC mass changes is largely restricted to the last few decades. Here we show the extent of 5327 Greenland GICs during Little Ice Age (LIA) termination (1900) and reveal that they have fragmented into 5467 glaciers in 2001, losing at least 587 km3 from their ablation areas, equating to 499 Gt at a rate of 4.34 Gt yr−1. We estimate that the long-term mean mass balance in glacier ablation areas has been at least −0.18 to −0.22 m w.e. yr−1 and note the rate between 2000 and 2019 has been three times that. Glaciers with ice-marginal lakes formed since the LIA termination have had the fastest changing mass balance. Considerable spatial variability in glacier changes suggest compounding regional and local factors present challenges for understanding glacier evolution.

The extent of the Southern Alps icefield in New Zealand is well-constrained chronologically for the last glacial cycle. The sediment-landform imprint of this glacial system, however, offers insight into ice-marginal processes that chronological control cannot. We present the first detailed investigation of sediments along the southwestern shores of Lake Tekapo, South Island. We identify seven lithofacies, from which a five-stage palaeoglaciological reconstruction of depositional and glaciotectonic events is proposed: (i) ice-marginal advance and deposition of outwash gravels in lithofacies (LF) 1; (ii) ice-marginal recession and the development of an ice-contact lake, manifest in rhythmite deposition and iceberg rafting of dropstones (LF 2), followed by a depositional hiatus; (iii) ice-marginal recession recorded in ice-proximal aggradation of glaciofluvial hyperconcentrated flows (LFs 3, 4); (iv) ice-marginal advance documented by glaciotectonic disturbance and localized hydrofracturing, coeval with the deposition of delta foresets and a subglacial diamicton/till (LFs 5, 6); (v) final stages of ice-marginal recession and deposition of outwash gravels in LF 7. Two infrared stimulated luminescence ages were obtained from the glaciolacustrine sediments and, whilst the dating has some limitations, the sediments pre-date both the global and local Last Glacial Maximum. Overall, this sequence, consistent with sediment fills recorded elsewhere across South Island, suggests recurrence of processes from different glacial advances and the role of topographic constraints on maintaining lake positions.

Ice-marginal lakes can affect glacier dynamics but are ignored in studies of the evolution of the Greenland ice sheet (GrIS) and of peripheral mountain glaciers and ice caps (PGICs). Here we show that lakes occupy 10% of the GrIS ice margin and occur on 5% of PGICs. Ice velocity at the GrIS margin is enhanced by ∼ 25% at lakes versus on land. Mean ice discharge into lakes is ∼4.9 Gt.yr, which is ∼1% of ice discharged through marine termini. We locate thousands of subglacial overdeepenings within which 7,404 km2 of future lakes could form, all of which will be ice-marginal at some time. Future lakes in the west and east will be restricted to the margin of the GrIS and within alpine valleys, respectively. This status and possible future leads us to contend that lakes should be incorporated into projections of Greenland ice loss.

Global climate change is evidently manifest in disappearing mountain glaciers and receding and thinning ice sheet margins. Concern about contemporaneous proglacial lake development has spurred an emerging area of research seeking to quantitatively understand lake - glacier interactions. This perspectives article identifies spatio-temporal disparity between the coverage of field data, remote sensing observations and numerical modeling efforts. Throughout, an overview of the physical effects of lakes on glaciers and on ice sheet margins is provided, drawing evidence together from very recent and high-impact studies of both modern glaciology and of the Quaternary record. We identify and discuss six challenges for numerical modeling of lake - glacier interactions, namely that there are meltwater exchanges between glaciers and ice-marginal lakes, lake bathymetry and glacier bed topography are often unknown, lake - glacier interactions affect the longitudinal stress regime of glaciers, lake water temperature affects glacier melting but is very poorly constrained, the interactions persist with considerable spatio-temporal variability and with boundary migration, and data for model parameterization and validation is extremely scarce. Overall, we contend that numerical modeling is a key frontier in the cryospheric sciences to deliver process understanding of lake - glacier interactions.

Abstract

Ice‐contact proglacial lakes are generally absent from numerical model simulations of glacier evolution, and their effects on ice dynamics and on rates of deglaciation remain poorly quantified. Using the BISICLES ice flow model, we analyzed the effects of an ice‐contact lake on the Pukaki Glacier, New Zealand, during recession from the Last Glacial Maximum. The ice‐contact lake produced a maximum effect on grounding line recession >4 times further and on ice velocities up to 8 times faster, compared to simulations of a land‐terminating glacier forced by the same climate. The lake contributed up to 82% of cumulative grounding line recession and 87% of ice velocity during the first 300 years of the simulations, but those values decreased to just 6% and 37%, respectively, after 5,000 years. Numerical models that ignore lake interactions will, therefore, misrepresent the rate of recession especially during the transition of a land‐terminating to a lake‐terminating environment.

Proglacial lakes can affect the stability of mountain glaciers and can partly disengage glacier behaviour from climatic perturbations. However, their role in controlling the onset and progression of deglaciation from the Last Glacial Maximum (LGM) remains poorly understood. This lack of understanding is partly because the evidence required to consistently and robustly identify the location and evolution of ice-contact lakes is not standardised. In this paper we therefore firstly present a new set of criteria for identifying the landform and sedimentary evidence that defines and characterises ice-marginal lakes. Secondly, we then apply these key criteria with the aid of high-resolution topographic mapping to produce the first holistic definition and assessment of major proglacial lake landforms and sediments pertaining to the end of the LGM across South Island, New Zealand. The major findings of this assessment can be grouped to include that: (i) The localised constraints to proglacial lake extent were topography, glacier size and meltwater/sediment fluxes, (ii) Lake damming was initiated by outwash fan-heads that interrupted water and sediment flows down-valley, and (iii) New Zealand LGM lakes were unequivocally in contact with a calving ice margin. These findings will be useful for reconstructing ice dynamics and landscape evolution in this region.

Quaternary glaciations have created impressive landform assemblages that can be used to understand palaeo-glacier extent, character and behaviour, and hence past global and local glacier forcings. However, in the southern hemisphere and especially in New Zealand, the Quaternary glacial landform record is relatively poorly investigated with regard to glaciological properties. In this study, a 1 m digital elevation model (DEM) was generated from airborne LiDAR data and supplemented with aerial imagery and field observations to analyse the exceptionally well-preserved glacial geomorphology surrounding Lake Tekapo, New Zealand. We describe a rich suite of Last Glacial Maximum (LGM) and recessional ice-marginal, subglacial, supraglacial, glaciofluvial and glaciolacustrine landform assemblages. These represent two landsystems comprising i) fluted till surfaces with low-relief push moraine ridges; and ii) crevasse-squeeze ridges, ‘zig-zag’ eskers and attenuated lineations. The former landsystem records the behaviour of an active temperate glacier and the latter landsystem, which is superimposed upon and inset within the former, strongly suggests intermittent surge phases. The two landsystem signatures indicate a sequential change in ice-marginal dynamics during recession that was likely to have been partially non-climatically driven. Overall, we present the first evidence of surge-type glacier behaviour in New Zealand.

The central sector of the British-Irish Ice Sheet during the last glaciation was characterised by complex ice-flow reflecting interacting ice streams and changing dominance of different ice dispersal centres. At Tunstall, east Yorkshire, two subglacial till units have been traditionally identified as the Late Devensian Skipsea and Withernsea tills, and thought to record two separate ice advances onto the Holderness coast, from divergent ice flow directions. Our study presents the first quantitative lithological, sedimentological and structural evaluation of glacial sediments at the site. The lithological composition of both till units suggests that ice extended southwards from southern Scotland, incorporating material from north-east England and the western margin of the North Sea Basin. Notably, the bulk lithological properties of both the Skipsea and Withernsea tills are very similar. Subtle variations in colour, texture and lithology that do occur simply appear to reflect spatial and temporal variability in subglacial entrainment along the flow path of the North Sea Lobe. The relative arrangements of the units plus the fracture sets also indicates phases of intra-till thrust-stacking and unloading (F2), consolidation and shrinkage (F1, F3) suggestive of cycles of ice re-advance (thrusting) and ice-marginal retreat (unloading and shrinkage) possibly relating to active recession. The findings from this study reveal a sedimentary and structural complexity that is not recognised by the current Late Devensian till stratigraphy of east Yorkshire.

Abstract

Rapid changes observed today in mountain glaciers need to be put into a longer-term context to understand global sea-level contributions, regional climate-glacier systems and local landscape evolution. In this study we determined volume changes for 400 mountain glaciers across the Southern Alps, New Zealand for three time periods; pre-industrial “Little Ice Age (LIA)” to 1978, 1978 to 2009 and 2009 to 2019. At least 60 km3 ± 12 km3 or between 41 and 62% of the LIA total ice volume has been lost. The rate of mass loss has nearly doubled from − 0.4 m w.e year−1 during 1,600 to 1978 to − 0.7 m w.e year−1 at present. In comparison Patagonia has lost just 11% of it’s LIA volume. Glacier ice in the Southern Alps has become restricted to higher elevations and to large debris-covered ablation tongues terminating in lakes. The accelerating rate of ice loss reflects regional-specific climate conditions and suggests that peak glacial meltwater production is imminent if not already passed, which has profound implications for water resources and riverine habitats.

Global glacier mass loss is causing expansion of proglacial landscapes and producing meltwater that can become impounded as lakes within natural topographic depressions or ‘overdeepenings’. It is important to understand the evolution of these proglacial landscapes for water resources, natural hazards and ecosystem services. In this study we (i) overview contemporary loss of glacier ice across the Southern Alps of New Zealand, (ii) analyse ice-marginal lake development since the 1980s, (iii) utilise modelled glacier ice thickness to suggest the position and size of future lakes, and (iv) employ a large-scale glacier evolution model to suggest the timing of future lake formation and future lake expansion rate. In recent decades, hundreds of Southern Alps glaciers have been lost and those remaining have fragmented both by separation of tributaries and by detachment of ablation zones. Glaciers with ice-contact margins in proglacial lakes (n > 0.1 km2 = 20 in 2020) have experienced the greatest terminus retreat and typically twice as negative mass balance compared to similar-sized land-terminating glaciers. Our analysis indicates a positive relationship between mean glacier mass balance and rate of lake growth (R2 = 0.34) and also with length of an ice-contact lake boundary (R2 = 0.44). We project sustained and relatively homogenous glacier volume loss for east-draining basins but in contrast a heterogeneous pattern of volume loss for west-draining basins. Our model results show that ice-marginal lakes will increase in combined size by ~150% towards 2050 and then decrease to 2100 as glaciers disconnect from them. Overall, our findings should inform (i) glacier evolution models into which ice-marginal lake effects need incorporating, (ii) studies of rapid landscape evolution and especially of meltwater and sediment delivery, and (iii) considerations of future meltwater supply and water quality.