Database of Ice-Rich Yedoma Permafrost (IRYP)

Vast portions of Arctic and sub-Arctic Siberia, Alaska and the Yukon Territory are covered by ice-rich silty to sandy deposits that are containing large ice wedges, resulting from syngenetic sedimentation and freezing. Accompanied by wedge-ice growth in polygonal landscapes, the sedimentation process was driven by cold continental climatic and environmental conditions in unglaciated regions during the late Pleistocene, inducing the accumulation of the unique Yedoma deposits up to >50 meters thick. Because of fast incorporation of organic material into syngenetic permafrost during its formation, Yedoma deposits include well-preserved organic matter. Ice-rich deposits like Yedoma are especially prone to degradation triggered by climate changes or human activity. When Yedoma deposits degrade, large amounts of sequestered organic carbon as well as other nutrients are released and become part of active biogeochemical cycling. This could be of global significance for future climate warming as increased permafrost thaw is likely to lead to a positive feedback through enhanced greenhouse gas fluxes. Therefore, a detailed assessment of the current Yedoma deposit coverage and its volume is of importance to estimate its potential response to future climate changes. We synthesized the map of the coverage and thickness estimation, which will provide critical data needed for further research. In particular, this preliminary Yedoma map is a great step forward to understand the spatial heterogeneity of Yedoma deposits and its regional coverage. There will be further applications in the context of reconstructing paleo-environmental dynamics and past ecosystems like the mammoth-steppe-tundra, or ground ice distribution including future thermokarst vulnerability. Moreover, the map will be a crucial improvement of the data basis needed to refine the present-day Yedoma permafrost organic carbon inventory, which is assumed to be between 83±12 (Strauss et al., 2013, doi:10.1002/2013GL058088) and 129±30 (Walter Anthony et al., 2014, doi:10.1038/nature13560) gigatonnes (Gt) of organic carbon in perennially-frozen archives. Hence, here we synthesize data on the circum-Arctic and sub-Arctic distribution and thickness of Yedoma for compiling a preliminary circum-polar Yedoma map. For compiling this map, we used (1) maps of the previous Yedoma coverage estimates, (2) included the digitized areas from Grosse et al. (2013) as well as extracted areas of potential Yedoma distribution from additional surface geological and Quaternary geological maps (1.: 1:500,000: Q-51-V,G; P-51-A,B; P-52-A,B; Q-52-V,G; P-52-V,G; Q-51-A,B; R-51-V,G; R-52-V,G; R-52-A,B; 2.: 1:1,000,000: P-50-51; P-52-53; P-58-59; Q-42-43; Q-44-45; Q-50-51; Q-52-53; Q-54-55; Q-56-57; Q-58-59; Q-60-1; R-(40)-42; R-43-(45); R-(45)-47; R-48-(50); R-51; R-53-(55); R-(55)-57; R-58-(60); S-44-46; S-47-49; S-50-52; S-53-55; 3.: 1:2,500,000: Quaternary map of the territory of Russian Federation, 4.: Alaska Permafrost Map). The digitalization was done using GIS techniques (ArcGIS) and vectorization of raster Images (Adobe Photoshop and Illustrator). Data on Yedoma thickness are obtained from boreholes and exposures reported in the scientific literature. The map and database are still preliminary and will have to undergo a technical and scientific vetting and review process. In their current form, we included a range of attributes for Yedoma area polygons based on lithological and stratigraphical information from the original source maps as well as a confidence level for our classification of an area as Yedoma (3 stages: confirmed, likely, or uncertain). In its current version, our database includes more than 365 boreholes and exposures and more than 2000 digitized Yedoma areas. We expect that the database will continue to grow. In this preliminary stage, we estimate the Northern Hemisphere Yedoma deposit area to cover approximately 625,000 km². We estimate that 53% of the total Yedoma area today is located in the tundra zone, 47% in the taiga zone. Separated from west to east, 29% of the Yedoma area is found in North America and 71 % in North Asia. The latter include 9% in West Siberia, 11% in Central Siberia, 44% in East Siberia and 7% in Far East Russia. Adding the recent maximum Yedoma region (including all Yedoma uplands, thermokarst lakes and basins, and river valleys) of 1.4 million km² (Strauss et al., 2013, doi:10.1002/2013GL058088) and postulating that Yedoma occupied up to 80% of the adjacent formerly exposed and now flooded Beringia shelves (1.9 million km², down to 125 m below modern sea level, between 105°E - 128°W and >68°N), we assume that the Last Glacial Maximum Yedoma region likely covered more than 3 million km² of Beringia. Acknowledgements: This project is part of the Action Group "The Yedoma Region: A Synthesis of Circum-Arctic Distribution and Thickness" (funded by the International Permafrost Association (IPA) to J. Strauss) and is embedded into the Permafrost Carbon Network (working group Yedoma Carbon Stocks). We acknowledge the support by the European Research Council (Starting Grant #338335), the German Federal Ministry of Education and Research (Grant 01DM12011 and "CarboPerm" (03G0836A)), the Initiative and Networking Fund of the Helmholtz Association (#ERC-0013) and the German Federal Environment Agency (UBA, project UFOPLAN FKZ 3712 41 106).

(Table 1) Radiocarbon and OSL ages of sediment cores C1-C5 drilled during the COAST expedition to Cape Mamontov Klyk in April 2005

The palaeoenvironmental development of the western Laptev Sea is understood primarily from investigations of exposed cliffs and surface sediment cores from the shelf. In 2005, a core transect was drilled between the Taymyr Peninsula and the Lena Delta, an area that was part of the westernmost region of the non-glaciated Beringian landmass during the late Quaternary. The transect of five cores, one terrestrial and four marine, taken near Cape Mamontov Klyk reached 12 km offshore and 77 m below sea level. A multiproxy approach combined cryolithological, sedimentological, geochronological (14C-AMS, OSL on quartz, IR-OSL on feldspars) and palaeoecological (pollen, diatoms) methods. Our interpretation of the proxies focuses on landscape history and the transition of terrestrial into subsea permafrost. Marine interglacial deposits overlain by relict terrestrial permafrost within the same offshore core were encountered in the western Laptev Sea. Moreover, the marine interglacial deposits lay unexpectedly deep at 64 m below modern sea level 12 km from the current coastline, while no marine deposits were encountered onshore. This implies that the position of the Eemian coastline presumably was similar to today's. The landscape reconstruction suggests Eemian coastal lagoons and thermokarst lakes, followed by Early to Middle Weichselian fluvially dominated terrestrial deposition. During the Late Weichselian, this fluvial landscape was transformed into a poorly drained accumulation plain, characterized by widespread and broad ice-wedge polygons. Finally, the shelf plain was flooded by the sea during the Holocene, resulting in the inundation and degradation of terrestrial permafrost and its transformation into subsea permafrost.

Inventory of thermo-erosional valleys and streams in three ice-rich permafrost lowlands adjacent to the Laptev Sea, with links to maps in ArcGIS format

This dataset provides an inventory of thermo-erosional landforms and streams in three lowland areas underlain by ice-rich permafrost of the Yedoma-type Ice Complex at the Siberian Laptev Sea coast. It consists of two shapefiles per study region: one shapefile for the digitized thermo-erosional landforms and streams, one for the study area extent. Thermo-erosional landforms were manually digitized from topographic maps and satellite data as line features and subsequently analyzed in a Geographic Information System (GIS) using ArcGIS 10.0. The mapping included in particular thermo-erosional gullies and valleys as well as streams and rivers, since development of all of these features potentially involved thermo-erosional processes. For the Cape Mamontov Klyk site, data from Grosse et al. [2006], which had been digitized from 1:100000 topographic map sheets, were clipped to the Ice Complex extent of Cape Mamontov Klyk, which excludes the hill range in the southwest with outcropping bedrock and rocky slope debris, coastal barrens, and a large sandy floodplain area in the southeast. The mapped features (streams, intermittent streams) were then visually compared with panchromatic Landsat-7 ETM+ satellite data (4 August 2000, 15 m spatial resolution) and panchromatic Hexagon data (14 July 1975, 10 m spatial resolution). Smaller valleys and gullies not captured in the maps were subsequently digitized from the satellite data. The criterion for the mapping of linear features as thermo-erosional valleys and gullies was their clear incision into the surface with visible slopes. Thermo-erosional features of the Lena Delta site were mapped on the basis of a Landsat-7 ETM+ image mosaic (2000 and 2001, 30 m ground resolution) [Schneider et al., 2009] and a Hexagon satellite image mosaic (1975, 10 m ground resolution) [G. Grosse, unpublished data] of the Lena River Delta within the extent of the Lena Delta Ice Complex [Morgenstern et al., 2011]. For the Buor Khaya Peninsula, data from Arcos [2012], which had been digitized based on RapidEye satellite data (8 August 2010, 6.5 m ground resolution), were completed for smaller thermo-erosional features using the same RapidEye scene as a mapping basis. The spatial resolution, acquisition date, time of the day, and viewing geometry of the satellite data used may have influenced the identification of thermo-erosional landforms in the images. For Cape Mamontov Klyk and the Lena Delta, thermo-erosional features were digitized using both Hexagon and Landsat data; Hexagon provided higher resolution and Landsat provided the modern extent of features. Allowance of up to decameters was made for the lateral expansion of features between Hexagon and Landsat acquisitions (between 1975 and 2000).

Bulk sediment parameter of three sediment cores from permafrost lowlands, Alaska Arctic Coastal Plain

Arctic lowland landscapes have been modified by thermokarst lake processes throughout the Holocene. Thermokarst lakes form as a result of ice-rich permafrost degradation and they may expand over time through thermal and mechanical shoreline erosion. We studied proximal and distal sedimentary records from a thermokarst lake located on the Arctic Coastal Plain of northern Alaska to reconstruct the impact of catchment dynamics and morphology on the lacustrine depositional environment and to quantify carbon accumulation in thermokarst lake sediments. Short cores were collected for analysis of pollen, sedimentological and geochemical proxies. Radiocarbon and Pb/Cs dating, as well as extrapolation of measured historic lake expansion rates, were applied to estimate a minimum lake age of ~ 1,400 calendar years BP. The pollen record is in agreement with the young lake age as it does not include evidence of the "alder high" that occurred in the region ~ 4.0 cal ka BP. The lake most likely initiated from a remnant pond in a drained thermokarst lake basin (DTLB) and deepened rapidly as evidenced by accumulation of laminated sediments. Increasing oxygenation of the water column as shown by higher Fe/Ti and Fe/S ratios in the sediment indicate shifts in ice regime with increasing water depth. More recently, the sediment source changed as the thermokarst lake expanded through lateral permafrost degradation, alternating from redeposited DTLB sediments, to increased amounts of sediment from eroding, older upland deposits, followed by a more balanced combination of both DTLB and upland sources. The characterizing shifts in sediment sources and depositional regimes in expanding thermokarst lakes were therefore archived in the thermokarst lake sedimentary record. This study also highlights the potential for Arctic lakes to recycle old carbon from thawing permafrost and thermokarst processes.

Characteristics of samples obtained during Pokhodsk 2012-2013 campaigns in the joint Russian-German POLYGON Project

Polygonal tundra, thermokarst basins and pingos are common and characteristic periglacial features of arctic lowlands underlain by permafrost in Northeast Siberia. Modern polygonal mires are in the focus of biogeochemical, biological, pedological, and cryolithological research with special attention to their carbon stocks and greenhouse-gas fluxes, their biodiversity and their dynamics and functioning under past, present and future climate scenarios. Within the frame of the joint German-Russian DFG-RFBR project Polygons in tundra wetlands: state and dynamics under climate variability in Polar Regions (POLYGON) field studies of recent and of late Quaternary environmental dynamics were carried out in the Indigirka lowland and in the Kolyma River Delta in summer 2012 and summer 2013. Using a multidisciplinary approach, several types of polygons and thermokarst lakes were studied in different landscapes units in the Kolyma Delta in 2012 around the small fishing settlement Pokhodsk. The floral and faunal associations of polygonal tundra were described during the fieldwork. Ecological, hydrological, meteorological, limnological, pedological and cryological features were studied in order to evaluate modern and past environmental conditions and their essential controlling parameters. The ecological monitoring and collection program of polygonal ponds were undertaken as in 2011 in the Indigirka lowland by a former POLYGON expedition (Schirrmeister et al. [eds.] 2012). Exposures, pits and drill cores in the Kolyma Delta were studied to understand the cryolithological structures of frozen ground and to collect samples for detailed paleoenvironmental research of the late Quaternary past. Dendrochronological and ecological studies were carried out in the tree line zone south of the Kolyma Delta. Based on previous work in the Indigirka lowland in 2011 (Schirrmeister et al. [eds.] 2012), the environmental monitoring around the Kytalyk research station was continued until the end of August 2012. In addition, a classical exposure of the late Pleistocene permafrost at the Achchaygy Allaikha River near Chokurdakh was studied. The ecological studies near Pokhodsk were continued in 2013 (chapter 13). Other fieldwork took place at the Pokhodsk-Yedoma-Island in the northwestern part of the Kolyma Delta.

Physical, chemical and biological characteristics of disturbed and undisturbed Arctic lakes, Northwest Territories

1. Global warming is predicted to cause changes in permafrost cover and stability in the Arctic. Zones of high ion concentration in regions of ice-rich permafrost are a reservoir of chemicals that can potentially be transferred to fresh waters during thawing. Consequently, input of enriched runoff from the thaw and sediment and vegetation from the landscape could alter lakes by affecting their geochemistry and biological production. 2. Three undisturbed lakes and five lakes disturbed by retrogressive permafrost thaw slumps were sampled during late summer of 2006 to assess the potential effects of thermokarst shoreline slumping on water and sediment chemistry, the underwater light regime, and benthic macrophyte biomass and community structure. 3. Undisturbed lakes had sediments rich in organic material and selected micronutrients, while disturbed lakes had sediments richer in calcium, magnesium and strontium, greater transparency of the water column, and a well-developed submerged macrophyte community. 4. It is postulated that enriched runoff chemistry may alter nutrient availability at the sediment-water interface and also the degradation of organic material, thus affecting lake transparency and submerged macrophytes. The results suggest that retrogressive permafrost slumping can significantly affect food webs in arctic tundra lakes through an increase in macrophyte biomass and development of a more complex benthic habitat.

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.

Description and chemical and mineral composition of iron oxide samples from two Finnish lakes

Fifteen iron oxide accumulations from the bottoms of two Finnish lakes ("lake ores") were found to contain as much as 50% Fe. Differential X-ray powder diffraction and selective dissolution by oxalate showed that the samples consisted of poorly crystallized goethite and ferrihydrite. The crust ores of one lake had higher ferrihydrite to goethite ratios than the nodular ores of the other lake. The higher ferrihydrite proportion was attributed to a higher rate of Fe2+ supply from the ground water and/or a higher rate of oxidation as a function of water depth and bottom-sediment permeability. Values of Al-for-Fe substitution of the goethites determined from unit-cell dimensions agreed with those obtained from chemical extraction if the unit-cell volume rather than the c dimension was used. In very small goethite crystals a slight expansion of the a unit-cell dimension is probaby compensated by a corresponding contraction of the c dimension, so that a contraction of the c dimension need not necessarily be caused by Al substitution. The goethites of the two lakes differed significantly in their Al-for-Fe substitutions and hence in their unit-cell sizes, OH-bending characteristics, dehydroxylation temperatures, dissolution kinetics, and Mössbauer parameters. The difference in Al substitution (0 vs. 7 mole %) is attributed to the Al-supplying power of the bottom sediments: the silty-clayey sediments in one lake appear to have supplied A1 during goethite formation, whereas the gravelly-sandy sediments in the other lake did not. The compositions of the goethites thus reflect their environments of formation.

Oxygen and carbon isotopic composition of shells from Radix sp. in lakes across the Tibetan Plateau

The Tibetan Plateau (TP), including its surrounding mountain ranges, represents the largest store of ice outside the polar regions. It hosts numerous lakes as well as the head waters of major Asian rivers, on which billions of people depend, and it is particularly sensitive to climate change. The moisture transport to the TP is controlled by the Indian and Pacific monsoon and the Westerlies. Understanding the evolution of the interaction of these circulation systems requires studies on climate archives in different spatial and temporal contexts. The objective of this study is to learn more about the interannual variability of precipitation patterns across the TP and how different hydrologic systems react to different climatic factors. Aragonite shells of the aquatic gastropod Radix, which is widely distributed in the region, may represent suitable archives for inferring hydrologic and climatic signals in particularly high resolution. Therefore, sclerochronological studies of d18O and d13C ratios in Radix shells from seven lakes were conducted, each representing a different hydrologic and climatic setting, on a transect from the Pamirs across the TP. The shell patterns exhibit an increasing influence of precipitation and a decreasing influence of evaporation on the isotope compositions from west to east. d18O values of shells from lakes on the eastern and central TP (Donggi Cona, Yamdrok Yumco, Tarab Co) mirror monsoon signals, indicated by more negative values and higher variabilities compared to the more western lakes (Karakul, Bangong/Nyak, Manasarovar). In Yadang Co, located on the central southern TP, the monsoon rains did not reach the lake in the sampling year, although it is located in a region which is usually affected by monsoon circulation. The d18O values are used to differentiate the annual hydrological cycle into ice cover period, melt water period, precipitation period and evaporation period. d13C compositions in the shells particularly depend on specific habitats, which vary in biological productivity and in carbon sources. d18O and d13C patterns show a positive covariance in shells originating from large closed basins. The results show that Radix shells mirror general climatic differences between the seven lake regions. These differences reflect both regional and local climate signals in sub-seasonal resolution, without noticeable dependence on the particular lake system.