Permafrost deposits in the Beringian Yedoma region store large amounts of organic carbon (OC). Walter Anthony et al. (2014) describe a previously unrecognized pool of 159 Pg OC accumulated in Holocene thermokarst sediments deposited in Yedoma region alases (thermokarst depressions). They claim that these alas sediments increase the previously recognized circumpolar permafrost peat OC pool by 50 %. It is stated that previous integrated studies of the permafrost OC pool have failed to account for these deposits because the Northern Circumpolar Soil Carbon Database (NCSCD) is biased towards non-alas field sites and that the soil maps used in the NCSCD underestimate coverage of organic permafrost soils. Here we evaluate these statements against a brief literature review, existing data sets on Yedoma region soil OC storage and independent field-based and geospatial data sets of peat soil distribution in the Siberian Yedoma region. Our findings are summarized in three main points. Firstly, the sediments described by Walter Anthony et al. (2014) are primarily mineral lake sediments and do not match widely used international scientific definitions of peat or organic soils. They can therefore not be considered an addition to the circumpolar peat carbon pool. We also emphasize that a clear distinction between mineral and organic soil types is important since they show very different vulnerability trajectories under climate change. Secondly, independent field data and geospatial analyses show that the Siberian Yedoma region is dominated by mineral soils, not peatlands. Thus, there is no evidence to suggest any systematic bias in the NCSCD field data or maps. Thirdly, there is spatial overlap between these Holocene thermokarst sediments and previous estimates of permafrost soil and sediment OC stocks. These carbon stocks were already accounted for by previous studies and they do not significantly increase the known circumpolar OC pool. We suggest that these inaccurate statements made in Walter Anthony et al. (2014) mainly resulted from misunderstandings caused by conflicting definitions and terminologies across different geoscientific disciplines. A careful cross-disciplinary review of terminologies would help future studies to appropriately harmonize definitions between different fields.
Soils and sediments of the northern permafrost region have accumulated large
stocks of organic carbon (OC) over millennia (Tarnocai et al., 2009). As the
global climate warms there is a concern that thawing permafrost will expose
soil organic matter (SOM) that was previously protected in permafrost to
decomposition, causing a positive permafrost-carbon feedback to climate
(Schuur et al., 2008, 2015). Hugelius et al. (2014) provide the most recent
integrated estimate of Northern circumpolar permafrost region soil and
sediment OC stocks with total stocks estimated at 1307 Pg and a 95 %
confidence interval of 1140–1476 Pg. Of this roughly 800 Pg is perennially
frozen with the remainder stored in active layer or talik deposits. A
substantial part of the perennially frozen OC is stored in the Beringian
Yedoma region with estimated permafrost deposit OC stocks of 213 Pg with an
uncertainty range of 164–267 Pg. Schirrmeister et al. (2013) provide an in
depth discussion and review on various aspects of these deposits. Schuur et
al. (2015), in a recent review of the permafrost carbon feedback, highlight
that there is considerable spread in estimates of Yedoma region permafrost
OC stocks. In a study describing the Holocene C dynamics of Siberian
thermokarst lakes Walter Anthony et al. (2014) estimate a pool of
Here we examine these important statements by evaluating the findings and data presented by Walter Anthony et al. (2014) against (1) a brief review of vulnerability to climatic changes and scientific definitions of peat, peatlands, organic soils and thermokarst sediments, (2) independent field data as well as independent geospatial databases showing the extent of organic soils and/or peatlands in the Siberian Yedoma region and (3) by analyzing the spatial overlap between these new estimates and existing data sets of Yedoma region soil and sediment OC storage.
Walter Anthony et al. (2014) claim that 159 Pg of OC accumulated in deep Holocene thermokarst deposits across the Yedoma region increase the previously recognized permafrost peat OC pool by 50 %. We argue that the use of imprecise terminology has caused misleading comparisons in relation to previous stock estimates. These Holocene thermokarst deposits do not meet the criteria of peat (or organic soils) used in any regional or circumpolar peat carbon stock study. Therefore they cannot be claimed to increase peat carbon stocks. They simply increase the stock of alas sediments known to be of Holocene age.
We emphasize that the properties of mineral and organic soil material are very different and the distinction is especially important in permafrost regions where studies have consistently shown that organic and mineral soils differ both in their vulnerability to thaw and in the potential post-thaw lability of soil OM. Correct classification of organic and mineral soil material is not a mere issue of semantics or putting a different label on something depending on your scientific background. For example, due to distinct differences in soil thermal properties, organic soils are much less vulnerable to active layer deepening under climate warming than mineral soils (Shur and Jorgenson, 2007; Jorgenson et al., 2010). High-resolution modeling of active layer dynamics from a Russian low-Arctic site showed that organic soil is projected to remain stable until the end of this century while near-surface permafrost degraded in mineral soils (Hugelius et al., 2011). Organic soils also show different vulnerabilities to thermokarst. Thick surface O-horizons can reduce lateral expansion rates of thermokarst (Jorgenson and Osterkamp, 2005) and modelling studies suggest that thermokarst lake taliks formed into organic soils are shallower than their mineral counterparts (West and Plug, 2008). Sjöberg et al. (2013) suggest that thermokarst lake formation and orientation in peatland terrain may partly be controlled by different processes than for mineral soil thermokarst. They also demonstrated that peat substrate thermokarst lake shorelines display more pronounced and heterogeneous erosion patterns than mineral substrate shorelines, both in shoreline morphology and lake geometry. Harden et al. (2006) also describe multiple feedbacks between the thickness of surface organic soil horizons and the vulnerability of ecosystems to combustion by wildfires, where deep organic layers could often preserve thermal and biological properties of soils through repeated fire cycles.
While organic soils are thus less vulnerable to permafrost thaw and
combustion than mineral soils, other studies have demonstrated that SOM in
organic soils is typically less decomposed than in mineral soils, and thus
assumed to be more vulnerable to microbial decomposition. Through
comprehensive analyses of Siberian permafrost sediments Strauss et al. (2015)
showed that high OC % content is associated with less degraded SOM,
as indicated by multiple geochemical proxies. In sub-Arctic tundra, SOM in
peatlands has been shown to be significantly less degraded than mineral soil
SOM (Hugelius et al., 2012; Routh et al., 2014). Incubation studies have
also confirmed these findings. In a circumpolar incubation synthesis, the
fractional loss of initial soil OC was a factor of 2 to 4 higher from organic
soils compared to mineral soils (over 50 incubation years at 5
In light of these studies showing clear differences in the properties and potential vulnerabilities of mineral and organic soils it is evident that clear definitions and distinctions are needed to properly assess and predict the response of these vulnerable landscapes under a changing climate. Below we provide a brief review of different definitions and classifications currently used in studies of periglacial terrain.
Across different scientific disciplines (and countries) the definition of what is peat varies. A commonly used definition states that peat is sedentarily accumulated material consisting of at least 30 % (dry weight) dead organic material while peatlands are areas (with or without vegetation) with a naturally accumulated peat layer (Joostens and Clark, 2002). Many studies have employed a minimum depth criterion of the surface peat layer to the definition of peatland, most frequently 30 cm (Kivinen and Pakarinen, 1981; Lappalainen 1996; Joostens and Clark, 2002). The Canadian definition of an organic wetland (or peatland) includes a depth of organic soil material (of 17 % OC or 30 % organic material) of at least 40 cm (National Wetlands Working Group, 1997).
Soil classification systems define organic soil material (or peat) based on organic carbon content, while the thickness of organic soil material in the upper soil column determines whether a soil is primarily considered to be a mineral soil or an organic soil. The U.S. soil taxonomy (Soil Survey Staff, 2010) and the World Reference Base for Soil Resources (IUSS Working Group WRB, 2007) defines waterlogged soil with more than 12–18 % OC (dry weight; range depending on clay content) as organic soil material while the Canadian System of Soil Classification (Soil Classification Working Group, 1998) defines soil with more than 17 % OC (or 30 % organic material; dry weight) as organic soil material. All these soil classification systems define a soil as an organic soil if there is 40 cm or more of accumulated organic soil material in the upper soil column (the Canadian system employs 60 cm for highly fibric moss-peat).
The literature describing sediments of thermokarst basins and lakes includes many different definitions of different facies or deposit types. These definitions are often not based on quantified physical or chemical properties of sediments, but rather reflect descriptive characteristics and the environments in which they formed. In addition to in situ peat, previous studies have described organic rich sedimentary thermokarst facies such as the following: (1) “detrital peat” described as layered organic deposits formed on beaches or in shallow waters (Murton, 1996) or as lee-shore deposits (Hopkins and Kid, 1988); (2) “organic rich silts” (or “lacustrine organic silt”) where primarily mineral lake sediments are interspersed with sedentary or allochthonous detrital organic sediments layers (Murton, 1996; Kanevskiy et al., 2014) and (3) “Mud/muddy peat” which differs from detrital peat based on a higher mud content. These deposits may contain blocks of peat or other materials and typically form thick sediment layers in deep water thermokarst lake environments by suspension settling of fine and/or low-density material (Hopkins and Kidd, 1988; Murton, 1996).
The bulk of the Holocene OC described by Walter Anthony et al. (2014) has accumulated in sediment facies these authors descriptively call “Stratified muddy peat”. The authors state that this facies corresponds to strata that previous authors have called “Mud/muddy peat”. These facies are described as deep-water lake sediments, predominantly of minerogenic origin and with an OC content of only 3–4 % by weight (Walter Anthony et al., 2014; Fig. 2 and extended data Table 2). Their use of terminology is in line with previous studies of thermokarst sedimentary facies. But the classification, origin and properties of these deposits are clearly very different from pedologically defined peat as being a primarily organic material, usually of terrestrial or shallow water origin.
Walter Anthony et al. (2014) describe a pool of Holocene OC which was previously unrecognized because earlier studies had not systematically accounted for alas deposits. Here we evaluate these statements and present independent evidence showing that their argumentation is based on flawed assumptions. We find no evidence of systematic biases in the data sets on which earlier studies were based.
The Beringian Yedoma region can be subdivided into areas of intact Yedoma
(ca. 30 % by area), areas that have been affected by thermokarst and
subsequently re-aggraded permafrost (56 %) and areas of open water
(14 %) which are commonly underlain by taliks (Strauss et al., 2013). The
study by Walter Anthony et al. (2014) uses an identical spatial subdivision
of this region but with different data and computational methods to estimate
the volume and OC stocks of the various sediments and deposit types in the
region. This includes a thermokarst-basin Holocene carbon pool (
We further evaluate these claims of systematic underrepresentation of organic soils in the NCSCD maps and pedon databases against independent inventories of geospatial data sets and field data. To provide independent estimates of Siberian Yedoma region peatland coverage four different geospatial data sets were used (Nilsson et al., 2002; Bartalev et al., 2003; Lehner and Döll, 2004; Arino et al., 2012). Thematic classes that corresponded to peatlands were identified and their respective coverage quantified. The independent field validation sites are all located in alases or thermoerosional gullies from across the Siberian Yedoma region and were classified and sampled using a transect-based semi-random approach during field campaigns in August (2010 and 2013). For detailed method descriptions and calculations we refer to the online supplementary materials.
Both the geospatial data sets and field inventory data show a limited extent
of organic soils in the Siberian Yedoma region (Fig. 1). The mapped Histel
coverage in the Siberian Yedoma region in the NCSCD is 9 % (Supplement Fig. S1).
This is comparable to peatland coverage estimated from independent
geospatial databases of 3–6 % (Fig. 1). It is notable that the degree of
overlap between independent data sets is limited, indicative of difficulties
with classifications and class definitions when mapping peatland extent
(Fig. S1). A spatial overlay analyses of regional land-cover and wetland
characterization maps (Nilsson et al., 2002; Stolbovoi, 2002) suggest that
Overview of field sites and estimated coverage of peatlands in the Siberian Yedoma region. Graduated colors within the region show coverage of peatlands in four global and/or regional map products that are independent from the NCSCD (see online supplementary material for detailed methods). The coverage is shown cumulatively so that the colors reflect how many of the four products that map peatlands in any given location. Points show locations of the Holocene alas profiles used by Walter Anthony et al. (2014) as well as independent soil profiles for validation (classified as mineral soils or peatlands). All of the independent validation points are known to be located in alases or thermoerosional gullies. Extent of the Siberian Yedoma region digitized from Grosse et al. (2013).
Conceptual diagram illustrating how organic soil/sediment C in
Yedoma region alases is described and estimated by
Our independent compilation of field sites located in alases or
thermoerosional gullies from across the Siberian Yedoma region reveals that
16 % of sites are peatlands (Fig. 1; 9 out of 49 sites). The surface peat
depth of these nine peatland sites was
The two main processes of global peatland formation (and expansion) are paludification or terrestrialization. Terrestrialization describes peatlands formed via gradual in-filling of water bodies. Paludification is the process by which peatlands expand into other established terrestrial ecosystems. Paludification is considered the most common form of high-latitude peatland formation (Charman, 2002; Kuhry and Turunen, 2006). The bulk of the Holocene alas deposits described by Walter Anthony et al. (2014) were formed through a terrestrialization process in combination with permafrost dynamics (sometimes causing rapid drainage). The characteristics of peatland sites from across the Yedoma region show that paludification of terrestrial alas ecosystems has also contributed to local peatland formation (e.g., facies F1 in Walter Anthony et al., 2014 see also Palmtag et al., 2015 and Weiss et al., 2015).
All of these combined lines of evidence support an interpretation that peatlands are locally present in alases of the Siberian Yedoma region, but rarely cover large surfaces. We recognize that the maps and pedon data set of the NCSCD are highly generalized, but find no support to the claim that they are systematically biased to non-alas soils. We conclude that the Siberian Yedoma region alases are dominated by mineral soils, often formed on parent material of reworked yedoma or lacustrine sediments. This interpretation is also supported by previous scientific studies from this region (e.g., Czudek and Demek, 1970; Veremeeva and Gubin, 2009; Wetterich et al., 2009; Schirrmeister et al., 2011; Morgenstern et al., 2013; Siewert et al., 2015).
Spatial overlap between different studies of soil carbon stocks may mislead data users and cause significant errors in estimates. Walter Anthony et al. (2014) claim that the pool of 159 Pg Holocene OC in Yedoma region alases increases the previously recognized circumpolar permafrost peat OC pool by 50 %. Here we show that these sediments were already accounted for by previous studies.
The calculations of overlap in soil carbon stocks between different estimates and data sets for the Siberian Yedoma region are based on data on soil and/or sediment carbon stocks from Tarnocai et al. (2009), Hugelius et al. (2013a, b, 2014) and Walter Anthony et al. (2014). By using the reported depth ranges and soil carbon densities of the different studies, the overlap between estimates has been calculated following the same methods used in the original studies. We refer to the online method section for more details on the calculations.
Previous integrated estimates of carbon stocks in the Beringian Yedoma
region (Tarnocai et al., 2009; Hugelius et al., 2014) are based on soil maps
linked to field-based soil data for the upper 3 m and generalized
estimates of Yedoma region deposits for deeper deposits (Zimov et al., 2006;
Strauss et al., 2013). The Holocene thermokarst deposits described by Walter
Anthony et al. (2014) overlap these previous estimates in space, but they
differ in their characterization of the sediment (Fig. 2). An important
difference compared to previous studies is that Walter Anthony et al. (2014)
include 24 Pg carbon in Holocene deposits assumed to occur in taliks
(perennially thawed ground) under present day lakes and rivers. We recognize
that these estimates are new, but they are also outside the scope of the
studies by Tarnocai et al. (2009) and Hugelius et al. (2014) as they are per
definition not terrestrial soils, nor are they permafrost deposits. Out of the 159 Pg of
Holocene alas carbon reported by Walter Anthony et al. (2014), this leaves 135 Pg
of Holocene carbon to be reconciled with previous estimates for
soil and/or sediment that occupy the same physical space. For the upper 3 m, Walter Anthony et al. (2014) estimate 76 Pg of Holocene carbon. This
overlaps soils from previous estimates with carbon stocks of 53–58 Pg in 0–3 m
soils from the NCSCD (range based on different versions of the NCSCD from
Tarnocai et al., 2009; Hugelius et al., 2013a, b, 2014). This comparison
would result in a
For alas deposits below 3 m, the estimate by Walter Anthony et al. (2014)
includes 60 Pg of Holocene OC and 155 Pg of Pleistocene OC which
overlaps with estimates of
We conclude that Holocene OC stocks in Siberian Yedoma region alases which overlap
estimates from previous studies are primarily stored in mineral soils and
lacustrine sediments rather than peat and do not increase estimates of
circumpolar permafrost peat carbon stocks. There is no evidence or reasoning
to suggest that these deposits increase the northern peatland pool or that
the NCSCD is systematically biased against upland soils. In fact, the
differences between the estimates of Hugelius et al. (2014) and Walter
Anthony et al. (2014) are rather small. If storage in taberites and
subaqueous sediments is accounted for, the difference in estimated Yedoma
region alas OC stocks is only
We emphasize that our concerns regarding use of terminology and spatial overlap of estimates in the discussed study in no way affects the validity of their other important findings regarding Holocene carbon dynamics of these ecosystems. It is relevant and important to contrast the Holocene accumulation of carbon in alas sediments to that estimated for peatlands. We attribute the misunderstandings to confusing overlap between terminologies in the respective fields of science that study soils and sediments in periglacial landscapes. We suggest that a careful and exhaustive review of these terminologies would help future studies to harmonize classifications and definitions. The need for reconciliation of terminologies is emphasized by accumulating evidence that the differing properties of mineral and organic soil affect their vulnerability under future climatic changes.
For more detailed descriptions of methods and the data sets used in this paper, please consult the supplementary materials.
The GIS database and pedon data for the NCSCDv2 (Hugelius et al., 2013b) is avaiable via:
We are grateful to K. M. Walter Anthony, S. A. Zimov, G. Grosse, M. C. Jones, P. M. Anthony, F. S. Chapin III, J. C. Finlay, M. C. Mack, S. Davydov, P. Frenzel and S. Frolking for constructive scientific correspondence regarding the issues discussed in this manuscript, both prior to our submission to this journal and during the open discussion stage. We are grateful to three anonymous reviewers for very constructive comments which helped us write a better and more interesting manuscript. We acknowledge J. Palmtag, N. Weiss and M. Siewert for contributing to soil classification of the validation field sites. G. Hugelius acknowledges funding from the Swedish Research Council (grant number E0689701) and the EU JPI COUP project. Edited by: R. Conant