skip to main content


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
NSF-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. The Amundsen Sea sector of Antarctica has long been considered the most vulnerable part of the West Antarctic Ice Sheet (WAIS) because of the great water depth at the grounding line, a subglacial bed seafloor deepening toward the interior of the continent, and the absence of substantial ice shelves. Glaciers in this configuration are thought to be susceptible to rapid or runaway retreat. Ice flowing into the Amundsen Sea Embayment is undergoing the most rapid changes of any sector of the Antarctic ice sheets outside the Antarctic Peninsula, including substantial grounding-line retreat over recent decades, as observed from satellite data. Recent models suggest that a threshold leading to the collapse of WAIS in this sector may have been already crossed and that much of the ice sheet could be lost even under relatively moderate greenhouse gas emission scenarios. Drill cores from the Amundsen Sea provide tests of several key questions about controls on ice sheet stability. The cores offer a direct offshore record of glacial history in a sector that is exclusively influenced by ice draining the WAIS, which allows clear comparisons between the WAIS history and low-latitude climate records. Today, relatively warm (modified) Circumpolar Deep Water (CDW) is impinging onto the Amundsen Sea shelf and causing melting under ice shelves and at the grounding line of the WAIS in most places. Reconstructions of past CDW intrusions can assess the ties between warm water upwelling and large-scale changes in past grounding-line positions. Carrying out these reconstructions offshore from the drainage basin that currently has the most substantial negative mass balance of ice anywhere in Antarctica is thus of prime interest to future predictions. The scientific objectives for this expedition are built on hypotheses about WAIS dynamics and related paleoenvironmental and paleoclimatic conditions. The main objectives are: 1. To test the hypothesis that WAIS collapses occurred during the Neogene and Quaternary and, if so, when and under which environmental conditions; 2. To obtain ice-proximal records of ice sheet dynamics in the Amundsen Sea that correlate with global records of ice-volume changes and proxy records for atmospheric and ocean temperatures; 3. To study the stability of a marine-based WAIS margin and how warm deepwater incursions control its position on the shelf; 4. To find evidence for the earliest major grounded WAIS advances onto the middle and outer shelf; 5. To test the hypothesis that the first major WAIS growth was related to the uplift of the Marie Byrd Land dome. International Ocean Discovery Program (IODP) Expedition 379 completed two very successful drill sites on the continental rise of the Amundsen Sea. Site U1532 is located on a large sediment drift, now called the Resolution Drift, and it penetrated to 794 m with 90% recovery. We collected almost-continuous cores from recent age through the Pleistocene and Pliocene and into the upper Miocene. At Site U1533, we drilled 383 m (70% recovery) into the more condensed sequence at the lower flank of the same sediment drift. The cores of both sites contain unique records that will enable study of the cyclicity of ice sheet advance and retreat processes as well as ocean-bottom water circulation and water mass changes. In particular, Site U1532 revealed a sequence of Pliocene sediments with an excellent paleomagnetic record for high-resolution climate change studies of the previously sparsely sampled Pacific sector of the West Antarctic margin. Despite the drilling success at these sites, the overall expedition experienced three unexpected difficulties that affected many of the scientific objectives: 1. The extensive sea ice on the continental shelf prevented us from drilling any of the proposed shelf sites. 2. The drill sites on the continental rise were in the path of numerous icebergs of various sizes that frequently forced us to pause drilling or leave the hole entirely as they approached the ship. The overall downtime caused by approaching icebergs was 50% of our time spent on site. 3. A medical evacuation cut the expedition short by 1 week. Recovery of core on the continental rise at Sites U1532 and U1533 cannot be used to indicate the extent of grounded ice on the shelf or, thus, of its retreat directly. However, the sediments contained in these cores offer a range of clues about past WAIS extent and retreat. At Sites U1532 and U1533, coarse-grained sediments interpreted to be ice-rafted debris (IRD) were identified throughout all recovered time periods. A dominant feature of the cores is recorded by lithofacies cyclicity, which is interpreted to represent relatively warmer periods variably characterized by sediments with higher microfossil abundance, greater bioturbation, and higher IRD concentrations alternating with colder periods characterized by dominantly gray laminated terrigenous muds. Initial comparison of these cycles to published late Quaternary records from the region suggests that the units interpreted to be records of warmer time intervals in the core tie to global interglacial periods and the units interpreted to be deposits of colder periods tie to global glacial periods. Cores from the two drill sites recovered sediments of dominantly terrigenous origin intercalated or mixed with pelagic or hemipelagic deposits. In particular, Site U1533, which is located near a deep-sea channel originating from the continental slope, contains graded silts, sands, and gravels transported downslope from the shelf to the rise. The channel is likely the pathway of these sediments transported by turbidity currents and other gravitational downslope processes. The association of lithologic facies at both sites predominantly reflects the interplay of downslope and contouritic sediment supply with occasional input of more pelagic sediment. Despite the lack of cores from the shelf, our records from the continental rise reveal the timing of glacial advances across the shelf and thus the existence of a continent-wide ice sheet in West Antarctica during longer time periods since at least the late Miocene. Cores from both sites contain abundant coarse-grained sediments and clasts of plutonic origin transported either by downslope processes or by ice rafting. If detailed provenance studies confirm our preliminary assessment that the origin of these samples is from the plutonic bedrock of Marie Byrd Land, their thermochronological record will potentially reveal timing and rates of denudation and erosion linked to crustal uplift. The chronostratigraphy of both sites enables the generation of a seismic sequence stratigraphy for the entire Amundsen Sea continental rise, spanning the area offshore from the Amundsen Sea Embayment westward along the Marie Byrd Land margin to the easternmost Ross Sea through a connecting network of seismic lines. 
    more » « less
  2. null (Ed.)
    The Amundsen Sea sector of Antarctica has long been considered the most vulnerable part of the West Antarctic Ice Sheet (WAIS) because of the great water depth at the grounding line and the absence of substantial ice shelves. Glaciers in this configuration are thought to be susceptible to rapid or runaway retreat. Ice flowing into the Amundsen Sea Embayment is undergoing the most rapid changes of any sector of the Antarctic Ice Sheet outside the Antarctic Peninsula, including changes caused by substantial grounding-line retreat over recent decades, as observed from satellite data. Recent models suggest that a threshold leading to the collapse of WAIS in this sector may have been already crossed and that much of the ice sheet could be lost even under relatively moderate greenhouse gas emission scenarios. Drill cores from the Amundsen Sea provide tests of several key questions about controls on ice sheet stability. The cores offer a direct record of glacial history offshore from a drainage basin that receives ice exclusively from the WAIS, which allows clear comparisons between the WAIS history and low-latitude climate records. Today, warm Circumpolar Deep Water (CDW) is impinging onto the Amundsen Sea shelf and causing melting of the underside of the WAIS in most places. Reconstructions of past CDW intrusions can assess the ties between warm water upwelling and large-scale changes in past grounding-line positions. Carrying out these reconstructions offshore from the drainage basin that currently has the most substantial negative mass balance of ice anywhere in Antarctica is thus of prime interest to future predictions. The scientific objectives for this expedition are built on hypotheses about WAIS dynamics and related paleoenvironmental and paleoclimatic conditions. The main objectives are 1. To test the hypothesis that WAIS collapses occurred during the Neogene and Quaternary and, if so, when and under which environmental conditions; 2. To obtain ice-proximal records of ice sheet dynamics in the Amundsen Sea that correlate with global records of ice-volume changes and proxy records for atmospheric and ocean temperatures; 3. To study the stability of a marine-based WAIS margin and how warm deep-water incursions control its position on the shelf; 4. To find evidence for earliest major grounded WAIS advances onto the middle and outer shelf; 5. To test the hypothesis that the first major WAIS growth was related to the uplift of the Marie Byrd Land dome. International Ocean Discovery Program (IODP) Expedition 379 completed two very successful drill sites on the continental rise of the Amundsen Sea. Site U1532 is located on a large sediment drift, now called Resolution Drift, and penetrated to 794 m with 90% recovery. We collected almost-continuous cores from the Pleistocene through the Pliocene and into the late Miocene. At Site U1533, we drilled 383 m (70% recovery) into the more condensed sequence at the lower flank of the same sediment drift. The cores of both sites contain unique records that will enable study of the cyclicity of ice sheet advance and retreat processes as well as bottom-water circulation and water mass changes. In particular, Site U1532 revealed a sequence of Pliocene sediments with an excellent paleomagnetic record for high-resolution climate change studies of the previously sparsely sampled Pacific sector of the West Antarctic margin. Despite the drilling success at these sites, the overall expedition experienced three unexpected difficulties that affected many of the scientific objectives: 1. The extensive sea ice on the continental shelf prevented us from drilling any of the proposed shelf sites. 2. The drill sites on the continental rise were in the path of numerous icebergs of various sizes that frequently forced us to pause drilling or leave the hole entirely as they approached the ship. The overall downtime caused by approaching icebergs was 50% of our time spent on site. 3. An unfortunate injury to a member of the ship's crew cut the expedition short by one week. Recovery of core on the continental rise at Sites U1532 and U1533 cannot be used to precisely indicate the position of ice or retreat of the ice sheet on the shelf. However, these sediments contained in the cores offer a range of clues about past WAIS extent and retreat. At Sites U1532 and U1533, coarse-grained sediments interpreted to be ice-rafted debris (IRD) were identified throughout all recovered time periods. A dominant feature of the cores is recorded by lithofacies cyclicity, which is interpreted to represent relatively warmer periods variably characterized by higher microfossil abundance, greater bioturbation, and higher counts of IRD alternating with colder periods characterized by dominantly gray laminated terrigenous muds. Initial comparison of these cycles to published records from the region suggests that the units interpreted as records of warmer time intervals in the core tie to interglacial periods and the units interpreted as deposits of colder periods tie to glacial periods. The cores from the two drill sites recovered sediments of purely terrigenous origin intercalated or mixed with pelagic or hemipelagic deposits. In particular, Site U1533, which is located near a deep-sea channel originating from the continental slope, contains graded sands and gravel transported downslope from the shelf to the abyssal plain. The channel is likely the path of such sediments transported downslope by turbidity currents or other sediment-gravity flows. The association of lithologic facies at both sites predominantly reflects the interplay of downslope and contouritic sediment supply with occasional input of more pelagic sediment. Despite the lack of cores from the shelf, our records from the continental rise reveal the timing of glacial advances across the shelf and thus the existence of a continent-wide ice sheet in West Antarctica at least during longer time periods since the late Miocene. Cores from both sites contain abundant coarse-grained sediments and clasts of plutonic origin transported either by downslope processes or by ice rafting. If detailed provenance studies confirm our preliminary assessment that the origin of these samples is from the plutonic bedrock of Marie Byrd Land, their thermochronological record will potentially reveal timing and rates of denudation and erosion linked to crustal uplift. The chronostratigraphy of both sites enables the generation of a seismic sequence stratigraphy not only for the Amundsen Sea rise but also for the western Amundsen Sea along the Marie Byrd Land margin through a connecting network of seismic lines. 
    more » « less
  3. 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
  4. Abstract

    Projections of Antarctica's contribution to future sea level rise are associated with significant uncertainty, in part because the observational record is too short to capture long‐term processes necessary to estimate ice mass changes over societally relevant timescales. Records of grounding line retreat from the geologic past offer an opportunity to extend our observations of these processes beyond the modern record and to gain a more comprehensive understanding of ice‐sheet change. Here, we present constraints on the timing and inland extent of deglacial grounding line retreat in the southern Ross Sea, Antarctica, obtained via direct sampling of a subglacial lake located 150 km inland from the modern grounding line and beneath >1 km of ice. Isotopic measurements of water and sediment from the lake enabled us to evaluate how the subglacial microbial community accessed radiocarbon‐bearing organic carbon for energy, as well as where it transferred carbon metabolically. Using radiocarbon as a natural tracer, we found that sedimentary organic carbon was microbially translocated to dissolved carbon pools in the subglacial hydrologic system during the 4.5‐year period of water accumulation prior to our sampling. This finding indicates that the grounding line along the Siple Coast of West Antarctica retreated more than 250 km inland during the mid‐Holocene (6.3 ± 1.0 ka), prior to re‐advancing to its modern position.

     
    more » « less
  5. 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