谷子不同輪作模式對(duì)土壤理化性質(zhì)及細(xì)菌群落的影響

李 丹， 李小霞， 李萬星， 曹晉軍， 靳鯤鵬， 韓文清， 胡丹珠， 劉 鑫， 田 崗， 黃學(xué)芳， 劉永忠

(1.山西農(nóng)業(yè)大學(xué)谷子研究所，山西長(zhǎng)治046000；2.山西農(nóng)業(yè)大學(xué)山西有機(jī)旱作農(nóng)業(yè)研究院，山西太原030000)

連作是指同一種植物或近緣植物在同一塊地上連續(xù)種植兩茬或以上的現(xiàn)象。連作可導(dǎo)致土壤中的某種養(yǎng)分過度消耗或產(chǎn)生自毒作用使得作物生長(zhǎng)發(fā)育遲緩，病蟲害發(fā)生概率增加，農(nóng)產(chǎn)品產(chǎn)量降低。谷子是不宜重茬的作物之一，近年來谷子種植面積不斷縮減，但是由于其豐富的營(yíng)養(yǎng)價(jià)值和保健作用，較高的經(jīng)濟(jì)效益，農(nóng)戶種植經(jīng)濟(jì)性較高，導(dǎo)致連作現(xiàn)象頻發(fā)[1]。

已有研究結(jié)果表明，避免連作障礙的最佳防治途徑是與不同作物輪作[2]。輪作可以合理利用土壤養(yǎng)分，調(diào)節(jié)土壤酸堿度，改善土壤結(jié)構(gòu)，減輕病蟲害發(fā)病率。孫倩等[3]發(fā)現(xiàn)，輪作能夠改善土壤酶活性，輪作模式下谷茬土壤中的過氧化氫酶、脲酶、堿性磷酸酶活性提高，土壤細(xì)菌群落豐富度指數(shù)降低。蔡晨等[4]研究發(fā)現(xiàn)，與水稻單作相比，稻蝦長(zhǎng)期輪作可使0～20.0 cm和20.1～40.0 cm土層容重降低，非毛管孔隙度、毛管孔隙度、總空隙度提高，同時(shí)增強(qiáng)土壤緩沖能力，提高耕層土壤pH。熊湖等[5]研究結(jié)果表明，連作可以改變馬鈴薯植株形態(tài)特征(株高、莖粗、總?cè)~面積降低，根長(zhǎng)增加)，通過輪作使土壤中速效養(yǎng)分轉(zhuǎn)化速率加快，可以緩解連作脅迫。

本研究擬以特色雜糧谷子為研究對(duì)象，探討不同作物與谷子輪作后對(duì)土壤的影響，包括物理結(jié)構(gòu)、土壤養(yǎng)分、細(xì)菌多樣性，篩選適用于晉東南武鄉(xiāng)地區(qū)的最佳谷子輪作模式，以期對(duì)谷子特色產(chǎn)業(yè)的可持續(xù)、綠色發(fā)展提供指導(dǎo)。

1 材料與方法

1.1 試驗(yàn)地概況及試驗(yàn)設(shè)計(jì)

大田試驗(yàn)于2017-2020年在山西省長(zhǎng)治市武鄉(xiāng)縣上司鎮(zhèn)鋪上村(E112°93′，N36°76′)進(jìn)行。2018年至2020年谷子生育期內(nèi)平均降水量為379.7 mm，為典型的雨養(yǎng)農(nóng)業(yè)區(qū)。土壤為典型的紅壤土，前茬為谷子，基礎(chǔ)肥力為：有機(jī)質(zhì)含量10.24 g/kg，全氮含量0.089%，速效磷含量29.14 mg/kg，速效鉀含量246.35 mg/kg，堿解氮含量46.96 mg/kg，pH值為8.28。試驗(yàn)設(shè)置3次重復(fù)，隨機(jī)區(qū)組排列，小區(qū)面積7.7 m×7.0 m，有機(jī)肥為羊糞，施入量為60 000 kg/hm2，不同小區(qū)管理方式一致。各處理設(shè)置見表1。

表1 試驗(yàn)設(shè)置

供試谷子為長(zhǎng)生13，玉米為晉單73，大豆為長(zhǎng)豆35號(hào)，高粱為晉雜22號(hào)，花生為長(zhǎng)治縣裕盛豐農(nóng)業(yè)專業(yè)合作社提供的黑美人，馬鈴薯為晉薯16號(hào)，甘薯為山西省長(zhǎng)治市長(zhǎng)子縣當(dāng)?shù)靥赜衅贩N，玉米(糯)為晉糯10號(hào)。

1.2 土壤樣品采集

用于測(cè)定土壤性質(zhì)及土壤團(tuán)聚體的土壤樣品于谷子成熟期采集：清除地表雜物后用環(huán)刀分別采集0～10.0 cm、10.1～20.0 cm的土壤裝入鋁盒內(nèi)帶回實(shí)驗(yàn)室并置于室內(nèi)風(fēng)干，運(yùn)輸過程中盡量避免碰撞擠壓。用于測(cè)定土壤肥力的根際土壤樣品采集：于收獲期采集每種作物的根際土壤，隨機(jī)選取5株作物，將作物連根拔起，抖落外側(cè)大型土塊，去除石子、根系等雜物，充分混勻后過篩，隨后帶回實(shí)驗(yàn)室自然風(fēng)干，備用。用于測(cè)定細(xì)菌多樣性的根際土壤樣品采集：于灌漿期分別采集不同種植模式的根際土壤，方法同上。每個(gè)處理樣品分成2份，一份放入凍存管中，液氮保存寄送至上海派森諾生物科技有限公司用于土壤真菌多樣性的檢測(cè)，一份放入-80 ℃?zhèn)溆谩?/p>

1.3 測(cè)定指標(biāo)與方法

1.3.1 土壤容重及土壤水穩(wěn)性團(tuán)聚體測(cè)定 于谷子收獲期采用環(huán)刀法在田間測(cè)定土壤容重。

土壤水穩(wěn)性團(tuán)聚體采用濕篩法[6-8]測(cè)定：將處理好的土壤置于團(tuán)聚體分析儀套篩上部篩網(wǎng)上，套篩孔徑依次為2.00 mm、1.00 mm、0.50 mm和0.25 mm，桶中加入蒸餾水使水沒過土樣，浸泡10 min后開始振蕩(振蕩10 min，上下移動(dòng)30次，振幅3.00 cm)。振蕩結(jié)束后將套篩取出，將每個(gè)篩子上面的土壤沖洗到鋁盒內(nèi)，靜置15 min后用塑料滴管去除上清液，烘箱烘干，稱量，得到粒徑<0.250 mm、0.251～0.500 mm、0.501～1.000 mm、1.001～2.000 mm、>2.000 mm的水穩(wěn)性團(tuán)聚體質(zhì)量。

1.3.2 土壤養(yǎng)分測(cè)定 土壤養(yǎng)分、酶活性的測(cè)定采用常規(guī)測(cè)定方法[9-11]。

1.3.3 土壤細(xì)菌多樣性測(cè)定 根據(jù)DNeasy PowerSoil Kit(QIAGEN)試劑盒說明書對(duì)土壤總DNA進(jìn)行提取，利用NanoDrop2000對(duì)提取的DNA質(zhì)量濃度和純度進(jìn)行檢測(cè)；對(duì)細(xì)菌16S區(qū)進(jìn)行擴(kuò)增。對(duì)構(gòu)建好的文庫(kù)進(jìn)行檢驗(yàn)、測(cè)序。對(duì)測(cè)序結(jié)果進(jìn)行操作單位聚類與分析、Alpha多樣性分析、Beta多樣性分析及不同分類水平的群落結(jié)構(gòu)分析。

1.4 數(shù)據(jù)計(jì)算與統(tǒng)計(jì)分析

采用Excel2000軟件進(jìn)行數(shù)據(jù)整理與分析，用IBM SPSS20.0軟件進(jìn)行單因素方差分析和主成分分析。

計(jì)算各粒級(jí)團(tuán)聚體含量、各粒級(jí)團(tuán)聚體比例、平均質(zhì)量直徑(MWD)，公式[12-14]如下：

各粒級(jí)團(tuán)聚體含量=各級(jí)團(tuán)聚體質(zhì)量/土樣樣品干質(zhì)量×100%

(1)

各粒級(jí)團(tuán)聚體比例=某級(jí)團(tuán)聚體質(zhì)量/團(tuán)聚體總質(zhì)量×100%

(2)

(3)

式中，MWD為平均質(zhì)量直徑；wi為第i級(jí)團(tuán)聚體的質(zhì)量占土壤樣品干質(zhì)量的質(zhì)量百分?jǐn)?shù)；Ri為第i級(jí)團(tuán)聚體平均直徑。本研究中各級(jí)團(tuán)聚體平均直徑分別為2.00 mm、1.50 mm、0.75 mm、0.38 mm、0.25 mm。

2 結(jié)果與分析

2.1 谷子不同輪作模式對(duì)土壤容重和水穩(wěn)性團(tuán)聚體的影響

2.1.1 土壤容重 圖1顯示，不同處理下土壤容重由大到小為：SSSi、SZGm>SZZm>SZSt、SZLb>SZSb>SZAh，SSSi處理下土壤容重為1.25 g/cm3，其他6個(gè)處理的土壤平均容重為1.19 g/cm3，較SSSi處理下的土壤容重降低了4.8%。其中SZAh、SZSb、SZZm、SZSt、SZLb處理下土壤容重較SSSi處理下的土壤容重降低了2.4%～8.8%，SZAh處理下土壤容重最低。

SSSi、SZGm、SZZm、SZSt、SZLb、SZSb、SZAh見表1。不同小寫字母表示不同處理間差異顯著(P<0.05)。圖1 不同處理下的土壤容重Fig.1 Soil bulk density under different treatments

2.1.2 土壤水穩(wěn)性團(tuán)聚體的粒徑分布 圖2顯示，0～10.0 cm和10.1～20.0 cm土層中各粒級(jí)團(tuán)聚體占比變化趨勢(shì)一致：占比隨團(tuán)聚體粒徑減小而降低。在0～10.0 cm土層，各處理下粒徑>2.000 mm土壤水穩(wěn)性團(tuán)聚體占比為42.98%～53.26%，在SZSb處理下的占比最低，在SZSt處理下的占比最高；粒徑為1.001～2.000 mm的土壤水穩(wěn)性團(tuán)聚體占比26.72%～32.72%，在SZSt處理下的占比最低，在SZSb處理下的占比最高；粒徑為0.501～1.000 mm的土壤水穩(wěn)性團(tuán)聚體占比9.90%～16.09%，在SZAh處理下的占比最低，在SZSb處理下的占比最高；粒徑為0.251～0.500 mm的土壤水穩(wěn)性團(tuán)聚體占比6.16%～7.86%，在SSSi處理下的占比最低，在SZZm處理下的占比最高；粒徑為0.100～0.250 mm的土壤水穩(wěn)性團(tuán)聚體占比0.22%～4.67%，在SSSi處理下的占比最低，在SZZm處理下的占比最高。在10.1～20.0 cm土層，各處理中粒徑>2.000 mm的土壤水穩(wěn)性團(tuán)聚體占比為41.78%～52.55%，在SZSb處理下的占比最低，在SZLb處理下的占比最高；SZSb處理中粒徑為1.001～2.000 mm、0.501～1.000 mm、0.251～0.500 mm的土壤水穩(wěn)性團(tuán)聚體占比均最高，分別為32.94%、15.17%、8.13%；粒徑為1.001～2.000 mm的土壤水穩(wěn)性團(tuán)聚體在SZSt處理中占比最低，為26.61%；粒徑為0.501～1.000 mm的土壤水穩(wěn)性團(tuán)聚體在SZLb處理中占比最低，為11.45%；粒徑為0.251～0.500 mm、0.100～0.250 mm的土壤水穩(wěn)性團(tuán)聚體在SSSi處理中占比最低，為6.70%、0.17%；粒徑為0.100～0.250 mm的土壤水穩(wěn)性團(tuán)聚體在SZGm處理中占比最高，為5.51%。

SSSi、SZGm、SZZm、SZSt、SZLb、SZSb、SZAh見表1。圖2 不同處理土壤團(tuán)聚體粒級(jí)分布Fig.2 Particle size distribution of soil aggregates under different treatments

2.1.3 土壤水穩(wěn)性團(tuán)聚體組成及穩(wěn)定性 可用平均質(zhì)量直徑(MWD)反映土壤水穩(wěn)性團(tuán)聚體穩(wěn)定性，MWD越大，代表土壤團(tuán)聚體的團(tuán)聚度越高，該土壤越穩(wěn)定[15-16]。表2顯示，同一處理10.1～20.0 cm土層的MWD均大于0～10.0 cm土層。0～10.0 cm土層中，MWD由大到小為：SZGm>SZZm>SSSi、SZSb>SZAh>SZLb、SZSt；10.1～20.0 cm土層中，MWD由大到小為：SZGm>SZZm>SZSb>SSSi>SZAh、SZSt>SZLb。在0～10.0 cm、10.1～20.0 cm分層土壤中，SZGm處理下的MWD值最大，且與其他處理差異顯著，說明SZGm處理下土壤穩(wěn)定性最優(yōu)，SZSt、SZLb處理下土壤穩(wěn)定性差，猜測(cè)塊莖、塊根類作物可降低土壤團(tuán)聚體穩(wěn)定性。

2.2 谷子不同輪作模式對(duì)土壤基礎(chǔ)肥力的影響

2.2.1 土壤養(yǎng)分 表3顯示，不同處理下土壤全氮含量、堿解氮含量、速效磷含量、速效鉀含量、有機(jī)質(zhì)含量、pH與初始土樣相比均發(fā)生改變。土壤全氮、有機(jī)質(zhì)含量均提高，增幅達(dá)30.3%～59.6%、96.4%～139.0%，SSSi處理下全氮、有機(jī)質(zhì)含量最高，SZGm處理下全氮含量最低，SZAh、SZLb處理下有機(jī)質(zhì)含量最低；不同處理下土壤均為堿性狀態(tài)，SSSi、SZZm處理下土壤pH與初始土壤相比呈降低趨勢(shì)；SZGm、SZLb、SZZm、SSSi處理下堿解氮含量升高，SZAh、SZSt處理下降低，SZZm處理下堿解氮含量最高，SZAh處理下堿解氮含量最低；SZGm、SSSi、SZSb、SZAh處理下速效磷含量升高，增幅為0.9%～56.5%，SZLb、SZSt、SZZm處理下速效磷含量降低，降幅為13.9%～26.2%；SSSi處理下的速效鉀含量變化不明顯，其他處理下的速效鉀含量均明顯降低，降幅為16.9%～26.8%，SZLb處理下的速效鉀含量最低。

表2 不同處理土壤水穩(wěn)性團(tuán)聚體組成特征

表3 不同處理土壤養(yǎng)分含量

2.2.2 土壤酶活性 不同種植模式會(huì)影響土壤酶活性。過氧化氫酶可催化土壤中的過氧化氫分解為水和氧氣，是一種與土壤呼吸作用相關(guān)的酶；同時(shí)可以分解土壤中的有害自由基，與自凈能力有關(guān)。表4顯示，不同處理下過氧化氫酶活性從大到小為：SZZm>SZGm>SZAh>SZSb>SSSi>SZLb>SZSt，說明SZZm處理下土壤呼吸強(qiáng)度高，自凈能力強(qiáng)，SZSt處理下最差。堿性磷酸酶在磷轉(zhuǎn)化方面發(fā)揮著重要作用，可以加速有機(jī)磷的脫磷速度，該酶主要由土壤細(xì)菌分泌，可將土壤中的有機(jī)磷水解為可供植物利用的無機(jī)磷[17]，不同處理下堿性磷酸酶活性由大到小為：SZSb>SZGm>SZLb、SZZm>SZSt>SSSi、SZAh。

脲酶可以分解土壤中的含氮化合物，是速效氮的主要來源，脲酶活性是評(píng)價(jià)土壤酶活性和肥力的重要指標(biāo)，不同處理下土壤中脲酶活性大小順序依次為：SZZm>SSSi>SZAh>SZLb>SZGm>SZSb>SZSt，SZZm處理下脲酶活性最高，比脲酶活性最低的SZSt處理高31.0%，其他5個(gè)處理間差異不顯著。

纖維素酶可由土壤真菌、細(xì)菌和放線菌分泌產(chǎn)生，能夠降解秸稈中的纖維素，將養(yǎng)分釋放到土壤中[18]。不同處理下土壤中纖維素酶活性由大到小為：SZSt>SZLb>SZAh>SSSi>SZZm>SZGm>SZSb，SZSt處理下纖維素酶活性最大，比酶活性最低的SZSb處理高47.8%。

轉(zhuǎn)化酶又稱蔗糖酶，能夠促進(jìn)蔗糖生成葡萄糖和果糖，可以促進(jìn)土壤熟化，增加土壤中的速效養(yǎng)分[19]。不同處理下土壤中蔗糖酶活性由大到小為：SZLb>SZZm、SZAh>SZSb>SSSi>SZSt、SZGm，SZLb處理下蔗糖酶活性最高，比活性最低的SZGm、SZSt處理高35.3%。

蛋白酶可將土壤中蛋白質(zhì)和含氮化合物分解為植物可以直接吸收利用的氨基酸或活性氮[20]。不同處理下土壤中蛋白酶活性由大到小為：SZLb>SZSt>SZZm>SSSi>SZGm>SZSb>SZAh，SZLb處理下蛋白酶活性最高，且與其他處理差異顯著，比酶活性最低的SZAh處理高174.2%。

表4 不同處理土壤酶活性

2.3 谷子不同輪作模式對(duì)細(xì)菌群落結(jié)構(gòu)的影響

2.3.1 土壤細(xì)菌組成 圖3顯示，不同處理的土壤具有相似的細(xì)菌組成，但是比例上存在一定差異。其中相對(duì)豐度大于1.00%的有6個(gè)門類，分別是Proteobacteria、Actinobacteria、Chloroflexi、Acidobacteria、Bacteroidetes和Gemmatimonadetes，6個(gè)門的豐度占總豐度的94.50%～96.70%。其中，Proteobacteria、Actinobacteria為優(yōu)勢(shì)菌門，在各處理中占比達(dá)64.20%～74.10%。與SSSi處理相比，SZLb、SZSt處理中Proteobacteria的相對(duì)豐度加，分別增加了2.30%、1.60%，其他處理中Proteobacteria的相對(duì)豐度降低，降低了3.60%～17.20%。與SSSi處理相比，SZSt、SZSb、SZAh、SZGm處理中Actinobacteria的相對(duì)豐度增加了3.98%～30.70%，SZLb、SZZm處理中Actinobacteria的相對(duì)豐度降低了7.50%、6.30%。另外，Nitrospirae、Verrucomicrobia、Patescibacteria、Rokubacteria是另外4個(gè)主要的門(至少在一個(gè)處理中的相對(duì)豐度大于0.50%)。

SSSi、SZGm、SZZm、SZSt、SZLb、SZSb、SZAh見表1。圖3 不同處理土壤細(xì)菌門水平的相對(duì)豐度Fig.3 Relative abundance of soil bacteria at phyla level in different treatments

不同處理土壤細(xì)菌屬水平群落結(jié)構(gòu)分析結(jié)果(圖4、圖5)表明，7種輪作模式聚類圖沒有明顯分支，說明土壤微生物結(jié)構(gòu)相似，但是比例上存在差異。相對(duì)豐度大于1.00%或至少在1個(gè)處理中的相對(duì)豐度大于1.00%的優(yōu)勢(shì)屬有27個(gè)，7個(gè)處理共有的優(yōu)勢(shì)屬有9種，分別為Subgroup_6占比4.5%～8.1%，BIrii41占比約1.9%～3.9%，KD4-96占比1.8%～3.4%，Sphingomonas占比1.8%～3.4%，MND1占比1.3%～3.6%，Blastococcus占比1.5%～2.4%，RB41占比1.0%～2.2%，Nocardioides占比1.0%～2.1%，Haliangium占比1.0%～1.3%。TRA3-20、IMCC26256為SZGm處理特有的屬，Acidibacter、CCD24為SZZm處理特有的屬，Micromonospora為SZAh處理特有的屬，剩余13屬為2個(gè)或多個(gè)處理的共有屬。聚類分析結(jié)果(圖5)表明，SZLb、SZSb處理與SSSi處理距離最遠(yuǎn)，說明這2個(gè)處理對(duì)于細(xì)菌群落的影響大于其他處理，SZAh處理與SSSi處理距離最近，說明這2個(gè)處理的細(xì)菌組成最相似。

SSSi、SZGm、SZZm、SZSt、SZLb、SZSb、SZAh見表1。圖4 不同處理土壤細(xì)菌屬水平相對(duì)豐度的比較(前20)Fig.4 Comparison of relative abundance of soil bacteria at genus level in different treatments (top 20)

SSSi、SZGm、SZZm、SZSt、SZLb、SZSb、SZAh見表1。圖5 不同輪作模式聚類分析Fig.5 Hierarchical clustering analysis of different crop rotation patterns

2.3.2 土壤細(xì)菌多樣性 各處理的Chao1指數(shù)(豐富度)間沒有顯著差異，與SSSi處理相比，其他6個(gè)處理的Chao1指數(shù)均有所提高；與SSSi處理相比，其他6個(gè)處理的香濃指數(shù)均降低；各處理的辛普森指數(shù)間沒有顯著差異(表5)。綜合3個(gè)指數(shù)可推斷，與連作相比，輪作會(huì)改變土壤微生物群落結(jié)構(gòu)：物種種類(豐富度)增多，而物種多樣性降低，更利于無優(yōu)勢(shì)種群的生存。

表5 不同處理土壤細(xì)菌多樣性指數(shù)

2.4 土壤養(yǎng)分與細(xì)菌優(yōu)勢(shì)屬的相關(guān)性分析

表6顯示，蛋白酶活性是影響土壤細(xì)菌群落最主要的因子。過氧化氫酶活性與Devosia呈顯著負(fù)相關(guān)(P<0.05)；堿性磷酸酶活性與Gitt-GS-136呈顯著正相關(guān)(P<0.05)；纖維素酶活性與Altererythrobacter呈顯著正相關(guān)(P<0.05)；轉(zhuǎn)化酶活性與Haliangium呈顯著正相關(guān)(P<0.05)；蛋白酶活性與Blastococcus呈極顯著負(fù)相關(guān)(P<0.01)，與Lysobacter呈極顯著正相關(guān)(P<0.01)；堿解氮含量與MND1、Nitrospira、TRA3-20、CCD24呈顯著正相關(guān)(P<0.05)，與Actinoplanes呈顯著負(fù)相關(guān)(P<0.05)；速效磷、速效鉀含量與Sphingomonas呈顯著正相關(guān)(P<0.05)。

2.5 谷子不同輪作模式對(duì)土壤物理、肥力指標(biāo)影響的綜合分析

2.5.1 土壤相關(guān)指標(biāo)選取與主成分分析結(jié)果 選取具有代表性的土壤相關(guān)因子，包括：物理結(jié)構(gòu)(土壤容重、10.1～20.0 cm的MWD)、土壤養(yǎng)分(全氮、堿解氮、速效磷、速效鉀、有機(jī)質(zhì))含量及土壤酶(過氧化氫酶、堿性磷酸酶、脲酶、纖維素酶、轉(zhuǎn)化酶、蛋白酶)活性，共13個(gè)指標(biāo)作為主成分分析的評(píng)價(jià)指標(biāo)。為使選取因子在數(shù)量級(jí)和量綱上達(dá)到一致，需對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，然后對(duì)土壤物理、養(yǎng)分指標(biāo)進(jìn)行相關(guān)性分析。

谷子土壤相關(guān)指標(biāo)主成分分析結(jié)果(表7)表明，從13個(gè)因子中提取出5個(gè)主成分且5個(gè)主成分累積貢獻(xiàn)率高達(dá)98.076%，即5個(gè)主成分可包含原有信息的98.076%，符合分析要求，所以這5個(gè)主成分可以代表谷子土壤相關(guān)指標(biāo)進(jìn)行分析。

表6 土壤肥力指標(biāo)與優(yōu)勢(shì)屬之間的相關(guān)性分析

2.5.2 不同處理谷子土壤物理、肥力綜合排名 綜合分析土壤理化性質(zhì)等相關(guān)指標(biāo)，得分最高的為SZZm處理(表8)，其次為SSSi處理，得分最低的為SZSt處理。認(rèn)為谷子-玉米-玉米(糯)-谷子輪作模式可以使土壤環(huán)境保持最優(yōu)。

3 討論

土壤肥力對(duì)農(nóng)作物生長(zhǎng)發(fā)育至關(guān)重要[21-26]，合理輪作可以提高土壤養(yǎng)分含量，改善土壤結(jié)構(gòu)[271]。本研究結(jié)果表明，合理的輪作可降低土壤容重，增強(qiáng)土壤團(tuán)聚體穩(wěn)定性(塊根、塊莖類作物除外)；谷子連作降低了土壤pH，使土壤趨于酸化，這與牛倩云等[28]、徐寧等[29]的研究結(jié)果一致；谷子連作模式下，土壤中速效磷、速效鉀、有機(jī)質(zhì)含量均居于首位，袁洪超等[30]在當(dāng)歸連作研究中得到類似結(jié)論。這可能與不同種植模式作物對(duì)養(yǎng)分的需求不同有關(guān)。

表7 主成分特征值與貢獻(xiàn)率

表8 不同輪作模式因子得分與綜合排名

土壤酶是作物根系和微生物產(chǎn)生的一類蛋白質(zhì)代謝產(chǎn)物，參與土壤的各種生化反應(yīng)。有研究認(rèn)為，合理輪作有利于提高土壤酶活性[31-32]。本研究中，馬鈴薯輪作模式的過氧化氫酶、脲酶活性最低；谷子連作模式、花生輪作模式堿性磷酸酶活性最低；高粱輪作模式纖維素酶活性最低?？傮w來講，合理的輪作有利于提高土壤酶活性，不合理的輪作會(huì)降低土壤酶活性。本試驗(yàn)中7種輪作模式的土壤酶活性沒有明顯規(guī)律，可能與試驗(yàn)地常年施用有機(jī)肥以及試驗(yàn)周期較短(4年)、連作障礙不明顯有關(guān)。

谷子連作模式的土壤細(xì)菌豐富度指數(shù)最低，但是其多樣性指數(shù)最高，推測(cè)與土壤酸化有關(guān)，該條件下更利于優(yōu)勢(shì)菌的生存。聚類分析結(jié)果表明，不同輪作模式細(xì)菌群落結(jié)構(gòu)組成基本相同，但在比例上有所差異，推測(cè)土壤細(xì)菌對(duì)于環(huán)境變化較不敏感，前人在番茄上得到類似結(jié)論[33]。

4 結(jié)論

為了篩選出適宜晉東南地區(qū)谷子最佳的輪作作物，使土壤保持健康生態(tài)環(huán)境,本研究通過主成分分析法可知，不同作物與谷子輪作的綜合得分順序?yàn)镾ZZm>SSSi>SZSb>SZGm>SZLb>SZAh>SZSt。其中，輪作玉米(糯)得分最高，其次依次為谷子連作、輪作高粱、輪作大豆，且四者得分均為正值，說明輪作玉米(糯)、高粱、大豆以及連作谷子可使谷子土壤保持較健康的水平。輪作馬鈴薯、花生、甘薯三者得分為負(fù)值，說明這3種輪作方式降低了土壤質(zhì)量。由此可知，不同輪作模式會(huì)影響土壤的質(zhì)量，且長(zhǎng)時(shí)間連作會(huì)產(chǎn)生連作障礙，玉米(糯)是較為適宜晉東南武鄉(xiāng)地區(qū)避免谷子連作障礙的理想作物。