Snow accumulation, albedo and melt patterns following road construction on permafrost, Inuvik–Tuktoyaktuk Highway, Canada

Roads constructed on permafrost can have a significant impact on the surrounding environment, potentially inducing permafrost degradation. These impacts arise from factors such as snow accumulation near the road, which affects the soil's thermal and hydrological regime, and road dust that decreases the snow's albedo, altering the timing of snowmelt. However, our current understanding of the magnitude and the spatial extent of these effects is limited. In this study we addressed this gap by using remote sensing techniques to assess the spatial effect of the Inuvik to Tuktoyaktuk Highway (ITH) in Northwest Territories, Canada, on snow accumulation, snow albedo and snowmelt patterns. With a new, high resolution snow depth raster from airborne laser scanning, we quantified the snow accumulation at road segments in the Trail Valley Creek area using digital elevation model differencing. We found increased snow accumulation up to 36 m from the road center. The magnitude of this snow accumulation was influenced by the prevailing wind direction and the embankment height. Furthermore, by analyzing 43 Sentinel-2 satellite images between February and May 2020, we observed reduced snow albedo values within 500 m of the road, resulting in a 12-days-earlier onset of snowmelt within 100 m from the road. We examined snowmelt patterns before, during and after the road construction using the normalized difference snow index from Landsat-7 and Landsat-8 imagery. Our analysis revealed that the road affected the snowmelt pattern up to 600 m from the road, even in areas which appeared undisturbed. In summary, our study improves our understanding of the spatial impact of gravel roads on permafrost due to enhanced snow accumulation, reduced snow albedo and earlier snowmelt. Our study underscores the important contribution that remote sensing can provide to improve our understanding of the effects of infrastructure development on permafrost environments.</p

Simulating ice segregation and thaw consolidation in permafrost environments with the CryoGrid community model

The ground ice content in cold environments influences the permafrost thermal regime and the thaw trajectories in a warming climate, especially for soils containing excess ice. Despite their importance, the amount and distribution of ground ice are often unknown due to lacking field observations. Hence, modeling the thawing of ice-rich permafrost soils and associated thermokarst is challenging as ground ice content has to be prescribed in the model setup. In this study, we present a model scheme, capable of simulating segregated ice formation during a model spinup together with associated ground heave. It provides the option to add a constant sedimentation rate throughout the simulation. Besides ice segregation, it can represent thaw consolidation processes and ground subsidence under a warming climate. The computation is based on soil mechanical processes, soil hydrology by the Richards equation and soil freezing characteristics. The code is implemented in the CryoGrid community model (version 1.0), a modular land surface model for simulations of the ground thermal regime. The simulation of ice segregation and thaw consolidation with the new model scheme allows us to analyze the evolution of ground ice content in both space and time. To do so, we use climate data from two contrasting permafrost sites to run the simulations. Several influencing factors are identified, which control the formation and thaw of segregated ice. (i) Model results show that high temperature gradients in the soil as well as moist conditions support the formation of segregated ice. (ii) We find that ice segregation increases in fine-grained soils and that especially organic-rich sediments enhance the process. (iii) Applying external loads suppresses ice segregation and speeds up thaw consolidation. (iv) Sedimentation leads to a rise of the ground surface and the formation of an ice-enriched layer whose thickness increases with sedimentation time. We conclude that the new model scheme is a step forward to improve the description of ground ice distributions in permafrost models and can contribute towards the understanding of ice segregation and thaw consolidation in permafrost environments under changing climatic conditions.</p

Observation and modelling of snow at a polygonal tundra permafrost site: spatial variability and thermal implications

The shortage of information on snow properties in high latitudes places a major limitation on permafrost and more generally climate modelling. A dedicated field program was therefore carried out to investigate snow properties and their spatial variability at a polygonal tundra permafrost site. Notably, snow samples were analysed for surface-normal thermal conductivity (Keff − z) based on X-ray microtomography. Also, the detailed snow model SNOWPACK was adapted to these Arctic conditions to enable relevant simulations of the ground thermal regime. Finally, the sensitivity of soil temperatures to snow spatial variability was analysed.Within a typical tundra snowpack composed of depth hoar overlain by wind slabs, depth hoar samples were found more conductive (Keff − z = 0.22±0.05&thinsp;W&thinsp;m−1&thinsp;K−1) than in most previously published studies, which could be explained by their high density and microstructural anisotropy. Spatial variations in the thermal properties of the snowpack were well explained by the microtopography and ground surface conditions of the polygonal tundra, which control depth hoar growth and snow accumulation. Our adaptations to SNOWPACK, phenomenologically taking into account the effects of wind compaction, basal vegetation, and water vapour flux, yielded realistic density and Keff − z profiles that greatly improved simulations of the ground thermal regime. Also, a density- and anisotropy-based parameterization for Keff − z lead to further slight improvements. Soil temperatures were found to be particularly sensitive to snow conditions during the early winter and polar night, highlighting the need for improved snow characterization and modelling over this period.</p

Thaw processes in ice-rich permafrost landscapes represented with laterally coupled tiles in a land surface model

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&thinsp;∘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.</p

Lena Delta hydrology and geochemistry: long-term hydrological data and recent field observations

The Lena River forms one of the largest deltas in the Arctic. We compare two sets of data to reveal new insights into the hydrological, hydrochemical, and geochemical processes within the delta: (i) long-term hydrometric observations at the Khabarova station at the head of the delta from 1951 to 2005; (ii) field hydrological and geochemical observations carried out within the delta since 2002. Periods with differing relative discharge and intensity of fluvial processes were identified from the long-term record of water and sediment discharge. Ice events during spring melt (high water) reconfigured branch channels and probably influenced sediment transport within the delta. Based on summer field measurements during 2005–2012 of discharge and sediment fluxes along main delta channels, both are increased between the apex and the front of the delta. This increase is to a great extent connected with an additional influx of water from tributaries, as well as an increase of suspended and dissolved material released from the ice complex. Summer concentrations of major ion and biogenic substances along the delta branches are partly explained by water sources within the delta, such as thawing ice complex waters, small Lena River branches and estuarine areas

Explicitly modelling microtopography in permafrost landscapes in a land surface model (JULES vn5.4_microtopography)

Microtopography can be a key driver of heterogeneity in the ground thermal and hydrological regime of permafrost landscapes. In turn, this heterogeneity can influence plant communities, methane fluxes, and the initiation of abrupt thaw processes. Here we have implemented a two-tile representation of microtopography in JULES (the Joint UK Land Environment Simulator), where tiles are representative of repeating patterns of elevation difference. Tiles are coupled by lateral flows of water, heat, and redistribution of snow, and a surface water store is added to represent ponding. Simulations are performed of two Siberian polygon sites, (Samoylov and Kytalyk) and two Scandinavian palsa sites (Stordalen and Iškoras). The model represents the observed differences between greater snow depth in hollows vs. raised areas well. The model also improves soil moisture for hollows vs. the non-tiled configuration (“standard JULES”) though the raised tile remains drier than observed. The modelled differences in snow depths and soil moisture between tiles result in the lower tile soil temperatures being warmer for palsa sites, as in reality. However, when comparing the soil temperatures for July at 20 cm depth, the difference in temperature between tiles, or “temperature splitting”, is smaller than observed (3.2 vs. 5.5 ∘C). Polygons display small (0.2 ∘C) to zero temperature splitting, in agreement with observations. Consequently, methane fluxes are near identical (+0 % to 9 %) to those for standard JULES for polygons, although they can be greater than standard JULES for palsa sites (+10 % to 49 %). Through a sensitivity analysis we quantify the relative importance of model processes with respect to soil moisture and temperatures, identifying which parameters result in the greatest uncertainty in modelled temperature. Varying the palsa elevation between 0.5 and 3 m has little effect on modelled soil temperatures, showing that using only two tiles can still be a valid representation of sites with a range of palsa elevations. Mire saturation is heavily dependent on landscape-scale drainage. Lateral conductive fluxes, while small, reduce the temperature splitting by ∼ 1 ∘C and correspond to the order of observed lateral degradation rates in peat plateau regions, indicating possible application in an area-based thaw model