Ice Cover, Subglacial Landscape, and Estimation of Bottom Melting of Mac. Robertson, Princess Elizabeth, Wilhelm II, and Western Queen Mary Lands, East Antarctica

This study demonstrates the results of Russian airborne radio-echo sounding (RES) investigations and also seismic reflection soundings carried out in 1971–2020 over a vast area of coastal part of East Antarctica. It is the first comprehensive summary mapping of these data. Field research, equipment, errors of initial RES data, and methods of gridding are discussed. Ice thickness, ice base elevation, and bedrock topography are presented. The ice thickness across the research area varies from a few meters to 3620 m, and is greatest in the local subglacial depressions. The average thickness is about 1220 m. The total volume of the ice is about 710,500 km3. The bedrock heights vary from 2860 m below sea level in the ocean bathyal zone to 2040 m above sea level in the Grove Mountains area (4900 m relief). The main directions of the bedrock orographic forms are concentrated mostly in three intervals: 345∘–30∘, 45∘–70∘, and 70∘–100∘. The bottom melting rate was estimated on the basis of the simple Zotikov model. Total annual melting under the study area is about 0.633 cubic meters. The total annual melting in the study area is approximately 1.5 mm/yr.

Basal melt of the southern Filchner Ice Shelf, Antarctica

Basal melt of ice shelves is a key factor governing discharge of ice from the Antarctic Ice Sheet as a result of its effects on buttressing. Here, we use radio echo sounding to determine the spatial variability of the basal melt rate of the southern Filchner Ice Shelf, Antarctica along the inflow of Support Force Glacier. We find moderate melt rates with a maximum of 1.13 m a−1 about 50 km downstream of the grounding line. The variability of the melt rates over distances of a few kilometres is low (all but one < 0.15 m a−1 at < 2 km distance), indicating that measurements on coarse observational grids are able to yield a representative melt rate distribution. A comparison with remote sensing based melt rates revealed that, for the study area, large differences were due to inaccuracies in the estimation of vertical strain rates from remote sensing velocity fields. These inaccuracies can be overcome by using modern velocity fields.

Development of a subglacial lake monitored with radio-echo sounding: case study from the eastern Skaftá cauldron in the Vatnajökull ice cap, Iceland

We present repeated radio-echo sounding (RES, 5 MHz) on a profile grid over the eastern Skaftá cauldron (ESC) in Vatnajökull ice cap, Iceland. The ESC is a ∼ 3 km wide and 50–150 m deep ice cauldron created and maintained by subglacial geothermal activity of ∼ 1 GW. Beneath the cauldron and 200–400 m thick ice, water accumulates in a subglacial lake and is released semi-regularly in jökulhlaups. The RES record consists of annual surveys conducted at the beginning of every summer during the period 2014–2020. Comparison of the RES surveys reveals variable lake area (0.5–4.1 km2) and enables traced reflections from the lake roof to be distinguished from bedrock reflections. This allows construction of a digital elevation model (DEM) of the bedrock in the area, further constrained by two borehole measurements at the cauldron centre. It also allows creation of lake thickness maps and an estimate of lake volume at the time of each survey, which we compare with lowering patterns and released water volumes obtained from pre- and post-jökulhlaup surface DEMs. The estimated lake volume was 250 GL (gigalitres = 106 m3) in June 2015, but 320 ± 20 GL drained from the ESC in October 2015. In June 2018, RES profiles revealed a lake volume of 185 GL, while 220 ± 30 GL were released in a jökulhlaup in August 2018. Considering the water accumulation over the periods between RES surveys and jökulhlaups, this indicates 10 %–20 % uncertainty in the RES-derived volumes at times when significant jökulhlaups may be expected.

AN AUTOMATIC EXTRACTION METHOD FOR ANTARCTIC SUBGLACIAL LAKE BASED ON RADIO ECHO SOUNDING DATA

Detection of subglacial lakes and interpretation their hydrological connectivity is of great importance to understanding the mass balance of Antarctic ice sheet. Over the past five decades, a large number of Radio Echo Sounding (RES) data has been collected in Antarctica. However, the identification of subglacial lakes based on RES data mainly relies on visual interpretation due to the lack of quantitative indicators for subglacial lakes distinguishing. To solve this problem, an automatic subglacial lakes detection method based on the reflective characteristics of B-scan echogram is proposed in this paper. Firstly, the echo strength of the original radar echogram is corrected for the dielectric attenuation and geometric spreading in the ice. Secondly, the original radar echogram is binarized to preserve the bright subglacial lake region, and the thickness of bright pixels is measured along the direction of A-scan. Thirdly, the thickness, the variance of the thickness in the neighbourhood and the corrected echo strength are used to obtain the response value for the subglacial lake. Finally, a threshold for the response value is determined to detect subglacial lakes. It is found that the proposed method can determine the location of the subglacial lake in Antarctica's Gamburtsev Province (AGAP) region with high accuracy.

Development of a subglacial lake monitored with radio-echo sounding: Case study from the Eastern Skaftá Cauldron in the Vatnajökull ice cap, Iceland

We present repeated radio-echo sounding (RES, 5 MHz) on a profile grid over the Eastern Skaftá Cauldron (ESC) in Vatnajökull ice cap, Iceland. The ESC is ~3 km wide and 50–150 m deep ice cauldron created and maintained by subglacial geothermal activity of ~1 GW. Beneath the cauldron and 200–400 m thick ice, water accumulates in a lake and is released semi-regularly in jökulhlaups. The RES record consists of annual surveys with 200–400 m between profiles in early summers of 2014–2020. Comparison of the RES surveys (2D migrated profiles) reveals variable lake area (0.5–4.1 km2) and enables traced reflections from the lake roof to be distinguished from bedrock reflections. This allows construction of a digital elevation model (DEM) of the bedrock in the area, further constrained by two borehole measurements at the cauldron centre. It also allows creation of lake thickness maps and an estimate of lake volume at the time of each survey, which we compare with lowering patterns and released water volumes obtained from surface DEMs obtained before and after jökulhlaups. The estimated lake volume is 250 × 106 m3 in June 2015 but 320 ± 20 × 106 m3 drained from the cauldron in October 2015. In June 2018, RES profiles reveal a lake volume of 185 × 106 m3 while 220 ± 30 × 106 m3 was released in a jökulhlaup in August 2018. Considering the water accumulation over the periods between RES surveys and jökulhlaups, this indicates 10–20 % uncertainty in the RES-derived volumes at times when significant jökulhlaups may be expected.

Radio-echo sounding and waveform modeling reveal abundant marine ice in former rifts and basal crevasses within Crary Ice Rise, Antarctica

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

Measuring and inferring the ice thickness distribution of four glaciers in the Tien Shan, Kyrgyzstan

Glaciers in the Tien Shan mountains contribute considerably to the fresh water used for irrigation, households and energy supply in the dry lowland areas of Kyrgyzstan and its neighbouring countries. To date, reconstructions of the current ice volume and ice thickness distribution remain scarce, and accurate data are largely lacking at the local scale. Here, we present a detailed ice thickness distribution of Ashu-Tor, Bordu, Golubin and Kara-Batkak glaciers derived from radio-echo sounding measurements and modelling. All the ice thickness measurements are used to calibrate three individual models to estimate the ice thickness in inaccessible areas. A cross-validation between modelled and measured ice thickness for a subset of the data is performed to attribute a weight to every model and to assemble a final composite ice thickness distribution for every glacier. Results reveal the thickest ice on Ashu-Tor glacier with values up to 201 ± 12 m. The ice thickness measurements and distributions are also compared with estimates composed without the use of in situ data. These estimates approach the total ice volume well, but local ice thicknesses vary substantially.

Subglacial topography and ice flux along the English Coast of Palmer Land, Antarctic Peninsula

Recent satellite data have revealed widespread grounding line retreat, glacier thinning, and associated mass loss along the Bellingshausen Sea sector, leading to increased concern for the stability of this region of Antarctica. While satellites have greatly improved our understanding of surface conditions, a lack of radio-echo sounding (RES) data in this region has restricted our analysis of subglacial topography, ice thickness, and ice flux. In this paper we analyse 3000 km of 150 MHz airborne RES data collected using the PASIN2 radar system (flown at 3–5 km line spacing) to investigate the subglacial controls on ice flow near the grounding lines of Ers, Envisat, Cryosat, Grace, Sentinel, Lidke, and Landsat ice streams as well as Hall and Nikitin glaciers. We find that each outlet is topographically controlled, and when ice thickness is combined with surface velocity data from MEaSUREs (Mouginot et al., 2019a), these outlets are found to discharge over 39.25 ± 0.79 Gt a−1 of ice to floating ice shelves and the Southern Ocean. Our RES measurements reveal that outlet flows are grounded more than 300 m below sea level and that there is limited topographic support for inland grounding line re-stabilization in a future retreating scenario, with several ice stream beds dipping inland at ∼ 5∘ km−1. These data reinforce the importance of accurate bed topography to model and understand the controls on inland ice flow and grounding line position as well as overall mass balance and sea level change estimates. RES data described in this paper are available through the UK Polar Data Centre: https://doi.org/10.5285/E07D62BF-D58C-4187-A019-59BE998939CC (Corr and Robinson, 2020).

Subglacial lakes and hydrology across the Ellsworth Subglacial Highlands, West Antarctica

Subglacial water plays an important role in ice sheet dynamics and stability. Subglacial lakes are often located at the onset of ice streams and have been hypothesised to enhance ice flow downstream by lubricating the ice–bed interface. The most recent subglacial-lake inventory of Antarctica mapped nearly 400 lakes, of which ∼ 14 % are found in West Antarctica. Despite the potential importance of subglacial water for ice dynamics, there is a lack of detailed subglacial-water characterisation in West Antarctica. Using radio-echo sounding data, we analyse the ice–bed interface to detect subglacial lakes. We report 33 previously uncharted subglacial lakes and present a systematic analysis of their physical properties. This represents a ∼ 40 % increase in subglacial lakes in West Antarctica. Additionally, a new digital elevation model of basal topography of the Ellsworth Subglacial Highlands was built and used to create a hydropotential model to simulate the subglacial hydrological network. This allows us to characterise basal hydrology, determine subglacial water catchments and assess their connectivity. We show that the simulated subglacial hydrological catchments of the Rutford Ice Stream, Pine Island Glacier and Thwaites Glacier do not correspond to their ice surface catchments.

Bed topography of Princess Elizabeth Land in East Antarctica

We present a topographic digital elevation model (DEM) for Princess Elizabeth Land (PEL), East Antarctica. The DEM covers an area of ∼900 000 km2 and was built from radio-echo sounding data collected during four campaigns since 2015. Previously, to generate the Bedmap2 topographic product, PEL's bed was characterized from low-resolution satellite gravity data across an otherwise large (>200 km wide) data-free zone. We use the mass conservation (MC) method to produce an ice thickness grid across faster flowing (>30 m yr−1) regions of the ice sheet and streamline diffusion in slower flowing areas. The resulting ice thickness model is integrated with an ice surface model to build the bed DEM. Together with BedMachine Antarctica and Bedmap2, this new bed DEM completes the first-order measurement of subglacial continental Antarctica – an international mission that began around 70 years ago. The ice thickness data and bed DEMs of PEL (resolved horizontally at 500 m relative to ice surface elevations obtained from the Reference Elevation Model of Antarctica – REMA) are accessible from https://doi.org/10.5281/zenodo.4023343 (Cui et al., 2020a) and https://doi.org/10.5281/zenodo.4023393 (Cui et al., 2020b).