High-resolution modelling of the seasonal evolution of surface water storage on the Greenland Ice Sheet

Seasonal meltwater lakes on the Greenland Ice Sheet form when surface runoff is temporarily trapped in surface topographic depressions. The development of such lakes affects both the surface energy balance and dynamics of the ice sheet. Although areal extents, depths, and lifespans of lakes can be inferred from satellite imagery, such observational studies have a limited temporal resolution. Here, we adopt a modelling-based strategy to estimate the seasonal evolution of surface water storage for the ~3600 km2 Paakitsoq region of W. Greenland. We use a high-resolution time dependent surface mass balance model to calculate surface melt, a supraglacial water routing model to calculate lake filling and a prescribed water-volume based threshold to predict rapid lake drainage events. This threshold assumes that drainage will occur through a fracture if V = Fa.H, where V is lake volume, H is the local ice thickness and Fa is the potential fracture area. The model shows good agreement between modelled lake locations and volumes and those observed in 9 Landsat 7 ETM+ images from 2001, 2002 and 2005. We use the model to investigate the lake water volume required to trigger drainage, and the impact that varying this threshold volume has on the proportion of meltwater that is stored in surface lakes and enters the subglacial drainage system. Model performance is maximised with values of Fa between 4000 and 7500 m2. For these thresholds, lakes transiently store <40% of available meltwater at the beginning of the melt season, decreasing to ~5 to 10% by the middle of the melt season; over the course of a melt-season, 40 to 50% of total meltwater production enters the subglacial drainage system through moulins at the bottom of drained lakes.This work was partially funded by the UK Natural Environment Research Council Doctoral Training grant to A. F. Banwell (LCAG/ 133) (CASE Studentship with The Geological Survey of Denmark and Greenland (GEUS))

Controls on rapid supraglacial lake drainage in West Greenland: An Exploratory Data Analysis approach

ABSTRACTThe controls on rapid surface lake drainage on the Greenland ice sheet (GrIS) remain uncertain, making it challenging to incorporate lake drainage into models of GrIS hydrology, and so to determine the ice-dynamic impact of meltwater reaching the ice-sheet bed. Here, we first use a lake area and volume tracking algorithm to identify rapidly draining lakes within West Greenland during summer 2014. Second, we derive hydrological, morphological, glaciological and surface-mass-balance data for various factors that may influence rapid lake drainage. Third, these factors are used within Exploratory Data Analysis to examine existing hypotheses for rapid lake drainage. This involves testing for statistical differences between the rapidly and non-rapidly draining lake types, as well as examining associations between lake size and the potential controlling factors. This study shows that the two lake types are statistically indistinguishable for almost all factors investigated, except lake area. Thus, we are unable to recommend an empirically supported, deterministic alternative to the fracture area threshold parameter for modelling rapid lake drainage within existing surface-hydrology models of the GrIS. However, if improved remotely sensed datasets (e.g. ice-velocity maps, climate model outputs) were included in future research, it may be possible to detect the causes of rapid drainage.</jats:p

Modeled Subglacial Water Flow Routing Supports Localized Intrusive Heating as a Possible Cause of Basal Melting of Mars' South Polar Ice Cap

The discovery of a ~20 km wide area of bright subsurface radar reflections, interpreted as liquid water, beneath the Martian south polar layered deposits (SPLD) in data from the Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) instrument, and the discovery of two geologically recent potential eskers (landforms produced by subglacial melt) associated with viscous flow features in Martian mid-latitudes, has suggested recent basal melting of Martian ice deposits may be feasible, possibly due to locally elevated geothermal heating. Locations of terrestrial subglacial lakes and major drainage axes have been successfully predicted from subglacial hydraulic potential surfaces calculated from surface topography and ice thickness. Here, we use surface topography from the Mars Orbiter Laser Altimeter and SPLD bed elevations derived from MARSIS data to calculate the subglacial hydraulic potential surface beneath the SPLD and determine whether the observed high reflectance area coincides with predicted subglacial lake locations. Given the sensitivity of terrestrial predictions of lake locations to basal topography, we derive over 1000 perturbed topographies (using noise statistics from the MARSIS data) to infer the most likely locations of possible subglacial water bodies and drainage axes. Our results show that the high reflectance area does not coincide with any substantial predicted lake locations; three nearby lake locations are robustly predicted however. We interpret this result as suggesting that the high reflectance area (assuming the interpretation as liquid is correct) is most likely a hydraulically-isolated patch of liquid confined by the surrounding cold-based ice, rather than a topographically-constrained subglacial lake

Dual-satellite (Sentinel-2 and Landsat 8) remote sensing of supraglacial lakes in Greenland

© 2018 All rights reserved. Remote sensing is commonly used to monitor supraglacial lakes on the Greenland Ice Sheet (GrIS); however, most satellite records must trade off higher spatial resolution for higher temporal resolution (e.g. MODIS) or vice versa (e.g. Landsat). Here, we overcome this issue by developing and applying a dual-sensor method that can monitor changes to lake areas and volumes at high spatial resolution (10-30&thinsp;m) with a frequent revisit time ( ~ 3 days). We achieve this by mosaicking imagery from the Landsat 8 Operational Land Imager (OLI) with imagery from the recently launched Sentinel-2 Multispectral Instrument (MSI) for a ~ 12 000 km2area of West Greenland in the 2016 melt season. First, we validate a physically based method for calculating lake depths with Sentinel-2 by comparing measurements against those derived from the available contemporaneous Landsat 8 imagery; we find close correspondence between the two sets of values (R2Combining double low line 0.841; RMSE Combining double low line 0.555 m). This provides us with the methodological basis for automatically calculating lake areas, depths, and volumes from all available Landsat 8 and Sentinel-2 images. These automatic methods are incorporated into an algorithm for Fully Automated Supraglacial lake Tracking at Enhanced Resolution (FASTER). The FASTER algorithm produces time series showing lake evolution during the 2016 melt season, including automated rapid ( ≤ 4 day) lake-drainage identification. With the dual Sentinel-2-Landsat 8 record, we identify 184 rapidly draining lakes, many more than identified with either imagery collection alone (93 with Sentinel-2; 66 with Landsat 8), due to their inferior temporal resolution, or would be possible with MODIS, due to its omission of small lakes &lt; 0.125 km2. Finally, we estimate the water volumes drained into the GrIS during rapid-lake-drainage events and, by analysing downscaled regional climate-model (RACMO2.3p2) run-off data, the water quantity that enters the GrIS via the moulins opened by such events. We find that during the lake-drainage events alone, the water drained by small lakes ( &lt; 0.125 km2) is only 5.1 % of the total water volume drained by all lakes. However, considering the total water volume entering the GrIS after lake drainage, the moulins opened by small lakes deliver 61.5 % of the total water volume delivered via the moulins opened by large and small lakes; this is because there are more small lakes, allowing more moulins to open, and because small lakes are found at lower elevations than large lakes, where run-off is higher. These findings suggest that small lakes should be included in future remote-sensing and modelling work.NER

Spatial, seasonal and interannual variability of supraglacial ponds in the Langtang Valley of Nepal, 1999–2013

Supraglacial ponds play a key role in absorbing atmospheric energy and directing it to the ice of debris-covered glaciers, but the spatial and temporal distribution of these features is not well documented. We analyse 172 Landsat TM/ETM+ scenes for the period 1999–2013 to identify thawed supraglacial ponds for the debris-covered tongues of five glaciers in the Langtang Valley of Nepal. We apply an advanced atmospheric correction routine (Landcor/6S) and use band ratio and image morphological techniques to identify ponds and validate our results with 2.5 m Cartosat-1 observations. We then characterize the spatial, seasonal and interannual patterns of ponds. We find high variability in pond incidence between glaciers (May–October means of 0.08–1.69% of debris area), with ponds most frequent in zones of low surface gradient and velocity. The ponds show pronounced seasonality, appearing in the pre-monsoon as snow melts, peaking at the monsoon onset at 2% of debris-covered area, then declining in the post-monsoon as ponds drain or freeze. Ponds are highly recurrent and persistent, with 40.5% of pond locations occurring for multiple years. Rather than a trend in pond cover over the study period, we find high interannual variability for each glacier after controlling for seasonality

Past water flow beneath Pine Island and Thwaites glaciers, West Antarctica

Abstract. Outburst floods from subglacial lakes beneath the Antarctic Ice Sheet modulate ice-flow velocities over periods of months to years. Although subglacial lake drainage events have been observed from satellite-altimetric data, little is known about their role in the long-term evolution of ice-sheet basal hydrology. Here, we systematically map and model past water flow through an extensive area containing over 1000 subglacial channels and 19 former lake basins exposed on over 19 000 km2 of seafloor by the retreat of Pine Island and Thwaites glaciers, West Antarctica. At 507 m wide and 43 m deep on average, the channels offshore of present-day Pine Island and Thwaites glaciers are approximately twice as deep, 3 times as wide, and cover an area over 400 times larger than the terrestrial meltwater channels comprising the Labyrinth in the Antarctic Dry Valleys. The channels incised into bedrock offshore of contemporary Pine Island and Thwaites glaciers would have been capable of accommodating discharges of up to 8.8×106 m3 s−1. We suggest that the channels were formed by episodic discharges from subglacial lakes trapped during ice-sheet advance and retreat over multiple glacial periods. Our results document the widespread influence of episodic subglacial drainage events during past glacial periods, in particular beneath large ice streams similar to those that continue to dominate contemporary ice-sheet discharge. UK Natural Environment Research Council’s iSTAR programme (grant nos. NE/J005703/1, NE/J005746/1, and NE/J005770/1). James D. Kirkham: Debenham Scholarship from the Scott Polar Research Institute, University of Cambridge, and a UK Natural Environment Research Council Ph.D. studentship awarded through the Cambridge Earth System Science Doctoral Training Partnership (grant no. NE/L002507/1

Recent Basal Melting of a Mid-Latitude Glacier on Mars

Evidence for past basal melting of young (late Amazonian), debris-covered glaciers in Mars’ mid-latitudes is extremely rare. Thus, it is widely thought that these viscous flow features (VFFs) have been perennially frozen to their beds. We identify an instance of recent, localized wet-based mid-latitude glaciation, evidenced by a candidate esker emerging from a VFF in a tectonic rift in Tempe Terra. Eskers are sedimentary ridges deposited in ice-walled meltwater conduits and are indicative of glacial melting. We compare the candidate esker to terrestrial analogues, present a geomorphic map of landforms in the rift, and develop a landsystem model to explain their formation. We propose that the candidate esker formed during a transient phase of wet-based glaciation. We then consider the similarity between the geologic setting of the new candidate esker and that of the only other candidate esker to be identified in association with an existing mid-latitude VFF; both are within tectonic graben/rifts proximal to volcanic provinces. Finally, we calculate potential basal temperatures for a range of VFF thicknesses, driving stresses, mean annual surface temperatures, and geothermal heat fluxes, which unlike previous studies, include the possible role of internal strain heating. Strain heating can form an important additional heat source, especially in flow convergence zones, or where ice is warmer due to elevated surface temperatures or geothermal heat flux. Elevated geothermal heat flux within rifts, perhaps combined with locally-elevated strain heating, may have permitted wet-based glaciation during the late Amazonian, when cold climates precluded more extensive wet-based glaciation on Mars

An investigation of the influence of supraglacial debris on glacier-hydrology

Abstract. The influence of supraglacial debris on the rate and spatial distribution of glacier surface melt is well established, but its potential impact on the structure and evolution of the drainage system of extensively debris-covered glaciers has not been previously investigated. Forty-eight dye injections were conducted on Miage Glacier, Italian Alps, throughout the 2010 and 2011 ablation seasons. An efficient conduit system emanates from moulins in the mid-part of the glacier, which are downstream of a high melt area of dirty ice and patchy debris. High melt rates and runoff concentration by intermoraine troughs encourages the early-season development of a channelized system downstream of this area. Conversely, the drainage system beneath the continuously debris-covered lower ablation area is generally inefficient, with multi-peaked traces suggesting a distributed network, which likely feeds into the conduit system fed by the upglacier moulins. Drainage efficiency from the debris-covered area increased over the season but trace flow velocity remained lower than from the upper glacier moulins. Low and less-peaked melt inputs combined with the hummocky topography of the debris-covered area inhibits the formation of an efficient drainage network. These findings are relevant to regions with extensive glacial debris cover and where debris cover is expanding.</jats:p