Hostname: page-component-848d4c4894-m9kch Total loading time: 0 Render date: 2024-05-23T09:49:44.460Z Has data issue: false hasContentIssue false

Topographic modulation of outlet glaciers in Greenland: a review

Published online by Cambridge University Press:  14 August 2023

Ginny Catania*
Institute for Geophysics, University of Texas at Austin, Austin, TX 78712, USA Department of Geological Sciences, University of Texas at Austin Austin, TX 78712, USA
Denis Felikson
Cryospheric Sciences Laboratory, NASA Goddard Space Flight Center Greenbelt, MD 20771, USA Goddard Earth Sciences Technology and Research II, Morgan State University Baltimore, MD 21251, USA
Corresponding author: Ginny Catania; Email:
Rights & Permissions [Opens in a new window]


Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change.

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 (, which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of The International Glaciological Society

1. Introduction

The Greenland Ice Sheet (GrIS) is presently the largest contributor of added mass to rising sea levels (10.9 mm; 22.3%), recently surpassing alpine glaciers and ice caps (7.5 mm; 15.4%), the Antarctic Ice Sheet (6.4 mm; 13.1%) and land water storage (7.2 mm; 14.8%) (IPCC, 2021). For comparison, sea level rise due to thermal expansion accounts for 16.7 mm (34.4%) (IPCC, 2021). For Greenland, most of the mass loss is concentrated around the periphery of the ice sheet (Velicogna and others, Reference Velicogna2020) due to the low elevations found in these regions (Pritchard and others, Reference Pritchard, Arthern, Vaughan and Edwards2009; Shepherd and others, Reference Shepherd2020), which promotes large negative surface mass balance (Fettweis and others, Reference Fettweis2020), abundant supraglacial melt that promotes fast flow at the ice sheet bed (Andrews and others, Reference Andrews2014; Nienow and others, Reference Nienow, Sole, Slater and Cowton2017), and the presence of nearly 300 fast-flowing outlet glaciers that act as conveyor belts draining ice from the interior toward the coast (Joughin and others, Reference Joughin, Smith, Howat, Scambos and Moon2010; Moon and others, Reference Moon, Gardner, Csatho, Parmuzin and Fahnestock2020). Outlet glaciers are impacted by numerous processes including external climate triggers (surface mass balance and oceanographic heat transfer), which in turn impact the boundary conditions of glaciers (Nick and others, Reference Nick, Vieli, Howat and Joughin2009; Howat and others, Reference Howat, Box, Ahn, Herrington and McFadden2010; Straneo and Heimbach, Reference Straneo and Heimbach2013; Catania and others, Reference Catania, Stearns, Moon, Enderlin and Jackson2020). Externally, topography of the ice sheet bed influences surface elevation and ice thickness gradients and thus exerts a first-order control on ice dynamics. This means that while outlet glaciers remain acutely sensitive to climate change (ocean and atmosphere), the magnitude of dynamic adjustment that they make in response to climate perturbations is to a large degree controlled by topography. This idea has been suggested previously by Pfeffer (Reference Pfeffer2007) and forms the premise of this brief review paper.

Our review of the current understanding of how topography controls GrIS mass loss focuses on outlet glacier dynamic response, which exhibits heterogeneity between individual glacier catchments in ice thickness (Csatho and others, Reference Csatho2014), velocity (Moon and others, Reference Moon, Gardner, Csatho, Parmuzin and Fahnestock2020), and terminus change (Murray and others, Reference Murray2015; Catania and others, Reference Catania2018). We first discuss early theories underlying our understanding of irreversible glacier retreat related to topography and how modern observations and modeling of subglacial topography may be used to explain the observed heterogeneous response of the GrIS outlet glaciers to recent climate change. We then discuss the role of topography in modulating the ability, timing and amount of terminus retreat of marine-terminating outlet glaciers as well as how topography far inland from the terminus influences the amount of mass loss resulting from a terminus perturbation. Finally, we will explore the commonly held assumption that topographic change with time is slow and not necessary to include in ice sheet models. We end with an outlook on future research priorities aimed at improving our understanding of topographic influence on the GrIS.

2. How topography controls ice sheet change

The bed topography of ice sheets is deeper in the interior than at the periphery of the ice sheet owing to lithospheric loading that depresses the crust where the overlying weight of ice is greatest (Fig. 1). As a result, Earth's remaining ice sheets both have basal topography that is largely characterized as retrograde - sloped inward to the ice sheet interior from the margins. Early research into the impact of ice sheets on retrograde bed slopes under climate warming led to the development of a theory now termed Marine Ice Sheet Instability (MISI) (e.g., Weertman, Reference Weertman1974; Mercer, Reference Mercer1978; Thomas and Bentley, Reference Thomas and Bentley1978). MISI suggests that glaciers and ice sheets resting on retrograde beds are more susceptible to rapid and irreversible mass loss as the climate warms. This is because the flux of ice through the grounding line increases with bed depth. Thus, ice sheet retreat on retrograde bed slopes results in greater flux of ice and additional retreat. This runaway process is suggested as a possible mechanism that could lead to rapid collapse of the West Antarctic Ice Sheet, with the potential for much faster rates of sea level rise than currently estimated (IPCC, 2021). Early theories regarding MISI laid the groundwork for modeling this process. This includes Schoof (Reference Schoof2007), who found that the presence of topographic overdeepenings results in hysteresis in terminus position; termini that sit on retrograde bed slopes have unstable grounding line states and can rapidly retreat given a climate perturbation, while termini that sit on prograde bed slopes are stable and can advance given a climate perturbation.

Figure 1. Topography of the Greenland Ice Sheet from BedMachine (Morlighem and others, Reference Morlighem2023) showing the location of knickpoints (black dots) and outlet glacier flow lines (black and red lines). Black lines are glaciers with knickpoints, red do not have knickpoints. All data from Felikson and others (Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020).

3. Topography observations for Greenland

While the theory and modeling of MISI demonstrated its importance, there has been little observational evidence of it in part due to the lack of adequately resolved bed topography. Bed topography became more readily resolved with the collection and publication of radar-derived bed topographic data and the use of these data via the production of BedMachine (Figure 1; Morlighem and others, Reference Morlighem2011), a data-constrained estimate of the bed elevation for the GrIS. BedMachine uses radar-derived ice thickness, surface velocity and surface mass balance data with the assumption of mass continuity to estimate bed elevation for fast-flowing portions of the ice sheet. This approach is an improvement over kriging because it provides a physically-realistic bed elevation in areas where radar data are unavailable. This technique is focused at the fast-flowing margin of the ice sheet and in the slower-flowing interior, BedMachine relies on kriging to interpolate between radar-derived ice thickness measurements to create a complete map of topography beneath the ice sheet. BedMachine has been improved over time with the addition of fjord bathymetry data (Fenty and others, Reference Fenty2020), providing seamless topography across outlet glacier termini (Morlighem and others, Reference Morlighem2023).

While BedMachine provides improvements to bed topographic estimates, the terminal regions of outlet glaciers remain poorly resolved because radar-derived ice thickness is more difficult to obtain in these wet, crevassed and deep regions. Further complicating improvements to outlet glacier topographic measurements is the fact that there are nearly 300 outlet glaciers, each with their own unique characteristics, making data collection onerous. The topography that shapes outlet glaciers partially arises because glaciers are natural erosive agents (Kessler and others, Reference Kessler, Anderson and Briner2008; Koppes and Montgomery, Reference Koppes and Montgomery2009; Love and others, Reference Love, Hallet, Pratt and O'Neel2016) capable of delivering large amounts of sediment to their termini where it often becomes visible as sediment plumes at the fjord surface (McGrath and others, Reference McGrath2010; Hudson and others, Reference Hudson2014). Glacier erosion is likely to be largest where ice is moving fastest and where there is efficient removal of sediment via an active subglacial system (Hallet and others, Reference Hallet, Hunter and Bogen1996; Cowton and others, Reference Cowton, Nienow, Bartholomew, Sole and Mair2012; Love and others, Reference Love, Hallet, Pratt and O'Neel2016; Brinkerhoff and others, Reference Brinkerhoff, Truffer and Aschwanden2017). Supraglacial water delivery to the bed downstream of the equilibrium line provides a steady water supply in summers, possibly explaining the much deeper beds observed in these locations (Fig. 2). In addition to erosion, because many outlet glaciers approach floatation toward their margins, sediment may preferentially deposit close to the terminus creating a moraine (sometimes called a sill in the literature), which also shapes the topography at termini (Figure 2; Batchelor and others, Reference Batchelor, Dowdeswell and Rignot2018). Comparisons of BedMachine to radar-derived topography in outlet glaciers reveals that there may be a systematic offset of the bed elevation between the two but that the radar-derived bed slopes are well-preserved in BedMachine (Morlighem and others, Reference Morlighem2017; Catania and others, Reference Catania2018; Narkevic and Anton, Reference Narkevic, Csatho and Anton2023). In addition, BedMachine is unrealistically smooth at small scales making it less useful for examining subglacial water routing and the influence of small-scale bed features on dynamics (MacKie and others, Reference MacKie, Schroeder, Zuo, Yin and Caers2021). Stochastic modeling of bed topography has emerged as an additional technique to simulate the most realistic bed topography using the statistics of the bed elevation uncertainty to drive the scope of simulations (e.g., Goff and Jordan, Reference Goff and Jordan1988; MacKie and others, Reference MacKie, Schroeder, Caers, Siegfried and Scheidt2020, Reference MacKie, Schroeder, Zuo, Yin and Caers2021).

Figure 2. Schematic of outlet glacier topography showing retrograde inland topography, the presence of a knickpoint at the location of the equilibrium line altitude (ELA), an overdeepening in the region of fast flow and the presence of moraines (both active and paleo) at the glacier terminus.

3.1 Topographic controls on retreat and inland thinning

Much of the dynamic mass loss for the GrIS is thought to have initiated by warm ocean temperatures and sufficient supraglacial melt to drive enhanced melt of the ice marginal region. This has been demonstrated as a consequence of 20th century ocean and atmospheric warming (Holland and others, Reference Holland, Thomas, Young, Ribergaard and Lyberth2008; Murray and others, Reference Murray2010; Straneo and Heimbach, Reference Straneo and Heimbach2013; Trusel and others, Reference Trusel2018; Wood and others, Reference Wood2021). However, while the overwhelming majority of GrIS outlet glaciers have experienced retreat (King and others, Reference King2020), acceleration (Moon and others, Reference Moon, Joughin, Smith and Howat2012) and thinning (Smith and others, Reference Smith2020), there is heterogeneity in the dynamic response of glaciers to climate (Csatho and others, Reference Csatho2014; Catania and others, Reference Catania2018; Hill and others, Reference Hill, Carr, Stokes and Gudmundsson2018; Moon and others, Reference Moon, Gardner, Csatho, Parmuzin and Fahnestock2020) with a small number of glaciers exhibiting stability over this time period. In part, this may be due to the inability for warm ocean waters to access the terminus where fjords are shallow or protected with shallow sills or paleo-moraines downstream of glacier termini (Bartholomaus and others, Reference Bartholomaus, Larsen and O'Neel2013; Carroll and others, Reference Carroll2016; Batchelor and others, Reference Batchelor, Dowdeswell, Rignot and Millan2019). However, the climate alone is not sufficient to force retreat of every glacier terminus. This was demonstrated by Christian and others (Reference Christian, Robel and Catania2022) who found that some bed bumps at glacier termini could cause persistent terminus stability even when climate change produces an increased probability of retreat. This is because ice flux reduces as ice flows up to a bed peak, flattening surface slopes and reducing the ability for climate to affect terminus position (Robel and others, Reference Robel, Pegler, Catania, Felikson and Simkins2022). The enhanced stability of glacier termini near the peaks of bed bumps may thus explain why there are a number of persistently stable glaciers despite widespread coincident warming of the ocean and atmosphere around Greenland since the late 1990s.

In addition to controlling the timing of glacier terminus retreat since the 1990s, topography also regulates how much retreat occurs. Using BedMachine data, Catania and others (Reference Catania2018) confirmed that glacier termini retreat from one bed bump to another further inland with retreat that appears insensitive to bed bumps smaller than the seasonal amplitude of the terminus. Carnahan and others (Reference Carnahan, Catania and Bartholomaus2022) examined this for neighboring glaciers Umiamiko Isbræ and Ingia Isbræ, both of which began to retreat in ~2001. While both glaciers retreated significantly, Umiamiko Isbræ restabilized on the prograde side of a large bed bump in 2010, while Ingia Isbræ has continued to retreat (Zhang and others, Reference Zhang, Catania and Trugman2023). The ongoing retreat of Ingia Isbræ was enabled, in part, because its fjord has flat bed topography, which permits low basal resistance to driving stress for several kilometers upstream of the terminus (Carnahan and others, Reference Carnahan, Catania and Bartholomaus2022). The ongoing retreat occurring for many GrIS glaciers in spite of widespread ocean cooling in ~2008 (Wood and others, Reference Wood2021) suggests that climate alone cannot sustain retreat, and that topography may permit the degree to which a climate trigger will influence future dynamics of individual glaciers.

While we have focused largely on bed topographic controls so far, we note that other variables can exert control on the pace and timing of retreat. In streaming ice, Greenwood and others (Reference Greenwood, Simkins, Winsborrow and Bjarnadóttir2021) found that along-flow bed slope was a poor predictor of retreat style with similar retreats occurring on all types of bed slopes. This suggests that bed topography alone may not be the dominant control to terminus retreat. For the relatively narrow outlet glaciers in Greenland, changes in fjord width also appear to exert control in model studies of glacier retreat (Enderlin and others, Reference Enderlin, Howat and Vieli2013; Akesson and others, Reference Akesson, Nisancioglu and Nick2018; Hill and others, Reference Hill, Carr, Stokes and Gudmundsson2018). Indeed, model simulations from Akesson and others (Reference Akesson, Nisancioglu and Nick2018) show that termini can retreat through fjord embayments even in the presence of bed topographic bumps in the bed. For the GrIS, it has been difficult to ascertain the degree to which changes in fjord width have impacted observed retreat rates because observed variations in width are small compared to variations in bed topography over the relatively short length scale of most retreats (Catania and others, Reference Catania2018).

Terminus retreat is thought to precede surface steepening, acceleration and inland thinning of outlet glaciers (Carnahan and others, Reference Carnahan, Catania and Bartholomaus2022). Inland thinning of ice is diffusive and leads to slow, long-term mass loss but represents the majority of future committed sea level rise (Price and others, Reference Price, Payne and Howat2011). The amount of inland thinning permitted is also controlled by bed topography (Felikson and others, Reference Felikson2017), with the presence of bed topographic ‘knickpoints’ - steep reaches where the bed rises from below sea level to above sea level far inland from the glacier terminus (Figures 1, 3; Felikson and others, Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020). These knickpoints also control the degree to which interior ice accelerates in response to thinning and retreat (Williams and others, Reference Williams, Gourmelen and Nienow2021). The proximity of knickpoints to a glacier's terminus and the steepness of that knickpoint is correlated with regional topographic steepness, which likely steers outlet glacier tributaries to converge enhancing bed erosion (Kessler and others, Reference Kessler, Anderson and Briner2008; Felikson and others, Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020). In regions of steeper terrain, there are steeper, more well-defined knickpoints, while in gentler terrain, there are no knickpoints or much gentler sloped knickpoints (Figures 1, 3; Felikson and others, Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020). Where they are present, knickpoint locations are roughly coincident with the location of the equilibrium line, suggesting that surface meltwater is required to sufficiently lubricate the bed enabling enhanced erosion downstream of knickpoints as opposed to frozen bed conditions upstream of knickpoints (Felikson and others, Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020). When retreat occurs on glaciers without well-defined knickpoints, inland diffusive thinning may be slower, but occurs for much longer than for glaciers with well-defined knickpoints. This suggests that glaciers without knickpoints may represent bottlenecks of mass loss; where ongoing long-term mass loss can continue well after the retreat has occurred.

Figure 3. Subglacial and surface topography along flowlines of two GrIS glaciers. Grey lines show smoothed bed topography along six individual flowlines for each glacier, black line shows the mean of all six flowlines. Blue shows ice surface topography. All topography data from (Felikson and others, Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020). Sea level is indicated at zero elevation with a red dotted line. The approximate equilibrium line elevation is ~1500 m (Noël and others, Reference Noël, Berg, Lhermitte and Broeke2019) and is indicated for each glacier. (a) Humbolt Glacier showing an overdeepened bed near the terminus but no presence of a strong knickpoint detected. (b) Helheim Glacier showing a strongly overdeepened bed topography near the terminus and a steep knickpoint at ~35 km where inland thinning would be limited according to Felikson and others (Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020).

Given the importance of bed topography, and our knowledge that glaciers are effective at erosion, it seems logical then to suspect that topography can co-evolve with glacier dynamics. Indeed, this was proposed to explain the ‘tidewater glacier cycle’ (Meier and Post, Reference Meier and Post1987; Nick and others, Reference Nick, Veen and Oerlemans2007), which is a long-term cycle of slow terminus advance and rapid retreat initially described as typical of Alaskan tidewater glaciers. The importance of coupling sediment to ice dynamics was demonstrated by Brinkerhoff and others (Reference Brinkerhoff, Truffer and Aschwanden2017) who showed that the tidewater glacier cycle could be reproduced within a steady climate simply through interactions between ice flow, glacier erosion and sediment transport. Glacier erosion is also responsible for the creation of overdeepenings (Patton and others, Reference Patton, Swift, Clark, Livingstone and Cook2016), which then feedback onto the ice dynamics including the rate of terminus retreat (Robel and others, Reference Robel, Pegler, Catania, Felikson and Simkins2022) and the ice flux (Hooke, Reference Hooke1991; Creyts and others, Reference Creyts, Clarke and Church2013). Differences in glacier bed erosion rates are also likely responsible for the formation of knickpoints (Kessler and others, Reference Kessler, Anderson and Briner2008; Felikson and others, Reference Felikson, Catania, Bartholomaus, Morlighem and Noël2020), which suggests that long-term erosion may be responsible for the heterogeneous dynamic response in inland thinning of the ice sheet that is observed today.

4. Future research priorities

While the role of topography in controlling outlet glacier dynamics is of clear importance, our ability to actually observe subglacial topography with reasonable accuracy is quite recent. For the GrIS, considerable data acquisition occurred through CReSIS (Gogineni and others, Reference Gogineni2001, Reference Gogineni2014), NASA's Operation IceBridge (Studinger and others, Reference Studinger, Koenig, Martin and Sonntag2010; MacGregor and others, Reference MacGregor2021) and Oceans, Melting, Greenland Missions (Fenty and others, Reference Fenty2020), which provided significant improvements in the spatial resolution and coverage of bed topography in Greenland. Despite the significant funding and effort that went to securing these data, we still lack adequate topographic data for many outlet glacier terminal regions in Greenland. Thick, warm, wet and steep-walled ice conditions here pose unique challenges for airborne radar data collection. Some progress is being made to counter these challenges using unmanned aircraft, which can house lower frequency radar systems than are typically used, permitting deeper penetration (Arnold and others, Reference Arnold2018). Additional dedicated funding is needed to fully map these parts of the ice sheet, perhaps focusing on those glaciers that are most susceptible to initiating large changes in mass loss.

Similarly, additional surveying of fjord topography is needed (Jakobsson and Mayer, Reference Jakobsson and Mayer2022), particularly for fjords that are persistently chocked with melange, making the terminus region much more difficult to access via ship (e.g., Helheim Glacier). For glaciers without persistent melange, coordination across nations can crowd-source data collection to happen during periods of time when accessibility is available. Perhaps this means that we make use of even single-beam sounding from fishing and expedition ships working in Greenland. For fjords that have persistent melange, a different approach is needed that allows remote sounding of the sea floor topography. This could occur via including remotely operated vehicles (Jakobsson and Mayer, Reference Jakobsson and Mayer2022), using novel remote-sensing techniques that make use of iceberg draft heights (Scheick and others, Reference Scheick, Enderlin, Miller and Hamilton2019) and the use of airborne gravity (Boghosian and others, Reference Boghosian2015; Tinto and others, Reference Tinto, Bell, Cochran and Münchow2015). While airborne gravity provides a coarser resolution bed topography estimate compared to multibeam data, it can be done more easily by plane and with complete coverage of the terminal zone (An and others, Reference An, Rignot, Millan, Tinto and Willis2019).

Given the importance of bed topography, particularly at the ice-ocean boundary, we must next understand the pace at which topographic change is possible. Glaciologists largely assume stable bed topography over time assuming that bed erosion and deposition rates are small compared to terminus change rates (Koppes and Montgomery, Reference Koppes and Montgomery2009). Yet, sedimentation has been observed to be critical to understanding the stability state of glacier termini (Alley and others, Reference Alley, Lawson, Larson, Evenson and Baker2003, Reference Alley, Anandakrishnan, Dupont, Parizek and Pollard2007). Indeed, glacier advance has been shown to be uniquely dependent on sediments infilling into bed lows in front of the advancing terminus (Nick and others, Reference Nick, Veen and Oerlemans2007). Further, sedimentation rates at Alaskan glaciers were recently measured to be on the order of several meters per year (Eidam and others, Reference Eidam2020), which is larger than the vertical motion of the solid Earth following deglaciation, now considered an important stabilizing feedback on ice loss (Barletta and others, Reference Barletta2018). Despite its clear importance, sedimentation is either entirely missing from sea level projecting models of ice sheets (Aschwanden and others, Reference Aschwanden2019) or is modeled without clear knowledge of the types and rates of processes that contribute to moraine building (Brinkerhoff and others, Reference Brinkerhoff, Truffer and Aschwanden2017). While there has been extensive research on glacial landforms from past glaciations, most of these studies are not able to produce a precise depiction of the coincident ice dynamics at the time of deposition. Thus, we lack a set of governing equations that describe how to couple ice and sediment dynamics in a way that is consistent with observations. We thus recommend new observations of sedimentation rates that can be paired with observed glacier dynamics so that we can build equations that describe moraine-building and erosion of overdeepenings that are consistent with observations.

Finally, there is a remaining need to consider outlet glacier dynamics holistically because there are multiple processes (both internal and external to the ice sheet) that impact dynamics making it difficult to tease apart cause and effect of glacier change. Poorly constrained boundary conditions for outlet glaciers exacerbates this (Malles and others, Reference Malles2023). New observations must therefore be coupled to focused modeling of ice sheet outlet glaciers to discern the processes that are most important for accurately estimating future sea level.

5. Conclusions

Ice sheet mass loss has direct implications to sea level rise for coastal communities who rely on accurate forecasts of sea level across a wide range of time scales (Larour and others, Reference Larour, Ivins and Adhikari2017; Ultee and others). To address this need, there have been increased efforts to coordinate and improve model predictions of ice sheets over the last decade (Nowicki and others, Reference Nowicki2016; Seroussi and others, Reference Seroussi2020). Capturing historical GrIS mass change in model simulations remains a challenge (Aschwanden and others, Reference Aschwanden, Bartholomaus, Brinkerhoff and Truffer2021), which is due to a range of uncertainties including the lack of understanding of processes that control ice sheet mass loss. Such uncertainties make accurate model prediction of sea level challenging. For example, recent modeled future mass loss of the GrIS suggests that it will contribute somewhere between 5–33 cm to sea level by 2100 (Aschwanden and others, Reference Aschwanden2019). Meanwhile, the most recent IPCC report for the first time included a ‘low-likelihood, high-impact storyline’ suggesting that sea level could be 0.5 m or more higher than anticipated by 2100 from ‘deeply uncertain processes related to ice sheet instability’ (IPCC, 2021). Such a large range in future sea level means the difference between a coastal city that remains largely untouched by sea level rise versus one that becomes submerged. The research community has the responsibility to improve ice sheet uncertainties. Within the focus of this review, we argue for improved radar-data coverage over the more difficult to access terminal regions of outlet glaciers that will improve the mass conserving bed solution. We also argue for observations of the role of sedimentation/erosion for sculpting and changing bed topography over time. To a first order we need a better understanding of the rates of topographic change that are possible and what controls such rates.


Funding for the work leading to these ideas comes from NASA (NNX12AP50G and 80NSSC18K1477).


Akesson, H, Nisancioglu, KH and Nick, FM (2018) Impact of fjord geometry on grounding line stability. Frontiers in Earth Science 6, 119216. doi:10.3389/feart.2018.00071Google Scholar
Alley, RB, Lawson, DE, Larson, GJ, Evenson, EB and Baker, GS (2003) Stabilizing feedbacks in glacier-bed erosion. Nature 424(6950), 758760. doi:10.1038/nature01839Google Scholar
Alley, RB, Anandakrishnan, S, Dupont, TK, Parizek, BR and Pollard, D (2007) Effect of sedimentation on ice-sheet grounding-line stability. Science 315(5), 18381841. doi:10.1126/science.1138396Google Scholar
An, L, Rignot, E, Millan, R, Tinto, K and Willis, J (2019) Bathymetry of northwest Greenland Using ‘Ocean Melting Greenland’ (OMG) high-resolution airborne gravity and other data. Remote Sensing 11(2), 131. doi:10.3390/rs11020131Google Scholar
Andrews, LC and 7 others (2014) Direct observations of evolving subglacial drainage beneath the Greenland Ice Sheet. Nature 514(7520), 8083. doi:10.1038/nature13796Google Scholar
Arnold, E and 13 others (2018) HF/VHF radar sounding of ice from manned and unmanned airborne platforms. Geosciences 8(5), 182. doi:10.3390/geosciences8050182Google Scholar
Aschwanden, A and 7 others (2019) Contribution of the Greenland Ice Sheet to sea level over the next millennium. Science Advances 5(6), eaav9396. doi:10.1126/sciadv.aav9396Google Scholar
Aschwanden, A, Bartholomaus, TC, Brinkerhoff, DJ and Truffer, M (2021) Brief communication: a roadmap towards credible projections of ice sheet contribution to sea-level. The Cryosphere Discussions 2021, 114. doi:10.5194/tc-2021-175Google Scholar
Barletta, VR and 16 others (2018) Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability. Science 360(6395), 13351339. doi:10.1126/science.aao1447Google Scholar
Bartholomaus, TC, Larsen, CF and O'Neel, S (2013) Does calving matter? evidence for significant submarine melt. Earth and Planetary Science Letters 380, 2130. doi:10.1016/j.epsl.2013.08.014Google Scholar
Batchelor, C, Dowdeswell, J and Rignot, E (2018) Submarine landforms reveal varying rates and styles of deglaciation in North-West Greenland fjords. Marine Geology 402, 6080. doi:10.1016/j.margeo.2017.08.003Google Scholar
Batchelor, CL, Dowdeswell, JA, Rignot, E and Millan, R (2019) Submarine moraines in southeast Greenland fjords reveal contrasting outlet–glacier behavior since the last glacial maximum. Geophysical Research Letters 46(6), 32793286. doi:10.1029/2019gl082556Google Scholar
Boghosian, A and 6 others (2015) Resolving bathymetry from airborne gravity along Greenland fjords. Journal of Geophysical Research Solid Earth 120(12), 85168533. doi:10.1002/2015jb012129Google Scholar
Brinkerhoff, D, Truffer, M and Aschwanden, A (2017) Sediment transport drives tidewater glacier periodicity. Nature Communications 8(1), 91. doi:10.1038/s41467-017-00095-5Google Scholar
Carnahan, E, Catania, G and Bartholomaus, TC (2022) Observed mechanism for sustained glacier retreat and acceleration in response to ocean warming around Greenland. The Cryosphere 16(10), 43054317. doi:10.5194/tc-16-4305-2022Google Scholar
Carroll, D and 11 others (2016) The impact of glacier geometry on meltwater plume structure and submarine melt in Greenland fjords. Geophysical Research Letters 43(18), 97399748. doi:10.1002/2016gl070170Google Scholar
Catania, GA and 7 others (2018) Geometric controls on tidewater glacier retreat in central western Greenland. Journal of Geophysical Research: Earth Surface 123, 20242038. doi:10.1029/2017jf004499Google Scholar
Catania, GA, Stearns, LA, Moon, T, Enderlin, EM and Jackson, RH (2020) Future evolution of Greenland's marine terminating outlet glaciers. Journal of Geophysical Research: Earth Surface 125(2), 128. doi:10.1029/2018jf004873Google Scholar
Christian, JE, Robel, AA and Catania, G (2022) A probabilistic framework for quantifying the role of anthropogenic climate change in marine-terminating glacier retreats. The Cryosphere 16(7), 27252743. doi:10.5194/tc-16-2725-2022Google Scholar
Cowton, TR, Nienow, PW, Bartholomew, ID, Sole, AJ and Mair, DWF (2012) Rapid erosion beneath the Greenland Ice Sheet. Geology 40(4), 343346. doi:10.1130/g32687.1Google Scholar
Creyts, TT, Clarke, GK and Church, M (2013) Evolution of subglacial overdeepenings in response to sediment redistribution and glaciohydraulic supercooling. Journal of Geophysical Research: Earth Surface 118(2), 423446. doi:10.1002/jgrf.20033Google Scholar
Csatho, BM and 9 others (2014) Laser altimetry reveals complex pattern of Greenland Ice Sheet dynamics. Proceedings of the National Academy of Sciences 111(52), 1847818483. doi:10.1073/pnas.1411680112Google Scholar
Eidam, EF and 5 others (2020) Morainal bank evolution and impact on terminus dynamics during a tidewater glacier stillstand. Journal of Geophysical Research: Earth Surface 125(11), e2019JF005359. doi:10.1029/2019jf005359Google Scholar
Enderlin, EM, Howat, IM and Vieli, A (2013) High sensitivity of tidewater outlet glacier dynamics to shape. The Cryosphere 7(3), 10071015. doi:10.5194/tc-7-1007-2013Google Scholar
Felikson, D and 11 others (2017) Inland thinning on the Greenland Ice Sheet controlled by outlet glacier geometry. Nature Geoscience 10(5), 366369. doi:10.1038/ngeo2934Google Scholar
Felikson, D, Catania, G, Bartholomaus, TC, Morlighem, M and Noël, BPY (2021) Steep glacier bed knickpoints mitigate inland thinning in Greenland. Geophysical Research Letters 48(2), e2020GL090112. doi:10.1029/2020gl090112Google Scholar
Fenty, IG and 16 others (2016) Oceans melting Greenland: early results from NASA's ocean-ice mission in Greenland. Oceanography 29(4), 7283. doi:10.5670/oceanog.2016.100Google Scholar
Fettweis, X and 40 others (2020) GrSMBMIP: intercomparison of the modelled 1980–2012 surface mass balance over the Greenland Ice Sheet. The Cryosphere 14(11), 39353958. doi:10.5194/tc-14-3935-2020Google Scholar
Goff, JA and Jordan, TH (1988) Stochastic modeling of seafloor morphology: inversion of sea beam data for second–order statistics. Journal of Geophysical Research: Solid Earth 93(B11), 1358913608. doi:10.1029/jb093ib11p13589Google Scholar
Gogineni, S and 9 others (2001) Coherent radar ice thickness measurements over the Greenland Ice Sheet. Journal of Geophysical Research: Atmospheres 106(D24), 3376133772. doi:10.1029/2001jd900183Google Scholar
Gogineni, SP and 10 others (2014) Bed topography of Jakobshavn Isbræ, Greenland, and Byrd Glacier, Antarctica. Journal of Glaciology 60(223), 813833. doi:10.3189/2014jog14j129Google Scholar
Greenwood, SL, Simkins, LM, Winsborrow, MCM and Bjarnadóttir, LR (2021) Exceptions to bed-controlled ice sheet flow and retreat from glaciated continental margins worldwide. Science Advances 7(3), eabb6291. doi:10.1126/sciadv.abb6291Google Scholar
Hallet, B, Hunter, L and Bogen, J (1996) Rates of erosion and sediment evacuation by glaciers: a review of field data and their implications. Global and Planetary Change 12(1-4), 213235. doi:10.1016/0921-8181(95)00021-6Google Scholar
Hill, EA, Carr, JR, Stokes, CR and Gudmundsson, GH (2018) Dynamic changes in outlet glaciers in northern Greenland from 1948 to 2015. The Cryosphere 12(10), 32433263. doi:10.5194/tc-12-3243-2018Google Scholar
Holland, DM, Thomas, RH, Young, BD, Ribergaard, MH and Lyberth, B (2008) Acceleration of Jakobshavn Isbræ triggered by warm subsurface ocean waters. Nature Geoscience 1, 659664. doi:10.1038/ngeo316Google Scholar
Hooke, RL (1991) Positive feedbacks associated with erosion of glacial cirques and overdeepenings. Bulletin of the Geological Society of America 103(8), 11041108.Google Scholar
Howat, IM, Box, JE, Ahn, Y, Herrington, A and McFadden, EM (2010) Seasonal variability in the dynamics of marine-terminating outlet glaciers in Greenland. Journal of Glaciology 56(198), 601613.Google Scholar
Hudson, B and 5 others (2014) MODIS observed increase in duration and spatial extent of sediment plumes in Greenland fjords. The Cryosphere 8(4), 11611176. doi:10.5194/tc-8-1161-2014Google Scholar
IPCC (2021) Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Technical report, Cambridge University Press.Google Scholar
Jakobsson, M and Mayer, LA (2022) Polar region bathymetry: critical knowledge for the prediction of global sea level rise. Frontiers in Marine Science 8, 788724. doi:10.3389/fmars.2021.788724Google Scholar
Joughin, IR, Smith, BE, Howat, IM, Scambos, TA and Moon, T (2010) Greenland flow variability from ice-sheet-wide velocity mapping. Journal of Glaciology 56(197), 415430.Google Scholar
Kessler, MA, Anderson, RS and Briner, JP (2008) Fjord insertion into continental margins driven by topographic steering of ice. Nature Geoscience 1(6), 365369. doi:10.1038/ngeo201Google Scholar
King, MD and 8 others (2020) Dynamic ice loss from the Greenland Ice Sheet driven by sustained glacier retreat. Nature Communications Earth and Environment 1(1), 1. doi:10.1038/s43247-020-0001-2Google Scholar
Koppes, MN and Montgomery, DR (2009) The relative efficacy of fluvial and glacial erosion over modern to orogenic timescales. Nature Geoscience 2(9), 644647. doi:10.1038/ngoe616Google Scholar
Larour, E, Ivins, ER and Adhikari, S (2017) Should coastal planners have concern over where land ice is melting?. Science Advances 3(11), e1700537. doi:10.1126/sciadv.1700537Google Scholar
Love, KB, Hallet, B, Pratt, TL and O'Neel, S (2016) Observations and modeling of fjord sedimentation during the 30 year retreat of Columbia Glacier, AK. Journal of Glaciology 62(234), 778793. doi:10.1017/jog.2016.67Google Scholar
MacGregor, JA and 45 others (2021) The scientific legacy of NASA's operation icebridge. Reviews of Geophysics 59(2), e2020RG000712. doi:10.1029/2020rg000712Google Scholar
MacKie, EJ, Schroeder, DM, Caers, J, Siegfried, MR and Scheidt, C (2020) Antarctic topographic realizations and geostatistical modeling used to map subglacial lakes. Journal of Geophysical Research: Earth Surface 125(3), e2019JF005420. doi:10.1029/2019jf005420Google Scholar
MacKie, EJ, Schroeder, DM, Zuo, C, Yin, Z and Caers, J (2021) Stochastic modeling of subglacial topography exposes uncertainty in water routing at Jakobshavn Glacier. Journal of Glaciology 67(261), 7583. doi:10.1017/jog.2020.84Google Scholar
Malles, JH and 5 others (2023) Exploring the impact of a frontal ablation parameterization on projected 21st-century mass change for Northern Hemisphere glaciers. Journal of Glaciology 116. doi:10.1017/jog.2023.19Google Scholar
McGrath, D and 5 others (2010) Sediment plumes as a proxy for local ice-sheet runoff in Kangerlussuaq Fjord, West Greenland. Journal of Glaciology 56(199), 813821.Google Scholar
Meier, MF and Post, A (1987) Fast tidewater glaciers. Journal of Geophysical Research 92(B9), 90519058.Google Scholar
Mercer, JH (1978) West Antarctic ice sheet and CO2 greenhouse effect: a threat of disaster. Nature 271(5643), 321325. doi:10.1038/271321a0Google Scholar
Moon, TA, Gardner, AS, Csatho, B, Parmuzin, I and Fahnestock, MA (2020) Rapid reconfiguration of the Greenland Ice Sheet coastal margin. Journal of Geophysical Research: Earth Surface 125(11), e2020JF005585. doi:10.1029/2020jf005585Google Scholar
Moon, T, Joughin, IR, Smith, BE and Howat, IM (2012) 21st-century evolution of Greenland outlet glacier velocities. Science 336(6081), 576578. doi:10.1126/science.1219985Google Scholar
Morlighem, M and 5 others (2011) A mass conservation approach for mapping glacier ice thickness. Geophysical Research Letters 38(19), 000. doi:10.1029/2011gl048659Google Scholar
Morlighem, M and 31 others (2017) BedMachine v3: complete bed topography and ocean bathymetry mapping of Greenland from multibeam echo sounding combined with mass conservation. Geophysical Research Letters 44, 1105111061. doi:10.1002/2017gl074954Google Scholar
Morlighem, WC and 30 others (2023) IceBridge BedMachine Greenland, Version 5.Google Scholar
Murray, T and 10 others (2010) Ocean regulation hypothesis for glacier dynamics in southeast Greenland and implications for ice sheet mass changes. Journal of Geophysical Research 115, F03026. doi:10.1029/2009jf001522Google Scholar
Murray, T and 14 others (2015) Extensive retreat of Greenland tidewater glaciers, 2000–2010. Arctic, Antarctic, and Alpine Research 47(3), 427447. doi:10.1657/aaar0014-049Google Scholar
Narkevic, A, Csatho, B and Anton, S (2023) Rapid basal channel growth beneath Greenland's longest floating ice shelf. ESS Open Archive, Feb 20.Google Scholar
Nick, FM, Veen, CJvd and Oerlemans, J (2007) Controls on advance of tidewater glaciers: results from numerical modeling applied to Columbia Glacier. Journal of Geophysical Research 112(F3), 33729. doi:10.1029/2006jf000551Google Scholar
Nick, FM, Vieli, A, Howat, IM and Joughin, IR (2009) Large-scale changes in Greenland outlet glacier dynamics triggered at the terminus. Nature 2(2), 110114. doi:10.1038/ngeo394Google Scholar
Nienow, PW, Sole, AJ, Slater, DA and Cowton, TR (2017) Recent advances in our understanding of the role of meltwater in the Greenland Ice Sheet system. Current Climate Change Reports 3(4), 115. doi:10.1007/s40641-017-0083-9Google Scholar
Noël, B, Berg, WJvd, Lhermitte, S and Broeke, MRvd (2019) Rapid ablation zone expansion amplifies north Greenland mass loss. Science Advances 5(9), eaaw0123. doi:10.1126/sciadv.aaw0123Google Scholar
Nowicki, SMJ and 8 others (2016) Ice sheet model intercomparison project (ISMIP6) contribution to CMIP6. Geoscientific Model Development 9(12), 45214545. doi:10.5194/gmd-9-4521-2016Google Scholar
Patton, H, Swift, DA, Clark, CD, Livingstone, SJ and Cook, SJ (2016) Distribution and characteristics of overdeepenings beneath the Greenland and Antarctic ice sheets: Implications for overdeepening origin and evolution. Quaternary Science Reviews 148(C), 128145. doi:10.1016/j.quascirev.2016.07.012Google Scholar
Pfeffer, WT (2007) A simple mechanism for irreversible tidewater glacier retreat. Journal of Geophysical Research 112(F03S25), 119. doi:10.1029/2006jf000590Google Scholar
Price, SF, Payne, A and Howat, IM (2011) Committed sea-level rise for the next century from {G}reenland ice sheet dynamics during the past decade. Proceedings of the National Academy of Sciences 108(22), 89788983.Google Scholar
Pritchard, HD, Arthern, RJ, Vaughan, DG and Edwards, LA (2009) Extensive dynamic thinning on the margins of the Greenland and Antarctic ice sheets. Nature 461(7266), 971975. doi:10.1038/nature08471Google Scholar
Robel, AA, Pegler, SS, Catania, G, Felikson, D and Simkins, LM (2022) Ambiguous stability of glaciers at bed peaks. Journal of Glaciology 68(272), 11771184. doi:10.1017/jog.2022.31Google Scholar
Scheick, J, Enderlin, EM, Miller, EK and Hamilton, GS (2019) First-order estimates of coastal bathymetry in Ilulissat and Naajarsuit fjords, Greenland, from remotely sensed iceberg observations. Remote Sensing 11(8), 935–19. doi:10.3390/rs11080935Google Scholar
Schoof, CG (2007) Ice sheet grounding line dynamics: Steady states, stability, and hysteresis. Journal of Geophysical Research 112(F3), 1720. doi:10.1029/2006jf000664Google Scholar
Seroussi, H and 46 others (2020) ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century. The Cryosphere 14(9), 30333070. doi:10.5194/tc-14-3033-2020Google Scholar
Shepherd, A and 87 others (2019) Mass balance of the Greenland Ice Sheet from 1992 to 2018. Nature 579(7798), 233239. doi:10.1038/s41586-019-1855-2Google Scholar
Smith, B and 14 others (2020) Pervasive ice sheet mass loss reflects competing ocean and atmosphere processes. Science 368(6496), 12391242. doi:10.1126/science.aaz5845Google Scholar
Straneo, F and Heimbach, P (2013) North Atlantic warming and the retreat of Greenland's outlet glaciers. Nature 504(7478), 3643. doi:10.1038/nature12854Google Scholar
Studinger, M, Koenig, LS, Martin, S and Sonntag, JG (2010) Operation Icebridge: Using instrumented aircraft to bridge the observational gap between ICESat and ICESat-2. 2010 IEEE International Geoscience and Remote Sensing Symposium. Honolulu, HI, USA, 19181919. doi:10.1109/IGARSS.2010.5650555Google Scholar
Thomas, RH and Bentley, CR (1978) A model for Holocene retreat of the West Antarctic ice sheet. Quaternary Research 10(2), 150170. doi:10.1016/0033-5894(78)90098-4Google Scholar
Tinto, KJ, Bell, RE, Cochran, JR and Münchow, A (2015) Bathymetry in Petermann fjord from operation icebridge aerogravity. Earth and Planetary Science Letters 422, 5866. doi:10.1016/j.epsl.2015.04.009Google Scholar
Trusel, LD and 8 others (2018) Nonlinear rise in Greenland runoff in response to post-industrial Arctic warming. Nature 564(7734), 118. doi:10.1038/s41586-018-0752-4Google Scholar
Ultee, L, Arnott, JC, Bassis, J and Lemos, MC (2018) From ice sheets to main streets: intermediaries connect climate scientists to coastal adaptation. Earth's Future 6(3), 299304. doi:10.1002/2018ef000827Google Scholar
Velicogna, I and 9 others (2020) Continuity of ice sheet mass loss in Greenland and Antarctica from the GRACE and GRACE follow–on missions. Geophysical Research Letters 47(8), e2020GL087291. doi:10.1029/2020gl087291Google Scholar
Weertman, J (1974) Stability of the junction of an ice sheet and an ice shelf. Journal of Glaciology 13(67), 311.Google Scholar
Williams, JJ, Gourmelen, N and Nienow, P (2021) Complex multi-decadal ice dynamical change inland of marine-terminating glaciers on the Greenland Ice Sheet. Journal of Glaciology 67(265), 833846. doi:10.1017/jog.2021.31Google Scholar
Wood, MH and 16 others (2021) Ocean forcing drives glacier retreat in Greenland. Science Advances 7(1), e2020GL087291. doi:10.1126/sciadv.aba7282Google Scholar
Zhang, E, Catania, G and Trugman, D (2023) AutoTerm: A ‘big data’ repository of Greenland glacier termini delineated using deep learning. The Cryosphere 2022, 134. doi:10.5194/egusphere-2022-1095Google Scholar
Figure 0

Figure 1. Topography of the Greenland Ice Sheet from BedMachine (Morlighem and others, 2023) showing the location of knickpoints (black dots) and outlet glacier flow lines (black and red lines). Black lines are glaciers with knickpoints, red do not have knickpoints. All data from Felikson and others (2020).

Figure 1

Figure 2. Schematic of outlet glacier topography showing retrograde inland topography, the presence of a knickpoint at the location of the equilibrium line altitude (ELA), an overdeepening in the region of fast flow and the presence of moraines (both active and paleo) at the glacier terminus.

Figure 2

Figure 3. Subglacial and surface topography along flowlines of two GrIS glaciers. Grey lines show smoothed bed topography along six individual flowlines for each glacier, black line shows the mean of all six flowlines. Blue shows ice surface topography. All topography data from (Felikson and others, 2020). Sea level is indicated at zero elevation with a red dotted line. The approximate equilibrium line elevation is ~1500 m (Noël and others, 2019) and is indicated for each glacier. (a) Humbolt Glacier showing an overdeepened bed near the terminus but no presence of a strong knickpoint detected. (b) Helheim Glacier showing a strongly overdeepened bed topography near the terminus and a steep knickpoint at ~35 km where inland thinning would be limited according to Felikson and others (2020).