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.

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.

Sedimentological, biogeomcial and geochronological data of thermokarst lake sediment core PG1967

Ice-rich permafrost landscapes are sensitive to climate and environmental change due to the melt-out of ground ice during thermokarst development. Thermokarst processes in the northern Yukon Territory are currently not well-documented. Lake sediments from Herschel Island (69°36'N; 139°04'W) in the western Canadian Arctic provide a record of thermokarst lake development since the early Holocene. A 727 cm long lake sediment core was analyzed for radiographic images, magnetic susceptibility, granulometry, and biogeochemical parameters (organic carbon, nitrogen, and stable carbon isotopes). Based on eight calibrated AMS radiocarbon dates, the sediment record covers the last ~ 11,500 years and was divided into four lithostratigraphic units (A to D) reflecting different thermokarst stages. Thermokarst initiation at the study area began ~ 11.5 cal ka BP. From ~ 11.5 to 10.0 cal ka BP, lake sediments of unit A started to accumulate in an initial lake basin created by melt-out of massive ground ice and thaw subsidence. Between 10.0 and 7.0 cal ka BP (unit B) the lake basin expanded in size and depth, attributed to talik formation during the Holocene thermal maximum. Higher-than-modern summer air temperatures led to increased lake productivity and widespread terrain disturbances in the lake's catchment. Thermokarst lake development between 7.0 and 1.8 cal ka BP (unit C) was characterized by a dynamic equilibrium, where lake basin and talik steadily expanded into ambient ice-rich terrain through shoreline erosion. Once lakes become deeper than the maximum winter lake ice thickness, thermokarst lake sediments show a great preservation potential. However, site-specific geomorphic factors such as episodic bank-shore erosion or sudden drainage through thermo-erosional valleys or coastal erosion breaching lake basins can disrupt continuous deposition. A hiatus in the record from 1.8 to 0.9 cal ka BP in Lake Herschel likely resulted from lake drainage or allochthonous slumping due to collapsing shore lines before continuous sedimentation of unit D recommenced during the last 900 years.

High-resolution digital elevation model (DEM) of a thermokarst depression in Siberian ice-rich permafrost deposits

This data set provides a high-resolution digital elevation model (DEM) of a thermokarst depression (~7 km²) on ice-complex deposits in the Arctic Lena Delta, Siberia. The DEM based on a geodetic field survey and was used for quantitative land surface analyses and detailed description of the thermokarst depression morphology. Detailed morphometrical analyses, volume calculations, and solar radiation modeling were performed and statistically analyzed by Ulrich et al. (2010) to investigate the asymmetrical thermokarst depression development and directed lake migration previously proposed by Morgenstern et al. (2008). Furthermore, the high-resolution DEM in combination with satellite data allowed detailed analyses of spatial and temporal landscape changes due to thermokarst development (Günther, 2009).

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.

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.

Age determinations from Bykovsky Peninsula, Siberia

A wide variety of environmental records is necessary for analysing and understanding the complex Late Quaternary dynamics of permafrost-dominated Arctic landscapes. A NE Siberian periglacial key region was studied in detail using sediment records, remote sensing data, and terrain modelling, all incorporated in a geographical information system (GIS). The study area consists of the Bykovsky Peninsula and the adjacent Khorogor Valley in the Kharaulakh Ridge situated a few kilometres southeast of the Lena Delta. In this study a comprehensive cryolithological database containing information from 176 sites was compiled. The information from these sites is based on the review of previously published borehole data, outcrop profiles, surface samples, and our own field data. These archives cover depositional records of three periods: from Pliocene to Early Pleistocene, the Late Pleistocene and the Holocene. The main sediment sequences on the Bykovsky Peninsula consist of up to 50 m thick ice-rich permafrost deposits (Ice Complex) that were accumulated during the Late Pleistocene. They were formed as a result of nival processes around extensive snowfields in the Kharaulakh Ridge, slope processes in these mountains (such as in the Khorogor Valley), and alluvial/proluvial sedimentation in a flat accumulation plain dominated by polygonal tundra in the mountain foreland (Bykovsky Peninsula). During the early to middle Holocene warming, a general landscape transformation occurred from an extensive Late Pleistocene accumulation plain to a strongly thermokarst-dominated relief dissected by numerous depressions. Thermokarst subsidence had an enormous influence on the periglacial hydrological patterns, the sediment deposition, and on the composition and distribution of habitats. Climate deterioration, lake drainage, and talik refreezing occurred during the middle to late Holocene. The investigated region was reached by the post-glacial sea level rise during the middle Holocene, triggering thermo-abrasion of ice-rich coasts and the marine inundation of thermokarst depressions.

Mid-Wisconsin to Holocene permafrost and landscape dynamics based on a drained lake basin core from the northern Seward Peninsula, Northwest Alaska

Permafrost-related processes drive regional landscape dynamics in the Arctic terrestrial system. A better understanding of past periods indicative of permafrost degradation and aggradation is important for predicting the future response of Arctic landscapes to climate change. Here, we used a multi-proxy approach to analyze a ~4 m long sediment core from a drained thermokarst lake basin on the northern Seward Peninsula in western Arctic Alaska (USA). Sedimentological, biogeochemistical, geochronological, micropaleontological (ostracoda, testate amoeba) and tephra analyses were used to determine the long-term environmental Early-Wisconsin to Holocene history preserved in our core for Central Beringia. Yedoma accumulation dominated throughout the Early to Late-Wisconsin but was interrupted by wetland formation from 44.5 to 41.5 ka BP. The latter was terminated by deposition of 1 m of volcanic tephra, most likely originating from the South Killeak Maar eruption at about 42 ka BP. Yedoma deposition continued until 22.5 ka BP and was followed by a depositional hiatus in the sediment core between 22.5 and 0.23 ka BP. We interpret this hiatus as due to intense thermokarst activity in the areas surrounding the site, which served as a sediment source during the Late-Wisconsin to Holocene climate transition. The lake forming the modern basin on the upland initiated around 0.23 ka BP, which drained catastrophically in spring 2005. The present study emphasizes that Arctic lake systems and periglacial landscapes are highly dynamic and permafrost formation as well as degradation in Central Beringia was controlled by regional to global climate patterns and as well as by local disturbances.

Pollen record with radiocarbon datings from Faddeyevskiy Island, Laptev Sea

Pollen, plant macrofossil, loss-on-ignition and radiocarbon analyses of a 1.4-m section in thermokarst topography from Faddeyevskiy Island (75°20'N, 143°50'E, 30 m elevation) provides new information on Late Pleistocene interstadial environmental history of this high Arctic region. Conventional radiocarbon dates (25,700 ± 1000, 32,780 ± 500, 35,200 ± 650 yr BP) and two AMS dates (29,950 ± 660 and 42,990 ± 1280 yr BP) indicate that the deposits accumulated during the Kargian (Boutellier) interval. Numerous mammoth (Mammuthus primigenius) remains that have been collected in vicinity of the site in this study were radio-carbon dated to 36,700-18,500 yr BP. Rare bison (Bison priscus) bones were dated to 32,200 ± 600 and 33,100 ± 320 yr BP. Poaceae, Cyperaceae, and Artemisia pollen dominate the spectra with some Ranunculaceae, Caryophyllaceae, Rosaceae, and Asteraceae. The pollen spectra reflect steppe-like (tundra-steppe) vegetation, which was dominant on the exposed shelf of the Arctic Ocean. Numerous Carex macrofossils suggest that the summer climate was at least 2°C warmer than today. The productivity of the local vegetation during the Kargian interstadial was high enough to feed the population of grazing mammals.