skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Variational formulation of marine ice-sheet and subglacial-lake grounding-line dynamics
Grounding lines exist where land-based glacial ice flows on to a body of water. Accurately modelling grounding-line migration at the ice–ocean interface is essential for estimating future ice-sheet mass change. On the interior of ice sheets, the shores of subglacial lakes are also grounding lines. Grounding-line positions are sensitive to water volume changes such as sea-level rise or subglacial-lake drainage. Here, we introduce numerical methods for simulating grounding-line dynamics in the marine ice sheet and subglacial-lake settings. Variational inequalities arise from contact conditions that relate normal stress, water pressure and velocity at the base. Existence and uniqueness of solutions to these problems are established using a minimisation argument. A penalty method is used to replace the variational inequalities with variational equations that are solved using a finite-element method. We illustrate the grounding-line response to tidal cycles in the marine ice-sheet problem and filling–draining cycles in the subglacial-lake problem. We introduce two computational benchmarks where the known lake volume change is used to measure the accuracy of the numerical method.  more » « less
Award ID(s):
1643970
PAR ID:
10322244
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of fluid mechanics
Volume:
919
Issue:
a23
ISSN:
1469-7645
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Knowledge of past ice sheet configurations is useful for informing projections of future ice sheet dynamics and for calibrating ice sheet models. The topology of grounding line retreat in the Ross Sea sector of Antarctica has been much debated, but it has generally been assumed that the modern ice sheet is as small as it has been for more than 100 000 years (Conway et al., 1999; Lee et al., 2017; Lowry et al., 2019; McKay et al., 2016; Scherer et al., 1998). Recent findings suggest that the West Antarctic Ice Sheet (WAIS) grounding line retreated beyond its current location earlier in the Holocene and subsequently readvanced to reach its modern position (Bradley et al., 2015; Kingslake et al., 2018). Here, we further constrain the post-LGM (Last Glacial Maximum) grounding line retreat and readvance in the Ross Sea sector using a two-phase model of radiocarbon input and decay in subglacial sediments from six sub-ice sampling locations. In addition, we reinterpret high basal temperature gradients, measured previously at three sites in this region (Engelhardt, 2004), which we explain as resulting from recent ice shelf re-grounding accompanying grounding line readvance. Atone location – Whillans Subglacial Lake (SLW) – for which a sedimentporewater chemistry profile is known, we estimate the grounding linereadvance by simulating ionic diffusion. Collectively, our analyses indicate that the grounding line retreated over SLW 4300-2500+1500 years ago, and over sites on Whillans Ice Stream (WIS), Kamb Ice Stream (KIS), and Bindschadler Ice Stream (BIS) 4700-2300+1500, 1800-700+2700, and 1700-600+2800 years ago, respectively. The grounding line only recently readvanced back over those sites 1100-100+200, 1500-200+500, 1000-300+200, and 800±100 years ago for SLW, WIS, KIS, and BIS, respectively. The timing of grounding line retreat coincided with a warm period in the mid-Holocene to late Holocene. Conversely, grounding line readvance is coincident with cooling climate in the last 1000–2000 years. Our estimates for the timing of grounding line retreat and readvance are also consistent with relatively low carbon-to-nitrogen ratios measured in our subglacial sediment samples (suggesting a marine source of organic matter) and with the lack of grounding zone wedges in front of modern grounding lines. Based on these results, we propose that the Siple Coast grounding line motions in the mid-Holocene to late Holocene were primarily driven by relatively modest changes in regional climate, rather than by ice sheet dynamics and glacioisostatic rebound, as hypothesized previously (Kingslake et al., 2018). 
    more » « less
  2. Abstract. Widespread existing geological records from above the modern ice sheet surface and outboard of the current ice margin show that the Antarctic IceSheet (AIS) was much more extensive at the Last Glacial Maximum (∼ 20 ka) than at present. However, whether it was ever smaller thanpresent during the last few millennia, and (if so) by how much, is known only for a few locations because direct evidence lies within or beneath theice sheet, which is challenging to access. Here, we describe how retreat and readvance (henceforth “readvance”) of AIS grounding lines during theHolocene could be detected and quantified using subglacial bedrock, subglacial sediments, marine sediment cores, relative sea-level (RSL) records,geodetic observations, radar data, and ice cores. Of these, only subglacial bedrock and subglacial sediments can provide direct evidence forreadvance. Marine archives are of limited utility because readvance commonly covers evidence of earlier retreat. Nevertheless, stratigraphictransitions documenting change in environment may provide support for direct evidence from subglacial records, as can the presence of transgressionsin RSL records, and isostatic subsidence. With independent age control, ice structure revealed by radar can be used to infer past changes in iceflow and geometry, and therefore potential readvance. Since ice cores capture changes in surface mass balance, elevation, and atmosphericand oceanic circulation that are known to drive grounding line migration, they also have potential for identifying readvance. A multidisciplinaryapproach is likely to provide the strongest evidence for or against a smaller-than-present AIS in the Holocene. 
    more » « less
  3. Antarctic subglacial lakes can play an important role in ice sheet dynamics, biology, geology, and oceanography, but it is difficult to definitively constrain their character and locations. Subglacial lake locations are related to factors including heat flux, ice surface slope, ice thickness, and bed topography, though these relationships are not fully quantified. Bed topography is particularly important for determining where water flows and accumulates, but digital elevation models of the ice sheet bed rely on interpolation and are unrealistically smooth, biasing estimates of subglacial lake location and surface area. To address this issue, we use geostatistical methods to simulate realistically rough bed topography. We use our simulated topography to predict subglacial lake distribution across the continent using a binomial logistic regression, which uses physical parameters and known lake locations to calculate the probabilities of lake occurrences. Our results suggest that topography models interpolated without appropriate geostatistics overestimate subglacial lake surface area and that total lake surface area is lower than previously predicted. We find that radar‐detected lakes are more likely to occur in the interior of East Antarctica, while altimetry‐detected (active) lakes are expected to be found in West Antarctica and near the grounding line. We observe that radar‐detected lakes have a high correlation with heat flux and ice thickness, while active lakes are associated with higher ice velocity. 
    more » « less
  4. The Antarctic Ice Sheet contains enough ice above floatation to contribute ∼57 m of sea-level rise. About a third of that ice rests on bedrock far below sea level, in deep subglacial basins. In places where the bedrock slopes downward, away from the coast (retrograde bed) or where groundling lines terminate in deep water, the ice-sheet margin is susceptible to dynamic instabilities (related to both viscous and brittle processes) that can cause rapid retreat. Here, we use a numerical ice sheet-shelf model that captures these processes to demonstrate that the current configuration of the Antarctic Ice Sheet, with marine-terminating grounding lines poised on the edges of deep subglacial basins, combined with current, early 21st century atmospheric and oceanic tempera- tures, maximizes the danger of rapid retreat and sea-level rise. The potential retreat rates from the modern state greatly exceed those starting from other past configurations such as glacial maxima and smaller terrestrial ice cover. We also explore the importance of various negative feedbacks on the pace of retreat, including the effects of mélange, meltwater-climate feedback, and bedrock rebound. We find that today’s Antarctic Ice Sheet is capable of contributing faster sea-level rise than at any other time in the ice sheet’s recent (three million-year) geologic history, and possibly its entire 34 million-year history. This potential for unprecedented sea-level rise from Antarctica, coinciding with the recent expansion of human population centers and infrastructure along low-lying coastlines maximizes the risk of a future sea-level catastrophe. 
    more » « less
  5. Abstract Adélie‐George V Land in East Antarctica, encompassing the vast Wilkes Subglacial Basin, has a configuration that could be prone to ice sheet instability: the basin's retrograde bed slope could make its marine terminating glaciers vulnerable to warm seawater intrusion and irreversible retreat under predicted climate forcing. However, future projections are uncertain, due in part to limited subglacial observations near the grounding zone. Here, we develop a novel statistical approach to characterize subglacial conditions from radar sounding observations. Our method reveals intermixed frozen and thawed bed within 100 km of the grounding‐zone near the Wilkes Subglacial Basin outflow, and enables comparisons to ice sheet model‐inferred thermal states. The signs of intermixed or near thawed conditions raises the possibility that changes in basal thermal state could impact the stability of Adélie‐George V Land, adding to the region's potentially vulnerable topographic configuration and sensitivity to ocean forcing driven grounding line retreat. 
    more » « less