ALBERT

Description: We present the results of mapping the limit of the tidal flexure (point F) and hydrostatic equilibrium (point H) of the grounding zone of Antarctic ice shelves from CryoSat-2 standard and swath elevation data. Overall we were able to map 31 % of the grounding zone of the Antarctic floating ice shelves and outlet glaciers. We obtain near-complete coverage of the Filchner–Ronne Ice Shelf. Here we manage to map areas of Support Force Glacier and the Doake Ice Rumples, which have previously only been mapped using break-in-slope methods. Over the Ross Ice Shelf, Dronning Maud Land and the Antarctic Peninsula, we obtained partial coverage, and we could not map a continuous grounding zone for the Amery Ice Shelf and the Amundsen Sea sector. Tidal amplitude and distance south (i.e. across-track spacing) are controlling factors in the quality of the coverage and performance of the approach. The location of the point F agrees well with previous observations that used differential satellite radar interferometry (DInSAR) and ICESat-1, with an average landward bias of 0.1 and 0.6 km and standard deviation of 1.1 and 1.5 km for DInSAR and ICESat measurements, respectively. We also compared the results directly with DInSAR interferograms from the Sentinel-1 satellites, acquired over the Evans Ice Stream and the Carlson Inlet (Ronne Ice Shelf), and found good agreement with the mapped points F and H. We also present the results of the spatial distribution of the grounding zone width (the distance between points F and H) and used a simple elastic beam model, along with ice thickness calculations, to calculate an effective Young modulus of ice of E=1.4±0.9 GPa.

Description: A major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry, and to a lesser extent satellite altimetry, is the poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA). In the past decade, much progress has been made in consistently modelling the ice sheet and solid Earth interactions; however, forward-modelling solutions of GIA in Antarctica remain uncertain due to the sparsity of constraints on the ice sheet evolution, as well as the Earth's rheological properties. An alternative approach towards estimating GIA is the joint inversion of multiple satellite data – namely, satellite gravimetry, satellite altimetry and GPS, which reflect, with different sensitivities, trends of recent glacial changes and GIA. Crucial to the success of this approach is the accuracy of the space-geodetic data sets. Here, we present reprocessed rates of surface-ice elevation change (Envisat/ICESat; 2003–2009), gravity field change (GRACE; 2003–2009) and bedrock uplift (GPS; 1995–2013.7). The data analysis is complemented by the forward-modelling of viscoelastic response functions to disc load forcing, allowing us to relate GIA-induced surface displacements with gravity changes for different rheological parameters of the solid Earth. The data and modelling results presented here are available in the Pangea archive; https://doi.pangaea.de/10.1594/PANGAEA.875745. The data sets are the input streams for the joint inversion estimate of present-day ice-mass change and GIA, focusing on Antarctica. However, the methods, code and data provided in this paper are applicable to solve other problems, such as volume balances of the Antarctic ice sheet, or to other geographical regions, in the case of the viscoelastic response functions. This paper presents the first of two contributions summarizing the work carried out within a European Space Agency funded study, REGINA.

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: Runoff from the Greenland Ice Sheet (GrIS) has increased in recent years due largely to declining albedo and enhanced surface melting. Some of the largest declines in GrIS albedo have occurred in the ablation zone of the south-west sector and are associated with the development of 'dark' ice surfaces. Field observations at local scales reveal that a variety of light-absorbing impurities (LAIs) can be present on the surface, ranging from inorganic particulates, to cryoconite materials and ice algae. Meanwhile, satellite observations show that the areal extent of dark ice has varied significantly between recent successive melt seasons. However, the processes that drive such large inter-annual variability in dark ice extent remain essentially unconstrained. At present we are therefore unable to project how the albedo of bare-ice sectors of the GrIS will evolve, causing uncertainty in the projected sea level contribution from the GrIS over the coming decades. Here we use MODIS satellite imagery to examine dark ice dynamics on the south-west GrIS each year from 2000 to 2016. We quantify dark ice in terms of its annual extent, duration, intensity and timing of first appearance. Not only does dark ice extent vary significantly between years, but so too does its duration (from 0 % to 〉 80 % of June–July–August, JJA), intensity and the timing of its first appearance. Comparison of dark ice dynamics with potential meteorological drivers from the regional climate model MAR reveals that the JJA sensible heat flux, the number of positive minimum-air-temperature days and the timing of bare ice appearance are significant inter-annual synoptic controls. We use these findings to identify the surface processes which are most likely to explain recent dark ice dynamics.We suggest that whilst the spatial distribution of dark ice is best explained by outcropping of particulates from ablating ice, these particulates alone do not drive dark ice dynamics. Instead, they may enable the growth of pigmented ice algal assemblages which cause visible surface darkening, but only when the climatological pre-requisites of liquid meltwater presence and sufficient photosynthetically-active radiation fluxes are met. Further field studies are required to fully constrain the processes by which ice algae growth proceeds and the apparent dependency of algae growth on melt-out particulates.

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: Runoff from the Greenland Ice Sheet (GrIS) has increased in recent years due largely to changes in atmospheric circulation and atmospheric warming. Albedo reductions resulting from these changes have amplified surface melting. Some of the largest declines in GrIS albedo have occurred in the ablation zone of the south-west sector and are associated with the development of dark ice surfaces. Field observations at local scales reveal that a variety of light-absorbing impurities (LAIs) can be present on the surface, ranging from inorganic particulates to cryoconite materials and ice algae. Meanwhile, satellite observations show that the areal extent of dark ice has varied significantly between recent successive melt seasons. However, the processes that drive such large interannual variability in dark ice extent remain essentially unconstrained. At present we are therefore unable to project how the albedo of bare ice sectors of the GrIS will evolve in the future, causing uncertainty in the projected sea level contribution from the GrIS over the coming decades. Here we use MODIS satellite imagery to examine dark ice dynamics on the south-west GrIS each year from 2000 to 2016. We quantify dark ice in terms of its annual extent, duration, intensity and timing of first appearance. Not only does dark ice extent vary significantly between years but so too does its duration (from 0 to 〉 80 % of June–July–August, JJA), intensity and the timing of its first appearance. Comparison of dark ice dynamics with potential meteorological drivers from the regional climate model MAR reveals that the JJA sensible heat flux, the number of positive minimum-air-temperature days and the timing of bare ice appearance are significant interannual synoptic controls. We use these findings to identify the surface processes which are most likely to explain recent dark ice dynamics. We suggest that whilst the spatial distribution of dark ice is best explained by outcropping of particulates from ablating ice, these particulates alone do not drive dark ice dynamics. Instead, they may enable the growth of pigmented ice algal assemblages which cause visible surface darkening, but only when the climatological prerequisites of liquid meltwater presence and sufficient photosynthetically active radiation fluxes are met. Further field studies are required to fully constrain the processes by which ice algae growth proceeds and the apparent dependency of algae growth on melt-out particulates.

Description: The subglacial environment of the Greenland Ice Sheet (GrIS) is poorly constrained, both in its bulk properties, for example geology, presence of sediment, and 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 two decades worth 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 of 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 bed, 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 bed/soft bed boundary constraint for ice-sheet models.

Description: There is widespread, but often indirect, evidence that a significant fraction of the bed beneath the Greenland Ice Sheet is thawed (at or above the pressure melting point for ice). This includes the beds of major outlet glaciers and their tributaries and a large area around the NorthGRIP borehole in the ice-sheet interior. The ice-sheet-scale distribution of basal water is, however, poorly constrained by existing observations. In principle, airborne radio-echo sounding (RES) enables the detection of basal water from bed-echo reflectivity, but unambiguous mapping is limited by uncertainty in signal attenuation within the ice. Here we introduce a new, RES diagnostic for basal water that is associated with wet–dry transitions in bed material: bed-echo reflectivity variability. This technique acts as a form of edge detector and is a sufficient, but not necessary, criteria for basal water. However, the technique has the advantage of being attenuation insensitive and suited to combined analysis of over a decade of Operation IceBridge survey data.The basal water predictions are compared with existing analyses of the basal thermal state (frozen and thawed beds) and geothermal heat flux. In addition to the outlet glaciers, we demonstrate widespread water storage in the northern and eastern interior. Notably, we observe a quasilinear corridor of basal water extending from NorthGRIP to Petermann Glacier that spatially correlates with elevated heat flux predicted by a recent magnetic model. Finally, with a general aim to stimulate regional- and process-specific investigations, the basal water predictions are compared with bed topography, subglacial flow paths and ice-sheet motion. The basal water distribution, and its relationship with the thermal state, provides a new constraint for numerical models.