Grain yield
In the early rice season, rice yield initially increased and then decreased with increasing substitution of N fertilizer. The yield of the N60G40 treatment is the highest, showing a 2.6% increase compared to the N100 treatment. However, when the substitution rate reached 80%, there was a significant decrease in crop yield. In the late rice season, yield gradually decreased with increasing N fertilizer substitution, although none of the reductions reached a significant level compared to the N100 treatment.
Table 2
Rice grain yield in different fertilization treatments
Treatment | Yield |
Early rice | Late rice | Annual |
CK | 2490 ± 161a) cb) | 4593 ± 233 b | 7083 ± 359 c |
N100 | 5027 ± 201 a | 7693 ± 291 a | 12720 ± 497 a |
N80G20 | 5047 ± 183 a | 7413 ± 462 a | 12460 ± 639 ab |
N60G40 | 5157 ± 175 a | 7077 ± 176 a | 12233 ± 130 ab |
N40G60 | 4723 ± 57 ab | 6934 ± 369 a | 11660 ± 346 ab |
N20G80 | 4407 ± 219 b | 6797 ± 181 a | 11203 ± 394 b |
-
Different lowercase letters represent the correlation was significant (P < 0.05).
-
Means ± standard errors (n = 3)
Distribution of soil aggregates and their stability indices
Under all treatments substituting partial N fertilizer with MV, the proportions of > 2mm and 1-2mm water-stable aggregates were the highest (68%-72%) and the lowest (4.2%-4.4%), respectively. The incorporation of MV increased the content of > 2mm water-stable aggregates across all treatments. The N60G40 treatment exhibited the highest content (72.41%), but as the amount of MV increased, there was a decreasing trend in the content of > 2mm water-stable aggregates. Analysis of the proportion of water-stable aggregates smaller than < 0.25mm showed similar trends for aggregates between 0.053-0.25mm and those smaller than 0.053mm across different treatments. In both cases, the N100 treatment had the highest proportion, while the N60G40 treatment had the lowest. Compared with the N100 treatment, the N60G40 treatment reduced the contents of 0.053-0.25mm and < 0.053mm water-stable aggregates by 15.44% and 18.67%, respectively. The subsequent reductions were observed in the N40G60 (13.04% and 10.00%) and N80G20 (9.60% and 9.24%) treatments. The changes in MWD and GMD were consistent across all treatments (Fig. 1c, Fig. 1d). Long-term fertilization reduced soil aggregate stability. Compared with the CK, the MWD and GMD of the N100 treatment decreased by 2.19% and 3.27%, respectively. Different substitution rates of MV improved soil aggregate stability. Compared with the N100 treatment, N80G20 increased MWD and GMD by 4.19% and 7.73%, respectively, while N60G40 increased MWD and GMD by 5.31% and 12.08%, respectively.
Figure 1 Effects of fertilizer-N substituted by MV on different proportions of soil aggregate (a, b), MWD(c) and GMD(d).
SOC density fractions
Incorporating MV positively impacted the content of soil active organic carbon (Fig. 2a). The N60G40 treatment showed the most significant improvement, followed by N40G60. The order of LOC content across various treatments was N60G40 > N40G60 > N80G20 > N20G80 > N100 > CK. Compared with the N100 treatment, N40G60 mainly increased the content of HLOC by 24.34%, while N60G40 primarily increased the content of MLOC by 16.40%. As shown in Fig. 2b, all treatments significantly increased soil AI, CPI, and CPMI, with increased ranging from 7.77–30.10%, 5.55–15.74%, and 18.91–44.96%, respectively. As the proportion of N fertilizer substitution increased, CPI and CPMI initially increased and then decreased. The N60G40 treatment showed the most significant improvement, with CPI and CPMI increasing by 15.74% and 44.96%, respectively, compared with the N100 treatment,
Figure 2 Effects of fertilizer-N substituted by MV on organic carbon content of whole soil.
Soil aggregate nutrients and its distribution
Substituting CF with MV increased the content of soil organic matter (SOM), total nitrogen (TN), total phosphorus (TP), and total potassium (TK) in the 1-2mm aggregate, with total nutrient content increasing alongside higher MV substitution rates (Fig. 3a). The SOM content initially increased and then decreased with decreasing aggregate particle size, peaking in the 1-2mm aggregate, followed by the > 2mm aggregate. Apart from the 1-2mm aggregate, the N80G20 treatment exhibited the highest SOM content in other particle size aggregates, ranging from 17.75–27.09 g kg− 1. The TN and TP contents of soil aggregate gradually decreased with decreasing particle size and were mainly concentrated in macro-aggregates (> 0.25mm). The highest TN content was observed in the > 2mm aggregate under the N80G20 treatment, at 3.19 g kg− 1. Compared with the N80G20 treatment, the N40G60 treatment significantly reduced the TN content in the 0.053-0.25mm and < 0.053mm aggregates. The N80G20 treatment also significantly increased TP content in the 1-2mm aggregate compared with the N100 treatment. Regarding TK content, the N60G40 treatment increased the content in the > 2mm and 1-2mm aggregates compared with the N100 treatment, although the increase was not statistically significant.
The alkali-hydrolysable N (AN) and available P (AP) in soil aggregates gradually decreased with decreasing aggregate particle size (Fig. 3b). Compared with the N100 treatment, the N80G20 treatment increased AN content in the 1-2mm and 0.25-0.5mm aggregates by 7.83% and 3.27%, respectively. Similarly, the N60G40 treatment followed this trend. Under different proportions of MV, the trend in AP content across various aggregate size was consistent with AN. The N80G20 treatment significantly increased AP content in the 1-2mm and 0.5-1mm aggregates by 29.23% and 27.79%, respectively, compared to the N100 treatment. Incorporating MV reduced the available potassium (AK) in various particle aggregates, except for the 1-2mm size, with the reduction increasing alongside higher MV substitution rates. Compared with the N100 treatment, the N20G80 treatment reduced the AK content in the 0.5-1mm aggregate significantly.
Under different substitution rates of N fertilizer, >2mm aggregates contributed more to SOM and TN contents, with their contribution rates initially increasing and then decreasing as the MV proportion increased (Fig. 4). Compared with the N100 treatment, all treatments substituting N fertilizer with MV increased the contribution rates of > 2mm aggregate to SOM and TN. The N40G60 and N60G40 treatments showed the greatest increase, with SOM contribution rates increasing by 3.41% and 4.07%, and TN contribution rates increasing by 6.00% and 8.17%, respectively. Meanwhile, these treatments reduced the contribution rates of < 0.25mm aggregate, with SOM contribution rates decreasing by 15.77% and 15.65%, and TN contribution rates decreasing by 5.17% and 7.55%, respectively. Substituting N fertilizer with MV also increased the contribution rates of > 2mm aggregate to soil TP and TK, with the N60G40 treatment showing the highest increases of 9.31% and 5.63%, respectively, compared with the N100 treatment. Conversely, it reduced the contribution rates of < 0.25mm aggregate to soil TP and TK, with reductions of 34.66% and 17.96%, respectively, compared with the N100 treatment.
Figure 3 Effects of N fertilizer substituted by MV on aggregate-associated SOM, TN, TP and TK; AN, AP and AK contents
Figure 4 Effects of N fertilizer substituted by MV on SOM, TN, TP and TK contribution rates in soil aggregates
SQI and rice grain yield
Soil chemical properties. Incorporating MV significantly increased the soil AN, AP and AK contents (Table 3), compared with CK, the AN in N40G60 treatment and the AP content in N60G40 treatment improved by 21.08% and 72.33%, respectively. Substitution N fertilizer with MV significantly increased SOM and TN contents, with the N60G40 treatment showing the greatest improvement among all treatments. Compared to the N100 treatment, SOM and TN contents improved by 10.16% and 9.52%, respectively.
Table 3
Chemical properties of the whole soil.
Treatments
|
pH
|
SOM
|
TN
|
TP
|
TK
|
AN
|
AP
|
AK
|
CK
|
7.58 ± 0.08 b
|
38.4 ± 1.31 c
|
2.67 ± 0.01 b
|
0.99 ± 0.03 c
|
22.2 ± 0.55 ab
|
204 ± 3.0 d
|
11.9 ± 3.45 b
|
55.1 ± 1.22 b
|
|
N100
|
7.67 ± 0.03 ab
|
43.3 ± 2.50 b
|
2.73 ± 0.19 b
|
1.16 ± 0.06 ab
|
21.4 ± 0.47 bc
|
224 ± 6.8 c
|
38.0 ± 7.36 a
|
74.6 ± 4.94 a
|
|
N80G20
|
7.65 ± 0.02 ab
|
43.7 ± 1.52 b
|
2.58 ± 0.14 b
|
1.09 ± 0.06 b
|
21.3 ± 0.33 bc
|
231 ± 2.7 bc
|
40.0 ± 3.60 a
|
69.2 ± 4.09 a
|
|
N60G40
|
7.63 ± 0.01 ab
|
47.7 ± 1.09 a
|
2.99 ± 0.01 a
|
1.24 ± 0.01 a
|
21.0 ± 0.05 c
|
236 ± 0.5 ab
|
43.0 ± 3.30 a
|
71.6 ± 1.63 a
|
|
N40G60
|
7.71 ± 0.05 a
|
46.8 ± 0.71 ab
|
2.72 ± 0.09 b
|
1.14 ± 0.02 b
|
22.5 ± 0.20 a
|
247 ± 10.8 a
|
36.0 ± 4.23 a
|
72.2 ± 1.69 a
|
|
N20G80
|
7.73 ± 0.04 a
|
47.3 ± 0.27 a
|
2.79 ± 0.13 ab
|
1.15 ± 0.04 ab
|
21.8 ± 0.67 abc
|
242 ± 8.5 ab
|
36.8 ± 1.48 a
|
67.8 ± 3.78 a
|
|
The relationship between SQI, grain yield and soil chemical properties.
Substituting N fertilizer with MV increased the SQI (Fig. 5a). The SQI showed an initial increase followed by a decreasing trend, with the N60G40 treatment exhibiting the highest SQI, which was 19.08% higher compared with the N100 treatment. Additionally, the annual rice yield significantly increased with the improvement of SQI (Fig. 5b). Random forest model analysis integrated the prediction results of soil chemical properties for SQI (Fig. 5c) and rice yield (Fig. 5d). The results revealed that AN, TP, AP, SOM and TN were the primary predictors influencing SQI, while AK, MLOC, AP and SQI were key factors for interpreting rice yield
Figure 5 Comprehensive analysis of SQI, grain yield and soil chemical properties
The SQI across different treatments (a) and its correlation with annual rice grain yield (b), the relative influence (%) of soil chemical properties on SQI (c) and rice grain yield (d) based on a relative important model. (*: P < 0.05; **: P < 0.01)