Refine
Document Type
- Article (32)
- Postprint (10)
- Doctoral Thesis (1)
- Other (1)
Language
- English (44)
Is part of the Bibliography
- yes (44)
Keywords
- permafrost (10)
- Siberia (4)
- Thermokarst (4)
- Permafrost (3)
- permafrost thaw (3)
- thermokarst (3)
- time series (3)
- Alaska (2)
- Arctic (2)
- Arctic lakes (2)
Sub-seasonal thaw slump mass wasting is not consistently energy limited at the landscape scale
(2018)
Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of timescales. On sub-seasonal timescales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellites. The estimated elevation changes at 12m resolution complement the commonly observed planimetric retreat rates by providing information on volume losses. Their high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km(2)) allow us to track mass wasting as drivers such as the available energy change during the summer of 2015 in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm day 1) despite the ample available energy, suggesting the widespread presence of early season insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm day 1), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August respectively. The observed temporal patterns are poorly related to slump characteristics like the headwall height. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.
Sub-seasonal thaw slump mass wasting is not consistently energy limited at the landscape scale
(2018)
Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of timescales. On sub-seasonal timescales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellites. The estimated elevation changes at 12m resolution complement the commonly observed planimetric retreat rates by providing information on volume losses. Their high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km(2)) allow us to track mass wasting as drivers such as the available energy change during the summer of 2015 in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm day 1) despite the ample available energy, suggesting the widespread presence of early season insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm day 1), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August respectively. The observed temporal patterns are poorly related to slump characteristics like the headwall height. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.
Glacial-interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (> 40 degrees N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.
Tundra be dammed
(2018)
Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.
Permafrost is a distinct feature of the terrestrial Arctic and is vulnerable to climate warming. Permafrost degrades in different ways, including deepening of a seasonally unfrozen surface and localized but rapid development of deep thaw features. Pleistocene ice-rich permafrost with syngenetic ice-wedges, termed Yedoma deposits, are widespread in Siberia, Alaska, and Yukon, Canada and may be especially prone to rapid-thaw processes. Freeze-locked organic matter in such deposits can be re-mobilized on short time-scales and contribute to a carbon-cycle climate feedback. Here we synthesize the characteristics and vulnerability of Yedoma deposits by synthesizing studies on the Yedoma origin and the associated organic carbon pool. We suggest that Yedoma deposits accumulated under periglacial weathering, transport, and deposition dynamics in non-glaciated regions during the late Pleistocene until the beginning of late glacial warming. The deposits formed due to a combination of aeolian, colluvial, nival, and alluvial deposition and simultaneous ground ice accumulation. We found up to 130 gigatons organic carbon in Yedoma, parts of which are well-preserved and available for fast decomposition after thaw. Based on incubation experiments, up to 10% of the Yedoma carbon is considered especially decomposable and may be released upon thaw. The substantial amount of ground ice in Yedoma makes it highly vulnerable to disturbances such as thermokarst and thermo-erosion processes. Mobilization of permafrost carbon is expected to increase under future climate warming. Our synthesis results underline the need of accounting for Yedoma carbon stocks in next generation Earth-System-Models for a more complete representation of the permafrost-carbon feedback.
Permafrost is warming globally which leads to widespread permafrost thaw. Particularly ice-rich permafrost is vulnerable to rapid thaw and erosion, impacting whole landscapes and ecosystems. Retrogressive thaw slumps (RTS) are abrupt permafrost disturbances that expand by several meters each year and lead to an increased soil organic carbon release. Local Remote Sensing studies identified increasing RTS activity in the last two decades by increasing number of RTS or heightened RTS growth rates. However, a large-scale assessment across diverse permafrost regions and at high temporal resolution allowing to further determine RTS thaw dynamics and its main drivers is still lacking. In this study we apply the disturbance detection algorithm LandTrendr for automated large-scale RTS mapping and high temporal thaw dynamic assessment to North Siberia (8.1 x 106km2). We adapted and parametrised the temporal segmentation algorithm for abrupt disturbance detection to incorporate Landsat+Sentinel-2 mosaics, conducted spectral filtering, spatial masking and filtering, and a binary machine-learning object classification of the disturbance output to separate between RTS and false positives (F1 score: 0.609). Ground truth data for calibration and validation of the workflow was collected from 9 known RTS cluster sites using very highresolution RapidEye and PlanetScope imagery. Our study presents the first automated detection and assessment of RTS and their temporal dynamics at largescale for 2001-2019. We identified 50,895 RTS and a steady increase in RTS-affected area from 2001 to 2019 across North Siberia, with a more abrupt increase from 2016 onward. Overall the RTS-affected area increased by 331% compared to 2000 (2000: 20,158 ha, 2001-2019: 66,699 ha). Contrary to this, 5 focus sites show spatiotemporal variability in their annual RTS dynamics, with alternating periods of increased and decreased RTS development, indicating a close relationship to thaw drivers. The majority of identified RTS was active from 2000 onward and only a small proportion initiated during the assessment period, indicating that the increase in RTS-affected area was mainly caused by enlarging existing RTS and not by new RTS. The detected increase in RTS dynamics suggests advancing permafrost thaw and underlines the importance of assessing abrupt permafrost disturbances with high spatial and temporal resolution at large-scales. Obtaining such consistent disturbance products will help to parametrise regional and global climate change models.
Permafrost is warming in the northern high latitudes, inducing highly dynamic thaw-related permafrost disturbances across the terrestrial Arctic. Monitoring and tracking of permafrost disturbances is important as they impact surrounding landscapes, ecosystems and infrastructure. Remote sensing provides the means to detect, map, and quantify these changes homogeneously across large regions and time scales. Existing Landsat-based algorithms assess different types of disturbances with similar spatiotemporal requirements. However, Landsat-based analyses are restricted in northern high latitudes due to the long repeat interval and frequent clouds, in particular at Arctic coastal sites. We therefore propose to combine Landsat and Sentinel-2 data for enhanced data coverage and present a combined annual mosaic workflow, expanding currently available algorithms, such as LandTrendr, to achieve more reliable time series analysis. We exemplary test the workflow for twelve sites across the northern high latitudes in Siberia. We assessed the number of images and cloud-free pixels, the spatial mosaic coverage and the mosaic quality with spectral comparisons. The number of available images increased steadily from 1999 to 2019 but especially from 2016 onward with the addition of Sentinel-2 images. Consequently, we have an increased number of cloud-free pixels even under challenging environmental conditions, which then serve as the input to the mosaicking process. In a comparison of annual mosaics, the Landsat+Sentinel-2 mosaics always fully covered the study areas (99.9–100 %), while Landsat-only mosaics contained data-gaps in the same years, only reaching coverage percentages of 27.2 %, 58.1 %, and 69.7 % for Sobo Sise, East Taymyr, and Kurungnakh in 2017, respectively. The spectral comparison of Landsat image, Sentinel-2 image, and Landsat+Sentinel-2 mosaic showed high correlation between the input images and mosaic bands (e.g., for Kurungnakh 0.91–0.97 between Landsat and Landsat+Sentinel-2 mosaic and 0.92–0.98 between Sentinel-2 and Landsat+Sentinel-2 mosaic) across all twelve study sites, testifying good quality mosaic results. Our results show that especially the results for northern, coastal areas was substantially improved with the Landsat+Sentinel-2 mosaics. By combining Landsat and Sentinel-2 data we accomplished to create reliably high spatial resolution input mosaics for time series analyses. Our approach allows to apply a high temporal continuous time series analysis to northern high latitude permafrost regions for the first time, overcoming substantial data gaps, and assess permafrost disturbance dynamics on an annual scale across large regions with algorithms such as LandTrendr by deriving the location, timing and progression of permafrost thaw disturbances
The Arctic-Boreal regions experience strong changes of air temperature and precipitation regimes, which affect the thermal state of the permafrost. This results in widespread permafrost-thaw disturbances, some unfolding slowly and over long periods, others occurring rapidly and abruptly. Despite optical remote sensing offering a variety of techniques to assess and monitor landscape changes, a persistent cloud cover decreases the amount of usable images considerably. However, combining data from multiple platforms promises to increase the number of images drastically. We therefore assess the comparability of Landsat-8 and Sentinel-2 imagery and the possibility to use both Landsat and Sentinel-2 images together in time series analyses, achieving a temporally-dense data coverage in Arctic-Boreal regions. We determined overlapping same-day acquisitions of Landsat-8 and Sentinel-2 images for three representative study sites in Eastern Siberia. We then compared the Landsat-8 and Sentinel-2 pixel-pairs, downscaled to 60 m, of corresponding bands and derived the ordinary least squares regression for every band combination. The acquired coefficients were used for spectral bandpass adjustment between the two sensors. The spectral band comparisons showed an overall good fit between Landsat-8 and Sentinel-2 images already. The ordinary least squares regression analyses underline the generally good spectral fit with intercept values between 0.0031 and 0.056 and slope values between 0.531 and 0.877. A spectral comparison after spectral bandpass adjustment of Sentinel-2 values to Landsat-8 shows a nearly perfect alignment between the same-day images. The spectral band adjustment succeeds in adjusting Sentinel-2 spectral values to Landsat-8 very well in Eastern Siberian Arctic-Boreal landscapes. After spectral adjustment, Landsat and Sentinel-2 data can be used to create temporally-dense time series and be applied to assess permafrost landscape changes in Eastern Siberia. Remaining differences between the sensors can be attributed to several factors including heterogeneous terrain, poor cloud and cloud shadow masking, and mixed pixels.
The Arctic-Boreal regions experience strong changes of air temperature and precipitation regimes, which affect the thermal state of the permafrost. This results in widespread permafrost-thaw disturbances, some unfolding slowly and over long periods, others occurring rapidly and abruptly. Despite optical remote sensing offering a variety of techniques to assess and monitor landscape changes, a persistent cloud cover decreases the amount of usable images considerably. However, combining data from multiple platforms promises to increase the number of images drastically. We therefore assess the comparability of Landsat-8 and Sentinel-2 imagery and the possibility to use both Landsat and Sentinel-2 images together in time series analyses, achieving a temporally-dense data coverage in Arctic-Boreal regions. We determined overlapping same-day acquisitions of Landsat-8 and Sentinel-2 images for three representative study sites in Eastern Siberia. We then compared the Landsat-8 and Sentinel-2 pixel-pairs, downscaled to 60 m, of corresponding bands and derived the ordinary least squares regression for every band combination. The acquired coefficients were used for spectral bandpass adjustment between the two sensors. The spectral band comparisons showed an overall good fit between Landsat-8 and Sentinel-2 images already. The ordinary least squares regression analyses underline the generally good spectral fit with intercept values between 0.0031 and 0.056 and slope values between 0.531 and 0.877. A spectral comparison after spectral bandpass adjustment of Sentinel-2 values to Landsat-8 shows a nearly perfect alignment between the same-day images. The spectral band adjustment succeeds in adjusting Sentinel-2 spectral values to Landsat-8 very well in Eastern Siberian Arctic-Boreal landscapes. After spectral adjustment, Landsat and Sentinel-2 data can be used to create temporally-dense time series and be applied to assess permafrost landscape changes in Eastern Siberia. Remaining differences between the sensors can be attributed to several factors including heterogeneous terrain, poor cloud and cloud shadow masking, and mixed pixels.
Permafrost is warming in the northern high latitudes, inducing highly dynamic thaw-related permafrost disturbances across the terrestrial Arctic. Monitoring and tracking of permafrost disturbances is important as they impact surrounding landscapes, ecosystems and infrastructure. Remote sensing provides the means to detect, map, and quantify these changes homogeneously across large regions and time scales. Existing Landsat-based algorithms assess different types of disturbances with similar spatiotemporal requirements. However, Landsat-based analyses are restricted in northern high latitudes due to the long repeat interval and frequent clouds, in particular at Arctic coastal sites. We therefore propose to combine Landsat and Sentinel-2 data for enhanced data coverage and present a combined annual mosaic workflow, expanding currently available algorithms, such as LandTrendr, to achieve more reliable time series analysis. We exemplary test the workflow for twelve sites across the northern high latitudes in Siberia. We assessed the number of images and cloud-free pixels, the spatial mosaic coverage and the mosaic quality with spectral comparisons. The number of available images increased steadily from 1999 to 2019 but especially from 2016 onward with the addition of Sentinel-2 images. Consequently, we have an increased number of cloud-free pixels even under challenging environmental conditions, which then serve as the input to the mosaicking process. In a comparison of annual mosaics, the Landsat+Sentinel-2 mosaics always fully covered the study areas (99.9–100 %), while Landsat-only mosaics contained data-gaps in the same years, only reaching coverage percentages of 27.2 %, 58.1 %, and 69.7 % for Sobo Sise, East Taymyr, and Kurungnakh in 2017, respectively. The spectral comparison of Landsat image, Sentinel-2 image, and Landsat+Sentinel-2 mosaic showed high correlation between the input images and mosaic bands (e.g., for Kurungnakh 0.91–0.97 between Landsat and Landsat+Sentinel-2 mosaic and 0.92–0.98 between Sentinel-2 and Landsat+Sentinel-2 mosaic) across all twelve study sites, testifying good quality mosaic results. Our results show that especially the results for northern, coastal areas was substantially improved with the Landsat+Sentinel-2 mosaics. By combining Landsat and Sentinel-2 data we accomplished to create reliably high spatial resolution input mosaics for time series analyses. Our approach allows to apply a high temporal continuous time series analysis to northern high latitude permafrost regions for the first time, overcoming substantial data gaps, and assess permafrost disturbance dynamics on an annual scale across large regions with algorithms such as LandTrendr by deriving the location, timing and progression of permafrost thaw disturbances
Vast areas of the terrestrial Subarctic and Arctic are underlain by permafrost. Landscape evolution is therefore largely controlled by climate-driven periglacial processes. The response of the frozen ground to late Quaternary warm and cold stages is preserved in permafrost sequences, and deducible by multi-proxy palaeoenvironmental approaches. Here, we analyse radiocarbon-dated mid-Wisconsin Interstadial and Holocene lacustrine deposits preserved in the Kit-1 pingo permafrost sequence combined with water and surface sediment samples from nine modern water bodies on Seward Peninsula (NW Alaska) to reconstruct thermokarst dynamics and determine major abiotic factors that controlled the aquatic ecosystem variability. Our methods comprise taxonomical diatom analyses as well as Detrended Correspondence Analysis (DCA) and Redundancy Analysis (RDA). Our results show, that the fossil diatom record reflects thermokarst lake succession since about 42 C-14 kyr BP. Different thermolcarst lake stages during the mid-Wisconsin Interstadial, the late Wisconsin and the early Holocene are mirrored by changes in diatom abundance, diversity, and ecology. We interpret the taxonomical changes in the fossil diatom assemblages in combination with both modern diatom data from surrounding ponds and existing micropalaeontological, sedimentological and mineralogical data from the pingo sequence. A diatom based quantitative reconstruction of lake water pH indicates changing lake environments during mid-Wisconsin to early Holocene stages. Mineralogical analyses indicate presence of tephra fallout and its impact on fossil diatom communities. Our comparison of modern and fossil diatom communities shows the highest floristic similarity of modern polygon ponds to the corresponding initial (shallow water) development stages of thermolcarst lakes. We conclude, that mid-Wisconsin thermokarst processes in the study area could establish during relatively warm interstadial climate conditions accompanied by increased precipitation due to approaching coasts, while still high continentality and hence high seasonal temperature gradients led to warm summers in the central part of Beringia. (C) 2017 Elsevier B.V. All rights reserved.
Local observations indicate that climate change and shifting disturbance regimes are causing permafrost degradation. However, the occurrence and distribution of permafrost region disturbances (PRDs) remain poorly resolved across the Arctic and Subarctic. Here we quantify the abundance and distribution of three primary PRDs using time-series analysis of 30-m resolution Landsat imagery from 1999 to 2014. Our dataset spans four continental-scale transects in North America and Eurasia, covering similar to 10% of the permafrost region. Lake area loss (-1.45%) dominated the study domain with enhanced losses occurring at the boundary between discontinuous and continuous permafrost regions. Fires were the most extensive PRD across boreal regions (6.59%), but in tundra regions (0.63%) limited to Alaska. Retrogressive thaw slumps were abundant but highly localized (< 10(-5)%). Our analysis synergizes the global-scale importance of PRDs. The findings highlight the need to include PRDs in next-generation land surface models to project the permafrost carbon feedback.
Local observations indicate that climate change and shifting disturbance regimes are causing permafrost degradation. However, the occurrence and distribution of permafrost region disturbances (PRDs) remain poorly resolved across the Arctic and Subarctic. Here we quantify the abundance and distribution of three primary PRDs using time-series analysis of 30-m resolution Landsat imagery from 1999 to 2014. Our dataset spans four continental-scale transects in North America and Eurasia, covering ~10% of the permafrost region. Lake area loss (−1.45%) dominated the study domain with enhanced losses occurring at the boundary between discontinuous and continuous permafrost regions. Fires were the most extensive PRD across boreal regions (6.59%), but in tundra regions (0.63%) limited to Alaska. Retrogressive thaw slumps were abundant but highly localized (<10−5%). Our analysis synergizes the global-scale importance of PRDs. The findings highlight the need to include PRDs in next-generation land surface models to project the permafrost carbon feedback.
Lakes are a ubiquitous landscape feature in northern permafrost regions. They have a strong impact on carbon, energy and water fluxes and can be quite responsive to climate change. The monitoring of lake change in northern high latitudes, at a sufficiently accurate spatial and temporal resolution, is crucial for understanding the underlying processes driving lake change. To date, lake change studies in permafrost regions were based on a variety of different sources, image acquisition periods and single snapshots, and localized analysis, which hinders the comparison of different regions. Here, we present a methodology based on machine-learning based classification of robust trends of multi-spectral indices of Landsat data (TM, ETM+, OLI) and object-based lake detection, to analyze and compare the individual, local and regional lake dynamics of four different study sites (Alaska North Slope, Western Alaska, Central Yakutia, Kolyma Lowland) in the northern permafrost zone from 1999 to 2014. Regional patterns of lake area change on the Alaska North Slope (-0.69%), Western Alaska (-2.82%), and Kolyma Lowland (-0.51%) largely include increases due to thermokarst lake expansion, but more dominant lake area losses due to catastrophic lake drainage events. In contrast, Central Yakutia showed a remarkable increase in lake area of 48.48%, likely resulting from warmer and wetter climate conditions over the latter half of the study period. Within all study regions, variability in lake dynamics was associated with differences in permafrost characteristics, landscape position (i.e., upland vs. lowland), and surface geology. With the global availability of Landsat data and a consistent methodology for processing the input data derived from robust trends of multi-spectral indices, we demonstrate a transferability, scalability and consistency of lake change analysis within the northern permafrost region.
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 similar to 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. (C) 2016 Elsevier Inc. All rights reserved.
Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (<5 m) circum-Arctic Permafrost Region Pond and Lake (PeRL) database recently becoming available. We have used this database to make the first consistent, high-resolution estimation of Arctic waterbody size distributions, with surface areas ranging from 0.0001 km(2) (100 m(2)) to 1 km(2). We found that the size distributions varied greatly across the thirty study regions investigated and that there was no single universal size distribution function (including power-law distribution functions) appropriate across all of the study regions. We did, however, find close relationships between the statistical moments (mean, variance, and skewness) of the waterbody size distributions from different study regions. Specifically, we found that the spatial variance increased linearly with mean waterbody size (R-2 = 0.97, p < 2.2e-16) and that the skewness decreased approximately hyperbolically. We have demonstrated that these relationships (1) hold across the 30 Arctic study regions covering a variety of (bio)climatic and permafrost zones, (2) hold over time in two of these study regions for which multi-decadal satellite imagery is available, and (3) can be reproduced by simulating rising water levels in a high-resolution digital elevation model. The consistent spatial and temporal relationships between the statistical moments of the waterbody size distributions underscore the dominance of topographic controls in lowland permafrost areas. These results provide motivation for further analyses of the factors involved in waterbody development and spatial distribution and for investigations into the possibility of using statistical moments to predict future hydrologic dynamics in the Arctic.
Thermal erosion is a major mechanism of permafrost degradation, resulting in characteristic landforms. We inventory thermo-erosional valleys in ice-rich coastal lowlands adjacent to the Siberian Laptev Sea based on remote sensing, Geographic Information System (GIS), and field investigations for a first regional assessment of their spatial distribution and characteristics. Three study areas with similar geological (Yedoma Ice Complex) but diverse geomorphological conditions vary in valley areal extent, incision depth, and branching geometry. The most extensive valley networks are incised deeply (up to 35 m) into the broad inclined lowland around Mamontov Klyk. The flat, low-lying plain forming the Buor Khaya Peninsula is more degraded by thermokarst and characterized by long valleys of lower depth with short tributaries. Small, isolated Yedoma Ice Complex remnants in the Lena River Delta predominantly exhibit shorter but deep valleys. Based on these hydrographical network and topography assessments, we discuss geomorphological and hydrological connections to erosion processes. Relative catchment size along with regional slope interact with other Holocene relief-forming processes such as thermokarst and neotectonics. Our findings suggest that thermo-erosional valleys are prominent, hitherto overlooked permafrost degradation landforms that add to impacts on biogeochemical cycling, sediment transport, and hydrology in the degrading Siberian Yedoma Ice Complex.
Ice-rich permafrost has been subject to abrupt thaw and thermokarst formation in the past and is vulnerable to current global warming. The ice-rich permafrost domain includes Yedoma sediments that have never thawed since deposition during the late Pleistocene and Alas sediments that were formed by previous thermokarst processes during the Lateglacial and Holocene warming. Permafrost thaw unlocks organic carbon (OC) and minerals from these deposits and exposes OC to mineralization. A portion of the OC can be associated with iron (Fe), a redox-sensitive element acting as a trap for OC. Post-depositional thaw processes may have induced changes in redox conditions in these deposits and thereby affected Fe distribution and interactions between OC and Fe, with knock-on effects on the role that Fe plays in mediating present day OC mineralization. To test this hypothesis, we measured Fe concentrations and proportion of Fe oxides and Fe complexed with OC in unthawed Yedoma and previously thawed Alas deposits. Total Fe concentrations were determined on 1,292 sediment samples from the Yedoma domain using portable X-ray fluorescence; these concentrations were corrected for trueness using a calibration based on a subset of 144 samples measured by inductively coupled plasma optical emission spectrometry after alkaline fusion (R (2) = 0.95). The total Fe concentration is stable with depth in Yedoma deposits, but we observe a depletion or accumulation of total Fe in Alas deposits, which experienced previous thaw and/or flooding events. Selective Fe extractions targeting reactive forms of Fe on unthawed and previously thawed deposits highlight that about 25% of the total Fe is present as reactive species, either as crystalline or amorphous oxides, or complexed with OC, with no significant difference in proportions of reactive Fe between Yedoma and Alas deposits. These results suggest that redox driven processes during past thermokarst formation impact the present-day distribution of total Fe, and thereby the total amount of reactive Fe in Alas versus Yedoma deposits. This study highlights that ongoing thermokarst lake formation and drainage dynamics in the Arctic influences reactive Fe distribution and thereby interactions between Fe and OC, OC mineralization rates, and greenhouse gas emissions.
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 toHolocene 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, C-13(org)) and stable water isotopes (O-18, D, dexcess) 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 (22800 +/- 280cal. a BP, Unit A), which has vertically subsided in places due to subsequent development of a deep thermokarst lake that initiated around 11800cal. a BP (Unit B). At about 9000cal. 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 1060cal. a BP, including post-drainage sediment freezing from the top down to 154cm 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 11800cal. 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.