ALBERT

Description: The study presented here is focused on the assessment of surface elevations derived from CryoSat-2 SARIn level 1b data over the Austfonna ice cap, Svalbard, in 2016. The processing chain that must be applied to the CryoSat-2 waveforms to derive heights is non-trivial, and consists of multiple steps, all requiring subjective choices of methods such as the choice of retracker, geo-relocation, and outlier rejection. Here, we compare six CryoSat-2 level-2 type data sets of surface elevations derived using different SARIn processing chains. These data sets are validated against surface elevation data collected from an airborne laser scanner, during a dedicated CryoSat validation experiment field campaign carried out in April 2016. The flight pattern of the airborne campaign was designed so that elevations were measured in a grid pattern rather than along single lines, as has previously been the standard procedure. The flight grid pattern was chosen to optimize the comparison with the CryoSat-2 SARIn elevation data, the location of which can deviate from nadir by several kilometers due to topography within the satellite footprint. The processing chains behind the six data sets include different outlier/error rejection approaches, and do not produce the same number of data points in our region of interest. To make a consistent analysis, we provide statistics from the validation of both the full data sets from each processing chain, and on only those data that all the six data sets provide a geo-located elevation estimate for. We find that the CryoSat-2 data sets that agree best with the validation data are those derived from dedicated land ice processing schemes. This study may serve as a benchmark for future CryoSat-2 retracker developments, and the evaluation software and data set are made publicly available.

Description: Mass balances of individual glaciers on ice sheets have been previously reported by forming a mass budget of discharged ice and modelled ice sheet surface mass balance or a complementary method which measures volume changes over the glaciated area that are subsequently converted to glacier mass change. On ice sheets, volume changes have been measured predominantly with radar and laser altimeters but InSAR DEM differencing has also been applied on smaller ice bodies. Here, we report for the first time on the synergistic use of volumetric measurements from the CryoSat-2 radar altimetry mission together with TanDEM-X DEM differencing and calculate the mass balance of the two major outlet glaciers of the Northeast Greenland Ice Stream: Zachariæ Isstrøm and Nioghalvfjerdsfjorden (79North). The glaciers lost 3.59 ± 1.15 G t a − 1 and 1.01 ± 0.95 G t a − 1 , respectively, between January 2011 and January 2014. Additionally, there has been substantial sub-aqueous mass loss on Zachariæ Isstrøm of more than 11 G t a − 1 . We attribute the mass changes on both glaciers to dynamic downwasting. The presented methodology now permits using TanDEM-X bistatic InSAR data in the context of geodetic mass balance investigations for large ice sheet outlet glaciers. In the future, this will allow monitoring the mass changes of dynamic outlet glaciers with high spatial resolution while the superior vertical accuracy of CryoSat-2 can be used for the vast accumulation zones in the ice sheet interior.

Description: Time-variable gravity measurements from the Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow On (GRACE-FO) missions and satellite altimetry measurements from CryoSat-2 enable independent mass balance estimates of the Earth’s glaciers and ice sheets. Both approaches vary in terms of their retrieval principles and signal-to-noise characteristics. GRACE/GRACE-FO recovers the gravity disturbance caused by changes in the mass of the entire ice sheet with a spatial resolution of 300 to 400 km. In contrast, CryoSat-2measures travel times of a radar signal reflected close to the ice sheet surface, allowing changes of the surface topography to be determined with about 5 km spatial resolution. Here, we present a method to combine observations from the both sensors, taking into account the different signal and noise characteristics of each satellite observation that are dependent on the spatial wavelength. We include uncertainties introduced by the processing and corrections, such as the choice of the re-tracking algorithm and the snow/ice volume density model for CryoSat-2, or the filtering of correlated errors and the correction for glacial-isostatic adjustment (GIA) for GRACE. We apply our method to the Antarctic ice sheet and the time period 2011–2017, in which GRACE and CryoSat-2 were simultaneously operational, obtaining a total ice mass loss of 178 ± 23 Gt yr−1. We present a map of the rate of mass change with a spatial resolution of 40 km that is evaluable across all spatial scales, and more precise than estimates based on a single satellite mission.

Description: Surface melt, driven by atmospheric temperatures and albedo, is a strong contribution of mass loss of the Greenland Ice Sheet. In the past, black carbon, algae and other light-absorbing impurities were suggested to govern albedo in Greenland’s ablation zone. Here we combine optical (MODIS/Sentinel-2) and radar (Sentinel-1) remote sensing data with airborne radar and laser scanner data, and engage firn modelling to identify the governing factors leading to dark glacier surfaces in Northeast Greenland. After the drainage of supraglacial lakes, the former lake ground is a clean surface represented by a high reflectance in Sentinel-2 data and aerial photography. These bright spots move with the ice flow and darken by more than 20% over only two years. In contrast, sites further inland do not exhibit this effect. This finding suggests that local deposition of dust, rather than black carbon or cryoconite formation, is the governing factor of albedo of fast-moving outlet glaciers. This is in agreement with a previous field study in the area which finds the mineralogical composition and grain size of the dust comparable with that of the surrounding soils.

Description: In response to increasing Arctic temperatures, ice-rich permafrost landscapes are undergoing rapid changes. In permafrost lowlands, polygonal ice wedges are especially prone to degradation. Melting of ice wedges results in deepening troughs and the transition from low-centered to high-centered ice-wedge polygons. This process has important implications for surface hydrology, as the connectivity of such troughs determines the rate of drainage for these lowland landscapes. In this study, we present a comprehensive, modular, and highly automated workflow to extract, to represent, and to analyze remotely sensed ice-wedge polygonal trough networks as a graph (i.e., network structure). With computer vision methods, we efficiently extract the trough locations as well as their geomorphometric information on trough depth and width from high-resolution digital elevation models and link these data within the graph. Further, we present and discuss the benefits of graph analysis algorithms for characterizing the erosional development of such thaw-affected landscapes. Based on our graph analysis, we show how thaw subsidence has progressed between 2009 and 2019 following burning at the Anaktuvuk River fire scar in northern Alaska, USA. We observed a considerable increase in the number of discernible troughs within the study area, while simultaneously the number of disconnected networks decreased from 54 small networks in 2009 to only six considerably larger disconnected networks in 2019. On average, the width of the troughs has increased by 13.86%, while the average depth has slightly decreased by 10.31%. Overall, our new automated approach allows for monitoring ice-wedge dynamics in unprecedented spatial detail, while simultaneously reducing the data to quantifiable geometric measures and spatial relationships.

Description: To understand the dynamics of ice shelves, a knowledge of their internal and basal structure is very important. As the capacity to perform local surveys is limited, remote sensing provides an opportunity to obtain the relevant information. We must prove, however, that the relevant information can be obtained from remote sensing of the surface. That is the aim of this study. The Jelbart Ice Shelf, Antarctica, exhibits a variety of surface structures appearing as stripe-like features in radar imagery. We performed an airborne geophysical survey across these features and compared the results to TerraSAR-X imagery. We find that the stripe-like structures indicate surface troughs coinciding with the location of basal channels and crevasse-like features, revealed by radio-echo sounding. HH and VV polarizations do not show different magnitude. In surface troughs, the local accumulation rate is larger than at the flat surface. Viscoelastic modelling is used to gain an understanding of the surface undulations and their origin. The surface displacement, computed with a Maxwell model, matches the observed surface reasonably well. Our simulations show that the surface troughs develop over decadal to centennial timescales.

Description: Using in situ data from 2011 and 2013, we evaluate the ability of CryoSat-2 (CS-2) to retrieve sea-ice freeboard over fast ice in McMurdo Sound. This provides the first systematic validation of CS-2 in the coastal Antarctic and offers insight into the assumptions currently used to process CS-2 data. European Space Agency Level 2 (ESAL2) data are compared with results of a Waveform Fitting (WfF) procedure and a Threshold-First-Maximum-Retracker-Algorithm employed at 40% (TFMRA40). A supervised freeboard retrieval procedure is used to reduce errors associated with sea surface height identification and radar velocity in snow. We find ESAL2 freeboards located between the ice and snow freeboard rather than the frequently assumed snow/ice interface. WfF is within 0.04m of the ice freeboard but is influenced by variable snow conditions causing increased radar backscatter from the air/snow interface. Given such snow conditions and additional uncertainties in sea surface height identification, a positive bias of 0.14m away from the ice freeboard is observed. TFMRA40 freeboards are within 0.03m of the snow freeboard. The separation of freeboard estimates is primarily driven by the different assumptions of each retracker, although waveform alteration by variations in snow properties and surface roughness is evident. Techniques are amended where necessary, and automatic freeboard retrieval procedures for ESAL2, WfF and TFMRA40 are presented. CS-2 detects annual fastice freeboard trends using all three automatic procedures that are in line with known sea-ice growth rates in the region.

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.