Soil physicochemical properties analysis
As shown in Table 1, the soil total organic carbon (TOC), available nitrogen (AN), available phosphorus (AP), available potassium (AK), and the total P (TP) contents decreased with increasing soil depth under three treatments. The pH value of the 0-20 cm layer was considerably lower in the WSC treatment than that in the WS treatment. Meanwhile, the AN and the total N (TN) contents of the 10-20 cm layer were significantly higher in the WSC treatment than those in the WS treatment.
Table 1
Soil physicochemical properties under different treatments. W, winter flooded fallow + no straw returning; WS, winter flooded fallow + straw returning; WSC, winter flooded fallow + straw returning + crayfish farming. Means with different letters for the same property in the same soil layer indicate significant differences at P < 0.05. Values are means ± standard errors.
Soil depths
|
|
pH
|
Total organic carbon
|
Available N
|
Available P
|
Available K
|
Total N
|
Total P
|
Total K
|
(cm)
|
|
|
(g·kg-1)
|
(mg·kg-1)
|
(mg·kg-1)
|
(mg·kg-1)
|
(g·kg-1)
|
(g·kg-1)
|
(g·kg-1)
|
0-10
|
W
|
6.81±0.06a
|
18.20±0.12a
|
114.88±4.42b
|
11.16±0.83a
|
197.05±4.88a
|
2.50±0.01a
|
0.67±0.02a
|
23.01±0.16a
|
WS
|
6.84±0.04a
|
18.49±0.51a
|
127.00±3.45ab
|
12.79±0.16a
|
205.52±5.08a
|
2.56±0.01a
|
0.68±0.00a
|
23.12±0.21a
|
WSC
|
6.42±0.03b
|
19.56±0.87a
|
136.04±4.24a
|
12.01±0.75a
|
194.09±12.25a
|
2.55±0.07a
|
0.67±0.02a
|
23.00±0.50a
|
10-20
|
W
|
6.98±0.01a
|
16.66±0.56a
|
105.12±4.95ab
|
10.61±0.24a
|
171.23±5.08a
|
2.18±0.03b
|
0.61±0.01a
|
23.18±0.30a
|
WS
|
6.91±0.03a
|
16.32±0.37a
|
89.43±2.07b
|
10.29±0.64a
|
178.00±2.96a
|
2.21±0.04b
|
0.62±0.01a
|
23.14±0.32a
|
WSC
|
6.73±0.05b
|
16.42±0.62a
|
108.93±5.44a
|
11.05±0.85a
|
185.62±1.53a
|
2.62±0.15a
|
0.64±0.01a
|
23.03±0.17a
|
Soil reducing substances properties analysis
Soil total reducing substances, Fe2+, and Mn2+contents increased with increasing soil depth (Table 2). Soil Fe2+ content in the 0-20 cm layer, Soil Mn2+ content in the 0-10 cm layer, and total reducing substances in the 0-20 cm layer was significantly lower in the WSC treatment than those in the WS treatment (P < 0.05). Meantime, the soil total reducing substances of the 0-10 cm layer was significantly higher in the WS treatment than that in the W treatment by 28.5% (P < 0.05).
Table 2
Soil reducing substances properties under different treatments. W, winter flooded fallow + no straw returning; WS, winter flooded fallow + straw returning; WSC, winter flooded fallow + straw returning + crayfish farming. Means with different letters for the same property in the same soil layer indicate significant differences at P < 0.05. Values are means ± standard errors.
Soil depths(cm)
|
|
Fe2+(cmol·kg-1)
|
Mn2+(cmol·kg-1)
|
Total Reducing Substances
(cmol·kg-1)
|
0-10
|
W
|
1.220±0.140ab
|
0.084±0.007ab
|
1.921±0.036b
|
WS
|
1.498±0.108a
|
0.095±0.002a
|
2.469±0.148a
|
WSC
|
0.778±0.033b
|
0.074±0.002b
|
1.881±0.121b
|
10-20
|
W
|
1.670±0.063a
|
0.085±0.001a
|
4.683±0.402a
|
WS
|
1.655±0.003a
|
0.092±0.005a
|
3.934±0.084a
|
WSC
|
1.055±0.086b
|
0.094±0.003a
|
2.519±0.155b
|
Soil enzyme activity analysis
Soil enzyme activity decreased with increasing soil depth (Table 3). Soil sucrase and acid phosphatase activities showed an increasing trend in the WSC and WS treatments compared with the W treatment in all of the layers examined. Soil sucrase activity was significantly higher in the WS treatment than that in the W treatment (P < 0.05). Meantime, soil sucrase and acid phosphatase activities in the 10-20 cm layer were higher in the WSC treatment than those in the WS treatment by 44.4% and 32.2%, respectively. However, there was no significant difference in soil urease activity among different treatments in 0-20 cm layer.
Table 3
Soil enzyme activity analysis under different treatments. W, winter flooded fallow + no straw returning; WS, winter flooded fallow + straw returning; WSC, winter flooded fallow + straw returning + crayfish farming. Means with different letters for the same property in the same soil layer indicate significant differences at P < 0.05. Values are means ± standard errors.
Soil depths(cm)
|
|
Acid phosphatase
(mg·phenol·g−1soil·24 h−1)
|
Urease
(mg·NH3-N·g−1 soil·24 h−1)
|
Sucrase
(mg·glucose·g−1soil·24 h−1)
|
0-10
|
W
|
47.38±3.19b
|
2.11±0.08a
|
131.60±6.16b
|
WS
|
54.73±2.72ab
|
2.12±0.03a
|
211.41±4.44a
|
WSC
|
61.71±2.84a
|
2.01±0.06a
|
222.16±6.79a
|
10-20
|
W
|
37.49±1.40b
|
2.11±0.08a
|
122.39±5.93c
|
WS
|
40.43±2.80b
|
2.12±0.06a
|
159.03±6.64b
|
WSC
|
53.44±1.98a
|
2.04±0.08a
|
229.71±3.06a
|
Soil functional diversity of microbial community analysis
The soil average well color development (AWCD) value is one of the indices to judge the total carbon utilization ability of microbial community, reflecting the soil microbial activity and the diversity of physiological functions of microbial community30. With the increasing duration of culture (Figure 1), the utilization degree of different carbon source increased gradually. In the 0-20 cm layer, the slope of the AWCD value was the highest at 24-120 h of culture, indicating that carbon source metabolic activity of soil microbes was the highest at this stage, and then a stable stage period was entered. During hours 72-192 in Fig. 1, the soil AWCD value in the WS treatment was higher than that in the W treatment, and the AWCD value of the WSC treatment was higher than that in the WS treatment in the 0-20 cm layer. This showed that the WSC and WS treatments increased microbial activity in the 0-20 cm layer.
The Shannon index reflects the diversity of bacterial community, with a higher index indicating the stronger the diversity of bacterial community31-32. The Simpson index reflects the changes in the population of each species, with a higher index indicating the position of dominant species more prominent33. The McIntosh index is an index to measure the community species consistency, with an increased degree of carbon source utilization, and consequently an increased McIntosh index34. Table 4 showed that at 96 h of culture in the Biolog ECO microplate, those indices were significant difference in the AWCD value in the 0-20 cm layer among the W, WS, and WSC treatments. In the 0-10 cm layer, the Shannon, Simpson, and McIntosh indices in the WS treatment were significantly higher than those in the W treatment. However, there were no significant difference in Shannon, Simpson, and McIntosh indices of the soil microbial community between the WSC and WS treatments. In the 10-20 cm layer, the Simpson index in the WS treatment had a significant increased compared with the W treatment. The Simpson, and McIntosh indices in the WSC treatment were significantly higher than those in the WS treatment. Therefore, the WS treatment increased significantly the diversity and species abundance of bacterial community in the 0-10 cm layer, the ability of utilizing carbon sources in the 0-10 cm layer, and the species abundance of bacterial community in the 10-20 cm layer. The species abundance of bacterial community, and the ability of utilizing carbon sources were increased significantly in the WSC treatment in the 10-20 cm layer.
Table 4
Functional diversity index of the soil microbial community under different treatments. W, winter flooded fallow + no straw returning; WS, winter flooded fallow + straw returning; WSC, winter flooded fallow + straw returning + crayfish farming. Means with different letters for the same property in the same soil layer indicate significant differences at P < 0.05. Values are means ± standard errors.
Soil depths (cm)
|
|
AWCD of 96 h
|
Shannon Index
|
Simpson Index
|
Mclntosh Index
|
0-10
|
W
|
0.87±0.06b
|
2.89±0.03b
|
0.939±0.002b
|
6.68±0.10b
|
WS
|
1.27±0.02a
|
3.11±0.05a
|
0.951±0.003a
|
8.71±0.08a
|
WSC
|
1.42±0.03a
|
3.21±0.02a
|
0.957±0.001a
|
9.08±0.04a
|
10-20
|
W
|
0.51±0.09b
|
2.50±0.20b
|
0.897±0.006c
|
4.81±0.48b
|
WS
|
0.69±0.10b
|
2.73±0.17ab
|
0.922±0.004b
|
5.40±0.65b
|
WSC
|
1.23±0.13a
|
3.16±0.04a
|
0.954±0.002a
|
8.10±0.69a
|
Analysis of soil microbe carbon source utilization
Based on the types of carbon sources in the Biolog ECO microplate, the 31 kinds of carbon sources were divided into carbohydrates, carboxylic acids, amino acids, phenols, amines and polymers. The changes in the AWCD values of microbial carbon source utilization at 96 h of three treatments were analyzed.
As shown in Fig. 2, with increasing soil layer depth, the soil microbial utilization rate of carbohydrates, carboxylic acids, amino acids, phenols, amines, and polymers showed a gradual decreasing trend among three treatments. Compared with the W treatment, there were a huge increase in the WS treatment among the utilization rates of carbohydrates, carboxylic acids, amino acids, amines and polymers in the 0-20 cm layers by 47.5-52.4%, 31.4-63.0%, 27.9-45.7%, 1.4-51.8%, and 51.8-65.0%, respectively. The WSC treatment also represented the same tendency in the 0-20 cm layers compared with the WS treatment, in which the utilization rates of phenols in the 0-10 cm layer and the utilization rates of carbohydrates, amino acids, and amines in the 10-20 cm layer had significant differences. Consequently, the straw returning significantly increased amino acids and polymers utilization rates in the 0-10 cm layer, carboxylic acids utilization rates in the 10-20 cm layer. The straw returning + crayfish farming significantly increased the utilization rates of phenols in the 0-10 cm layer and the utilization rates of carbohydrates, amino acids, and amines in the 10-20 cm layer.
Interaction of physicochemical and reducing substances properties with soil enzyme activity and microbial diversity indices
As shown in table 5, the pH value, Fe2+ content, and total reducing substances had a highly significant negative correlation with the AWCD of 96 h, acid phosphatase, sucrase, Shannon index, Simpson index, and McIntosh index. The TOC, AK, AN, and TP contents had significant positive relationship with the AWCD of 96 h, acid phosphatase, Shannon index, Simpson index, and McIntosh index. The TN content had also significant positive relationship with the AWCD of 96 h, acid phosphatase, sucrase, Shannon index, Simpson index, and McIntosh index. In addition, the AP content with the AWCD of 96 h and McIntosh index showed a significant direct correlation. Thus, it can be seen that the pH value and TN content had a closer relationship with soil enzymes; the pH value, AK, TN, and Fe2+ contents and total reducing substances had a closer relationship with the soil microbial diversity indices.
Table 5
Correlation coefficient of soil physicochemical properties, soil reducing substances properties, soil enzyme activity and microbial diversity index. * significant at P < 0.05. ** significant at P < 0.01.
|
AWCD of 96 h
|
Acid phosphatase
|
Urease
|
Sucrase
|
McIntosh index
|
Simpson index
|
Shannon index
|
pH
|
-0.719**
|
-0.812**
|
0.317
|
-0.644**
|
-0.697**
|
-0.702**
|
-0.654**
|
TOC
|
0.503*
|
0.507*
|
-0.31
|
0.256
|
0.498*
|
0.483*
|
0.514*
|
AP
|
0.537*
|
0.466
|
0.013
|
0.429
|
0.565*
|
0.442
|
0.421
|
AK
|
0.712**
|
0.654**
|
0.006
|
0.435
|
0.752**
|
0.778**
|
0.685**
|
AN
|
0.721**
|
0.641**
|
-0.321
|
0.454
|
0.727**
|
0.588*
|
0.583*
|
TK
|
-0.259
|
-0.07
|
0.014
|
-0.081
|
-0.258
|
-0.221
|
-0.346
|
TN
|
0.748**
|
0.644**
|
-0.416
|
0.637**
|
0.762**
|
0.789**
|
0.668**
|
TP
|
0.616**
|
0.572*
|
-0.151
|
0.404
|
0.629**
|
0.651**
|
0.501*
|
Fe2+
|
-0.766**
|
-0.723**
|
0.373
|
-0.556*
|
-0.767**
|
-0.781**
|
-0.688**
|
Mn2+
|
-0.132
|
-0.212
|
0.093
|
-0.011
|
-0.148
|
-0.028
|
-0.034
|
Total Reducing Substances
|
-0.772**
|
-0.780**
|
0.147
|
-0.540*
|
-0.770**
|
-0.890**
|
-0.773**
|