Correlation dispersion as a measure to better estimate uncertainty in remotely sensed glacier displacements

In recent years a vast amount of glacier surface velocity data from satellite imagery has emerged based on correlation between repeat images. Thereby, much emphasis has been put on the fast processing of large data volumes and products with complete spatial coverage. The metadata of such measurements are often highly simplified when the measurement precision is lumped into a single number for the whole dataset, although the error budget of image matching is in reality neither isotropic nor constant over the whole velocity field. The spread of the correlation peak of individual image offset measurements is dependent on the image structure and the non-uniform flow of the ice and is used here to extract a proxy for measurement uncertainty. A quantification of estimation error or dispersion for each individual velocity measurement can be important for the inversion of, for instance, rheology, ice thickness and/or bedrock friction. Errors in the velocity data can propagate into derived results in a complex and exaggerating way, making the outcomes very sensitive to velocity noise and outliers. Here, we present a computationally fast method to estimate the matching precision of individual displacement measurements from repeat imaging data, focusing on satellite data. The approach is based upon Gaussian fitting directly on the correlation peak and is formulated as a linear least-squares estimation, making its implementation into current pipelines straightforward. The methodology is demonstrated for Sermeq Kujalleq (Jakobshavn Isbræ), Greenland, a glacier with regions of strong shear flow and with clearly oriented crevasses, and Malaspina Glacier, Alaska. Directionality within an image seems to be the dominant factor influencing the correlation dispersion. In our cases these are crevasses and moraine bands, while a relation to differential flow, such as shear, is less pronounced on the correlation spread.</p

Circum-Arctic changes in the flow of glaciers and ice caps from satellite SAR data between the 1990s and 2017

We computed circum-Arctic surface velocitymaps of glaciers and ice caps over the Canadian Arctic, Svalbard and the Russian Arctic for at least two times between the 1990s and 2017 using satellite SAR data. Our analyses are mainly performed with offset-tracking of ALOS-1 PALSAR-1 (2007–2011) and Sentinel-1 (2015–2017) data. In certain cases JERS-1 SAR (1994–1998), TerraSAR-X (2008–2012), Radarsat-2 (2009–2016) and ALOS-2 PALSAR-2 (2015–2016) data were used to fill-in spatial or temporal gaps. Validation of the latest Sentinel-1 results was accomplished by means of SAR data at higher spatial resolution (Radarsat-2Wide Ultra Fine) and ground-basedmeasurements. In general, we observe a deceleration of flow velocities for the major tidewater glaciers in the Canadian Arctic and an increase in frontal velocity along with a retreat of frontal positions over Svalbard and the Russian Arctic. However, all regions have strong accelerations for selected glaciers. The latter developments can be well traced based on the very high temporal sampling of Sentinel-1 acquisitions since 2015, revealing new insights in glacier dynamics. For example, surges on Spitsbergen (e.g., Negribreen, Nathorsbreen, Penckbreen and Strongbreen) have a different characteristic and timing than those over Eastern Austfonna and Edgeoya (e.g., Basin 3, Basin 2 and Stonebreen). Events similar to those ongoing on Eastern Austofonna were also observed over the Vavilov Ice Cap on Severnaya Zemlya and possibly Simony Glacier on Franz-Josef Land. Collectively, there seems to be a recently increasing number of glaciers with frontal destabilization over Eastern Svalbard and the Russian Arctic compared to the 1990s

Sensitivity of glacier volume change estimation to DEM void interpolation

Glacier mass balance has been estimated on individual glacier and regional scales using repeat digital elevation models (DEMs). DEMs often have gaps in coverage (“voids”), the properties of which depend on the nature of the sensor used and the surface being measured. The way that these voids are accounted for has a direct impact on the estimate of geodetic glacier mass balance, though a systematic comparison of different proposed methods has been heretofore lacking. In this study, we determine the impact and sensitivity of void interpolation methods on estimates of volume change. Using two spatially complete, high-resolution DEMs over southeast Alaska, USA, we artificially generate voids in one of the DEMs using correlation values derived from photogrammetric processing of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) scenes. We then compare 11 different void interpolation methods on a glacier-by-glacier and regional basis. We find that a few methods introduce biases of up to 20&thinsp;% in the regional results, while other methods give results very close (&lt;1&thinsp;% difference) to the true, non-voided volume change estimates. By comparing results from a few of the best-performing methods, an estimate of the uncertainty introduced by interpolating voids can be obtained. Finally, by increasing the number of voids, we show that with these best-performing methods, reliable estimates of glacier-wide volume change can be obtained, even with sparse DEM coverage.</p

MMASTER: improved ASTER DEMs for elevation change monitoring

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) system on board the Terra (EOS AM-1) satellite has been a source of stereoscopic images covering the whole globe at 15-m resolution with consistent quality for over 16 years. The potential of these data in terms of geomorphological analysis and change detection in three dimensions is unrivaled and should be exploited more. Due to uncorrected errors in the image geometry due to sensor motion (“jitter”), however, the quality of the DEMs and orthoimages currently available is often insufficient for a number of applications, including surface change detection. We have therefore developed a series of algorithms packaged under the name MicMac ASTER (MMASTER). It is composed of a tool to compute Rational Polynomial Coefficient (RPC) models from the ASTER metadata, a method that improves the quality of the matching by identifying and correcting jitter-induced cross-track parallax errors and a correction for along-track jitter when computing differences between DEMs (either with another MMASTER DEM or with another data source). Our method outputs more precise DEMs with less unmatched areas and reduced overall noise compared to NASA’s standard AST14DMO product. The algorithms were implemented in the open source photogrammetric library and software suite MicMac. Here, we briefly examine the potential of MMASTER-produced DEMs to investigate a variety of geomorphological changes, including river erosion, seismic deformation, changes in biomass, volcanic deformation and glacier mass balance

From high friction zone to frontal collapse: dynamics of an ongoing tidewater glacier surge, Negribreen, Svalbard

Abstract Negribreen, a tidewater glacier located in central eastern Svalbard, began actively surging after it experienced an initial collapse in summer 2016. The surge resulted in horizontal surface velocities of more than 25 m d −1 , making it one of the fastest-flowing glaciers in the archipelago. The last surge of Negribreen likely occurred in the 1930s, but due to a long quiescent phase, investigations of this glacier have been limited. As Negribreen is part of the Negribreen Glacier System, one of the largest glacier systems in Svalbard, investigating its current surge event provides important information on surge behaviour among tidewater glaciers within the region. Here, we demonstrate the surge development and discuss triggering mechanisms using time series of digital elevation models (1969–2018), surface velocities (1995–2018), crevasse patterns and glacier extents from various data sources. We find that the active surge results from a four-stage process. Stage 1 (quiescent phase) involves a long-term, gradual geometry change due to high subglacial friction towards the terminus. These changes allow the onset of Stage 2, an accelerating frontal destabilization, which ultimately results in the collapse (Stage 3) and active surge (Stage 4)

Glacier–permafrost relations in a high-mountain environment: 5 decades of kinematic monitoring at the Gruben site, Swiss Alps

Digitized aerial images were used to monitor the evolution of perennially frozen debris and polythermal glacier ice at the intensely investigated Gruben site in the Swiss Alps over a period of about 50 years. The photogrammetric analysis allowed for a compilation of detailed spatio-temporal information on flow velocities and thickness changes. In addition, high-resolution GNSS (global navigation satellite system) and ground surface temperature measurements were included in the analysis to provide insight into short-term changes. Over time, extremely contrasting developments and landform responses are documented. Viscous flow within the warming and already near-temperate rock glacier permafrost continued at a constant average but seasonally variable speed of typically decimetres per year, with average surface lowering limited to centimetres to a few decimetres per year. This constant flow causes the continued advance of the characteristic convex, lava-stream-like rock glacier with its oversteepened fronts. Thawing rates of ice-rich perennially frozen ground to strong climate forcing are very low (centimetres per year) and the dynamic response strongly delayed (timescale of decades to centuries). The adjacent cold debris-covered glacier tongue remained an essentially concave landform with diffuse margins, predominantly chaotic surface structure, intermediate thickness losses (decimetres per year), and clear signs of down-wasting and decreasing flow velocity. The former contact zone between the cold glacier margin and the upper part of the rock glacier with disappearing remains of buried glacier ice embedded on top of frozen debris exhibits complex phenomena of thermokarst in massive ice and backflow towards the topographic depression produced by the retreating glacier tongue. As is typical for glaciers in the Alps, the largely debris-free glacier part shows a rapid response (timescale of years) to strong climatic forcing with spectacular retreat (>10 m a−1) and mass loss (up to >1 m w.e. specific mass loss per year). The system of periglacial lakes shows a correspondingly dynamic evolution and had to be controlled by engineering work for hazard protection

Recent glacier and lake changes in High Mountain Asia and their relation to precipitation changes

We present an updated, spatially resolved estimate of 2003–2008 glacier surface elevation changes for the entire region of High Mountain Asia (HMA) from ICESat laser altimetry data. The results reveal a diverse pattern that is caused by spatially greatly varying glacier sensitivity, in particular to precipitation availability and changes. We introduce a spatially resolved zonation where ICESat samples are grouped into units of similar glacier behaviour, glacier type and topographic settings. In several regions, our new zonation reveals local differences and anomalies that have not been described previously. Glaciers in the Eastern Pamirs, Kunlun Shan and central TP were thickening by 0.1–0.7 m a−1, and the thickening anomaly has a crisp boundary in the Eastern Pamirs that continues just north of the central Karakoram. Glaciers in the south and east of the TP were thinning, with increasing rates towards southeast. We attribute the glacier thickening signal to a stepwise increase in precipitation around ∼1997–2000 on the Tibetan Plateau (TP). The precipitation change is reflected by growth of endorheic lakes in particular in the northern and eastern TP. We estimate lake volume changes through a combination of repeat lake extents from Landsat data and shoreline elevations from ICESat and the Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) for over 1300 lakes. The rise in water volume contained in the lakes corresponds to 4–25 mm a−1, when distributed over entire catchments, for the areas where we see glacier thickening. The precipitation increase is also visible in sparse in situ measurements and MERRA-2 climate reanalysis data but less visible in ERA-Interim reanalysis data. Taking into account evaporation loss, the difference between average annual precipitation during the 1990s and 2000s suggested by these datasets is 34–100 mm a−1, depending on region, which can fully explain both lake growth and glacier thickening (Kunlun Shan) or glacier geometry changes such as thinning tongues while upper glacier areas were thickening or stable (eastern TP). The precipitation increase reflected in these glacier changes possibly extended to the northern slopes of the Tarim Basin, where glaciers were nearly in balance in 2003–2008. Along the entire Himalaya, glaciers on the first orographic ridge, which are exposed to abundant precipitation, were thinning less than glaciers in the dryer climate of the inner ranges. Thinning rates in the Tien Shan vary spatially but are rather stronger than in other parts of HMA

Mass-balance reconstruction for Glacier No. 354, Tien Shan, from 2003 to 2014

This study presents a reconstruction of the seasonal mass balance of Glacier No. 354, located in the Akshiirak range, Kyrgyzstan, from 2003 to 2014. We use a distributed accumulation and temperature-index melt model driven by daily air temperature and precipitation from a nearby meteorological station. The model is calibrated with in situ measurements of the annual mass balance collected from 2011 to 2014. The snow-cover depletion pattern observed using satellite imagery provides additional information on the dynamics of mass change throughout the melting season. Two digital elevation models derived from high-resolution satellite stereo images acquired in 2003 and 2012 are used to calculate glacier volume change for the corresponding period. The geodetic mass change thus derived is used to validate the modelled cumulative glacier-wide balance. For the period 2003–12 we find a cumulative mass balance of –0.40±10mw.e.a-1. This result agrees well with the geodetic balance of –0.48±0.07mw.e.a-1over the same period

Three different glacier surges at a spot: what satellites observe and what not

In the Karakoram, dozens of glacier surges occurred in the past 2 decades, making the region a global hotspot. Detailed analyses of dense time series from optical and radar satellite images revealed a wide range of surge behaviour in this region: from slow advances longer than a decade at low flow velocities to short, pulse-like advances over 1 or 2 years with high velocities. In this study, we present an analysis of three currently surging glaciers in the central Karakoram: North and South Chongtar Glaciers and an unnamed glacier referred to as NN9. All three glaciers flow towards the same small region but differ strongly in surge behaviour. A full suite of satellites (e.g. Landsat, Sentinel-1 and 2, Planet, TerraSAR-X, ICESat-2) and digital elevation models (DEMs) from different sources (e.g. Shuttle Radar Topography Mission, SRTM; Satellite Pour l’Observation de la Terre, SPOT; High Mountain Asia DEM, HMA DEM) are used to (a) obtain comprehensive information about the evolution of the surges from 2000 to 2021 and (b) to compare and evaluate capabilities and limitations of the different satellite sensors for monitoring surges of relatively small glaciers in steep terrain. A strongly contrasting evolution of advance rates and flow velocities is found, though the elevation change pattern is more similar. For example, South Chongtar Glacier had short-lived advance rates above 10 km yr−1, velocities up to 30 m d−1, and surface elevations increasing by 170 m. In contrast, the neighbouring and 3-times-smaller North Chongtar Glacier had a slow and near-linear increase in advance rates (up to 500 m yr−1), flow velocities below 1 m d−1 and elevation increases up to 100 m. The even smaller glacier NN9 changed from a slow advance to a full surge within a year, reaching advance rates higher than 1 km yr−1. It seems that, despite a similar climatic setting, different surge mechanisms are at play, and a transition from one mechanism to another can occur during a single surge. The sensor inter-comparison revealed a high agreement across sensors for deriving flow velocities, but limitations are found on small and narrow glaciers in steep terrain, in particular for Sentinel-1. All investigated DEMs have the required accuracy to clearly show the volume changes during the surges, and elevations from ICESat-2 ATL03 data fit neatly to the other DEMs. We conclude that the available satellite data allow for a comprehensive observation of glacier surges from space when combining different sensors to determine the temporal evolution of length, elevation and velocity changes