• Research dataclose
• Other research productsclose
• NO close
• English 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.

Acquiring data to analyse change in topography is often a costly endeavour requiring either extensive, potentially risky, fieldwork and/or expensive equipment or commercial data. Bringing the cost down while keeping the precision and accuracy has been a focus in geoscience in recent years. Structure from motion (SfM) photogrammetric techniques are emerging as powerful tools for surveying, with modern algorithm and large computing power allowing for the production of accurate and detailed data from low-cost, informal surveys. The high spatial and temporal resolution permits the monitoring of geomorphological features undergoing relatively rapid change, such as glaciers, moraines, or landslides. We present a method that takes advantage of light-transport flights conducting other missions to opportunistically collect imagery for geomorphological analysis. We test and validate an approach in which we attach a consumer-grade camera and a simple code-based Global Navigation Satellite System (GNSS) receiver to a helicopter to collect data when the flight path covers an area of interest. Our method is based and builds upon Welty et al. (2013), showing the ability to link GNSS data to images without a complex physical or electronic link, even with imprecise camera clocks and irregular time lapses. As a proof of concept, we conducted two test surveys, in September 2014 and 2015, over the glacier Midtre Lovénbreen and its forefield, in northwestern Svalbard. We were able to derive elevation change estimates comparable to in situ mass balance stake measurements. The accuracy and precision of our DEMs allow detection and analysis of a number of processes in the proglacial area, including the presence of thermokarst and the evolution of water channels.

Permafrost is a sensitive element of the cryosphere, but operational monitoring of the ground thermal conditions on large spatial scales is still lacking. Here, we demonstrate a remote-sensing-based scheme that is capable of estimating the transient evolution of ground temperatures and active layer thickness by means of the ground thermal model CryoGrid 2. The scheme is applied to an area of approximately 16 000 km2 in the Lena River delta (LRD) in NE Siberia for a period of 14 years. The forcing data sets at 1 km spatial and weekly temporal resolution are synthesized from satellite products and fields of meteorological variables from the ERA-Interim reanalysis. To assign spatially distributed ground thermal properties, a stratigraphic classification based on geomorphological observations and mapping is constructed, which accounts for the large-scale patterns of sediment types, ground ice and surface properties in the Lena River delta. A comparison of the model forcing to in situ measurements on Samoylov Island in the southern part of the study area yields an acceptable agreement for the purpose of ground thermal modeling, for surface temperature, snow depth, and timing of the onset and termination of the winter snow cover. The model results are compared to observations of ground temperatures and thaw depths at nine sites in the Lena River delta, suggesting that thaw depths are in most cases reproduced to within 0.1 m or less and multi-year averages of ground temperatures within 1–2 °C. Comparison of monthly average temperatures at depths of 2–3 m in five boreholes yielded an RMSE of 1.1 °C and a bias of −0.9 °C for the model results. The highest ground temperatures are calculated for grid cells close to the main river channels in the south as well as areas with sandy sediments and low organic and ice contents in the central delta, where also the largest thaw depths occur. On the other hand, the lowest temperatures are modeled for the eastern part, which is an area with low surface temperatures and snow depths. The lowest thaw depths are modeled for Yedoma permafrost featuring very high ground ice and soil organic contents in the southern parts of the delta. The comparison to in situ observations indicates that transient ground temperature modeling forced by remote-sensing data is generally capable of estimating the thermal state of permafrost (TSP) and its time evolution in the Lena River delta. The approach could hence be a first step towards remote detection of ground thermal conditions and active layer thickness in permafrost areas.

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 unknown glacier detachments have been discovered and described, and new ones occurred. In the current contribution, we compile, compare and discuss 19 actual or possible large-volume detachments of low-angle mountain glaciers at nine different sites in the Caucasus, the Pamirs, Tibet, Alaska’s St. Elias mountains, and the Southern Andes. Many of the detachments reached volumes in the order of 10–100 million m3. Commonalities and differences between the cases investigated suggest that a set of different conditions drives a transient combination of factors related to low basal friction, high driving stress, concentration of shear stress, and low resistance to exceed stability thresholds. Particularly, soft bedrocks below the detached glaciers seem to be a common condition among the observed events, as they offer smooth contact areas between the glacier and its substratum while being prone to till-strength weakening and eventually basal failure under high pore-water pressure. Surface slopes of the detached glaciers range between around 10° and 20°, possibly on the one hand low enough to enable development of thick and thus large-volume glaciers, and on the other hand steep enough to allow critical basal stresses to build up. Most of the ice-rock avalanches resulting from the detachments in this study have a particularly low angle of reach, down to around 0.1 (apparent friction angle), likely due to their high ice content and connected liquefaction potential, the ready availability of soft basal slurries and large amounts of basal water, and 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 or their geographic proximity to other surge-type glaciers suggests the glacier detachments investigated can be interpreted as end-members of the continuum of surge-like glacier instabilities. Though rare, glacier detachments appear more frequent than previously 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.

In this paper we describe a waves-in-ice model (WIM), which calculates ice breakage and the wave radiation stress (WRS). This WIM is then coupled to the new sea-ice model neXtSIM, which is based on the elasto-brittle (EB) rheology. We highlight some numerical issues involved in the coupling and investigate the impact of the WRS, and of modifying the EB rheology to lower the stiffness of the ice in the area where the ice has broken up (the marginal ice zone or MIZ). In experiments in the absence of wind, we find that wind waves can produce noticeable movement of the ice edge in loose ice (concentration around 70 %) – up to 36 km, depending on the material parameters of the ice that are used and the dynamical model used for the broken ice. The ice edge position is unaffected by the WRS if the initial concentration is higher (≳ 0.9). Swell waves (monochromatic waves with low frequency) do not affect the ice edge location (even for loose ice), as they are attenuated much less than the higher-frequency components of a wind wave spectrum, and so consequently produce a much lower WRS (by about an order of magnitude at least).In the presence of wind, we find that the wind stress dominates the WRS, which, while large near the ice edge, decays exponentially away from it. This is in contrast to the wind stress, which is applied over a much larger ice area. In this case (when wind is present) the dynamical model for the MIZ has more impact than the WRS, although that effect too is relatively modest. When the stiffness in the MIZ is lowered due to ice breakage, we find that on-ice winds produce more compression in the MIZ than in the pack, while off-ice winds can cause the MIZ to be separated from the pack ice.

The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting heat flux anomalies with various locations and intensities. In our model experiment, a minimum heat flux value of 970 mW m−2 located close to the East Greenland Ice-core Project (EGRIP) is required locally to reproduce the observed NEGIS velocities, giving basal melt rates consistent with previous estimates. The value cannot be attributed to geothermal heat flux alone and we suggest hydrothermal circulation as a potential explanation for the high local heat flux. By including high heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.

Effects of the short-term temporal variability of meteorological variables on soil temperature in northern high-latitude regions have been investigated. For this, a process-oriented land surface model has been driven using an artificially manipulated climate dataset. Short-term climate variability mainly impacts snow depth, and the thermal diffusivity of lichens and bryophytes. These impacts of climate variability on insulating surface layers together substantially alter the heat exchange between atmosphere and soil. As a result, soil temperature is 0.1 to 0.8 ∘C higher when climate variability is reduced. Earth system models project warming of the Arctic region but also increasing variability of meteorological variables and more often extreme meteorological events. Therefore, our results show that projected future increases in permafrost temperature and active-layer thickness in response to climate change will be lower (i) when taking into account future changes in short-term variability of meteorological variables and (ii) when representing dynamic snow and lichen and bryophyte functions in land surface models.