ALBERT

Description: Digital elevation models of Antarctic bed topography are smoothed and interpolated onto low-resolution ( 〉 1 km) grids as current observed topography data are generally sparsely and unevenly sampled. This issue has potential implications for numerical simulations of ice-sheet dynamics, especially in regions prone to instability where detailed knowledge of the topography, including fine-scale roughness, is required. Here, we present a high-resolution (100 m) synthetic bed elevation terrain for Antarctica, encompassing the continent, continental shelf, and seas south of 60° S. Although not identically matching observations, the synthetic bed surface – denoted as HRES – preserves topographic roughness characteristics of airborne and ground-based ice-penetrating radar data measured by the ICECAP (Investigating the Cryospheric Evolution of the Central Antarctic Plate) consortium or used to create the Bedmap1 compilation. Broad-scale ( 〉 5 km resolution) features of the Antarctic landscape are incorporated using a low-pass filter of the Bedmap2 bed elevation data. HRES has applicability in high-resolution ice-sheet modelling studies, including investigations of the interaction between topography, ice-sheet dynamics, and hydrology, where processes are highly sensitive to bed elevations and fine-scale roughness. The data are available for download from the Australian Antarctic Data Centre (doi:10.4225/15/57464ADE22F50).

Description: We present a new bed elevation digital elevation model (DEM), with a 1 km spatial resolution, for the Weddell Sea sector of the West Antarctic Ice Sheet. The DEM has a total area of ~125,000 km2 covering the Institute, Möller and Foundation ice streams and the Bungenstock ice rise. In comparison with the Bedmap2 product, our DEM includes several new aerogeophysical datasets acquired by the Center for Remote Sensing of Ice Sheets (CReSIS) through the NASA Operation IceBridge (OIB) program in 2012, 2014 and 2016. We also update bed elevation information from the single largest existing dataset in the region, collected by the British Antarctic Survey (BAS) Polarimetric Airborne Survey Instrument (PASIN) in 2010-11, as BEDMAP2 included only relatively crude ice thickness measurements determined in the field for quality control purposes. This have resulted in the deep parts of the topography not being visible in the fieldwork non-SAR processed radargrams. While the gross form of the new DEM is similar to Bedmap2, there are some notable differences. For example, the position and size of a deep trough (~ 2 km below sea level) between the ice sheet interior and the grounding line of Foundation ice stream has been redefined. From the revised DEM, we are able to better derive the expected routing of basal water at the ice-bed interface, and by comparison with that calculated using Bedmap2 we are able to assess regions where hydraulic flow is sensitive to change. Given the sensitivity of this sector of the ice sheet to ocean-induced melting at the grounding line, especially in light of improved definition of the Foundation ice stream trough, our revised DEM will be of value to ice-sheet modelling in efforts to quantify future glaciological changes in the region, and therefore the potential impact on global sea level. The new 1 km bed elevation product of the Weddell Sea sector, West Antarctica can be found in the http://doi.org/10.5281/zenodo.1035488.

Description: Subglacial roughness can be determined at variety of length scales from radio-echo sounding (RES) data; either via statistical analysis of along-track topography, or inferred from basal radar scattering. Past studies have demonstrated that subglacial terrain exhibits self-affine (fractal) scaling behaviour, where vertical roughness has a power-law relationship with the horizontal length scale. A self-affine statistical framework, which enables a consistent integration of topographic roughness and radar scattering, has yet to be applied to RES. Here we do this for recent RES data from northern Greenland, and demonstrate that subglacial topography exhibits pronounced spatial variation in the Hurst (roughness power-law) exponent. A radar scattering model then enables us to explain how the Hurst exponent exerts strong topographic control upon radar scattering, which we map using the waveform abruptness (pulse peakiness) parameter. Notably, lower abruptness (associated with diffuse scattering) occurs for regions with a higher Hurst exponent, and higher abruptness (associated with specular reflections) occurs for regions with a lower Hurst exponent. Finally, we compare the RES-derived data with an independent prediction for the subglacial thermal state of northern Greenland. This analysis shows that the majority of predicted thawed regions do not have the specular RES scattering signature of deep subglacial lakes, and instead have a diffuse scattering signature.

Description: The subglacial environment of the Greenland Ice Sheet (GrIS) is poorly constrained both in its bulk properties, for example geology, the presence of sediment, and the presence of water, and interfacial conditions, such as roughness and bed rheology. There is, therefore, limited understanding of how spatially heterogeneous subglacial properties relate to ice-sheet motion. Here, via analysis of 2 decades of radio-echo sounding data, we present a new systematic analysis of subglacial roughness beneath the GrIS. We use two independent methods to quantify subglacial roughness: first, the variability in along-track topography – enabling an assessment of roughness anisotropy from pairs of orthogonal transects aligned perpendicular and parallel to ice flow and, second, from bed-echo scattering – enabling assessment of fine-scale bed characteristics. We establish the spatial distribution of subglacial roughness and quantify its relationship with ice flow speed and direction. Overall, the beds of fast-flowing regions are observed to be rougher than the slow-flowing interior. Topographic roughness exhibits an exponential scaling relationship with ice surface velocity parallel, but not perpendicular, to flow direction in fast-flowing regions, and the degree of anisotropy is correlated with ice surface speed. In many slow-flowing regions both roughness methods indicate spatially coherent regions of smooth beds, which, through combination with analyses of underlying geology, we conclude is likely due to the presence of a hard flat bed. Consequently, the study provides scope for a spatially variable hard- or soft-bed boundary constraint for ice-sheet models.

Description: Digital elevation models of Antarctic bed topography are heavily smoothed and interpolated onto low-resolution (〉 1 km) grids as our current observed topography data are generally sparsely and unevenly sampled. This issue has potential implications for numerical simulations of ice-sheet dynamics, especially in regions prone to instability where detailed knowledge of the topography, including fine-scale roughness, is required. Here, we present a high-resolution (100 m) synthetic bed elevation terrain for the whole Antarctic continent. The synthetic bed surface preserves topographic roughness characteristics of airborne and ground-based ice-penetrating radar data from the Bedmap1 compilation and the ICECAP consortium. Broad-scale features of the Antarctic landscape are incorporated using a low-pass filter of the Bedmap2 bed-elevation data. Although not intended as a substitute for Bedmap2, the simulated bed elevation terrain has applicability in high-resolution ice-sheet modelling studies, including investigations of the interaction between topography, ice-sheet dynamics, and hydrology, where processes are highly sensitive to bed elevations. The data are available for download at the Australian Antarctic Data Centre (doi:10.4225/15/57464ADE22F50).

Description: Ice cores provide some of the best-dated and most comprehensive proxy records, as they yield a vast and growing array of proxy indicators. Selecting a site for ice core drilling is nonetheless challenging, as the assessment of potential new sites needs to consider a variety of factors. Here, we demonstrate a systematic approach to site selection for a new East Antarctic high-resolution ice core record. Specifically, seven criteria are considered: (1) 2000-year-old ice at 300 m depth; (2) above 1000 m elevation; (3) a minimum accumulation rate of 250 mm years−1 IE (ice equivalent); (4) minimal surface reworking to preserve the deposited climate signal; (5) a site with minimal displacement or elevation change in ice at 300 m depth; (6) a strong teleconnection to midlatitude climate; and (7) an appropriately complementary relationship to the existing Law Dome record (a high-resolution record in East Antarctica). Once assessment of these physical characteristics identified promising regions, logistical considerations (for site access and ice core retrieval) were briefly considered. We use Antarctic surface mass balance syntheses, along with ground-truthing of satellite data by airborne radar surveys to produce all-of-Antarctica maps of surface roughness, age at specified depth, elevation and displacement change, and surface air temperature correlations to pinpoint promising locations. We also use the European Centre for Medium-Range Weather Forecast ERA 20th Century reanalysis (ERA-20C) to ensure that a site complementary to the Law Dome record is selected. We find three promising sites in the Indian Ocean sector of East Antarctica in the coastal zone from Enderby Land to the Ingrid Christensen Coast (50–100° E). Although we focus on East Antarctica for a new ice core site, the methodology is more generally applicable, and we include key parameters for all of Antarctica which may be useful for ice core site selection elsewhere and/or for other purposes.

Description: Ice sheets provide exceptional archives of past changes in polar climate, regional environment and global atmospheric composition. The oldest dated deep ice core drilled in Antarctica has been retrieved at EPICA Dome C (EDC), reaching ~ 800,000 years. Obtaining an older paleoclimatic record from Antarctica is one of the greatest challenges of the ice core community. Here, we use internal isochrones, identified from airborne radar coupled to ice-flow modelling to estimate the age of basal ice along transects in the Dome C area. Three glaciological properties are inverted from isochrones: surface accumulation rate; geothermal flux; and the exponent of the Lliboutry velocity profile. We find that old ice (〉 1 Myr, 1 million years) likely exists in two regions: one ~ 40 km south-west of Dome C along the ice divide to Vostok, close to a secondary dome that we name "Little Dome C" (LDC); and a second region named "North Patch" (NP) located 10–30 km north-east of Dome C, in a region where the geothermal flux is apparently relatively low. Our work demonstrates the value of combining radar observations with ice flow modelling to accurately represent the true nature of ice flow, and the formation of ice-sheet architecture, in the centre of large ice sheets.

Description: Bed topography is a critical boundary for the numerical modelling of ice sheets and ice–ocean interactions. A persistent issue with existing topography products for the bed of the Greenland Ice Sheet and surrounding sea floor is the poor representation of coastal bathymetry, especially in regions of floating ice and near the grounding line. Sparse data coverage, and the resultant coarse resolution at the ice–ocean boundary, poses issues in our ability to model ice flow advance and retreat from the present position. In addition, as fjord bathymetry is known to exert strong control on ocean circulation and ice–ocean forcing, the lack of bed data leads to an inability to model these processes adequately. Since the release of the last complete Greenland bed topography–bathymetry product, new observational bathymetry data have become available. These data can be used to constrain bathymetry, but many fjords remain completely unsampled and therefore poorly resolved. Here, as part of the development of the next generation of Greenland bed topography products, we present a new method for constraining the bathymetry of fjord systems in regions where data coverage is sparse. For these cases, we generate synthetic fjord geometries using a method conditioned by surveys of terrestrial glacial valleys as well as existing sinuous feature interpolation schemes. Our approach enables the capture of the general bathymetry profile of a fjord in north-west Greenland close to Cape York, when compared to observational data. We validate our synthetic approach by demonstrating reduced overestimation of depths compared to past attempts to constrain fjord bathymetry. We also present an analysis of the spectral characteristics of fjord centrelines using recently acquired bathymetric observations, demonstrating how a stochastic model of fjord bathymetry could be parameterised and used to create different realisations.

Description: Subglacial roughness can be determined at a variety of length scales from radio-echo sounding (RES) data either via statistical analysis of topography or inferred from basal radar scattering. Past studies have demonstrated that subglacial terrain exhibits self-affine (power law) roughness scaling behaviour, but existing radar scattering models do not take this into account. Here, using RES data from northern Greenland, we introduce a self-affine statistical framework that enables a consistent integration of topographic-scale roughness with the electromagnetic theory of radar scattering. We demonstrate that the degree of radar scattering, quantified using the waveform abruptness (pulse peakiness), is topographically controlled by the Hurst (roughness power law) exponent. Notably, specular bed reflections are associated with a lower Hurst exponent, with diffuse scattering associated with a higher Hurst exponent. Abrupt waveforms (specular reflections) have previously been used as a RES diagnostic for basal water, and to test this assumption we compare our radar scattering map with a recent prediction for the basal thermal state. We demonstrate that the majority of thawed regions (above pressure melting point) exhibit a diffuse scattering signature, which is in contradiction to the prior approach. Self-affine statistics provide a generalised model for subglacial terrain and can improve our understanding of the relationship between basal properties and ice-sheet dynamics.

Description: Ice sheets provide exceptional archives of past changes in polar climate, regional environment and global atmospheric composition. The oldest dated deep ice core drilled in Antarctica has been retrieved at EPICA Dome C (EDC), reaching ∼ 800 000 years. Obtaining an older paleoclimatic record from Antarctica is one of the greatest challenges of the ice core community. Here, we use internal isochrones, identified from airborne radar coupled to ice-flow modelling to estimate the age of basal ice along transects in the Dome C area. Three glaciological properties are inferred from isochrones: surface accumulation rate, geothermal flux and the exponent of the Lliboutry velocity profile. We find that old ice (〉 1.5 Myr, 1.5 million years) likely exists in two regions: one ∼ 40 km south-west of Dome C along the ice divide to Vostok, close to a secondary dome that we name Little Dome C (LDC), and a second region named North Patch (NP) located 10–30 km north-east of Dome C, in a region where the geothermal flux is apparently relatively low. Our work demonstrates the value of combining radar observations with ice flow modelling to accurately represent the true nature of ice flow, and understand the formation of ice-sheet architecture, in the centre of large ice sheets.