The decomposition of thawing permafrost organic matter (OM) to the greenhouse gases (GHG) carbon dioxide (CO2) and methane forms a positive feedback to global climate change. Data on in situ GHG fluxes from thawing permafrost OM are scarce and OM degradability is largely unknown, causing high uncertainties in the permafrost-carbon climate feedback. We combined in situ CO2 and methane flux measurements at an abrupt permafrost thaw feature with laboratory incubations and dynamic modeling to quantify annual CO2 release from thawing permafrost OM, estimate its in situ degradability and evaluate the explanatory power of incubation experiments. In July 2016 and 2019, CO2 fluxes ranged between 0.24 and 2.6 g CO2-C m(-2) d(-1). Methane fluxes were low, which coincided with the absence of active methanogens in the Pleistocene permafrost. CO2 fluxes were lower three years after initial thaw after normalizing these fluxes to thawed carbon, indicating the depletion of labile carbon. Higher CO2 fluxes from thawing Pleistocene permafrost than from Holocene permafrost indicate OM preservation for millennia and give evidence that microbial activity in the permafrost was not substantial. Short-term incubations overestimated in situ CO2 fluxes but underestimated methane fluxes. Two independent models simulated median annual CO2 fluxes of 160 and 184 g CO2-C m(-2) from the thaw slump, which include 25%-31% CO2 emissions during winter. Annual CO2 fluxes represent 0.8% of the carbon pool thawed in the surface soil. Our results demonstrate the potential of abrupt thaw processes to transform the tundra from carbon neutral into a substantial GHG source.
Arctic warming is leading to substantial changes to permafrost including rapid degradation of ice and ice-rich coasts and riverbanks. In this study, we present and evaluate a high spatiotemporal resolution three-year time series of X-Band microwave satellite data from the TerraSAR-X (TSX) satellite to quantify cliff-top erosion (CTE) of an ice-rich permafrost riverbank in the central Lena Delta. We apply a threshold on TSX backscatter images and automatically extract cliff-top lines to derive intra- and inter-annual CTE. In order to examine the drivers of erosion we statistically compare CTE with climatic baseline data using linear mixed models and analysis of variance (ANOVA). Our evaluation of TSX-derived CTE against annual optical-derived CTE and seasonal in situ measurements showed good agreement between all three datasets. We observed continuous erosion from June to September in 2014 and 2015 with no significant seasonality across the thawing season. We found the highest net annual cliff-top erosion of 6.9 m in 2014, in accordance with above-average mean temperatures and thawing degree days as well as low precipitation. We found high net annual erosion and erosion variability in 2015 associated with moderate mean temperatures but above average precipitation. According to linear mixed models, climate parameters alone could not explain intra-seasonal erosional patterns and additional factors such as ground ice content likely drive the observed erosion. Finally, mean backscatter intensity on the cliff surface decreased from -5.29 to -6.69 dB from 2013 to 2015, respectively, likely resulting from changes in surface geometry and properties that could be connected to partial slope stabilization. Overall, we conclude that X-Band backscatter time series can successfully be used to complement optical remote sensing and in situ monitoring of rapid tundra permafrost erosion at riverbanks and coasts by reliably providing information about intra-seasonal dynamics.
Warming of the arctic landscape results in permafrost thaw, which causes ground subsidence or thermokarst. Thermokarst formation on hillslopes leads to the formation of thermal erosion features that dramatically alter soil properties and likely affect soil carbon emissions, but such features have received little study in this regard. In order to assess the magnitude and persistence of altered emissions, we use a space-for-time substitution (thaw slump chronosequence) to quantify and compare peak growing season soil carbon dioxide (CO2) fluxes from undisturbed tundra, active, and stabilized thermal erosion features over two seasons. Measurements of soil temperature and moisture, soil organic matter, and bulk density are used to evaluate the factors controlling soil CO2 emissions from each of the three chronosequence stages. Soil CO2 efflux from the active slump is consistently less than half that observed in the undisturbed tundra or stabilized slump (1.8 versus 5.2 g CO2-C m(-2) d(-1) in 2011; 0.9 versus 3.2 g CO2-C m(-2) d(-1) in 2012), despite soil temperatures on the floor of the active slump that are 10-15 degrees C warmer than the tundra and stabilized slump. Environmental factors such as soil temperature and moisture do not exert a strong control on CO2 efflux, rather, local soil physical and chemical properties such as soil organic matter and bulk density, are strongly and inversely related among these chronosequence stages (r(2) = 0.97), and explain similar to 50% of the variation in soil CO2 efflux. Thus, despite profound soil warming and rapid exposure of buried carbon in the active slump, the low organic matter content, lack of stable vegetation, and large increases in the bulk densities in the uppermost portion of active slump soils (up to similar to 2.2 g(-1) cm(-3)) appear to limit CO2 efflux from the active slump. Future studies should assess seasonal fluxes across these features and determine whether soil CO2 fluxes from active features with high organic content are similarly low.