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

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)

Modeling past and future surface mass balance of the Northern Patagonia Icefield

Glaciers are strongly retreating and thinning in Patagonia. We present new inferences about the climatic situation and the surface mass balance on the Northern Patagonia Icefield in the past and the future using a combined modeling approach. The simulations are driven by NCAR/NCEP Reanalysis and ECHAM5 data, which were physically downscaled using the Weather Research and Forecasting regional climate model and simple sub-grid parameterizations. The surface mass balance model was calibrated with geodetic mass balance data of three large non-calving glaciers and with point mass balance measurements. An increase of accumulation on the Northern Patagonia Icefield was detected from 1990–2011 as compared to 1975–1990. Using geodetic mass balance data, calving losses from the Northern Patagonia Icefield could be inferred, which doubled in 2000–2009 as compared to 1975–2000. The 21st century projection of future mass balance of the Northern Patagonia Icefield shows a strong increase in ablation from 2050 and a reduction of solid precipitation from 2080, both due to higher temperatures. The total mass loss in the 21st century is estimated to be 592±50 Gt with strongly increasing rates towards the end of the century. The prediction of the future mass balance of the Northern Patagonia Icefield includes several additional sources of errors due to uncertainties in the prediction of future climate and due to possible variations in ice dynamics, which might modify the geometry of the icefield and change the rate of mass losses due to calving

Estimating the volume of glaciers in the Himalayan–Karakoram region using different methods

Ice volume estimates are crucial for assessing water reserves stored in glaciers. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan–Karakoram (HK) region. In this study, different existing methodologies are used to estimate the ice reserves: three area–volume relations, one slope-dependent volume estimation method, and two ice-thickness distribution models are applied to a recent, detailed, and complete glacier inventory of the HK region, spanning over the period 2000–2010 and revealing an ice coverage of 40 775 km2. An uncertainty and sensitivity assessment is performed to investigate the influence of the observed glacier area and important model parameters on the resulting total ice volume. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 to 4737 km3, depending on the approach. This range is lower than most previous estimates. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses. However, total volume estimates from area-related relations are larger than those from other approaches. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation

Re-establishing glacier monitoring in Kyrgyzstan and Uzbekistan, Central Asia

Glacier mass loss is among the clearest indicators of atmospheric warming. The observation of these changes is one of the major objectives of the international climate monitoring strategy developed by the Global Climate Observing System (GCOS). Long-term glacier mass balance measurements are furthermore the basis for calibrating and validating models simulating future runoff of glacierised catchments. This is essential for Central Asia, which is one of the driest continental regions of the Northern Hemisphere. In the highly populated regions, water shortage due to decreased glacierisation potentially leads to pronounced political instability, drastic ecological changes and endangered food security. As a consequence of the collapse of the former Soviet Union, however, many valuable glacier monitoring sites in the Tien Shan and Pamir Mountains were abandoned. In recent years, multinational actors have re-established a set of important in situ measuring sites to continue the invaluable long-term data series. This paper introduces the applied monitoring strategy for selected glaciers in the Kyrgyz and Uzbek Tien Shan and Pamir, highlights the existing and the new measurements on these glaciers, and presents an example for how the old and new data can be combined to establish multi-decadal mass balance time series. This is crucial for understanding the impact of climate change on glaciers in this region

Reconstructions of the 1900–2015 Greenland ice sheet surface mass balance using the regional climate MAR model

With the aim of studying the recent Greenland ice sheet (GrIS) surface mass balance (SMB) decrease relative to the last century, we have forced the regional climate MAR (Modèle Atmosphérique Régional; version 3.5.2) model with the ERA-Interim (ECMWF Interim Re-Analysis; 1979–2015), ERA-40 (1958–2001), NCEP–NCARv1 (National Centers for Environmental Prediction–National Center for Atmospheric Research Reanalysis version 1; 1948–2015), NCEP–NCARv2 (1979–2015), JRA-55 (Japanese 55-year Reanalysis; 1958–2014), 20CRv2(c) (Twentieth Century Reanalysis version 2; 1900–2014) and ERA-20C (1900–2010) reanalyses. While all these forcing products are reanalyses that are assumed to represent the same climate, they produce significant differences in the MAR-simulated SMB over their common period. A temperature adjustment of +1 °C (respectively −1 °C) was, for example, needed at the MAR boundaries with ERA-20C (20CRv2) reanalysis, given that ERA-20C (20CRv2) is ∼ 1 °C colder (warmer) than ERA-Interim over Greenland during the period 1980–2010. Comparisons with daily PROMICE (Programme for Monitoring of the Greenland Ice Sheet) near-surface observations support these adjustments. Comparisons with SMB measurements, ice cores and satellite-derived melt extent reveal the most accurate forcing datasets for the simulation of the GrIS SMB to be ERA-Interim and NCEP–NCARv1. However, some biases remain in MAR, suggesting that some improvements are still needed in its cloudiness and radiative schemes as well as in the representation of the bare ice albedo. Results from all MAR simulations indicate that (i) the period 1961–1990, commonly chosen as a stable reference period for Greenland SMB and ice dynamics, is actually a period of anomalously positive SMB (∼ +40 Gt yr−1) compared to 1900–2010; (ii) SMB has decreased significantly after this reference period due to increasing and unprecedented melt reaching the highest rates in the 120-year common period; (iii) before 1960, both ERA-20C and 20CRv2-forced MAR simulations suggest a significant precipitation increase over 1900–1950, but this increase could be the result of an artefact in the reanalyses that are not well-enough constrained by observations during this period and (iv) since the 1980s, snowfall is quite stable after having reached a maximum in the 1970s. These MAR-based SMB and accumulation reconstructions are, however, quite similar to those from Box (2013) after 1930 and confirm that SMB was quite stable from the 1940s to the 1990s. Finally, only the ERA-20C-forced simulation suggests that SMB during the 1920–1930 warm period over Greenland was comparable to the SMB of the 2000s, due to both higher melt and lower precipitation than normal

Changes in Greenland’s peripheral glaciers linked to the North Atlantic Oscillation

Glaciers and ice caps peripheral to the main Greenland Ice Sheet contribute markedly to sea-level rise1,2,3. Their changes and variability, however, have been difficult to quantify on multi-decadal timescales due to an absence of long-term data4. Here, using historical aerial surveys, expedition photographs, spy satellite imagery and new remote-sensing products, we map glacier length fluctuations of approximately 350 peripheral glaciers and ice caps in East and West Greenland since 1890. Peripheral glaciers are found to have recently undergone a widespread and significant retreat at rates of 12.2 m per year and 16.6 m per year in East and West Greenland, respectively; these changes are exceeded in severity only by the early twentieth century post-Little-Ice-Age retreat. Regional changes in ice volume, as reflected by glacier length, are further shown to be related to changes in precipitation associated with the North Atlantic Oscillation (NAO), with a distinct east–west asymmetry; positive phases of the NAO increase accumulation, and thereby glacier growth, in the eastern periphery, whereas opposite effects are observed in the western periphery. Thus, with projected trends towards positive NAO in the future5,6, eastern peripheral glaciers may remain relatively stable, while western peripheral glaciers will continue to diminish

The future sea-level rise contribution of Greenland’s glaciers and ice caps

We calculate the future sea-level rise contribution from the surface mass balance of all of Greenland's glaciers and ice caps (GICs, ~90 000 km2) using a simplified energy balance model which is driven by three future climate scenarios from the regional climate models HIRHAM5, RACMO2 and MAR. Glacier extent and surface elevation are modified during the mass balance model runs according to a glacier retreat parameterization. Mass balance and glacier surface change are both calculated on a 250 m resolution digital elevation model yielding a high level of detail and ensuring that important feedback mechanisms are considered. The mass loss of all GICs by 2098 is calculated to be 2016 ± 129 Gt (HIRHAM5 forcing), 2584 ± 109 Gt (RACMO2) and 3907 ± 108 Gt (MAR). This corresponds to a total contribution to sea-level rise of 5.8 ± 0.4, 7.4 ± 0.3 and 11.2 ± 0.3 mm, respectively. Sensitivity experiments suggest that mass loss could be higher by 20–30% if a strong lowering of the surface albedo were to take place in the future. It is shown that the sea-level rise contribution from the north-easterly regions of Greenland is reduced by increasing precipitation while mass loss in the southern half of Greenland is dominated by steadily decreasing summer mass balances. In addition we observe glaciers in the north-eastern part of Greenland changing their characteristics towards greater activity and mass turnover

Firn data compilation reveals widespread decrease of firn air content in western Greenland

The perennial snow, or firn, on the Greenland ice sheet each summer stores part of the meltwater formed at the surface, buffering the ice sheet’s contribution to sea level. We gathered observations of firn air content, indicative of the space available in the firn to retain meltwater, and find that this air content remained stable in cold regions of the firn over the last 65 years but recently decreased significantly in western Greenland