Description: Ice shelves control the dynamic mass loss of ice sheets through buttressing and their integrity depends on the spatial variability of their basal mass balance (BMB), i.e. the difference between refreezing and melting. Here, we present an improved technique – based on satellite observations – to capture the small-scale variability in the BMB of ice shelves. As a case study, we apply the methodology to the Roi Baudouin Ice Shelf, Dronning Maud Land, East Antarctica, and derive its yearly averaged BMB at 10 m horizontal gridding. We use mass conservation in a Lagrangian framework based on high-resolution surface velocities, atmospheric-model surface mass balance and hydrostatic ice-thickness fields (derived from TanDEM-X surface elevation). Spatial derivatives are implemented using the total-variation differentiation, which preserves abrupt changes in flow velocities and their spatial gradients. Such changes may reflect a dynamic response to localized basal melting and should be included in the mass budget. Our BMB field exhibits much spatial detail and ranges from −14.7 to 8.6 m a−1 ice equivalent. Highest melt rates are found close to the grounding line where the pressure melting point is high, and the ice shelf slope is steep. The BMB field agrees well with on-site measurements from phase-sensitive radar, although independent radar profiling indicates unresolved spatial variations in firn density. We show that an elliptical surface depression (10 m deep and with an extent of 0.7 km × 1.3 km) lowers by 0.5 to 1.4 m a−1, which we tentatively attribute to a transient adaptation to hydrostatic equilibrium. We find evidence for elevated melting beneath ice shelf channels (with melting being concentrated on the channel's flanks). However, farther downstream from the grounding line, the majority of ice shelf channels advect passively (i.e. no melting nor refreezing) toward the ice shelf front. Although the absolute, satellite-based BMB values remain uncertain, we have high confidence in the spatial variability on sub-kilometre scales. This study highlights expected challenges for a full coupling between ice and ocean models.

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: Coastal infrastructure, cultural, and archeological sites are increasingly vulnerable to erosion and flooding along permafrost coasts. Amplified warming of the Arctic, sea level rise, lengthening of the open water period, and a predicted increase in frequency of major storms compound these threats. Mitigation necessitates decision-making tools at an appropriate scale. We present a study of coastal erosion combining it with a flooding risk assessment for the culturally important historic settlement on Herschel Island, a UNESCO World Heritage candidate site. The resulting map may help local stakeholders devise management strategies to cope with rapidly changing environmental conditions. We analyzed shoreline movement using the Digital Shoreline Analysis System (DSAS) after digitizing shorelines from 1952, 1970, and 2011. Using these data, forecasts of shoreline positions were made for 20 and 50 years into the future. Flooding risk was assessed using a cost-distance map based on a high-resolution Light Detection and Ranging (LiDAR) dataset and current Intergovernmental Panel on Climate Change sea level estimates. Widespread erosion characterizes the study area. The rate of shoreline movement for different periods of the study ranges from -5.5 to 2.7 m·a-1 (mean -0.6 m·a-1). Mean coastal retreat decreased from -0.6 m·a-1 to -0.5 m·a-1, for 1952-1970 and 1970-2000, respectively, and increased to -1.3 m·a-1 in the period 2000-2011. Ice-rich coastal sections, and coastal sections most exposed to wave attack exhibited the highest rates of coastal retreat. The geohazard map resulting from shoreline projections and flood risk analysis indicates that most of the area occupied by the historic settlement is at extreme or very high risk of flooding, and some buildings are vulnerable to coastal erosion. The results of this study indicate a greater threat by coastal flooding than erosion. Our assessment may be applied in other locations where limited data are available.

Description: Arctic coastal infrastructure, cultural, and archeological sites are increasingly vulnerable to erosion and flooding due to amplified warming of the Arctic, sea level rise, lengthening of open water periods, and a predicted increase in frequency of major storms. Mitigating these hazards necessitates decision-making tools at an appropriate scale. The objectives of this paper are to provide such a tool by assessing potential erosion and flood hazards at Herschel Island, a UNESCO World Heritage candidate site. This study focused on Simpson Point and the adjacent coastal sections, because of their archeological, historical, and cultural significance. Shoreline movement was analyzed using the Digital Shoreline Analysis System (DSAS) after digitizing shorelines from 1952, 1970, 2000, and 2011. For purposes of this analysis, the coast was divided in seven coastal reaches (CRs) reflecting different morphologies and/or exposures. Using linear regression rates obtained from these data, projections of shoreline position were made for 20 and 50 years into the future. Flood hazard was assessed using a least cost-path analysis based on a high-resolution Light Detection and Ranging (LiDAR) dataset and current Intergovernmental Panel on Climate Change sea level estimates. Widespread erosion characterizes the study area. The rate of shoreline movement in different periods of the study ranges from -5.5 to 2.7 m·a-1 (mean -0.6 m·a-1). Mean coastal retreat decreased from -0.6 m·a-1 to -0.5 m·a-1, for 1952-1970 and 1970-2000, respectively, and increased to -1.3 m·a-1 in the period 2000-2011. Ice-rich coastal sections most exposed to wave attack exhibited the highest rates of coastal retreat. The geohazard map combines shoreline projections and flood hazard analyses to show that most of the spit area has extreme or very high flood hazard potential, and some buildings are vulnerable to coastal erosion. This study demonstrates that transgressive forcing may provide ample sediment for the expansion of depositional landforms, while growing more susceptible to overwash and flooding.

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.