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.
Thermokarst lakes are prevalent in Arctic coastal lowland regions and sublake permafrost degradation and talik development contributes to greenhouse gas emissions by tapping the large permafrost carbon pool. Whereas lateral thermokarst lake expansion is readily apparent through remote sensing and shoreline measurements, sublake thawed sediment conditions and talik growth are difficult to measure. Here we combine transient electromagnetic surveys with thermal modeling, backed up by measured permafrost properties and radiocarbon ages, to reveal closed-talik geometry associated with a thermokarst lake in continuous permafrost. To improve access to talik geometry data, we conducted surveys along three transient electromagnetic transects perpendicular to lakeshores with different decadal-scale expansion rates of 0.16, 0.38, and 0.58m/year. We modeled thermal development of the talik using boundary conditions based on field data from the lake, surrounding permafrost and a borehole, independent of the transient electromagnetics. A talik depth of 91m was determined from analysis of the transient electromagnetic surveys. Using a lake initiation age of 1400years before present and available subsurface properties the results from thermal modeling of the lake center arrived at a best estimate talk depth of 80m, which is on the same order of magnitude as the results from the transient electromagnetic survey. Our approach has provided a noninvasive estimate of talik geometry suitable for comparable settings throughout circum-Arctic coastal lowland regions.
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.
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
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
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.
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.
Arctic tundra ecosystems have experienced unprecedented change associated with climate warming over recent decades. Across the Pan-Arctic, vegetation productivity and surface greenness have trended positively over the period of satellite observation. However, since 2011 these trends have slowed considerably, showing signs of browning in many regions. It is unclear what factors are driving this change and which regions/landforms will be most sensitive to future browning. Here we provide evidence linking decadal patterns in arctic greening and browning with regional climate change and local permafrost-driven landscape heterogeneity. We analyzed the spatial variability of decadal-scale trends in surface greenness across the Arctic Coastal Plain of northern Alaska (similar to 60,000 km(2)) using the Landsat archive (1999-2014), in combination with novel 30 m classifications of polygonal tundra and regional watersheds, finding landscape heterogeneity and regional climate change to be the most important factors controlling historical greenness trends. Browning was linked to increased temperature and precipitation, with the exception of young landforms (developed following lake drainage), which will likely continue to green. Spatiotemporal model forecasting suggests carbon uptake potential to be reduced in response to warmer and/or wetter climatic conditions, potentially increasing the net loss of carbon to the atmosphere, at a greater degree than previously expected.