Next Article in Journal
Profiling the Instantaneous Power Consumption of Electric Machinery in Agricultural Environments: An Algebraic Approach
Previous Article in Journal
Application of Fuzzy Optimization Model Based on Entropy Weight Method in Atmospheric Quality Evaluation: A Case Study of Zhejiang Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variations and Indications of δ13CSOC and δ15NSON in Soil Profiles in Karst Critical Zone Observatory (CZO), Southwest China

1
Institute of Earth Sciences, China University of Geosciences, Beijing 100083, China
2
School of Water Resources and Environment, China University of Geosciences, Beijing 100083, China
3
Institute of Surface-Earth System Science, Tianjin University, Tianjin 300072, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(7), 2144; https://0-doi-org.brum.beds.ac.uk/10.3390/su11072144
Submission received: 16 March 2019 / Revised: 7 April 2019 / Accepted: 9 April 2019 / Published: 10 April 2019

Abstract

:
Soil carbon and nitrogen storage and stabilization are the key to solving the problems of mitigation of global warming and maintaining of crop productivity. In this study, the contents of soil organic carbon (SOC) and soil organic nitrogen (SON) and their stable isotope compositions (δ13CSOC and δ15NSON) in soil profiles were determined in two agricultural lands (including a farmland and an abandoned farmland) and four non-agricultural lands (including two shrub-grass lands and two shrub lands) in the karst critical zone observatory (CZO), Southwest China. The contents of SOC and SON were used for research on the effects of land use on SOC and SON storage, and the change of δ13CSOC and δ15NSON values in soil profiles were used to indicate SOC and SON stabilization. The results showed that agricultural activities reduced SOC and SON storage in the whole soil layers of farmland compared to non-agricultural lands, and farmland abandonment slightly increased SOC and SON storage. Crop rotation between peanut (C3) and corn (C4) affected the δ13CSOC in surface soils of agricultural lands (−21.6‰), which were intermediate between shrub lands (−22.7‰) and shrub-grass lands (−19.6‰). 15N-depleted SON in surface soils in farmland compared to those soil in other lands possibly associated with synthetic N fertilizer application. In soil layers below 30 cm depth the δ13CSOC deceased with depth, while the δ15NSON displayed irregular fluctuation. The change in δ13CSOC and δ15NSON through soil profiles in karst soils were more intensive than those in semiarid grassland soils indicating the less stabilization of SOC and SON in karst soils.

1. Introduction

Global population growth in the 21st century leads to an increasing pressure on agricultural productivity, such as food and fiber [1]. Soil organic carbon (SOC) and soil organic nitrogen (SON) contents are important indexes of soil fertility controlling agroecosystem productivity [2]. Previous studies have reported that agricultural activities, such as application of organic and inorganic fertilizers, irrigation, fallow, and intercultivation, affect SOC and SON storage, greenhouse gas emissions, loss of inorganic N, and soil structure [3,4,5,6,7]. However, these agricultural activities cause different kinds of feedback on SOC and SON dynamic due to the differences in climate, topography, and soil parent material in some studies [8,9,10]. Hence, it is necessary to further understand the soil C and N cycles associated with mitigation of global warming and maintaining of crop productivity [2].
Stable carbon isotope signature of soil is widely used to indicate the sources and turnover of SOC in agricultural ecosystem where there is a shift between C3 and C4 crop [11,12,13,14,15]. The 13C compositions in surface soils inherit that of vegetation, and there were marked δ13C discrepancies between C3 (from −20‰ to −33‰) [16] and C4 (from −17‰ to −9‰) vegetation [17]. However, the 13C fractionation by soil microorganisms in the soil organic matter (SOM) decomposition process confuses the contribution to SOM from C3 or C4 plants, and further disturbs accurate indication of SOC sources and estimation of its turnover [18]. Since the vertical change of δ13CSOC with soil, depth is widely used for the indication of the rate of SOM decomposition [19]; the research of variation of δ13CSOC in soil profile is beneficial for more accurate interpretation of δ13CSOC as an indicator of C sources based on further understanding the 13C fractionation in SOM decomposition process.
Synthetic N fertilizer can be fixed into SON through microbial immobilization and crop uptake, transformation and returning into soil as litter [20]. Therefore, the δ15NSON in cropland commonly decreases with the application of 15N-depleted inorganic N fertilizers [21]. However, overuse of N fertilizer causes many forms of N loss, such us leaching of NO3 derived from nitrification, releasing of N2 and nitrogen oxides (NOX) after denitrification, and ammonia (NH3) volatilization, and these N losses generally lead to 15N enrichment in the remaining SON [21]. Thus, the δ15NSON is generally used as a coarse indicator of the N sources and loss processes in the agroecosystem [22,23]. Improved understanding of soil N processes is important for the coordination between fertilization availability and environmental influence in agricultural lands [24].
Soils are the most key part of Earth’s critical zone, associated intimately with the sustainable development of mankind [25]. In the karst critical zone observatory (CZO) which is located in Puding county, Southwest China, we investigate the karst soil production and erosion processes, and the integrated geophysical–geochemical–ecological responses to anthropogenic perturbations and global climate change [26]. The karst soils are characterized by small capacity, uneven distribution, serious soil loss, and depletion, and these strongly threaten agricultural production [27]. Meanwhile, intensive agricultural activities lead to many environmental problems, such as karst rocky desertification and loss of water and soil [27]. SOC and SON play important roles in the maintenance of soil fertility and soil structure; therefore, the research of their biogeochemical cycling can provide guidance to solve the coordination problems between soil productivity and environmental sustainability in the karst CZO. The SOC and SON distribution, δ13CSOC and δ15NSON in soil profiles, were measured in agricultural region, karst CZO. The objectives of this study are: (1) To estimate the SOC and SON storage in agricultural lands and non-agricultural lands; (2) to identify the sources of SOC and SON in agricultural lands by using the SOC and SON contents and their stable isotope compositions; and (3) to illustrate the fractionation of C and N isotope and SOC and SON stabilization in karst soil according to the δ13CSOC and δ15NSON values in soil profiles. This research can perhaps provide fundamental information of soil nutrient element (C and N) dynamics to support the maintenance of agroecosystem productivity and protection of fragile soil resources in the karst CZO, Southwest China.

2. Materials and Methods

2.1. Study Area

The study area (Figure 1) is located in a small karst watershed (26°15.5′ N, 105°46.7′ E) which is one of the karst CZO sites, with an area of 1.54 km2, Guizhou province, Southwest China. The land in the study area underwent serious soil degradation due to intensive tillage, and then the sloped arable lands were abandoned and recovered under the “Grain for Green Programme” (GGP) [26]. This study area is beneficial for research on SOC and SON dynamics response to land management. This region has a sub-tropical monsoonal climate; the average annual air temperature is 15.1 °C; the mean annual precipitation is 1315 mm, more than 80% of which occurs during the wet season (April to September) [28]. The altitude of this region ranges from 1310 m to 1524 m. The small karst watershed is surrounded by hills, the slope of which is generally more than 30° (~60%) [28]. Quaternary deposits are mainly located in the center of the depression in the watershed. Limestones are widely distributed on the hillsides, and many of them are exposed on the surface [28]. The calcareous soils developed from limestone are classified as Mollic Inceptisols according to United States Department of Agriculture (USDA) soil taxonomy [29]. The spatial distribution of soils is significantly discontinuous, with the thickness of soil layers ranging from 10 cm to 160 cm (average thickness is 30 cm) [30]. Ploughing using the moldboard plow in the 0–30 cm soil layer depths, crop rotation between peanut (C3) and corn (C4), fertilizer (urea and compound fertilizer) and pig manure application, and straw non-recycling are managed in farmland. Most of the farmlands are located in the low-lying center of the watershed, where the deposition region of erosive soil is derived from that of hillsides, thus the thickness of soil layers in farmlands is generally over 70 cm, which is enough for ploughing by moldboard plow. Paddy land accounted for 14.39% of total watershed area, 55.65% was dry land, 23.35% was shrubland, and 6.61% was secondary forest land in this study area before 2008 [31]; subsequently, many farmlands were abandoned and even evolved into shrublands or shrub-grass lands.

2.2. Soil Sampling

Agricultural lands including farmland and abandoned farmland are located in center of the depression; non-agricultural lands including shrubland and shrub-grass land are mainly located in the middle slope of hills in study area (Figure 1). Six soil profiles were classified as farmland (FL), abandoned farmland (AFL), shrub-grass land (SGL1, SGL2), and shrubland (SL1, SL2). Three soil profiles of FL, SGL1, and SL2 were located in the north of the watershed; the other three soil profiles of AFL, SGL2, and SL1 were located in the middle of the watershed. The pictures of the six soil profiles are shown in Figure 2, and visible characteristics of these soil profiles are described carefully in Table 1. All the soil profiles were chosen with a thickness over 70 cm, in order to compare the changes of SOC and SON contents and their stable isotope compositions with depth. Generally, SOM content in the soil layers profile is significantly changed in the soil layers of 0–30 cm soil depth, whereas it is slightly varied below 30 cm soil depth [30]. Thus, soil samples were collected with 10 cm interval in the soil layers 0–30 cm deep, and with 20 cm intervals in the soil layers below 30 cm deep in July, 2016. Descriptions about six soil profiles, thickness, dominant vegetations (or main crops), and δ13C values of their leaf are given in Table 1.

2.3. Soil Analysis

Soil samples were dried with air (25 °C) and coarse roots and stones were removed. Subsamples separated by quartering were ground and passed through a 100-mesh sieve (<150 μm) for analysis of chemical properties. Soil pH (soil:water = 1:2.5) was measured using a pH meter [35], with a precision of ±0.05. Soil particle distribution was measured by a particle size analyzer (Malvern, Mastersizer 2000, England) and the results were shown as proportion of equivalent volume with a precision of ±1%. Carbonates in the powder samples (<150 μm) were removed using 0.5 mol L−1 diluted hydrochloric acid (HCl) for 24 h [36]. Inorganic N (including NH4+-N and NO3-N) were removed using 2 mol L−1 potassium chloride (KCl) for 24 h [37]. The treated samples were washed using deionized water until the supernatant liquid was neutral, dried at 55 °C until a constant weigh, ground, and passed through a 100-mesh sieve, used for analysis of the SOC and SON contents and their stable isotope ratios. The SOC and SON contents were calibrated due to loss of carbonate and inorganic N.
The SOC and SON contents were analyzed using an elemental analyzer (Elementar, Vario TOC cube, Germany) with a precision of C ± 0.1% and N ± 0.02%, monitored with standard samples (low organic content soil OAS, CatNo B2152, C: 1.55% ± 0.04%; N: 0.13% ± 0.01%) in the Laboratory of Surficial Environment Geochemistry, China University of Geosciences (Beijing).
The stable carbon isotope ratio (13C/12C) of SOC and stable nitrogen isotope ratio (15N/14N) of SON were measured utilizing an isotope mass spectrometer (Thermo, MAT-253, USA) in the Center Laboratory for Physical and Chemical Analysis, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. Measurements were normalized based on the measured values of standards material (Urea and L-glutamic acid, GBW04494, δ13CVPDB: 45.6‰ ± 0.08‰; δ15NAir: 0.24‰ ± 0.13‰), and expressed on standard δ (δ13C and δ15N) notation (‰) relative to Vienna Pee Dee Belemnite (VPDB) and air, respectively, where:
δ13C(‰) = [(Rsample − RVPDB)/RVPDB] × 1000, where R = 13C/12C
δ15N(‰) = [(Rsample − Rair)/Rair] × 1000, where R = 15N/14N.
Reproducibility as determined through replicate measurements was better than 0.1‰ for δ13C, and better than 0.2‰ for δ15N.

2.4. Statistical Analysis

Statistical analysis was performed by SPSS 18.0 (SPSS Inc., Chicago, IL, USA) and SigmaPlot 12.5 (Systat Software GmbH, Erkrath, Germany) software package. The δ13CSOC values in the same layer were reported as the means ± standard errors for the middle sites (AFL, SGL2, and SL1) and the northern sites (FL, SGL1, and SL2). Least significant difference (LSD) test (P < 0.05) was used to examine the significance of δ13CSOC values in the same layer between middle sites and northern sites. The δ15NSON values in the same layer were reported as the means ± standard errors for agricultural lands (FL, AFL) and non-agricultural lands (SGL1, SGL2, SL1, and SL2). Least significant difference (LSD) test (P < 0.05) was used to examine the significance of δ15NSON values in the same layer between agricultural lands and non-agricultural lands.

3. Results

3.1. Soil Particle Distribution and Soil pH

The soils in all sampling sites were a silt loamy texture, with sand size proportion of 1.6–8.6% (mean: 4.1%) in the whole particle mass; silt size proportion of 70.4–79.5% (mean: 75.7%); and clay size proportion of 16.2–24.6% (mean: 20.2%) (Table 2).
Soil pH values in all samples ranged from 6.4 to 7.7 (mean: 7.2), which showed there were neutral and alkaline soils in the study area (Table 2).

3.2. SOC and SON Content and C/N Ratio

The SOC contents in all soil layers of six profiles ranged from 4.58 g kg−1 to 43.30 g kg−1, and the SON contents ranged from 0.92 g kg−1 to 4.18 g kg−1 (Figure 3a,b). The SOC and SON contents in farmland soil were the lowest compared to non-agricultural land-use soils in the same-depth soil layer. Both the SOC and SON contents in six profiles generally decreased with increasing of soil depth. The SOC and SON contents in the most soil layers of abandoned farmland slightly increased compared to that in farmland in the same-depth soil layer, but those were still lower than in shrublands and shrub-grass lands.
The C/N ratios of SOM in the 0–10-cm-deep soil layer of six profiles varied from 8.22 to 10.52 (Figure 3c), and most of them obviously decreased with increases in soil depth in soil profiles. The C/N ratios in all soil layers of farmland were lower than those soils in other land uses in the same depth, especially the significantly low value (5.00–5.60) in the 10–50 cm soil layer depths.

3.3. δ13CSOC and δ15NSON in Soil Profiles

In the 0–10 cm soil layer depth, the δ13CSOC was near −21.6‰ in farmland (FL) and abandoned farmland (AFL) soils, near −19.6‰ in shrub-grass land (SGL1 and SGL2) soils, and near −22.7‰ in shrubland (SL1 and SL2) soils (Figure 4a). In the 10–30 cm soil layer depth, the δ13CSOC focused on −22.0‰ in soil profiles from northern sites (FL, SGL1, and SL2), while close to −19.1‰ in soil profiles from middle sites (AFL, SGL2, and SL1). The δ13CSOC values in six profiles all decreased by 1.1–2.4‰ from the 30 cm depth to the bottom of the soil profiles. In the 10–90 cm soil layer depth, the mean δ13CSOC in soil profiles from middle sites was 2.6–3.4‰ higher than the mean δ13CSOC in soil profiles from northern sites in the same-depth soil layer (Figure 5a). Mean δ13CSOC values in soil profiles from both northern and middle sites decreased by 1.6‰ from the 10 cm depth to the bottom in the soil profiles, contrasting with increases in δ13CSOC (increased by 0.4–1.0‰) as increasing soil depth in northern semiarid grassland (the data were from Reference [38]).
In the 0–10 cm soil layer depth, the δ15NSON values from six soil profiles were slightly varied from 2.9‰ to 4.3‰ (Figure 4b). The δ15NSON in non-agricultural land remarkably increased in the 10–20 cm soil layer depth compared to that in 0–10 cm depth (increased by 2.3–3.6‰). In most soil layer depths, the δ15NSON values in farmland were lower than those soils in other lands. Especially, in the 10–20 cm soil layer depth, the δ15NSON values in agricultural lands (~3.9‰) were significantly lower than in non-agricultural lands (5.5–7.9‰). The δ15NSON through soil profiles fluctuated more strongly in karst soils (2.4–4.1‰) compared to that in northern semiarid grassland (1.0–1.6‰) (Figure 5b).

4. Discussion

4.1. Indication of Agricultural Activity, C and N Sources in Agricultural Land

The SOC content and its spatial distribution mainly depend on input of surface vegetation including litters, secretions, and residues of roots [39,40]. In the study area, straw non-recycling led to decreased input of organic matter in farmland [41]. Thus, the SOC contents in agricultural land soils were lower than in non-agricultural land in all soil layers, as shown in Figure 3a. Furthermore, intensive tillage in farmland increases the rate of SOC decomposition associated with disturbance of protected-SOC within macro-aggregates [14,42], since aggregate protection for SOC is generally benefit for SOC storage through reduction in the SOM turnover [13]. Continuous input of new C provides abundant SOM and strongly affects the δ13CSOC in the topsoil, while resistant old C mostly determines the δ13CSOC in deeper soil layers [43]. The δ13CSOC value in topsoil is inherited from the stable C isotope composition of surface vegetation [44]. There was no doubt that the δ13CSOC values in the 0–10 cm soil layer depth were reformed by present vegetation composition (or crops), in which they had an alteration of land use. In the 0–10 cm soil layer depth, the δ13CSOC in shrub-grass land (−19.6‰) was higher than in shrubland (−22.7‰), which contributed to additional input from C4 vegetation (Table 1); while in the 0–10 cm soil layer depth of agricultural land, the δ13CSOC (−21.6‰) reflected the rotation between peanut (C3) and corn (C4). Although the abandoned farmland (AFL) had been covered by C3 weed, the δ13CSOC in the 0–10 cm soil layer depth was similar to that in farmland, resulting from short abandonment time (two years) and decomposition of new C [45]. In the 10–30 cm soil layer depth, the δ13CSOC showed obvious difference between northern sites (mean: −22.0‰) and middle sites (mean: –19.1‰), indicating the δ13CSOC in subsurface soils were mainly determined by old C the since effect of present vegetation should not reach this depth in a few years. Furthermore, we speculated that the historical vegetation composition in middle sites and in northern sites had a significant difference in δ13C. Since the δ13CSOC in middle sites remained 2.9–3.4‰ higher than that in northern sites in the 10–70 cm soil layer depth (Figure 5a), we also speculated that the historical vegetation composition did not chang for a long time.
The SON content was significantly correlated with the SOC content due to primary source SOM [46]; thus, the SON distribution was also affected by straw non-recycling and tillage practices. Extraneous inorganic N, synthetic N fertilizer, and atmospheric N deposition can also be transformed into SON through uptake and transformation by plants and microbes [20]. In farmland (FL), markedly lower C/N rations in the 10–50-cm-deep soil layer likely resulted from biological immobilization of nitrogenous fertilizer. Since SOM decomposition by microbes loses C as CO2 while N is reserved, the C/N ratios gradually deceased in this process. The decreased C/N ratios with increases in soil depth indicated the accumulation of old organic matter in Figure 3c. However, significantly lower ratios might have indicated extraneous inorganic N input. The δ15N values of synthetic N fertilizers (~0‰) are lower than that of organic N derived from the litters and microbes, due to negligible fractionation in chemical fixation of N2 [21]. Choi et al. [21] reported that various synthetic N fertilizer (mean δ15N: −0.3 ± 0.2‰) were significantly more 15N-depleted than raw or composted livestock manure (mean δ15N: 7.8 ± 0.6‰). The synthetic N fertilizer in agroecosystem is the largest N source [47]. The low δ15NSON in the 0–30 cm soil layer depth, which is strongly affected by agricultural activity (Figure 4b), possibly resulted from application of 15N-depleted N fertilizer. In non-agricultural land, the pathway that plant organic debris enter the soil is of great importance for influencing the δ15NSON in the surface soils [48]. The δ15N of plants is also affected by soil δ15N, resulting in intimate feedbacks between soil N and vegetation N in the ecosystem [44]. The δ15NSON values in the 0–10 cm soil layer depth were higher in karst soil (2.9–4.3‰) than that in semiarid grassland soil (−0.1–1.6‰) (Figure 5b). Firstly, the difference in δ15NSON may have resulted from the discrepancy of foliar δ15N of vegetations in the two regions, which depend on many factors, for example precipitation, temperature, species, foliar N concentration, N availability, and degree of N2 fixation [49]. Secondly, the δ15NSON values in the two regions were affected by atmospheric deposition, for example wet deposition (δ15N of NO3 is 2‰, δ15N of NH4+ is −12‰) in Southwest China [50] and dry deposition (10–15‰) in England [51]. The climatic difference between southwestern karst region and northern semiarid region results in the different types of atmospheric N deposition.

4.2. Fractionation of C and N Isotopic in Decomposition, Transformation, and Translocation Processes

The δ13CSOC values through soil profiles commonly increase with increases in soil depth [19], resulting from the following processes: (a) The Suess effect (13C-depleted CO2 in modern atmosphere since the industrial revolution); (b) the change of environment factors, such as water and light, affect the efficiency of CO2 conservation in photosynthesis; (c) preferential utilization of 13C-depleted plant compounds and accumulation of 13C-enriched microbial biomass; (d) downward translocation of 13C-enriched dissolved organic carbon (DOC) through profiles [18]. However, the δ13CSOC can decline with increases in soil depth when (e) recalcitrant and 13C-depleted lignin, lipid, and cellulose are accumulated at depth in the decomposition process of C4 plant debris [52]. One or more of these processes are dominant in different ecosystems, resulting in various changes in δ13CSOC through soil profiles [44]. The δ13CSOC values decreased from the intermediate layer to the bottom in the soil profiles (Figure 4a) and analogues were widely reported in other regions of Southwest China [53,54]. Accumulation of 13C-depleted lignin can likely respond to declined δ13CSOC with depth in karst soils [53]. But van Bergen et al. [55] stated that lignin accounts for a significant fraction of SOM in the topsoil, while markedly decreased in the deeper soils. Furthermore, the δ13CSOC values in the soil layer below 30-cm deep decreased from −22.5‰ to −23.6‰ with increases in soil depth in northern sites, and decreased from −19.1‰ to −20.9‰ in middle sites; both of them exceeded the range of C4 vegetation (from −17‰ to −9‰). Thus, it was disputable that accumulation of 13C-depleted lignin in the decomposition process of C4 plant debris led to the decrease in δ13CSOC with depth. Remaining SOC was 13C-depleted at depth, for which 13C-enriched carbonaceous matter could left from original materials; for example, CO2 released in SOM decomposition or CH4 released in methane fermentation processes and translocated DOC. However, soil-released CO2 and CH4 were commonly 13C-depleted compared to organic substrate [56]. In this small karst watershed, the permeable soils with high permeation are convenient for DOC leaching in a well-drained watershed [57]. A more reasonable assumption, that 13C-enriched DOC was translocated from its original site into underground water, was proposed in this paper. Liu et al. [58] reported that DOC/SOC ratios increased from 0.2% to 1.1% with depth, and δ13CDOC increased by 3–6‰ with depth in some soil profiles of karst soils. Anyhow, the 13C fractionation of SOC at depth should be further researched in karst soils.
In non-agricultural lands, the largest N source for plant and microbe uptake and loss from system derive from mineralization of SON [59]. Rapid decomposition of SON leads to 15N enrichment in the remaining SOM and mycorrhizal fungi [60]. Therefore, the SON in non-agricultural lands showed an obvious 15N enrichment in the 10–20 cm soil layer depth compared to that in the 0–10 cm soil layer depth (Figure 4b). The SON content generally decreases with increases in depth with downward water translocation, while the proportion of the inorganic N in total N increases [20]. Therefore, soil N processes are commonly dominated by SON ammonification and mineralization in the surface layers, while dominated by inorganic N transformation processes, such as nitrification and denitrification, in the intermediate and bottom soil layers [20]. The leaching of 15N-depleted NO3 derived from nitrification and releasing of 15N-depleted N2 and N2O derived from denitrification result in 15N enrichment in the remaining SON, which are commonly used for explanation of abnormally high δ15NSON in the intermediate soil layer [60]. The δ15NSON in the soil layer below 30 cm deep strongly fluctuated with depth in karst watershed (Figure 4b), which likely responded to intensive nitrification and denitrification. We also assumed that soil water processes controlled the change in δ15NSON through soil profiles. On the one hand, loss of the dissolved organic nitrogen (DON) as water transported affected δ15NSON in soil profile [20]. On the other hand, soil water processes affected redox environment, which was of great importance to control the occurrence and scope of nitrification and denitrification [61]. Furthermore, in the small karst watershed, special environmental conditions, such as monsoon climate, large terrain gradient, and uneven distribution of soil, led to uneven distribution of soil water in temporal and spatial [62]. Further, the nonuniformity of soil water complicated the processes of nitrification and denitrification. For example, Hefting et al. [63] found that the N transformation process was dominated by ammonification, where the water table level was less than 10 cm, by denitrification where it was 10–30 cm, and by nitrification where it was over than 30 cm, respectively. Sogbedji et al. [64,65] found that the paths of soil N loss were dominated by denitrification and gaseous emission in waterlogged soils, while by nitrification and NO3 leaching under well-drained aerobic condition. Thus, soil water processes may control the strongly fluctuated δ15NSON through soil profiles in the karst soils.

4.3. The Impacts of the Soil Organic C and N Isotopic Fractionation on Source Identification in the Karst Soil

The variations of δ13CSOC and δ15NSON through soil profiles in karst soils compared with that in semiarid grassland soils of northern China [38], as shown in Figure 5. In the soil layers below 10 cm deep, the change of δ13CSOC in karst soils (~1.6‰) were larger than that in semiarid grassland soils (0.4%–1.0‰) and showed an opposite trend with depth. The δ13CSOC increased with depth through soil profiles in semiarid grassland soils, which was the same as that in other studies [18,19,44]. Although the 13C fractionation of SOC at depth was not clear in karst soils, the direction and magnitude of change in δ13CSOC with depth were significantly different compared to that in semiarid grassland soils. The fractionation of 13C of SOC in subsurface soils should be considered adequately, which also indicate the lower stabilization of SOC in karst soils.
The δ15NSON in the soil organic layer can indicate the relative rate of soil N cycling due to significant correlation between δ15N and mineralization and nitrification rate [66]. The δ15NSON values in the 0–10 cm soil layer depth in karst soils (2.4–4.1‰) were higher than those in semiarid grassland soils (−0.1–1.6‰) (Figure 5b), possibly indicating more rapid soil N cycling in karst soils. In non-agricultural land, the δ15NSON in the 10–20 cm soil layer depth was 2.9‰ larger than that in the 0–10 cm soil layer depth in karst soils, while increased by 0.6–0.9‰ in 10–20 cm depth layers of semiarid grassland soils (Figure 5b), showing more 15N enrichment in remaining SON and mycorrhizal fungi in karst soils. In karst soils, the more 15N-enriched SON in the 0–10 and 10–20 cm soil layer depths indicated the stronger microbial activity and SOM decomposition, attributed to intensive microbial activities in a wetter and hotter environment. The low terrain gradient and dry climate in semiarid grassland determined the more stable soil water processes. The δ15NSON values in intermediate layers of soil profiles strongly fluctuated in karst soils compared to that in semiarid grassland soil (Figure 5b), resulting from strong spatial dynamics of soil water processes and their controlled redox environment associated with nitrification and denitrification. On the other hand, intensive soil water dynamics in the karst region intensify the N losses into other ecosystems. Inorganic N migrates more easily than organic N; therefore, application of organic fertilizer is beneficial for the controlling of agricultural nitrogen pollution in the agroecosystem, especially soils.

5. Conclusions

The vertical distribution of the SOC and SON contents, δ13CSOC and δ15NSON values, from six soil profiles revealed storage and stabilization of SOC and SON under agricultural lands and non-agricultural lands in the karst CZO, Southwest China. SOC and SON storage reduced in all soil layers of agricultural lands compared to those in non-agricultural lands. Agricultural activities, such as crop rotation and fertilizer application, affected δ13CSOC and δ15NSON in the surface soils. The more intensive changes in δ13CSOC and δ15NSON through soil profiles in karst soils compared to those in semiarid grassland soils indicated the lower stabilization of SOC and SON in karst soils. This study suggests the application of organic fertilizer is of great importance to coordinate the relationship between soil productivity and environmental sustainability in the karst CZO.

Author Contributions

Conceptualization, G.H., Z.S., and M.L.; methodology, M.L.; formal analysis, M.L. and Q.Z.; writing—original draft preparation, M.L.; writing—review and editing, M.L., G.H., and Z.S.; Project administration, G.H.; Supervision, G.H.

Funding

This work was funded by National Natural Science Foundation of China (No. 41325010; 41661144029).

Acknowledgments

The authors thank Taoze Liu from Institute of Geochemistry, Chinese Academy of Sciences and Yuntao Wu from Tianjin University for invaluable support on field sampling. The authors would like to thank the reviewers for their valuable suggestions and comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qadir, M.; Noble, A.D.; Schubert, S.; Thomas, R.J.; Arslan, A. Sodicity-induced land degradation and its sustainable management: Problems and prospects. Land Degrad. Dev. 2010, 17, 661–676. [Google Scholar] [CrossRef]
  2. Lehmann, J.; Kleber, M. The contentious nature of soil organic matter. Nature 2015, 528, 60–68. [Google Scholar] [CrossRef]
  3. Anaya, C.A.; Huber-Sannwald, E. Long-term soil organic carbon and nitrogen dynamics after conversion of tropical forest to traditional sugarcane agriculture in East Mexico. Soil Till. Res. 2015, 147, 20–29. [Google Scholar] [CrossRef]
  4. Congreves, K.A.; Hooker, D.C.; Hayes, A.; Verhallen, E.A.; Eerd, L.L.V. Interaction of long-term nitrogen fertilizer application, crop rotation, and tillage system on soil carbon and nitrogen dynamics. Plant Soil 2017, 410, 1–15. [Google Scholar] [CrossRef]
  5. Rasmussen, P.E.; Collins, H.P. Long-term impacts of tillage, fertilizer, and crop residue on soil organic matter in temperate semiarid regions. Adv. Agron. 1991, 45, 93–134. [Google Scholar]
  6. Sanz-Cobena, A.; Lassaletta, L.; Garnier, J.; Smith, P. Mitigation and quantification of greenhouse gas emissions in Mediterranean cropping systems. Agric. Ecosyst. Environ. 2017, 238, 1–4. [Google Scholar] [CrossRef]
  7. Shelton, R.E.; Jacobsen, K.L.; Mcculley, R.L. Cover crops and fertilization alter nitrogen loss in organic and conventional conservation agriculture systems. Front. Plant Sci. 2017, 8, 2260. [Google Scholar] [CrossRef]
  8. Araujo, M.A.; Zinn, Y.L.; Lal, R. Soil parent material, texture and oxide contents have little effect on soil organic carbon retention in tropical highlands. Geoderma 2017, 300, 1–10. [Google Scholar] [CrossRef]
  9. Fissore, C.; Dalzell, B.J.; Berhe, A.A.; Voegtle, M.; Evans, M.; Wu, A. Influence of topography on soil organic carbon dynamics in a Southern California grassland. Catena 2017, 149, 140–149. [Google Scholar] [CrossRef]
  10. Reddy, K.R.; Patrick, W.H. Effect of alternate aerobic and anaerobic conditions on redox potential, organic matter decomposition and nitrogen loss in a flooded soil. Soil Biol. Biochem. 1975, 7, 87–94. [Google Scholar] [CrossRef]
  11. Balesdent, J.; Mariotti, A.; Boisgontier, D. Effect of tillage on soil organic carbon mineralization estimated from 13C abundance in maize fields. Eur. J. Soil Sci. 1990, 41, 587–596. [Google Scholar] [CrossRef]
  12. Biggs, T.H.; Quade, J.; Webb, R.H. δ13C values of soil organic matter in semiarid grassland with mesquite (Prosopis) encroachment in southeastern Arizona. Geoderma 2002, 110, 109–130. [Google Scholar] [CrossRef]
  13. Blagodatskaya, E.; Yuyukina, T.; Blagodatsky, S.; Kuzyakov, Y. Turnover of soil organic matter and of microbial biomass under C3–C4 vegetation change: Consideration of 13C fractionation and preferential substrate utilization. Soil Biol. Biochem. 2011, 43, 159–166. [Google Scholar] [CrossRef]
  14. Six, J.; Elliott, E.T.; Paustian, K. Soil macroaggregate turnover and microaggregate formation: A mechanism for C sequestration under no-tillage agriculture. Soil Biol. Biochem. 2000, 32, 2099–2103. [Google Scholar] [CrossRef]
  15. West, J.B.; Bowen, G.J.; Cerling, T.E.; Ehleringer, J.R. Stable isotopes as one of nature’s ecological recorders. Trends Ecol. Evol. 2006, 21, 408–414. [Google Scholar] [CrossRef]
  16. Vagen, T.G.; Walsh, M.G.; Shepherd, K.D. Stable isotopes for characterisation of trends in soil carbon following deforestation and land use change in the highlands of Madagascar. Geoderma 2006, 135, 133–139. [Google Scholar] [CrossRef]
  17. Farquhar, G.D.; Ehleringer, J.R.; Hubick, K.T. Carbon isotope discrimination and photosynthesis. Annu. Rev. Plant Biol. 1989, 40, 503–537. [Google Scholar] [CrossRef]
  18. Krull, E.S.; Skjemstad, J.O. δ13C and δ15N profiles in 14C-dated Oxisol and Vertisols as a function of soil chemistry and mineralogy. Geoderma 2003, 112, 1–29. [Google Scholar] [CrossRef]
  19. Chen, Q.Q.; Shen, C.D.; Sun, Y.M.; Peng, S.L.; Yi, W.X.; Li, Z.A.; Jiang, M.T. Spatial and temporal distribution of carbon isotopes in soil organic matter at the Dinghushan Biosphere Reserve, South China. Plant Soil 2005, 273, 115–128. [Google Scholar] [CrossRef]
  20. Zhu, Q.; Castellano, M.J.; Yang, G. Coupling soil water processes and nitrogen cycle across spatial scales: Potentials, bottlenecks and solutions. Earth-Sci. Rev. 2018, 187, 248–258. [Google Scholar] [CrossRef]
  21. Choi, W.J.; Kwak, J.H.; Lim, S.S.; Park, H.J.; Chang, S.X.; Lee, S.M.; Arshad, M.A.; Yun, S.I.; Kim, H.Y. Synthetic fertilizer and livestock manure differently affect δ15N in the agricultural landscape: A review. Agric. Ecosyst. Environ. 2017, 237, 1–15. [Google Scholar] [CrossRef]
  22. Lim, S.S.; Kwak, J.H.; Lee, K.S.; Chang, S.X.; Yoon, K.S.; Kim, H.Y.; Choi, W.J. Soil and plant nitrogen pools in paddy and upland ecosystems have contrasting δ15N. Biol. Fertil. Soils 2015, 51, 231–239. [Google Scholar] [CrossRef]
  23. Robinson, D. δ15N as an integrator of the nitrogen cycle. Trends Ecol. Evol. 2001, 16, 153–162. [Google Scholar] [CrossRef]
  24. Denk, T.R.A.; Mohn, J.; Decock, C.; Lewicka-Szczebak, D.; Harris, E.; Butterbach-Bahl, K.; Kiese, R.; Wolf, B. The nitrogen cycle: A review of isotope effects and isotope modeling approaches. Soil Biol. Biochem. 2017, 105, 121–137. [Google Scholar] [CrossRef]
  25. Mobley, M.L. Monitoring Earth’s critical zone. Science 2009, 326, 1067–1068. [Google Scholar]
  26. Quine, T.; Guo, D.; Green, S.M.; Tu, C.; Hartley, I.; Zhang, X.; Dungait, J.; Wen, X.; Song, Z.; Liu, H. Ecosystem service delivery in Karst landscapes: Anthropogenic perturbation and recovery. Acta Geochim. 2017, 36, 416–420. [Google Scholar] [CrossRef]
  27. Wang, S.J.; Liu, Q.M.; Zhang, D.F. Karst rocky desertification in southwestern China: Geomorphology, landuse, impact and rehabilitation. Land Degrad. Dev. 2004, 15, 115–121. [Google Scholar] [CrossRef]
  28. Zhao, M.; Zeng, C.; Liu, Z.H.; Wang, S.J. Effect of different land use/land cover on karst hydrogeochemistry: A paired catchment study of Chenqi and Dengzhanhe, Puding, Guizhou, SW China. J. Hydrol. 2010, 388, 121–130. [Google Scholar] [CrossRef]
  29. Soil Survey Staff. Keys to Soil Taxonomy, 12th ed.; USDA Natural Resources Conservation Service: Washington, DC, USA, 2014; pp. 161–195.
  30. Han, G.L.; Li, F.S.; Tang, Y. Organic matter impact on distribution of rare earth elements in soil under different land uses. Clean Soil Air Water 2017, 45, 1600235. [Google Scholar] [CrossRef]
  31. Yang, R.; Zhao, M.; Zeng, C.; Chen, B.; Liu, Z. Spatiotemporal variations of Soil CO2 in Chenqi, Puding, SW China: The effects of weather and LUCC. In Hydrogeological and Environmental Investigations in Karst Systems; Andreo, B., Carrasco, F., Durán, J.J., Jiménez, P., LaMoreaux, J.W., Eds.; Springer: Berlin, Germany, 2015; pp. 191–205. [Google Scholar]
  32. Piao, H.C.; Liu, Q.M.; Deng, Y.U.; Guo, J.H.; Ran, J.C. Origins of soil organic carbon with the method of natural 13C abundance in maize fields. Acta Ecol. Sin. 2001, 21, 434–439. (In Chinese) [Google Scholar]
  33. Piao, H.C.; Zhu, J.M.; Yu, D.L.; Ran, J.C. The contrlling factors of C4-grass C/N ratios and their relationships with soil organic carbon accumulation. Quat. Sci. 2004, 24, 621–629. (In Chinese) [Google Scholar]
  34. Du, X.L.; Wang, S.J.; Luo, X.Q. Effects of different soil types on the foliar δ13C values of common local plant species in karst rocky desertification area in central Guizhou Province. Environ. Sci. 2014, 35, 3587–3594. (In Chinese) [Google Scholar]
  35. Liu, G.S.; Jiang, N.H.; Zhang, L.D.; Liu, Z.L. Soil Physical and Chemical Analysis and Description of Soil Profiles; China Standard Methods Press: Beijing, China, 1996; pp. 24–25. [Google Scholar]
  36. Midwood, A.J.; Boutton, T.W. Soil carbonate decomposition by acid has little effect on δ13C of organic matter. Soil Biol. Biochem. 1998, 30, 1301–1307. [Google Scholar] [CrossRef]
  37. Meng, L.; Ding, W.; Cai, Z. Long-term application of organic manure and nitrogen fertilizer on N2O emissions, soil quality and crop production in a sandy loam soil. Soil Biol. Biochem. 2005, 37, 2037–2045. [Google Scholar] [CrossRef]
  38. Qiu, L.P.; Wei, X.R.; Ma, T.; Wei, Y.C.; Horton, R.; Zhang, X.C.; Cheng, J.M. Effects of land-use change on soil organic carbon and nitrogen in density fractions and soil δ13C and δ15N in semiarid grasslands. Plant Soil 2015, 390, 419–430. [Google Scholar] [CrossRef]
  39. Jackson, R.B.; Canadell, J.; Ehleringer, J.R.; Mooney, H.A.; Sala, O.E.; Schulze, E.D. A global analysis of root distributions for terrestrial biomes. Oecologia 1996, 108, 389–411. [Google Scholar] [CrossRef]
  40. Jobbagy, E.G.; Jackson, R.B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 2000, 10, 423–436. [Google Scholar] [CrossRef]
  41. Liu, M.; Han, G.L.; Li, Z.C.; Liu, T.Z.; Yang, X.M.; Wu, Y.T.; Song, Z.L. Effects of slope position and land use on the stability of aggregate–associated organic carbon in calcareous soils. Acta Geochim. 2017, 36, 456–461. [Google Scholar] [CrossRef]
  42. Schjønning, P.; Thomsen, I.K. Shallow tillage effects on soil properties for temperate-region hard-setting soils. Soil Till. Res. 2013, 132, 12–20. [Google Scholar] [CrossRef]
  43. Ellert, B.H.; Janzen, H.H. Long-term biogeochemical cycling in agroecosystems inferred from 13C, 14C and 15N. J. Geochem. Explor. 2016, 88, 198–201. [Google Scholar] [CrossRef]
  44. Nel, J.A.; Craine, J.M.; Cramer, M.D. Correspondence between δ13C and δ15N in soils suggests coordinated fractionation processes for soil C and N. Plant Soil 2018, 423, 1–15. [Google Scholar] [CrossRef]
  45. Wedin, D.A.; Tieszen, L.L.; Dewey, B.; Pastor, J. Carbon isotope dynamics during grass decomposition and soil organic matter formation. Ecology 1995, 76, 1383–1392. [Google Scholar] [CrossRef]
  46. Liu, W.G.; Wang, Z. Nitrogen isotopic composition of plant-soil in the Loess Plateau and its responding to environmental change. Chin. Sci. Bull. 2009, 54, 272–279. [Google Scholar] [CrossRef]
  47. Schlesinger, W.H. On the fate of anthropogenic nitrogen. Proc. Natl. Acad. Sci. USA 2009, 106, 203–208. [Google Scholar] [CrossRef]
  48. Eissfeller, V.; Beyer, F.; Valtanen, K.; Hertel, D.; Maraun, M.; Polle, A.; Scheu, S. Incorporation of plant carbon and microbial nitrogen into the rhizosphere food web of beech and ash. Soil Biol. Biochem. 2013, 62, 76–81. [Google Scholar] [CrossRef]
  49. Craine, J.M.; Elmore, A.J.; Aidar, M.P.M.; Mercedes, B.; Dawson, T.E.; Hobbie, E.A.; Ansgar, K.; Mack, M.C.; Mclauchlan, K.K.; Anders, M. Global patterns of foliar nitrogen isotopes and their relationships with climate, mycorrhizal fungi, foliar nutrient concentrations, and nitrogen availability. New Phytol. 2010, 183, 980–992. [Google Scholar] [CrossRef]
  50. Xiao, H.Y.; Liu, C.Q. Sources of nitrogen and sulfur in wet deposition at Guiyang, Southwest China. Atmos. Environ. 2002, 36, 5121–5130. [Google Scholar] [CrossRef]
  51. Yeatman, S.G.; Spokes, L.J.; Dennis, P.F.; Jickells, T.D. Comparisons of aerosol nitrogen isotopic composition at two polluted coastal sites. Atmos. Environ. 2001, 35, 1307–1320. [Google Scholar] [CrossRef]
  52. Wynn, J.G.; Harden, J.W.; Fries, T.L. Stable carbon isotope depth profiles and soil organic carbon dynamics in the lower Mississippi Basin. Geoderma 2006, 131, 89–109. [Google Scholar] [CrossRef]
  53. Zhu, S.F.; Liu, C.Q. Vertical patterns of stable carbon isotope in soils and particle-size fractions of karst areas, Southwest China. Environ. Geol. 2006, 50, 1119–1127. [Google Scholar] [CrossRef]
  54. Han, G.L.; Li, F.S.; Tang, Y. Variations in soil organic carbon contents and isotopic compositions under different land uses in a typical karst area in Southwest China. Geochem. J. 2015, 49, 63–71. [Google Scholar] [CrossRef]
  55. Van Bergen, P.F.V.; Nott, C.J.; Bull, I.D.; Poulton, P.R.; Evershed, R.P. Organic geochemical studies of soils from the Rothamsted Classical Experiments—IV. Preliminary results from a study of the effect of soil pH on organic matter decay. Org. Geochem. 1998, 29, 1779–1795. [Google Scholar] [CrossRef]
  56. Kim, J.H.; Torres, M.E.; Choi, J.; Bahk, J.J.; Park, M.H.; Hong, W.L. Inferences on gas transport based on molecular and isotopic signatures of gases at acoustic chimneys and background sites in the Ulleung Basin. Org. Geochem. 2012, 43, 26–38. [Google Scholar] [CrossRef]
  57. Bonacci, O.; Pipan, T.; Culver, D.C. A framework for karst ecohydrology. Environ. Geol. 2009, 56, 891–900. [Google Scholar] [CrossRef]
  58. Liu, T.Z.; Liu, C.Q.; Lang, Y.C. Dissolved organic carbon and its carbon isotope compositions in hill slope soils of the karst area of southwest China: Implications for carbon dynamics in limestone soil. Geochem. J. 2014, 48, 277–285. [Google Scholar] [CrossRef] [Green Version]
  59. Poffenbarger, H.J.; Sawyer, J.E.; Barker, D.W.; Olk, D.C.; Six, J.; Castellano, M.J. Legacy effects of long-term nitrogen fertilizer application on the fate of nitrogen fertilizer inputs in continuous maize. Agric. Ecosyst. Environ. 2018, 265, 544–555. [Google Scholar] [CrossRef]
  60. Hobbie, E.A.; Ouimette, A.P. Controls of nitrogen isotope patterns in soil profiles. Biogeochemistry 2009, 95, 355–371. [Google Scholar] [CrossRef]
  61. Stewart, K.J.; Grogan, P.; Coxson, D.S.; Siciliano, S.D. Topography as a key factor driving atmospheric nitrogen exchanges in arctic terrestrial ecosystems. Soil Biol. Biochem. 2014, 70, 96–112. [Google Scholar] [CrossRef]
  62. Tokumoto, I.; Heilman, J.L.; Schwinning, S.; Mcinnes, K.J.; Litvak, M.E.; Morgan, C.L.S.; Kamps, R.H. Small-scale variability in water storage and plant available water in shallow, rocky soils. Plant Soil 2014, 385, 193–204. [Google Scholar] [CrossRef]
  63. Hefting, M.; Clement, J.C.; Dowrick, D.; Cosandey, A.C.; Bernal, S.; Cimpian, C.; Tatur, A.; Burt, T.P.; Pinay, G. Water table elevation controls on soil nitrogen cycling in riparian wetlands along a European climate gradient. Biogeochemistry 2004, 67, 113–134. [Google Scholar] [CrossRef]
  64. Sogbedji, J.M.; Es, H.M.V.; Klausner, S.D.; Bouldin, D.R.; Cox, W.J. Spatial and temporal processes affecting nitrogen availability at the landscape scale. Soil Till. Res. 2001, 58, 233–244. [Google Scholar] [CrossRef]
  65. Sogbedji, J.M.; Es, H.M.V.; Yang, C.L.; Geohring, L.D.; Magdoff, F.R. Nitrate leaching and nitrogen budget as affected by maize nitrogen rate and soil type. J. Environ. Qual. 2000, 29, 1813–1820. [Google Scholar] [CrossRef]
  66. Templer, P.H.; Arthur, M.A.; Lovett, G.M.; Weathers, K.C. Plant and soil natural abundance δ15N: Indicators of relative rates of nitrogen cycling in temperate forest ecosystems. Oecologia 2007, 153, 399–406. [Google Scholar] [CrossRef]
Figure 1. Location of study area and sampling sites.
Figure 1. Location of study area and sampling sites.
Sustainability 11 02144 g001
Figure 2. The pictures of the six soil sampling profiles. FL is farmland (cultivation over 50 years); AFL is abandoned farmland (farmland have be abandoned for two years, covered by weeds); SGL is shrub-grass land (SGL1: Farmland that has been abandoned for five years, and evolved to shrub-grass land; SGL2: Shrub-grass land with non-disturbance); SL is shrubland (SL1: Farmland that has been abandoned for eight years, and evolved to shrubland; SL2: Pear orchard that has been abandoned for eight years, and evolved to shrubland).
Figure 2. The pictures of the six soil sampling profiles. FL is farmland (cultivation over 50 years); AFL is abandoned farmland (farmland have be abandoned for two years, covered by weeds); SGL is shrub-grass land (SGL1: Farmland that has been abandoned for five years, and evolved to shrub-grass land; SGL2: Shrub-grass land with non-disturbance); SL is shrubland (SL1: Farmland that has been abandoned for eight years, and evolved to shrubland; SL2: Pear orchard that has been abandoned for eight years, and evolved to shrubland).
Sustainability 11 02144 g002
Figure 3. Variation of soil organic carbon (SOC) in and soil organic nitrogen (SON) content and C/N ratio in soil profiles. (a) SOC in soil profiles; (b) SON in soil profiles and (c) C/N ratio in soil profiles. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Figure 3. Variation of soil organic carbon (SOC) in and soil organic nitrogen (SON) content and C/N ratio in soil profiles. (a) SOC in soil profiles; (b) SON in soil profiles and (c) C/N ratio in soil profiles. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Sustainability 11 02144 g003
Figure 4. Variation of δ13CSOC, δ15NSON in soil profiles. (a) δ13CSOC in soil profiles and (b) δ15NSON in soil profiles. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Figure 4. Variation of δ13CSOC, δ15NSON in soil profiles. (a) δ13CSOC in soil profiles and (b) δ15NSON in soil profiles. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Sustainability 11 02144 g004
Figure 5. Contrastive variations of δ13CSOC, δ15NSON in soil profiles between southwestern karst region and northern semiarid grassland in China. (a) Contrastive δ13CSOC in soil profiles and (b) contrastive δ15NSON in soil profiles. Different lowercase letters indicate significant differences of δ13CSOC in same soil layer between middle sites (AFL, SGL2, and SL1) and northern sites (FL, SGL1, and SL2) in (a), and significant differences of δ15NSON in same soil layer between agricultural land and non-agricultural land in (b), at P < 0.05 based on the least significant difference (LSD) test. The δ13CSOC and δ15NSON data in grassland, shrubland, and cropland in northern semiarid grassland are from Reference [38]. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Figure 5. Contrastive variations of δ13CSOC, δ15NSON in soil profiles between southwestern karst region and northern semiarid grassland in China. (a) Contrastive δ13CSOC in soil profiles and (b) contrastive δ15NSON in soil profiles. Different lowercase letters indicate significant differences of δ13CSOC in same soil layer between middle sites (AFL, SGL2, and SL1) and northern sites (FL, SGL1, and SL2) in (a), and significant differences of δ15NSON in same soil layer between agricultural land and non-agricultural land in (b), at P < 0.05 based on the least significant difference (LSD) test. The δ13CSOC and δ15NSON data in grassland, shrubland, and cropland in northern semiarid grassland are from Reference [38]. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Sustainability 11 02144 g005
Table 1. Soil profile, profile thickness, visible characteristics, dominant vegetation species, and its δ13C value in study area.
Table 1. Soil profile, profile thickness, visible characteristics, dominant vegetation species, and its δ13C value in study area.
ProfileThickness (cm)Dominant vegetation speciesδ13C of leaf (‰)Visible characteristics
Agricultural land
FL70 cmZea mays (C4) −11.4‰ a0–20 cm: brawn, block structure, tight, abundant plant roots and debris
Arachis hypogaea (C3) 20–70 cm: yellow, block structure, tight, no rootlet
AFL70 cmArtemisia carvifolia (C3)−28.4‰ c0–6 cm: brawn, block structure, tight, abundant plant roots and debris
6–70 cm: yellow, block structure, tight, no rootlet
Non-agricultural land
SGL190 cmMiscanthus floridulus (C4)−12.3‰ b0–22 cm: brawn, fine grained, loose, abundant plant roots
Rubus biflorus (C3)−28.7‰ c24–42 cm, Brawn, block structure, tight, few plant roots
42–90 cm: yellow, block structure, tight, no rootlet, merges to weathered crust below
SGL290 cmMiscanthus floridulus (C4)−12.3‰ b0–10 cm: black humus layer, fine grained, loose, abundant plant roots
Pyracantha fortuneana (C3)−28.5‰ c10–35 cm: black brawn, block structure, tight, medium amount of plant roots
35–90 cm: red brawn, block structure, tight, no rootlet
SL190 cmRubus biflorus (C3)−28.7‰ c0–15 cm: black humus layer, granular structure, loose, abundant plant roots and debris
15–40 cm: red, block structure, tight, few plant roots
40–90 cm: brawn, clayey, tight, no rootlet
SL270 cmRhamnus davurica (C3)−28.4‰ c0–10 cm: black humus layer, fine grained, loose, abundant plant roots
Rubus biflorus (C3)−28.7‰ c10–50 cm: brawn, block structure, tight, medium amount of plant roots
50–70 cm: red, clayey, tight, no rootlet
a Data are from Reference [32]; b Data are from Reference [33]; c Data are from Reference [34]. FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.
Table 2. Soil particle distribution and soil pH in different soil depth.
Table 2. Soil particle distribution and soil pH in different soil depth.
Sampling siteDepth (cm)Clay (%)Silt (%)Sand (%)pH
FL0–1020.973.85.37.5
10–2020.673.65.87.4
20–3020.071.48.67.5
30–5019.273.17.87.3
50–7018.775.26.17.3
AFL0–1023.071.85.27.2
10–2024.670.45.07.4
20–3021.474.93.87.5
30–5021.173.95.07.4
50–7023.271.75.17.4
SGL10–1021.974.23.96.8
10–2021.474.73.96.7
20–3020.276.73.16.7
30–5021.974.33.86.4
50–7020.675.33.86.4
70–9019.377.03.86.4
SGL20–1019.876.93.37.7
10–2019.177.23.77.5
20–3019.477.43.17.7
30–5017.779.52.97.5
50–7020.278.31.57.4
70–9022.576.31.37.4
SL10–1018.077.34.76.7
10–2019.776.93.46.8
20–3020.977.61.56.8
30–5020.378.21.56.8
50–7023.174.62.37.2
70–9023.674.42.07.5
SL20–1018.077.34.77.0
10–2017.677.74.77.2
20–3016.579.14.47.3
30–5016.279.34.57.3
50–7016.678.64.87.4
Maximum 24.679.58.67.7
Minimum 16.270.41.36.4
Mean 20.275.74.17.2
Note: FL is farmland, AFL is abandoned farmland, SGL is shrub-grass land, and SL is shrubland.

Share and Cite

MDPI and ACS Style

Liu, M.; Han, G.; Zhang, Q.; Song, Z. Variations and Indications of δ13CSOC and δ15NSON in Soil Profiles in Karst Critical Zone Observatory (CZO), Southwest China. Sustainability 2019, 11, 2144. https://0-doi-org.brum.beds.ac.uk/10.3390/su11072144

AMA Style

Liu M, Han G, Zhang Q, Song Z. Variations and Indications of δ13CSOC and δ15NSON in Soil Profiles in Karst Critical Zone Observatory (CZO), Southwest China. Sustainability. 2019; 11(7):2144. https://0-doi-org.brum.beds.ac.uk/10.3390/su11072144

Chicago/Turabian Style

Liu, Man, Guilin Han, Qian Zhang, and Zhaoliang Song. 2019. "Variations and Indications of δ13CSOC and δ15NSON in Soil Profiles in Karst Critical Zone Observatory (CZO), Southwest China" Sustainability 11, no. 7: 2144. https://0-doi-org.brum.beds.ac.uk/10.3390/su11072144

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop