An official website of the United States government
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
Modeling oscillations in connected glacial lakes
Mountain glaciers and ice sheets often host marginal and subglacial lakes that are hydraulically connected through subglacial drainage systems. These lakes exhibit complex dynamics that have been the subject of models for decades. Here we introduce and analyze a model for the evolution of glacial lakes connected by subglacial channels. Subglacial channel equations are supplied with effective pressure boundary conditions that are determined by a simple lake model. While the model can describe an arbitrary number of lakes, we solve it numerically with a finite element method for the case of two connected lakes. We examine the effect of relative lake size and spacing on the oscillations. Complex oscillations in the downstream lake are driven by discharge out of the upstream lake. These include multi-peaked and anti-phase filling–draining events. Similar filling–draining cycles have been observed on the Kennicott Glacier in Alaska and at the confluence of the Whillans and Mercer ice streams in West Antarctica. We further construct a simplified ordinary differential equation model that displays the same qualitative behavior as the full, spatially-dependent model. We analyze this model using dynamical systems theory to explain the appearance of filling–draining cycles as the meltwater supply varies.
more »
« less
- Award ID(s):
- 1643970
- PAR ID:
- 10144205
- Date Published:
- Journal Name:
- Journal of Glaciology
- Volume:
- 65
- Issue:
- 253
- ISSN:
- 1727-5652
- Page Range / eLocation ID:
- 745-758
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. The Greenland Ice Sheet (GrIS) is losing mass as the climate warms through both increased meltwater runoff and ice discharge at marine-terminating sectors. At the ice sheet surface, meltwater runoff forms a dynamic supraglacial hydrological system which includes stream and river networks and large supraglacial lakes (SGLs). Streams and rivers can route water into crevasses or into supraglacial lakes with crevasses underneath, both of which can then hydrofracture to the ice sheet base, providing a mechanism for the surface meltwater to access the bed. Understanding where, when, and how much meltwater is transferred to the bed is important because variability in meltwater supply to the bed can increase ice flow speeds, potentially impacting the hypsometry of the ice sheet in grounded sectors, and iceberg discharge to the ocean. Here we present a new, physically based, supraglacial hydrology model for the GrIS that is able to simulate (a) surface meltwater routing and SGL filling; (b) rapid meltwater drainage to the ice sheet bed via the hydrofracture of surface crevasses both in and outside of SGLs; (c) slow SGL drainage via overflow in supraglacial meltwater channels; and, by offline coupling with a second model, (d) the freezing and unfreezing of SGLs from autumn to spring. We call the model the Supraglacial Hydrology Evolution and Drainage (or SHED) model. We apply the model to three study regions in southwest Greenland between 2015 and 2019 (inclusive) and evaluate its performance with respect to observed supraglacial lake extents and proglacial discharge measurements. We show that the model reproduces 80 % of observed lake locations and provides good agreement with observations in terms of the temporal evolution of lake extent. Modelled moulin density values are in keeping with those previously published, and seasonal and inter-annual variability in proglacial discharge agrees well with that which is observed, though the observations lag the model by a few days since they include transit time through the subglacial system, while the model does not. Our simulations suggest that lake drainage behaviours may be more complex than traditional models suggest, with lakes in our model draining through a combination of both overflow and hydrofracture and with some lakes draining only partially and then refreezing. This suggests that, in order to simulate the evolution of Greenland's surface hydrological system with fidelity, a model that includes all of these processes needs to be used. In future work, we will couple our model to a subglacial model and an ice flow model and thus use our estimates of where, when, and how much meltwater gets to the bed to understand the consequences for ice flow.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
-
Morlighem, Mathieu (Ed.)Abstract Subglacial lakes have been mapped across Antarctica with two methods, radio‐echo sounding (RES) and ice‐surface deformation. At sites where both are coincident, these methods typically provide conflicting interpretations about the ice‐bed interface. With a single exception,activesubglacial lakes identified by surface deformation do not display the expected flat, bright, and specular bed reflection in RES data, characteristic ofnon‐activelakes. This observational conundrum suggests that our understanding of Antarctic subglacial hydrology, especially beneath important fast‐moving ice streams, remains incomplete. Here, we use an airborne RES campaign that surveyed a well‐characterized group of active subglacial lakes on lower Mercer and Whillans ice streams, West Antarctica, to explore inconsistency between the two observational techniques. We test hypotheses of increased scattering and attenuation due to the presence of an active subglacial lake system that could suppress reflected bed‐echo power for RES observations in these locations, finding that entrained water is most plausible.more » « less
-
Abstract Over the past 50 years, the discovery and initial investigation of subglacial lakes in Antarctica have highlighted the paleoglaciological information that may be recorded in sediments at their beds. In December 2018, we accessed Mercer Subglacial Lake, West Antarctica, and recovered the first in situ subglacial lake-sediment record—120 mm of finely laminated mud. We combined geophysical observations, image analysis, and quantitative stratigraphy techniques to estimate long-term mean lake sedimentation rates (SRs) between 0.49 ± 0.12 mm a–1 and 2.3 ± 0.2 mm a–1, with a most likely SR of 0.68 ± 0.08 mm a–1. These estimates suggest that this lake formed between 53 and 260 a before core recovery (BCR), with a most likely age of 180 ± 20 a BCR—coincident with the stagnation of the nearby Kamb Ice Stream. Our work demonstrates that interconnected subglacial lake systems are fundamentally linked to larger-scale ice dynamics and highlights that subglacial sediment archives contain powerful, century-scale records of ice history and provide a modern process-based analogue for interpreting paleo–subglacial lake facies.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/https://doi.org/10.1017/
