skip to main content


Title: Spatial heterogeneity and environmental predictors of permafrost region soil organic carbon stocks
Large stocks of soil organic carbon (SOC) have accumulated in the Northern Hemisphere permafrost region, but their current amounts and future fate remain uncertain. By analyzing dataset combining >2700 soil profiles with environmental variables in a geospatial framework, we generated spatially explicit estimates of permafrost-region SOC stocks, quantified spatial heterogeneity, and identified key environmental predictors. We estimated that 1014 − 175 + 194 Pg C are stored in the top 3 m of permafrost region soils. The greatest uncertainties occurred in circumpolar toe-slope positions and in flat areas of the Tibetan region. We found that soil wetness index and elevation are the dominant topographic controllers and surface air temperature (circumpolar region) and precipitation (Tibetan region) are significant climatic controllers of SOC stocks. Our results provide first high-resolution geospatial assessment of permafrost region SOC stocks and their relationships with environmental factors, which are crucial for modeling the response of permafrost affected soils to changing climate.  more » « less
Award ID(s):
1931333
NSF-PAR ID:
10291956
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; « less
Date Published:
Journal Name:
Science Advances
Volume:
7
Issue:
9
ISSN:
2375-2548
Page Range / eLocation ID:
eaaz5236
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Permafrost soils store between 1330 and 1580 Pg carbon (C), which is 3 times the amount of C in global vegetation, almost twice the amount of C in the atmosphere, and half of the global soil organic C pool. Despite the massive amount of C in permafrost, estimates of soil C storage in the high-latitude permafrost region are highly uncertain, primarily due to undersampling at all spatial scales; circumpolar soil C estimates lack sufficient continental spatial diversity, regional intensity, and replication at the field-site level. Siberian forests are particularly undersampled, yet the larch forests that dominate this region may store more than twice as much soil C as all other boreal forest types in the continuous permafrost zone combined. Here we present above- and belowground C stocks from 20 sites representing a gradient of stand age and structure in a larch watershed of the Kolyma River, near Chersky, Sakha Republic, Russia. We found that the majority of C stored in the top 1 m of the watershed was stored belowground (92 %), with 19 % in the top 10 cm of soil and 40 % in the top 30 cm. Carbon was more variable in surface soils (10 cm; coefficient of variation (CV)  =  0.35 between stands) than in the top 30 cm (CV  =  0.14) or soil profile to 1 m (CV  =  0.20). Combined active-layer and deep frozen deposits (surface – 15 m) contained 205 kg C m−2 (yedoma, non-ice wedge) and 331 kg C m−2 (alas), which, even when accounting for landscape-level ice content, is an order of magnitude more C than that stored in the top meter of soil and 2 orders of magnitude more C than in aboveground biomass. Aboveground biomass was composed of primarily larch (53 %) but also included understory vegetation (30 %), woody debris (11 %) and snag (6 %) biomass. While aboveground biomass contained relatively little (8 %) of the C stocks in the watershed, aboveground processes were linked to thaw depth and belowground C storage. Thaw depth was negatively related to stand age, and soil C density (top 10 cm) was positively related to soil moisture and negatively related to moss and lichen cover. These results suggest that, as the climate warms, changes in stand age and structure may be as important as direct climate effects on belowground environmental conditions and permafrost C vulnerability. 
    more » « less
  2. Abstract

    Permafrost thaw is typically measured with active layer thickness, or the maximum seasonal thaw measured from the ground surface. However, previous work has shown that this measurement alone fails to account for ground subsidence and therefore underestimates permafrost thaw. To determine the impact of subsidence on observed permafrost thaw and thawed soil carbon stocks, we quantified subsidence using high‐accuracy GPS and identified its environmental drivers in a permafrost warming experiment near the southern limit of permafrost in Alaska. With permafrost temperatures near 0°C, 10.8 cm of subsidence was observed in control plots over 9 years. Experimental air and soil warming increased subsidence by five times and created inundated microsites. Across treatments, ice and soil loss drove 85–91% and 9–15% of subsidence, respectively. Accounting for subsidence, permafrost thawed between 19% (control) and 49% (warming) deeper than active layer thickness indicated, and the amount of newly thawed carbon within the active layer was between 37% (control) and 113% (warming) greater. As additional carbon thaws as the active layer deepens, carbon fluxes to the atmosphere and lateral transport of carbon in groundwater could increase. The magnitude of this impact is uncertain at the landscape scale, though, due to limited subsidence measurements. Therefore, to determine the full extent of permafrost thaw across the circumpolar region and its feedback on the carbon cycle, it is necessary to quantify subsidence more broadly across the circumpolar region.

     
    more » « less
  3. Abstract

    The magnitude of future emissions of greenhouse gases from the northern permafrost region depends crucially on the mineralization of soil organic carbon (SOC) that has accumulated over millennia in these perennially frozen soils. Many recent studies have used radiocarbon (14C) to quantify the release of this “old” SOC as CO2or CH4to the atmosphere or as dissolved and particulate organic carbon (DOC and POC) to surface waters. We compiled ~1,90014C measurements from 51 sites in the northern permafrost region to assess the vulnerability of thawing SOC in tundra, forest, peatland, lake, and river ecosystems. We found that growing season soil14C‐CO2emissions generally had a modern (post‐1950s) signature, but that well‐drained, oxic soils had increased CO2emissions derived from older sources following recent thaw. The age of CO2and CH4emitted from lakes depended primarily on the age and quantity of SOC in sediments and on the mode of emission, and indicated substantial losses of previously frozen SOC from actively expanding thermokarst lakes. Increased fluvial export of aged DOC and POC occurred from sites where permafrost thaw caused soil thermal erosion. There was limited evidence supporting release of previously frozen SOC as CO2, CH4, and DOC from thawing peatlands with anoxic soils. This synthesis thus suggests widespread but not universal release of permafrost SOC following thaw. We show that different definitions of “old” sources among studies hamper the comparison of vulnerability of permafrost SOC across ecosystems and disturbances. We also highlight opportunities for future14C studies in the permafrost region.

     
    more » « less
  4. Conifer forests historically have been resilient to wildfires in part due to thick organic soil layers that regulate combustion and post-fire moisture and vegetation change. However, recent shifts in fire activity in western North America may be overwhelming these resilience mechanisms with potential impacts for energy and carbon exchange. Here, we quantify the long-term recovery of the organic soil layer and its carbon pools across 511 forested plots. Our plots span ~ 140,000 km2 across two ecozones of the Northwest Territories, Canada, and allowed us to investigate the impacts of time-after-fire, site moisture class, and dominant canopy type on soil organic layer thickness and associated carbon stocks. Despite thinner soil organic layers in xeric plots immediately after fire, these drier stands supported faster post-fire recovery of the soil organic layer than in mesic plots. Unlike xeric or mesic stands, post-fire soil carbon accumulation rates in hydric plots were negligible despite wetter forested plots having greater soil organic carbon stocks immediately post-fire compared to other stands. While permafrost and high-water tables inhibit combustion and maintain thick organic soils immediately after fire, our results suggest that these wet stands are not recovering their pre-fire carbon stocks on a century timescale. We show that canopy conversion from black spruce to jack pine or deciduous dominance could reduce organic soil carbon stocks by 60–80% depending on stand age. Our two main findings—decreasing organic soil carbon storage with increasing deciduous cover and the lack of post-fire SOL recovery in hydric sites—have implications for the turnover time of carbon stocks in the western boreal forest region and also will impact energy fluxes by controlling albedo and surface soil moisture. 
    more » « less
  5. Abstract

    Estimates of soil organic carbon (SOC) stocks are essential for many environmental applications. However, significant inconsistencies exist in SOC stock estimates for the U.S. across current SOC maps. We propose a framework that combines unsupervised multivariate geographic clustering (MGC) and supervised Random Forests regression, improving SOC maps by capturing heterogeneous relationships with SOC drivers. We first used MGC to divide the U.S. into 20 SOC regions based on the similarity of covariates (soil biogeochemical, bioclimatic, biological, and physiographic variables). Subsequently, separate Random Forests models were trained for each SOC region, utilizing environmental covariates and SOC observations. Our estimated SOC stocks for the U.S. (52.6 ± 3.2 Pg for 0–30 cm and 108.3 ± 8.2 Pg for 0–100 cm depth) were within the range estimated by existing products like Harmonized World Soil Database, HWSD (46.7 Pg for 0–30 cm and 90.7 Pg for 0–100 cm depth) and SoilGrids 2.0 (45.7 Pg for 0–30 cm and 133.0 Pg for 0–100 cm depth). However, independent validation with soil profile data from the National Ecological Observatory Network showed that our approach (R2 = 0.51) outperformed the estimates obtained from Harmonized World Soil Database (R2 = 0.23) and SoilGrids 2.0 (R2 = 0.39) for the topsoil (0–30 cm). Uncertainty analysis (e.g., low representativeness and high coefficients of variation) identified regions requiring more measurements, such as Alaska and the deserts of the U.S. Southwest. Our approach effectively captures the heterogeneous relationships between widely available predictors and the current SOC baseline across regions, offering reliable SOC estimates at 1 km resolution for benchmarking Earth system models.

     
    more » « less