Next Article in Journal
Viability of Entomopathogenic Fungi in Oil Suspensions and Their Effectiveness against the Agave Pest Scyphophorus acupunctatus under Laboratory Conditions
Previous Article in Journal
The Chemical Capping Regulation Mechanism of Cotton Main Stem Growth
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long-Term Chemical and Organic Fertilization Differently Affect Soil Aggregates and Associated Carbon and Nitrogen in the Loess Plateau of China

1
Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, College of Urban and Environmental Science, Northwest University, Xi’an 710127, China
2
Shaanxi Key Laboratory for Carbon Neutral Technology, Northwest University, Xi’an 710127, China
3
Northern Plains Agricultural Research Laboratory, USDA-ARS, 1500 North Central Avenue, Sidney, MT 59270, USA
4
College of Resource and Environmental Sciences, Hebei Agricultural University, Baoding 071001, China
5
State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Xianyang 712100, China
*
Author to whom correspondence should be addressed.
Submission received: 26 April 2023 / Revised: 15 May 2023 / Accepted: 23 May 2023 / Published: 25 May 2023
(This article belongs to the Section Soil and Plant Nutrition)

Abstract

:
Fertilizer sources may have variable effects on soil aggregation, aggregated-associated C and N, and wheat yield. A 34-year field experiment was performed to evaluate the influences of chemical and organic fertilization on soil aggregates and associated carbon and nitrogen under winter wheat in a Cumulic Haplustoll of the Loess Plateau, China. Treatments included unfertilized control (CK), inorganic N fertilizer (NF), inorganic P fertilizer (PF), inorganic N and P fertilizer (NP), organic manure (M), inorganic N fertilizer plus manure (NM), inorganic P fertilizer plus manure (PM), and inorganic N plus P fertilizers plus manure (NPM). Compared to CK, long-term fertilization significantly increased the proportion of soil macro-aggregates, mean weight diameter (MWD), and geometric mean diameter (GMD), but decreased the proportion of micro-aggregates and fractal dimension, especially fertilizer plus manure. Compared to CK, manure treatments (M, NM, PM, and NPM) had a better improvement on soil organic carbon (SOC), soil total nitrogen (STN), particle organic C, and microbial biomass C in all aggregates than the fertilizer alone. The SOC in different aggregates increased with the increased aggregate size, which was because the larger aggregates formed by the binding of the smaller aggregates and organic matter. PON increased in NM and NPM, and MBN was more sensitive to N fertilizer. The C/N ratio in bulk soil and aggregates decreased with fertilization, especially in fertilizer with manure and in macro-aggregates. The improved soil structure was related to the increased SOC and STN, which was proved by the positive correlations among SOC and STN with macro-aggregates and MWD. A correlation analysis also showed that the contribution rate of SOC and STN in macro-aggregates was positively associated with the macro-aggregate and stability. Therefore, the sequestration of C and N in soil was related to aggregate size and was mainly affected by larger aggregates. The results demonstrated that fertilizer with manure improved the soil structure and fertility better than fertilizer alone, thus increasing crop yield.

1. Introduction

Maintaining soil good structure and high organic matter is essential for the continuous supply of soil nutrients, plant growth, and ecosystem balance [1,2]. Soil aggregation promotes root growth and microbial activity by optimizing the soil structure [3]. Materials secreted by roots and microorganisms that act as adhesives can promote soil aggregations [4,5]. Soil organic carbon (C) and nitrogen (N) are the core substances of nutrient cycling, and the sequestration of C and N in soil is considered to be a potential mechanism for reducing the emissions of greenhouse gases in agricultural systems [6,7]. It has been reported that nearly 90% of soil organic C is in the aggregates [3,8]. Aggregates protect soil C and N from mineralization, and soil organic C and N provides the adhesives for the formation of aggregates [9]. Therefore, the changes in soil aggregates affect the sequestration of C and N in soil, and C and N pools that are protected by soil aggregates are more stable than that in the unprotected pools [10].
The storage and transformation of C and N in aggregates vary with the aggregate size [11]. Soil organic C (SOC) and total N (STN) can be greater in macro-aggregates (>2 mm) than in micro-aggregates (<0.25 mm) because fresh crop residues are incorporated first in macro-aggregates [9,12,13]. However, Udom et al. [14] reported that SOC can be higher in micro-aggregates because of the more stable structure and slower decomposition. Large aggregates provide an ecological niche for storing labile organic C while small aggregates contain more stabilized C [15,16,17]. Therefore, C and N in large aggregates have greater dispersion and turnover rates under the influence of external conditions than those in small aggregates [18,19,20]. Particle organic C (POC) and N (PON) are intermediate fractions of C and N that are more labile than SOC and STN, which provide substrates for soil microorganisms and contribute to aggregate stability [21,22]. Soil microbial biomass C (MBC) and N (MBN) are labile C and N fractions that are stored in the microbial body and change rapidly within a growing season [11,16,23]. Therefore, it is necessary to know how aggregate stability and turnover under management practices affect C and N storage, mineralization, and availability for plant uptake.
Fertilization, as an important management to improve crop yield, has a great impact on soil aggregates and associated C and N fractions [17,24]. Previous studies have shown that frequent use of chemical fertilizers deteriorates soil aggregation and reduces the formation and stability of large aggregates, thereby exposing previously protected organic C [25]. In contrast, the application of organic fertilizer is beneficial to the formation of large aggregates that are rich in organic C and increases the retention of organic C [16,24,26]. Wang et al. [27] found that the combined application of chemical and organic fertilizer had a positive effect on improving aggregate stability and nutrient accumulation in dryland greenhouses. In addition, Duan et al. [28] found that the long-term application of manure led to an increase in macro-aggregates (>250 μm) and the contents of plant available N in aggregates, thus contributing to high yield. However, Das et al. [29] and Hati et al. [30] reported the long-term application of chemical fertilizers improved macro-aggregation and stability of aggregates. Therefore, the response of soil aggregates’ associated C and N fractions to fertilization practices may vary with the fertilizer types and soil properties [2,31].
In the Loess Plateau of China, different fertilization practices, mainly with chemical fertilizers, have been used since the 1970s to improve soil quality and enhance winter wheat (Triticum aestivum L.) yield [32]. This has degraded the soil quality by reducing aggregation, resulting in a decline in the crop yield [32,33]. Improved management strategies are needed to enhance soil aggregation, quality, and crop yields. This study evaluated a 34-year field experiment that received chemical fertilizers, organic manure, and a combination of both on soil aggregation and associated C and N in a Cumulic Haplustoll. The objectives were to determine the long-term effect of chemical fertilizers and organic manure on soil aggregation and stability, aggregated-associated C and N, and winter wheat yields in the Loess Plateau of China. We hypothesized that organic manure, with or without chemical fertilizers, would enhance soil aggregation and stability, the associated C and N, and the winter wheat yield compared to chemical fertilizers alone or no fertilization.

2. Materials and Methods

2.1. Site Description

A field experiment was established in 1984 that evaluated the effects of chemical fertilizers and organic manure on winter wheat yield at the Shaanxi Changwu Agro-Ecological Station of China (107°44.70′ E, 35°12.79′ N) [32,34]. The site (elevation 1200 m) had a sub-humid temperate climate with an average (30 year) annual precipitation of 584 mm, an average annual air temperature of 9.1 °C, and a frost-free period of 171 days. The soil was classified as an aridic and loamy, Cumulic Haplustoll, according to USDA Soil Taxonomy [35], with sand, silt, and clay concentrations of 35, 656, and 309 g kg−1, respectively, pH 8.4, soil organic C (SOC) 6.5 g kg−1, total N (STN) 0.80 g kg−1, total phosphorus (P) 1.26 g kg−1, and CaCO3 105 g kg−1 at the 0–30 cm depth at the initiation of the experiment.

2.2. Experimental Design

Treatments were (1) unfertilized control (CK), (2) inorganic N fertilizer (NF), (3) inorganic P fertilizer (PF), (4) inorganic N and P fertilizer (NP), (5) organic manure (M), (6) inorganic N fertilizer plus organic manure (NM), (7) inorganic P fertilizer plus organic manure (PM), and (8) inorganic N and P fertilizers plus organic manure (NPM). All treatments were arranged in a randomized block design with three replications. The plot size was 10.3 m × 6.5 m, with a 0.5 m strip between plots. Following the local practices, N fertilizer as urea (46% N), P fertilizer as triple superphosphate (45% P), and organic manure were applied, and the doses of N and P were at a rate of 120 and 90 kg ha−1, respectively, in 75 t ha−1 of organic nutrients. The organic manure had an average C, N, and P concentrations of 800, 87, and 44 kg ha−1, respectively. Because the soil was rich in K concentration, no K fertilizer was applied. However, organic manure also supplied K at 30 kg K ha−1.
Before planting, chemical fertilizers and organic manure were broadcast and incorporated at a depth of 20 cm with a rotary tiller in late September each year [33]. Seed beds were prepared by cultivating soil at a depth of 20 cm using a rotary tiller. Immediately after tillage operation, winter wheat (cv. Changwu 131) was sown at 195 kg ha−1 with a 20 cm row-spacing. No irrigation was applied. Weeds were controlled manually and pesticides were applied as needed. Winter wheat was harvested in June of the following year and grain yield was determined by cutting all plants at 2 cm above the ground with a sickle within the plot, threshing plants on the ground, and oven drying at 70 °C for 7 days.

2.3. Soil Sampling

Soil samples were collected from 0–20 cm after wheat harvest in late June of 2018. In each plot, five cores were collected randomly with a hand probe (5 cm inside diameter), and composited within a plot. After air-drying, one part was gently broken along the natural destruction surface to obtain aggregates of <10 mm for the aggregate analysis. Another part was sieved to 2 mm to perform C and N analysis.

2.4. Soil Aggregate Analysis

Dry sieving is an appropriate method that separates aggregates in dryland cropping systems with limited precipitation where aggregates are disrupted more by wind than water erosion [7]. Briefly, we prepared a set of upper and lower sets of 2.00 mm and 0.25 mm, respectively. A total of 100 g of air-dried soil samples was dry sieved through a nest of sieves containing 2 mm and 0.25 mm sieves and vibrated mechanically for 5 min in a shaker. Aggregates that remained in each sieve after sieving (10–2 mm, 2–0.25 mm, and <0.25 mm) were collected and weighed. The stability index of soil aggregates included aggregate content >0.25 mm (R0.25), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (D). The MWD (mm), GMD (mm), and D were estimated using the following equations [36]:
MWD = i = 1 3 m i d i i = 1 3 m i
GMD = E x p i = 1 3 ( m i ln d i ¯ ) i = 1 3 m i
m d < d i i = 1 3 m i = d i d m a x 3 D
where mi is the mass of aggregate fraction i (g); di is the mean diameter of the aggregate fraction i (mm); m(d<di) is the mass of aggregate fraction with diameter of <di (g); and dmax = 10 mm.

2.5. Soil C and N Fractions Analyses

Soil organic carbon (SOC) and total nitrogen (STN) in bulk soil and aggregates were analyzed using a C and N analyzer (Euro Vector EA3000, Manzoni, Italy) after further grinding the bulk soil and aggregates to 0.15 mm and pretreating with 1.0 mol L−1 HCl to remove inorganic C. The POC and PON in aggregates were determined using the method shown by Cambardella and Elliott [37]. A total of 10 g of aggregates was mixed with 30 mL of 5 g L−1 sodium hexametaphosphate solution and shaken in a reciprocating oscillator for 16 h. The dispersed solution was passed through a 0.053 mm sieve, dried at 50 °C, and analyzed for C and N concentrations with a C and N analyzer, as described above. The POC and PON concentrations were calculated as the difference between the SOC or STN concentrations in the bulk soil and the particle concentration that passed through the sieve after correction for sand concentration. Soil microbial biomass C and N (MBC and MBN) were determined using the chloroform fumigation leaching method [38]. The contribution rate (CR) and enrichment coefficient (EC) of SOC and STN in soil aggregates were calculated by the following equations [39]:
CR = S O C i ( or   S T N i ) × M i SOC ( or   STN ) × 100 %
EC = S O C i ( or   S T N i ) SOC ( or   STN ) × 100 %
where SOCi (or STNi) is the content of SOC and STN in aggregate (g kg−1); Mi is the proportion of the aggregate (%); and SOC (or STN) is the content of SOC (or STN) in the bulk soil (g kg−1).

2.6. Statistical Analysis

The SPSS 19.0 (SPSS, Inc., Chicago, IL, USA) was used for the statistical analysis of data. A one-way analysis of variance (ANOVA) was carried out to determine significant differences between treatments using LSD’s test with a significance level of p < 0.05. Treatment was considered as the main plot and aggregate size was considered as split-plot treatment for data analysis. Correlation heatmap was constructed using the pheatmap-package and R package vegan, based on the Spearman correlation matrix. Co-occurrence networks among soil aggregates, the aggregate stability, and aggregate-associated C and N fractions were constructed according to the Pearson’s correlation analysis and p value (p < 0.01, Bonferroni-corrected). The software Gephi 0.9.2 was used to visualize the networks of co-occurring relationships [40].

3. Results

3.1. Soil’s Basic Properties

Compared to CK, fertilizer with manure decreased the soil bulk density (BD) and PM decreased by 7.6% (p < 0.05) (Table 1). Inversely, the soil porosity increased by 3.6–6.6% in fertilizer with manure and PM was significant. Compared to CK, all fertilized treatments increased WHC and fertilizer with manure showed a more obvious effect than the fertilizer alone. M, NM, PM, and NPM increased WHC by 36, 39, 38, and 47% (p < 0.05), respectively. All fertilized treatments, except NF, had a greater SOC than CK, and the SOC was 19–54% greater in fertilizer with manure than fertilizer alone. Similarly, compared to CK, all fertilized treatments had a greater STN, which was 46–67% greater in fertilizer with manure than fertilizer alone. As a result, the soil C/N ratio decreased by 17–37% in fertilized treatments compared to CK and NP, and the fertilizer with manure had more obvious changes.

3.2. Aggregate Distribution and Stability

Compared to CK, the proportion of macro-aggregates (10–2 mm) increased by 28, 13, 35, 45, and 61% in NP, M, NM, PM, and NPM (p < 0.05, Figure 1), respectively. The proportion of medium aggregates (2–0.25 mm) decreased in all fertilized treatments, which decreased by 9.2, 11, 19, and 22% (p < 0.05) in N, NM, PM, and NPM, respectively. The proportion of micro-aggregates (<0.25 mm) was opposite to the changes in macro-aggregates, which decreased by 37, 22, 34, 29, and 49% (p < 0.05) in NP, M, NM, PM, and NPM, respectively. Table 2 showed that the R0.25 increased by 7.8, 4.5, 7.1, 6.1, and 10% (p < 0.05) in NP, M, NM, PM, and NPM, respectively, compared to CK. Meanwhile, NP, M, NM, PM, and NPM significantly increased MWD and GMD by 12–51% and 7.8–29%, respectively. The fractal dimension (D) decreased by 4.2, 4.9, and 7.2% (p < 0.05) in NM, PM, and NPM, respectively. These results indicated that soil aggregate stability improved better in fertilizer with manure than in fertilizer alone, especially in NPM.

3.3. C and N Fractions in Aggregates

Compared to CK, the SOC in M, NM, PM, and NPM increased by 3.1–4.7% and 3.3–4.8% in macro- and medium-aggregates, respectively, and all fertilized treatments increased SOC by 1.1–5.1% in micro-aggregates (p < 0.05, Figure 2a). The POC in M, NM, PM, and NPM increased by 85–115%, 30–67%, and 3.3–4.8% in macro-, medium-, and micro-aggregates, respectively, while NF, PF, and NPF only increased POC in macro- and micro-aggregates (p < 0.05, Figure 2b). Compared to CK, MBC in NM, PM, and NPM increased by 15–40%, 6.7–33%, and 12–54% in macro-, medium-, and micro-aggregates, respectively, and M only increased MBC in macro-aggregates (p < 0.05, Figure 2c). However, NF, PF, and NPF decreased MBC in different aggregates. Compared to CK, NP, M, NM, PM, and NPM had significant improvements on aggregated-associated STN, which increased by 39–66%, 36–60%, and 47–76% in macro-, medium-, and micro-aggregates, respectively (p < 0.05, Figure 2d). The PON in different aggregates significantly increased in NM and NPM, while other treatments even decreased in macro- and medium-aggregates (Figure 2e). MBN showed the highest value in NPM and increased by 49, 48, and 59% in macro-, medium-, and micro-aggregates, respectively (p < 0.05, Figure 2f).

3.4. C and N Distributions in Aggregates

Compared to CK, all fertilized treatments decreased the C/N ratio in different aggregates (Table 3). The fertilizer with manure had a more significant effect than fertilizer alone, which reached the lowest value in NPM and decreased by 37, 35, and 41% (p < 0.05) in macro-, medium-, and micro-aggregates, respectively. Compared to CK, CR of SOC (C-CR) increased in macro-aggregates and decreased in medium- and micro-aggregates, especially in fertilizer with manure (Table 4). However, the EC of SOC (C-EC) in fertilizer with manure had no significant changes, nor did it decrease. In the macro-aggregates, CR of STN (N-CR) increased by 54, 38, 29, 42, and 66% (p < 0.05) in NP, M, NM, PM, and NPM, respectively. However, N-CR in medium-aggregates only increased in NP and M, while it decreased by 18–23% (p < 0.05) in NM, PM, and NPM. Fertilizer with manure also decreased N-CR in the micro-aggregates by 25–44% (p < 0.05). Consistent with N-CR in the medium-aggregates, the EC of STN (N-EC) significantly increased in NP and M in different aggregates. Meanwhile, N-EC also significantly increased in the micro-aggregates in PM and NPM.

3.5. Correlation between Aggregate Stability and the C and N Fractions

Figure 3 shows the correlations between the soil’s basic properties and the aggregate parameters. Soil porosity and WHC were positively correlated with the macro-aggregates, R0.25, MWD, and GWD, while they were negatively correlated with the micro-aggregates and fractal dimension (D). However, BD showed opposite relations with the aggregate parameters compared to porosity and WHC. SOC and STN were positively correlated with WHC, macro-aggregates, MWD, and GWD, while they were negatively correlated with the fractal dimension (D). Meanwhile, SOC and STN were negatively correlated with the medium- and micro-aggregates.
The correlations between the aggregate size and stability of bulk soil and the C and N fractions in different aggregates are shown in Figure 4. In macro-aggregates, SOC and MBC were positively correlated with the SOC of bulk soil, with coefficients of 1.00 and 0.77, respectively. Meanwhile, SOC, POC, MBC, and C-CR in macro-aggregates were positively correlated with MWD, with coefficients of 0.71, 0.54, 0.80, and 0.72, respectively. SOC in medium- and micro-aggregates and POC in micro-aggregates showed the highest positive correlations with the STN of bulk soil, with coefficients of 0.77 and 0.79, respectively. Similarly to macro-aggregates, MBC in medium- and micro-aggregates showed positive correlations with the SOC of bulk soil. However, C-CR in the medium- and micro-aggregates was negatively correlated with MWD, with coefficients of 0.84 and 0.82, respectively. STN in different aggregates was positively correlated with MWD (0.82, 0.83, and 0.83) and GMD (0.85, 0.86, and 0.86). The PON in aggregates was positively correlated with the SOC of bulk soil, with coefficients of 0.51, 0.48, and 0.46 in macro-, medium-, and micro-aggregates, respectively. N-CR was positively correlated with GWD (0.90) in macro-aggregates while negatively correlated with GWD (0.87) in macro-aggregates. The C/N ration in different aggregates was positively correlated with D (0.80, 0.82, and 0.81) and negatively with GMD (0.81, 0.82, and 0.81).

3.6. Average Wheat Yield

The winter wheat yield, averaged across years, had more obvious improvements in fertilizer with manure than the fertilizer alone (Figure 5). Compared to CK, N and NP increased the yield by 84 and 140% (p < 0.05), respectively, while P had no significant effect. M, NM, PM, and NPM increased the yield by 99, 96, 166, and 147%, respectively, compared to CK.

4. Discussion

4.1. Effects of Long-Term Fertilization on SOC and STN in Bulk Soil

The manure treatments increased SOC and STN while reducing the C/N ratio more than the chemical fertilizer, especially the NPM (Table 1). The total SOC mainly depended on the balance of the input and degradation of organic C [34]. The chemical fertilizer plus manure could increase the soil microbial activity on the one hand, thus promoting crop growth and biomass, which could increase the return amount of root stubble and secretion [41]. On the other hand, manure directly promoted increases in SOC by the direct input of C [6,42]. However, the study reported that the input of SOM by chemical fertilizer was not enough to make up for the loss of mineralization, and N fertilizer provided a soil microbial N source effectively, which could improve the soil microbial activity and strengthen the decomposition of organic C by physical protection [30]. Thus, the chemical fertilizer alone had no significant changes in SOC, or the promotion was poorer (Table 1). STN increased significantly with the application of N fertilizer. The soil C/N ratio decreased after fertilization, especially with chemical fertilizer plus manure, which was closely related to the improvement of organic matter decomposition and the application of N fertilizer.

4.2. Effects of Long-Term Fertilization on Aggregate Proportion and Stability

The greater proportion of macro-aggregates as well as increased MWD and GWD for NPM than other treatments (Figure 1 and Table 2) were likely due to higher C and N concentrations (Figure 2). These results were consistent with previous findings that the application of manure increased soil aggregation and stability compared to chemical fertilization [43,44]. The increased soil aggregation and stability could be related to the increased SOM, which acted as a binding agent for soil particles to form aggregates [4,20,45,46], as proved by the positive relationships between SOC and macro-aggregates, MWD, and GWD (Figure 3). Although chemical fertilizer could increase cementing substances indirectly by improving the crop biomass, this effect was far less obvious than with organic fertilizer [47]. Organic fertilizer could not only improve crop biomass by providing various nutrients, but it also contained a lot of organic cementing substances [6,48,49]. The increased microbial biomass and fungal mycelia, as shown by greater POC and MBC, may increase the binding agents and promote soil aggregation for NPM [42]. Manure alone also increased macro-aggregates compared to chemical fertilization; however, the greater macro-aggregate proportion in NPM was also likely due to the increased winter wheat yield (Figure 3) that returned a greater crop residue to the soil, which provided more C inputs along with that provided by organic manure. This subsequently reduced the aggregate proportion for NPM compared to other treatments in medium- and micro-aggregates. Therefore, supplementing exogenous organic matter was an important measure to improve soil structure and maintain soil fertility.

4.3. Effects of Long-Term Fertilization on Aggregate-Associated Carbon and Nitrogen

The manure treatments increased the SOC, POC, and MBC in all aggregates, while chemical fertilizer only increased the POC in macro- and micro-aggregates and decreased the MBC in all aggregates (Figure 3). In addition to N, P, and K, organic manure also supplied 800 kg C ha−1 annually. As a result, more C may have been stored in aggregates that were in the treatments receiving manure. The SOC in different aggregates increased with the increased aggregate size. Larger aggregates were formed by the binding of small, free aggregates to temporary cementing agents, and the structure of larger aggregates could better protect the organic C and its components; thus, the larger aggregates contained more organic C [18,50,51]. Larger soil aggregates were also formed by the cementation of organic residues and mycelia, and newly imported organic C that was easily mineralized and decomposed mainly existed in large aggregates and was physically protected, while small aggregates were formed by the cementation of polysaccharides or inorganic colloid [52,53]. However, the chemical fertilizer resulted in a low content of newly input organic C, which reduced cementation material and contents of large aggregates, thus increasing the contribution rate of small aggregates to organic C [25]. At the same time, the effect of chemical fertilizer on POC was inferior to that of manure. On the one hand, POC lacked the physical protection due to poor soil aggregate structure and was easily used by microorganisms. On the other hand, chemical fertilizer could provide direct and available nutrients for microorganisms to enhance microbial activity and accelerate the conversion of POC to humic compounds [21]. In addition, larger aggregates could maintain a higher nutrient level than micro-aggregates and powder clay particles, which was conducive to microbial colonization [54]; thus, the MBC content in large aggregates was higher under the application of manure. A correlation analysis also showed that R0.25 was positively correlated with the contribution rate of SOC (C-CR) in macro-aggregates, while it was negatively correlated with C-CR in medium- and micro-aggregates (Figure 4), which confirmed the physical protection theory of aggregates to a certain extent.
The changing trends of STN in bulk soil and different aggregates were basically consistent with those of SOC, which was due to the fact that 95% N in soil existed in organic matter in the form of organic N; thus, the changes in C in organic matter would inevitably lead to changes in organic N [11,55] (Figure 3). STN in different aggregates significantly increased in NP, M, NM, PM, and NPM. The application of manure not only directly increased the input of N, but also indirectly increased STN through the increase in roots and exudates. MBN was more sensitive to the addition of N fertilizer, which was due to the microorganisms that were no longer limited by the N element after N input, thus increasing the transformation rate and improving the absorption and utilization of N [30]. The fertilizer plus manure was more effective in reducing the C/N ratio than the fertilizer and manure alone, especially in macro-aggregates (Table 3). With the input of exogenous nutrients, although SOC increased, the soil microbial activity also improved, which accelerated the decomposition of organic matter, thus resulting in the accumulation of N in soil faster than that of organic C [23,56]. Li et al. [57] showed that SOC and STN in subtropical paddy soil reserved in large aggregates (≥0.25 mm) accounted for 64–81% and 54–82% of the total soil, respectively, which was because C and N in large aggregates were not easily degraded by physical action, microorganisms, or enzymes; thus, it is more difficult to be mineralized and detained in soil [48]. A correlation analysis showed that N-CR was positively correlated with R0.25, MWD, and GMD in macro-aggregates, while negatively correlated with R0.25, MWD, and GMD in micro-aggregates (Figure 4).

4.4. Effects of Long-Term Fertilization on Winter Wheat Yield

A balance of nutrient availability throughout the growing season as well as improved soil aggregation and stability and associated C and N may have increased the winter wheat yield for PM and NPM compared to other treatments, except NP (Figure 3). It is likely that chemical fertilizers supplied nutrients at the beginning of wheat growth, while organic manure supplied them at the later stage of growth, resulting in a consistent availability of nutrients throughout the growing season and an enhanced wheat yield for PM and NPM [58,59]. The improved soil aggregation and stability and associated C and N (Figure 1 and Figure 2) may have also increased water and nutrient movements, root growth, and microbial activity due to greater energy and substrate availability, thereby increasing the wheat yield. As crops require both N and P for enhanced growth and production, the increased wheat yield for NP compared to CK, PF, M, and NM was probably due to the increased availability of both N and P. Soil aggregation and stability and associated C and N were also greater with NP than CK, PF, M, and NM (Figure 1 and Figure 2, Table 2), which may have had a positive effect on the wheat yield, as described above. A decline in the winter wheat yield along with the degradation of soil aggregation due to the continuous application of chemical fertilizers in the Loess Plateau of China have been reported by several researchers [33,34]. A lack of fertilization and manure applications probably reduced the wheat yield for CK. Similarly, the absence of N fertilization reduced the wheat yield for PF.
The decreased C/N ratio in the soil and aggregates indicated that the increased rate of C was lower than that of N, which would accelerate the decomposition of microorganisms and the mineralization rate of N, and was not conducive to the sequestration of C [9]. Weil and Kroontje [60] found that when the C/N ratio was close to 20, most of the organic matter was plant residue in the early stages of decomposition. When the C/N ratio was close to 10, most of the organic matter was humus (the refractory component lignin and the aromatic substances in small aggregates that bind closely to the clay). In this study, the C/N ratio was all close to or lower than 10, indicating that most soil organic matter components were in the later stage of decomposition under long-term fertilization, which converted to humus and existed in soil and aggregates. This result was also verified by the conclusion that the POC was greater in manure (Figure 2). Only the contribution rates of SOC and STN in macro-aggregates were positively correlated with R0.25 (Figure 4), indicating that the aggregate of >0.25 mm was more sensitive to fertilization management and was a good indicator of changes in organic C under agricultural management. Therefore, manure could improve soil fertility and physical properties with the increase in SOC, thus improving the soil’s ability to supply crop nutrients and increasing crop yield. A reasonable combination of chemical fertilizers and organic manure is an effective fertilization model for the sustainable development of agricultural production.

5. Conclusions

The long-term (34 years) application of chemical fertilizers and organic manure had variable effects on soil aggregation and stability, aggregated-associated C and N, and the winter wheat yield. Soil aggregate stability and the proportion of macro-aggregates increased, but that of micro-aggregates decreased. The manure treatments better improved the contents of SOC, STN, POC, and MBC in different aggregates. The content of SOC in aggregates increased with increased aggregate size, which was because the larger aggregates were formed by the binding of the smaller aggregates and organic matter. PON significantly increased in NM and NPM, and MBN was more sensitive to the addition of N fertilizer. However, the C/N ratio in bulk soil and aggregates decreased with fertilization. The improved soil structure was mainly due to the increased SOC and STN, which was proved by the positive correlation among SOC and STN with macro-aggregates and MWD. A correlation analysis also showed that the contribution rate of SOC and STN in macro-aggregates was positively correlated with aggregate size and stability, while it showed a negative relationship with micro-aggregates. Therefore, the sequestration of C and N fractions in soil was closely related to the aggregate size and was mainly affected by the larger aggregates. The chemical fertilizer plus manure could enhance the ability of the soil to supply nutrients and increase crop yield better than the chemical fertilizer by improving soil structure and fertility.

Author Contributions

Conceptualization, J.W.; methodology and software, C.L.; validation, C.Y.; formal analysis, X.F.; investigation and resources, C.L.; data curation, writing—original draft preparation, C.Y.; writing—review and editing, U.M.S. and F.Z.; visualization, C.Y.; supervision, project administration and funding acquisition, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 42277322, 31570440), the Key International Scientific and Technological Cooperation and Exchange Project of Shaanxi Province, China (2020KWZ-010), and the Shaanxi Agricultural Science and Technology Innovation-Driven Project (NYKJ-2021-XA-005, NYKY-2022-XA-004).

Data Availability Statement

Data available on request from the corresponding author.

Acknowledgments

Special thanks to the editor and reviewers for their valuable comments. In addition, the authors gratefully acknowledge Jun Wang and every teacher, classmate, and friend who helped the authors with their experiment and writing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sarker, J.R.; Singh, B.P.; Cowie, A.L.; Fang, Y.; Collins, D.; Badgery, W.; Dalal, R.C. Agricultural management practices impacted carbon and nutrient concentrations in soil aggregates, with minimal influence on aggregate stability and total carbon and nutrient stocks in contrasting soils. Soil Tillage Res. 2018, 178, 209–223. [Google Scholar] [CrossRef]
  2. Okolo, C.C.; Gebresamuel, G.; Zenebe, A.; Haile, M.; Eze, P.N. Accumulation of organic carbon in various soil aggregate sizes under different land use systems in a semi-arid environment. Agric. Ecosyst. Environ. 2020, 297, 106924. [Google Scholar] [CrossRef]
  3. Zhang, W.; Munkholm, L.J.; Liu, X.; An, T.; Xu, Y.; Ge, Z.; Xie, N.; Li, A.; Dong, Y.; Peng, C.; et al. Soil aggregate microstructure and microbial community structure mediate soil organic carbon accumulation: Evidence from one-year field experiment. Geoderma 2023, 430, 116324. [Google Scholar] [CrossRef]
  4. Feng, M.; Li, T.; He, B.; Bai, L.; Wu, X. Soil water repellency and micro-aggregate fractions in response to crop growth periods in a sloping cropland subjecting to long-term fertilization management. Sci. Total Environ. 2022, 838, 156350. [Google Scholar] [CrossRef]
  5. Rodrigues, L.A.T.; Giacomini, S.J.; Aita, C.; Lourenzi, C.R.; Brunetto, G.; Bacca, A.; Ceretta, C.A. Short-and long-term effects of animal manures and mineral fertilizer on carbon stocks in subtropical soil under no-tillage. Geoderma 2021, 386, 114913. [Google Scholar] [CrossRef]
  6. Pan, G.; Smith, P.; Pan, W. The role of soil organic matter in maintaining the productivity and yield stability of cereals in China. Agric. Ecosyst. Environ. 2009, 129, 344–348. [Google Scholar] [CrossRef]
  7. Sainju, U.M.; Caesar-TonThat, T.; Jabro, J.D. Carbon and nitrogen fractions in dryland soil aggregates affected by long-term tillage and cropping sequence. Soil Sci. Soc. Am. J. 2009, 73, 1488–1495. [Google Scholar] [CrossRef]
  8. Rabot, E.; Wiesmeier, M.; Schlüter, S.; Vogel, H.J. Soil structure as an indicator of soil functions: A review. Geoderma 2018, 314, 122–137. [Google Scholar] [CrossRef]
  9. Blanco-Canqui, H.; Lal, R. Mechanisms of carbon sequestration in soil aggregates. Crit. Rev. Plant Sci. 2004, 23, 481–504. [Google Scholar] [CrossRef]
  10. Jin, V.L.; Wienhold, B.J.; Mikha, M.M.; Schmer, M.R. Cropping system partially offsets tillage-related degradation of soil organic carbon and aggregate properties in a 30-yr rainfed agroecosystem. Soil Tillage Res. 2021, 209, 104968. [Google Scholar] [CrossRef]
  11. Fu, X.; Wang, J.; Sainju, U.M.; Liu, W. Aggregate size distribution and associated carbon and nitrogen in mulched winter wheat and spring corn. Can. J. Soil Sci. 2019, 99, 367–379. [Google Scholar] [CrossRef]
  12. Elliott, E.T. Aggregate structure and carbon, nitrogen, and phosphorus in native and cultivated soils. Soil Sci. Soc. Am. J. 1986, 50, 627–633. [Google Scholar] [CrossRef]
  13. Xue, B.; Huang, L.; Li, X.; Lu, J.; Gao, R.; Kamran, M. Straw residue incorporation and potassium fertilization enhances soil aggregate stability by altering soil content of iron oxide and organic carbon in a rice–rape cropping system. Land Degrad. Dev. 2022, 33, 2567–2584. [Google Scholar] [CrossRef]
  14. Udom, B.E.; Nuga, B.O.; Adesodun, J.K. Water-stable aggregates and aggregate-associated organic carbon and nitrogen after three annual applications of poultry manure and spent mushroom wastes. Appl. Soil Ecol. 2016, 101, 5–10. [Google Scholar] [CrossRef]
  15. Zhang, H.; Ding, W.; He, X.; Yu, H.; Fan, J.; Liu, D. Influence of 20-year organic and inorganic fertilization on organic carbon accumulation and microbial community structure of aggregates in an intensively cultivated sandy loam soil. PLoS ONE 2014, 9, e92733. [Google Scholar] [CrossRef]
  16. Dhaliwal, S.S.; Naresh, R.K.; Gupta, R.K.; Panwar, A.S.; Mahajan, N.C.; Singh, R.; Mandal, A. Effect of tillage and straw return on carbon footprints, soil organic carbon fractions and soil microbial community in different textured soils under rice–wheat rotation: A review. Rev. Environ. Sci. Bio/Technol. 2020, 19, 103–115. [Google Scholar] [CrossRef]
  17. Zhao, Z.; Gao, S.; Lu, C.; Li, X.; Li, F.; Wang, T. Effects of different tillage and fertilization management practices on soil organic carbon and aggregates under the rice–wheat rotation system. Soil Tillage Res. 2021, 212, 105071. [Google Scholar] [CrossRef]
  18. Six, J.; Paustian, K.; Elliott, E.T.; Combrink, C. Soil structure and organic matter I. Distribution of aggregate-size classes and aggregate-associated carbon. Soil Sci. Soc. Am. J. 2020, 64, 681–689. [Google Scholar] [CrossRef]
  19. Mikha, M.M.; Rice, C.W. Tillage and manure effects on soil and aggregate-associated carbon and nitrogen. Soil Sci. Soc. Am. J. 2004, 68, 809–816. [Google Scholar] [CrossRef]
  20. Yang, C.D.; Lu, S.G. Effects of five different biochars on aggregation, water retention and mechanical properties of paddy soil: A field experiment of three-season crops. Soil Tillage Res. 2021, 205, 104798. [Google Scholar] [CrossRef]
  21. Huang, R.; Tian, D.; Liu, J.; Lv, S.; He, X.; Gao, M. Responses of soil carbon pool and soil aggregates associated organic carbon to straw and straw-derived biochar addition in a dryland cropping mesocosm system. Agric. Ecosyst. Environ. 2018, 265, 576–586. [Google Scholar] [CrossRef]
  22. Galantini, J.A.; Senesi, N.; Brunetti, G.; Rosell, R. Influence of texture on organic matter distribution and quality and nitrogen and sulphur status in semiarid Pampean grassland soils of Argentina. Geoderma 2004, 123, 143–152. [Google Scholar] [CrossRef]
  23. Yilmaz, E.; Sönmez, M. The role of organic/bio–fertilizer amendment on aggregate stability and organic carbon content in different aggregate scales. Soil Tillage Res. 2017, 168, 118–124. [Google Scholar] [CrossRef]
  24. Li, T.; Zhang, Y.; Bei, S.; Li, X.; Reinsch, S.; Zhang, H.; Zhang, J. Contrasting impacts of manure and inorganic fertilizer applications for nine years on soil organic carbon and its labile fractions in bulk soil and soil aggregates. Catena 2020, 194, 104739. [Google Scholar] [CrossRef]
  25. Chen, Z.; Ti, J.S.; Chen, F. Soil aggregates response to tillage and residue management in a double paddy rice soil of the Southern China. Nutr. Cycl. Agroecosyst. 2017, 109, 103–114. [Google Scholar] [CrossRef]
  26. Mondal, S.; Naik, S.K.; Haris, A.A.; Mishra, J.S.; Mukherjee, J.; Rao, K.K.; Bhatt, B.P. Effect of conservation tillage and rice-based cropping systems on soil aggregation characteristics and carbon dynamics in Eastern Indo-Gangetic Plain. Paddy Water Environ. 2020, 18, 573–586. [Google Scholar] [CrossRef]
  27. Wang, J.; Dai, W.; Fang, K.; Gao, H.; Sha, Z.; Cao, L. Nutrient Characterization in Soil Aggregate Fractions with Different Fertilizer Treatments in Greenhouse Vegetable Cultivation. Agriculture 2022, 12, 440. [Google Scholar] [CrossRef]
  28. Duan, Y.; Yang, H.; Shi, T.; Zhang, W.; Xu, M.; Gao, S. Long-term manure application to improve soil macroaggregates and plant-available nitrogen in a Mollisol. Soil Tillage Res. 2021, 211, 105035. [Google Scholar] [CrossRef]
  29. Das, B.; Chakraborty, D.; Singh, V.K.; Aggarwal, P.; Singh, R.; Dwivedi, B.S.; Mishra, R.P. Effect of integrated nutrient management practice on soil aggregate properties, its stability and aggregate-associated carbon content in an intensive rice–wheat system. Soil Tillage Res. 2014, 136, 9–18. [Google Scholar] [CrossRef]
  30. Hati, K.M.; Jha, P.; Dalal, R.C.; Jayaraman, S.; Dang, Y.P.; Kopittke, P.M.; Kirchhof, G.; Menzies, N.W. 50 years of continuous no-tillage, stubble retention and nitrogen fertilization enhanced macro-aggregate formation and stabilisation in a Vertisol. Soil Tillage Res. 2021, 214, 105163. [Google Scholar] [CrossRef]
  31. Guo, Z.C.; Zhang, Z.B.; Zhou, H.; Rahman, M.T.; Wang, D.Z.; Guo, X.S.; Li, L.J.; Peng, X.H. Long-term animal manure application promoted biological binding agents but not soil aggregation in a Vertisol. Soil Tillage Res. 2018, 180, 232–237. [Google Scholar] [CrossRef]
  32. Wang, J.; Fu, X.; Ghimire, R.; Sainju, U.M.; Jia, Y.; Zhao, F. Responses of soil bacterial community and enzyme activity to organic matter components under long-term fertilization on the loess plateau of China. Appl. Soil Ecol. 2021, 166, 103992. [Google Scholar] [CrossRef]
  33. Xie, J.; Hou, M.; Zhou, Y.; Wang, R.; Zhang, S.; Yang, X.; Sun, B. Carbon sequestration and mineralization of aggregate-associated carbon in an intensively cultivated Anthrosol in north China as affected by long term fertilization. Geoderma 2017, 296, 1–9. [Google Scholar] [CrossRef]
  34. Wang, Y.; Ji, H.; Wang, R.; Guo, S.; Gao, C. Impact of root diversity upon coupling between soil C and N accumulation and bacterial community dynamics and activity: Result of a 30-year rotation experiment. Geoderma 2017, 292, 87–95. [Google Scholar] [CrossRef]
  35. Soil Survey Staff. Keys to Soil Taxonomy, 11th ed.; USDA Natural Resources Conservation Service: Washington, DC, USA, 2010.
  36. Lal, R.; Shukla, M.J. Principles of Soil Physics; Marcel Dekker, Inc.: New York, NY, USA, 2004; pp. 93–145. [Google Scholar]
  37. Cambardella, C.A.; Elliott, E.T. Particulate soil organic-matter changes across a grassland cultivation sequence. Soil Sci. Soc. Am. J. 1992, 56, 777–783. [Google Scholar] [CrossRef]
  38. Brookes, P.C.; Powlson, D.S.; Jenkinson, D.S. Measurement of microbial biomass phosphorus in soil. Soil Biol. Biochem. 1982, 14, 319–329. [Google Scholar] [CrossRef]
  39. Cao, H.; Xie, J.; Ciren, Q.; Guo, L.; Hong, J.; Jing, Y.; Meng, H. Effects of fertilization regimes on carbon and nitrogen contents of aggregates and maize yield in reclaimed soils. Trans. Chin. Soc. Agric. Eng. 2020, 36, 135–143. [Google Scholar]
  40. Bastian, M.; Heymann, S.; Jacomy, M. Gephi: An open source software for exploring and manipulating networks. In Proceedings of the Third International AAAI Conference on Weblogs and Social Media, San Jose, CA, USA, 17–20 May 2009; AAAI Publications: San Jose, CA, USA, 2009; pp. 361–362. [Google Scholar] [CrossRef]
  41. Zhao, B.; Chen, J.; Zhang, J.; Qin, S. Soil microbial biomass and activity response to repeated drying–rewetting cycles along a soil fertility gradient modified by long-term fertilization management practices. Geoderma 2010, 160, 218–224. [Google Scholar] [CrossRef]
  42. 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]
  43. Zhou, H.; Peng, X.; Peth, S.; Xiao, T.Q. Effects of vegetation restoration on soil aggregate microstructure quantified with synchrotron-based micro-computed tomography. Soil Tillage Res. 2012, 124, 17–23. [Google Scholar] [CrossRef]
  44. Mustafa, A.; Hu, X.; Shah, S.A.A.; Abrar, M.M.; Maitlo, A.A.; Kubar, K.A.; Saeed, Q.; Kamran, M.; Naveed, M.; Boren, W.; et al. Long-term fertilization alters chemical composition and stability of aggregate-associated organic carbon in a Chinese red soil: Evidence from aggregate fractionation, C mineralization, and 13C NMR analyses. J. Soils Sediments 2021, 21, 2483–2496. [Google Scholar] [CrossRef]
  45. Álvaro-Fuentes, J.; López, M.V.; Cantero-Martínez, C.; Arrúe, J.L. Tillage effects on soil organic carbon fractions in Mediterranean dryland agroecosystems. Soil Sci. Soc. Am. J. 2008, 72, 541–547. [Google Scholar] [CrossRef]
  46. Li, J.; Wu, X.; Gebremikael, M.T.; Wu, H.; Cai, D.; Wang, B.; Li, B.; Zhang, J.; Li, Y.; Xi, J. Response of soil organic carbon fractions, microbial community composition and carbon mineralization to high-input fertilizer practices under an intensive agricultural system. PLoS ONE 2018, 13, e0195144. [Google Scholar] [CrossRef] [PubMed]
  47. Yang, Z.; Singh, B.R.; Hansen, S. Aggregate associated carbon, nitrogen and sulfur and their ratios in long-term fertilized soils. Soil Tillage Res. 2007, 95, 161–171. [Google Scholar] [CrossRef]
  48. Bandyopadhyay, P.K.; Saha, S.; Mani, P.K.; Mandal, B. Effect of organic inputs on aggregate associated organic carbon concentration under long-term rice–wheat cropping system. Geoderma 2010, 154, 379–386. [Google Scholar] [CrossRef]
  49. Ding, X.; Han, X.; Liang, Y.; Qiao, Y.; Li, L.; Li, N. Changes in soil organic carbon pools after 10 years of continuous manuring combined with chemical fertilizer in a Mollisol in China. Soil Tillage Res. 2012, 122, 36–41. [Google Scholar] [CrossRef]
  50. Six, J.; Elliott, E.T.; Paustian, K.; Doran, J.W. Aggregation and soil organic matter accumulation in cultivated and native grassland soils. Soil Sci. Soc. Am. J. 1998, 62, 1367–1377. [Google Scholar] [CrossRef]
  51. Green, V.S.; Cavigelli, M.A.; Dao, T.H.; Flanagan, D.C. Soil physical properties and aggregate-associated C, N, and P distributions in organic and conventional cropping systems. Soil Sci. 2005, 170, 822–831. [Google Scholar] [CrossRef]
  52. Zhang, Y.; Shengzhe, E.; Wang, Y.; Su, S.; Bai, L.; Wu, C.; Zeng, X. Long-term manure application enhances the stability of aggregates and aggregate-associated carbon by regulating soil physicochemical characteristics. Catena 2021, 203, 105342. [Google Scholar] [CrossRef]
  53. Totsche, K.U.; Amelung, W.; Gerzabek, M.H.; Guggenberger, G.; Klumpp, E.; Knief, C.; Lehndorff, E.; Mikutta, R.; Peth, S.; Prechtel, A.; et al. Microaggregates in soils. J. Plant Nutr. Soil Sci. 2018, 181, 104–136. [Google Scholar] [CrossRef]
  54. Liao, H.; Zhang, Y.; Wang, K.; Hao, X.; Chen, W.; Huang, Q. Complexity of bacterial and fungal network increases with soil aggregate size in an agricultural Inceptisol. Appl. Soil Ecol. 2020, 154, 103640. [Google Scholar] [CrossRef]
  55. Banger, K.; Kukal, S.S.; Toor, G.; Sudhir, K.; Hanumanthraju, T.H. Impact of long-term additions of chemical fertilizers and farm yard manure on carbon and nitrogen sequestration under rice-cowpea cropping system in semi-arid tropics. Plant Soil 2009, 318, 27–35. [Google Scholar] [CrossRef]
  56. Wang, L.M.; Huang, D.F.; Fang, Y.; Wang, F.; Li, F.L.; Liao, M. Soil fungal communities in tea plantation after 10 years of chemical vs. integrated fertilization. Chil. J. Agric. Res. 2017, 77, 355–364. [Google Scholar] [CrossRef]
  57. Li, W.T.; Liu, M.; Jiang, C.Y.; Wu, M.; Chen, X.F.; Ma, X.Y.; Li, Z.P. Changes in soil aggregate-associated enzyme activities and nutrients under long-term chemical fertilizer applications in a phosphorus-limited paddy soil. Soil Use Manag. 2017, 33, 25–33. [Google Scholar] [CrossRef]
  58. Li, S.X.; Wang, Z.H.; Miao, Y.F.; Li, S.Q. Soil organic nitrogen and its contribution to crop production. J. Integr. Agric. 2014, 13, 2061–2080. [Google Scholar] [CrossRef]
  59. Fan, J.; Xiao, J.; Liu, D.; Ye, G.; Luo, J.; Houlbrooke, D.; Laurenson, S.; Yan, J.; Chen, L.; Tian, J.; et al. Effect of application of dairy manure, effluent and inorganic fertilizer on nitrogen leaching in clayey fluvo-aquic soil: A lysimeter study. Sci. Total Environ. 2017, 592, 206–214. [Google Scholar] [CrossRef]
  60. Weil, R.R.; Kroontje, W. Physical condition of a Davidson clay loam after five years of heavy poultry manure applications. J. Environ. Qual. 1979, 8, 387–392. [Google Scholar] [CrossRef]
Figure 1. Soil aggregate size distribution under different fertilization treatments. Error bars indicate ±standard deviation. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Different letters above error bars indicate significant difference between treatments (p < 0.05, LSD test).
Figure 1. Soil aggregate size distribution under different fertilization treatments. Error bars indicate ±standard deviation. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Different letters above error bars indicate significant difference between treatments (p < 0.05, LSD test).
Agronomy 13 01466 g001
Figure 2. Soil organic C (SOC), particulate organic C (POC), microbial biomass C (MBC), soil total N (STN), particulate organic N (PON), and microbial biomass N (MBN) in different aggregates as affected by chemical fertilization and organic manure. The subfigures (af) represent SOC, POC, MBC, STN, PON, and MBN, respectively. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Error bars indicate ±standard deviation. Different letters above error bars indicate significant difference between treatments (p < 0.05, LSD test).
Figure 2. Soil organic C (SOC), particulate organic C (POC), microbial biomass C (MBC), soil total N (STN), particulate organic N (PON), and microbial biomass N (MBN) in different aggregates as affected by chemical fertilization and organic manure. The subfigures (af) represent SOC, POC, MBC, STN, PON, and MBN, respectively. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Error bars indicate ±standard deviation. Different letters above error bars indicate significant difference between treatments (p < 0.05, LSD test).
Agronomy 13 01466 g002aAgronomy 13 01466 g002b
Figure 3. Correlation heatmap between soil’s basic properties and aggregate parameters. The intensity of color indicates the correlation coefficient (R) (blue and red indicate a positive and negative correlation, respectively). *, **, and *** indicate p < 0.05, 0.01, and 0.001, respectively.
Figure 3. Correlation heatmap between soil’s basic properties and aggregate parameters. The intensity of color indicates the correlation coefficient (R) (blue and red indicate a positive and negative correlation, respectively). *, **, and *** indicate p < 0.05, 0.01, and 0.001, respectively.
Agronomy 13 01466 g003
Figure 4. Co-occurrence networks of soil aggregate C/N fractions and the aggregate parameters of bulk soil (BS) under different fertilization based on correlation analysis. A red edge indicates a positive relationship, while a green edge indicates a negative relationship between two individual nodes. An orange node indicates the aggregate size and stability. A green node indicates the C/N fractions of bulk soil. A purple node indicates the soil aggregate C/N fractions.
Figure 4. Co-occurrence networks of soil aggregate C/N fractions and the aggregate parameters of bulk soil (BS) under different fertilization based on correlation analysis. A red edge indicates a positive relationship, while a green edge indicates a negative relationship between two individual nodes. An orange node indicates the aggregate size and stability. A green node indicates the C/N fractions of bulk soil. A purple node indicates the soil aggregate C/N fractions.
Agronomy 13 01466 g004
Figure 5. The average wheat yield across 34 years after long-term fertilization. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Different letters above error bars indicate significant difference between treatments (p < 0.05, LSD test).
Figure 5. The average wheat yield across 34 years after long-term fertilization. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Different letters above error bars indicate significant difference between treatments (p < 0.05, LSD test).
Agronomy 13 01466 g005
Table 1. Soil’s basic properties under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
Table 1. Soil’s basic properties under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
TreatmentsBD (g cm−3)Porosity (%)WHC (%)SOC (g kg−1)STN (g kg−1)C/N Ratio
CK1.19 ± 0.04 b55.00 ± 1.57 b20.25 ± 0.61 f9.29 ± 0.36 ef0.98 ± 0.07 d11.09 ± 0.51 a
NF1.22 ± 0.11 ab54.04 ± 4.02 bc24.52 ± 0.72 e10.11 ± 0.55 de1.39 ± 0.06 b7.27 ± 0.09 cd
PF1.26 ± 0.11 a52.62 ± 4.17 c25.62 ± 0.25 de10.43 ± 0.37 cd1.02 ± 0.03 cd9.10 ± 0.06 b
NP1.15 ± 0.04 b56.64 ± 1.64 ab26.30 ± 0.41 cd10.85 ± 0.32 cd1.15 ± 0.07 c10.90 ± 0.02 a
M1.13 ± 0.06 bc57.23 ± 2.35 ab27.54 ± 0.61 bc12.55 ± 0.71 b1.12 ± 0.06 cd9.33 ± 0.15 b
NM1.12 ± 0.08 bc57.81 ± 2.93 ab28.19 ± 0.98 b11.45 ± 0.35 c1.64 ± 0.09 a6.98 ± 0.17 d
PM1.10 ± 0.07 c58.62 ± 2.56 a27.84 ± 0.69 b11.07 ± 0.69 cd1.43 ± 0.07 b7.74 ± 0.13 c
NPM1.14 ± 0.05 bc57.00 ± 1.83 ab29.85 ± 0.45 a14.27 ± 0.18 a1.58 ± 0.06 a9.05 ± 0.21 b
Table 2. Soil aggregate stability is measured as the content of aggregates larger than 0.25 mm (R0.25), mean weight diameter (MWD), geometric mean diameter (GMD), and the fractal dimension (D) under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
Table 2. Soil aggregate stability is measured as the content of aggregates larger than 0.25 mm (R0.25), mean weight diameter (MWD), geometric mean diameter (GMD), and the fractal dimension (D) under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
TreatmentsR0.25 (%)MWD (mm)GMD (mm)D
CK82.58 ± 1.25 d2.27 ± 0.23 e1.03 ± 0.04 e2.65 ± 0.03 a
NF79.09 ± 1.17 e2.35 ± 0.14 e1.01 ± 0.03 e2.66 ± 0.02 a
PF83.58 ± 2.53 d2.32 ± 0.19 e1.04 ± 0.05 e2.64 ± 0.03 a
NP89.01 ± 1.50 b2.82 ± 0.11 c1.18 ± 0.03 c2.55 ± 0.02 c
M86.32 ± 1.76 c2.53 ± 0.13 d1.11 ± 0.03 d2.60 ± 0.02 a
NM88.47 ± 0.88 b2.95 ± 0.17 c1.21 ± 0.04 bc2.54 ± 0.02 cd
PM87.58 ± 1.68 bc3.17 ± 0.13 b1.24 ± 0.04 b2.52 ± 0.03 d
NPM91.10 ± 0.95 a3.42 ± 0.25 a1.33 ± 0.07 a2.46 ± 0.05 e
Table 3. Soil C/N ratio in aggregates under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
Table 3. Soil C/N ratio in aggregates under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
TreatmentsMacro-Aggregates
(>2 mm)
Medium-Aggregates
(2–0.25 mm)
Micro-Aggregates
(<0.25 mm)
CK10.76 ± 0.07 a10.33 ± 0.45 a11.13 ± 0.92 a
NF10.22 ± 0.18 a10.04 ± 0.48 a10.18 ± 0.42 a
PF10.80 ± 0.84 a10.30 ± 0.65 a10.37 ± 0.59 a
NP7.68 ± 0.43 bc7.56 ± 0.18 bc7.48 ± 0.46 b
M8.01 ± 0.42 b7.84 ± 0.35 b7.85 ± 0.36 b
NM6.98 ± 0.28 bc6.99 ± 0.18 bc6.94 ± 0.30 b
PM7.82 ± 0.48 bc7.66 ± 0.49 bc7.59 ± 0.18 b
NPM6.81 ± 0.37 c6.76 ± 0.19 c6.62 ± 0.37 b
Table 4. The contribution rate (CR) and enrichment coefficient (EC) of soil aggregate organic C and total N under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
Table 4. The contribution rate (CR) and enrichment coefficient (EC) of soil aggregate organic C and total N under different fertilization treatments. Treatments were CK, unfertilized control; NF, inorganic N fertilizer; PF, inorganic P fertilizer; NP, inorganic N plus P fertilizer; M, organic manure; NM, inorganic N fertilizer plus manure; PM, inorganic P fertilizer plus manure; and NPM, inorganic N plus P fertilizers plus manure. Values are means ± standard errors (n = 3). Different letters indicate significant difference between treatments (p < 0.05, LSD test).
ItemTreatmentsMacro-Aggregates (>2 mm)Medium-Aggregates (2–0.25 mm)Micro-Aggregates (<0.25 mm)
CR (%)ECCR (%)ECCR (%)EC
SOCCK32.04 ± 3.10 c0.99 ± 0.02 b47.17 ± 2.92 bc0.99 ± 0.02 b17.42 ± 1.74 b0.94 ± 0.02 cd
NF35.79 ± 3.69 bc1.07 ± 0.05 b46.74 ± 2.69 bc1.07 ± 0.05 b20.83 ± 1.46 a1.02 ± 0.05 b
PF38.65 ± 5.05 bc1.16 ± 0.06 a55.92 ± 1.75 a1.16 ± 0.06 a16.92 ± 1.63 bc1.11 ± 0.04 a
NP35.73 ± 1.88 bc0.86 ± 0.05 c39.50 ± 2.29 d0.86 ± 0.05 c7.98 ± 0.82 fg0.83 ± 0.05 e
M38.92 ± 0.86 bc1.06 ± 0.03 b50.52 ± 2.92 ab1.06 ± 0.03 b14.22 ± 1.25 cd1.02 ± 0.02 bc
NM42.48 ± 3.30 ab0.98 ± 0.02 b41.73 ± 2.31 cd0.98 ± 0.02 b10.62 ± 0.19 ef0.93 ± 0.02 d
PM47.28 ± 1.93 a1.01 ± 0.05 b39.32 ± 3.36 d1.01 ± 0.05 b11.74 ± 1.75 de0.96 ± 0.05 c
NPM41.12 ± 3.70 ab0.79 ± 0.01 c29.25 ± 3.81 e0.79 ± 0.01 c6.71 ± 1.00 g0.75 ± 0.01 f
STNCK33.03 ± 3.62 d1.02 ± 0.06 bc53.96 ± 3.37 b1.07 ± 0.03 b16.31 ± 0.31 a0.94 ± 0.01 d
NF25.41 ± 2.17 e0.76 ± 0.02 d35.84 ± 1.50 d0.78 ± 0.01 c15.26 ± 0.64 ab0.73 ± 0.00 e
PF32.43 ± 1.49 d0.98 ± 0.03 c52.54 ± 2.06 b1.04 ± 0.03 b16.06 ± 2.75 a0.97 ± 0.03 c
NP50.69 ± 0.17 ab1.23 ± 0.01 a60.10 ± 1.77 ab1.26 ± 0.03 a13.26 ± 0.54 bc1.21 ± 0.01 a
M45.45 ± 1.26 bc1.24 ± 0.01 a63.62 ± 2.32 a1.28 ± 0.01 a16.52 ± 1.37 a1.21 ± 0.00 a
NM42.49 ± 2.76 c0.98 ± 0.01 c44.16 ± 3.21 c0.98 ± 0.07 b10.78 ± 0.65 cd0.94 ± 0.01 d
PM46.90 ± 1.77 bc1.00 ± 0.02 bc42.23 ± 5.69 cd1.03 ± 0.11 b12.18 ± 0.99 c0.98 ± 0.02 c
NPM54.86 ± 6.50 a1.05 ± 0.02 b41.52 ± 5.22 cd1.06 ± 0.00 b9.11 ± 0.12 d1.02 ± 0.03 b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, C.; Sainju, U.M.; Li, C.; Fu, X.; Zhao, F.; Wang, J. Long-Term Chemical and Organic Fertilization Differently Affect Soil Aggregates and Associated Carbon and Nitrogen in the Loess Plateau of China. Agronomy 2023, 13, 1466. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13061466

AMA Style

Yang C, Sainju UM, Li C, Fu X, Zhao F, Wang J. Long-Term Chemical and Organic Fertilization Differently Affect Soil Aggregates and Associated Carbon and Nitrogen in the Loess Plateau of China. Agronomy. 2023; 13(6):1466. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13061466

Chicago/Turabian Style

Yang, Caidi, Upendra M. Sainju, Chao Li, Xin Fu, Fazhu Zhao, and Jun Wang. 2023. "Long-Term Chemical and Organic Fertilization Differently Affect Soil Aggregates and Associated Carbon and Nitrogen in the Loess Plateau of China" Agronomy 13, no. 6: 1466. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13061466

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