Soil organic carbon (SOC) in Arctic coastal polygonal tundra is vulnerable to climate change, especially in soils with occurrence of large amounts of ground ice. Pan-arctic studies of mapping SOC exist, yet they fail to describe the high spatial variability of SOC storage in permafrost landscapes. An important factor is the landscape history which determines landform development and consequently the spatial variability of SOC. Our aim was to map SOC stocks, and which environmental variables that determine SOC, in two adjacent coastal areas along Canadian Beaufort Sea coast with different glacial history. We used the machine learning technique random forest and environmental variables to map the spatial distribution of SOC stocks down to 1 m depth at a spatial resolution of 2 m for depth increments of 0-5, 5-15, 15-30, 30-60 and 60-100 cm. The results show that the two study areas had large differences in SOC stocks in the depth 60-100 cm due to high amounts of ground ice in one of the study areas. There are also differences in variable importance of the explanatory variables between the two areas. The area low in ground ice content had with 66.6 kg C/m(-2) more stored SOC than the area rich in ground ice content with 40.0 kg C/m(-2). However, this SOC stock could be potentially more vulnerable to climate change if ground ice melts and the ground subsides. The average N stock of the area low in ground ice is 3.77 kg m(-2) and of the area rich in ground ice is 3.83 kg m(-2). These findings support that there is a strong correlation between ground ice and SOC, with less SOC in ice-rich layers on a small scale. In addition to small scale studies of SOC mapping, detailed maps of ground ice content and distribution are needed for a validation of large-scale quantifications of SOC stocks and transferability of models.

Ice-wedge polygon landscapes make up a substantial part of high-latitude permafrost landscapes. The hydrological conditions shape how these landscapes store and release organic carbon. However, their coupled water‑carbon dynamics are poorly understood as field measurements are sparse in smaller catchments and coupled hydrology-dissolved organic carbon (DOC) models are not tailored for these landscapes. Here we present a model that simulates the hydrology and associated DOC export of high-centered and low-centered ice-wedge polygons and apply the model to a small catchment with abundant polygon coverage along the Yukon Coast, Canada. The modeled seasonal pattern of water and carbon fluxes aligns with sparse field data. These modeled seasonal patterns indicate that early-season runoff is mostly surficial and generated by low-centered polygons and snow trapped in troughs of high-centered polygons. High-centered polygons show potential for deeper subsurface flow under future climate conditions. This suggests that high-centered polygons will be responsible for an increasing proportion of annual DOC export compared to low-centered polygons. Warming likely shifts low-centered polygons to high-centered polygons, and our model shows that this shift will cause a deepening of the active layer and a lengthening of the thawing season. This, in turn, intensifies seasonal runoff and DOC flux, mainly through its duration. Our model provides a physical hypothesis that can be used to further quantify and refine our understanding of hydrology and DOC export of arctic ice-wedge polygon terrain.