The sensitivity of soil carbon and nitrogen to warming is a major uncertainty in projections of climate. However, previous studies about soil organic carbon (SOC) stocks and potential emission predominantly concentrated on the shallow soil layer in high latitude ecosystems. In this study, we analyzed the SOC, total nitrogen (TN) and soil inorganic carbon (SIC) stocks, C/N ratios, and stable carbon isotope (δ13C) in the active layer and permafrost layer on the Eboling Mountain in the upper reach of Heihe River basin, northwestern China. Our results showed that the average stocks of SOC, TN, and SIC in permafrost layer above soil parent materials (71.7 kg m-2, 8.0 kg m-2, 34.7 kg m-2) were much higher than those in the active layer (44.3 kg m-2, 5.3 kg m-2, 12.2 kg m-2). The δ13C pattern in the soil profiles indicated that historical drainage conditions and pedogenesis were important factors in determining soil organic matter (SOM) stocks in this permafrost region. The δ13C and C/N ratios of the transient layer and some layers of permafrost implied that the degradation of SOM was different. These results highlight that carbon and nitrogen in permafrost regions with Alpine Kobresia meadow could make significant contribution to China's terrestrial carbon cycle.
Introduction
Recently, terrestrial ecosystem carbon and nitrogen have received increased attention worldwide because emissions of greenhouse gases into the atmosphere play an important role in driving global warming (Jobbágy and Jackson, 2000; Wynn et al., 2006; Fu et al., 2010). The northern permafrost regions contain approximately 1672 Pg C, which is roughly double the amount of carbon currently in the atmosphere (Tarnocai et al., 2009). It had been proposed that decrease of soil organic carbon (SOC) in permafrost can lead to an increase of greenhouse gases input to the atmosphere (Schuur et al., 2009; Dorrepaal et al., 2009; Smith and Fang, 2010; Schaefer et al., 2011). However, most studies of SOC stocks and its potential emission have focused on soils and upper permafrost (Ping et al., 2008a, 2008b; Zimov et al., 2006), while the deep permafrost carbon (>3 m) remains largely unknown. Especially, most distribution of SOC was studied in high latitude ecosystems, however, magnitude and spatial patterns of carbon and nitrogen in the mid-latitude regions remain largely unknown. Therefore, accurate estimates of carbon stocks, distribution, and its potential emission in deep permafrost of the mid-latitude ecosystems are critical for predicting feedbacks between soil carbon and global warming.
As the largest carbon pool in the terrestrial biosphere, soil carbon stocks consist of organic and inorganic components. However, previous studies dominantly concentrated on SOC stocks, and little is known about the magnitude and patterns of soil inorganic carbon (SIC). The SIC stocks in the top 1 m of soil depth are about 930–1738 Pg (Schlesinger, 2002; Wang et al., 2012). SIC stocks in the Tibetan alpine grasslands are about 2.1 times the corresponding SOC stocks (Yang et al., 2010). Total SIC stocks in China are approximately 60 Pg C, representing 1/20 of the global SIC pool (Wu et al., 2009). The SIC was mostly stored below the 1 m of soil depth (Mikhailova and Post, 2006). Small changes in SOC and SIC may radically alter the carbon balance (Li et al., 2007; Gleixner et al., 2009). SOC and SIC can interact with each other—for example, SOC accumulation with carbonates may induce the dissolution of SIC (Duan et al., 1999)—while respiratory CO2 from SOM decomposition and root respiration can be used in secondary SIC formation (Entry et al., 2004). Therefore, a concurrent study on the changes in the SIC and SOC in permafrost is important for carbon balance study.
Not only are carbon stocks affected by permafrost degradation, but also soil nitrogen content can play an important role in permafrost carbon cycle (Harden et al., 2012). The rate of SOC turnover is strongly linked to nitrogen availability. However, little work has been conducted on nitrogen stocks within the soil profile in comparison to carbon stocks (Crowe et al., 2004). It has long been recognized that nitrogen limitations often constrain carbon accumulations in high latitude ecosystems (Tamm et al., 1982). High nitrogen availability may also strongly increase microbial growth, thus the rate of SOM decomposition (Sistla et al., 2012). Several modeling studies and empirical experiments have identified nitrogen as a potential regulator of SOM decomposition (Churchland et al., 2010; Lavoie et al., 2011). Vertical distribution of nitrogen contents and C/N ratios at different depth in permafrost plays a key role in how soils respond to environmental changes (Callesen et al., 2007). Soil C/N ratios usually are used as a scalar parameter to predict nitrous oxide and methane emissions (Andersson et al., 2012).
Stable carbon isotope analyses of peat and other geological materials is an important tool for gaining insight into biogeochemical processes involved in SOM degradation and preservation (Kracht and Gleixner, 2000), and also in paleoclimate studies (Jones et al., 2010). Detailed sampling of the entire soil profile is crucial to any of these interpretations because δ13C depth profiles reflect systematic fractionation processes of SOC (Becker-Heidmann and Scharpenseel, 1986).
The Eboling Mountain is located in the Qilian Mountains, northwestern China, and it has experienced severe degradation in recent years. During the International Polar Year (ca. 2007–2008), mean annual ground temperature at depth of 6 m ranged from -3.2 °C to 0.2 °C (Zhao et al., 2010). Few studies about carbon in permafrost were carried out in the alpine permafrost area in the middle latitudes (Wu et al., 2012). Thus it is important to study the deep permafrost here. The aim of this study is to understand the characteristics of soil carbon and nitrogen, especially the dynamics of SOC and SIC in deep permafrost of mid-latitude ecosystem by determining the distribution of SOC, total N, C/N ratios, and δ13C with depth, and their possible relationships with soil properties. It is important in predicting the response of alpine permafrost in the middle latitudes to global warming.
Field Sampling and Laboratory Analyses
SITE DESCRIPTION AND PERMAFROST SAMPLING
Sampling sites were located in the permafrost zone on the Eboling Mountain in the upper reach of Heihe River basin in the Eastern Qilian Mountain, Northwestern China (Site A: 37°59.872′N, 100°54.977′E, 3700 m), (Site B: 38°00.196′N, 100°54.411′E, 3615 m). The slope of Site A was 11° (north by west 5°) and that of Site B was 7° (north by west 20°) (Mu et al., 2013). The physiography was mountains. Geomorphic positions were the top of slope at Site A and gentle slope at Site B, with the microtopography of turf hummocks. The dominant plant is Kobresia tibetic. This watershed is characterized by alpine subhumid climate with the mean annual precipitation (MAP) of 433 mm and the mean annual evaporation of 1080 mm (Wang et al., 2008). It should be noted that about 90% of precipitation falls between May and September in the form of rainfall, and the other 10% occurs from mid-September to mid-May as snow (Liu et al., 2005). Mean annual air temperature of Heihe River basin increased by an average of about 0.6–1 °C during the period 1950 to 2006 (from the daily data of 12 weather stations in the upper reach of Heihe River basin shown in Fig. 1) (Peng et al., 2013). Mean annual ground temperature at depth of 6 m at Site A and B were approximately -0.5 °C and -0.8 °C. Wetlands are prevalent in this study area due to the poor drainage caused by permafrost, and typical periglacial phenomena such as thaw slumping and thermokarst depression are common (Fig. 2) (Mu et al., 2013).
The two permafrost cores were collected by drilling on the north slope of Eboling Mountain from 1 to 14 March 2012. Site A was sampled to depth of 20.4 m and Site B to 11.7 m. Each core was cut into smaller segments, then photographed, wrapped, labeled, and stored in a freezer in the truck. Upon returning to the laboratory of Lanzhou University, the samples were transferred into the ultra-low temperature freezer.
LABORATORY ANALYSES
The samples were stored in the freezer at -50 °C in order to keep the samples in frozen state before processing. Then we conducted samples cutting and analyses of soil physical and chemical characteristics. The subsampling depths for Sites A and B were 600 cm and 470 cm, respectively.
Sample Cutting
All permafrost cores were cut in half lengthwise. Half of each core was cut into about 3–5 cm slices. Every fifth sample was split in duplicate for quality assurance and quality control.
Total Water Content (TWC) and Soil pH
TWC was determined by drying soils at 60 °C for about 48 h and measuring the weight before and after drying (Ahuja et al., 1985). It is expressed as the ratio of the volume of wet soil and the loss mass. The pH of the soil suspension (1:5 soil-water ratio) was measured by an acidity meter (PHS-3C, Shanghai, China).
Bulk Density (BD)
BD was determined by measuring the volume (length, width, height) of a section of frozen core, and then drying the segment at 60 °C and determining its mass.
SOC, SIC, and Elemental Analysis of C and N
SOC, C, and N of homogenized samples were quantified by dry combustion on a Vario EL elemental analyzer (Elemental, Hanau, Germany). For SOC measurement, 0.5 g dry soil samples were pretreated by HCl (10 mL 1 mol L-1) for 24 h to remove carbonate. Then the samples were rinsed with deionized water (millipore) until pH = 7. SIC contents were estimated using the following equation: SIC = Total C - SOC.
Stable Carbon Isotope Analyses
Approximately 2–3 µg of each sample was measured out into tin cups and weighted using an Electronic Precision and Analytical Balance (Sartorius, Germany). Samples were analyzed using OI Analytical Analyzer (Picarro, U.S.A.). Results are based on the mean of three replicates of each sample and expressed as δ-values relative to the Vienna Pee Dee Belemnite (VPDB) standard for δ13C. The δ-values are defined as
where Rsample and Rstandard are 13C/12C ratios of the samples and standards, respectively.CALCULATING SOC, TN, AND SIC STOCKS
The calculation of total SOC stock (SOCS) (kg m-2), soil N stock (TN) (kg m-2), and total SIC stock (SICS) (kg m-2) within soil profile was obtained using Equations 2, 3, and 4.
where SOCi = organic carbon (wt%), SICi = inorganic carbon (wt%), BDi = bulk density (g cm-3), Ni = total nitrogen concentration (wt%), hi = soil thickness (cm), and Ci = percentage of the rock fragments fraction (>2 mm) at layer i, respectively. The C/N ratios were calculated using mass ratios between SOC and TN contents.Results
SOIL BULK DENSITY, TOTAL WATER CONTENT, AND PH
The general trend for total water content increased with depth gradually with a step increase at depth of about 113 cm of Site A and the depth of about 134 cm of Site B (Figs. 3 and 4). Total water content decreased below the depth of 350 cm at Site B, where a clay layer appeared. There was a good correlation between total water content and bulk density (P < 0.001).
Soil bulk density increased with depth and differed at Sites A and B (Figs. 3 and 4). Soil bulk density of Site B was higher than that of Site A. Soil bulk density had a step increase at the depth from 60 cm to 113 cm of Site A and from 80 cm to 134 cm of Site B. Soil pH values increased with depth at Sites A and B. Soil pH value of Site A ranged from 6.08 to 7.03; the average pH value of the active layer was 6.24, and the average value of the permafrost layer was 6.67. Soil pH value of Site B ranged from 5.95 to 8.46; the average value of the active layer was 6.61, and the average value of the permafrost layer was 7.73. Overall, the pH value in permafrost was higher than that in the active layer at both sites.
SOC, SIC, TN, δ13C, AND C/N RATIOS
The SIC concentration of Site A ranged from 0.45% to 12.8%, with an average value of 5.8% (Fig. 5), while the SIC concentration ranged from 0.03% to 17.3% at Site B, with an average value of 4.7% (Fig. 6). SIC decreased with soil depth, and its distribution had the same trend for change with SOC within the soil profile. The average N concentrations were 0.89% and 0.52% at Site A and Site B, respectively. The distribution of TN with depth had the same trend for change with SOC, and the correlation coefficient was 0.95 (P < 0.001).
The δ13C values for SOM at Site A ranged from -27.3‰ to -24.9‰, and the average value was -25.8‰ (Fig. 5). The δ13C values below the permafrost table were more negative, ranging from -25.3‰ to -26.7‰, where soil water content was high. The δ13C values of the depth from 350 cm to 460 cm have wider range from -24.9‰ to -27.3‰. The δ13C values for SOM from the top layer to 350 cm at Site B ranged from -28.0‰ to -23.5‰, and the average value was -25.7‰. The δ13C values below the depth of 350 cm were larger than -10.0‰ (Fig. 6).
The soils at Site A had lower C/N ratios and the average value was 7.7 (Fig. 5). The C/N ratios have smaller changes in upper permafrost, reflecting the continuous proportional decreases with depth in the concentrations of SOC and N. The C/N ratios decreased slowly from the active layer to permafrost, then increased toward the bottom of the profile. As for Site B, the C/N ratios ranged from 3.0 to 10.0 within the soil profile (Fig. 6). The C/N ratios below 350 cm at Site B are low because there was excess of inorganic N. The C/N ratios from the depth of 80 cm to 280 cm were higher than those of the depth from 300 cm to 350 cm.
SOC, SIC, AND TN STOCKS
The SOC, TN, and SIC stocks in the active layer at Site A were 40.7 kg m-2, 4.9 kg m-2, and 10.1 kg m-2, respectively, and those in permafrost layer were 88.9 kg m-2, 12.2 kg m-2, and 23.4 kg m-2. At Site B, the SOC, TN, and SIC stocks in the active layer were 47.8 kg m-2, 5.7 kg m-2, and 14.2 kg m-2, and those in permafrost layer were 54.5 kg m-2, 3.8 kg m-2, and 45.9 kg m-2 (Fig. 7). The average stocks of SOC, TN, and SIC in permafrost layer (71.7 kg m-2, 8.0 kg m-2, 34.7 kg m-2) were probably much higher than those in the active layer (44.3 kg m-2, 5.3 kg m-2, 12.2 kg m-2) because the two permafrost cores could not represent the large area since the pedogenesis was quite complex.
Discussion
At the depth of approximately 113 cm of Site A and 134 cm of Site B, ataxitic cryostructure was well developed and total water content was near 1.0 g cm-3. It was reasonable to assume that the transient layer was at ∼113–150 cm of Site A and ∼134–160 cm of Site B, based on a relatively consistent trend of water content, bulk density, and pH values (Osterkamp et al., 2009; French and Shur, 2010; Murton and French, 1994). There is a wider range in soil bulk density and total water content in the transient layer (Figs. 3 and 4), which is caused by the migration of water in the freeze-thaw process (French, 1999). Studies showed that the transient layer is particularly sensitive to thermokarst due to the high aggradational ice contents (Shur and Jorgenson, 2007; Osterkamp et al., 2009). Thus the thaw slumping and thermokarst depression on the Eboling Mountain are mainly caused by thawing of these ice-rich layers.
The distribution of SOC, N, and SIC have the same trends for all the soil layers at Site B. This result implied that studies of soil carbon balance should include the three components or else an underestimation of the soil carbon sequestration or depletion will occur. The average SIC stocks were 0.6 times of the SOC at the two sites. It is lower than that at depth of 1 m in the Tibetan alpine grasslands, which is about 2.1 times the corresponding SOC stocks (Yang et al., 2010). One explanation is that the moisture of our study area is higher than that of most regions in the Tibetan grassland. High soil moisture can lead to lower SIC content. The extra SIC in the active layer of Site B may come from slope transport and deposits.
High C/N ratios probably reflect better preservation of SOM owing to permafrost or to wetter and anoxic conditions in the peat before the onset of permafrost (Andersson et al., 2012). Soil pH values and C/N ratios had negative correlation in permafrost (P < 0.001), which implies alkaline soil environment is favorable to SOM decomposition. The negative relationship between SIC and soil pH (P < 0.001) is because SIC can regulate soil acidity (Bronick and Lal, 2005).
Our results support the suggestion by Skrzypek et al. (2010) that δ13C values of SOM were not significantly impacted by diagenesis. The δ13C values were more positive in lower layer, which may indicate that the vascular plants existed here. Fluctuation in δ13C values was greater in the transient layer and the depth from 350 cm to 430 cm relative to the whole profile at Site A, but all δ13C values were within the range generally associated with C3 peat-forming plants (Menot and Burns, 2001). The overall trend of δ13C values displayed a little change, different with the values reported by Oksanen et al. (2001) for Rogovaya River peat plateau in the Russian Arctic, and the values in Andersson et al. (2012) for peat deposit from East European Russian Arctic.
Isotopic depth trends in modern soils that are independent of vegetation changes can be divided into three distinct categories: (1) a δ13C depth trend toward slightly lower values, (2) a uniform δ13C depth trend, or (3) a δ13C depth trend toward higher values (Alewell et al., 2011). The δ13C values of Site A belonged to the first categories, implying the soils were constantly waterlogged (Stout and Rafter, 1978). The anaerobic conditions apparently produced the small variations in δ13C values with depth and caused little microbial fractionation of original SOM (Stout et al., 1981). Total water contents, bulk densities, SOC contents, and δ13C values at Site A can indicate the syngenesis, with the intermittent soil formation due to the aggradation of permafrost caused by periodic addition of new material. The δ13C values of Site B belonged to the second categories, implying the soils were relatively young and poorly drained soils that underwent more fluctuating, less waterlogged conditions than Site A (Stout and Rafter, 1978). Thus the two sites can represent the different drainage conditions. The uniform δ13C depth profiles are mainly the result of little time for soil formation, and thus limited decomposable fractionation, which is further impeded by partial waterlogging. Pronounced δ13C increases are the norm in mature well-drained modern soils or rice land soils (Balesdent and Mariotti, 1996). In addition to aerobic conditions within the soil favoring selective loss of 12CO2 accompanying microbial degradation and fractionation (Martin and Haider, 1986), the abundant clays in subsurface were also known to selectively bind 13C (Becker-Heidmann and Scharpenseel, 1986). This is the main reason for higher δ13C values below the depth of 350 cm at Site B.
The SOC and TN stocks were greater than those of northern China's grasslands (Wu et al., 2003; Yang et al., 2010), and those in Shulehe River basin of Qinghai-Tibetan Plateau, where the average stocks of SOC and TN at depth of 1 m were about 7.72 kg m-2 and 0.93 kg m-2 (Liu et al., 2012). The SOC and SIC stocks were also greater than those in the region of the Larix gmelinii plantations in northeastern China (15.34 kg m-2, 83.38 g m-2) (Wang et al., 2012), and those in the inner Mongolian grasslands (Shi et al., 2012). Though the sampling sites in this study are limited, the two sites can show the potential for distribution and stocks of SOC, TN, and SIC in deep permafrost. This could be because climatic factors and pedogenesis are more favorable to the accumulation of soil carbon and nitrogen in this area (Baumann et al., 2009; Yang et al., 2010).
Conclusions
The transient layers in the permafrost region were at depth of ∼113–150 cm of Site A and ∼134–160 cm of Site B. The thaw slumping and thermokarst depression on the Eboling Mountain was caused by thawing of the ice-rich layers. The distribution of SIC and TN with depth had the same trend for change with SOC, therefore studies of the soil carbon balance should include three components or else an underestimation of the soil carbon sequestration/depletion will occur. Soil pH value in permafrost layer was higher than that in the active layer. The two sites have different types of drainage based on the isotopic depth trends. The lower δ13C and C/N values of the transient layer and the depth from 350 cm to 430 cm implied that the SOM was highly degraded at Site A.
The average stocks of SOC, TN, and SIC in permafrost layer above Crf layers (71.7 kg m-2, 8.0 kg m-2, 34.7 kg m-2) were higher than those in the active layer (44.3 kg m-2, 5.3 kg m-2, 12.2 kg m-2). These results highlight that SOC, TN, and SIC in permafrost regions with Alpine Kobresia meadow could make significant contribution to interaction between SOC and SIC in response to climate change in the alpine and plateau of China.
Acknowledgments
This work was supported by the National Key Scientific Research Project (grant 2013CBA01802), the National Natural Science Foundation of China (grants 91325202 and 41330634), the Open Foundations of State Key Laboratory of Cryospheric Sciences (grant SKLCS-OP-2014-08), and the State Key Laboratory of Frozen Soil Engineering (grant SKLFSE201408).