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.

Abstract. Permafrost degradation in steep rock walls and associated slope destabilization have been studied increasingly in recent years. While most studies focus on mountainous and sub-Arctic regions, the occurring thermo-mechanical processes also play an important role in the high Arctic. A more precise understanding is required to assess the risk of natural hazards enhanced by permafrost warming in high-Arctic rock walls. This study presents one of the first comprehensive datasets of rock surface temperature measurements of steep rock walls in the high Arctic, comparing coastal and near-coastal settings. We applied the surface energy balance model CryoGrid 3 for evaluation, including adjusted radiative forcing to account for vertical rock walls. Our measurements comprise 4 years of rock surface temperature data from summer 2016 to summer 2020. Mean annual rock surface temperatures ranged from −0.6 in a coastal rock wall in 2017/18 to −4.3 ∘C in a near-coastal rock wall in 2019/20. Our measurements and model results indicate that rock surface temperatures at coastal cliffs are up to 1.5 ∘C higher than at near-coastal rock walls when the fjord is ice-free in winter, resulting from additional energy input due to higher air temperatures at the coast and radiative warming by relatively warm seawater. An ice layer on the fjord counteracts this effect, leading to similar rock surface temperatures to those in near-coastal settings. Our results include a simulated surface energy balance with shortwave radiation as the dominant energy source during spring and summer with net average seasonal values of up to 100 W m−2 and longwave radiation being the main energy loss with net seasonal averages between 16 and 39 W m−2. While sensible heat fluxes can both warm and cool the surface, latent heat fluxes are mostly insignificant. Simulations for future climate conditions result in a warming of rock surface temperatures and a deepening of active layer thickness for both coastal and near-coastal rock walls. Our field data present a unique dataset of rock surface temperatures in steep high-Arctic rock walls, while our model can contribute towards the understanding of factors influencing coastal and near-coastal settings and the associated surface energy balance.

Abstract. Ground temperatures in coarse, blocky deposits such as mountain blockfields and rock glaciers have long been observed to be lower in comparison with other (sub)surface material. One of the reasons for this negative temperature anomaly is the lower soil moisture content in blocky terrain, which decreases the duration of the zero curtain in autumn. Here we used the CryoGrid community model to simulate the effect of drainage on the ground thermal regime and ground ice in blocky terrain permafrost at two sites in Norway. The model set-up is based on a one-dimensional model domain and features a surface energy balance, heat conduction and advection, as well as a bucket water scheme with adjustable lateral drainage. We used three idealized subsurface stratigraphies, blocks only, blocks with sediment and sediment only, which can be either drained (i.e. with strong lateral subsurface drainage) or undrained (i.e. without drainage), resulting in six scenarios. The main difference between the three stratigraphies is their ability to retain water against drainage: while the blocks only stratigraphy can only hold small amounts of water, much more water is retained within the sediment phase of the two other stratigraphies, which critically modifies the freeze–thaw behaviour. The simulation results show markedly lower ground temperatures in the blocks only, drained scenario compared to other scenarios, with a negative thermal anomaly of up to 2.2 ∘C. For this scenario, the model can in particular simulate the time evolution of ground ice, with build-up during and after snowmelt and spring and gradual lowering of the ice table in the course of the summer season. The thermal anomaly increases with larger amounts of snowfall, showing that well-drained blocky deposits are less sensitive to insulation by snow than other soils. We simulate stable permafrost conditions at the location of a rock glacier in northern Norway with a mean annual ground surface temperature of 2.0–2.5 ∘C in the blocks only, drained simulations. Finally, transient simulations since 1951 at the rock glacier site (starting with permafrost conditions for all stratigraphies) showed a complete loss of perennial ground ice in the upper 5 m of the ground in the blocks with sediment, drained run; a 1.6 m lowering of the ground ice table in the sediment only, drained run; and only 0.1 m lowering in the blocks only, drained run. The interplay between the subsurface water–ice balance and ground freezing/thawing driven by heat conduction can at least partly explain the occurrence of permafrost in coarse blocky terrain below the elevational limit of permafrost in non-blocky sediments. It is thus important to consider the subsurface water–ice balance in blocky terrain in future efforts in permafrost distribution mapping in mountainous areas. Furthermore, an accurate prediction of the evolution of the ground ice table in a future climate can have implications for slope stability, as well as water resources in arid environments.

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 parameterizations, 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.

Abstract. Ground temperatures in coarse, blocky deposits such as mountain blockfields and rock glaciers have long been observed to be lower in comparison with other (sub)surface material. One of the reasons for this negative temperature anomaly is the lower soil moisture content in blocky terrain, which decreases the duration of the zero curtain in autumn. Here we used the CryoGrid community model to simulate the effect of drainage in blocky terrain permafrost at two sites in Norway. The model setup features a surface energy balance, heat conduction and advection, as well as a bucket water scheme with adjustable lateral drainage. We used three idealized subsurface stratigraphies, denoted blocks only, blocks with sediment and sediment only and are either drained or undrained of water, resulting in six ‘scenarios’. The main difference between the three stratigraphies is their ability to retain water against drainage: while the blocks only stratigraphy can only hold small amounts of water, much more water is retained within the sediment phase of the two other stratigraphies, which critically modifies the freeze-thaw behaviour. The simulation results show markedly lower ground temperatures in the blocks only, drained scenario compared to other scenarios, with negative thermal anomaly of up to 1.8–2.2 °C. For this scenario, the model can in particular simulate the time evolution of ground ice, with build-up during and after snow melt and spring and gradual lowering of the ice table in the course of the summer season. We simulate stable permafrost conditions at the location of a rock glacier in northern Norway with a mean annual ground surface temperature of 2.0–2.5 °C in the blocks only, drained simulations. Finally, transient simulations at the rock glacier site showed a complete or partial lowering of the ground ice table since 1951 for all simulations except the blocks only, drained run. The interplay between the subsurface water/ice balance and ground freezing/thawing driven by heat conduction can at least partly explain the occurrence of permafrost in coarse blocky terrain below the assumed elevational limit of permafrost. It is thus important to consider this effect in future efforts on permafrost distribution mapping in mountainous areas. Furthermore, an accurate prediction of the evolution of the ground ice table in a future climate can have implications for slope stability, as well as water resources in arid environments.