Fractions of soil organic and inorganic carbon and their relationships in typical loess cropland of Fenggu Basin

There is evidence of connections between soil organic carbon (SOC) and inorganic carbon (SIC) in dryland of north China. However, fractions of SOC and SIC and the relationship are not well understood in the Loess Plateau that undergoes profound erosion and redeposition. A study was conducted in low-elevation cropland of Loess Plateau across two distinctive basins: Lingfeng basin (LFB) with lower soil pH (< 8.4) and subject to erosion-redeposition, and Yuncheng basin (YCB) with higher soil pH (> 8.6) and under the inuence of the Yellow River. Soil samples were collected from 30 sites over 100 cm. We determined SOC, SIC, dissolved organic carbon (DOC) and other properties. Above 100 cm, SOC stock is signicantly higher in LFB (10.0 ± 2.6 kg C m − 2 ) than in YCB (6.9 ± 1.5 kg C m − 2 ), but SIC lower in LFB (14.0 ± 2.5 kg C m − 2 ) than in YCB (17.0 ± 5.7 kg C m − 2 ). We nd a signicantly negative correlation between SOC and SIC stocks in LFB, but no clear relationship in YCB. On average, DOC:SOC ratio is signicantly higher below 40 cm in YCB (1.9%) than LFB (1.2%), indicating stronger DOC desorption in YCB that has stronger hydrological process due to the inuence of the Yellow River. Overall, SOC has a negative correlation with SIC and soil pH, and DOC:SOC ratio has a signicantly positive correlation with soil pH. Our analyses suggest that erosion/re-deposition of topsoil is partly responsible for the negative SIC-SOC relationship in LFB, and high soil pH and stronger hydrological processes are attributable to relatively lower levels of SOC in YCB. This study highlights that soil carbon fractions in the lowland of Loess Plateau are inuenced by many drivers, which leads to complex relationships between major soil carbon pools. There is evidence that precipitation has large impacts on both SOC and SIC in arid/semi-arid lands Raheb particularly in the Loess Plateau (Han et al. 2018). The YCB has a large spatial variation in both precipitation (500-750 mm) and elevation (380-820 m a.s.l.), which could have large impacts on the distributions of SOC and SIC, leading to alterations of the SIC-SOC relationship. Our analyses indicate that there is an increasing trend (from 1.0 to 5.8) in SIC:SOC ratio with increasing elevation in the YCB. In addition, the distance from the Yellow River is probably another factor that can alter the SIC-SOC relationship through the inuences of hydrological processes on either SIC (Shi et al. 2017b) or SOC (Shi et al. 2017a). We further discuss the potential inuence of hydrological processes on SOC below. over 0-100 cm in low-elevation cropland of the Loess Plateau, i.e., LFB and YCB. Our data show that SIC stock is signicantly higher than SOC stock over 0-100 cm in the whole study area; SOC stock is negatively correlated with soil pH, but positively correlated with EC. SOC and SIC stocks reveal a signicantly negative correlation in LFB, which is partly related to erosion and re-deposition/redistribution of topsoil. We nd no clear relationship between SOC and SIC stocks in YCB, which may reect the large spatial variations of elevation and precipitation. Our data show a signicantly positive correlation between DOC and SOC in both topsoil and subsoil of YCB, but only in topsoil of LFB. We nd that DOC:SOC ratio is signicantly higher below 40 cm in YCB than LFB, and DOC:SOC ratio is positively correlated with soil pH, and negatively correlated with EC. Our analyses suggest that high soil pH and stronger hydrological processes are attributable to the relatively lower levels of SOC in YCB. Further studies are needed to investigate how soil properties and environmental conditions interplay to regulate the dynamics of SOC and SIC in the Loess Plateau.


Introduction
The storage (~ 1500 Pg) of soil organic carbon (SOC) is greater than the sum of carbon stock in the atmosphere (~ 750 Pg) and terrestrial biosphere (~ 560 Pg), acting as both sources and sinks and contributing huge share to the regional and/or global carbon budgets (Amundson 2001, Jobbágy & Jackson 2000. The pool of soil inorganic carbon (SIC) is comparable to that of SOC in both the global land (Lal 2004) and China's land ). Although there are many studies that demonstrate some kinds of connection between SOC and SIC (Mehra et al. 2019, Monger et al. 2015, Zamanian et al. 2016, the relationship between SOC and SIC has not been well understood. There are some studies addressing SIC-SOC relationship, which report inconsistent ndings. While there is evidence of a positive correlation between SIC and SOC stocks in the cropland of north China (Guo et al. 2016, Shi et al. 2017b, Wang et al. 2015b, studies also report a negative SIC-SOC relationship under various land uses in the north China (Li et al. 2010, Zhao et al. 2016. In particular, Zhao et al. (2016) found that SOC had a negative relationship with SIC under mixed land uses (i.e., forest, grass, shrub lands) in the Loess Plateau that undergoes profound erosion.
The differences in SIC-SOC relationship may re ect the differences in the responses to changes of environmental conditions between SOC and SIC. For instance, there is evidence that high soil pH often leads to low levels of SOC but high levels of SIC (Oste et al. 2002, Tavakkoli et al. 2015, Wang et al. 2015a. Recent studies suggest that hydrologic process may have in uences on various fractions of soil carbon, including dissolved organic (DOC) and inorganic carbon (DIC) that result from the desorption of organic matter and the dissolution of inorganic carbon, respectively (Shi et al. 2017b, Zhang et al. 2020. There is evidence that higher levels of CaCO 3 or Ca 2+ /Mg 2+ are bene cial for SOC stabilization (Tavakkoli et al. 2015, Virto et al. 2011) because of enhanced formation of soil aggregates (Rowley et al. 2018).
There have been some studies on soil carbon dynamics in the Loess Plateau, showing large variability in both SOC and SIC, which is in uenced by many environmental factors such as climate (e.g., precipitation, temperature), and vegetation types (Han et al. 2018. The Loess Plateau has undergone severe soil erosions due to the poor structure of loess (Fu et al. 2011). Numerous studies have reported that erosion associated processes, including detachment, transport and deposition of soil materials, have large impacts on SOC distribution (Schiettecatte et al. 2008, Zhu et al. 2014).
There is a distinctive difference in vertical distribution between SOC and SIC in the north China's cropland: i.e., a sharp decrease in SOC but a general increase in SIC with depth (Shi et al. 2017b, Zhang et al. 2015. Erosion and redeposition of topsoil from highlands to lowlands could have in uences on the storages of both SOC and SIC, thus altering the SOC-SIC relationship in the Loess Plateau. The Fengu Basin, consisting of Linfeng Basin (LFB) and Yuncheng Basin (YCB), is located at the southeast edge of the Loess Plateau, which is subject to erosion-redeposition (in LFB) and under the in uence of Yellow River (In YCB). We hypothesize that the relationship between SIC and SOC is more complex in the cropland of Fenggu Basin relative to other regions of north China (e.g., the North China Plain). The objectives of this work are to study the dynamics of SOC and SIC in the low-elevation cropland of the Loess Plateau, to assess the relationship between major soil carbon fractions, and to explore the underlying mechanism responsible for the variability of SOC and SIC in semi-arid region.

Characteristics of the study area
Our study area (110°24¢12²-111°42¢29²E, 34°54¢43²-36°15¢52²N, ~4000 km 2 ) is located at the southwest of Shanxi Province on the Loess Plateau in China (Fig. 1), which spans almost the entire Fengu Basin's cropland. The area is dominated by semiarid continental monsoon climate with noticeable seasonal changes: a severely dry spring, hot and rainy summer, and a cold and dry winter. Main soil type is identi ed as Calcareous cinnamon (Huang et al. 2007, Shi et al. 2006. The region has an agricultural history of thousands of years of a double cropping system, i.e., winter wheat and summer maize rotation, with similar fertilization and tillage management. Irrigation is often applied using waters from the rivers.
Despite of the similarity in climate and agricultural practice, there are some differences in other aspects between the LFB and YCB.
LFB is a canyon basin with elevation of 389-554 m and Fenhe river owing from the north to the south, which is surrounded by mountains (> 1000 m a.s.l.). The parent material in the region is redeposited loess. Annual mean temperature is 8.9-12.9 C. Annual mean precipitation ranges from 420 at lower altitude to 550 mm at higher altitude, with 70% during June-September. Annual average evaporation is ~1660 mm.
The parent material in the YCB is the alluvial loess. The elevation varies from 380 m to 820 m. Annual mean temperature is ~13 C, and annual mean evaporation is ~1800 mm. Annual precipitation has a range of ~500-750 mm, showing an increasing trend with an increase in elevation, and 69% during the period of June-September. The majority of YCB is in uenced by the Yellow River, and the northern part is also in uenced by Fenhe river.

Soil sampling and analyses
In August 2017, during mature stage of maize, soil samples were collected from 30 sites ( Figure 1) using a soil auger to a depth of 100 cm, at 20-cm intervals. At each site, four soil cores were randomly taken, and soils for each layer were mixed. Soil bulk density was determined for the 0-20, 20-40 and 40-60 cm at several representative sites. Soils were air-dried, crushed, mixed thoroughly and passed 2-mm screen. A portion of soil samples were crushed to <0.25 mm, which were used for the measurements of total carbon (TC), SOC and total nitrogen (TN). We prepared soil-water mixtures (1:5) using 2-mm soils for measurements of soil pH, electrical conductivity (EC), and water extractable Ca 2+ and Mg 2+ (by using an Atomic Absorption Spectrophotometer). Contents of TC and TN were measured using a CNHS-O analyzer (Model EuroEA3000). SOC content was determined by K 2 Cr 2 O 7 oxidation titration (Walkley & Black 1934). Content of SIC was obtained by subtracting SOC from TC.
For the measurements of DOC and dissolved inorganic carbon (DIC), 6 g 2-mm soil was treated with 24 ml 0.05 M K 2 SO 4 solution for 4 hours at 25 C, the mixture was shaken for one hour, then centrifugation. The supernatant was ltered through a 0.45-μm membrane, followed by analyses of DOC and DIC by a TOC analyzer (TOC-VCPH, Shimadzu).

Data calculation and Statistical analyses
Stocks of SOC, SIC and DOC were calculated using the following equations: where X i is carbon content for layer i, D i soil layer thickness (cm), B i bulk density (g cm -3 ), and n the number of soil layers.
The normality test of our data on organic and inorganic carbon shows a normal distribution. Two-way analyses of variance followed by least signi cant differences (LSD) were performed to evaluate the differences of various soil indices between basins and layers. Linear regression analyses were employed to analyze the SOC-SIC relationship, and the relationship between DOC:SOC ratio and other soil properties. Statistical analysis was executed using SPSS 19.0.

Results
Basic soil properties and spatial variations of water-soluble Ca+Mg density As presented in Fig. 1, soil pH had a large range in both LFB (7.9~8.9) and YCB (8.4~9.1); mean soil pH was signi cantly lower in the former (<8.4) than in the latter (>8.6) ( Table 1). But EC value was observably high in the LFB (0.49-0.72 mS/cm) than in the YCB (0.19-0.29 mS/cm), with higher values in the subsoil. Soil C:N ratio was similar between LFB (9.2-13.4) and YCB (10.2-12.4) above 80 cm, but signi cantly lower in LFB (7.3) than in YCB (9.3) over 80-100 cm. Fig. 2 showed that water-soluble Ca+Mg (the sum of water-soluble Ca 2+ and Mg 2+ ) was generally lower in topsoil (i.e., 98 to 341 g m -3 ) than in subsoils (i.e., 133-450 g m -3 over 20-40 cm, and 134-465 g m -3 over 40-100 cm). Overall, water-soluble Ca+Mg was low at sites in the north of LFB, and at the sites close to the Yellow River in the YCB above 40 cm. Interestingly, water-soluble Ca+Mg showed similar spatial pattern and magnitudes over 40-100 cm to those over 20-40 cm in the YCB, but overall higher levels in the LFB. There were no signi cant differences in water-soluble Ca+Mg between two basins or over depth (Table 1).
Spatial variations of SOC, DOC density and DOC:SOC ratio Fig. 3 shows SOC had a large range in topsoil (5.9-23.0 kg C m -3 ) than in subsoils (4.3-16.7 kg C m -3 over 20-40 cm, and 3.1-12.1 kg C m -3 over 40-100 cm). Overall, SOC density was greater in LFB than in YCB, especially in subsoils, although there were a couple of sites (along the Fenhe river and Yellow River) showing high levels of SOC over 20-40 cm in YCB. It appeared that SOC density was lower at high-altitude sites, in particular over 0-20 cm.
DOC varied from 64 to 147 g C m -3 in LFB, and 45 to 118 g C m -3 in YCB over 0-20 cm (Fig. 3b), showing a similar spatial distribution with SOC (Fig. 3a). However, the spatial pattern of DOC was somehow different to that of SOC in subsoils. Overall, DOC level in subsoil was higher in the YCB than in the LFB. For example, DOC over 40-100 cm varied from 46 to 106 g C m -3 in LFB, but 48-144 g C m -3 in the YCB.
Ratio of DOC:SOC ranged from 0.37 to 0.98 over 0-20 cm, 0.49 to 1.87 over 20-40 cm, and 0.71 to 2.52 over 40-100 cm (Fig. 4). Despite of a large spatial variability, DOC:SOC ratio showed general higher values in YCB than in LFB in particular below 40 cm. As showed in Table 2, mean DOC:SOC ratio was slightly higher in the YCB (0.76-1.09) than in the LFB (0.61-0.95) above 40 cm, signi cantly higher in the YCB (1.24) than in the LFB (1.88) above 40 cm. It appeared that lower ratio of DOC:SOC was found at sites in high-altitude and far away from the Yellow River.
Spatial variations of SIC and DIC density SIC also showed large spatial variation, with a range from 1.5 to 37.8 kg C m -3 above 20 cm, 7.7 to 25.5 kg C m -3 over 20-40 cm, and 5.3 to 28.2 kg C m -3 over 40-100 cm (Fig. 5ace). Overall, lower levels of SIC were found in the north LFB; however, the lowest SIC (< 5 kg C m -3 ) was at high-altitude sites in the YCB over 40-100 cm. Nevertheless, the highest levels of SIC (> 37 kg C m -3 ) were found in the YCB for all soil layers. It appeared that the spatial variability of SIC was larger over 20-40 and 40-100 cm than over 0-20 cm.
The spatial pattern of DIC was markedly different between the three layers although the range was similar, i.e., from 34 to 85 g C m -3 over 0-20 cm, 37 to 96 g C m -3 over 20-40 cm, and 38 to 98 g C m -3 over 40-100 cm (Fig. 5bdf). Overall, the spatial distribution of DIC showed no similarity to that of SIC (Fig. 5). For the 0-20 cm layer, high DIC (> 50 g C m -3 ) was found at sites in the north and south part of LFB and the southwest of YCB. However, DIC below 40 cm displayed high values (>60 g C m -3 ) at sites in the east of YCB. In general, low values of DIC were found at sites near the Yellow River below 20 cm (Fig. 5df).
The relationship between carbon fractions Fig. 6 shows a signi cantly negative correlation between SIC and SOC stocks using entire dataset (30 sites) over both 0-20 cm (P < 0.05) and 0-100 cm (P < 0.01). On a regional scale, there was a signi cantly negative correlation between SIC and SOC stocks in LFB over the 0-20 cm (P < 0.01) and 0-100 cm (P < 0.05), but weak or no signi cant correlation in YCB. Our analysis also showed a strong positive relationship (P < 0.01) between DOC and SOC over 0-20 cm in both LFB and YCB. However, for the 20-100 cm layer, DOC had a signi cantly positive correlation with SOC only in the YCB (Fig. 7).

Discussion
Variations of SOC and SIC in loess cropland On average, SOC content is signi cantly higher in whole soil pro le in the LFB (4.35-15.1 g kg -1 ) relative to YCB (3.14-10.3 g kg -1 ) ( Table 2). Previous studies have demonstrated that erosion can transport topsoil from highlands to lowlands which leads to depositions of SOM/SOC in valleys and basins with lower elevation , Zheng et al. 2005, Zhong & Xu 2009). The LFB, as a valley, may have accumulated a thick layer of redeposited topsoil from surrounding highlands that contain high levels of SOC .
SOC content of the YCB (from 10.3 g kg -1 over 0-20 cm to 3.1 g kg -1 over 80-100 cm) is similar to those in loess cropland, e.g., from 9.7 to 2.7 g kg -1 in the Hebei Plain , and from 11.0 to 3 g kg -1 in the west part of Loess Plateau , Zhang et al. 2015. The relative low levels of SOC in those loess croplands may be largely related to the texture of loess that has less clay thus less protection for SOM (Li et al. 2019, Li et al. 2017).
Our study showed that SIC was signi cantly lower in LFB (9.2-10.6 g kg -1 ) than in YCB (11.4-13.4 g kg -1 ) above 80 cm, which corresponded with signi cantly higher DIC:SIC ratio in the LFB (0.91%) than in the YCB (0.39%) in the topsoil (Table 2). There were some studies on SIC dynamics in the croplands of north China that have the same or similar parent materials, such as in the Loess Plateau (Zhang et al. 2015) and the North China Plain (Guo et al. 2016, Shi et al. 2017b. Clearly, SIC content in the LFB (with redeposited topsoil) is much lower than that in the western highland (>15 g kg -1 ) of Loess Plateau (Zhang et al. 2015) that has much drier climate and higher soil pH (Table 3). Interestingly, SIC levels in the YCB are close to those (10.5-12.7 g kg -1 ) in the upper YRD (Guo et al. 2016), but modest higher than those (7.1-11.4 g kg -1 ) in the North China Plain (Shi et al. 2017b) though all soils had same parent materials and similar climatic conditions (Table 3). Both the YCB and upper YRD are close to the Yellow River, thus experience strong in uences of hydrological processes that could supply extra Ca and Mg ions over the past, which is bene cial for SIC formation (Guo et al. 2016).

The relationship between SIC and SOC in loess
Our analysis shows no clear relationship between SIC and SOC stocks in YCB, but a signi cant negative correlation between SIC and SOC stocks in the LFB. The overall negative correlation of SIC and SOC in this study area disagrees with our previous ndings of a positive correlation in the loess cropland of North China Plain (Guo et al. 2016, Shi et al. 2017b, and under various land uses in northwest China , Wang et al. 2015b). There is also evidence of a negative SIC-SOC relationship in the cropland of Hebei Plain (Li et al. 2010) and under various land use types in the Loess Plateau, north China (Zhao et al. 2016).
We found a negative SIC-SOC relationship in the LFB, which was consistent with a report for another part of the Loess Plateau (Zhao et al. 2016). Soil erosion and re-deposition/redistribution is profound in the Loess Plateau (Zheng et al. 2005), which would move topsoil from highland to lowland, leading to enhanced SOC storage (but with lower SIC stock) in upper 100 cm (Table 3). In addition, elevated SOC level in the soil pro le of LFB could also result in more CO 2 production thus lower pH (Table 1), which causes dissolution of soil carbonate (Chang et al. 2012, Raheb et al. 2017).
The non-signi cantly SIC-SOC relationship in YCB may re ect the complex in uences of multi drivers. There is evidence that precipitation has large impacts on both SOC and SIC in arid/semi-arid lands , Raheb et al. 2017, Wu et al. 2009), particularly in the Loess Plateau (Han et al. 2018). The YCB has a large spatial variation in both precipitation (500-750 mm) and elevation (380-820 m a.s.l.), which could have large impacts on the distributions of SOC and SIC, leading to alterations of the SIC-SOC relationship. Our analyses indicate that there is an increasing trend (from 1.0 to 5.8) in SIC:SOC ratio with increasing elevation in the YCB. In addition, the distance from the Yellow River is probably another factor that can alter the SIC-SOC relationship through the in uences of hydrological processes on either SIC (Shi et al. 2017b) or SOC (Shi et al. 2017a). We further discuss the potential in uence of hydrological processes on SOC below.

Regulating factors for SOC
The variability of SOC is large under the same cropping system in the loess of north China (Fig. 3 and Shi et al. 2017b), implying that "outputs" of SOC/SOM (e.g., decomposition of SOM), rather than "inputs", primarily regulate the dynamics of SOC, which is in uenced by environmental conditions , Shirale et al. 2019. Given that subsoils are mainly subject to decomposition, we evaluate the relationships between soil carbon indices and other variables over 40-100 cm. Our analyses show that SOC has a positive correlation with EC, and a negative correlation with soil pH (Fig. 6). There is evidence that high EC in saline soils can cause occulation of clay particles into aggregates, which restricts substrate availability for microbial thus retards decomposition of SOM (Wong et al. 2010). Additionally, soils with high pH often have poor physical-chemical conditions that are harmful to crop growth and root system development, resulting in less organic carbon inputs into soil (Kemmitt et al. 2006, Wong et al. 2010).
In addition to decomposition process, desorption could also lead to lower SOC (Kalbitz et al. 2000, Mavi et al. 2012. Our analyses show that DOC:SOC ratio (representing desorption potential) is positively correlated with pH, and negatively correlated with EC in subsoils (Fig. 8). There are studies demonstrating that high soil pH can destroy soil aggregates, thus reduce the protection of SOM against degradation (Kalbitz et al. 2000, Tavakkoli et al. 2015 whereas high EC may be bene cial for the formation and stabilization of soil aggregates (Rahimi et al. 2000).
Our analyses show that DOC:SOC ratio is higher in YCB (0.76-1.88 %) than in LFB (0.61-1.24 %), with a signi cant difference below 40 cm (Table 2), which indicates stronger desorption of SOC in the subsoils of YCB. Similarly, Zhang et al. (2020) reported greater DOC:SOC ratio close to the Yellow River over 60-100 cm, relative to those in other parts of North China Plain. Apparently, the YCB and other regions that have shorter distance to the Yellow River are in uenced by stronger hydrological processes (such as water movements), which would result in more desorption and remove of DOC thus lower levels of SOC in subsoils (Zhang et al. 2020).

Conclusions
We evaluate the spatial variation, relationship and driving factors of soil carbon fractions over 0-100 cm in low-elevation cropland of the Loess Plateau, i.e., LFB and YCB. Our data show that SIC stock is signi cantly higher than SOC stock over 0-100 cm in the whole study area; SOC stock is negatively correlated with soil pH, but positively correlated with EC. SOC and SIC stocks reveal a signi cantly negative correlation in LFB, which is partly related to erosion and re-deposition/redistribution of topsoil. We nd no clear relationship between SOC and SIC stocks in YCB, which may re ect the large spatial variations of elevation and precipitation.
Our data show a signi cantly positive correlation between DOC and SOC in both topsoil and subsoil of YCB, but only in topsoil of LFB. We nd that DOC:SOC ratio is signi cantly higher below 40 cm in YCB than LFB, and DOC:SOC ratio is positively correlated with soil pH, and negatively correlated with EC. Our analyses suggest that high soil pH and stronger hydrological processes are attributable to the relatively lower levels of SOC in YCB. Further studies are needed to investigate how soil properties and environmental conditions interplay to regulate the dynamics of SOC and SIC in the Loess Plateau.

Declarations
Availability of data and materials Xiujun Wang provided supervision and nancial support for this study, and checked/corrected all the versions of the manuscript.
Tongping Lu collected soils, conducted laboratory measurements and data analyses, and prepared for the manuscript. Wenxi Zhang helped with sampling and analyses, and commented on earlier versions of the manuscript. All authors read and approved the nal manuscript.  Values followed by the same letter (capital letter between two basins and lowercase letter between layers) are not signi cantly different at P < 0.05 based on LSD test.   Water soluble Ca+Mg density (g m-3) over (a) 0-20 cm, (b) 20-40 cm, and (c) 40-100 cm in the YCB (red) and LFB (green). The gure was generated using ArcMap 10.5 (http://www.esri.com/).  Soil inorganic carbon (SIC) (left panel) and dissolved inorganic carbon (DIC) (right panel) over (a, b) 0-20 cm, (c, d) 20-40 cm, and (e, f) 40-100 cm in the YCB (red) and LFB (green). The gure was generated using ArcMap 10.5 (http://www.esri.com/).

Figure 6
Correlation between SIC and SOC stocks in LFB, YCB and combined data (black dash lines) over (a) 0-20 cm, (b) 0-100 cm. One asterisk indicates signi cance at P < 0.05, and two asterisks at P < 0.01.

Figure 7
Correlation between DOC and SOC content in the LFB, YCB and combined data (black dash lines) over (a) 0-20 cm and (b) 20-100 cm. One asterisk indicates signi cance at P < 0.05, and two asterisks at P < 0.01.