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.

Borehole temperature, submarine permafrost methane concentrations and pore water chemistry measured on borehole BK2, central Laptev Sea shelf

Submarine permafrost degradation has been invoked as a cause for recent observations of methane emissions from the seabed to the water column and atmosphere of the East Siberian shelf. Sediment drilled 52 m down from the sea ice in Buor Khaya Bay, central Laptev Sea revealed unfrozen sediment overlying ice-bonded permafrost. Methane concentrations in the overlying unfrozen sediment were low (mean 20 µM) but higher in the underlying ice-bonded submarine permafrost (mean 380 µM). In contrast, sulfate concentrations were substantially higher in the unfrozen sediment (mean 2.5 mM) than in the underlying submarine permafrost (mean 0.1 mM). Using deduced permafrost degradation rates, we calculate potential mean methane efflux from degrading permafrost of 120 mg/m**2 per year at this site. However, a drop of methane concentrations from 190 µM to 19 µM and a concomitant increase of methane d13C from -63 per mil to -35 per mil directly above the ice-bonded permafrost suggest that methane is effectively oxidized within the overlying unfrozen sediment before it reaches the water column. High rates of methane ebullition into the water column observed elsewhere are thus unlikely to have ice-bonded permafrost as their source.

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.

Dissolved organic carbon (DOC) in Arctic ground ice, from northwest Canada, east Siberia, and Alaska

Thermal permafrost degradation and coastal erosion in the Arctic remobilize substantial amounts of organic carbon (OC) and nutrients which have accumulated in late Pleistocene and Holocene unconsolidated deposits. Permafrost vulnerability to thaw subsidence, collapsing coastlines and irreversible landscape change are largely due to the presence of large amounts of massive ground ice such as ice wedges. However, ground ice has not, until now, been considered to be a source of dissolved organic carbon (DOC), dissolved inorganic carbon (DIC) and other elements which are important for ecosystems and carbon cycling. Here we show, using biogeochemical data from a large number of different ice bodies throughout the Arctic, that ice wedges have the greatest potential for DOC storage, with a maximum of 28.6 mg/L (mean: 9.6 mg/L). Variation in DOC concentration is positively correlated with and explained by the concentrations and relative amounts of typically terrestrial cations such as Mg2+ and K+. DOC sequestration into ground ice was more effective during the late Pleistocene than during the Holocene, which can be explained by rapid sediment and OC accumulation, the prevalence of more easily degradable vegetation and immediate incorporation into permafrost. We assume that pristine snowmelt is able to leach considerable amounts of well-preserved and highly bioavailable DOC as well as other elements from surface sediments, which are rapidly frozen and stored in ground ice, especially in ice wedges, even before further degradation. We found that ice wedges in the Yedoma region represent a significant DOC (45.2 Tg) and DIC (33.6 Tg) pool in permafrost areas and a freshwater reservoir of 4200 km**3. This study underlines the need to discriminate between particulate OC and DOC to assess the availability and vulnerability of the permafrost carbon pool for ecosystems and climate feedback upon mobilization.

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.

Radiocarbon ages, dose rates, water isotopic composition and hydrochemistry of samples from Komakuk Beach and Herschel Island

Terrestrial permafrost archives along the Yukon Coastal Plain (northwest Canada) have recorded landscape development and environmental change since the Late Wisconsinan at the interface of unglaciated Beringia (i.e. Komakuk Beach) and the northwestern limit of the Laurentide Ice Sheet (i.e. Herschel Island). The objective of this paper is to compare the late glacial and Holocene landscape development on both sides of the former ice margin based on permafrost sequences and ground ice. Analyses at these sites involved a multi-proxy approach including: sedimentology, cryostratigraphy, palaeoecology of ostracods, stable water isotopes in ground ice, hydrochemistry, and AMS radiocarbon and infrared stimulated luminescence (IRSL) dating. AMS and IRSL age determinations yielded full glacial ages at Komakuk Beach that is the northeastern limit of ice-free Beringia. Herschel Island to the east marks the Late Wisconsinan limit of the northwest Laurentide Ice Sheet and is composed of ice-thrust sediments containing plant detritus as young as 16.2 cal ka BP that might provide a maximum age on ice arrival. Late Wisconsinan ice wedges with sediment-rich fillings on Herschel Island are depleted in heavy oxygen isotopes (mean d18O of -29.1 per mil); this, together with low d-excess values, indicates colder-than-modern winter temperatures and probably reduced snow depths. Grain-size distribution and fossil ostracod assemblages indicate that deglaciation of the Herschel Island ice-thrust moraine was accompanied by alluvial, proluvial, and eolian sedimentation on the adjacent unglaciated Yukon Coastal Plain until ~11 cal ka BP during a period of low glacio-eustatic sea level. The late glacial-Holocene transition was marked by higher-than-modern summer temperatures leading to permafrost degradation that began no later than 11.2 cal ka BP and caused a regional thaw unconformity. Cryostructures and ice wedges were truncated while organic matter was incorporated and soluble ions were leached in the thaw zone. Thermokarst activity led to the formation of ice-wedge casts and deposition of thermokarst lake sediments. These were subsequently covered by rapidly accumulating peat during the early Holocene Thermal Maximum. A rising permafrost table, reduced peat accumulation, and extensive ice-wedge growth resulted from climate cooling starting in the middle Holocene until the late 20th century. The reconstruction of palaeolandscape dynamics on the Yukon Coastal Plain and the eastern Beringian edge contributes to unraveling the linkages between ice sheet, ocean, and permafrost that have existed since the Late Wisconsinan.

High-resolution dynamical downscaling products for the North Slope of Alaska and surrounding areas, links to model result files

Climatic changes are most pronounced in northern high latitude regions. Yet, there is a paucity of observational data, both spatially and temporally, such that regional-scale dynamics are not fully captured, limiting our ability to make reliable projections. In this study, a group of dynamical downscaling products were created for the period 1950 to 2100 to better understand climate change and its impacts on hydrology, permafrost, and ecosystems at a resolution suitable for northern Alaska. An ERA-interim reanalysis dataset and the Community Earth System Model (CESM) served as the forcing mechanisms in this dynamical downscaling framework, and the Weather Research & Forecast (WRF) model, embedded with an optimization for the Arctic (Polar WRF), served as the Regional Climate Model (RCM). This downscaled output consists of multiple climatic variables (precipitation, temperature, wind speed, dew point temperature, and surface air pressure) for a 10 km grid spacing at three-hour intervals. The modeling products were evaluated and calibrated using a bias-correction approach. The ERA-interim forced WRF (ERA-WRF) produced reasonable climatic variables as a result, yielding a more closely correlated temperature field than precipitation field when long-term monthly climatology was compared with its forcing and observational data. A linear scaling method then further corrected the bias, based on ERA-interim monthly climatology, and bias-corrected ERA-WRF fields were applied as a reference for calibration of both the historical and the projected CESM forced WRF (CESM-WRF) products. Biases, such as, a cold temperature bias during summer and a warm temperature bias during winter as well as a wet bias for annual precipitation that CESM holds over northern Alaska persisted in CESM-WRF runs. The linear scaling of CESM-WRF eventually produced high-resolution downscaling products for the Alaskan North Slope for hydrological and ecological research, together with the calibrated ERA-WRF run, and its capability extends far beyond that. Other climatic research has been proposed, including exploration of historical and projected climatic extreme events and their possible connections to low-frequency sea-atmospheric oscillations, as well as near-surface permafrost degradation and ice regime shifts of lakes. These dynamically downscaled, bias corrected climatic datasets provide improved spatial and temporal resolution data necessary for ongoing modeling efforts in northern Alaska focused on reconstructing and projecting hydrologic changes, ecosystem processes and responses, and permafrost thermal regimes. The dynamical downscaling methods presented in this study can also be used to create more suitable model input datasets for other sub-regions of the Arctic.

Ages and isotopic signatures of peat cores from a boreal landscape, Interior Alaska

To determine the influence of fire and thermokarst in a boreal landscape, we investigated peat cores within and adjacent to a permafrost collapse feature on the Tanana River Floodplain of Interior Alaska. Radioisotope dating, diatom assemblages, plant macrofossils, charcoal fragments, and carbon and nitrogen content of the peat profile indicate ~600 years of vegetation succession with a transition from a terrestrial forest to a sedge-dominated wetland over 100 years ago, and to a Sphagnum-dominated peatland in approximately 1970. The shift from sedge to Sphagnum, and a decrease in the detrended tree-ring width index of black spruce trees adjacent to the collapse coincided with an increase in the growing season temperature record from Fairbanks. This concurrent wetland succession and reduced growth of black spruce trees indicates a step-wise ecosystem-level response to a change in regional climate. In 2001, fire was observed coincident with permafrost collapse and resulted in lateral expansion of the peatland. These observations and the peat profile suggest that future warming and/or increased fire disturbance could promote permafrost degradation, peatland expansion, and increase carbon storage across this landscape; however, the development of drought conditions could reduce the success of both black spruce and Sphagnum, and potentially decrease the long-term ecosystem carbon storage.

Single non-normalized data of electron probe analyses of all glass shard samples from the Seward Peninsula and the Lipari obsidian reference standard

Permafrost degradation influences the morphology, biogeochemical cycling and hydrology of Arctic landscapes over a range of time scales. To reconstruct temporal patterns of early to late Holocene permafrost and thermokarst dynamics, site-specific palaeo-records are needed. Here we present a multi-proxy study of a 350-cm-long permafrost core from a drained lake basin on the northern Seward Peninsula, Alaska, revealing Lateglacial to Holocene thermokarst lake dynamics in a central location of Beringia. Use of radiocarbon dating, micropalaeontology (ostracods and testaceans), sedimentology (grain-size analyses, magnetic susceptibility, tephra analyses), geochemistry (total nitrogen and carbon, total organic carbon, d13Corg) and stable water isotopes (d18O, dD, d excess) of ground ice allowed the reconstruction of several distinct thermokarst lake phases. These include a pre-lacustrine environment at the base of the core characterized by the Devil Mountain Maar tephra (22 800±280 cal. a BP, Unit A), which has vertically subsided in places due to subsequent development of a deep thermokarst lake that initiated around 11 800 cal. a BP (Unit B). At about 9000 cal. a BP this lake transitioned from a stable depositional environment to a very dynamic lake system (Unit C) characterized by fluctuating lake levels, potentially intermediate wetland development, and expansion and erosion of shore deposits. Complete drainage of this lake occurred at 1060 cal. a BP, including post-drainage sediment freezing from the top down to 154 cm and gradual accumulation of terrestrial peat (Unit D), as well as uniform upward talik refreezing. This core-based reconstruction of multiple thermokarst lake generations since 11 800 cal. a BP improves our understanding of the temporal scales of thermokarst lake development from initiation to drainage, demonstrates complex landscape evolution in the ice-rich permafrost regions of Central Beringia during the Lateglacial and Holocene, and enhances our understanding of biogeochemical cycles in thermokarst-affected regions of the Arctic.

Ecological classification of Herschel Island based on supervised classification of multispectral satellite imagery

Permafrost landscapes experience different disturbances and store large amounts of organic matter, which may become a source of greenhouse gases upon permafrost degradation. We analysed the influence of terrain and geomorphic disturbances (e.g. soil creep, active-layer detachment, gullying, thaw slumping, accumulation of fluvial deposits) on soil organic carbon (SOC) and total nitrogen (TN) storage using 11 permafrost cores from Herschel Island, western Canadian Arctic. Our results indicate a strong correlation between SOC storage and the topographic wetness index. Undisturbed sites stored the majority of SOC and TN in the upper 70 cm of soil. Sites characterised by mass wasting showed significant SOC depletion and soil compaction, whereas sites characterised by the accumulation of peat and fluvial deposits store SOC and TN along the whole core. We upscaled SOC and TN to estimate total stocks using the ecological units determined from vegetation composition, slope angle and the geomorphic disturbance regime. The ecological units were delineated with a supervised classification based on RapidEye multispectral satellite imagery and slope angle. Mean SOC and TN storage for the uppermost 1?m of soil on Herschel Island are 34.8 kg C/m**2 and 3.4 kg N/m**2, respectively.

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).

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.