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.

Landform characterization using geophysics—Recent advances, applications, and emerging tools

This paper presents an overview of the strengths and limitations of existing and emerging geophysical tools for landform studies. The objectives are to discuss recent technical developments and to provide a review of relevant recent literature, with a focus on propagating field methods with terrestrial applications. For various methods in this category, including ground-penetrating radar (GPR), electrical resistivity (ER), seismics, and electromagnetic (EM) induction, the technical backgrounds are introduced, followed by section on novel developments relevant to landform characterization. For several decades, GPR has been popular for characterization of the shallow subsurface and in particular sedimentary systems. Novel developments in GPR include the use of multi-offset systems to improve signal-to-noise ratios and data collection efficiency, amongst others, and the increased use of 3D data. Multi-electrode ER systems have become popular in recent years as they allow for relatively fast and detailed mapping. Novel developments include time-lapse monitoring of dynamic processes as well as the use of capacitively-coupled systems for fast, non-invasive surveys. EM induction methods are especially popular for fast mapping of spatial variation, but can also be used to obtain information on the vertical variation in subsurface electrical conductivity. In recent years several examples of the use of plane wave EM for characterization of landforms have been published. Seismic methods for landform characterization include seismic reflection and refraction techniques and the use of surface waves. A recent development is the use of passive sensing approaches. The use of multiple geophysical methods, which can benefit from the sensitivity to different subsurface parameters, is becoming more common. Strategies for coupled and joint inversion of complementary datasets will, once more widely available, benefit the geophysical study of landforms.Three cases studies are presented on the use of electrical and GPR methods for characterization of landforms in the range of meters to 100. s of meters in dimension. In a study of polygonal patterned ground in the Saginaw Lowlands, Michigan, USA, electrical resistivity tomography was used to characterize differences in subsurface texture and water content associated with polygon-swale topography. Also, a sand-filled thermokarst feature was identified using electrical resistivity data. The second example is on the use of constant spread traversing (CST) for characterization of large-scale glaciotectonic deformation in the Ludington Ridge, Michigan. Multiple CST surveys parallel to an ~. 60. m high cliff, where broad (~. 100. m) synclines and narrow clay-rich anticlines are visible, illustrated that at least one of the narrow structures extended inland. A third case study discusses internal structures of an eolian dune on a coastal spit in New Zealand. Both 35 and 200. MHz GPR data, which clearly identified a paleosol and internal sedimentary structures of the dune, were used to improve understanding of the development of the dune, which may shed light on paleo-wind directions.

Characterization and Mapping of Patterned Ground in the Saginaw Lowlands, Michigan: Possible Evidence for Late-Wisconsin Permafrost

We identified, mapped, and characterized a widespread area (gt;1,020 km2) of patterned ground in the Saginaw Lowlands of Michigan, a wet, flat plain composed of waterlain tills, lacustrine deposits, or both. The polygonal patterned ground is interpreted as a possible relict permafrost feature, formed in the Late Wisconsin when this area was proximal to the Laurentide ice sheet. Cold-air drainage off the ice sheet might have pooled in the Saginaw Lowlands, which sloped toward the ice margin, possibly creating widespread but short-lived permafrost on this glacial lake plain. The majority of the polygons occur between the Glacial Lake Warren strandline (~14.8 cal. ka) and the shoreline of Glacial Lake Elkton (~14.3 cal. ka), providing a relative age bracket for the patterned ground. Most of the polygons formed in dense, wet, silt loam soils on flat-lying sites and take the form of reticulate nets with polygon long axes of 150 to 160 m and short axes of 60 to 90 m. Interpolygon swales, often shown as dark curvilinears on aerial photographs, are typically slightly lower than are the polygon centers they bound. Some portions of these interpolygon swales are infilled with gravel-free, sandy loam sediments. The subtle morphology and sedimentological characteristics of the patterned ground in the Saginaw Lowlands suggest that thermokarst erosion, rather than ice-wedge replacement, was the dominant geomorphic process associated with the degradation of the Late-Wisconsin permafrost in the study area and, therefore, was primarily responsible for the soil patterns seen there today.

Permafrost response to last interglacial warming: field evidence from non-glaciated Yukon and Alaska

We present stratigraphic observations from three sites in eastern Beringia - Ch'ijee's Bluff in northern Yukon and nearby exposures on the Old Crow River, the Palisades on the Yukon River in Alaska, and placer mining exposures at Thistle Creek in west-central Yukon - which provide insight into the response of permafrost to regional warming during the last interglaciation. Chronology is based on the presence of Old Crow tephra, an important regional stratigraphic marker that dates to late Marine Isotope Stage 6, supplemented by paleoecology and non-finite C ages on wood-rich organic silts. Old Crow tephra overlies several relict ice wedges at the Palisades and Thistle Creek, indicating that permafrost at these sites did not thaw completely during the last interglaciation. Prominent deposits of last interglacial wood-rich organic silt are present at multiple sites in eastern Beringia, and probably represent accumulations of reworked forest vegetation due to thaw slumping or deposition into thermokarst ponds or depressions. Consistent stratigraphic relations between these deposits, Old Crow tephra, and ice wedge pseudomorphs at our three study sites, and at least six other sites in eastern Beringia, suggest that thaw of shallow permafrost was widespread during the last interglaciation. Limited stratigraphic evidence suggests that thaw was probably on the order of meters, rather than 10s of meters. The ubiquity of shallow permafrost degradation during the last interglaciation suggests that current ground warming may foreshadow widespread near-surface thaw under even modest future warming scenarios. However, the persistence of relict pre-last interglacial ice wedges highlights the potential for the regional antiquity of discontinuous permafrost, and provides compelling field evidence for the long-term resilience of deep permafrost during sustained periods of warmer-than-present climate.