• Other research productsclose
• Open Access close
• NO close
• The Cryosphere (TC) close

With dense SAR satellite data time series it is possible to map surface and subsurface glacier properties that vary in time. On Sentinel-1A and RADARSAT-2 backscatter time series images over mainland Norway and Svalbard, we outline how to map glaciers using descriptive methods. We present five application scenarios. The first shows potential for tracking transient snow lines with SAR backscatter time series and correlates with both optical satellite images (Sentinel-2A and Landsat 8) and equilibrium line altitudes derived from in situ surface mass balance data. In the second application scenario, time series representation of glacier facies corresponding to SAR glacier zones shows potential for a more accurate delineation of the zones and how they change in time. The third application scenario investigates the firn evolution using dense SAR backscatter time series together with a coupled energy balance and multilayer firn model. We find strong correlation between backscatter signals with both the modeled firn air content and modeled wetness in the firn. In the fourth application scenario, we highlight how winter rain events can be detected in SAR time series, revealing important information about the area extent of internal accumulation. In the last application scenario, averaged summer SAR images were found to have potential in assisting the process of mapping glaciers outlines, especially in the presence of seasonal snow. Altogether we present examples of how to map glaciers and to further understand glaciological processes using the existing and future massive amount of multi-sensor time series data.

Soils in the northern high latitudes are a key component in the global carbon cycle, with potential feedback on climate. This study aims to improve the previous soil organic carbon (SOC) and total nitrogen (TN) storage estimates for the Zackenberg area (NE Greenland) that were based on a land cover classification (LCC) approach, by using geomorphological upscaling. In addition, novel organic carbon (OC) estimates for deeper alluvial and deltaic deposits (down to 300 cm depth) are presented. We hypothesise that landforms will better represent the long-term slope and depositional processes that result in deep SOC burial in this type of mountain permafrost environments. The updated mean SOC storage for the 0–100 cm soil depth is 4.8 kg C m−2, which is 42 % lower than the previous estimate of 8.3 kg C m−2 based on land cover upscaling. Similarly, the mean soil TN storage in the 0–100 cm depth decreased with 44 % from 0.50 kg (± 0.1 CI) to 0.28 (±0.1 CI) kg TN m−2. We ascribe the differences to a previous areal overestimate of SOC- and TN-rich vegetated land cover classes. The landform-based approach more correctly constrains the depositional areas in alluvial fans and deltas with high SOC and TN storage. These are also areas of deep carbon storage with an additional 2.4 kg C m−2 in the 100–300 cm depth interval. This research emphasises the need to consider geomorphology when assessing SOC pools in mountain permafrost landscapes.

The detachment of large parts of low-angle mountain glaciers resulting in massive ice–rock avalanches have so far been believed to be a unique type of event, made known to the global scientific community first for the 2002 Kolka Glacier detachment, Caucasus Mountains, and then for the 2016 collapses of two glaciers in the Aru range, Tibet. Since 2016, several so-far unrecognized low-angle glacier detachments have been recognized and described, and new ones have occurred. In the current contribution, we compile, compare, and discuss 20 actual or suspected large-volume detachments of low-angle mountain glaciers at 10 different sites in the Caucasus, the Pamirs, Tibet, Altai, the North American Cordillera, and the Southern Andes. Many of the detachments reached volumes in the order of 10–100 million m3. The similarities and differences between the presented cases indicate that glacier detachments often involve a coincidental combination of factors related to the lowering of basal friction, high or increasing driving stresses, concentration of shear stress, or low resistance to exceed stability thresholds. Particularly soft glacier beds seem to be a common condition among the observed events as they offer smooth contact areas between the glacier and the underlying substrate and are prone to till-strength weakening and eventually basal failure under high pore-water pressure. Partially or fully thawed glacier bed conditions and the presence of liquid water could thus play an important role in the detachments. Surface slopes of the detached glaciers range between around 10∘ and 20∘. This may be low enough to enable the development of thick and thus large-volume glaciers while also being steep enough to allow critical driving stresses to build up. We construct a simple slab model to estimate ranges of glacier slope and width above which a glacier may be able to detach when extensively losing basal resistance. From this model we estimate that all the detachments described in this study occurred due to a basal shear stress reduction of more than 50 %. Most of the ice–rock avalanches resulting from the detachments in this study have a particularly low angle of reach, down to around 5∘, likely due to their high ice content and connected liquefaction potential, the availability of soft basal slurries, and large amounts of basal water, as well as the smooth topographic setting typical for glacial valleys. Low-angle glacier detachments combine elements and likely also physical processes of glacier surges and ice break-offs from steep glaciers. The surge-like temporal evolution ahead of several detachments and their geographic proximity to other surge-type glaciers indicate the glacier detachments investigated can be interpreted as endmembers of the continuum of surge-like glacier instabilities. Though rare, glacier detachments appear to be more frequent than commonly thought and disclose, despite local differences in conditions and precursory evolutions, the fundamental and critical potential of low-angle soft glacier beds to fail catastrophically.

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

We present an updated, spatially resolved estimate of 2003–2008 glacier volume changes for entire High Mountain Asia (HMA) from ICESat laser altimetry data. The results reveal a diverse pattern that is driven 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. A step-increase in precipitation around 1997–2000 on the Tibetan Plateau (TP) caused thickening of glaciers in the Eastern Pamirs, Kunlun Shan and central TP by 0.1–0.7 m a−1. The thickening anomaly has a crisp boundary in the Eastern Pamir that continues just north of the central Karakoram. Glaciers in the south and east of the TP were thinning, with increasing rates towards southeast. The precipitation increase 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 SRTM DEM for over 1300 lakes. The rise in water volume contained in the lakes corresponds to 4–25 m 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 well in ERA Interim reanalysis data. Considering evaporation loss, the difference between average annual precipitation during the 1990s and 2000s suggested by these datasets is 34–100 m a−1, depending on region, which can fully explain both lake growth, and glacier thickening (Kunlun Shan) or glacier geometry changes (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, are 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.

Spatio-temporal patterns of rock glacier creep have rarely been studied outside the densely populated European Alps. This study investigates the spatial and temporal variability of rock glacier motion in the Ile Alatau and Kungöy Ala-Too mountain ranges, northern Tien Shan. Over the study region of more than 3000 km2, an inventory of slope movements is constructed using a large number of radar interferograms and high-resolution optical imagery. The inventory includes more than 900 landforms, of which around 550 are interpreted as rock glaciers. Out of the active rock glaciers, 45 are characterised by a rate of motion exceeding 1 m/a. From these fast rock glaciers we select six and study them in more detail (Gorodetzky, Morenny, Archaly, Ordzhonikidze, Karakoram and Kugalan Tash rock glaciers) using offset tracking between airphotos, and historical and modern very high resolution optical satellite data. Most of them show an overall increase of decadal surface velocities from the 1950s onwards with speeds being roughly two to three times higher in recent years compared to the 1950s and 1960s. This development indicates a possible significant increase in sediment and ice fluxes through rock glaciers and implies that – when compared to glacier shrinkage – periglacial sediment transport in the region seems to gain importance relative to glacial sediment transport. Those rock glacier fronts reaching the valley floors show a strongly compressive flow regime, and changes in speeds further upstream affect them only in a damped way. The only rock glacier investigated in detail that does not exhibit an overall increase in speed since the 1950s is Gorodetsky where glacier retreat and dead-ice degradation seem to have decoupled the rock glacier from its supply by glacial sediments and ice.

Snow consists of non-spherical grains of various shapes and sizes. Still, in radiative transfer calculations, snow grains are often treated as spherical. This also applies to the computation of snow albedo in the Snow, Ice, and Aerosol Radiation (SNICAR) model and in the Los Alamos sea ice model, version 4 (CICE4), both of which are employed in the Community Earth System Model and in the Norwegian Earth System Model (NorESM). In this study, we evaluate the effect of snow grain shape on climate simulated by NorESM in a slab ocean configuration of the model. An experiment with spherical snow grains (SPH) is compared with another (NONSPH) in which the snow shortwave single-scattering properties are based on a combination of three non-spherical snow grain shapes optimized using measurements of angular scattering by blowing snow. The key difference between these treatments is that the asymmetry parameter is smaller in the non-spherical case (0.77–0.78 in the visible region) than in the spherical case ( ≈ 0.89). Therefore, for the same effective snow grain size (or equivalently, the same specific projected area), the snow broadband albedo is higher when assuming non-spherical rather than spherical snow grains, typically by 0.02–0.03. Considering the spherical case as the baseline, this results in an instantaneous negative change in net shortwave radiation with a global-mean top-of-the-model value of ca. −0.22 W m−2. Although this global-mean radiative effect is rather modest, the impacts on the climate simulated by NorESM are substantial. The global annual-mean 2 m air temperature in NONSPH is 1.17 K lower than in SPH, with substantially larger differences at high latitudes. The climatic response is amplified by strong snow and sea ice feedbacks. It is further demonstrated that the effect of snow grain shape could be largely offset by adjusting the snow grain size. When assuming non-spherical snow grains with the parameterized grain size increased by ca. 70 %, the climatic differences to the SPH experiment become very small. Finally, the impact of assumed snow grain shape on the radiative effects of absorbing aerosols in snow is discussed.

In cold and arid climate, small glaciers with cold accumulation zone are often thought to be entirely cold based. However, Ground Penetrating Radar (GPR) measurements on Rikha Samba Glacier in the Nepal Himalaya reveal a large amount of temperate ice that seems to be influenced by the presence of crevassed areas. We used a coupled thermo-mechanical model forced by a firn model accounting for firn heating to interpret the observed thermal regime. We show that the addition of water percolation and refreezing in crevassed areas using a simple energy conservative approach is able to explain the observations. Model experiment shows that both steady and transient thermal regimes are significantly affected by latent heat release in crevassed areas. It makes half of the glacier base temperate, resulting in the dynamics mainly controlled by basal friction instead of ice deformation. Timescale of thermal regime change in response to warming climate is also greatly diminished with a potential switch from cold to temperate basal ice in 50–60 years in the upper part of the glacier while it would take 100–150 years without the crevasses effect. This study highlights the crucial role of water percolation through the crevasses on the thermal regime of glaciers and validates a simple method to take it into account in glacier thermo-mechanical models.

Spatio-temporal patterns related to the viscous creep in perennially frozen sediments of rock glaciers in cold mountains have rarely been studied outside the densely populated European Alps. This study investigates the spatial and temporal variability of rock glacier movement in the Ile Alatau and Kungöy Ala-Too mountain ranges, northern Tien Shan, a region with particularly large and fast rock glaciers. Over the study region of more than 3000 km2, an inventory of slope movements was constructed using a large number of radar interferograms and high-resolution optical imagery. The inventory includes more than 900 landforms, of which around 550 were interpreted as rock glaciers. Out of the active rock glaciers inventoried, 45 are characterized by a rate of motion exceeding 100 cm/a. From these fast rock glaciers we selected six (Gorodetzky, Morenny, Archaly, Ordzhonikidze, Karakoram, and Kugalan Tash) and studied them in more detail using offset tracking between repeat aerial images and historical and modern high-resolution optical satellite data. Two of these rock glaciers showed a steady increase in decadal surface velocities from the 1950s onwards, with speeds being roughly 2 to 4 times higher in recent years compared to the 1950s and 1960s. Three rock glaciers showed similar accelerations over the last 1 to 2 decades but also phases of increased speeds in the 1960s. This development indicates a possible significant increase in current sediment and ice fluxes through rock glaciers and implies that their material transport in the region might gain geomorphodynamic importance relative to material transport by glaciers, assuming the latter decreases together with the regional glacier shrinkage. The study demonstrates how air and satellite image archives are exploited to construct one of the longest decennial times series of rock glacier speeds currently available. Our results are in line with findings from Europe about rock glacier speeds increasing with atmospheric warming and underline local variability of such an overall response.

Recent estimates of the contribution of glaciers to sea-level rise during the 20th century are strongly divergent. Advances in data availability have allowed revisions of some of these published estimates. Here we show that outside of Antarctica, the global estimates of glacier mass change obtained from glacier-length-based reconstructions and from a glacier model driven by gridded climate observations are now consistent with each other, and also with an estimate for the years 2003–2009 that is mostly based on remotely sensed data. This consistency is found throughout the entire common periods of the respective data sets. Inconsistencies of reconstructions and observations persist in estimates on regional scales.