Soil bulk density and porosity
The soil bulk density (BD) with mulching treatments were lower than pre-treatment (PS) value after 8-year continuous film mulching, average decreased by 3.85%, while CK was only decreased by 0.15% (Fig. 4). After maize harvest in 2020, the BD value of each treatment was ordered as CK > P > R > D > S, which was lowest in 20–30 cm, and highest in 0–10 cm depth. Compared with CK, the average BD under film mulching (P, R, D, and S) were significant decreased by 3.68% (P < 0.05). Only S was significantly lower than P and R by 1.83% and 1.65%, respectively, and no significant difference among other film treatments.
Contrary to BD, the soil porosity (SP) with mulching treatments were all higher than PS value, average increased by 3.27% in 0–40 cm depth, and the difference increase gradually with the deepening of the soil layer (Fig. 4). The order of the value in each soil depth was S > D > R > P > CK, the film mulching treatments were all significantly (P < 0.05) higher than CK, average increased by 7.76%, 5.70%, 1.72%, and 2.01% under S, D, R, and P, respectively, but no significant difference among the four film treatments.
Figure 4.
Soil mechanical and water stable aggregates
After the 8-year continuous plastic film mulching, only > 5 mm size fraction significant increased with the soil layer depth. And the proportion of the soil mechanical stable aggregates in < 0.25 mm size fraction was highest in 0–40 cm depth, then 1-0.5, 2 − 1 and 5 − 2 mm, and the 0.5 − 0.25 mm was the lowest. Compared with pre-treatment (PS), each mulching treatments was increased the macroaggregate (> 0.25 mm) content by 8.64% in average, and also was 15.68% (P < 0.05) higher than CK (Fig. 5). In 0–10 cm layer, the S was significantly increased each aggregate fraction-size than that with other mulching treatments (D, R, and P), the macroaggregate (> 0.25 mm) content average increased by 10.38% (P < 0.05). The difference between mulching treatments increased gradually as the soil layer deepened. In 10–20 cm, S or D treatment was significant (P < 0.05) increased the proportions of aggregates in the size fractions 5 − 2, 1-0.5 and 0.5 − 0.25 mm by 18.80%, 14.35%, 16.82%, respectively, compared to P treatment. Meanwhile, the S or D treatment increased (P < 0.05) the percentages of aggregates in the > 5 and 5–2 mm fractions to R or P treatment in 20–30 cm. In 30–40 cm, the S or D significant increased (P < 0.05) in the > 5 mm fraction, but decreased in other fractions by 8.10% or 15.97% in average to R or P treatment.
Figure 5.
The proportion of the soil water stable aggregates in < 0.25 mm size fraction was highest in 0–40 cm depth, and the five treatments increased the macroaggregate content (> 0.25 mm) by 34.63% to pre-treatment in average (Fig. 6). In 0–40 cm layer, every mulching treatment increased the proportions of macroaggregate by 24.67%, 17.99%, 11.85%, and 8.49% with S, D, R, and P, respectively, compared to CK. The proportions of aggregates in the size fractions under mulching treatments was ordered as 1-0.5 mm > 0.5 − 0.25 mm > 2 − 1 mm > 5 − 2 mm, however, the 0.5 − 0.25 mm size fraction was highest and 5 − 2 mm lowest under CK. The small variation of soil water stable aggregate fraction-size distribution was found between mulching treatments in different soil layers, except the 5 − 2 mm fraction significant increased (P < 0.05) by 16.10% and 7.20% in average in 0–10 and 20–30 cm depth under S or D compared to R or P treatment, respectively.
Figure 6.
Mean-weight diameter (MWD) and geometric mean diameter (GMD) of soil aggregates
The variation of each treatment on soil MWD was different after 8-year continuous mulching (Table 1). The MWD of the mechanical-stable aggregates (MMWD) was increased gradually with the layer depth. Compared with pre-treatment (PS), only S and D were increased by 16.04% and 11.76% in average, while CK was 12.83% decreased in 0–40 cm layer. In 0–40 cm layer, very mulching treatments significant (P < 0.05) increased the MMWD by 33.13%, 28.22%, 19.63%, and 9.20% with S, D, R, and P, respectively, compared to CK. Contrary to the MMWD, the MWD of the water-stable aggregates (WMWD) was decreased gradually with the layer depth, and the five treatments, which was ordered as S > D > R > P > CK in each depth in 0–40 cm, increased by 10.53% to pre-treatment in average (Table 1). Compared with CK, mulching treatments average increased by 6.25%, and nearly no variation of was found between mulching treatments in different soil layers, except S and D significant (P < 0.05) increased by 11.11% and 6.67% in average in 0–10 cm depth compared to P treatment, respectively.
Table 1
Mean weight diameter (MWD) of soil mechanical-stable and water-stable aggregates in 0–40 cm depth under different treatments after 8-year continuous film mulching.
Treatment | Mechanical-stable aggregates (MMWD) | | Water-stable aggregates (WMWD) |
0–10 cm | 10–20 cm | 20–30 cm | 30–40 cm | 0–40 cm AVG | | 0–10 cm | 10–20 cm | 20–30 cm | 30–40 cm | 0–40 cm AVG |
PS | 1.38 | 1.65 | 2.09 | 2.34 | 1.87 | | 0.39 | 0.37 | 0.39 | 0.38 | 0.38 |
S | 1.39a | 1.82a | 2.71a | 2.77a | 2.17a | | 0.50a | 0.44a | 0.42a | 0.39a | 0.44a |
D | 1.36b | 1.77a | 2.57a | 2.67a | 2.09a | | 0.48a | 0.43a | 0.42a | 0.39a | 0.43a |
R | 1.35b | 1.72b | 2.17b | 2.55a | 1.95a | | 0.47ab | 0.42ab | 0.41ab | 0.39ab | 0.42ab |
P | 1.31b | 1.66b | 1.98b | 2.16bc | 1.78bc | | 0.45b | 0.41ab | 0.41ab | 0.39ab | 0.41ab |
CK | 1.26c | 1.56c | 1.67c | 2.04c | 1.63c | | 0.42c | 0.40b | 0.39b | 0.38b | 0.40b |
Note: PS, pre-sowing; S, planting furrows separated by consecutive plastic film-mulched ridges; R, traditional ridge-furrow rainwater harvesting planting; D, double ridge-furrow with full plastic film mulching; P, alternate mulching and bare rows without ridges; CK, flat planting without mulching. Different lower-case letters indicate significant differences at P < 0.05. |
Table 1
The GMD of the mechanical-stable aggregates (MGMD) was far greater than that in water-stable aggregates (WGMD) after continuous 8-year mulching (Table 2). Compared with pre-treatment (PS), the MGMD under S, D, and R treatments was average increased by 28.26%, 21.74%, and 9.78%, respectively, however, CK was 13.04% lower. The MGMD under four mulching treatments was ordered as S > D > R > P in each depth in 0–40 cm, and average increased by 47.50%, 40.00%, 26.25%, and 11.25%, respectively, compared to CK treatment. Small variation was found between mulching treatments, and only S or D significantly higher than P in 20–30 or 30–40 cm depth by 48.7% and 35.40% in average. Different from MGMD, the WGMD of the five treatments were ordered as S > D > R > P > CK in each depth in 0–40 cm, increased by 7.59% to pre-treatment in average (Table 2). While each mulching treatments were all higher than CK, and average increased by 5.01%, however, no significant difference among four mulching treatments.
Table 2
Geometric mean diameter (GMD) of soil mechanical-stable and water-stable aggregates in 0–40 cm depth under different treatments after 8-year continuous film mulching.
Treatment | Mechanical-stable aggregates (MGMD) | | Water-stable aggregates (WGMD) |
0–10 cm | 10–20 cm | 20–30 cm | 30–40 cm | 0–40 cm AVG | | 0–10 cm | 10–20 cm | 20–30 cm | 30–40 cm | 0–40 cm AVG |
PS | 0.66 | 0.78 | 1.02 | 1.22 | 0.92 | | 0.30 | 0.29 | 0.29 | 0.29 | 0.29 |
S | 0.72a | 0.92a | 1.51a | 1.57a | 1.18a | | 0.34a | 0.32a | 0.31a | 0.30a | 0.32a |
D | 0.69ab | 0.88a | 1.40a | 1.49a | 1.12a | | 0.33a | 0.32a | 0.31a | 0.30a | 0.32a |
R | 0.67ab | 0.85a | 1.10b | 1.41a | 1.01ab | | 0.33ab | 0.31ab | 0.31ab | 0.30a | 0.31ab |
P | 0.65ab | 0.80ab | 0.98b | 1.13b | 0.89bc | | 0.32ab | 0.31ab | 0.31ab | 0.30a | 0.31ab |
CK | 0.61b | 0.72b | 0.82c | 1.04c | 0.80c | | 0.31b | 0.31b | 0.30b | 0.29a | 0.30b |
Note: PS, pre-sowing; S, planting furrows separated by consecutive plastic film-mulched ridges; R, traditional ridge-furrow rainwater harvesting planting; D, double ridge-furrow with full plastic film mulching; P, alternate mulching and bare rows without ridges; CK, flat planting without mulching. Different lower-case letters indicate significant differences at P < 0.05. |
Table 2
FD Fractal dimension (FD) of soil Aggregate
The fractal dimension (FD) with mechanical- or water-stable aggregates have a range of values from 2.64–2.80 or 2.91–2.96, and each mulching treatments decreased FD value in 0–40 cm layer, compared to PS (Fig. 7). The FD with mechanical-stable aggregates (MFD) of five treatments was ordered as CK > P > R > D > S, mulching treatments decreased by 2.10% in average than CK in 0–40 cm layers. The variation between four mulching treatments was gradually decreased with the depth in 0–40 cm layer, in 0–10 cm, the S was significantly lower than D, R, or P by 1.25%, 1.89%, and 2.46%, respectively. Contrary to the MFD, the FD with water-stable aggregates (WFD) was gradually increase with soil depth. Mulching treatments was ordered as P > R > D > S in 0–40 cm depth, which were lower than CK by 0.36% in average, and no significant differences were observed among the four mulching treatments.
Figure 7. Fractal dimension (FD) of soil Aggregate at 0–40 cm depth under per-treatment and different film mulching after 8-year continuous film mulching.
Percentage of aggregate destruction (PAD) and unstable aggregates index (ELT)
Continuous plastic film mulching decreased the percentage of aggregate destruction (PAD) in 0–40 cm layer by 6.39%, compared to pre-treatment (PS) (Table 3). Different variation was found in different soil depth. In 0–20 cm, each mulching treatment was higher than that in CK, and the greatest increase showed in 10–20 cm by 2.81%. Together with the soil layer's deepening, mulching treatments (except R) gradually lower than CK, and average decreased by 0.64% in 30–40 cm. Overall, nearly no difference observed among different four mulching treatments in different soil depth.
Table 3
Percentage of aggregate destruction (PAD) and Unstable aggregates index (ELT) in 0–40 cm soil depth under different treatments after 8-years continuous film mulching.
Treatment | Percentage of aggregate destruction (PAD) | | Unstable aggregates index (ELT) |
0–10 cm | 10–20 cm | 20–30 cm | 30–40 cm | 0–40 cm AVG | | 0–10 cm | 10–20 cm | 20–30 cm | 30–40 cm | 0–40 cm AVG |
PS | 71.78 | 78.50 | 81.67 | 84.93 | 78.92 | | 85.28 | 87.77 | 88.42 | 89.43 | 87.23 |
S | 63.34a | 72.01ab | 77.71a | 81.58a | 74.33a | | 77.91b | 81.41b | 83.15b | 85.71b | 82.05b |
D | 63.31a | 71.77ab | 77.33a | 81.57a | 74.29a | | 79.22ab | 82.02ab | 83.22b | 85.90b | 82.59ab |
R | 63.56a | 72.30a | 76.57a | 82.11a | 74.43a | | 80.13ab | 83.12ab | 84.19b | 86.55ab | 83.50ab |
P | 64.10a | 71.81ab | 76.25a | 80.56a | 73.81a | | 81.12ab | 83.68ab | 84.64ab | 86.53ab | 83.99ab |
CK | 63.08a | 70.01b | 77.60a | 81.97a | 74.07a | | 82.35a | 84.25a | 86.31a | 88.06a | 85.24a |
Note: PS, pre-sowing; S, planting furrows separated by consecutive plastic film-mulched ridges; R, traditional ridge-furrow rainwater harvesting planting; D, double ridge-furrow with full plastic film mulching; P, alternate mulching and bare rows without ridges; CK, flat planting without mulching. Different lower-case letters indicate significant differences at P < 0.05. |
Similar with PAD, the unstable aggregates index (ELT) of each treatment in 0–40 cm soil layer was average decreased by 5.09%, compared to pre-treatment (PS) (Table 3). Four mulching treatments was ordered as S < D < R < P in 0–40 cm depth, which decreased the ELT by 2.74% (P < 0.05), 2.51% (P < 0.05), 1.74%, and 1.77% with S, D, R, and P, respectively, compared to CK. And no significant difference (P > 0.05) observed among different four mulching treatments in different soil depth.
Table 3
Root system characteristic parameter
Maize roots system gradually increased with the growth process, and reached the maximum in 120 days after planting (DAP) of each treatments. During 2019–2020, film mulching treatments significantly (P < 0.05) increased root length (RL), root diameter (RD), root surface area (RS), and root volume (RV) compared with CK, and both of these in S and D were better than that with R and P (Fig. 8). In 2019, the RL of D treatment was 12.02%, 25.82% and 21.06% higher than that under S, R and P in early growth stage (0–30 DAP), respectively, while there were no significant difference in RD and RV between mulching treatments. After sowing 60 days, the RL of each mulching treatments increased to 1996.07-2284.03 cm with the growth process moving forward, and no significant difference between mulching treatments; the RD of S treatment was 11.11% (P < 0.05) and 18.86% (P < 0.05) higher than R and P treatment, respectively. And the RS of S and D were significantly higher than that under R and P by 32.20% and 20.72%; while the RV of S treatment significantly (P < 0.05) increased by 16.40% than R treatment. After sowing 60–120 days, the root growth under S treatment gradually better than others, i.e. the RL was significantly (P < 0.05) increased by 6.30%, 13.30% and 15.34% compared with D, R and P treatment, respectively, and RD, RS and RV average increased by 15.84%, 32.25%, and 24.50%, respectively. The effects of mulching treatment on root characteristics in 2020 was similar to 2019.
Figure 8.
Root biomass (RB), grain and biomass yield (GY and BY)
The film mulching significantly affected maize grain yield (GY) and biomass yield (BY) during experimental years (Fig. 9). Compared with CK (GY = 7.48 t ha− 1 and BY = 17.92 t ha− 1), the GY of each film mulching treatments (S, D, R, and P) were significantly (P < 0.05) increased by 63.96%, 61.81%, 30.99%, and 29.98%, and BY were increased by 44.31%, 44.22%, 22.14%, and 19.26%, respectively. The GY order of mulching treatments was S > D > R > P. The variation coefficient of GY with S and D were lower than R and P, and the S were significant (P < 0.05) increased by 25.17% and 26.14%, and D were 23.53% and 24.49%, respectively. Similar with GY, the BY of S and D were significantly higher than R and P, average increased by 18.12% (P < 0.05) and 20.97% (P < 0.05).
Figure 9.
Film mulching enhanced root dry mass during 2013–2020, and the interaction between the mulching and experimental year was significant (Fig. 10). The mean values with film mulching treatments were all higher than CK, the difference increased gradually with the increase of mulching years, which averaged increased by 47.25% (P < 0.05), 43.59% (P < 0.05), 25.15%, and 22.87% under S, D, R and P, respectively. In each experimental year (except 2013), the order of mulching treatments was S > D > R > P. S and D were significantly higher than R and P in most experimental years (except 2013 and 2018), average increased by 18.13% and 23.39%, and the variation coefficient was also significantly lower (28.61% in average).
Figure 10.
Correlativity
The structural equation modeling (SEM) indicated that all predictor variables explained 94% of variations in crop yield. Film mulching method showed significantly effects on root biomass and crop yield (0.76 and − 0.26, respectively) (Fig. 11A). Root biomass showed significantly positive effects on macroaggregate content, which showed significantly positive effects on soil aggregate stability and crop yield, and negative effects on soil bulk density. And soil bulk density suppressed crop yield (Fig. 11A). In general, film mulching method had positive total effects on root biomass, soil aggregate stability, macroaggregate content and crop yield, and negative total effects on soil bulk density (Fig. 11B).
Figure 11.