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
‘Stable’ and ‘unstable’ are not useful descriptions of marine ice sheets in the Earth's climate system
Investigations of the time-dependent behavior of marine ice sheets and their sensitivity to basal conditions require numerical models because existing theoretical analyses focus only on steady-state configurations primarily with a power-law basal shear stress. Numerical results indicate that the choice of the sliding law strongly affects ice-sheet dynamic behavior. Although observed or simulated grounding-line retreat is typically interpreted as an indication of marine ice sheet instability introduced by Weertman (1974), this (in)stability is a characteristic of the ice sheet's steady states – not time-variant behavior. To bridge the gap between theoretical and numerical results, we develop a framework to investigate grounding line dynamics with generalized basal and lateral stresses (i.e. the functional dependencies are not specified). Motivated by observations of internal variability of the Southern Ocean conditions we explore the grounding-line response to stochastic variability. We find that adding stochastic variability to submarine melt rates that produced stable steady-state configurations leads to intermittently advancing and retreating grounding lines. They can also retreat in an unstoppable manner on time-scales significantly longer than the stochastic correlation time-scales. These results suggest that at any given time of their evolution, the transient behavior of marine ice sheets cannot be described in terms of ‘stable’ or ‘unstable’.
more »
« less
- Award ID(s):
- 2218463
- PAR ID:
- 10436317
- Date Published:
- Journal Name:
- Journal of Glaciology
- ISSN:
- 0022-1430
- Page Range / eLocation ID:
- 1 to 17
- 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 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
-
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
-
Marine ice stream dynamics are sensitive to conditions at the grounding line and basal shear stress. Variations in subglacial hydrology have been implicated in ice stream speed-up and shutdown. To investigate the interplay between marine ice stream flow and subglacial hydrology, we couple models of a marine ice stream and a subglacial drainage system. The coupled system evolves dynamically due to a positive feedback between ice flow, heat dissipation at the ice stream bed, and basal lubrication. Our results show that, depending on the hydraulic conductivity of the bed, distinct dynamic regimes can be identified. These regimes include steady streaming, hydraulically controlled oscillations, and thermally controlled oscillations. Periodic fast flow can be initiated by activation waves travelling upstream or downstream or quasi-simultaneously everywhere. Different dynamical regimes are characterized by large differences in grounded ice volume, even under modest changes of grounding line positions. These results imply a strong dependence of marine ice sheet dynamics on evolving hydrological conditions at the bed and highlight the importance of a better understanding of subglacial hydrology.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1017/jog.2023.40
