Hostname: page-component-848d4c4894-pftt2 Total loading time: 0 Render date: 2024-05-31T20:28:04.428Z Has data issue: false hasContentIssue false
  • English
  • Français

GHOSTly flute music: drumlins, moats and the bed of Thwaites Glacier

Published online by Cambridge University Press:  23 August 2023

Richard B. Alley Open the ORCID record for Richard B. Alley [Opens in a new window] ,
Nick Holschuh Open the ORCID record for Nick Holschuh [Opens in a new window] ,
Byron Parizek ,
Lucas K. Zoet Open the ORCID record for Lucas K. Zoet [Opens in a new window] ,
Kiya Riverman Open the ORCID record for Kiya Riverman [Opens in a new window] ,
Atsuhiro Muto Open the ORCID record for Atsuhiro Muto [Opens in a new window] ,
Knut Christianson ,
Elisabeth Clyne ,
Sridhar Anandakrishnan  and
Nathan T. Stevens Open the ORCID record for Nathan T. Stevens [Opens in a new window]
Richard B. Alley*
Affiliation:
Department of Geosciences, and Earth and Environmental Systems Institute, Pennsylvania State University, University Park, PA, USA
Nick Holschuh
Affiliation:
Department of Geology, Amherst College, Amherst, MA, USA
Byron Parizek
Affiliation:
Department of Geosciences, and Earth and Environmental Systems Institute, Pennsylvania State University, University Park, PA, USA Mathematics and Geoscience, Pennsylvania State University, Dubois, PA, USA
Lucas K. Zoet
Affiliation:
Department of Geoscience, University of Wisconsin, Madison, WI, USA
Kiya Riverman
Affiliation:
Department of Environmental Studies, University of Portland, Portland, OR, USA
Atsuhiro Muto
Affiliation:
Department of Earth and Environmental Science, Temple University, Philadelphia, PA, USA
Knut Christianson
Affiliation:
Department of Earth and Space Sciences, University of Washington, Seattle, WA, USA
Elisabeth Clyne
Affiliation:
Department of Geosciences, and Earth and Environmental Systems Institute, Pennsylvania State University, University Park, PA, USA
Sridhar Anandakrishnan
Affiliation:
Department of Geosciences, and Earth and Environmental Systems Institute, Pennsylvania State University, University Park, PA, USA
Nathan T. Stevens
Affiliation:
Department of Geosciences, and Earth and Environmental Systems Institute, Pennsylvania State University, University Park, PA, USA
*
Corresponding author: Richard B. Alley; Email: rba6@psu.edu
Rights & Permissions [Opens in a new window]

Abstract

Glacier-bed characteristics that are poorly known and modeled are important in projected sea-level rise from ice-sheet changes under strong warming, especially in the Thwaites Glacier drainage of West Antarctica. Ocean warming may induce ice-shelf thinning or loss, or thinning of ice in estuarine zones, reducing backstress on grounded ice. Models indicate that, in response, more-nearly-plastic beds favor faster ice loss by causing larger flow acceleration, but more-nearly-viscous beds favor localized near-coastal thinning that could speed grounding-zone retreat into interior basins where marine-ice-sheet instability or cliff instability could develop and cause very rapid ice loss. Interpretation of available data indicates that the bed is spatially mosaicked, with both viscous and plastic regions. Flow against bedrock topography removes plastic lubricating tills, exposing bedrock that is eroded on up-glacier sides of obstacles to form moats with exposed bedrock tails extending downglacier adjacent to lee-side soft-till bedforms. Flow against topography also generates high-ice-pressure zones that prevent inflow of lubricating water over distances that scale with the obstacle size. Extending existing observations to sufficiently large regions, and developing models assimilating such data at the appropriate scale, present large, important research challenges that must be met to reliably project future forced sea-level rise.

Keywords

Glacial tillsglacier modellingsubglacial processes
Type
Letter
Information
Annals of Glaciology , Volume 63 , Issue 87-89 , September 2022 , pp. 153 - 157
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - ND
This is an Open Access article, distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives licence (http://creativecommons.org/licenses/by-nc-nd/4.0), which permits non-commercial re-use, distribution, and reproduction in any medium, provided that no alterations are made and the original article is properly cited. The written permission of Cambridge University Press must be obtained prior to any commercial use and/or adaptation of the article.
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of The International Glaciological Society

Introduction

Thwaites Glacier is the most-likely path from modern ice sheets for large, rapid loss of grounded ice leading to sea-level rise (e.g. Scambos and others, Reference Scambos2017; DeConto and others, Reference DeConto2021). At present, ice in the grounding zone is typically ~500 m thick, but retreat into interior basins could increase the grounding-zone thickness several-fold, allowing much faster ice discharge, with further retreat likely through the marine-ice-sheet instability. Thwaites Glacier flows through an outlet that is exceptionally wide among ice streams, but that is relatively narrow compared to the several-fold wider front that could develop during retreat. An outlet both much wider and thicker would have much greater ice discharge leading to sea-level rise. Furthermore, if Thwaites Glacier loses its ice shelf and retreats into interior basins, the evolving calving front would become higher and wider than any now existing on Earth. Highly uniform, fine-grained ice could support tall cliffs (e.g. Hanson and Hooke, Reference Hanson and Hooke2003; Clerc and others, Reference Clerc, Minchew and Behn2019), but such ice does not exist in large volumes in actual ice sheets, where many types of inhomogeneities and damage are likely to limit subaerial cliffs to values not much taller than 100 m (Bassis and Walker, Reference Bassis and Walker2012; Parizek and others, Reference Parizek2019; see review by Alley and others, Reference Alley2023), especially if there is meltwater wedging open surface crevasses (Pollard and others, Reference Pollard, DeConto and Alley2015). Cliff failure could then drive sea-level rise at unprecedented rates.

Large, rapid sea-level rise might be avoided if ice shelves persist or if the now-retreating grounding zone stabilizes before retreating from the modern ‘bottleneck’ with its relatively narrow, shallow outlet on a bed that is not generally deepening inland (Fig. 1). Ice-shelf persistence or loss involves a wide range of important issues (reviewed by, e.g. Scambos and others, Reference Scambos2017 and Alley and others, Reference Alley2023), some of which affect the grounding-zone stress balance and thus the likelihood or rate of grounding-zone retreat. Retreat also depends on many stabilizing and destabilizing factors; forcing by reduction or loss of ice-shelf backstress and melting from warmer waters circulating in grounding-zone estuaries (Horgan and others, Reference Horgan2013) may be especially important.

Figure 1. Key features of the bed of Thwaites Glacier. (a). Bed elevation from BedMachine Antarctica version 2 (Morlighem, Reference Morlighem2020), with elevation indicated by color scale, ice-shelf front (gray), and grounding line (black), with survey location and seismic lines in red. Thwaites discharges primarily between ~109 and 104 W longitude, across a relatively shallow grounding zone. (b). Bed of the up-glacier part of Thwaites Glacier (Holschuh and others, Reference Holschuh, Christianson, Paden, Alley and Anandakrishnan2020), with elevation (left color bar, m; note that scale differs from 1a, as indicated) from airborne swath radar, and acoustic impedance (right color bar, kg m−2 s−1) from seismic surveys (Muto and others, Reference Muto2019a, Reference Muto, Alley, Parizek and Anandakrishnan2019b). Lower acoustic impedance indicates softer bed (deforming tills), and higher acoustic impedance indicates harder bed (bedrock). Holschuh and other (Reference Holschuh, Christianson, Paden, Alley and Anandakrishnan2020) and Alley and others (Reference Alley2021) provide additional, annotated data on the bed character. Radar data and gridded topographies are accessible through the University of Washington ResearchWorks Archive (http://hdl.handle.net/1773/44950).

Basal ‘sliding law’

Even if ice-shelf and grounding-zone processes were well-known and modeled accurately, great uncertainty is attached to the internal and basal ice-flow response. Changes in internal deformation of the ice are important, with significant uncertainties, but the largest issues arise from uncertainties about changing basal motion of the ice.

Basal motion is often represented in models by a ‘sliding law’. Importantly, this is neither a law, nor truly sliding in many cases (e.g. Cuffey and Paterson, Reference Cuffey and Paterson2010, ch. 7). The ‘law’ is a parameterization of important contributions to velocity from processes in debris-laden basal ice, between that ice and bedrock or till, and within till. The ‘law’ must represent ice motion over a wide range of topographic features and a mosaic of bed types involving both till and bedrock. Many of these features are large enough to be resolved in geophysical surveys and models, but smaller features require parameterizations (Hoffman and others, Reference Hoffman2022). We will refer to sliding laws here for simplicity.

How the basal processes are parameterized in ice-flow models fundamentally affects the model projections. For example, in a flowline model of Thwaites Glacier, Parizek and others (Reference Parizek2013) showed that under strong forcing, whether or not the grounding zone retreated into interior basins depended on the sliding law.

Uncertainty in the sliding law impacts most or all projections of ice-sheet and glacier change (e.g. Åkesson and others, Reference Åkesson, Morlighem, O'Regan and Jakobsson2021), and is illustrated clearly for the Antarctic ice sheet including Thwaites Glacier by the results of the ABUMIP model intercomparison (Antarctic BUttressing Model Intercomparison Project; Sun and others, Reference Sun2020). A total of 15 models of various types, with a wide range of basal sliding laws, subjected their modeled Antarctic ice sheet to either very high sub-ice-shelf melt rates or sustained ice-shelf loss, although without the additional cliff calving that can drive the fastest sea-level rise. The forcing was applied at the start of the runs and sustained throughout, and thus notably exceeds realistic forcing, but provides an excellent way to assess the effects of model parameterizations. For the large sub-ice-shelf melt rate, response over the first model century ranged from slight ice-sheet growth to mass loss of almost 7 m of sea-level equivalent. For specified ice-shelf loss, all models showed ice-sheet shrinkage, but the range of results spanned roughly an order of magnitude after a century, again with almost 7 m of sea-level equivalent for the fastest loss. Many differences between the models were important in the disparate results, but strong dependence on the level of plasticity of the sliding law was clear (see their Figure 5 comparing sea-level rise from different models, ordered by sliding law, and Figure 6 showing results from one model with different sliding laws; Sun and others, Reference Sun2020). Additional insights into the importance of the sliding law for ice-sheet response to forcing are found in many other studies, including Brondex and others (Reference Brondex, Gillet-Chaulet and Gagliardini2019), Joughin and others (Reference Joughin, Smith and Schoof2019), Schwans and others (Reference Schwans, Parizek, Alley, Anandakrishnan and Morlighem2023), and Kazmierczak and others (Reference Kazmierczak, Sun, Coulon and Pattyn2022).

Sliding laws are often formulated such that basal velocity increases with basal shear stress raised to some power, which may range from 1 (linear-viscous) to approaching infinite (more-nearly plastic, perhaps with zero motion specified for shear stress below some shear strength), or perhaps some combination such as linear-viscous motion for low shear stresses but plastic or frictional motion for stress above some ‘Iken limit’, the maximum shear stress that the bed can support (see, e.g. Budd and others, Reference Budd, Keage and Blundy1979; Iken, Reference Iken1981; Cuffey and Paterson, Reference Cuffey and Paterson2010; Helanow and others, Reference Helanow, Iverson, Woodard and Zoet2021). At low stresses, viscous laws tend to produce faster sliding than plastic laws, but deformation in the ice column generally dominates mass flux for both parameterizations in central regions of ice sheets with low stresses. (For simplicity, we will refer to the end-member viscous and plastic laws.) Where sliding is dominant, plastic laws give larger velocity changes for shear-stress perturbations; this means that reduction in ice-shelf buttressing causes significant flow acceleration in the grounding zone but minimal thinning there, with the faster flow propagating well inland (Parizek and others, Reference Parizek2013). In contrast, for viscous behavior, most of the response to loss of ice-shelf buttressing is localized in the grounding zone, giving relatively little acceleration well inland but favoring grounding-zone thinning that speeds retreat from the stabilizing sill. For models lacking cliff calving, the plastic behavior tends to give larger sea-level rise in response to reduced buttressing as in the ABUMIP experiments; however, if retreat from the sill is followed by ice-shelf loss and then very rapid cliff calving, a viscous bed might drive larger sea-level rise (Pollard and others, Reference Pollard, DeConto and Alley2015; DeConto and others, Reference DeConto2021).

Sliding laws often include dependence on the basal water system. Frozen-bedded glaciers lacking a water system have very low basal velocity. For melted beds with basal water, the sliding velocity tends to increase with water supply up to some limit, beyond which efficient channelization of flow can partially drain the bed and cause decrease of basal velocity with increasing water supply (Cuffey and Paterson, Reference Cuffey and Paterson2010). In the inefficient drainage regime, which likely applies to most of the area of the beds of ice sheets (e.g. Alley and others, Reference Alley1997), additional water supply tends to raise the water pressure, growing water-filled cavities that separate ice from bed and thus facilitate sliding, and weakening tills by supporting more of the overburden pressure that otherwise tends to interlock grains so they cannot move past each other. Some ice-flow models include, or can include, such a dependence of basal velocity on the basal water system, although only one of the ABUMIP models included a water-pressure-dependent sliding law (Kazmierczak and others (Reference Kazmierczak, Sun, Coulon and Pattyn2022) provide additional insights to the importance of the water system to the ABUMIP-type experiments).

Extensive evidence (reviewed in, e.g. Alley and others, Reference Alley1997) shows that thawed till beds generally deform, and exhibit more-nearly-plastic behavior. Hard beds almost surely have a more-nearly-viscous behavior under at least many conditions, although the recent emphasis on the possibility that ice deformation depends on the fourth power of stress rather than the third (Millstein and others, Reference Millstein, Minchew and Pegler2022) may recommend testing a higher power in models, with potential consequences (see Sun and others, Reference Sun2020, Fig. 5). Models that experience an ‘Iken limit’ typically have been formulated assuming that the ice interacting with rigid bedrock lacks debris, and thus ignoring the dynamic (skin) drag from the widespread to ubiquitous occurrence of debris in basal ice, which may raise the limit or avoid it entirely (see, e.g. Zoet and others, Reference Zoet2020; Helanow and others, Reference Helanow, Iverson, Woodard and Zoet2021).

Importantly, extensive observations of deglaciated beds and beneath modern glaciers and ice sheets, including beneath Thwaites Glacier as summarized briefly below, show that the basal character is highly variable spatially in ways that will cause spatial variations in the sliding law. Models, however, almost never use a spatially varying basal sliding law that is registered to the real bed of the glacier to reflect the observed variations in bed type. Omitting plastic regions may induce only small error if they are always exceptionally well lubricated and thus unimportant in the force balance, but otherwise, models lacking this spatial information must be missing aspects of the actual behavior. Koellner and others (Reference Koellner, Parizek, Alley, Muto and Holschuh2019) conducted idealized flowline experiments forced from the downglacier end for domains in which spatially varying sliding laws were draped over basal topography, and found that while under many conditions the mixed bed produced behaviors that were intermediate between the end-member beds, under some conditions the mixed bed produced more-extreme behavior than either end-member. This is just one set of model experiments, but suggests the need to accurately represent the true bed with its appropriate sliding laws and their spatial extent, a task that will not be easy (see below).

Insights to the bed

Three-dimensional maps of large areas of ice-sheet beds generally have neither sufficient detail to resolve the main bedforms (e.g. the moats and the drumlins or flutes observed beneath Thwaites Glacier; see below), nor supporting geophysical data to characterize the nature of the bed (whether bedrock, till or other materials) or to characterize the water system. Some data with such resolution are available. Here, we focus especially on results of CReSIS surveys of the inland regions of Thwaites Glacier (Fig. 1a) (Muto and others, Reference Muto2019a, Reference Muto, Alley, Parizek and Anandakrishnan2019b; Clyne and others, Reference Clyne, Anandakrishnan, Muto, Alley and Voigt2020; Holschuh and others, Reference Holschuh, Christianson, Paden, Alley and Anandakrishnan2020, as summarized by Alley and others, Reference Alley2021), together with data from Rutford Ice Stream (Smith and others, Reference Smith2007) and Whillans Ice Stream (Barcheck and others, Reference Barcheck, Schwartz and Tulaczyk2020).

The map in Figure 1b shows the bed of Thwaites Glacier as surveyed by airborne swath radar (Holschuh and others, Reference Holschuh, Christianson, Paden, Alley and Anandakrishnan2020), with seismic data draped over it. As reviewed by Alley and others (Reference Alley2021), much of the West Antarctic ice sheet occupies the extended crust of the West Antarctic Rift System, with the Transantarctic Mountains as one rift flank, and Marie Byrd Land and other highlands as the other rift flank. Flow is preferentially directed along the fault-block topography, as on the Siple Coast, but the distance from the center of the ice sheet to the front of the Ross Ice Shelf is quite long, and was even longer when the ice was advanced at glacial maxima. This caused ice to build up in central regions, driving flow across the tectonic fabric and out the gap at Thwaites Glacier. While ice flow across topography occurs frequently in many environments, such flow is probably better developed for Thwaites Glacier than for most or all other regions.

Although additional geophysical confirmation would be helpful, the linear topographic features shown in Figure 1 oriented across flow are likely the normal faults of the rift system, with upthrown blocks on the downglacier sides. The ~100-m-high flutes in the lee of the higher parts of the upthrown blocks (the tails of crag-and-tail features) are seismically thick, soft till, which thus is at least a few meters thick and perhaps much more, as seismically observed internal structures are scarce or absent. The up-glacier-facing parts of these flutes, and of shorter features that might be called drumlins (Holschuh and others, Reference Holschuh, Christianson, Paden, Alley and Anandakrishnan2020), are seismically more rigid (or ‘hard’), probably poorly or moderately lithified sedimentary rocks. Seismically hard regions extend into the moats eroded along the up-glacier sides of the bedrock obstacles, and then downglacier in the troughs between the raised flutes.

This distribution of features is to some extent repeated beneath Rutford Ice Stream, and perhaps also beneath Whillans Ice Stream. Beneath Rutford, the seismically harder regions of the bed are on up-glacier-facing noses of bedforms and in troughs between flutes. Microseismicity shows that important tangential drag (dynamic drag, as opposed to the form drag of the ice needing to get around the bumps) is supported in these hard-bedded regions (Kufner and others, Reference Kufner2021), likely arising from interactions between bedrock and clasts carried in basal ice (e.g. Zoet and others, Reference Zoet2013). Microseismicity was also observed in trough locations beneath Whillans Ice Stream (Barcheck and others, Reference Barcheck, Schwartz and Tulaczyk2020), suggesting similar conditions there.

Occurrence in the Thwaites Glacier data of the large, bedrock-floored moats, incised as much as 50 m or more beneath the more-or-less planar surface up-glacier, was a special focus of the modeling in Alley and others (Reference Alley2021). Briefly, as shown first by Stokes (Reference Stokes1851) and subsequently in all similar studies including those of Weertman (Reference Weertman1957)-type sliding, the form drag of an obstacle creates a high-pressure zone on the upstream side to drive the flow divergence around the obstacle, with a coupled low-pressure zone on the downstream side. For ice flow over a bedrock knob, the high pressure is also exerted across the ice-bed interface, over a horizontal distance that scales with the obstacle size. For relatively rapid ice flow, hundreds of meters per year in the survey grid, the magnitude of this pressure perturbation caused by obstacles on their up-glacier sides is large compared to the pressure difference between the ice-overburden and the water pressure well away from the obstacle, and compared to changes in water pressure or ice-overburden pressure over similar length scales away from obstacles.

The soft beds of Thwaites Glacier are rather clearly in a regime for which increasing water pressure closer to the ice-overburden pressure tends to float the ice off the bed and reduce subglacial till deformation (Muto and others, Reference Muto, Alley, Parizek and Anandakrishnan2019b; Hansen and Zoet, Reference Hansen and Zoet2022). The zone of high ice pressure on the up-glacier side of an obstacle tends to prevent inflow of through-going subglacial water from up-glacier, which instead is routed far around the obstacle. Lack of this water then allows tighter coupling of the ice into the till in the high-pressure zone, increasing till deformational flux, and causing flux divergence that removes the till, exposing bedrock to erosion; this process acting over time likely explains the moats. The exclusion zone for throughgoing water extends to the side of an obstacle over a similar length scale, explaining the bedrock in troughs. Till transported over obstacles and deposited in lee-side positions then forms the large lee-side drumlins or flutes.

This stoss-side high-pressure behavior applies to all scales of obstacles. As described in Alley and others (Reference Alley2021), it likely lowers the melting point of ice up-glacier of abrading clasts in basal ice, causing heat-flow convergence that drives melting at the ice-bed interface and thus increases downward ice motion, resulting in an increase in the stress between the bed and abrading clasts within the basal ice. This stoss-side high pressure also occurs for all fixed obstacles in the bed, and contributes to deviations in water routing from paths that otherwise would be expected. Notably, though, the bed slope into a moat tends to capture basal water flow; learning whether moats extend far enough up-glacier and have steep enough bed slopes to offset the diversion of water away from the high-pressure region will require high-quality data and careful modeling.

Research priorities

After COVID-induced loss of field seasons and ongoing disruptions, the GHOST project (Geophysical Habitat Of Subglacial Thwaites; https://thwaitesglacier.org/projects/ghost) deployed a limited field season in 2022–2023 to collect new data to extend the high-resolution geophysical datasets noted above along and across larger regions of Thwaites Glacier, and to assess spatial occurrence of microearthquakes, although with major COVID-forced reductions in the planned program. One goal of this work is to test the hypothesis that this template is widespread: that bedrock is exposed on stoss sides of topography, and in moats with trailing troughs eroded around up-glacier sides of topographic features, with soft till in lee-side positions. Microearthquake data can be combined with results from laboratory experiments to understand the sliding law in the bedrock regions (e.g. Zoet and others, Reference Zoet2013, Reference Zoet2020). Careful surveys with targeted modeling can be used to calculate flow paths for subglacial water, and learn which downstream regions are, or are not, lubricated by water flow from up-glacier. As discussed in Alley and others (Reference Alley2021), the bedforms of Thwaites Glacier have experienced more-nearly steady conditions for longer times than the bedforms of most deglaciated regions, so similarities and differences may provide additional insights to formation mechanisms of bedforms.

To improve accuracy of ice-sheet projections, though, new data from Thwaites Glacier and elsewhere will need to be used in prognostic models. As noted above, whole-ice-sheet models generally are not resolving such topographic features, and so they are not applying the correct sliding laws in spatially accurate ways. Furthermore, whole-ice-sheet models generally are not routing subglacial water, with its effects on lubrication, in ways that include the pressure variations generated by flow against bedrock topography. Also, the Iken-limit sliding models may not be accurately parameterizing the frictional behavior when the effects of abrading clasts are also considered. Such detailed modeling likely is most important on the sill where Thwaites Glacier currently ends, and in the basin immediately up-glacier (e.g. Alley and others, Reference Alley, Li, Parizek and Zhang2019), but extending accurate modeling farther inland and ultimately to the whole ice sheet may be necessary. Modeling should address the scales at which bed types and processes need to be resolved to provide accurate projections. All of this is a huge challenge, and various intermediate models may provide guidance on the effects of including more of these processes.

Acknowledgements

We thank Doug MacAyeal, who was essential in the research but who is an editor of this volume and chose not to be a coauthor to avoid perceived conflict of interest. We thank Andy Smith and other colleagues in GHOST and other parts of the International Thwaites Glacier Collaboration (ITGC), many of whom were in the field as this was being written. We also thank colleagues at the Center for Remote Sensing of Ice Sheets (CReSIS). This research was supported by the U.S. National Science Foundation under grants AGS-1338832 and NSF-NERC-OPP-1738934, NASA under grant 80NSSC22K0384, and the Heising-Simons Foundation under grant 2018-0769. The seismic and radar data for this study were published with the papers cited. No new data are presented here; data are available through the papers cited.

References

Åkesson, H, Morlighem, M, O'Regan, M and Jakobsson, M (2021) Future projections of Petermann Glacier under ocean warming depend strongly on friction law. Journal of Geophysical Research: Earth Surface 126, e2020JF005921. doi: 10.1029/2020JF005921.Google Scholar
Alley, RB and 5 others (1997) How glaciers entrain and transport basal sediment: physical constraints. Quaternary Science Reviews 16, 1017–1038. doi: 10.1016/S0277-3791(97)00034-6.Google Scholar
Alley, RB and 11 others (2021) Bedforms of Thwaites Glacier, West Antarctica: character and origin. Journal of Geophysical Research – Earth Surface 126, e2021JF006339. doi: 10.1029/2021JF006339.Google Scholar
Alley, RB and 8 others (2023) Iceberg calving: regimes and transitions. Annual Reviews of Earth and Planetary Sciences 51, 189–215. doi: 10.1146/annurev-earth-032320-11.Google Scholar
Alley, RB, Li, W, Parizek, BR and Zhang, F (2019) Evaluation of ice-stream model sensitivities for parameter estimation. Earth and Planetary Science Letters 516, 49–55. doi: 10.1016/j.epsl.2019.03.035.Google Scholar
Barcheck, CG, Schwartz, SY and Tulaczyk, S (2020) Icequake streaks linked to potential mega-scale glacial lineations beneath an Antarctic ice stream. Geology 48, 99–102. doi: 10.1130/G46626.1.Google Scholar
Bassis, JN and Walker, CC (2012) Upper and lower limits on the stability of calving glaciers from the yield strength envelope of ice. Proceedings of the Royal Society of London A 468, 913–931. doi: 10.1098/rspa.2011.0422.Google Scholar
Brondex, J, Gillet-Chaulet, F and Gagliardini, O (2019) Sensitivity of centennial mass loss projections of the Amundsen basin to the friction law. The Cryosphere 13, 177–195. doi: 10.5194/tc-13-177-2019.Google Scholar
Budd, W, Keage, P and Blundy, N (1979) Empirical studies of ice sliding. Journal of Glaciology 23(89), 157–170. doi: 10.3189/S0022143000029804.Google Scholar
Clerc, F, Minchew, BM and Behn, MD (2019) Marine ice cliff instability mitigated by slow removal of ice shelves. Geophysical Research Letters 46, 12108–12116. doi: 10.1029/2019GL084183.Google Scholar
Clyne, ER, Anandakrishnan, S, Muto, A, Alley, RB and Voigt, DE (2020) Reflection seismic interpretation of topography and acoustic impedance beneath Thwaites Glacier, West Antarctica. Earth and Planetary Science Letters 550, 116543. doi: 10.1016/j.epsl.2020.116543.Google Scholar
Cuffey, KM and Paterson, WSB (2010) The Physics of Glaciers, 4th Edn. Oxford: Butterworth-Heinemann.Google Scholar
DeConto, RM and 12 others (2021) The Paris climate agreement and future sea-level rise from Antarctica. Nature 593, 83–89. doi: 10.1038/s41586-021-03427-0.Google Scholar
Hansen, DD and Zoet, LK (2022) Characterizing sediment flux of deforming glacier beds. Journal of Geophysical Research Earth Surface 127, e2021JF006544. doi: 10.1029/2021JF006544.Google Scholar
Hanson, B and Hooke, R (2003) Buckling rate and overhang development at a calving face. Journal of Glaciology 49, 577–586. doi: 10.3189/172756503781830476.Google Scholar
Helanow, C, Iverson, NR, Woodard, JB and Zoet, LK (2021) A slip law for hard-bedded glaciers derived from observed bed topography. Science Advances 7(20), eabe7798. doi: 10.1126/sciadv.abe7798.Google Scholar
Hoffman, AO and 5 others (2022) The impact of basal roughness on inland Thwaites Glacier sliding. Geophysical Research Letters 49, e2021GL096564. doi: 10.1029/2021GL096564.Google Scholar
Holschuh, N, Christianson, K, Paden, J, Alley, RB and Anandakrishnan, S (2020) Linking postglacial landscapes to glacier dynamics using swath radar at Thwaites Glacier. Geology 48, 268–272. doi: 10.1130/G46772.1.Google Scholar
Horgan, HJ and 7 others (2013) Estuaries beneath ice sheets. Geology 41(11), 1159–1162. doi: 10.1130/G34654.1.Google Scholar
Iken, A (1981) The effect of the subglacial water pressure on the sliding velocity of a glacier in an idealized numerical model. Journal of Glaciology 27, 407–421. doi: 10.3189/S0022143000011448.Google Scholar
Joughin, I, Smith, BE and Schoof, CG (2019) Regularized Coulomb friction laws for ice sheet sliding: application to Pine Island Glacier, Antarctica. Geophysical Research Letters 46, 4764–4771. doi: 10.1029/2019GL082526.Google Scholar
Kazmierczak, E, Sun, S, Coulon, V and Pattyn, F (2022) Subglacial hydrology modulates basal sliding response of the Antarctic ice sheet to climate forcing. The Cryosphere 16, 4537–4552. doi: 10.5194/tc-16-4537-2022.Google Scholar
Koellner, S, Parizek, BR, Alley, RB, Muto, A and Holschuh, N (2019) The impact of spatially-variable basal properties on outlet glacier flow. Earth and Planetary Science Letters 515, 200–208. doi: 10.1016/j.epsl.2019.03.026.Google Scholar
Kufner, S-K and 8 others (2021) Not all icequakes are created equal: basal icequakes suggest diverse bed deformation mechanisms at Rutford Ice Stream, West Antarctica. Journal of Geophysical Research: Earth Surface 126, e2020JF006001. doi: 10.1029/2020JF006001.Google Scholar
Millstein, JD, Minchew, BM and Pegler, SS (2022) Ice viscosity is more sensitive to stress than commonly assumed. Communications Earth and Environment 3, 57. doi: 10.1038/s43247-022-00385-x.Google Scholar
Morlighem, M (2020) MEaSUREs BedMachine Antarctica, Version 2 [Bed Elevation]. Boulder, CO, USA: NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: 10.5067/E1QL9HFQ7A8M (Accessed 3 February 2021).Google Scholar
Muto, A and 7 others (2019 a) Relating bed character and subglacial morphology using seismic data from Thwaites Glacier, West Antarctica. Earth and Planetary Science Letters 507, 199–206. doi: 10.1016/j.epsl.2018.12.008.Google Scholar
Muto, A, Alley, RB, Parizek, BR and Anandakrishnan, S (2019 b) Bed-type variability and till (dis)continuity beneath Thwaites Glacier, West Antarctica. Annals of Glaciology 60(80), 1–9. doi: 10.1017/aog.2019.32.Google Scholar
Parizek, BR and 10 others (2013) Dynamic (In)stability of Thwaites Glacier, West Antarctica. Journal of Geophysical Research-Earth Surface 118, 1–18. doi: 10.1002/jgrf.20044.Google Scholar
Parizek, BR and 7 others (2019) Ice-cliff failure via retrogressive slumping. Geology 47, 449–452. doi: 10.1130/G45880.1.CrossRefGoogle Scholar
Pollard, D, DeConto, RM and Alley, RB (2015) Potential Antarctic ice sheet retreat driven by hydrofracturing and ice cliff failure. Earth and Planetary Science Letters 412, 112–121. doi: 10.1016/j.epsl.2014.12.035.Google Scholar
Scambos, T and 22 others (2017) How much, how fast?: a review and science plan for research on the instability of Antarctica's Thwaites Glacier in the 21st century. Global and Planetary Change 153, 16–34. doi: 10.1016/j.gloplacha.2017.04.008.Google Scholar
Schwans, E, Parizek, BR, Alley, RB, Anandakrishnan, S and Morlighem, MM (2023) Model output for model insights into bed control on retreat of Thwaites Glacier, West Antarctica. Journal of Glaciology, 1–19. doi: 10.1017/jog.2023.13.Google Scholar
Smith, AM and 6 others (2007) Rapid erosion, drumlin formation, and changing hydrology beneath an Antarctic ice stream. Geology 35(2), 127–130. doi: 10.1130/G23036A.1.Google Scholar
Stokes, GG (1851) On the effect of internal friction of fluids on the motion of pendulums. Transactions of the Cambridge Philosophical Society 9(part ii), 8–106. doi: 10.1017/CBO9780511702266.002.Google Scholar
Sun, S and 28 others (2020) Antarctic ice sheet response to sudden and sustained ice-shelf collapse (ABUMIP). Journal of Glaciology 66, 891–904. doi: 10.1017/jog.2020.67.Google Scholar
Weertman, J (1957) On the sliding of glaciers. Journal of Glaciology 3(21), 33–38. doi: 10.3189/S0022143000024709.Google Scholar
Zoet, LK and 6 others (2013) The effects of entrained debris on the basal sliding stability of a glacier. Journal of Geophysical Research, Earth Surface 118, 656–666. doi: 10.1002/jgrf.20052.Google Scholar
Zoet, LK and 6 others (2020) Application of constitutive friction laws to glacier seismicity. Geophysical Research Letters 47, e2020GL088964. doi: 10.1029/2020GL088964.Google Scholar
Figure 0

Figure 1. Key features of the bed of Thwaites Glacier. (a). Bed elevation from BedMachine Antarctica version 2 (Morlighem, 2020), with elevation indicated by color scale, ice-shelf front (gray), and grounding line (black), with survey location and seismic lines in red. Thwaites discharges primarily between ~109 and 104 W longitude, across a relatively shallow grounding zone. (b). Bed of the up-glacier part of Thwaites Glacier (Holschuh and others, 2020), with elevation (left color bar, m; note that scale differs from 1a, as indicated) from airborne swath radar, and acoustic impedance (right color bar, kg m−2 s−1) from seismic surveys (Muto and others, 2019a, 2019b). Lower acoustic impedance indicates softer bed (deforming tills), and higher acoustic impedance indicates harder bed (bedrock). Holschuh and other (2020) and Alley and others (2021) provide additional, annotated data on the bed character. Radar data and gridded topographies are accessible through the University of Washington ResearchWorks Archive (http://hdl.handle.net/1773/44950).