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: Antarctic Topographic Realizations and Geostatistical Modeling Used to Map Subglacial Lakes
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
Award ID(s):
1745137
PAR ID:
10375171
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
125
Issue:
3
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. null (Ed.)
    Abstract Subglacial topography is an important feature in numerous ice-sheet analyses and can drive the routing of water at the bed. Bed topography is primarily measured with ice-penetrating radar. Significant gaps, however, remain in data coverage that require interpolation. Topographic interpolations are typically made with kriging, as well as with mass conservation, where ice flow dynamics are used to constrain bed geometry. However, these techniques generate bed topography that is unrealistically smooth at small scales, which biases subglacial water flowpath models and makes it difficult to rigorously quantify uncertainty in subglacial drainage patterns. To address this challenge, we adapt a geostatistical simulation method with probabilistic modeling to stochastically simulate bed topography such that the interpolated topography retains the spatial statistics of the ice-penetrating radar data. We use this method to simulate subglacial topography using mass conservation topography as a secondary constraint. We apply a water routing model to each of these realizations. Our results show that many of the flowpaths significantly change with each topographic realization, demonstrating that geostatistical simulation can be useful for assessing confidence in subglacial flowpaths. 
    more » « less
  3. 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
  4. Hercules Dome, Antarctica, has long been identified as a prospective deep ice core site due to the undisturbed internal layering, climatic setting and potential to obtain proxy records from the Last Interglacial (LIG) period when the West Antarctic ice sheet may have collapsed. We performed a geophysical survey using multiple ice-penetrating radar systems to identify potential locations for a deep ice core at Hercules Dome. The surface topography, as revealed with recent satellite observations, is more complex than previously recognized. The most prominent dome, which we term ‘West Dome’, is the most promising region for a deep ice core for the following reasons: (1) bed-conformal radar reflections indicate minimal layer disturbance and extend to within tens of meters of the ice bottom; (2) the bed is likely frozen, as evidenced by both the shape of the measured vertical ice velocity profiles beneath the divide and modeled ice temperature using three remotely sensed estimates of geothermal flux and (3) models of layer thinning have 132 ka old ice at 45–90 m above the bed with an annual layer thickness of ~1 mm, satisfying the resolution and preservation needed for detailed analysis of the LIG period. 
    more » « less
  5. Abstract Subglacial lakes require a thawed bed either now or in the past; thus, their presence and stability have implications for current and past basal conditions, ice dynamics, and climate. Here, we present the most extensive geophysical exploration to date of a subglacial lake near the geographic South Pole, including radar‐imaged stratigraphy, surface velocities, and englacial vertical velocities. We use a 1.5‐dimensional temperature model, optimized with our geophysical data set and nearby temperature measurements, to estimate past basal‐melt rates. The ice geometry, reflected bed‐echo power, surface and vertical velocities, and temperature model indicate that the ice‐bed interface is regionally thawed, contradicting prior studies. Together with an earlier active‐source seismic study, which showed a 32‐m deep lake underlain by 150 m of sediment, our results suggest that the lake has been thermodynamically stable through at least the last 120,000 years and possibly much longer, making it a promising prospective site for sediment coring. 
    more » « less