Search

Due to adverse snow and cloud conditions, only a few inventories are available for the maritime glaciers in New Zealand. These are difficult to compare as different approaches and baseline data have been used to create them. In consequence, glacier fluctuations in New Zealand over the past two decades are only known for a few glaciers based on field observations. Here we present the results of a new inventory for the ‘year 2000’ (some scenes are from 2001 and 2002) that is based on glacier outlines from a recently published inventory for the year 2016 and allowed consistent change assessment for nearly 3000 glaciers over this period. The year 2000 inventory was created by manual on-screen digitizing using Landsat ETM+ satellite imagery (15 m panchromatic band) in the background and the year 2016 outlines as a starting point. Major challenges faced were late and early seasonal snow, clouds and shadow, the geo-location mismatch between Landsat and Sentinel-2 as well as the correct interpretation of ice patches and ice under debris cover. In total, we re-mapped 2967 glaciers covering an area of 885.5 km$^{2}$ in 2000, which is 91.7 km$^{2}$ (or 10.4%) more than the 793.8 km$^{2}$ mapped in 2016. Area change rates (mean rate −0.65% a$^{−1}$) increase towards smaller glaciers. Strongest area loss from 2000 to 2016 occurred at elevations ~1900 m but the highest relative loss was found below 800 m a.s.l. In total, 109 glaciers split into two or more entities and 264 had wasted away by 2016.

Glacier monitoring has been internationally coordinated for more than 125 years. Despite this long history, there is no authoritative answer to the popular question: ‘Which glaciers are the largest in the world?’ Here, we present the first systematic assessment of this question and identify the largest glaciers in the world – distinct from the two ice sheets in Greenland and Antarctica but including the glaciers on the Antarctic Peninsula. We identify the largest glaciers in two domains: on each of the seven geographical continents and in the 19 first-order glacier regions defined by the Global Terrestrial Network for Glaciers. Ranking glaciers by area is non-trivial. It depends on how a glacier is defined and mapped and also requires differentiating between a glacier and a glacier complex, i.e. glaciers that meet at ice divides such as ice caps and icefields. It also depends on the availability of a homogenized global glacier inventory. Using separate rankings for glaciers and glacier complexes, we find that the largest glacier complexes have areas on the order of tens of thousands of square kilometers whereas the largest glaciers are several thousands of square kilometers. The world's largest glaciers and glacier complexes are located in the Antarctic, Arctic and Patagonia.

Climate change, increasing population and changes in land use are all rapidly driving the need to be able to better understand surface water dynamics. The targets set by the United Nations under Sustainable Development Goal 6 in relation to freshwater ecosystems also make accurate surface water monitoring increasingly vital. However, the last decades have seen a steady decline in in situ hydrological monitoring and the availability of the growing volume of environmental data from free and open satellite systems is increasingly being recognized as an essential tool for largescale monitoring of water resources. The scientific literature holds many promising studies on satellite-based surface-water mapping, but a systematic evaluation has been lacking. Therefore, a round robin exercise was organized to conduct an intercomparison of 14 different satellite-based approaches for monitoring inland surface dynamics with Sentinel-1, Sentinel-2, and Landsat 8 imagery. The objective was to achieve a better understanding of the pros and cons of different sensors and models for surface water detection and monitoring. Results indicate that, while using a single sensor approach (applying either optical or radar satellite data) can provide comprehensive results for very specific localities, a dual sensor approach (combining data from both optical and radar satellites) is the most effective way to undertake largescale national and regional surface water mapping across bioclimatic gradients.

Meltwater from the cryosphere contributes a significant fraction of the freshwater resources in the countries receiving water from the Third Pole. Within the ESA-MOST Dragon 4 project, we addressed in particular changes of glaciers and proglacial lakes and their interaction. In addition, we investigated rock glaciers in permafrost environments. Here, we focus on the detailed investigations which have been performed in the Poiqu River Basin, central Himalaya. We used in particular multi-temporal stereo satellite imagery, including high-resolution 1960/70s Corona and Hexagon spy images and contemporary Pleiades data. Sentinel-2 data was applied to assess the glacier flow. The results reveal that glacier mass loss continuously increased with a mass budget of −0.42 ± 0.11 m w.e.a−1 for the period 2004–2018. The mass loss has been primarily driven by an increase in summer temperature and is further accelerated by proglacial lakes, which have become abundant. The glacial lake area more than doubled between 1964 and 2017. The termini of glaciers that flow into lakes moved on average twice as fast as glaciers terminating on land, indicating that dynamical thinning plays an important role. Rock glaciers are abundant, covering approximately 21 km2, which was more than 10% of the glacier area (approximately 190 km2) in 2015. With ongoing glacier wastage, rock glaciers can become an increasingly important water resource.

A number of glaciological observations on debris-covered glaciers around the globe have shown a delayed length and mass adjustment in relation to climate variability, a behavior normally attributed to the ice insulation effect of thick debris layers. Dynamic interactions between debris cover, geometry and surface topography of debris-covered glaciers can nevertheless govern glacier velocities and mass changes over time, with many glaciers exhibiting high thinning rates in spite of thick debris cover. Such interactions are progressively being incorporated into glacier evolution research. In this paper we reconstruct changes in debris-covered area, surface velocities and surface features of three glaciers in the Patagonian Andes over the 1958–2020 period, based on satellite and aerial imagery and Digital Elevation Models. Our results show that debris cover has increased from 40 ± 0.6 to 50 ± 6.7% of the total glacier area since 1958, whilst glacier slope has slightly decreased. The gently sloping tongues have allowed surface flow velocities to remain relatively low (<60 m a−1) for the last two decades, preventing evacuation of surface debris, and contributing to the formation and rise of the ice cliff zone upper boundary. In addition, mapping of end of summer snowline altitudes for the last two decades suggests an increase in the Equilibrium Line Altitudes, which promotes earlier melt out of englacial debris and further increases debris-covered ice area. The strongly negative mass budget of the three investigated glaciers throughout the study period, together with the increases in debris cover extent and ice cliff zones up-glacier, and the low velocities, shows a strong linkage between debris cover, mass balance evolution, surface velocities and topography. Interestingly, the presence of thicker debris layers on the lowermost portions of the glaciers has not lowered thinning rates in these ice areas, indicating that the mass budget is mainly driven by climate variabil

The ongoing glacier shrinkage in the Alps requires frequent updates of glacier outlines to provide an accurate database for monitoring, modelling purposes (e.g. determination of run-off, mass balance, or future glacier extent), and other applications. With the launch of the first Sentinel-2 (S2) satellite in 2015, it became possible to create a consistent, Alpine-wide glacier inventory with an unprecedented spatial resolution of 10 m. The first S2 images from August 2015 already provided excellent mapping conditions for most glacierized regions in the Alps and were used as a base for the compilation of a new Alpine-wide glacier inventory in a collaborative team effort. In all countries, glacier outlines from the latest national inventories have been used as a guide to compile an update consistent with the respective previous interpretation. The automated mapping of clean glacier ice was straightforward using the band ratio method, but the numerous debris-covered glaciers required intense manual editing. Cloud cover over many glaciers in Italy required also including S2 scenes from 2016. The outline uncertainty was determined with digitizing of 14 glaciers several times by all participants. Topographic information for all glaciers was obtained from the ALOS AW3D30 digital elevation model (DEM). Overall, we derived a total glacier area of 1806±60 km2 when considering 4395 glaciers >0.01 km2. This is 14 % (−1.2 % a−1) less than the 2100 km2 derived from Landsat in 2003 and indicates an unabated continuation of glacier shrinkage in the Alps since the mid-1980s. It is a lower-bound estimate, as due to the higher spatial resolution of S2 many small glaciers were additionally mapped or increased in size compared to 2003. Median elevations peak around 3000 m a.s.l., with a high variability that depends on location and aspect. The uncertainty assessment revealed locally strong differences in interpretation of debris-covered glaciers, resulting in limitations for change asses

Glaciers in the Alps and several other regions in the world have experienced strong negative mass balances over the past few decades. Some of them are disappearing, undergoing exceptionally negative mass balances that impact the mean regional value, and require replacement. In this study, we analyse the geomorphometric characteristics of 46 mass-balance glaciers in the Alps and the long-term mass-balance time series for a subset of nine reference glaciers. We identify regime shifts in the mass-balance time series (when non-climatic controls started impacting) and develop a glacier vulnerability index (GVI) as a proxy for their possible future development, based on criteria such as hypsometric index, breaks in slope, thickness distribution and elevation change pattern. We found that the subset of 46 mass-balance glaciers reflects the characteristics of the total glacier sample very well and identified a region-specific variability of the mass balance. As the GVI is strongly related to cumulative glacier mass balances, it can be used as a pre-selector of future mass-balance glaciers. We conclude that measurements on rapidly shrinking glaciers should be continued as long as possible to identify regime shifts in hind-cast and better understand the impacts of climatic variability on such glaciers.

Rock glaciers are an important component of the cryosphere and are one of the most visible manifestations of permafrost. While the significance of rock glacier contribution to streamflow remains uncertain, the contribution is likely to be important for certain parts of the world. High-resolution remote sensing data has permitted the creation of rock glacier inventories for large regions. However, due to the spectral similarity between rock glaciers and the surrounding material, the creation of such inventories is typically conducted based on manual interpretation, which is both time consuming and subjective. Here, we present a novel method that combines deep learning (convolutional neural networks or CNNs) and object-based image analysis (OBIA) into one workflow based on freely available Sentinel-2 optical imagery (10 m spatial resolution), Sentinel-1 interferometric coherence data, and a digital elevation model (DEM). CNNs identify recurring patterns and textures and produce a prediction raster, or heatmap where each pixel indicates the probability that it belongs to a certain class (i.e. rock glacier) or not. By using OBIA we can segment the datasets and classify objects based on their heatmap value as well as morphological and spatial characteristics. We analysed two distinct catchments, the La Laguna catchment in the Chilean semi-arid Andes and the Poiqu catchment in the central Himalaya. In total, our method mapped 108 of the 120 rock glaciers across both catchments with a mean overestimation of 28%. Individual rock glacier polygons howevercontained false positives that are texturally similar, such as debris-flows, avalanche deposits, or fluvial material causing the user's accuracy to be moderate (63.9–68.9%) even if the producer's accuracy was higher (75.0–75.4%). We repeated our method on very-high-resolution Pléiades satellite imagery and a corresponding DEM (at 2 m resolution) for a subset of the Poiqu catchment to ascertain what difference image resolution

Mapping snow cover (SC) on glaciers at the end of the ablation period provides a possibility to rapidly obtain a proxy for their equilibrium line altitude (ELA) which in turn is a metric for the mass balance. Satellite determination of glacier snow cover, derived over large regions, can reveal its spatial variability and temporal trends. Accordingly, snow mapping on glaciers has been widely applied using several satellite sensors. However, as glacier ice originates from compressed snow, both have very similar spectral properties and standard methods to map snow struggle to distinguish snow on glaciers. Hence, most studies applied manual delineation of snow extent on glaciers. Here we present an automated tool, named ‘ASMAG’ (automated snow mapping on glaciers), to map SC on glaciers and derive the related snow line altitude (SLA) for individual glaciers using multi-temporal Landsat satellite imagery and a digital elevation model (DEM). The method has been developed using the example of the Ötztal Alps, where an evaluation of the method is possible using field-based observations of the annual equilibrium line altitude (ELA) and the accumulation area ratio (AAR) measured for three glaciers for more than 30 years. The tool automatically selects a threshold to map snow on glaciers and robustly calculates the SLA based on the frequency distribution of elevation bins with more than 50% SC. The accuracy of the SC mapping was about 90% and the SLA was determined successfully in 80% of all cases with a mean uncertainty of ±19 m. When cloud-free scenes close to the date of the highest snowline are available, a good to very good agreement of SC ratios (SCR)/SLA with field data of AAR/ELA are obtained, otherwise values are systematically higher/lower as useful images were often acquired too early in the summer season. However, glacier specific differences are still well captured. Snow mapping on glaciers is impeded by clouds and their shadows or when fresh snow is covering the

Accurately identifying the extent of surface snow cover on glaciers is important for extrapolating end of year mass balance measurements, constraining the glacier surface radiative energy balance and evaluating model simulations of snow cover. Here, we use auxiliary information from Riegl VZ-6000 Terrestrial Laser Scanner (TLS) return signals to accurately map the snow cover over a glacier throughout an ablation season. Three classification systems were compared, and we find that supervised classification based on TLS signal intensity alone is outperformed by a rule-based classification employing intensity, surface roughness and an associated optical image, which achieves classification accuracy of 68–100%. The TLS intensity signal shows no meaningful relationship with surface or bulk snow density. Finally, we have also compared our Snow Line Altitude (SLA) derived from TLS with SLA derived from the model output, as well as one Landsat image. The results of the model output track the SLA from TLS well, however with a positive bias. In contrast, automatic Landsat-derived SLA slightly underestimates the SLA from TLS. To conclude, we demonstrate that the snow cover extent can be mapped successfully using TLS, although the snow mass remains elusive.

Monitoring glacier changes in remote Arctic regions are strongly facilitated by satellite data. This is especially true for the Russian Arctic where recently increased optical and SAR satellite imagery (Landsat 8 OLI, Sentinel 1/2), and digital elevation models (TanDEM-X, ArcticDEM) are becoming available. These datasets offer new possibilities to create high-quality glacier inventories. Here, we present a new glacier inventory derived from a fusion of multi-source satellite data for Novaya Zemlya in the Russian Arctic. We mainly used Landsat 8 OLI data to automatically map glaciers with the band ratio method. Missing debris-covered glacier parts and misclassified lakes were manually corrected. Whereas perennial snow fields were a major obstacle in glacier identification, seasonal snow was identified and removed using Landsat 5 TM scenes from the year 1998. Drainage basins were derived semi-automatically using the ArcticDEM (gap-filled by the ASTER GDEM V2) and manually corrected using fringes from ALOS PALSAR. The new glacier inventory gives a glacierized area of 22,379 ± 246.16 km2 with 1474 glacier entities >0.05 km2. The region is dominated by large glaciers, as 909 glaciers <0.5 km2 (62% by number) cover only 156 ± 1.7 km2 or 0.7% of the area, whereas 49 glaciers >100 km2 (3.3% by number) cover 18,724 ± 205.9 km2 or 84%. In total, 41 glaciers are marine terminating covering an area of 16,063.7 ± 118.8 km2. The mean elevation is 596 m for all glaciers in the study region (528 m in the northern part, 641 in the southern part). South-east (north-west) facing glaciers cover >35% (20%) of the area. For the smaller glaciers in the southern part we calculated an area loss of ~5% (52.5 ± 4.5 km2) from 2001 to 2016.

Despite renewed efforts to better understand glacier change and recognize glacier change trends in the Andes, relatively large areas in the Andes of Argentina and Chile are still not investigated. In this study, we report on glacier elevation and mass changes in the outer region of the Northern and Southern Patagonian Icefields in the Southern Patagonian Andes. A newly-compiled Landsat ETM+ derived glacier inventory (consisting of 2253 glaciers and ~1314 ± 66 km2 of ice area) and differencing of the SRTM and SPOT5 DEMs were used to derive glacier-specific elevation changes over the 2000–12 period. The investigated glaciers showed a volume change of −0.71 ± 0.55 km3 a−1, yielding a surface lowering of 0.52 ± 0.35 m a−1 on average and an overall mass loss of 0.46 ± 0.37 m w.e. a−1. Highly variable individual glacier responses were observed and interestingly, they were less negative than previously reported for the neighboring Patagonian Icefields.