Characteristics and changes of thermokarst lakes in Canada, Russia and Sweden

Peat plateaus are widespread at high northern latitudes and are important soil organic carbon reservoirs. A warming climate can cause either increased ground subsidence (thermokarst) resulting in lake formation or increased drainage as the permafrost thaws. A better understanding of spatiotemporal variations in these landforms in relation to climate change is important for predicting the future thawing permafrost carbon feedback. In this study, dynamics in thermokarst lake extent during the last 35-50 years has been quantified through time series analysis of aerial photographs and high-resolution satellite images (IKONOS/QuickBird) in three peat plateau complexes, spread out across the northern circumpolar region along a climatic and permafrost gradient. From the mid-1970s until the mid-2000s there has been an increase in mean annual air temperature, winter precipitation, and ground temperature in all three study areas. The two peat plateaus located in the continuous and discontinuous permafrost zones, respectively, where mean annual air temperatures are below -5°C and ground temperatures are -2°C or colder, have experienced small changes in thermokarst lake extent. In the peat plateau located in the sporadic permafrost zone where the mean annual air temperature is around -3°C, and the ground temperature is close to 0°C, lake drainage and infilling with fen vegetation has been extensive and many new thermokarst lakes have formed. In a future progressively warmer and wetter climate permafrost degradation can cause significant impacts on landscape composition and greenhouse gas exchange also in areas with extensive peat plateaus, which presently still experience stable permafrost conditions.

Temperature, water level and bathymetry of thermokarst lakes in the continuous permafrost zone of northern Siberia - Lena River Delta, Siberia

Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009-2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths > 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W/m**2 in summer and with heat released back into the water column at a rate of less than 1 W/m**2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Spatial distribution and morphometric analysis of thermokarst lakes and other water bodies : Case study from the sporadic permafrost region, Tavvavuoma, Sweden

Projected air and ground temperatures are expected to be higher in Arctic and sub-Arcticlatitudes and with temperatures already close to the limit where permafrost can exist,resistance against degradation is low. With thawing permafrost, the landscape is modifiedwith depression in which thermokarst lakes emerge. In permafrost soils a considerableamount of soil organic carbon is stored, with the potential of altering climate even furtherif expansion and formation of new thermokarst lakes emerge, as decay releasesgreenhouse gases (C02 and CH4) to the atmosphere. Analyzing the spatial distribution andmorphometry over time of thermokarst lakes and other water bodies, is of importance inaccurately predict carbon budget and feedback mechanisms, as well as to assess futurelandscape layout and these features interaction. Different types of high-spatial resolutionaerial and satellite imageries from 1963, 1975, 2003, 2010 and 2015, were used in bothpre- and post-classification change detection analyses. Using object oriented segmentationin eCognition combined with manual adjustments, resulted in digitalized water bodies>28m2 from which direction of change and morphometric values were extracted. Thequantity of thermokarst lakes and other water bodies was in 1963 n=92, with succeedingyears as a trend decreased in numbers, until 2010-2015 when eleven water bodies wereadded in 2015 (n=74 to n=85). In 1963-2003, area of these water bodies decreased with50 651m2 (189 446-138 795m2) and continued to decrease in 2003-2015 ending at 129337m2. Limnicity decreased from 19.9% in 1963 to 14.6% in 2003 (-5.3%). In 2010 and2015 13.7-13.6%. The late increase in water bodies differs from an earlier hypothesis thatsporadic permafrost regions experience decrease in both area and quantity of thermokarstlakes and water bodies. During 1963-2015, land gain has been in dominance of the ratiobetween the two competing processes of expansion and drainage. In 1963-1975, 55/45%,followed by 90/10% in 1975-2003. After major drainage events, land loss increased to62/38% in 2010-2015. Drainage and infilling rates, calculated for 15 shorelines werevaried across both landscape and parts of shorelines, with in average 0.17/0.15/0.14m/yr.Except for 1963-1975 when rate of change in average was in opposite direction (-0.09m/yr.), likely due to evident expansion of a large thermokarst lake. Using a squaregrid, distribution of water bodies was determined, with an indistinct cluster located in NEand central parts. Especially for water bodies <250m2, which is the dominant area classthroughout 1963-2015 ranging from n=39-51. With a heterogeneous composition of bothsmall and large thermokarst lakes, and with both expansion and drainage altering thelandscape in Tavvavuoma, both positive and negative climate feedback mechanisms are inplay - given that sporadic permafrost still exist.

Evidence for past variations in methane availability in a Siberian thermokarst lake based on δ13C of chitinous invertebrate remains

Understanding past methane dynamics in arctic wetlands and lakes is crucial for estimating future methane release. Methane fluxes from lake ecosystems have increasingly been studied, yet only few reconstructions of past methane emissions from lakes are available. In this study, we develop an approach to assess changes in methane availability in lakes based on δ13C of chitinous invertebrate remains and apply this to a sediment record from a Siberian thermokarst lake. Diffusive methane fluxes from the surface of ten newly sampled Siberian lakes and seven previously studied Swedish lakes were compared to taxon-specific δ13C values of invertebrate remains from lake surface sediments to investigate whether these invertebrates assimilated 13C-depleted carbon typical for methane. Remains of chironomid larvae of the tribe Orthocladiinae that, in the study lakes, mainly assimilate plant-derived carbon had higher δ13C than other invertebrate groups. δ13C of other invertebrates such as several chironomid groups (Chironomus, Chironomini, Tanytarsini, and Tanypodinae), cladocerans (Daphnia), and ostracods were generally lower. δ13C of Chironomini and Daphnia, and to a lesser extent Tanytarsini was variable in the lakes and lower at sites with higher diffusive methane fluxes. δ13C of Chironomini, Tanytarsini, and Daphnia were correlated significantly with diffusive methane flux in the combined Siberian and Swedish dataset (r = −0.72, p = 0.001, r = −0.53, p = 0.03, and r = −0.81, p < 0.001, respectively), suggesting that δ13C in these invertebrates was affected by methane availability. In a second step, we measured δ13C of invertebrate remains from a sediment record of Lake S1, a shallow thermokarst lake in northeast Siberia. In this record, covering the past ca 1000 years, δ13C of taxa most sensitive to methane availability (Chironomini, Tanytarsini, and Daphnia) was lowest in sediments deposited from ca AD 1250 to ca AD 1500, and after AD 1970, coinciding with warmer climate as indicated by an independent local temperature record. As a consequence the offset in δ13C between methane-sensitive taxa and bulk organic matter was higher in these sections than in other parts of the core. In contrast, δ13C of other invertebrate taxa did not show this trend. Our results suggest higher methane availability in the study lake during warmer periods and that thermokarst lakes can respond dynamically in their methane output to changing environmental conditions.

Isotopic composition of ground ice, ebullition gases and thermokarst lake water, North America

Thermokarst lakes are thought to have been an important source of methane (CH4) during the last deglaciation when atmospheric CH4 concentrations increased rapidly. Here we demonstrate that meltwater from permafrost ice serves as an H source to CH4 production in thermokarst lakes, allowing for region-specific reconstructions of dD-CH4 emissions from Siberian and North American lakes. dD CH4 reflects regionally varying dD values of precipitation incorporated into ground ice at the time of its formation. Late Pleistocene-aged permafrost ground ice was the dominant H source to CH4 production in primary thermokarst lakes, whereas Holocene-aged permafrost ground ice contributed H to CH4 production in later generation lakes. We found that Alaskan thermokarst lake dD-CH4 was higher (-334 ± 17 per mil) than Siberian lake dD-CH4 (-381 ± 18 per mil). Weighted mean dD CH4 values for Beringian lakes ranged from -385 per mil to -382 per mil over the deglacial period. Bottom-up estimates suggest that Beringian thermokarst lakes contributed 15 ± 4 Tg CH4 /yr to the atmosphere during the Younger Dryas and 25 ± 5 Tg CH4 /yr during the Preboreal period. These estimates are supported by independent, top-down isotope mass balance calculations based on ice core dD-CH4 and d13C-CH4 records. Both approaches suggest that thermokarst lakes and boreal wetlands together were important sources of deglacial CH4.

Evolution of thermokarst in East Siberian ice-rich permafrost

Thermokarst lakes and basins are major components of ice-rich permafrost landscapes in East Siberian coastal lowlands and are regarded as indicators of regional climatic changes. We investigate the temporal and spatial dynamics of a 7.5 km**2, partly drained thermokarst basin (alas) using field investigations, remote sensing, Geographic Information Systems (GIS), and sediment analyses. The evolution of the thermokarst basin proceeded in two phases. The first phase started at the Pleistocene/Holocene transition (13 to 12 ka BP) with the initiation of a primary thermokarst lake on the Ice Complex surface. The lake expanded and persisted throughout the early Holocene before it drained abruptly about 5.7 ka BP, thereby creating a > 20 m deep alas with residual lakes. The second phase (5.7 ka BP to present) is characterized by alternating stages of lower and higher thermokarst intensity within the alas that were mainly controlled by local hydrological and relief conditions and accompanied by permafrost aggradation and degradation. It included diverse concurrent processes like lake expansion and stepwise drainage, polygonal ice-wedge growth, and the formation of drainage channels and a pingo, which occurred in different parts of the alas. This more dynamic thermokarst evolution resulted in a complex modern thermokarst landscape. However, on the regional scale, the changes during the second evolutionary phase after drainage of the initial thermokarst lakes were less intense than the early Holocene extensive thermokarst development in East Siberian coastal lowlands as a result of a significant regional change to warmer and wetter climate conditions.

Physico-chemical parameters of sediment core PG1972-1 (09-Tik-03) from a Tundra Lake, Northeastern Siberia

Thermokarst lakes in the Siberian Arctic contain sediment archives that can be used for paleoenvironmental inference. Until now, however, there has been no study from the inner Lena River Delta with a focus on diatoms. The objective of this study was to investigate how the diatom community in a thermokarst lake responded to past limnogeological changes and what specific factors drove variations in the diatom assemblage. We analysed fossil diatom species, organic content, grain-size distribution and elemental composition in a sediment core retrieved in 2009 from a shallow thermokarst lake in the Arga Complex, western Lena River Delta. The core contains a 3,000-year record of sediment accumulation. Shifts in the predominantly benthic and epiphytic diatom species composition parallel changes in sediment characteristics. Paleoenvironmental and limnogeological development, inferred from multiple biological and sedimentological variables, are discussed in the context of four diatom zones, and indicate a strong relation between changes in the diatom assemblage and thermokarst processes. We conclude that limnogeological and thermokarst processes such as lake drainage, rather than direct climate forcing, were the main factors that altered the aquatic ecosystem by influencing, for example, habitat availability, hydrochemistry, and water level.

Age determination, sedimentological, biogeochemical, and plant macrofossils analyses of outcrop PG2038-1

Thermokarst lakes are a widespread feature of the Arctic tundra, in which highly dynamic processes are closely connected with current and past climate changes. We investigated late Quaternary sediment dynamics, basin and shoreline evolution, and environmental interrelations of Lake El'gene-Kyuele in the NE Siberian Arctic (latitude 71°17'N, longitude 125°34'E). The water-body displays thaw-lake characteristics cutting into both Pleistocene Ice Complex and Holocene alas sediments. Our methods are based on grain size distribution, mineralogical composition, TOC/N ratio, stable carbon isotopes and the analysis of plant macrofossils from a 3.5-m sediment profile at the modern eastern lake shore. Our results show two main sources for sediments in the lake basin: terrigenous diamicton supplied from thermokarst slopes and the lake shore, and lacustrine detritus that has mainly settled in the deep lake basin. The lake and its adjacent thermokarst basin rapidly expanded during the early Holocene. This climatically warmer than today period was characterized by forest or forest tundra vegetation composed of larches, birch trees and shrubs. Woodlands of both the HTM and the Late Pleistocene were affected by fire, which potentially triggered the initiation of thermokarst processes resulting later in lake formation and expansion. The maximum lake depth at the study site and the lowest limnic bioproductivity occurred during the longest time interval of ~7 ka starting in the Holocene Thermal Maximum and lasting throughout the progressively cooler Neoglacial, whereas partial drainage and an extensive shift of the lake shoreline occurred ~0.9 cal. ka BP. Correspondingly, this study discusses different climatic and environmental drivers for the dynamics of a thermokarst basin.