The length of the world’s glaciers – a new approach for the global calculation of center lines

Glacier length is an important measure of glacier geometry. Nevertheless, global glacier inventories are mostly lacking length data. Only recently semi-automated approaches to measure glacier length have been developed and applied regionally. Here we present a first global assessment of glacier length using an automated method that relies on glacier surface slope, distance to the glacier margins and a set of trade-off functions. The method is developed for East Greenland, evaluated for East Greenland as well as for Alaska and eventually applied to all ~ 200 000 glaciers around the globe. The evaluation highlights accurately calculated glacier length where digital elevation model (DEM) quality is high (East Greenland) and limited accuracy on low-quality DEMs (parts of Alaska). Measured length of very small glaciers is subject to a certain level of ambiguity. The global calculation shows that only about 1.5% of all glaciers are longer than 10 km, with Bering Glacier (Alaska/Canada) being the longest glacier in the world at a length of 196 km. Based on the output of our algorithm we derive global and regional area–length scaling laws. Differences among regional scaling parameters appear to be related to characteristics of topography and glacier mass balance. The present study adds glacier length as a key parameter to global glacier inventories. Global and regional scaling laws might prove beneficial in conceptual glacier models

Assessment of evolution and risks of glacier lake outbursts in the Djungarskiy Alatau, Central Asia, using Landsat imagery and glacier bed topography modelling

Changes in the abundance and area of mountain lakes in the Djungarskiy (Jetysu) Alatau between 2002 and 2014 were investigated using Landsat imagery. The number of lakes increased by 6.2 % from 599 to 636 with a growth rate of 0.51 % a−1. The combined areas were 16.26 ± 0.85 to 17.35 ± 0.92 km2 respectively and the overall change was within the uncertainty of measurements. Fifty lakes, whose potential outburst can damage existing infrastructure, were identified. The glacier bed topography version 2 (GlabTop2) model was applied to simulate ice thickness and subglacial topography using glacier outlines for 2000 and SRTM DEM (Shuttle Radar Topography Mission digital elevation model) as input data achieving realistic patterns of ice thickness. A total of 513 overdeepenings in the modelled glacier beds, presenting potential sites for the development of lakes, were identified with a combined area of 14.7 km2. Morphometric parameters of the modelled overdeepenings were close to those of the existing lakes. A comparison of locations of the overdeepenings and newly formed lakes in the areas de-glacierized in 2000–2014 showed that 67 % of the lakes developed at the sites of the overdeepenings. The rates of increase in areas of new lakes correlated with areas of modelled overdeepenings. Locations where hazardous lakes may develop in the future were identified. The GlabTop2 approach is shown to be a useful tool in hazard management providing data on the potential evolution of future lakes

A daily, 1 km resolution data set of downscaled Greenland ice sheet surface mass balance (1958–2015)

This study presents a data set of daily, 1 km resolution Greenland ice sheet (GrIS) surface mass balance (SMB) covering the period 1958–2015. Applying corrections for elevation, bare ice albedo and accumulation bias, the high-resolution product is statistically downscaled from the native daily output of the polar regional climate model RACMO2.3 at 11 km. The data set includes all individual SMB components projected to a down-sampled version of the Greenland Ice Mapping Project (GIMP) digital elevation model and ice mask. The 1 km mask better resolves narrow ablation zones, valley glaciers, fjords and disconnected ice caps. Relative to the 11 km product, the more detailed representation of isolated glaciated areas leads to increased precipitation over the southeastern GrIS. In addition, the downscaled product shows a significant increase in runoff owing to better resolved low-lying marginal glaciated regions. The combined corrections for elevation and bare ice albedo markedly improve model agreement with a newly compiled data set of ablation measurements

The length of the glaciers in the world:a straightforward method for the automated calculation of glacier center lines

Glacier length is an important measure of glacier geometry but global glacier inventories are mostly lacking length data. Only recently semi-automated approaches to measure glacier length have been developed and applied regionally. Here we present a first global assessment of glacier length using a fully automated method based on glacier surface slope, distance to the glacier margins and a set of trade-off functions. The method is developed for East Greenland, evaluated for the same area as well as for Alaska, and eventually applied to all ∼200 000 glaciers around the globe. The evaluation highlights accurately calculated glacier length where DEM quality is good (East Greenland) and limited precision on low quality DEMs (parts of Alaska). Measured length of very small glaciers is subject to a certain level of ambiguity. The global calculation shows that only about 1.5% of all glaciers are longer than 10 km with Bering Glacier (Alaska/Canada) being the longest glacier in the world at a length of 196 km. Based on model output we derive global and regional area-length scaling laws. Differences among regional scaling parameters appear to be related to characteristics of topography and glacier mass balance. The present study adds glacier length as a central parameter to global glacier inventories. Global and regional scaling laws might proof beneficial in conceptual glacier models

Quantifying mass balance processes on the Southern Patagonia Icefield

Artículo de publicación ISIWe present surface mass balance simulations of the Southern Patagonia Icefield (SPI) driven by downscaled reanalysis data. The simulations were evaluated and interpreted using geodetic mass balances, measured point balances and a complete velocity field of the icefield for spring 2004. The high measured accumulation of snow of up to 15.4 m w.e. yr−1 (meters water equivalent per year) as well as the high measured ablation of up to 11 m w.e. yr−1 is reproduced by the model. The overall modeled surface mass balance was positive and increasing during 1975–2011. Subtracting the surface mass balance from geodetic balances, calving fluxes were inferred. Mass losses of the SPI due to calving were strongly increasing from 1975–2000 to 2000– 2011 and higher than losses due to surface melt. Calving fluxes were inferred for the individual glacier catchments and compared to fluxes estimated from velocity data. Measurements of ice thickness and flow velocities at the glaciers’ front and spatially distributed accumulation measurements can help to reduce the uncertainties of the different terms in the mass balance of the Southern Patagonia Icefield.FONDECYT 3140135 European Union 22637

Long-term firn and mass balance modelling for Abramov Glacier in the data-scarce Pamir Alay

Several studies identified heterogeneous glacier mass changes in western High Mountain Asia over the last decades. Causes for these mass change patterns are still not fully understood. Modelling the physical interactions between glacier surface and atmosphere over several decades can provide insight into relevant processes. Such model applications, however, have data needs which are usually not met in these data-scarce regions. Exceptionally detailed glaciological and meteorological data exist for the Abramov Glacier in the Pamir Alay range. In this study, we use weather station measurements in combination with downscaled reanalysis data to force a coupled surface energy balance–multilayer subsurface model for Abramov Glacier for 52 years. Available in situ data are used for model calibration and validation. We find an overall negative mass balance of −0.27 mw.e.a-1 for 1968/1969–2019/2020 and a loss of firn pore space causing a reduction of internal accumulation. Despite increasing air temperatures, we do not find an acceleration of glacier-wide mass loss over time. Such an acceleration is compensated for by increasing precipitation rates (+0.0022 mw.e.a-1, significant at a 90 % confidence level). Our results indicate a significant correlation between annual mass balance and precipitation (R2 = 0.72).</p

Mass-balance parameters derived from a synthetic network of mass-balance glaciers

Glacier mass-balance parameters such as the equilibrium-line altitude (ELA) play an important role when working with large glacier samples. While the number of observational mass-balance series to derive such parameters is limited, more and more modeled data are becoming available. Here we explore the possibilities of analyzing such ‘synthetic’ mass-balance data with respect to mass-balance parameters. A simplified energy-balance model is driven by bias-corrected regional climate model output to model mass-balance distributions for 94 glaciers in the Swiss Alps over 15 years. The modeling results in realistic interannual variability and mean cumulative mass balance. Subsequently model output is analyzed with respect to 18 topographic and mass-balance parameters and a correlation analysis is performed. Well-known correlations such as for ELA and median elevation are confirmed from the synthetic data. Furthermore, previously unreported parameter relationships are found such as a correlation of the balance rate at the tongue with the accumulation–area ratio (AAR) and of the glacier elevation range with the AAR. Analyzing modeled data complements in situ observations and highlights their importance: the small number of accurate mass-balance observations available for validation is a major challenge for the presented approach

Modelling glacier-bed overdeepenings and possible future lakes for the glaciers in the Himalaya—Karakoram region

Surface digital elevation models (DEMs) and slope-related estimates of glacier thickness enable modelling of glacier-bed topographies over large ice-covered areas. Due to the erosive power of glaciers, such bed topographies can contain numerous overdeepenings, which when exposed following glacier retreat may fill with water and form new lakes. In this study, the bed overdeepenings for ~28 000 glaciers (40 775 km²) of the Himalaya-Karakoram region are modelled using GlabTop2 (Glacier Bed Topography model version 2), in which ice thickness is inferred from surface slope by parameterizing basal shear stress as a function of elevation range for each glacier. The modelled ice thicknesses are uncertain (±30%), but spatial patterns of ice thickness and bed elevation primarily depend on surface slopes as derived from the DEM and, hence, are more robust. About 16 000 overdeepenings larger than 10⁴m² were detected in the modelled glacier beds, covering an area of ~2200 km² and having a volume of ~120km³ (3-4% of present-day glacier volume). About 5000 of these overdeepenings (1800 km²) have a volume larger than 10⁶m³. The results presented here are useful for anticipating landscape evolution and potential future lake formation with associated opportunities (tourism, hydropower) and risks (lake outbursts)