Description: Although ponds make up roughly half of the total area of surface water in permafrost landscapes, their relevance to carbon dioxide emissions on a landscape scale has, to date, remained largely unknown. We have therefore investigated the inflows and outflows of dissolved organic and inorganic carbon from lakes, ponds, and outlets on Samoylov Island, in the Lena Delta of northeastern Siberia in September 2008, together with their carbon dioxide emissions. Outgassing of carbon dioxide (CO2) from these ponds and lakes, which cover 25% of Samoylov Island, was found to account for between 74 and 81% of the calculated net landscape-scale CO2 emissions of 0.2–1.1 g C m�2 d�1 during September 2008, of which 28–43% was from ponds and 27–46% from lakes. The lateral export of dissolved carbon was negligible compared to the gaseous emissions due to the small volumes of runoff. The concentrations of dissolved inorganic carbon in the ponds were found to triple during freezeback, highlighting their importance for temporary carbon storage between the time of carbon production and its emission as CO2. If ponds are ignored the total summer emissions of CO2-C from water bodies of the islands within the entire Lena Delta (0.7–1.3 Tg) are underestimated by between 35 and 62%.

Description: We conducted eddy covariance measurements from April to August 2014 on a Siberian thermokarst lake. The study site is located in the Lena River Delta and characterized as a floating ice lake. Heat fluxes differed in magnitudes, directions and temporal patterns depending on the lake surface conditions (“frozen” ice cover, ice cover melt, and open water). Significant heat release during frozen ice cover conditions highlighted the importance of lakes for the landscape heat budget and water balance. The energy balance was nearly closed during the open water period and highlighted the impact of melting energy on its closure during the ice cover period. Sensible and latent heat dynamics were driven by temperature and water vapor gradients scaled by wind speed, respectively. We calculated bulk aerodynamics transfer coefficients and evaluated the performance of the derived in situ and three independent heat flux parameterization schemes. We found that bulk transfer models perform moderately to poorly for the different lake surface conditions. During the open water period small‐scale temporal variability could not be represented by the models, particularly in case of latent heat flux. The model results were less sensitive to the specific model type than to the accuracy of the surface water temperature measurement, which is dependent on a well‐thought‐out measurement design. Our study stresses considerations that are crucial for similar campaigns in the future, in order to face the measurement challenges encountered on arctic lakes especially during the ice cover period.

Description: In high‐latitude and mountain regions, local processes such as redistribution by wind, snow metamorphism and percolation of water, produce a complex spatial distribution of snow depths and snow densities. With its strong control on the ground thermal regime, this snow distribution has pronounced effects on ground temperatures at small spatial scales which are typically not resolved by land surface models (LSMs). This limits our ability to simulate the local impacts of climate change on for example vegetation and permafrost. Here, we present a tiling approach combining the CryoGrid permafrost model with snow microphysics parametrizations from the CROCUS snow scheme to account for sub‐grid lateral exchange of snow and water in a process‐based way. We demonstrate that a simple setup with three coupled tiles, each representing a different snow accumulation class with a specific topographic setting, can reproduce the observed spread of winter‐time ground surface temperatures (GST) and end‐of‐season snow distribution for a high‐Arctic site on Svalbard. For the three‐year study period, the three‐tile simulations showed substantial improvement compared to traditional single‐tile simulations which naturally cannot account for sub‐grid variability. Amongst others, the representation of the warmest and coldest 5% of the observed GST distribution was improved by 1‐2°C, while still capturing the average of the distribution. The simulations also reveal positive mean annual GSTs at the locations receiving the greatest snow cover. This could be an indication for the onset of localized permafrost degradation which would be obscured in single‐tile simulations.