An official website of the United States government
Abstract. Ice-sheet models used to predict sea-level rise often neglect subglacial hydrology. However, theory and observations suggest that ice flow and subglacial water flow are bidirectionally coupled: ice geometry affects hydraulic potential, hydraulic potential modulates basal shear stress via the basal water pressure, and ice flow advects the subglacial drainage system. This coupling could impact rates of ice mass change but remains poorly understood. We develop a coupled ice–subglacial-hydrology model to investigate the effects of coupling on the long-term evolution of marine-terminating ice sheets. We combine a one-dimensional channelized subglacial hydrology model with a depth-integrated marine-ice-sheet model, incorporating each component of the coupling listed above, yielding a set of differential equations that we solve using a finite-difference, implicit time-stepping approach. We conduct a series of experiments with this model, using either bidirectional or unidirectional coupling. These experiments generate profiles of channel cross-sectional area, channel flow rate, channel effective pressure, ice thickness, and ice velocity. We discuss how the profiles shape one another, resulting in the effective pressure reaching a local maximum in a region near the grounding line. We also describe the impact of bidirectional coupling on the transient retreat of ice sheets through a comparison of our coupled model with ice-flow models that have imposed static basal conditions. We find that including coupled subglacial hydrology leads to grounding-line retreat that is virtually absent when static basal conditions are assumed. This work highlights the role time-evolving subglacial drainage may have in ice-sheet change and informs efforts to include it in ice-sheet models. This work also supplies a physical basis for a commonly used parameterization which assumes that the subglacial water pressure is set by the bed's depth beneath the sea surface.
more »
« less
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
- 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
-
-
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
-
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
-
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
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.394
