Abstract. The CryoGrid community model is a flexible toolbox for simulating the ground thermal regime and the ice/water balance for permafrost and glaciers, extending a well-established suite of permafrost models (CryoGrid 1, 2 and 3). The CryoGrid community model can accommodate a wide variety of application scenarios, which is achieved by fully modular structures through object-oriented programming. Different model components, characterized by their process representations and parametrizations, are realized as classes (i.e. objects) in CryoGrid. Standardized communication protocols between these classes ensure that they can be stacked vertically. For example, the CryoGrid community model features several classes with different complexity for the seasonal snow cover which can be flexibly combined with a range of classes representing subsurface materials, each with their own set of process representations (e.g. soil with and without water balance, glacier ice). We present the CryoGrid architecture as well as the model physics and defining equations for the different model classes, focusing on one-dimensional model configurations which can also interact with external heat and water reservoirs. We illustrate the wide variety of simulation capabilities for a site on Svalbard, with point-scale permafrost simulations using e.g. different soil freezing characteristics, drainage regimes and snow representations, as well as simulations for glacier mass balance and a shallow water body. The CryoGrid community model is not intended as a static model framework, but aims to provide developers with a flexible platform for efficient model development. In this study, we document both basic and advanced model functionalities to provide a baseline for the future development of novel cryosphere models.

No abstract

Abstract. Understanding the future evolution of permafrost requires a better understanding of its climatological past. This requires permafrost models to efficiently simulate the thermal dynamics of permafrost over the past centuries to millennia, taking into account highly uncertain soil and snow properties. In this study, we present a computationally efficient numerical permafrost model which satisfactorily reproduces the current thermal state of permafrost in the Arctic and its recent trend over the last decade. Also, the active layer dynamics and its trend is realistically captured. The performed simulations provide insights into the evolution of permafrost since the 18th century and show that permafrost on the North American continent is subject to early degradation, while permafrost on the Eurasian continent is relatively stable over the investigated 300-year period. Permafrost warming since industrialization has occurred primarily in three "hotspot" regions in northeastern Canada, northern Alaska, and, to a lesser extent, western Siberia. The extent of near-surface permafrost has changed substantially since the 18th century. In particular, loss of continuous permafrost has accelerated from low (−0.10 × 105 km2 dec−1) to moderate (−0.77 × 105 km2 dec−1) rates for the 18th and 19th centuries, respectively. In the 20th century, the loss rate nearly doubled (−1.36 × 105 km2 dec−1), with the highest near-surface permafrost losses occurring in the last 50 years. Our simulations further indicate that climate disturbances due to large volcanic eruptions in the Northern Hemisphere, can only counteract near-surface permafrost loss for a relatively short period of a few decades. Despite some limitations, the presented model shows great potential for further investigation of the climatological past of permafrost, especially in conjunction with paleoclimate modeling.

Abstract. Long-term measurements of permafrost temperatures do not provide a complete picture of the Arctic subsurface thermal regime. Regions with warmer permafrost often show little to no long-term change in ground temperature due to the uptake and release of latent heat during freezing and thawing. Thus, regions where the least warming is observed may also be the most vulnerable to permafrost degradation. Since direct measurements of ice and liquid water contents in the permafrost layer are not widely available, thermal modeling of the subsurface plays a crucial role in understanding how permafrost responds to changes in the local energy balance. In this work, we first analyze trends in observed air and permafrost temperatures at four sites within the continuous permafrost zone, where we find substantial variation in the apparent relationship between long-term changes in permafrost temperatures (0.02 K yr−1 to 0.16 K yr−1) and air temperature (0.09 K yr−1 to 0.11 K yr−1). We then apply recently developed Bayesian inversion methods to link observed changes in borehole temperatures to unobserved changes in latent heat and thaw depth using a transient model of heat conduction with phase change. Our results suggest that the degree to which recent warming trends correlate with permafrost thaw and variations in latent heat is heavily dependent on both local soil properties as well as historical climatology. At the warmest site, a nine meter borehole near Ny-Ålesund, Svalbard, modeled annual maximum thaw depth increases by an average of (12 ± 1) cm K−1 rise in mean annual ground temperature. In stark contrast, modeled thaw rates for a borehole on Samoylov Island in the Lena River Delta (northeastern Siberia) appear far less sensitive to temperature change, with an almost negligible increase of (1 ± 1) cm K−1. Although our study is limited to just four sites, the results urge caution in the interpretation and comparison of warming trends in Arctic boreholes, indicating substantial uncertainty in their implications for the current and future thermal state of permafrost.