Next Article in Journal
Spatial Distribution and Development of Sequential Sampling Plans for Diaphorina citri Kuwayama (Hemiptera: Liviidae)
Previous Article in Journal
Exogenous Application of Thiourea for Improving the Productivity and Nutritional Quality of Bread Wheat (Triticum aestivum L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Organic Carbon Changes for Croplands across China from 1991 to 2012

1
Institute of Plant Nutrition and Resources, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
2
Key Laboratory of Plant Nutrition and Fertilizer, Ministry of Agriculture and Rural Affairs/Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, China
3
Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, 960 Carling Avenue, Ottawa, ON K1A 0C6, Canada
4
Harrow Research and Development Centre, Agriculture and Agri-Food Canada, 2585 County Road 20, Harrow, ON N0R 1G0, Canada
*
Authors to whom correspondence should be addressed.
Agronomy 2021, 11(7), 1433; https://doi.org/10.3390/agronomy11071433
Submission received: 17 June 2021 / Revised: 15 July 2021 / Accepted: 16 July 2021 / Published: 19 July 2021

Abstract

:
Accurate estimates of soil organic carbon (SOC) are critical for evaluating the impacts of crop and nutrient management practices on soil sustainability and global climate change. Temporal and spatial variations in topsoil (0–0.20 m) SOC were analyzed using 43,743 soil samples in China’s croplands. The soil database in our study was collected from the International Plant Nutrition Institute (IPNI) China Program. The results showed an increasing trend in SOC density (SOCD) for both grain and cash crops from 1991 to 2012. The average SOCD increased by 16.8, 17.4, 11.8 and 8.7% in the north central, northwest, southeast and southwest regions, respectively, whereas average SOCD decreased by 1.3% for the northeast region between the 1991–2001 and 2002–2012 periods. For both grain and cash crops, the SOCD frequency distribution (%) increased in the ranges of 10–20, 20–30 and 30–40 Mg C ha−1 from the 1991–2001 to the 2002–2012, but decreased in the ranges of 0–10 and 50–60 Mg C ha−1. Additionally, SOCD increased in most major soil types across China’s cropland regions, except in phaeozems, chernozems and umbrisols, where it decreased by 8.6–18.7% mainly due to water runoff, soil erosion, and low C input. The overall SOC stock (SOCS) in China’s cropland increased by 260 Tg C (23.7 Tg C yr−1) from 1991–2001 to 2002–2012, which was partially due to the increased crop residue return, improved fertilization and adopted conservation tillage over the period. This SOC increase represents a potential offset in carbon dioxide (CO2) emissions that could help reduce the overall net CO2 emissions in China.

1. Introduction

Globally, soil organic carbon (SOC), estimated to be approximately 1500 Pg C within 1 m depth [1], is the largest organic C pool in terrestrial ecosystems, being two or three-fold greater than the C in the atmosphere or in vegetation [2,3]. As soil stores vast quantities of C, the increase in net soil C storage by even a few percent represents a substantial C sink potential [4]. Thus, enhancing SOC sequestration could be viewed as a potential approach to mitigate global warming induced by carbon dioxide (CO2) emissions [5]. Soil organic carbon plays an important role in agricultural productivity [5,6,7]. A quantitative model was developed at the global level indicating that the yield of maize and wheat increased with increase in SOC concentrations, but yield increase levelled off when SOC concentration was approximately 2.0% [6]. The increased SOC in depleted/degraded soils could have a positive effect on crop yield [7]. They further indicated that it was difficult to establish a direct cause–effect relationship between crop yield and SOC content due to the numerous factors including biotic and abiotic stresses. Soil organic carbon is also an important indicator of soil fertility and soil productivity, which could affect soil physical properties and biochemical processes [8,9], such as water infiltration ability, soil water retention, nutrient availability, and the biological activity of microorganisms [7,10]. Hence, understanding changes of SOC storage in cropland is useful to develop management strategies for increasing agricultural productivity and mitigating climate change.
Cropland SOC is affected by management practices (e.g., fertilizer, tillage) and environmental conditions (e.g., precipitation, temperature), which can lead to changes in soil physical and chemical properties as well impacting carbon cycling [11,12]. Fertilization is an important agricultural management practice that could directly or indirectly affect the SOC inputs and thereby change the availability of nutrients and soil carbon turnover [11,13]. Compared with traditional tillage, conservation tillage can help to reduce or prevent soil carbon losses and increase overall topsoil SOC stock (SOCS) [14,15,16]. Other management practices such as crop rotation, straw return and irrigation have also been shown to significantly impact SOC content [3,7,17,18]. Changes in temperatures and precipitation, wet–dry and freeze–thaw cycles affect microbial and biotic activities, which can regulate SOM decomposition and SOC accumulation [19,20]. Increases in soil temperature have been associated with higher respiration and decomposition rates, but decreases in soil moisture slow down microbial activity, which may reduce SOC decomposition [21]. Drought restricts diffusion of extracellular enzymes and soluble organic C substrates, and flooding slows oxygen diffusion to decomposition reaction sites, which all could reduce substrate availability for C decomposition [22]. In cold regions, soil C decomposition is extremely low under frozen soils in winter due to slowly diffused substrates and low activity of extracellular enzymes, but the low substrate constraints to decomposition are removed when the soil thaws in the spring [22]. Thus, management practices and environmental factors could directly or indirectly affect SOC stock by changing C input and affecting C decomposition process. Furthermore, land use changes could also be one of the important factors that affects soil carbon stock over time by influencing C input and decomposition rates [23,24].
Exploring temporal and spatial changes in SOC is importance as SOC is widely used to evaluate the impacts of past and current agricultural practices on soil quality and C cycling [17,25,26,27,28]. Globally, a number of studies on SOC change were conducted under various soil, climate and cropping systems. For example, VandenBygaart et al. [17] reported the increased rate of SOC stock ranged from 0.10 to 0.33 Mg C ha−1 yr−1 when considering the changes from conventional tillage to no-tillage, annual to perennial cropping and conversion of crop-fallow to continuous cropping in Canada using the Century model from 1990 to 2004. Minasny et al. [25] found that SOC increased at an average of 0.3 Tg C yr−1 (1970–2010) for agricultural soils in South Korea, which could be partially attributed to an increase in biomass production resulting from increased fertilization and improved field management practices. Li et al. [27] demonstrated that the total SOC decreased by 7% from 1980 to 1995 and then rose by 6.2% again by 2012 in the Midwest croplands of the United States using a modelling approach. The SOC losses were primarily attributed to continuous cropping and intensive tillage practices, whereas the increases were due to the increased level of cropland production, adoption of conservation tillage and improved crop species. Zhao et al. [28] indicated that the average SOC stock in the topsoil (0–0.20 m) increased by 0.14 Mg C ha−1 yr−1 from 1980 to 2011 in China’s croplands, which was largely attributed to increased organic inputs driven by economics and policy. Overall, the changes of global SOC stock in croplands were mainly attributed to the changes in C inputs (e.g., residue removal vs. return), fertilizer application (inorganic vs. organic fertilizer), cropping systems (monoculture vs. rotation), tillage practices (e.g., conventional vs. reduced tillage), land use change (e.g., cropland vs. forest), and climate conditions (e.g., cool vs. warm).
China has an extensive terrestrial cropland area with a diverse set of soils, cropping systems, and climatic regions, which accounts for approximately 10% of the global cultivated area [29], but it provides food for nearly 20% of the world’s population [29,30]. Thus, accurate estimation of the temporal and spatial changes of SOC in China’s croplands is essential to assessing the sustainability of current agricultural productivity, evaluating global CO2 emissions, and examining and developing strategies to address C losses. Huang and Sun [31] and Xie et al. [32], using published data, reported an overall increasing trend of SOC in topsoil between 1980 and 2000 in China. Pan et al. [33] reported a similar trend using observations from soil monitoring studies between 1985 and 2006. However, there were some limitations in the previously conducted studies that focused on estimating spatial and temporal changes in SOC for cropland soils, including (1) the changes in SOC under different cropping systems (e.g., grain and cash crops) were ignored and (2) the insufficient field measurements or large modelling uncertainty in upscaling the results from individual sites to regional scales [34,35]. Thus, it is important to include and assess this more recent research documenting the changes in SOC density (SOCD) and stock at the national and regional scales in China over the last two decades. The International Plant Nutrition Institute (IPNI) China program has initiated a national-wide nutrient management survey in China since early 1990s and has accumulated large soil testing datasets, which provide a systematic estimation of SOC changes across the country. Therefore, the objectives of this study were (1) to evaluate the changes of SOCD for both grain and cash crops in China’s croplands from 1991 to 2012; (2) to compare spatial variations of SOCD and SOCS at regional scale between the 1991–2001 and the 2002–2012; and (3) to estimate the temporal and spatial changes of SOCD under different soil types over the periods.

2. Materials and Methods

2.1. Data Source

The datasets used for analyzing SOC changes originate mainly from the International Plant Nutrition Institute (IPNI) China Program database from 1991 to 2012. The soil samplings were collected from farmers’ fields, which were conducted by the IPNI Program in coordination with the local research institutes. The sampling sites were distributed across climatic zones encompassing all major agricultural regions of China (Figure 1 and Figure S1). At each site, three or four replicates of soil samples were collected and measured separately, and then the average over the replicates was calculated as the single observation. For each replicate of soil sample, three soil subsamples (0.20 × 0.20 m square) were collected and mixed using a soil-sampling spade. Approximate 1 kg composite sample was obtained and placed into a ziplock bag for laboratory analysis. The SOC content was converted from soil organic matter (SOC = SOM × 0.58), which was extracted and determined using sodium hydroxide, sodium EDTA and methanol based on ASI method [36,37]. Bulk density was determined from undisturbed soil cores collected in 100 cm3 cylinders after drying at 105 °C for 24 h. In total, 43,743 SOC topsoil (0–0.20 m) samples were collected and analyzed to assess the temporal and spatial variations in SOC for major agricultural regions across China’s croplands. The dataset was separated into two time periods to compare the SOC changes including the 1991–2001 and the 2002–2012. There were 15,526 records (2099 sites) from the 1991–2001 and 28,217 records (3512 sites) from the 2002–2012 which included 1227 resampling sites between the two periods (Table 1). The increasing number of records in the 2002–2012 was mainly from increased sampling replicates and density in the original monitored sites. In addition, soil classification was based on World Reference Base for Soil Resources (WRB) in our study.
Five agricultural regions were grouped based on geographical locations and China’s administrative divisions to evaluate regional variation of SOC in these soils. The five regions considered in this analysis are the northeast (NE), north central (NC), northwest (NW), southeast (SE) and southwest (SW). In addition, each agricultural region was further divided into two sub-groups based on grain crops (e.g., wheat, maize, rice, soybean) and cash crops (e.g., vegetables, sunflower, rapeseed, fruit trees, cotton) [34]. The geographical samples for different agricultural regions and site information are listed in Table 1 and Table 2.

2.2. Soil Organic Carbon Density and Stock

The SOC density (SOCD) of a single layer of soil was calculated using the following equation:
SOCD = SOC × t × BD × ( 1 Rm ) × 10 1
where SOCD is expressed in Mg C ha−1; SOC is the C concentration in g C kg−1; t and BD refer to the corresponding soil sampling thickness (cm) and soil bulk density (g cm−3), respectively; and Rm is the gravimetric fraction of rock fragments (>2 mm) in soils [38,39]. In this study, the minimum value, maximum value, standard deviation and coefficient variation were used to estimate the changes of SOCD. Statistical analysis of the SOCD across geographic regions and soil types at the resampling sites between the 1991–2001 and the 2002–2012 was conducted using group T test at the 0.05 probability level (p < 0.05) in the SPSS 20.0 package. The missing data of topsoil bulk density in China’s croplands was estimated with the topsoil SOC content using the regression equation (R2 = 0.7868, n = 410) based on the paired data of BD and SOC content (Figure S2) from our study and the data from previous publication [40,41] as follows:
BD = 1.4229   exp ( 0.0062 × SOC )
where BD represents topsoil bulk density (g cm−3) and SOC is the SOC concentration in g C kg−1.
The SOC stock (SOCS) of the topsoil (0–0.20 m) was calculated as follows:
SOCS = SOCD × Area
where SOCS is expressed in Tg C. SOCD represents the SOC density in the 0–0.20 m depth, and Area represents cropland area (ha). The annual cropland area for grain and cash crops was used to calculate the SOC stock based on statistical data from 1991 to 2012 [30].

3. Results

3.1. Changes of SOCD in Croplands from 1991 to 2012

It was found that the SOCD in the topsoil (0–0.20 m) generally increased across China’s croplands from 1991 to 2012 for both grain and cash crops (Figure 2). There was a similar trend of increased SOCD when the results were separated into grain crops and cash crops, however, the rate of increase was a little greater for cash crops than for grain crops (Figure 2b,c). The estimated changes of SOCD based on resampling sites between the 1991–2001 and the 2002–2012 also showed an increasing trend in SOCD (Table S1 and Figure S3). The SOCD significantly increased by 7.7%, 9.0% and 5.9% from the 1991–2001 to the 2002–2012 for total crops (19.4 to 20.9 Mg C ha−1), grain crops (18.8 to 20.5 Mg C ha−1) and cash crops (20.3 to 21.5 Mg C ha−1), respectively (Table S1).

3.2. Temporal and Spatial Changes of SOCD and SOCS

The frequency distribution of topsoil SOCD in China’s croplands in the 1991–2001 and the 2002–2012 is shown in Figure 3. The most common occurrence of SOCD was in the 10–20 Mg C ha−1 range. For total crops, the SOCD frequency distribution in the range of 0–10 Mg C ha−1 decreased by 41.0% from the 1991–2001 to the 2002–2012, whereas the SOCD frequency in the ranges of 10–20, 20–30 and 30–40 Mg C ha−1 increased by 10.4%, 2.9% and 36.0%, respectively. Furthermore, a similar shift towards higher SOCD in the 2002–2012 relative to the 1991–2001 was observed for grain and cash crops. The frequency of SOCD occurring in the 0–10 Mg C ha−1 range was decreased by 39.3% and 41.9% from the 1991–2001 to the 2002–2012 for grain and cash crops, respectively, whereas there was an increase in the higher ranges, except for the 50–60 Mg C ha−1 range.
Estimates of the topsoil SOCD and SOCS for the entire country and for the five geographical regions were conducted to compare their variations between the 1991–2001 and 2002–2012 (Table 3 and Table 4, Figure 4). For total crops, the results indicated that the topsoil SOCD in the 2002–2012 slightly decreased, relative to the 1991–2001, in the NE region (1.2%, 0.5 Mg C ha−1), but increased in the NC (16.8%, 2.0 Mg C ha−1), NW (17.4%, 2.2 Mg C ha−1), SE (11.8%, 1.9 Mg C ha−1) and SW (8.7%, 2.0 Mg C ha−1) regions. The SOCD from resampling sites significantly increased by 16.8%, 11.8%, 12.2% and 11.8% for total crops in the NC, NW, SE and SW regions from the 1991–2002 to the 2002–2012, respectively (Table S1). The results also indicated significant increases in SOCD for both grain and cash crops in most regions, except in the NW region for grain crops and in the SE and SW regions for cash crops. Although the SOCD decreased in the NE region, there was no significant difference between the two periods.
In this study, the topsoil SOCS increased by 260 Tg C (23.7 Tg C yr−1) from the 1991–2001 to the 2002–2012 across all the sampling regions (Table 4). For total crops, the annual changes of SOCS were −0.9, 6.7, 3.9 7.4 and 6.5 Tg C yr−1 between the 1991–2001 and the 2002–2012 in the NE, NC, NW, SE and SW regions, respectively. In addition, there were significant differences in the temporal and spatial changes of SOCS between grain crops and cash crops from the 1991–2001 to the 2002–2012. For grain crops, compared with the 1991–2001, SOCS in the 2002–2012 increased by 6.0%, 17.0%, 14.5% and 3.0% in the NC, NW, SE and SW regions, respectively. However, SOCS decreased by 0.9% in the NE. The SOCS for cash crops in the 2002–2012 increased by 33.7%, 13.6%, 5.3%, 4.0% in the NC, NW, SE and SW regions, respectively, compared to those in the 1991–2001.

3.3. Changes of SOCD under Different Soil Types

At the regional scale, the magnitude of the SOCD changes varied greatly across the main soil types (Figure 5, Table S2). For total crops, compared with the 1991–2001, the SOCD in the 2002–2012 increased by 37.0 and 10.2% in anthrosols and regosols in the NE, 18.5, 28.3 and 19.1% in fluvisols, regosols and luvisols in the NC, 17.2, 17.3 and 30.4% in cambisols, fluvisols and calcisols in the NW, 23.1, 20.1 and 28.0% in cambisols, fluvisols and ferralsols in the SE, and 9.3, 28.4, 6.8 and 44.2% in ferralsols, humic acrisols, anthrosols and cambisols in the SW region, respectively. Nevertheless, the SOCD decreased by 8.6, 18.7 and 8.9% in phaeozems, umbrisols and chernozems in the NE region, 2.9% in kastanozems in the NW region and 0.9% in anthrosols in the SE region. The same trend of SOCD in grain and cash crops were observed among most of the soil types, but a few reverse trends were found. For example, the SOCD in ferralsols increased in grain crops, but decreased in cash crops in the SW. At the resampling sites, the SOCD significantly increased by 13.1%, 25.4% and 17.2% for total crops, 12.1%, 28.3% and 21.0% for grain crops, and 14.5%, 23.1% and 14.4% for cash crops in fluvisols, regosols and luvisols in the NC region from the 1991–2001 to the 2002–2012, respectively. In addition, there were significant increases in SOCD for total crops in cambisols (21.1%), fluvisols (10.7%) and ferralsols (18.3%) in the SE region and in cambisols (21.5%) in the NW region from the 1991–2001 to the 2002–2012. For grain crops, the SOCD significantly increased by 25.4%,16.0% and 26.5% in cambisols, fluvisols and ferralsols, respectively, in the SE region and 20.5% in cambisols in the SW region between the two periods (Figure S4).

4. Discussion

4.1. Changes of SOCD and SOCS at National Scale

The increasing trend in SOC, as determined in our study, which might be caused by higher rates of increased fertilizer application and crop yields with greater rates of increased crop C inputs in the 1991–2001 relative to the 2002–2012 [42] (Table 5). The overall trend of increased SOCD in China might be caused by the improved farm management practices. These included adoption of balanced fertilization, increased organic fertilizer (human and livestock waste) inputs and higher amounts of crop residues returned to soils for different geographical regions in China’s croplands between the 1991–2001 and the 2002–2012. Similar to other countries [17,43], some field management implemented by farmers in China have contributed to SOC increase, including the increase of organic fertilizer input [44,45], improvement in crop cultivars [46] and conversion of conventional tillage to minimum or no tillage [15,47]. Previous studies indicated that the SOC stock was positively correlated with crop residue C, N fertilizer and organic fertilizer inputs [18,28,48]. The SOC sequestration and CO2 emission reduction could be achieved by optimizing management practices (e.g., residue return combined with NPK fertilizer, organic manure application, and mixed application of organic and inorganic fertilizers) [49,50]. It has been reported that the SOC stock would increase by ≥ 25.0 Tg C yr−1 or 0.63% yr−1 through the implementation of the recommended management practices in China’s croplands, compensating for CO2 emissions by 1.0% [51]. In addition, Zhao et al. [28] further indicated that application of N fertilizer increased SOC at low SOC cropland, which was attributed to the increased crop dry matter production, resulting in more C inputs to soils. However, excessive chemical N fertilizer may constrain SOC sequestration through stimulating decomposition of cellulose-dominant crop residues, reduce C retention efficiency by favoring bacteria over fungi, and even limit root growth by stimulating soil acidification [28].
Our results overall fell into the range of SOCS changes (9.6–27.6 Tg C yr−1) in topsoil in China’s croplands reported in previous studies (Table S3). Variations between studies were mainly due to different data source, monitoring periods, cultivation history, management practices and cropland area [28,31,32,52,53,54]. For example, Xie et al. [32] reported that the SOC stocks in China’s croplands (156 Mha) increased by 472 Tg C with an average rate of 23.6 Tg C yr−1 from the 1980s to the 2000s based on the Second National Soil Survey and published data. Zhao et al. [28] indicated that the SOC stock in topsoil increased by 560 Tg C with a rate of 18.1 Tg C yr−1 in topsoil croplands (130 Mha) from 1980–2011. Comparatively, the SOCS increase rate was 19.9 Tg C yr−1 in this study when we convert the area to 130 Mha, which was comparable to the SOCS increase rate in previous studies (Table S3). A previous study reported that SOC in the topsoil of farmland in China (no data for Hainan, Tibet and Yunnan) increased in 53–59%, decreased in 30–31% and stabilized in 4–6% from the early 1980s to the late 1990s [31]. Yu et al. [52], in a study using the Agro-C model, revealed that SOC stock (0–0.30 m) ranged from 11.0–36.5 Tg C yr−1 and increased in approximately 81% of Chinese croplands from 1980 to 2009, which was due to the increased crop production and retention of crop residues. However, a study using the DNDC (DeNitrification-DeComposition) model simulated decreasing trends in SOC stock with annual losses of 1.6–2.0% [55,56]. This was mainly a result of underestimated crop residues being returned to the soil over time. In addition, the changes of SOCD in this study were closely related to the variations of crop yields for grain and cash crops across five regions during the same period (1991–2012) [34]. Overall, the increased SOC sequestration in China’s croplands is identified over two decades as the result of the increased organic C input with greater returned crop residue, appropriate inorganic and organic fertilizer application and improved conservation tillage [28,51]. The net increase of SOCS in China’s croplands would help contribute to offsetting CO2 emissions and be an important factor in mitigating potential climate change impacts during the 2002–2012. Furthermore, previous reports indicated that there were great potentials for C sequestration under improved agricultural management strategies in China’s croplands [39,51]. The SOC will increase asymptotically with additional C input and reach an equilibrium C level, especially for low SOC soil over a certain time, but the rate of increase of equilibrium SOC declines gradually without continuous C input indefinitely [57]. The increasing C input no longer results in additional soil carbon when SOC at steady level reaches the ceiling [39,57]. Thus, it is worth considering that the less of C budget would remain for future exploitation by improving management practices, especially for high SOC soils. These improved management practices efforts to increase SOC sequestration and to restore productivity should be further concentrated on the sites with lower SOC.

4.2. Changes of SOCD and SOCS at Regional Scale

Similar trends of the SOCD were observed for grain crops and cash crops in different regions between 1991–2001 and 2002–2012. The increased SOCD was detected across all regions in China’s croplands, except for Northeast China (approximately a 0.5 Mg C ha−1 of SOC loss), which was consistent with a previous study, which reported that the about 0.41 Mg C ha−1 decreased in the NE region [28]. The slight decrease in SOCD in the NE region was partially due to the negative impacts of the relative low straw return rate, high initial SOC, soil erosion and degradation as well as climate change (e.g., drought, temperature stress) [28,58,59]. In addition, different agricultural management practices in each region, such as the use of chemical fertilizer, manure and the retention of crop straw (Table 5), would affect the distribution of SOCD. Statistical data could indicate the changes trend of nutrient input from chemical fertilizer, manure and straw return at the regional scale. Although the resolution of nutrient input data from statistical data did not perfectly match measurement data from actual field in our study, both datasets were collected based on farmers’ croplands. Thus, regional nutrient input data was a supplement of actual field measurements and could reflect the changes of famers’ field practices over time which provided an explanation for regional SOC changes. In addition, regional climate change impacts can be beneficial or harmful to agricultural systems depending on local climatic conditions and management practices. For example, increasing precipitation contributes to an increase in SOCD in northern China from the arid to semi-humid zone due to the increase in vegetation productivity and the increase of SOCD from tropical to cold-temperate zone in the eastern part of the country, which was likely attributed to decreased mineralization caused by decreased temperature from the south to the north [41]. Moreover, the changes in the spatial pattern of SOC can be attributed partly to the effects of soil types and texture, which have significant impacts on the soil water holding capacity and the decomposition rate of SOC [60]. Furthermore, the increased SOC stock was observed in our study across most agricultural regions in China’s croplands, indicating that the existing SOC may still has a potential to increase mainly due to that the SOC stock has not reached the upper limit of the saturation in these regions [39]. However, the decreased SOC stock in the NE region was detected, possibly because the new and lower equilibriums have not been reached yet due to a relative short cultivation history after conversion from forests and grasslands with high initial SOC [61,62]. It may reflect a great challenge of how to alleviate soil degradation in this region and achieve the C sequestration and emission reduction through effectively improved field management practices. Although CO2 leakage from carbon capture and storage was not analyzed in our study, recent study has shown that radial diffusion of CO2 leakage caused farmland degradation, degraded the morphological characteristics of spring wheat, and the impacts of CO2 leakage on horizontal topsoil CO2 concentrations was 2.67 times that of the source depth [63]. The SOC saturation, reversibility and leakage risk under different SOC levels at the regional scale should be considered in the future investigations.

4.3. Changes of SOCD under Different Soil Types

In our study, the decreased SOCD in phaeozems in the NE region from the 1991–2001 to the 2002–2012 was consistent with Xie et al. [32] and Wang et al. [64], which was partially attributed to the soil erosion by intense tillage in phaeozems region. Tillage erosion has been reported as an important erosion process contributing to SOC redistribution, which could result in SOC depletion, especially in the top and lower section of slopes [65]. Reducing tillage depth and size of implements, adopting tillage direction to upslope–downslope tillage, and contouring tillage patterns should be considered as effective strategies to reduce the negative impacts on tillage erosion [65]. Huang and Sun [31] reported the significant decrease of SOC in phaeozems and umbrisols in the NE region and Yan et al. [61] indicated the decreased SOCD in phaeozems in this region mainly due to water runoff, soil erosion and low C input. In addition, Wang et al. [64] indicated that the SOC loss in phaeozems region mainly derived from the decreases of slow and resistant SOC fractions. The increase of SOCD in anthrosols was observed in the NE, being attributed to the greater dry matter production and C sequestration efficiency than dry upland [66,67]. Increases of SOCD in anthrosols were also found in the SE and SW regions for grain crops; however, the SOCD decreased in anthrosols for cash crops because of the increased C mineralization by improving soil aeration [68]. The increased SOCD in fluvisols was found in the NC region in this study, which was consistent with Huang and Sun [31] and Xie et al. [32]. High rates of straw returned for major crops to cropland [69,70], application of combined organic manure and inorganic fertilizer [71], and development of no-tillage and reduced-tillage methods [15] helped contribute to the higher overall crop production, which resulted in increased root exudates and crop C inputs to soil and ultimately to higher SOC contents. The SOCD increase in the SE and SW regions, especially for grain crops, where ferralsols are mainly distributed, was associated with its high C saturation deficit due to the low initial SOC level [72] and organic fertilizer application [73,74].

4.4. Uncertainty and Limitations

Although we compared the SOC changes under different soil and crop types over two decades using large datasets, the uncertainties in estimating regional changes of SOC could be caused by the lack of more paired resampling data over time, detailed field management practices (e.g., fertilizer, crop residue, tillage, and irrigation) and precise cultivation history. The change of land use patterns in the past decades with the rapid development of urbanization and economic could also result in the uncertainty of resampling in the same site over time. The same method/procedure was used to reduce the deviation caused by the laboratory analysis errors, but huge soil samples were tested in different cooperating institutes under the IPNI program which could lead to uncertainties. In addition, the unavailable BD values were calculated by the regression equation based on the measurements in this study and the collected data from previous publications from 1991 to 2012. However, the changes of BD values could be affected by the different land use patterns, cultivation history, vegetation, and soil types, which may lead to the underestimation or overestimation. Furthermore, there is no accurate data in the changes of cultivated land area, especially for area changes of grain crops and cash crops every year, thus the SOCD would be affected by the new increased or reduced cropland area. It is important to implement more long-term monitoring studies to improve our understanding of the impacts of management practices, climate, and soils on SOC changes in China.

5. Conclusions

Spatial and temporal variation in topsoil organic carbon dynamics was estimated based on large soil testing datasets in China’s croplands. The increased SOC sequestration was identified between the 1991–2001 and 2002–2012 periods as a result of the increased organic C input with greater returned crop residue, appropriate inorganic and organic fertilizer application and improved conservation tillage. Spatial variation of topsoil SOC differed significantly among five agricultural regions mainly due to the differences in climate conditions, land use and soil properties. Additionally, the SOC increased in most major soil types across regions in China’s croplands, but it decreased in phaeozems, chernozems and umbrisols being attributed to water runoff, soil erosion and low C input. Since SOC levels in soils are known to be highly variable spatially, the extensive array of soil monitoring stations used in this study were valuable for understanding and encapsulating the differences between climatic zones, soil types and crops grown. In regions that were identified as having decreased SOC levels and degraded soils, alternative management practices could be implemented to combat this trend in SOC loss and soil degradation. Although we compared the changes in SOC under different soil and crop types over two decades using large datasets, the lack of more paired resampling data over time, detailed field management practices (e.g., fertilizer, crop residue, tillage and irrigation) and precise cultivation history led to the uncertainties in estimating the regional changes of SOC. In order to better understand soil carbon cycling and to identify management practices that benefit C sequestration, it is important that combination of process-based models and long-term monitoring sites are more effective to estimate the response of SOC stock to different soil types and cropping systems under climate variability. Additionally, future research should focus on exploring the effects of various field management practices with potential to build C in deeper soil depths.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/agronomy11071433/s1, Figure S1: Number of soil samples from 1991 to 2012 in China’s croplands. Figure S2: Correlation of bulk density with soil organic carbon (0–0.20 m depth) in China’s croplands. Figure S3. Soil organic carbon density (SOCD) for (a) total crops, (b) grain crops and (c) cash crops for different geographical regions in China’s croplands between the 1991–2001 and the 2002–2012. NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest). Figure S4. Changes of soil organic carbon density (SOCD) for resampling sites under different soil types for (a) total crops, (b) grain crops and (c) cash crops in China’s croplands between the 1991–2001 and the 2002–2012. NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest). Table S1: Estimated changes of soil organic carbon density (SOCD) for resampling sites between the 1991–2001 and the 2002–2012 from different geographical regions in China’s croplands. Table S2: Estimated changes of soil organic carbon density (SOCD) under different soil types between the 1991–2001 and the 2002–2012 from different geographical regions in China’s croplands. Table S3: Comparisons of the estimated topsoil soil organic carbon stock (SOCS) in China’s croplands.

Author Contributions

W.H., conceptualization, methodology, formal analysis, writing—original draft preparation; P.H., conceptualization, resources, supervision, project administration and funding acquisition, writing—review and editing; R.J., conceptualization, methodology, validation, writing—original draft; J.Y. and C.F.D., methodology and writing—review and editing; W.N.S. and B.B.G., methodology, software, writing—review and editing; W.Z., conceptualization, resources and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (No. 31972515), China Agricultural Research System (CARS-09-P31), and the International Plant Nutrition Institute China Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Relevant data applicable to this research are within the paper.

Acknowledgments

We thank the cooperators from 31 provinces in China for taking soil samples and managing the field experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Batjes, N.H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 1996, 47, 151–163. [Google Scholar] [CrossRef]
  2. Ewaran, H.E.; Berg, V.D.; Reich, P. Organic carbon in soil of the world. Soil Sci. Soc. Am. J. 1993, 57, 192–194. [Google Scholar] [CrossRef] [Green Version]
  3. Lal, R. Soil carbon sequestration impacts on global climate change and food security. Science 2004, 304, 1623–1627. [Google Scholar] [CrossRef] [Green Version]
  4. Paustian, K.; Lehmann, J.; Ogle, S.; Reay, D.; Robertson, G.P.; Smith, P. Climate-smart soils. Nature 2016, 532, 49–57. [Google Scholar] [CrossRef] [Green Version]
  5. Rumpel, C.; Amiraslani, F.; Koutika, L.S.; Smith, P.; Whitehead, D.; Wollenberg, E. Put more carbon in soils to meet Paris climate pledges. Nature 2018, 564, 32–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Oldfield, E.E.; Bradford, M.A.; Wood, S.A. Global meta-analysis of the relationship between soil organic matter and crop yields. Soil 2019, 5, 5–32. [Google Scholar] [CrossRef] [Green Version]
  7. Lal, R. Soil organic matter content and crop yield. J. Soil Water Conserv. 2020, 75, 27A–32A. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, X.; Xin, X.; Zhu, A.; Yang, W.; Zhang, J.; Ding, S.; Mu, L.; Shao, L. Linking macroaggregation to soil microbial community and organic carbon accumulation under different tillage and residue managements. Soil Tillage Res. 2018, 178, 99–107. [Google Scholar] [CrossRef]
  9. Yost, J.L.; Hartemink, A.E. Chapter Four-Soil organic carbon in sandy soils: A review. Adv. Agron. 2019, 158, 217–310. [Google Scholar]
  10. Rawls, W.J.; Pachepsky, Y.A.; Ritchie, J.C.; Sobecki, T.M.; Bloodworth, H. Effect of soil organic carbon on soil water retention. Geoderma 2003, 116, 61–76. [Google Scholar] [CrossRef]
  11. Malhi, S.S.; Nyborg, M.; Solberg, E.D.; Mcconkey, B.G.; Dyck, M.; Puurveen, D. Long-term straw management and N fertilizer rate effects on quantity and quality of organic C and N and some chemical properties in two contrasting soils in Western Canada. Biol. Fertil. Soils 2011, 47, 785–800. [Google Scholar] [CrossRef]
  12. Tang, X.; Zhao, X.; Bai, Y.; Tang, Z.; Wang, W.; Zhao, Y.; Wan, H.; Xie, Z.; Shi, X.; Wu, B.; et al. Carbon pools in China’s terrestrial ecosystems: New estimates based on an intensive field survey. Proc. Natl. Acad. Sci. USA 2018, 115, 4021–4026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tian, K.; Zhao, Y.; Xu, X.; Hai, N.; Huang, B.; Deng, W. Effects of long-term fertilization and residue management on soil organic carbon changes in paddy soils of China: A meta-analysis. Agric. Ecosyst. Environ. 2015, 204, 40–50. [Google Scholar] [CrossRef]
  14. Haddaway, N.R.; Hedlund, K.; Jackson, L.E.; Kätterer, T.; Lugato, E.; Thomsen, I.K.; Jørgense, H.B.; Isberg, P.E. How does tillage intensity affect soil organic carbon? A systematic review. Environ. Evid. 2017, 6, 30. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, Y.; Lia, X.; Gregorich, E.G.; McLaughlin, N.B.; Zhang, X.; Guo, Y.; Liang, A.; Fand, R.; Sun, B. No-tillage with continuous maize cropping enhances soil aggregation and organic carbon storage in Northeast China. Geoderma 2018, 330, 204–211. [Google Scholar] [CrossRef]
  16. Ogle, S.M.; Alsaker, C.; Baldock, F.; Bernoux, M.; Breidt, F.J.; McConkey, B.; Regina, K. Vazquez-Amabile, G. Climate and soil characteristics determine where no-till management can store carbon in soils and mitigate greenhouse gas emissions. Sci. Rep. 2019, 9, 11665. [Google Scholar] [CrossRef] [PubMed]
  17. VandenBygaart, A.J.; McConkey, B.G.; Angers, D.A.; Smith, W.; de Gooijer, H.; Bentham, M.; Martin, T. Soil carbon change factors for the Canadian agriculture national greenhouse gas inventory. Can. J. Soil Sci. 2008, 88, 671–680. [Google Scholar] [CrossRef]
  18. Berhane, M.; Xu, M.; Liang, Z.; Shi, J.; Wei, G.; Tian, X. Effects of long-term straw return on soil organic carbon storage and sequestration rate in North China upland crops: A meta-analysis. Glob. Chang. Biol. 2020, 26, 2686–2701. [Google Scholar] [CrossRef]
  19. Kirschbaum, M.U.F. The temperature dependence of soil organic matter decomposition, and the effect of global warming on soil organic storage. Soil Bid Biochem. 1995, 27, 753–760. [Google Scholar] [CrossRef]
  20. Choudhury, B.U.; Fiyaz, A.R.; Mohapatra, K.P.; Ngachan, S. Impact of land uses, agrophysical variables and altitudinal gradient on soil organic carbon concentration of north-eastern Himalayan region of India. Land Degrad. Dev. 2016, 27, 1163–1174. [Google Scholar] [CrossRef]
  21. Bronick, C.J.; Lal, R. Soil structure and management: A review. Geoderma 2005, 124, 3–22. [Google Scholar] [CrossRef]
  22. Davidson, E.A.; Janssens, I.A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 2006, 440, 165–173. [Google Scholar] [CrossRef] [PubMed]
  23. Guo, L.B.; Gifford, R.M. Soil carbon stocks and land use change: A meta analysis. Glob. Chang. Biol. 2002, 8, 345–360. [Google Scholar] [CrossRef]
  24. Xia, X.; Yang, Z.; Xue, Y.; Shao, X.; Yu, T.; Hou, Q. Spatial analysis of land use change effect on soil organic carbon stocks in the eastern regions of China between 1980 and 2000. Geosci. Front. 2016, 8, 597–603. [Google Scholar] [CrossRef] [Green Version]
  25. Minasny, B.; McBratney, A.B.; Hong, S.Y.; Sulaeman, Y.; Kim, M.S.; Zhang, Y.S.; Kim, Y.H.; Han, K.H. Continuous rice cropping has been sequestering carbon in soils in Java and South Korea for the past 30 years. Glob. Biogeochem. Cycles 2012, 26, GB3027. [Google Scholar] [CrossRef]
  26. Meersmans, J.; Martin, P.M.; Lacarce, E.; De Baets, S.; Jolivet, C.; Boulonne, L.; Lehmann, S.; Saby, N.P.A.; Bispo, A.; Arrouays, D. A high resolution map of French soil organic carbon. Agron. Sustain. Dev. 2012, 32, 841–851. [Google Scholar] [CrossRef]
  27. Li, Z.; Liu, S.; Tan, Z.; Sohl, T.L.; Wu, Y. Simulating the effects of management practices on cropland soil organic carbon changes in the Temperate Prairies Ecoregion of the United States from 1980 to 2012. Ecol. Model. 2017, 365, 68–79. [Google Scholar] [CrossRef]
  28. Zhao, Y.; Wang, M.; Hu, S.; Zhang, X.; Ouyang, Z.; Zhang, G.; Huang, B.; Zhao, S.; Wu, J.; Xie, D.; et al. Economics- and policy-driven organic carbon input enhancement dominates soil organic carbon accumulation in Chinese croplands. Proc. Natl. Acad. Sci. USA 2018, 115, 4045–4050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Food and Agriculture Organization of the United Nations. FAOSTAT Statistical Database; FAO: Rome, Italy, 2020. [Google Scholar]
  30. China Statistical Yearbook; Editorial Board of National Bureau of Statistics of China; China Statistics Press: Beijing, China, 2020.
  31. Huang, Y.; Sun, W. Changes in topsoil organic carbon of croplands in mainland China over the last two decades. Chin. Sci. Bull. 2006, 51, 1785–1803. [Google Scholar] [CrossRef]
  32. Xie, Z.B.; Zhu, J.G.; Liu, G.; Cadisch, G.; Hasegawa, T.; Chen, C.M.; Sun, H.F.; Tang, H.Y.; Zeng, Q. Soil organic carbon stocks in China and changes from 1980s to 2000s. Glob. Chang. Biol. 2007, 13, 1989–2007. [Google Scholar] [CrossRef]
  33. Pan, G.; Xu, X.; Smith, O.; Pan, W.; Rattan, L. An increase in topsoil SOC stock of China’s croplands between 1985 and 2006 revealed by soil monitoring. Agric. Ecosyst. Environ. 2010, 136, 133–138. [Google Scholar] [CrossRef]
  34. He, P.; Yang, L.; Xu, X.; Zhao, S.; Chen, F.; Li, S.; Tu, S.; Jin, J.; Johnston, A.M. Temporal and spatial variation of soil available potassium in China (1990–2012). Field Crop. Res. 2015, 173, 49–56. [Google Scholar] [CrossRef]
  35. Ma, J.; He, P.; Xu, X.; He, W.; Liu, Y.; Yang, F.; Chen, F.; Li, S.; Tu, S.; Jin, J.; et al. Temporal and spatial changes in soil available phosphorus in China (1990–2012). Field Crop. Res. 2016, 192, 13–20. [Google Scholar] [CrossRef]
  36. Portch, S.; Hunter, A. A Systematic Approach to Soil Fertility Evaluation and Improvement; Special Publ. No. 5; PPIC China Program: Hong Kong, China, 2002. [Google Scholar]
  37. Huang, X.; Du, L.; Hong, J.; Chen, G.; Wang, S.; Lian, Z.; Zhang, G.; Juang, L.; Zhang, L.; Ye, L. Study on the correlation between potassium dichromate external heating method and ASI for soil organic matter determination. Hubei Agric. Sci. 2002, 59, 122–125. [Google Scholar]
  38. Sun, W.; Huang, Y.; Zhang, W.; Yu, Y. Carbon sequestration and its potential in agricultural soils of China. Glob. Biogeochem. Cycles 2010, 24, GB3001. [Google Scholar] [CrossRef]
  39. Qin, Z.; Huang, Y.; Zhuang, Q. Soil organic carbon sequestration potential of cropland in China. Glob. Biogeochem. Cycles 2013, 27, 1–12. [Google Scholar] [CrossRef]
  40. Song, G.H.; Li, L.Q.; Pan, G. Topsoil organic carbon storage of China and its loss by cultivation. Biogeochemistry 2005, 74, 47–62. [Google Scholar] [CrossRef]
  41. Yang, Y.; Mohammat, A.; Feng, J.; Zhou, R.; Fang, J. Storage, patterns and environmental controls of soil organic carbon in China. Biogechemistry 2007, 84, 131–141. [Google Scholar] [CrossRef]
  42. He, W.; Jiang, R.; He, P.; Yang, J.Y.; Zhou, W.; Ma, J.; Liu, Y. Estimating soil nitrogen balance at regional scale in China’s croplands from 1984 to 2014. Agric. Syst. 2018, 167, 125–135. [Google Scholar] [CrossRef]
  43. Gubler, A.; Wächter, D.; Schwab, P.; Müller, M.; Keller, A. Twenty-five years of observations of soil organic carbon in Swiss croplands showing stability overall but with some divergent trends. Environ. Monit. Assess. 2019, 191, 277. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, C.A.; Zhou, L.M. Soil organic carbon sequestration and fertility response to newly-built terraces with organic manure and mineral fertilizer in a semi-arid environment. Soil Tillage Res. 2017, 172, 39–47. [Google Scholar] [CrossRef]
  45. Barreto, M.S.C.; Ramlogan, M.; Oliveir, D.M.S.; Verburg, E.E.J.; Elzinga, E.J.; Rouff, A.A.; Jemo, M.; Alleoni, L.R.F. Thermal stability of soil organic carbon after long-term manure application across land uses and tillage systems in an oxisol. CATENA 2021, 200, 105164. [Google Scholar] [CrossRef]
  46. Yin, X.; Huang, M.; Zou, Y. Changes in rice yield stability in southern China from 1949 to 2015. Agric. Environ. Lett. 2018, 3, 170038. [Google Scholar] [CrossRef] [Green Version]
  47. Wang, H.; Wang, S.; Yu, Q.; Zhang, Y.; Wang, R.; Li, J.; Wang, X. No tillage increases soil organic carbon storage and decreases carbon dioxide emission in the crop residue returned farming system. J. Environ. Manag. 2020, 261, 110261. [Google Scholar] [CrossRef]
  48. Oladele, S.O.; Adetunji, A.T. Agro-residue biochar and N fertilizer addition mitigates CO2-C emission and stabilized soil organic carbon pools in a rain-fed agricultural cropland. Int. Soil Water Conserv. Res. 2021, 9, 76–86. [Google Scholar] [CrossRef]
  49. Liu, B.; Wang, X.; Ma, L.; Chadwick, D.; Chen, X. Combined applications of organic and synthetic nitrogen fertilizers for improving crop yield and reducing reactive nitrogen losses from China’s vegetable systems: A meta-analysis. Environ. Pollut. 2021, 269, 116143. [Google Scholar] [CrossRef]
  50. Zhang, P.; Xu, S.; Zhang, G.; Pu, X.; Wang, J.; Zhang, W. Carbon cycle in response to residue management and fertilizer application in a cotton field in arid Northwest China. J. Integr. Agric. 2019, 18, 1103–1119. [Google Scholar] [CrossRef]
  51. Tao, F.; Palosuo, T.; Valkama, E.; Mäkipää, R. Cropland soils in China have a large potential for carbon sequestration based on literature survey. Soil Tillage Res. 2019, 186, 70–78. [Google Scholar] [CrossRef]
  52. Yu, Y.; Huang, Y.; Zhang, W. Modeling soil organic carbon in croplands of China, 1980-2009. Glob. Planet. Chang. 2012, 82–83, 115–128. [Google Scholar] [CrossRef]
  53. Liao, Q.; Zhang, X.H.; Li, Z.P.; Pan, G.; Smith, P.; Jin, Y.; Wu, X. Increase in soil organic carbon stock over the last two decades in China’s Jiangsu Province. Glob. Chang. Biol. 2009, 15, 861–875. [Google Scholar] [CrossRef]
  54. Lal, R. Offsetting China’s CO2 emissions by soil carbon sequestration. Clim. Chang. 2004, 65, 263–275. [Google Scholar] [CrossRef]
  55. Li, C.S.; Zhuang, Y.H.; Frolking, S.; Galloway, J.; Harriss, R.; Moore, I.I.I.B.; Schimel, D.; Wang, X. Modeling soil organic carbon change in croplands of China. Ecol. Appl. 2003, 13, 327–336. [Google Scholar] [CrossRef] [Green Version]
  56. Tang, H.J.; Qiu, J.H.; Ranst, E.V. Estimations of soil organic carbon storage in cropland of China based on DNDC model. Geoderma 2006, 134, 200–206. [Google Scholar] [CrossRef]
  57. Qin, Z.; Huang, Y. Quantification of soil organic carbon sequestration potential in cropland: A model approach. Sci. China Life Sci. 2010, 53, 868–884. [Google Scholar] [CrossRef] [PubMed]
  58. Zhang, Y.; Wu, Y.; Liu, B.; Zheng, Q.; Yin, J. Characteristics and factors controlling the development of ephemeral gullies in cultivated catchments of black soil region, Northeast China. Soil Tillage Res. 2007, 96, 28–41. [Google Scholar] [CrossRef]
  59. Ou, Y.; Rousseau, A.N.; Wang, L.; Yan, B. Spatio-temporal patterns of soil organic carbon and pH in relation to environmental factors–A case study of the black soil region of Northeastern China. Agric. Ecosyst. Environ. 2017, 245, 22–31. [Google Scholar] [CrossRef]
  60. Chaplot, V.; Bouahom, B.; Valentin, C. Soil organic carbon stocks in Laos: Spatial variations and controlling factors. Glob. Chang. Biol. 2010, 16, 1380–1393. [Google Scholar] [CrossRef]
  61. Yan, X.Y.; Cai, Z.C.; Wang, S.W.; Smith, P. Direct measurement of soil organic carbon content change in the croplands of China. Glob. Chang. Biol. 2011, 17, 1487–1496. [Google Scholar] [CrossRef]
  62. Zhou, Y.; Hartemink, A.E.; Zhou, S.; Liang, Z.; Lu, Y. Land use and climate change effects on soil organic carbon in North and Northeast China. Sci. Total Environ. 2019, 647, 1230–1238. [Google Scholar] [CrossRef]
  63. Ma, X.; Zhang, X.; Tian, D. Farmland degradation caused by radial diffusion of CO2 leakage from carbon capture and storage. J. Clean. Prod. 2020, 255, 20059. [Google Scholar] [CrossRef]
  64. Wang, X.; Yu, D.; Wang, C.; Pan, Y.; Pan, J.; Shi, X. Variations in cropland soil organic carbon fractions in the black soil region of China. Soil Tillage Res. 2018, 184, 93–99. [Google Scholar] [CrossRef]
  65. Zhao, P.; Li, S.; Wang, E.; Chen, X.; Deng, J.; Zhao, Y. Tillage erosion and its effect on spatial variations of soil organic carbon in the black soil region of China. Soil Tillage Res. 2018, 178, 72–81. [Google Scholar] [CrossRef]
  66. Liu, D.; Wang, Z.; Zhang, B.; Song, K.; Li, X.; Li, J.; Li, F.; Duan, H. Spatial distribution of soil organic carbon and analysis of related factors in croplands of the black soil region, Northeast China. Agric. Ecosyt. Environ. 2006, 113, 73–81. [Google Scholar] [CrossRef]
  67. Yan, X.; Zhou, H.; Zhu, Q.H.; Wang, X.F.; Zhang, Y.Z.; Yu, X.C.; Peng, X. Carbon sequestration efficiency in paddy soil and upland soil under long-term fertilization in southern China. Soil Tillage Res. 2013, 130, 42–51. [Google Scholar] [CrossRef]
  68. Yang, C.; Yang, L.; Ouyang, Z. Organic carbon and its fractions in paddy soil as affected by different nutrient and water regimes. Geoderma 2005, 124, 133–142. [Google Scholar] [CrossRef]
  69. Liu, X.; Li, S. Temporal and spatial distribution characteristics of crop straw nutrient resources and returning to farmland in China. TCSAE 2017, 33, 1–19. [Google Scholar]
  70. Jiang, G.; Xu, M.; He, X.; Zhang, W.; Huang, M.; Yang, X.; Liu, H.; Peng, C.; Shirato, Y.; Lizumi, T.; et al. Soil organic carbon sequestration in upland soils of northern China under variable fertilizer management and climate change scenarios. Glob. Biogeochem. Cycles 2014, 28, 319–333. [Google Scholar] [CrossRef]
  71. Xin, X.; Zhang, J.; Zhu, A.; Zhang, C. Effects of long-term (23 years) mineral fertilizer and compost application on physical properties of fluvo-aquic soil in the North China Plain. Soil Tillage Res. 2016, 156, 166–172. [Google Scholar] [CrossRef]
  72. Li, Z.T.; Li, X.G.; Li, M.; Yang, J.Y.; Turner, N.C.; Wang, X.Y.; Li, F.M. County-Scale changes in soil organic carbon of croplands in southeastern Gansu province of China from the 1980s to the Mid-2000s. Soil Sci. Soc. Am. J. 2013, 77, 2111–2121. [Google Scholar] [CrossRef]
  73. Liu, Z.; Xu, M.; Duan, Y.; Zhang, L.; Zhang, Y. Distribution of nitrogen in aggregates of black soil and red soil under long-term fertilization. J. Plant. Nutr. Fertil. 2013, 19, 1386–1392. [Google Scholar]
  74. Zhang, J.; Zhang, L.; Wang, P.; Huang, Q.; Yu, G.; Li, D.; Shen, Q.; Ran, W. The role of non-crystalline Fe in the increase of SOC after long-term organic manure application to the red soil of southern China. Eur J. Soil Sci. 2013, 64, 797–804. [Google Scholar] [CrossRef]
Figure 1. Distribution of experimental sites across five geographical regions in China’s croplands from 1991 to 2012. The green and blue dots represent grain and cash crops, respectively.
Figure 1. Distribution of experimental sites across five geographical regions in China’s croplands from 1991 to 2012. The green and blue dots represent grain and cash crops, respectively.
Agronomy 11 01433 g001
Figure 2. Trends of soil organic carbon density (SOCD) for (a1,a2) total crops, (b1,b2) grain crops and (c1,c2) cash crops in China’s croplands from 1991 to 2012.
Figure 2. Trends of soil organic carbon density (SOCD) for (a1,a2) total crops, (b1,b2) grain crops and (c1,c2) cash crops in China’s croplands from 1991 to 2012.
Agronomy 11 01433 g002
Figure 3. Frequency distribution of soil organic carbon density (SOCD) in the term of percentage of number of soil samples for (a) total crops, (b) grain crops and (c) cash crops in China’s croplands between the 1991–2001 and the 2002–2012.
Figure 3. Frequency distribution of soil organic carbon density (SOCD) in the term of percentage of number of soil samples for (a) total crops, (b) grain crops and (c) cash crops in China’s croplands between the 1991–2001 and the 2002–2012.
Agronomy 11 01433 g003
Figure 4. Geographical variation of soil organic carbon density (SOCD) for (a1,a2) total crops, (b1,b2) grain crops and (c1,c2) cash crops in China’s croplands between 1991–2001 and 2002–2012.
Figure 4. Geographical variation of soil organic carbon density (SOCD) for (a1,a2) total crops, (b1,b2) grain crops and (c1,c2) cash crops in China’s croplands between 1991–2001 and 2002–2012.
Agronomy 11 01433 g004
Figure 5. Changes of soil organic carbon density (SOCD) under different soil types for (a) total crops, (b) grain crops and (c) cash crops in China’s croplands between the 1991–2001 and the 2002–2012. NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest).
Figure 5. Changes of soil organic carbon density (SOCD) under different soil types for (a) total crops, (b) grain crops and (c) cash crops in China’s croplands between the 1991–2001 and the 2002–2012. NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest).
Agronomy 11 01433 g005
Table 1. Numbers of experimental observation across five geographical regions and two periods of the 1991–2001 and the 2002–2012 in China’s croplands.
Table 1. Numbers of experimental observation across five geographical regions and two periods of the 1991–2001 and the 2002–2012 in China’s croplands.
ItemRegion 1Total CropsGrain CropsCash Crops
1991–20012002–20121991–20012002–20121991–20012002–2012
Soil TestNE (199) 21826445713443620482837
NC (352)381770952751409710662998
NW (122)1548334460612019422143
SE (407)489897853320617415783611
SW (187)343735361851202615861510
China (1227)15,52628,217987217,118565411,099
1 Regions used in this study are NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest). 2 The values in brackets represent the number of the resampling sites between two time periods at different regions.
Table 2. Summary of experimental sites from 1991 to 2012 for different geographical regions in China’s croplands.
Table 2. Summary of experimental sites from 1991 to 2012 for different geographical regions in China’s croplands.
Region 1ProvinceMain CropMain Soil Type
(WRB Reference Soil Group)
Sample Number (n)Soil pHPrecipitation
(mm)
Latitude (°N)Longitude (°E)
NEJilin, Liaoning,Maize, rice, soybean, tomato Phaeozems, Umbrisols, Regosols,62833.7–9.5400–100037.74–53.53118.86–135.07
HeilongjiangCabbage, cucumber, fluxChernozems, Anthrosols
NCBeijing, Tianjin,Wheat, maize, cotton, peanutFluvisols, Regosols, 10,9123.4–10.0350–90031.41–42.67111.25–122.63
Hebei, Henan,Eggplant, cabbage, cucumberLuvisols, Solonchaks
Shandong, ShanxiTomato, cauliflower, pumpkin
NWShaanxi, Ningxia,Maize, wheat, potato Cambisols, Kastanozems,48925.0–9.9100–60027.23–53.3573.45–126.04
Gansu, Xinjiang,Cotton, carrot, cabbageFluvisols, Calcisols
Inner Mongolia,Spinach, pepper, cucumber
Qinghai, TibetOnion, tomato, rapeseed
SEHubei, Hunan,Wheat, maize, rice Cambisols, Fluvisols,14,6833.6–8.8700–160023.58–28.28108.38–122.20
Jiangsu, Anhui,Beans, cotton, sugarcaneFerralsols, Anthrosols
Shanghai, Jiangxi,Cabbage, citrus, banana
Zhejiang, FujianRapeseed, sesame
SWChongqing, Maize, wheat, rice, teaFerralsols, Humic Acrisols,69733.4–8.5600–200018.17–34.3097.39–117.06
Guizhou, Yunnan,Tomato, sugarcane, rapeAnthrosols, Cambisols
Guangxi, Hainan,Rapeseed, banana, cassava
Guangdong, SichuanPepper, pineapple
1 Regions used in this study are NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest).
Table 3. Estimated changes of soil organic carbon density (SOCD) between the 1991–2001 and the 2002–2012 from different geographical regions in China’s croplands.
Table 3. Estimated changes of soil organic carbon density (SOCD) between the 1991–2001 and the 2002–2012 from different geographical regions in China’s croplands.
Crop TypeItemNE 1 NC NW SE SW China
1991–20012002–20121991–20012002–20121991–20012002–20121991–20012002–20121991–20012002–20121991–20012002–2012
Total Min. 27.97.56.37.16.47.26.56.96.68.66.37.1
cropsMedian39.739.310.611.711.712.515.015.521.823.014.315.7
Mean41.140.611.713.612.714.915.817.623.425.518.620.9
Max.96.4106.033.852.269.361.147.460.268.369.496.4106.0
SD20.921.15.25.19.46.96.87.511.510.913.614.8
CV (%)51.451.242.638.963.947.140.843.348.343.773.170.8
Grain Min.7.97.56.37.16.47.26.56.96.68.66.37.1
cropsMedian35.636.310.511.211.712.214.914.921.420.414.015.5
Mean38.938.511.912.612.414.615.017.222.823.518.120.7
Max.80.1106.033.852.253.247.147.455.168.369.480.1106.0
SD19.220.15.24.69.56.27.16.910.710.312.914.9
CV (%)51.652.842.936.662.743.443.741.547.046.271.271.9
Cash Min.10.87.96.37.26.97.66.78.07.09.46.37.2
cropsMedian47.946.610.712.611.712.917.716.326.126.215.416.3
Mean49.248.911.315.113.315.117.318.325.926.919.821.2
Max.96.4104.731.848.769.361.142.260.261.465.996.4104.7
SD22.321.45.05.69.37.36.28.213.310.915.114.7
CV (%)44.142.941.540.264.748.535.344.249.040.576.169.3
1 Regions used in this study are NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest). 2 Min., minimum; Max., maximum; SD, standard deviation; CV, coefficient variation.
Table 4. Estimated changes of soil organic carbon stock (SOCS) between the 1991–2001 and the 2002–2012 from different geographical regions in China’s croplands.
Table 4. Estimated changes of soil organic carbon stock (SOCS) between the 1991–2001 and the 2002–2012 from different geographical regions in China’s croplands.
Crop TypeRegion 1Crop Area 2SOCS (Tg C)Total Changes Annual Changes
(Mha)1991–20012002–2012(Tg C)(%)(Tg C yr−1)
Total cropsNE18.5759749−9.8−1.3−0.9
NC37.743951373.916.86.7
NW19.725029343.317.43.9
SE43.869177281.611.87.4
SW35.182389471.38.76.5
Grain cropsNE16.0622617−5.7−0.9−0.5
NC27.032033919.26.01.7
NW13.917320329.417.02.7
SE28.342648861.714.55.6
SW23.553555116.03.01.5
Cash cropsNE2.5121120−0.9−0.7−0.1
NC10.712116240.833.73.7
NW5.8778710.413.60.9
SE15.526828214.15.31.3
SW11.630131312.04.01.1
TopsoilChina154.7296132212608.823.7
1 Regions used in this study are NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest). 2 The values were collected from the China Statistical Yearbook (1991–2012) [30].
Table 5. Changes in chemical and organic fertilizer inputs, straw returned to soils and crop production for different geographical regions in China’s croplands between the 1991–2001 and the 2002–2012.
Table 5. Changes in chemical and organic fertilizer inputs, straw returned to soils and crop production for different geographical regions in China’s croplands between the 1991–2001 and the 2002–2012.
Region 1Chemical Fertilizer 2
(Tg)
Organic Manure 3
(Tg)
Straw Returned 4
(Tg)
Total Crop Production 5
(Tg)
1991–20012002–20121991–20012002–20121991–20012002–20121991–20012002–2012
NE3.2 *4.7174.7254.97.313.4103.9142.3
NC10.514.9549.7672.763.8102.7274.2441.5
NW3.25.7230.2268.65.710.576.4138.7
SE12.315.3537.5701.033.545.3287.0368.7
SW7.410.4598.9747.925.534.6248.3353.0
1 Regions used in this study are NE (Northeast), NC (North Central), NW (Northwest), SE (Southeast) and SW (Southwest). 2 Chemical fertilizer indicates total N, P, and K nutrients from inorganic fertilizers. 3 Organic manure indicates total amount of human and livestock waste applied to soils. 4 Straw returned indicates the amount of crop residues returned to soils. 5 Total crop production indicates total production of grain and cash crops. * The values were collected from the China Statistical Yearbook (1991–2012) [30,42].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

He, W.; He, P.; Jiang, R.; Yang, J.; Drury, C.F.; Smith, W.N.; Grant, B.B.; Zhou, W. Soil Organic Carbon Changes for Croplands across China from 1991 to 2012. Agronomy 2021, 11, 1433. https://doi.org/10.3390/agronomy11071433

AMA Style

He W, He P, Jiang R, Yang J, Drury CF, Smith WN, Grant BB, Zhou W. Soil Organic Carbon Changes for Croplands across China from 1991 to 2012. Agronomy. 2021; 11(7):1433. https://doi.org/10.3390/agronomy11071433

Chicago/Turabian Style

He, Wentian, Ping He, Rong Jiang, Jingyi Yang, Craig F. Drury, Ward N. Smith, Brian B. Grant, and Wei Zhou. 2021. "Soil Organic Carbon Changes for Croplands across China from 1991 to 2012" Agronomy 11, no. 7: 1433. https://doi.org/10.3390/agronomy11071433

APA Style

He, W., He, P., Jiang, R., Yang, J., Drury, C. F., Smith, W. N., Grant, B. B., & Zhou, W. (2021). Soil Organic Carbon Changes for Croplands across China from 1991 to 2012. Agronomy, 11(7), 1433. https://doi.org/10.3390/agronomy11071433

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