CREp is able to calculate a large number of ages in a reasonable time (typically < 30 s for 50 samples). The user may export a summary table of the computed ages and the density probability function associated with each age (in the form of a spreadsheet).

Ice-wedge polygons are common features of lowland tundra in the continuous permafrost zone and prone to rapid degradation through melting of ground ice. There are many interrelated processes involved in ice-wedge thermokarst and it is a major challenge to quantify their influence on the stability of the permafrost underlying the landscape. In this study we used a numerical modelling approach to investigate the degradation of ice wedges with a focus on the influence of hydrological conditions. Our study area was Samoylov Island in the Lena River delta of northern Siberia, for which we had in situ measurements to evaluate the model. The tailored version of the CryoGrid 3 land surface model was capable of simulating the changing microtopography of polygonal tundra and also regarded lateral fluxes of heat, water, and snow. We demonstrated that the approach is capable of simulating ice-wedge degradation and the associated transition from a low-centred to a highcentred polygonal microtopography. The model simulations showed ice-wedge degradation under recent climatic conditions of the study area, irrespective of hydrological conditions. However, we found that wetter conditions lead to an earlier onset of degradation and cause more rapid ground subsidence. We set our findings in correspondence to observed types of ice-wedge polygons in the study area and hypothesized on remaining discrepancies between modelled and observed ice-wedge thermokarst activity. Our quantitative approach provides a valuable complement to previous, more qualitative and conceptual, descriptions of the possible pathways of ice-wedge polygon evolution. We concluded that our study is a blueprint for investigating thermokarst landforms and marks a step forward in understanding the complex interrelationships between various processes shaping ice-rich permafrost landscapes.

The ice-and organic-rich permafrost of the northeast Siberian Arctic lowlands (NESAL) has been projected to remain stable beyond 2100, even under pessimistic climate warming scenarios. However, the numerical models used for these projections lack processes which induce widespread landscape change termed thermokarst, precluding realistic simulation of permafrost thaw in such ice-rich terrain. Here, we consider thermokarst-inducing processes in a numerical model and show that substantial permafrost degradation, involving widespread landscape collapse, is projected for the NESAL under strong warming (RCP8.5), while thawing is moderated by stabilizing feedbacks under moderate warming (RCP4.5). We estimate that by 2100 thaw-affected carbon could be up to threefold (twelve-fold) under RCP4.5 (RCP8.5), of what is projected if thermokarst-inducing processes are ignored. Our study provides progress towards robust assessments of the global permafrost carbon-climate feedback by Earth system models, and underlines the importance of mitigating climate change to limit its impacts on permafrost ecosystems.

The Younger Dryas stadial, a cold event spanning 12,800 to 11,500 years ago, during the last deglaciation, is thought to coincide with the last major glacial re-advance in the tropical Andes. This interpretation relies mainly on cosmic-ray exposure dating of glacial deposits. Recent studies, however, have established new production rates for cosmogenic (10)Be and (3)He, which make it necessary to update all chronologies in this region and revise our understanding of cryospheric responses to climate variability. Here we present a new (10)Be moraine chronology in Colombia showing that glaciers in the northern tropical Andes expanded to a larger extent during the Antarctic cold reversal (14,500 to 12,900 years ago) than during the Younger Dryas. On the basis of a homogenized chronology of all (10)Be and (3)He moraine ages across the tropical Andes, we show that this behaviour was common to the northern and southern tropical Andes. Transient simulations with a coupled global climate model suggest that the common glacier behaviour was the result of Atlantic meridional overturning circulation variability superimposed on a deglacial increase in the atmospheric carbon dioxide concentration. During the Antarctic cold reversal, glaciers advanced primarily in response to cold sea surface temperatures over much of the Southern Hemisphere. During the Younger Dryas, however, northern tropical Andes glaciers retreated owing to abrupt regional warming in response to reduced precipitation and land-surface feedbacks triggered by a weakened Atlantic meridional overturning circulation. Conversely, glacier retreat during the Younger Dryas in the southern tropical Andes occurred as a result of progressive warming, probably influenced by an increase in atmospheric carbon dioxide. Considered with evidence from mid-latitude Andean glaciers, our results argue for a common glacier response to cold conditions in the Antarctic cold reversal exceeding that of the Younger Dryas.

Our paleoprecipitation map of the Altiplano reveals the atmospheric processes underlying a major hydroclimatic change 16 ka ago.

Abstract. Earth system models (ESMs) are our primary tool for projecting future climate change, but their ability to represent small-scale land surface processes is currently limited. This is especially true for permafrost landscapes in which melting of excess ground ice and subsequent subsidence affect lateral processes which can substantially alter soil conditions and fluxes of heat, water, and carbon to the atmosphere. Here we demonstrate that dynamically changing microtopography and related lateral fluxes of snow, water, and heat can be represented through a tiling approach suitable for implementation in large-scale models, and we investigate which of these lateral processes are important to reproduce observed landscape evolution. Combining existing methods for representing excess ground ice, snow redistribution, and lateral water and energy fluxes in two coupled tiles, we show that the model approach can simulate observed degradation processes in two very different permafrost landscapes. We are able to simulate the transition from low-centered to high-centered polygons, when applied to polygonal tundra in the cold, continuous permafrost zone, which results in (i) a more realistic representation of soil conditions through drying of elevated features and wetting of lowered features with related changes in energy fluxes, (ii) up to 2 ∘C reduced average permafrost temperatures in the current (2000–2009) climate, (iii) delayed permafrost degradation in the future RCP4.5 scenario by several decades, and (iv) more rapid degradation through snow and soil water feedback mechanisms once subsidence starts. Applied to peat plateaus in the sporadic permafrost zone, the same two-tile system can represent an elevated peat plateau underlain by permafrost in a surrounding permafrost-free fen and its degradation in the future following a moderate warming scenario. These results demonstrate the importance of representing lateral fluxes to realistically simulate both the current permafrost state and its degradation trajectories as the climate continues to warm. Implementing laterally coupled tiles in ESMs could improve the representation of a range of permafrost processes, which is likely to impact the simulated magnitude and timing of the permafrost–carbon feedback.

Peat plateaus and palsas are characteristic morphologies of sporadic permafrost, and the transition from permafrost to permafrost‐free ground typically occurs on spatial scales of meters. They are particularly vulnerable to climate change and are currently degrading in Fennoscandia. Here we present a spatially distributed data set of ground surface temperatures for two peat plateau sites in northern Norway for the year 2015–2016. Based on these data and thermal modeling, we investigate how the snow depth and water balance modulate the climate signal in the ground. We find that mean annual ground surface temperatures are centered around 2 to 2.5 °C for stable permafrost locations and 3.5 to 4.5 °C for permafrost‐free locations. The surface freezing degree days are characterized by a noticeable threshold around 200 °C.day, with most permafrost‐free locations ranging below this value and most stable permafrost ones above it. Freezing degree day values are well correlated to the March snow cover, although some variability is observed and attributed to the ground moisture level. Indeed, a zero curtain effect is observed on temperature time series for saturated soils during winter, while drained peat plateaus show early freezing surface temperatures. Complementarily, modeling experiments allow identifying a drainage effect that can modify 1‐m ground temperatures by up to 2 °C between drained and water accumulating simulations for the same snow cover. This effect can set favorable or unfavorable conditions for permafrost stability under the same climate forcing.

No abstract