ALBERT

Description: Lakes beneath the Antarctic Ice Sheet are known to decrease traction at the ice base and therefore can have a great impact on ice dynamics. However, the total extent of Antarctic subglacial lakes is still unknown. We address this issue by combining modeling and remote-sensing strategies to predict potential lake locations using the general hydraulic potential equation. We are able to reproduce the majority of known lakes, as well as predict the existence of many new and so far undetected potential lakes. To validate our predictions, we analyzed ice-penetrating radar profiles from radio-echo sounding flights acquired over 1994-2013 in Dronning Maud Land, East Antarctica, and this led to the identification of 31 new subglacial lakes. Based on these findings, we estimate the total number of Antarctic subglacial lakes to be similar to 1300, a factor of three higher than the total number of lakes discovered to date. We estimate that only similar to 30% of all Antarctic subglacial lakes and similar to 65% of the total estimated lake-covered area have been discovered, and that lakes account for 0.6% of the Antarctic ice/bed interface.

Description: The internal stratigraphy of snow and ice as imaged by ground-penetrating radar may serve as a source of information on past accumulation. This study presents results from two ground-based radar surveys conducted in Greenland in 2007 and 2015, respectively. The first survey was conducted during the traverse from the ice-core station NGRIP (North Greenland Ice Core Project) to the ice-core station NEEM (North Greenland Eemian Ice Drilling). The second survey was carried out during the traverse from NEEM to the ice-core station EGRIP (East Greenland Ice Core Project) and then onwards to Summit Station. The total length of the radar profiles is 1427 km. From the radar data, we retrieve the large-scale spatial variation of the accumulation rates in the interior of the ice sheet. The accumulation rates range from 0.11 to 0.26 m a−1 ice equivalent with the lowest values found in the northeastern sector towards EGRIP. We find no evidence of temporal or spatial changes in accumulation rates when comparing the 150-year average accumulation rates with the 321-year average accumulation rates. Comparisons with regional climate models reveal that the models underestimate accumulation rates by up to 35% in northeastern Greenland. Our results serve as a robust baseline to detect present changes in either surface accumulation rates or patterns.

Description: Basal motion of ice sheets depends in part on the roughness and material properties of the subglacial bed and the occurrence of water. To date, basal motion represents one of the largest uncertainties in ice-flow models. It is that component of the total flow velocity that can change most rapidly and can, therefore, facilitate rapid variations in dynamic behaviour. In this study, we investigate the subglacial properties of the East Antarctic Ice Sheet by statistically analysing the roughness of the bed topography, inferred from radio-echo sounding measurements. We analyse two sets of roughness parameters, one derived in the spatial and the other in the spectral domain, with two roughness parameters each. This enables us to compare the suitability of the four roughness parameters to classify the subglacial landscapes below the ice sheet. We further investigate the relationship of the roughness parameters with observed surface flow velocity and modelled basal temperatures of the ice sheet. We find that one of the roughness parameters, the Hurst exponent derived in the spatial domain, coincides with the thermal condition at the base of the ice sheet for slow flow velocities and varies with flow velocity.

Description: The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO_2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the International Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry surrounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79° N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities. The data set is available in full and in regional subsets in NetCDF format from the PANGAEA database at https://doi.pangaea.de/10.1594/PANGAEA.856844.

Description: The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

Description: This paper presents a new compilation of gridded datasets for three-dimensional modelling of the Antarctic ice sheet. These are for surface elevation, ice thickness, bedrock elevation and accumulation rate as interpolated on a 281 × 281 mesh with 20 km spacing, and encompass all the ice sheet and surrounding continental shelf. Data sources include the Bamber digital-elevation model from ERS-1 radar-altimeter data, a redigitization of available ice-thickness data, the Giovinetto accumulation data, recent ice-thickness data from British and German expeditions as well as accumulation data from German and Norwegian expeditions. In particular, new data were incorporated for the Filchner-Ronne Ice Shelf and for Dronning Maud Land, Antarctica, arising from the EPICA pre-site survey. Special attention was devoted to matching the various data sources carefully, both among themselves and across the grounding line and below the ice shelves, to enable ice-sheet expansion and retreat in dynamic situations. As an application, the balance flow is calculated over the entire ice sheet using a two-dimensional finite-difference scheme and compared with a previous assessment. This brought to light the existence of ice-streaming features extending well inland. A detailed zoom over Dronning Maud Land exhibits the general flow characteristics of interest for locating a future deep-drilling site. As a by-product, an updated value of 26.4 × 106km3 was obtained for the total volume of the ice sheet and ice shelves, or equivalent to 61.1 m of global sea-level rise after removal of the ice sheet and subsequent oceanic invasion and isostatic rebound. The total accumulation over the grounded ice sheet, including the Antarctic Peninsula, is 1924 Gta−1, or between 5 and 20% higher than earlier estimates. Including all the ice shelves, the value is 2344 Gt a−1.

Description: During the summer of 2003, a ground-penetrating radar survey around the North Greenland Icecore Project (NorthGRIP) deep ice-core drilling site (75˚06’N, 42˚20’W; 2957ma.s.l.) was carried out using a shielded 250 MHz radar system. The drill site is located on an ice divide, roughly 300 km north-northwest of the summit of the Greenland ice sheet. More than 430 km of profiles were measured, covering a 10 km by 10 km area, with a grid centered on the drilling location, and eight profiles extending beyond this grid. Seven internal horizons within the upper 120 m of the ice sheet were continuously tracked, containing the last 400 years of accumulation history. Based on the age-depth and density-depth distribution of the deep core, the internal layers have been dated and the regional and temporal distribution of accumulation rate in the vicinity of NorthGRIP has been derived. The distribution of accumulation shows a relatively smoothly increasing trend from east to west from 145 kgm–2a–1 to 200 kg m–2 a -1 over a distance of 50 km across the ice divide. The general trend is overlain by small-scale variations on the order of 2.5 kgm–2a-1 km- 1 , i.e. around 1.5% of the accumulation mean. The temporal variations of the seven periods defined by the seven tracked isochrones are on the order of ± 4% of the mean of the last 400 years, i.e. at NorthGRIP ± 7 kg m–2 a-1. If the regional accumulation pattern has been stable for the last several thousand years during the Holocene, and ice flow has been comparable to today, advective effects along the particle trajectory upstream of NorthGRIP do not have a significant effect on the interpretation of climatically induced changes in accumulation rates derived from the deep ice core over the last 10 kyr.

Description: The East Antarctic Ice Sheet contains a wealth of information that can be extracted from its internal architecture such as distribution of age, past flow features and surface and basal properties. Airborne radar surveys can sample this stratigraphic archive across broad areas. Here, we identify and trace key horizons across several radar surveys to obtain the stratigraphic information. We transfer the age-depth scales from ice cores to intersecting radar data. We then propagate these age scales across the ice sheet using the high fidelity continuity of the radar horizons. In Dronning Maud Land, including Dome Fuji, we mapped isochrones with ages of 38 ka and 74 ka. In the central region of East Antarctica around Dome Concordia, Vostok, and Dome Argus, we use isochrone ages of 38 ka, 48 ka, 90 ka, and 161 ka. Taking together both regions, we provide isochrone depths traced along a combined profile length of more than 40,000 km and discuss uncertainties of the obtained stratigraphy, as well as factors important to consider for further expansion. This dataset the most extensive distribution of internal horizons in East Antarctica to date. The isochrone depths are available on PANGAEA: https://doi.pangaea.de/10.1594/PANGAEA.895528.

Description: This study aims to demonstrate that deep ice cores can be synchronized using internal horizons in the ice between the drill sites revealed by airborne radio-echo sounding (RES) over a distance of 〉1000km, despite significant variations in glaciological parameters, such as accumulation rate between the sites. In 2002/03 a profile between the Kohnen station and Dome Fuji deep ice-core drill sites, Antarctica, was completed using airborne RES. The survey reveals several continuous internal horizons in the RES section over a length of 1217 km. The layers allow direct comparison of the deep ice cores drilled at the two stations. In particular, the counterpart of a visible layer observed in the Kohnen station (EDML) ice core at 1054 m depth has been identified in the Dome Fuji ice core at 575 m depth using internal RES horizons. Thus the two ice cores can be synchronized, i.e. the ice at 1560 m depth (at the bottom of the 2003 EDML drilling) is ∼49ka old according to the Dome Fuji age/depth scale, using the traced internal layers presented in this study.

Description: We used internal ice layers from a radio-echo sounding profile between the Kohnen and Dome Fuji deep drilling sites to infer the spatio-temporal pattern of accumulation rate in this sector of Dronning Maud Land, East Antarctica. Continuous internal reflection horizons can be traced to about half of the ice thickness and have a maximum age of approximately 72.7 ka BP. To infer palaeo-accumulation rates from the dated layers, we derived the thinning functions from a flow calculation with a high-resolution higher-order model of Dronning Maud Land embedded into a three-dimensional thermomechanical model of the Antarctic ice sheet. The method takes into account complex ice-flow dynamics and advection effects that cannot be dealt with using traditional local approaches. We selected seven time intervals over which we determine the average accumulation rate and average surface temperature at the place and time of origin of the layer particles. Our results show lower accumulation rates along eastern parts of the profile for the late Holocene (0–5 ka BP) than are shown by existing maps, which had no surface control points. During the last glacial period we find a substantially lower accumulation rate than predicted by the usual approach linking palaeo-accumulation rates to the condensation temperature above the surface inversion layer. These findings were used to fine-tune the relation between accumulation rate and temperature.