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.

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.

Robust trends of landscape dynamics in the Arctic Lena Delta with temporally dense Landsat time-series stacks, with links to GeoTIFFs

Arctic permafrost landscapes are among the most vulnerable and dynamic landscapes globally, but due to their extent and remoteness most of the landscape changes remain unnoticed. In order to detect disturbances in these areas we developed an automated processing chain for the calculation and analysis of robust trends of key land surface indicators based on the full record of available Landsat TM, ETM +, and OLI data. The methodology was applied to the ~ 29,000 km**2 Lena Delta in Northeast Siberia, where robust trend parameters (slope, confidence intervals of the slope, and intercept) were calculated for Tasseled Cap Greenness, Wetness and Brightness, NDVI, and NDWI, and NDMI based on 204 Landsat scenes for the observation period between 1999 and 2014. The resulting datasets revealed regional greening trends within the Lena Delta with several localized hot-spots of change, particularly in the vicinity of the main river channels. With a 30-m spatial resolution various permafrost-thaw related processes and disturbances, such as thermokarst lake expansion and drainage, fluvial erosion, and coastal changes were detected within the Lena Delta region, many of which have not been noticed or described before. Such hotspots of permafrost change exhibit significantly different trend parameters compared to non-disturbed areas. The processed dataset, which is made freely available through the data archive PANGAEA, will be a useful resource for further process specific analysis by researchers and land managers. With the high level of automation and the use of the freely available Landsat archive data, the workflow is scalable and transferrable to other regions, which should enable the comparison of land surface changes in different permafrost affected regions and help to understand and quantify permafrost landscape dynamics.