蘇堯，葉蘇梅，魯夢醒，馬躍，王玉寶，王珊珊，柴如山，葉新新，張震，馬超*

（1. 安徽省綠色磷肥智能制造與高效利用工程研究中心， 農(nóng)田生態(tài)保育與污染防控安徽省重點實驗室， 安徽農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院， 安徽 合肥 230036；2. 肥西縣嚴(yán)店鄉(xiāng)植保服務(wù)專業(yè)合作社， 安徽 合肥 231232；3. 安徽省文勝生物工程股份有限公司， 安徽 阜陽 236000）

雜草發(fā)生狀況是影響農(nóng)田生態(tài)系統(tǒng)穩(wěn)定性的重要因素之一， 草害防治一直是農(nóng)田生態(tài)管理的重要工作[1］。單施除草劑因其對作物引起的藥害頻發(fā)， 以及雜草抗藥性增強而難以滿足農(nóng)田除草需求[2-3］。秸稈作為農(nóng)田生態(tài)系統(tǒng)中的重要資源， 已有大量研究表明其富含的養(yǎng)分可有效改善土壤環(huán)境， 有利于作物的優(yōu)質(zhì)增產(chǎn)， 然而關(guān)于其在不同環(huán)境條件和管理措施下還田對雜草發(fā)生的影響仍不明確[4-6］。因此， 研究不同環(huán)境條件和管理措施下秸稈還田對草害發(fā)生的影響， 有利于秸稈還田條件下的草害科學(xué)防控， 以減輕過度依賴化學(xué)除草帶來的問題。

目前關(guān)于不同環(huán)境條件和管理措施下秸稈還田對農(nóng)田雜草的影響尚無一致結(jié)論。張大偉等[7］在江蘇水稻（Oryza sativa）田上的試驗得出，與免耕不還田相比， 免耕條件下將小麥（Triticum aestivum）秸稈以2250 kg·hm-2的量還田后雜草總株數(shù)顯著降低了97.3%。王慧敏等[8］在內(nèi)蒙古春小麥田上的研究結(jié)果則表明， 深翻條件下秸稈還田比不還田春小麥成熟期雜草密度增加了22.3%。張震等[9］在安徽巢湖流域小麥田上研究發(fā)現(xiàn)秸稈還田可提高雜草多樣性， 而陳浩等[10］在陜西的試驗卻表明常規(guī)施肥條件下秸稈還于稻田與不還田相比， 雜草Shannon指數(shù)顯著降低了61.5%。此外， 農(nóng)田雜草發(fā)生對秸稈還田地響應(yīng)在不同條件下有所差別。前期相關(guān)研究較為關(guān)注不同作物秸稈類型、還田方式等條件下秸稈還田對草害發(fā)生的影響， 而對于氣候條件、土壤性質(zhì)等相關(guān)因子的影響并未給予足夠的重視。桂浩然等[11］在河南開封開展的研究發(fā)現(xiàn)， 雙倍秸稈還田于大豆（Glycine max）地后雜草Shannon 指數(shù)增加了26.0%， 而還于玉米（Zea mays）地后雜草Shannon 指數(shù)顯著降低了71.6%。李貴等[12］在江蘇沿江地區(qū)沙壤土上的研究表明， 秸稈全量覆蓋還田的抑草效果要顯著優(yōu)于全量深埋還田。秸稈還田對農(nóng)田草害影響的關(guān)聯(lián)因素有很多， 如作物類型、耕作方式、秸稈還田量等， 現(xiàn)有研究大多側(cè)重于單一因素對秸稈還田抑草效應(yīng)的影響分析， 而不同影響因子在此過程中的相對重要性仍不明確。

因此， 本研究整理全國范圍內(nèi)的相關(guān)獨立試驗數(shù)據(jù)， 運用整合分析方法， 綜合分析不同環(huán)境條件和管理措施下秸稈還田對我國農(nóng)田雜草多度和多樣性的影響， 評估作物類型、秸稈類型、秸稈還田量、年均氣溫、年均降水量等因子對秸稈還田影響草害發(fā)生過程的關(guān)聯(lián)度， 并比較分析了雜草多度和多樣性對不同條件下秸稈還田的不同響應(yīng)機制及其原因， 以期為秸稈還田條件下優(yōu)化田間管理措施， 以便于有效抑制農(nóng)田草害的發(fā)生提供科學(xué)指導(dǎo)。

1 材料與方法

1.1 數(shù)據(jù)來源

本研究基于中國知網(wǎng)（CNKI）和Web of Science （WoS）等文獻數(shù)據(jù)庫， 以“秸稈還田”“雜草數(shù)量”“雜草密度”“雜 草 發(fā) 生 規(guī) 律”“雜 草 多 樣 性”和“straw returning”“weed abundance”“weed diversity”等 為 關(guān) 鍵 詞 檢 索 于2000-2022 年間發(fā)表的相關(guān)論文。對檢索出的文獻按照以下標(biāo)準(zhǔn)進行篩選： 1）田間試驗范圍限于中國； 2）同一試驗需包含相應(yīng)的對照組（秸稈不還田）和處理組（秸稈還田），且對照組與處理組除秸稈是否還田外， 其余試驗條件（包括環(huán)境條件、農(nóng)藝措施等）應(yīng)完全一致； 3）試驗處理的重復(fù)數(shù)≥3； 4）從結(jié)果中可直接或通過計算間接獲得雜草密度或雜草多樣性數(shù)據(jù)。為了消除試驗指標(biāo)差異所帶來的影響， 本研究選擇雜草密度反映雜草多度， 選擇Shannon 指數(shù)反映雜草多樣性[13］。本研究提取文獻中對照組和處理組的雜草密度、密度標(biāo)準(zhǔn)差及雜草Shannon多樣性指數(shù)、Shannon 多樣性指數(shù)標(biāo)準(zhǔn)差、處理重復(fù)數(shù)等信息。直接提取表格中的相關(guān)數(shù)據(jù)， 使用GetData Graph Digitizer 軟件提取以圖片形式展示的數(shù)據(jù)， 然后利用R studio 軟件進行標(biāo)準(zhǔn)差數(shù)據(jù)的補充[14］， 最終共獲得426 組試驗數(shù)據(jù)， 其中包括359 組雜草密度數(shù)據(jù)及67 組雜草多樣性數(shù)據(jù)。

1.2 數(shù)據(jù)分類

考慮到秸稈還田對草害發(fā)生的影響可能會受到其他條件因子的控制， 因此對收集到的相關(guān)因子信息進行歸納分組， 整理得到以下12 種影響因素（表1）： 氣候類型、年平均降水量、年平均氣溫、土壤質(zhì)地、土壤有機質(zhì)、土壤酸堿性、秸稈類型、還田方式、耕作方式、秸稈還田量、還田條件、作物類型， 并通過Meta 亞組分析比較不同因素下秸稈還田對草害發(fā)生的影響程度。

表1 試驗相關(guān)數(shù)據(jù)分類Table 1 Classification of experiment data

1.3 整合分析

本研究采用隨機效應(yīng)模型評估秸稈還田對我國農(nóng)田雜草的影響， 以處理組（秸稈還田）農(nóng)田的草害參數(shù)Xt與對照組（秸稈不還田）農(nóng)田的草害參數(shù)Xc的比值為響應(yīng)比， 并以響應(yīng)比（R）的自然對數(shù)作為Meta 分析的效應(yīng)值（lnR）， 具體公式如下[15］：

式中：Xt和Xc代表秸稈還田處理和秸稈不還田處理下雜草密度（株·m-2）和Shannon 多樣性指數(shù)參數(shù)。

通過Meta 分析得到綜合效應(yīng)值， 為便于理解和展示， 通過公式（2）計算秸稈還田對我國主要糧食作物草害效應(yīng)的變化百分?jǐn)?shù)：

式中：lnRR為處理組相對于對照組草害增加或減少的變化百分?jǐn)?shù)。

使用R 語言中的“metafor”軟件包進行發(fā)表偏倚檢驗。采用Egger test 法或羅森塔爾失安全系數(shù)法， 若P＜0.05 或Nfs 系數(shù)＞5n+10 （n為數(shù)據(jù)量），則認(rèn)為不存在發(fā)表偏倚[16-17］。總體效應(yīng)大小是用隨機效應(yīng)模型的加權(quán)重采樣方法計算的。如果95%的置信區(qū)間在0.05 的水平上不與0 重合， 則效果顯著為正或負(fù)。如果類別中95%的置信區(qū)間不重疊， 則兩類效應(yīng)大小的差異是顯著的。如果文獻的重復(fù)數(shù)大于3， 但是原始報告中沒有給出標(biāo)準(zhǔn)差（standard deviation， SD）或標(biāo)準(zhǔn)誤（standard error， SE）時， 用R 語言（version 3.6.1）的“metagear”數(shù)據(jù)包計算出標(biāo)準(zhǔn)差。對于每個屬性， 總異質(zhì)性（QT）分為組內(nèi)變異（QW）和組間變異（QB）。QB的顯著性表示分類亞組之間的平均效應(yīng)量的差異顯著水平。

1.4 統(tǒng)計分析

整合分析研究使用MetaWin 2.1 軟件對隨機效應(yīng)模型進行加權(quán)重采樣， 計算總體效應(yīng)大小， 并針對不同環(huán)境條件及農(nóng)業(yè)管理措施進行亞組效應(yīng)值分析來評估不同條件下秸稈還田對草害發(fā)生的影響程度[18-19］。試驗點分布圖、羅森塔爾失安全系數(shù)檢驗及相對重要性分析于R 4.1.2 中實現(xiàn)， 森林圖及線性回歸分析圖在Origin 9.0 中完成。

2 結(jié)果與分析

2.1 數(shù)據(jù)收集情況及秸稈還田對農(nóng)田草害的總體影響

本研究共收集到符合要求的論文41 篇， 獲得426 組觀測數(shù)據(jù)， 所選研究涉及11個省份， 大多位于我國東部地區(qū)（圖1a）。經(jīng)Egger tests 和Nfs 分析得到： 本研究羅森塔爾失安全系數(shù)分別為518952690 和2175， 其統(tǒng)計學(xué)意義為至少需要518952690 組具有發(fā)表偏倚性的雜草多度數(shù)據(jù)和2175 組具有發(fā)表偏倚性的雜草多樣性數(shù)據(jù)才能改變本試驗的研究結(jié)論， 即本研究所用數(shù)據(jù)不存在偏倚。由圖1b 可知， 與秸稈不還田（CK）相比， 秸稈還田可顯著降低農(nóng)田的雜草多度和雜草多樣性， 抑制效應(yīng)分別為-42.4%和-8.5%（P＜0.05）。

圖1 秸稈還田對農(nóng)田草害的整體影響Fig.1 Overall effects of straw returning on farmland weeds

2.2 不同條件下秸稈還田對農(nóng)田雜草多度的影響

由圖2 可知： 在亞熱帶季風(fēng)氣候區(qū)和溫帶季風(fēng)氣候區(qū)秸稈還田均可顯著降低雜草多度（P＜0.05）， 而溫帶大陸性氣候區(qū)的抑制效應(yīng)不顯著； 試驗地土壤的有機質(zhì)含量和酸堿性也會影響秸稈還田對田間雜草多度的抑制效果， 土壤有機質(zhì)含量為15～25 g·kg-1或土壤呈堿性時秸稈還田可顯著降低雜草多度（P＜0.05）。從秸稈類型上看， 麥稈還田可顯著降低雜草多度（P＜0.05）， 而稻稈、玉米稈及油菜稈還田對雜草多度的抑制效應(yīng)并不顯著；比較不同還田量對雜草多度的抑制效果發(fā)現(xiàn)， 以＞7000 kg·hm-2的秸稈量還田對雜草多度的抑制效果最好（P＜0.05）。秸稈無論是還于水田還是旱地均可顯著降低雜草多度（P＜0.05）， 而且水田環(huán)境下的抑制效果要明顯優(yōu)于旱地； 從作物類型上看， 種植水稻、玉米或油菜后秸稈還田能顯著降低雜草多度， 其中秸稈還于稻田對雜草的抑制效果最好（P＜0.05）。

圖2 不同條件下秸稈還田對農(nóng)田雜草多度的影響Fig.2 Effects of straw returning on farmland weed abundance under different conditions

2.3 不同條件下秸稈還田對農(nóng)田雜草多樣性的影響

相較于對雜草多度的影響， 秸稈還田對雜草多樣性的抑制效應(yīng)較弱， 其效應(yīng)值的顯著性和強度也會因氣候類型、土壤質(zhì)地等條件的不同而有所差別（圖3）。綜合不同條件下的效應(yīng)值可以發(fā)現(xiàn)： 將秸稈還田于溫帶季風(fēng)氣候區(qū)的砂質(zhì)土壤中可顯著降低雜草多樣性（P＜0.05）。此外， 將麥稈以＞7000 kg·hm-2的還田量翻埋還田后可顯著降低雜草多樣性（P＜0.05）。作物類型的不同也會影響秸稈還田對雜草多樣性的抑制效果， 秸稈還于玉米田或油菜田后雜草多樣性會顯著降低（P＜0.05）， 而還于水田或小麥田后的抑制效應(yīng)不顯著。另外， 通過比較不同耕作方式下秸稈還田對雜草多樣性的影響可以發(fā)現(xiàn)， 翻耕處理可顯著降低雜草多樣性（P＜0.05）， 而旋耕或免耕處理下秸稈還田對雜草多樣性并沒有顯著影響。

圖3 不同條件下秸稈還田對我國農(nóng)田雜草多樣性的影響Fig.3 Effects of straw returning on farmland weed diversity under different conditions

2.4 秸稈還田下農(nóng)田草害發(fā)生的主要變化及其影響因子

秸稈還田后農(nóng)田雜草的多度和多樣性會顯著降低， 相對重要性分析結(jié)果發(fā)現(xiàn)： 秸稈還田影響雜草多度的最關(guān)鍵因子是作物類型， 其次是秸稈還田量、年平均氣溫、年平均降水量、土壤pH、還田條件（圖4a）。秸稈還于不同作物農(nóng)田后對雜草多度的影響存在顯著差異， 其中秸稈還于稻田的雜草多度抑制效果最好， 而還于小麥田后卻不存在顯著抑制效果（圖2）。線性擬合結(jié)果顯示： 雜草多度效應(yīng)值與年平均氣溫、年平均降水量、土壤pH 均呈極顯著線性負(fù)相關(guān)（P＜0.01， 圖5a～c）， 即隨著年平均氣溫、年平均降水量及土壤pH 的增加， 秸稈還田對雜草多度的效應(yīng)值會逐漸降低。秸稈還田影響雜草多樣性的最關(guān)鍵因子是秸稈還田量， 其次是作物類型和土壤質(zhì)地（圖4b）， 其中雜草多樣性效應(yīng)值與秸稈還田量呈顯著線性負(fù)相關(guān)關(guān)系（R2=0.021，P＜0.05， 圖5d）， 即隨著秸稈還田量的增加， 秸稈還田對雜草多樣性的效應(yīng)值是逐漸降低的。秸稈還于不同作物農(nóng)田后對雜草多樣性的影響存在顯著差異， 其中秸稈還于油菜田的雜草多樣性抑制效果最好， 而還于水稻田或小麥田的抑制效應(yīng)均不顯著。此外， 秸稈還于砂質(zhì)土壤可顯著降低雜草多樣性， 而還于壤土或黏土后雜草多樣性并不能顯著降低（圖3）。

圖4 影響秸稈還田對農(nóng)田草害發(fā)生的預(yù)測因子的相對重要性分析Fig.4 Model-averaged importance of the predictors for the effects of straw returning on farmland weeds

圖5 關(guān)鍵因子與秸稈還田抑草效應(yīng)的關(guān)系Fig.5 The relationship between significant factors and effects of straw returning on weed suppression

3 討論

整合分析結(jié)果顯示秸稈還田可顯著降低農(nóng)田雜草多度， 但在不同條件下秸稈還田對雜草多度的抑制效應(yīng)存在差異， 本研究表明作物類型是影響該抑制效應(yīng)的最重要因子。已有研究指出種植不同作物農(nóng)田的雜草發(fā)生規(guī)律存在顯著差異[20］， 本研究結(jié)果與其一致。作物類型對抑草效應(yīng)的影響可能歸因于以下幾點： 一是不同作物產(chǎn)生的根系分泌物不同， 導(dǎo)致根際微生物群落結(jié)構(gòu)存在差異， 進而影響作物對相鄰雜草的抑制程度[21-22］； 二是種植不同作物農(nóng)田的土壤呼吸存在差異， 土壤溫度是引起上述差異的主要因素之一[23］。小麥生長季氣溫較低， 因此相較于其他作物， 小麥田的土壤呼吸速率較低， 秸稈的腐解效率下降， 進而降低了抑草效應(yīng)。除此之外， 不同作物農(nóng)田所使用的還田秸稈類型也是不同的， 而不同秸稈還田后所釋放的化感物質(zhì)及其腐解產(chǎn)物不同， 使得其對農(nóng)田雜草發(fā)生的影響也存在差異[24-25］。

除了作物類型之外， 秸稈還田量、年平均氣溫、年平均降水量對秸稈還田的雜草多度抑制效應(yīng)也存在顯著影響。秸稈還田量的差異不會改變雜草多度抑制效應(yīng)的顯著性， 但會影響其作用強度。本研究表明＞7000 kg·hm-2秸稈還田量的抑草效果最好， 這可能是由于較高的秸稈還田量可有效改善土壤養(yǎng)分和微生物群落結(jié)構(gòu)， 更有利于秸稈的腐解[26-27］。氣溫和降水量是影響土壤有機質(zhì)含量的主要因子[28］。相對較高的溫度和較大的降水量可能會從一定程度上改良微生物群落結(jié)構(gòu)并增加其生物活性[29］， 以加速秸稈的腐解與化感物質(zhì)的釋放[30］， 進而提高秸稈對雜草的抑制效應(yīng)。干燥少雨的環(huán)境下由于水熱不足， 可能會導(dǎo)致秸稈養(yǎng)分的凈釋放量降低， 最終使得抑草效應(yīng)有所降低[31］。

本研究結(jié)果表明秸稈還田也可以顯著抑制雜草多樣性， 而不同條件下的抑制效果存在差異， 其中秸稈還田量是最主要的影響因子。線性擬合結(jié)果表明， 隨著秸稈還田量的增加， 秸稈還田對雜草多樣性的抑制效應(yīng)值是逐漸降低的。秸稈還田量對該抑草效應(yīng)的影響可能與土壤碳氮比有關(guān)， 桂浩然等[11］研究發(fā)現(xiàn)玉米地秸稈還田量和土壤碳氮比呈正相關(guān)關(guān)系， 而隨著秸稈還田量的增加， 雜草多樣性逐漸降低。這可能是由于過高的碳氮比不能為具有較高多樣性的農(nóng)田雜草區(qū)系提供足夠的氮素。除了秸稈還田量之外， 作物類型和土壤質(zhì)地也是影響該抑制效應(yīng)的關(guān)鍵因子。徐琳等[32］研究發(fā)現(xiàn)， 種植不同作物耕地土壤的理化性質(zhì)及微生物群落結(jié)構(gòu)存在差異， 因此玉米田的雜草多樣性抑制效果最好可能是由于土壤有效養(yǎng)分含量更高， 微生物代謝活性更高， 進而導(dǎo)致秸稈腐解效率更高。然而秸稈還于水稻田后雜草多樣性并未發(fā)生顯著降低， 這可能是由于秸稈還田作為一項田間管理措施， 其產(chǎn)生的生態(tài)環(huán)境變化可能會抑制稻田淹水環(huán)境下原有水生雜草優(yōu)勢種群的生長[33］， 使得雜草群落結(jié)構(gòu)組成趨于平衡， 進而不利于秸稈還田對雜草多樣性的抑制。土壤質(zhì)地作為土壤重要的物理特性之一， 嚴(yán)重影響土壤的水肥氣熱狀況[34］。相較于壤土與黏土， 砂質(zhì)土壤水熱條件較差， 不利于雜草種子庫的存活[35］。

盡管秸稈還田對雜草多度和多樣性均有顯著的抑制效果， 但對比二者總體效應(yīng)值可以發(fā)現(xiàn)： 秸稈還田對前者的抑制效應(yīng)要顯著高于后者， 這在某種程度上說明： 雜草多度與多樣性對秸稈還田的響應(yīng)機制是不同的[36］。因此， 單獨秸稈還田防除雜草缺乏持續(xù)穩(wěn)定的控制作用， 需結(jié)合其他物化因素對農(nóng)田草害進行綜合防治[37］?？傮w效應(yīng)值存在較大差異可能是由于雜草多度和多樣性對各種條件因子的響應(yīng)機制不同。例如， 在亞熱帶季風(fēng)氣候條件下秸稈還田對雜草多度的抑制效果最好， 而雜草多樣性的抑制效應(yīng)在溫帶季風(fēng)氣候區(qū)最為顯著， 這可能是由于更溫暖的氣候下雜草物種的豐富度會更高[38］。因此， 相較于亞熱帶季風(fēng)氣候， 溫帶季風(fēng)氣候更有利于田間雜草多樣性的降低。另外， 秸稈還于水田對雜草多度的抑制效應(yīng)要明顯高于旱地， 而水田環(huán)境下秸稈還田對雜草多樣性的抑制效應(yīng)卻并不顯著。將秸稈覆蓋還于壤質(zhì)土壤雖然可以顯著降低雜草多度， 但不能使雜草多樣性顯著降低。此外， 不同耕作方式可通過擾動土壤結(jié)構(gòu)來改變土壤理化性質(zhì)及微生物特性， 進而影響農(nóng)田雜草的生長發(fā)育進程[39-41］。本研究主要涉及免耕、旋耕及翻耕3 種耕作方式， 其中免耕與旋耕翻動的土壤深度相對較淺， 從一定程度上限制了深土層雜草種子的出苗[42］， 減小了可出苗雜草的土層分布范圍， 進而提高了秸稈還田對雜草的集中防除效果。然而， 這兩種耕作方式下雜草多樣性并未顯著降低。這可能是由于保護性耕作對土壤的擾動較小， 有利于雜草群落內(nèi)不同生態(tài)位的保持， 因此該條件下農(nóng)田雜草多樣性會有所提升[43-44］。綜合上述差異可以看出， 雜草多度的顯著降低可能并不會伴隨著雜草多樣性的顯著降低， 這一研究結(jié)果說明在草害防治工作中應(yīng)注意雜草多度和多樣性的控制平衡。雜草多度在草害防治工作中是必須降低的指標(biāo)， 而保持一定的田間雜草多樣性或許有利于農(nóng)田養(yǎng)分的循環(huán)與資源的多級利用[45］。

4 結(jié)論

秸稈還田可有效抑制農(nóng)田草害的發(fā)生， 但是當(dāng)作物類型、秸稈還田量、氣候因子、土壤質(zhì)地等因素不同時抑制效應(yīng)會有所差別?？傮w來看： 在亞熱帶季風(fēng)氣候區(qū)將秸稈以＞7000 kg·hm-2的量還于稻田對雜草多度的抑制效果最好， 而在溫帶季風(fēng)氣候區(qū)以同等量的秸稈還于油菜田時可以更好地抑制雜草多樣性。值得注意的是， 秸稈還田抑草效應(yīng)的影響因素應(yīng)不僅限于本研究中所討論的， 可能的影響因素還有雜草類型（如單子葉和雙子葉）、土壤初始碳氮比、還田年限等。由于現(xiàn)階段相關(guān)文獻數(shù)據(jù)的缺乏， 目前暫時無法對其他因素進行整合分析研究， 因此進一步的研究可以嘗試從更多的角度去開展深入探討。