長期翻壓綠肥對提高灌漠土土壤肥力和作物產量的貢獻

吳科生，車宗賢*，包興國，張久東，曹衛(wèi)東，盧秉林，楊新強，楊蕊菊，崔 恒

（1 甘肅省農業(yè)科學院土壤肥料與節(jié)水農業(yè)研究所，甘肅蘭州 730070；2 農業(yè)農村部甘肅耕地保育與農業(yè)環(huán)境科學觀測實驗站，甘肅蘭州 733070；3 中國農業(yè)科學院農業(yè)資源與農業(yè)區(qū)劃研究所，北京 100081）

大量施用化肥是我國作物高產的重要支撐[1]，然而，長期大量施用化肥也帶來了土壤有機質下降，資源環(huán)境壓力加大，面源污染嚴重，土壤質量退化[2]等諸多問題，嚴重制約了土地和現代農業(yè)綠色可持續(xù)發(fā)展。通過輪作、間作、套種綠肥，翻壓綠肥作為農作物生長和培肥地力在我國有悠久的歷史，綠肥是將其植株翻壓，促進農作物生長和改善土壤理化性狀的綠色植物[3]。具有改良土壤[4]，降低土壤容重[5]，促進土壤團聚體形成和累積[6]，促進農作物生長發(fā)育以及觀賞、飼用等多方面的作用。因此，綠肥是農業(yè)綠色發(fā)展的有效物質支撐[7]?？茖W施用綠肥能減少化肥施用量，有效提高土壤肥力、改善生態(tài)環(huán)境[8]，對保障我國糧食安全和現代農業(yè)的綠色可持續(xù)發(fā)展意義重大。前人在施用綠肥能有效提高作物產量和肥料利用率[9–10]、提升農產品品質[11]、提高水肥利用效率[12–13]、改善土壤微生物群落結構[14]、施用綠肥替代或減施化肥[15–16]以及減少溫室氣體排放[17]等方面闡述了綠肥的施用效果和作用機理，但以上研究大多是基于綠肥的短期研究結果。灌漠土是河西綠洲灌區(qū)主要的土壤類型之一，灌漠土長期施用綠肥的土壤肥力演變規(guī)律和作物產量穩(wěn)定性等研究鮮見報道。本研究利用1988年在武威設置的灌漠土培肥長期定位田間試驗，從較長時間尺度上探索灌漠土長期施用綠肥條件下，土壤有機質、全氮、堿解氮、全磷、有效磷、全鉀、速效鉀等土壤肥力指標31年間的演變規(guī)律和綠肥肥料貢獻率，以及對作物產量的變化過程、穩(wěn)定性和可持續(xù)性進行分析，以期更科學地揭示長期施用綠肥對提高土壤肥力和增加作物產量的貢獻，為科學施用綠肥和綠肥培肥土壤提供理論依據。

1 材料與方法

1.1 研究區(qū)域概況

研究在甘肅省武威市涼州區(qū)永昌鎮(zhèn)白云村(38°37′N，102°40′E)，海拔 1504 m，年均溫 7.7℃，無霜期150天左右，降水量150 mm，蒸發(fā)量2021 mm。日照時數3028 h，相對濕度53%。一年中≥10℃的有效年積溫3016℃，年太陽輻射總量580～660 kJ/cm2，麥收后≥10℃的有效年積溫為1350℃[18]，屬典型的兩季不足，一季有余的自然生態(tài)區(qū)。區(qū)內降水量少蒸發(fā)量大，主要集中在6―9月，降水在農業(yè)生產中的可利用率較低。灌溉方式以地下水(井灌)灌溉為主。試驗地土壤為石灰性灌漠土，表土質地為輕壤，試驗開始時耕層(0—20 cm)土壤理化性質為：pH 8.8、有機質 16.4g/kg、全氮 1.06 g/kg、全磷 1.50 g/kg、全鉀 17.6 g/kg、堿解氮 64.5 mg/kg、有效磷 13.1 mg/kg、速效鉀 180 mg/kg、土壤容重 1.40 g/cm3[18]。

1.2 試驗設計與施肥方法

選擇長期定位試驗中的4個處理進行了分析：1) 綠肥磷肥配施 (綠肥 45000 kg/hm2+P2O5150 kg/hm2，GP)；2)綠肥氮磷配施 [綠肥 22500 kg/hm2+N 187.5 kg/hm2+P2O5150 kg/hm2，1/2 (G+N) P]；3)氮磷配施 (N 375 kg/hm2+P2O5150 kg/hm2，NP)；4)空白對照(CK)。小區(qū)長6.9 m，寬4.5 m，小區(qū)面積31.05 m2，各處理肥料施用量為每年總用量，所用綠肥為箭筈豌豆 (Viciasativa)，其有機碳 455 g/kg、N 36.9 g/kg、P 7.0 g/kg、K 30.3 g/kg，品種為隴碗 2 號，在播種前將綠肥粉碎翻壓在30 cm的耕作層內。小麥收獲前最后一次灌溉時在小麥行撒播套種箭筈豌豆，翻壓量中一部分是在小麥行種植的箭筈豌豆，不足部分用異地種植的相同品種綠肥補足，所用綠肥全部作基肥。氮肥為尿素，其中20%為底肥，在小麥播種前結合整地翻入土壤耕層；20%為小麥拔節(jié)期追肥，追施在小麥行；30%為玉米拔節(jié)期追肥，追施在玉米行；30%為玉米灌漿期追肥，追施在玉米行。磷肥用過磷酸鈣，施肥量150 kg/hm2，對照不施磷肥，其余處理作基肥一次性施入耕層[19]。

1.3 種植制度與栽培管理

試驗開始于1988年，栽培作物為小麥和玉米。種植模式是從1988年開始，以小麥/玉米間作―小麥單作―玉米單作為一個輪作周期種植。1999年以后至今，采用小麥/玉米間作的種植模式。小麥/玉米間作采用傳統的6∶2高產栽培模式，帶幅寬1.5 m，其中小麥帶寬70 cm，種6行小麥，玉米帶寬80 cm，種2行玉米。每小區(qū)設3個小麥玉米組合種植帶，中間一個組合種植帶為取樣和測產帶。試驗栽培作物為小麥和玉米，小麥品種為武春121或隴春系列，3月中旬播種，行距12 cm，播種量為600 kg/hm2，4月中上旬出苗，7月中旬收獲；玉米品種為中單系列或武科2號，4月下旬覆膜播種，株距20 cm，行距 40 cm，保苗 67500 株/hm2，5 月初出苗，10月上旬收獲。試驗全生育期共灌水6～7次[19]。

1.4 樣品采集與測定方法

每年作物收獲前在中間組合帶選取有代表性植株作為考種樣品，每個小區(qū)小麥取20株，玉米取10株，風干后測定其農藝性狀。取樣后小區(qū)剩余部分單打單收，將收獲的果穗晾曬風干后脫粒稱重。用烘干法測定作物籽粒含水量，采用標準含水量計算作物產量，測產面積為 10.35 m2(1.5 m×6.9 m)。每茬作物收獲后，在各小區(qū)用土鉆按“S”形路線取0―20 cm的耕作層土樣，土壤樣品在室內風干后按測定項目要求研磨過篩，進行土壤理化性質測定分析。pH采用電位法(土水比為1∶2.5)，有機質含量采用重鉻酸鉀氧化―油浴加熱法，全氮含量采用凱氏定氮法，全磷含量采取氫氧化鈉熔融―鉬銻抗比色法，全鉀含量采用氫氧化鈉熔融―火焰光度法，堿解氮含量采用堿解擴散法，有效磷采用Olsen法，速效鉀采用乙酸銨浸提―火焰光度法測定[20]。

1.5 計算方法與數據處理

產量穩(wěn)定性以統計學變異系數 (coefficient of variation，CV) 表示，用于衡量作物產量年際間的變異程度，CV 越大，說明產量穩(wěn)定性越低[21]。產量可持續(xù)性以可持續(xù)產量指數 (sustainable yield index，SYI)表示，用于表征產量可持續(xù)性，指數越高，說明該系統的可持續(xù)性越好[22]。

式中：σ 為標準差 (kg/hm2)；為平均產量 (kg/hm2)；Ymax為所有年份中產量的最大值(kg/hm2)。

肥料貢獻率 (fertilizer contribution rate，FCR)指施用肥料增加的作物產量占總產量的百分比，反映了年投入肥料的生產能力[23]。

式中：Y為施肥區(qū)作物產量；Y0為不施肥區(qū)作物產量。

采用 Microsoft Excel 2016 軟件對數據進行處理和繪圖，并用SAS 8.0軟件進行統計分析。

2 結果與分析

2.1 施用綠肥對土壤肥力的影響

2.1.1 土壤有機質含量和pH 施用綠肥提升了土壤有機質含量，隨綠肥施用年限的延長土壤有機質含量逐年增加(圖1)。各處理土壤有機質含量大小順序是，GP>1/2 (G+N) P>NP>CK。GP、1/2 (G+N) P、NP、CK歷年平均值較試驗初始值分別增加28.1%、21.1%、11.1%、4.8%。GP、1/2 (G+N) P較 NP 處理分別增加 15.3%、9.0%。2018 年 GP、1/2 (G+N) P、NP、CK比試驗初始數值分別提高68.8%、56.6%、30.9%、11.3%。GP、1/2 (G+N) P 較 NP 分別增加29.0%、19.6%。經方差分析，各處理歷年平均值GP、1/2 (G+N) P、NP 之間差異顯著，NP 與 CK 差異不顯著。

圖1 不同處理土壤有機質含量隨試驗年限的變化Fig.1 Variation of soil organic matter with fertilization years

2018年測定的土壤pH結果見圖2。1988年試驗初始值 pH 為 8.8，2018 年 GP、1/2 (G+N) P、NP、CK各處理pH分別為8.41、8.36、8.41、8.53，較試驗初始值分別降低了0.39、0.44、0.39、0.27。CK處理的pH顯著高于3個施肥處理，而綠肥氮磷配施土壤pH最小，相比初始值(8.8)降低幅度為5.0%。

圖2 不同施肥處理土壤pHFig.2 Soil pH under different fertilization treatments

2.1.2 土壤全氮和堿解氮含量 土壤全氮含量隨綠肥施用量的增加而增加，隨著施用年限的延長而上升 (圖3)。GP、1/2 (G+N) P、NP、CK 逐年平均值較試驗初始值分別增加30.4%、25.3%、10.3%、12.1%，GP、1/2 (G+N) P 較 NP 分別增加 18.2%、13.6%。2018 年，GP、1/2 (G+N) P、NP、CK 較試驗初始數值分別提高75.8%、56.3%、22.6%、46.2%，GP、1/2 (G+N) P 較 NP 分別增加 43.3%、27.4%。各處理歷年平均值GP、1/2 (G+N) P之間差異不顯著，NP與CK差異不顯著，前兩者與后兩者差異顯著。土壤堿解氮含量年際間變異較大(圖3)。GP、1/2 (G+N)P、NP、CK逐年平均值較試驗初始值分別增加36.4%、25.2%、14.9%、10.7%，GP、1/2 (G+N)P較NP分別增加18.8%、9.0%。2018年，GP、1/2(G+N) P、NP、CK比試驗初始值分別增加73.5%、53.2%、14.5%、7.1%，GP、1/2 (G+N) P 較 NP 分別增加51.5%、33.8%。各處理歷年平均值GP、1/2 (G+N)P、NP之間差異顯著，NP與CK差異不顯著。

圖3 不同施肥處理全氮和堿解氮含量隨試驗年限的變化Fig.3 Variations of soil total N and available N contents with fertilization years

2.1.3 土壤全磷和有效磷含量 土壤全磷含量隨種植年限的延長下降趨勢明顯(圖4)。與試驗初始值相比，只有GP處理全磷含量31年平均增加了0.9%，GP、1/2 (G+N) P 較 NP 分別增加了 4.9%、1.5%，在第 31 年 (2018 年)，GP、1/2 (G+N) P、NP、CK 處理的全磷含量比初始值分別降低了20.3%、27.8%、32.7%、48.4%，GP、1/2 (G+N) P 較 NP 分別增加了18.4%、7.2%。經方差分析，各處理歷年平均值GP、1/2 (G+N) P 之間差異不顯著，GP、NP、CK 之間差異顯著。土壤有效磷含量不同處理隨試驗年限的變化差異較大(圖4)。GP處理呈現上升趨勢，1/2(G+N) P處理穩(wěn)中有升，NP處理始終低于翻壓綠肥的兩個處理，而CK則呈現下降趨勢，GP、1/2 (G+N)P較NP分別增加了142.0%、54.1%。各處理歷年平均值GP、1/2 (G+N) P、NP處理的土壤有效磷含量較試驗初始值分別增加176.4%、76.1%、14.2%，CK降低了9.4%。以2018年與初始值比較，GP、1/2 (G+N) P、NP 分別增加 464.6%、260.3%、102.5%，CK 降低 70.0%，GP、1/2 (G+N) P 較 NP 增加178.8%、77.9%。經方差分析，各處理歷年平均值 GP、1/2 (G+N) P、CK 之間差異顯著，1/2 (G+N)P、NP差異不顯著。

圖4 不同施肥處理土壤全磷和有效磷含量隨試驗年限的變化Fig.4 Variations of soil total P and available P contents with fertilization years

2.1.4 土壤全鉀和速效鉀含量 不同施肥處理的土壤全鉀含量變化規(guī)律基本一致，前10年為逐年增加，到1997年達到最高值，之后的7年(1997―2004 年)快速下降，保持相對穩(wěn)定(2004―2012年)，之后又降低。31年土壤平均全鉀含量GP、1/2(G+N) P、NP、CK 處理較試驗初始值分別增加了34.9%、33.1%、31.3%、34.2%。GP、1/2 (G+N) P 較NP分別增加了2.8%、1.4%。經方差分析，各處理歷年平均值 GP、1/2 (G+N) P、NP、CK 之間差異不顯著。2018 年，GP、1/2 (G+N) P、NP、CK 較試驗初始數值分別提高16.1%、18.0%、14.2%、16.1%，GP、1/2 (G+N) P 較 NP 分別增加 1.7%、3.3%。不同處理的土壤速效鉀含量在年際間變化規(guī)律不明顯(圖5)。31年平均值中，只有GP較試驗初始值增加 2.8%，1/2 (G+N) P、NP、CK 分別降低 13.7%、31.3%、17.9%，GP、1/2 (G+N) P 較 NP 處理分別增加49.6%、25.7%。2018年速效鉀含量與初始值比較均有所下降，GP、1/2 (G+N) P、NP、CK 分別降低15.0%、37.2%、52.0%、29.8%，GP、1/2 (G+N) P較NP分別增加了77.2%、30.9%。經方差分析，各處理歷年平均值 GP、1/2 (G+N) P、NP、CK 之間差異顯著。

圖5 不同施肥處理土壤全鉀和速效鉀含量隨試驗年份的變化Fig.5 Variations of soil total K and available K contents with fertilization years

2.2 施用綠肥對小麥/玉米間作體系作物產量的影響

長期施用綠肥小麥產量相對其他處理較高(表1)。GP、1/2 (G+N) P、NP、CK 連續(xù) 28 年小麥平均產量分別為 4109、4052、3804、1713 kg/hm2。GP 和 1/2(G+N) P差異不顯著，但二者與NP差異顯著，分別較NP增加305、248 kg/hm2，增幅分別為8.3%、6.8%。相關分析表明，小麥產量與土壤全氮(R2=0.577，P<0.05)、土壤全磷(R2=0.838，P<0.01)和土壤有效磷(R2=0.715，P<0.05)呈顯著或極顯著正相關(表2)。

表1 不同處理小麥/玉米間作體系歷年平均產量Table 1 Average yield of wheat and maize intercropping system with different treatments over the years

施用綠肥對玉米產量影響較大(表1)。GP、1/2(G+N) P、NP、CK連續(xù)27年玉米平均產量分別為7578、9081、9116、3762 kg/hm2。NP 和 1/2 (G+N)P差異不顯著，但二者與GP差異顯著，分別較GP增加1538、1503 kg/hm2，增幅分別為20.3%、19.8%。相關分析表明，玉米產量與土壤全磷(R2=0.772，P<0.01)和土壤有效磷(R2=0.612，P<0.05)呈顯著正相關 (表2)。

長期綠肥化肥配施對小麥/玉米間作體系產量影響最大 (表1)。GP、1/2 (G+N) P、NP、CK 連續(xù)24年小麥/玉米間作體系平均產量分別為10938、12690、12363、5032 kg/hm2。NP 和 1/2 (G+N) P 差異不顯著，但二者與GP差異顯著，分別較GP增加1425、1752 kg/hm2，增幅分別為13.0%、16.0%。相關分析表明，間作體系總產量與土壤全磷(R2=0.812，P<0.01)、土壤有效磷(R2=0.659，P<0.05)呈極顯著或顯著正相關(表2)。

表2 土壤養(yǎng)分特征與作物產量的相關關系(R2）Table 2 Correlation coefficients of soil nutrient characteristics and yield of crops

2.3 施用綠肥的作物產量穩(wěn)定性、可持續(xù)性和肥料貢獻率

不同作物和種植方式下，作物產量的變異系數(CV)差異較大(表3)。小麥產量的變異系數大小順序為 CK>NP>1/2 (G+N) P>GP，兩個綠肥配施化肥(GP、1/2 (G+N) P) 處理的小麥產量穩(wěn)定性好于NP處理，又大大好于CK。玉米產量變異系數大小順序為 CK>GP>1/2 (G+N) P> NP，不施氮肥 (GP 和CK處理)大大降低了玉米產量的穩(wěn)定性。權衡小麥/玉米間作體系產量，其CV大小為CK>GP>NP>1/2(G+N) P，說明整個體系在綠肥配合氮磷肥下獲得的產量最穩(wěn)定，單施氮磷肥獲得的產量穩(wěn)定性好于不施氮肥的GP處理。

表3 施肥對作物產量穩(wěn)定性(CV)、可持續(xù)性(SYI)和肥料貢獻率(FCR)的影響Table 3 Effects of fertilization on crop yield stability (CV), sustainability (SYI) and fertilizer contribution rate (FCR)

作物的可持續(xù)性指數(SYI)不論小麥、玉米還是小麥/玉米間作體系，均以1/2 (G+N) P處理最大，不施肥CK的可持續(xù)性最低。就小麥來講，肥料的貢獻率以GP和1/2 (G+N) P處理較高，換句話說，翻壓綠肥提高了肥料對小麥產量的貢獻率；而就玉米或者小麥/玉米間作體系來講，施用氮肥提高肥料貢獻率的作用更重要。

3 討論

3.1 施用綠肥對土壤肥力的影響

張久東等[24]在甘肅河西綠洲灌區(qū)的研究結果顯示，長期減施30%化肥配合復種毛葉苕子，土壤有機質和堿解氮含量分別提高了1.8 g/kg和10.5 mg/kg。李繼明等[25]在紅壤上26年的定位試驗結果表明，綠肥與化肥長期配合施用可減少化肥用量，化肥氮、磷、鉀養(yǎng)分的農學利用效率高達60%以上。王勁松等[26]在紅壤上也證明了綠肥和化肥配施能顯著增加土壤有機質含量，改良土壤結構。在本研究中，小麥玉米間作體系下長期復種并翻壓箭筈豌豆能顯著提高土壤有機質、全氮、全鉀含量。其中綠肥磷肥配施土壤有機質含量歷年平均值和2018年末值，較試驗初始值分別提高28.1%和68.8%，較氮磷配施分別提高15.3%和29.0%。說明綠肥磷肥配施能顯著提高土壤有機質含量，而且隨著施用年限延長呈逐年連續(xù)增加趨勢。綠肥氮磷配施歷年平均值和2018年末值較試驗初始值、氮磷配施分別提高21.1%和56.5%、9.0%和19.6%?？梢钥闯鼍G肥氮磷配施和綠肥磷肥配施均能顯著提高土壤有機質含量。由此，長期施用綠肥對提高土壤肥力的貢獻可歸納為3個方面。一是施用綠肥土壤有機質含量顯著提高；二是隨綠肥施用年限延長土壤肥力持續(xù)增加；三是隨著綠肥施用量增加培肥土壤的貢獻增大。長期氮磷配施土壤有機質含量略有提高，空白對照土壤有機質含量呈下降或基本保持穩(wěn)定狀態(tài)。原因是長期施用化肥和不施肥土壤有機碳主要來源于作物枯枝落葉、根茬、以及土壤微生物和小動物，沒有其他外源有機物料輸入，長期施用化肥的作物生物量較大，留在土壤中的枯枝落葉和根茬相對較多，不施肥對照則較少，因此土壤有機質含量降低。綠肥磷肥配施和綠肥氮磷肥配施土壤全氮和堿解氮含量顯著高于氮磷配施和空白對照處理，原因是化學氮肥多為無機態(tài)，施入土壤后容易發(fā)生揮發(fā)和硝酸鹽淋失[27]，而有機態(tài)綠肥養(yǎng)分釋放較慢，肥料利用率高。

綠肥與化肥長期配合施用對土壤全磷有累積效應[25]。本研究結果表明，土壤全磷含量逐年降低。主要原因是1988年以前的種植制度是單作小麥，之后以小麥玉米間作為主，間作相對于單作作物生長周期長、產量高，磷素攜出量高，而施磷量一直不變，導致土壤全磷含量逐年降低。綠肥磷肥配施和綠肥氮磷配施土壤全磷含量歷年平均比氮磷配施分別增加4.9%和1.5%，2018年末值較氮磷配施處理分別增加18.4%和7.2%。說明施用綠肥能提高土壤全磷含量，并隨綠肥施用量的增加而增加。另外，土壤有效磷含量顯著提高。原因是石灰性土壤中CaCO3含量較高，有效磷施入土壤后易被Ca2+固定而轉化為緩效態(tài)磷[28]。施用綠肥后增加了土壤有機質，在其轉化過程中能產生有機酸等物質，可減少土壤對速效磷酸鹽的固定，促進了無機磷的溶解[29–30]。

鄧小華等[31]、陳曉波等[32]研究表明，綠肥翻壓還田能提高土壤全鉀含量1.24%～5.15%，能顯著提高土壤速效鉀含量335%以上。王勁松等[26]研究表明，綠肥和化肥配施能明顯增加土壤中速效鉀含量，維持土壤鉀素平衡。本研究中全鉀含量先增加后降低。其原因是施用綠肥補充了鉀素，導致土壤鉀素增加。而隨著試驗種植年限的延長，間作體系鉀素攜出量大，導致土壤全鉀含量逐年降低。但歷年平均值和2018年末值都比試驗初始值高，綠肥磷肥配施高于綠肥氮磷肥配施或者氮磷配施，說明土壤中鉀素含量仍處于較高水平。綠肥磷肥配施和綠肥氮磷肥配施土壤速效鉀含量較氮磷配施增加25.7%～77.2%。原因是不施用外源化學鉀肥條件下，綠肥施入土壤后會帶入一部分鉀素，而且不同處理作物的生物量和經濟產量不同，導致鉀素的攜出量不同，進而影響了土壤鉀含量。

3.2 施用綠肥對作物產量的影響

作物產量是衡量肥料效果的一項重要指標。有研究表明，長期綠肥與化肥配施能增加作物產量，提高土壤微生物量及酶活性[33–34]。本研究結果表明，在本試驗條件下，施用綠肥對小麥產量增產貢獻最大，對玉米產量增產貢獻小于化學氮肥。小麥/玉米間作體系總產量的歷年平均值以綠肥氮磷配施最高，綠肥磷肥配施較低。原因是小麥生物量比玉米小，需要吸收的土壤養(yǎng)分相對較少。在本研究綠肥施用量下，能夠滿足小麥全生育期的養(yǎng)分需求，加上綠肥顯著的培肥累積效應疊加，小麥產量顯著高于氮磷配施處理。相反玉米生物量大，僅施綠肥不能充分滿足玉米全生育期對氮素的需求，玉米生長后期出現脫肥現象，導致產量低于氮磷配施。為探索長期施用綠肥對土壤肥力和作物產量的累積效應，進一步對2018年土壤養(yǎng)分和作物產量進行相關分析，結果表明，小麥產量與土壤全氮、全磷、有效磷呈顯著正相關，玉米產量、間作體系總產量與土壤全磷和有效磷呈顯著正相關。可見，本試驗的綠肥施用量對小麥/玉米間作體系作物產量影響較大。在生產中可以考慮適當減少綠肥用量，增加氮磷施用量的綠肥化肥協同配施提高作物產量。

3.3 施用綠肥作物產量的穩(wěn)定性、可持續(xù)性和肥料貢獻率

作物產量穩(wěn)定性是判斷農田生態(tài)系統質量好壞的重要指標[35]。本研究表明，小麥產量穩(wěn)定性最好的是綠肥磷肥處理，其次是綠肥氮磷肥配施。玉米產量穩(wěn)定性最好是氮磷配施，綠肥氮磷配施次之。這與王婷等[36]的研究結果一致。小麥/玉米間作體系中各處理穩(wěn)定性均較好，其中綠肥氮磷配施作物產量穩(wěn)定性好于綠肥磷肥配施，說明綠肥化肥配施是保持作物產量穩(wěn)定性的更好措施。

可持續(xù)性(SYI)是衡量農田生態(tài)系統是否能持續(xù)生產的一個重要參數[37]。SYI值大于0.55時產量可持續(xù)性為較好，0.45～0.55為一般，小于0.45為較差[38]。在相同降水年型，氮肥+磷肥、秸稈+氮肥+磷肥冬小麥26年種植期產量可持續(xù)性均高于CK[36]。本研究得出了相似的結果，綠肥氮磷配施小麥、玉米產量及總產量的SYI值分別為0.36、0.49、0.62，綠肥磷肥配施分別為0.33、0.37、0.54，空白對照僅為0.14、0.28、0.48。說明小麥、玉米在綠肥磷肥配施、綠肥氮磷配施下產量可持續(xù)性一般，而小米/玉米間作體系總產量的可持續(xù)性均較好，尤其綠肥氮磷配施的產量可持續(xù)性最好，SYI為0.62。原因是不同處理對小麥和玉米產量影響不同。綠肥磷肥配施小麥產量最高，玉米產量較低，氮磷配施玉米產量最高小麥產量較低，小麥/玉米間作體系總產量是小麥玉米產量之和，對體系總產量而言，綠肥氮磷配施和氮磷配施處理都表現出較好的產量可持續(xù)性，而且綠肥氮磷配施的可持續(xù)性好于綠肥磷肥配施。

綠肥磷肥配施對小麥、小麥/玉米間作體系的肥料貢獻率較高，分別為0.58、0.53。這也再次證明了綠肥磷肥配施對小麥生產的貢獻主要來自于綠肥的高肥料貢獻率。綠肥氮磷配施對小麥、玉米、小麥/玉米間作體系的肥料貢獻率分別為0.58、0.57、0.59，氮磷配施分別為0.55、0.58、0.59?？梢?，綠肥氮磷配施的肥料貢獻率最大。這與樊廷錄等[39]在黃土高原黑壚土進行的24年定位試驗結果一致，原因是綠肥養(yǎng)分全面且釋放慢，而化肥養(yǎng)分單一而見效快，綠肥氮磷配施能適時適量和長時間為作物提供養(yǎng)分，更適宜促進作物的健康生長和產量形成。

4 結論

在河西灌漠土區(qū)，小麥、玉米產量與土壤有機質、全氮和有效磷含量顯著相關。與單施化肥相比，長期施用綠肥可顯著提升土壤有機質、全氮及有效磷含量，顯示出較好的土壤持續(xù)培肥效果。綠肥磷肥或者綠肥氮磷肥配施提高間作體系中小麥的產量但降低玉米的產量，肥料貢獻率、產量穩(wěn)定性和可持續(xù)性的效果好于單施化肥，綠肥與氮磷肥配合施用的增產貢獻、產量穩(wěn)定性、可持續(xù)性最好。從提高土壤有機質、氮、磷、鉀含量，增加作物產量等方面綜合考慮，綠肥和化肥配合施用對提升土壤肥力和增加作物產量的貢獻最大。因此，在本試驗綠肥施用量的基礎上，適當減少綠肥用量、增加氮磷施用量的綠肥化肥協同配施模式是在甘肅河西綠洲灌區(qū)灌漠土上實現培肥土壤和增加產量雙贏的施肥模式。