Repeat optical satellite images reveal widespread and long term decrease in land-terminating glacier speeds

By matching of repeat optical satellite images it is now possible to investigate glacier dynamics within large regions of the world and also between regions to improve knowledge about glacier dynamics in space and time. In this study we investigate whether the negative glacier mass balance seen over large parts of the world has caused the glaciers to change their speeds. The studied regions are Pamir, Caucasus, Penny Ice Cap, Alaska Range and Patagonia. In addition we derive speed changes for Karakoram, a region assumed to have positive mass balance and that contains many surge-type glaciers. We find that the mapped glaciers in the five regions with negative mass balance have over the last decades decreased their velocity at an average rate per decade of: 43 % in the Pamir, 8 % in the Caucasus, 25 % on Penny Ice Cap, 11 % in the Alaska Range and 20 % in Patagonia. Glaciers in Karakoram have generally increased their speeds, but surging glaciers and glaciers with flow instabilities are most prominent in this area. Therefore the calculated average speed change is not representative for this area

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

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

Geochemical characterization of supraglacial debris via in situ and optical remote sensing methods: a case study in Khumbu Himalaya, Nepal

Surface glacier debris samples and field spectra were collected from the ablation zones of Nepal Himalaya Ngozumpa and Khumbu glaciers in November and December 2009. Geochemical and mineral compositions of supraglacial debris were determined by X-ray diffraction and X-ray fluorescence spectroscopy. This composition data was used as ground truth in evaluating field spectra and satellite supraglacial debris composition and mapping methods. Satellite remote sensing methods for characterizing glacial surface debris include visible to thermal infrared hyper- and multispectral reflectance and emission signature identification, semi-quantitative mineral abundance indicies and spectral image composites. Satellite derived supraglacial debris mineral maps displayed the predominance of layered silicates, hydroxyl-bearing and calcite minerals on Khumbu Himalayan glaciers. Supraglacial mineral maps compared with satellite thermal data revealed correlations between glacier surface composition and glacier surface temperature. Glacier velocity displacement fields and shortwave, thermal infrared false color composites indicated the magnitude of mass flux at glacier confluences. The supraglacial debris mapping methods presented in this study can be used on a broader scale to improve, supplement and potentially reduce errors associated with glacier debris radiative property, composition, areal extent and mass flux quantifications

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

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

Cold-regions river flow observed from space

Knowledge of water-surface velocities in rivers is useful for understanding a wide range of lotic processes and systems, such as water and ice fluxes and forces, mixing, solute and sediment transport, bed and bank stability, aquatic and riparian ecology, and extreme hydrologic events. In cold regions, river-ice break up and the associated downstream transport of ice debris is often the most important hydrological event of the year, producing flood levels that commonly exceed those for the open-water period and dramatic consequences for river infrastructure and ecology. Quantification of river surface velocity and currents has relied mostly on very scarce in situ measurements or particle tracking in laboratory models, with few attempts to cover entire river reaches. Accurate and complete surface-velocity fields on rivers have rarely been produced. Here, we use river-ice debris as an index of surface water velocity, and track it over a time period of about one minute, which is the typical time lapse between the two or more images that form a stereo data set in spaceborne, alongtrack optical-stereo mapping. In this way, we measure and visualize for the first time, the almost complete surface velocity field of a river. Examples are used from approximately 80 km and 40 km long reaches of the St. Lawrence and Mackenzie rivers, respectively. The methodology and results will be valuable to a number of disciplines requiring detailed information about river flow, such as hydraulics, hydrology, river ecology and natural-hazard management

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)

Stagnation and mass loss on a Himalayan debris-covered glacier: processes, patterns and rates

This research was supported financially by the University Centre in Svalbard (UNIS), National Geographic Society GRANT #W135-10, The Natural Environmental Research Council and the European Commission FP7-MC-IEF.The ablation areas of debris-covered glaciers typically consist of a complex mosaic of surface features with contrasting processes and rates of mass loss. This greatly complicates glacier response to climate change, and increases the uncertainty of predictive models. In this paper we present a series of high-resolution DEMs and repeat lake bathymetric surveys on Ngozumpa Glacier, Nepal, to study processes and patterns of mass loss on a Himalayan debris-covered glacier in unprecedented detail. Most mass loss occurs by melt below supraglacial debris, and melt and calving of ice cliffs (backwasting). Although ice cliffs cover only ∼5% of the area of the lower tongue, they account for 40% of the ablation. The surface debris layer is subject to frequent re-distribution by slope processes, resulting in large spatial and temporal differences in debris-layer thickness, enhancing or inhibiting local ablation rates and encouraging continuous topographic inversion. A moraine-dammed lake on the lower glacier tongue (Spillway Lake) underwent a period of rapid expansion from 2001 to 2009, but later experienced a reduction of area and volume as a result of lake level lowering and sediment redistribution. Rapid lake growth will likely resume in the near future, and may eventually become up to 7 km long.Publisher PDFPeer reviewe

Sedimentological characterization of Antarctic moraines using UAVs and Structure-from-Motion photogrammetry

In glacial environments particle-size analysis of moraines provides insights into clast origin, transport history, depositional mechanism and processes of reworking. Traditional methods for grain-size classification are labour-intensive, physically intrusive and are limited to patch-scale (1m2) observation. We develop emerging, high-resolution ground- and unmanned aerial vehicle-based ‘Structure-from-Motion’ (UAV-SfM) photogrammetry to recover grain-size information across an moraine surface in the Heritage Range, Antarctica. SfM data products were benchmarked against equivalent datasets acquired using terrestrial laser scanning, and were found to be accurate to within 1.7 and 50mm for patch- and site-scale modelling, respectively. Grain-size distributions were obtained through digital grain classification, or ‘photo-sieving’, of patch-scale SfM orthoimagery. Photo-sieved distributions were accurate to <2mm compared to control distributions derived from dry sieving. A relationship between patch-scale median grain size and the standard deviation of local surface elevations was applied to a site-scale UAV-SfM model to facilitate upscaling and the production of a spatially continuous map of the median grain size across a 0.3 km2 area of moraine. This highly automated workflow for site scale sedimentological characterization eliminates much of the subjectivity associated with traditional methods and forms a sound basis for subsequent glaciological process interpretation and analysis