晉西黃土區核桃－谷子間作系統水肥影響域

王丹，畢華興，2，3，許華森，常譯方，蔡智才

（1.北京林業大學水土保持學院，北京100083；2.山西吉縣森林生態系統國家野外科學觀測研究站，北京100083；3.北京林果業生態環境功能提升協同創新中心，北京102206；4.中國農業大學資源與環境學院，北京100094）

晉西黃土區核桃－谷子間作系統水肥影響域

王丹1，畢華興1，2，3，許華森4，常譯方1，蔡智才1

（1.北京林業大學水土保持學院，北京100083；2.山西吉縣森林生態系統國家野外科學觀測研究站，北京100083；3.北京林果業生態環境功能提升協同創新中心，北京102206；4.中國農業大學資源與環境學院，北京100094）

為了定量判定晉西黃土區果農間作系統水肥主要競爭區域，以核桃－谷子間作系統為研究對象，對土壤水分和土壤養分進行定位監測，采用移動窗口法定量判定間作系統中土壤水分和土壤養分的影響域。結果顯示：當窗口寬度為2時，移動窗口法能有效地判定出土壤水分和土壤養分的影響域。40～100 cm土層，距樹行不同距離處的土壤水分和土壤養分含量相異性較小，難以判定影響域；綜合0～40 cm、0～100 cm土層的土壤水分狀況，土壤水分的影響域范圍是距樹南側1.75m到距樹北側1.25m，影響寬度為3m，即主要競爭區域為距樹南側1.75 m到距樹北側1.25m；不同的土壤養分指標影響域的范圍不同，綜合考慮土壤養分的四個不同指標（有機質、全氮、速效磷和速效鉀）及其在0～40 cm、0～100 cm兩個土層的分布情況，土壤養分的影響域是距樹南側2.25m到距樹北側2.25m，影響寬度為4.5m，即主要競爭區域為距樹南側2.25 m到距樹北側2.25 m。建議在土壤水分和土壤養分的主要競爭區域，采取改變壟溝、起壟覆膜、選擇性施肥等措施，緩解核桃、谷子之間的競爭，提高資源利用效率，增加當地農民經濟收入。

核桃－谷子間作系統；土壤水分；土壤養分；影響域；主要競爭區域；晉西黃土區

農林復合是生態農業的主要內容之一。果農間作作為一種典型的農林復合模式，因其具有較高的生態效益、經濟效益和社會效益［1－2］，受到廣泛的推廣種植。果樹和農作物在同一土地單元內共享自然資源，會引起一定程度的競爭［3－4］，尤其是資源不足以滿足果樹和農作物的共同生長需求時。這種資源競爭主要包括地上部分光的競爭［5－7］以及地下部分土壤水分［8－10］和土壤養分的競爭［9－13］。雖然現階段國內外大量研究學者開展了對農林復合的相關研究，取得了一定的成果，但是大部分集中在光、土壤水分、土壤養分、根系［14－16］等空間分布特征的研究，缺乏對果農間作系統中地下部分競爭水平變化的定量判定。由于谷子根重90%以上集中在0～40 cm土層［17］，核桃根系主要分布在0～100 cm土層［18－21］，本研究將土層分為0～40、40～100 cm和0～100 cm三個層次對核桃－谷子間作系統中土壤水分和土壤養分進行定位監測，利用移動窗口法，從水分、養分競爭的角度分析核桃與谷子之間的關系，定量判定土壤養分、土壤水分的影響域，劃分核桃－谷子間作系統的主要競爭區域，以期為農林復合系統的高效可持續經營提供理論依據。

1 研究區概況

研究區位于山西省吉縣（35°53′10″～36°21′02″N，110°27′30″～111°07′20″E），地處黃河中游東岸，屬于典型的黃土高原殘垣溝壑區。吉縣土壤以褐土為主，黃土母質，土層深厚。年均降水量522.8 mm，降水具有季節變化，4—10月降水量占全年降水量的90%以上。吉縣屬于典型的暖溫帶大陸性季風氣候，日照時數2 538 h，大于10℃的有效積溫3 361.5℃。無霜期年平均172 d，年均氣溫10.2℃，年均日較差11.5℃。常見的有木本植物193種，49科；草本植物180種，44科（不包括農作物）；各類中草藥約141種。果農間作系統中主要經濟林樹種包括蘋果（Malus pumila）、核桃（Juglans regia）等，主要農作物包括大豆（Glycine max）、花生（Arachis hypogaea）、谷子（Setaria）、苜蓿（Medicago sativa）等。

2 研究方法

2.1 研究對象

根據晉西黃土區果農間作的經營特點以及實際情況，在山西省吉縣東城鄉柏東村建立核桃－谷子間作系統試驗示范區。試驗區核桃品種為晉龍1號，種植于2009年，樹齡為5 a，株行距為6m×6m，樹行走向為東西向；谷子品種為晉谷21，種植于2014年4月7日，距樹0.4 m開始種植作物，株距0.3m×0.4 m。在試驗區，以中央核桃帶為中心從南至北布設樣線，樣線與核桃樹行向垂直。以距核桃樹行南側3 m處為第一個樣點，每隔0.5 m取一個樣點，直到距樹北側3m，每條樣線布設13個樣點（圖1），共布設3條重復樣線。

圖1 核桃－谷子間作土壤水分、土壤養分定位監測試驗布設圖Fig.1 Experiment layout of monitoring soil moisture and nutrients in walnut－millet intercropping systems

2014年8月，在谷子開花期，選擇晴朗的典型天氣，在3條樣線的各監測點采用土鉆法取土，取土深度為100 cm，每20 cm為一層。利用烘干稱重法確定土壤水分的含量。土壤養分按照中華人民共和國林業行業標準進行測定［22］，具體方法如下：土壤有機質采用濃硫酸消煮－重鉻酸鉀氧化容量法；全氮采用濃硫酸消煮－凱氏定氮法；速效鉀采用醋酸銨浸提－火焰光度計法；速效磷采用碳酸氫鈉浸提－鉬銻抗比色法。測定3次，取其平均值。

2.2 數據處理與分析

本文的數據處理方法主要有移動窗口法［23］和Microsoft Excel 2010。采用移動窗口法判定核桃－谷子間作的影響域范圍的原理如圖2所示：將一個窗口（含偶數個取樣點）平均分割成兩個半窗口，利用歐氏距離的平方（SED）計算比較兩個半窗口的相異性，按順序將窗口向右側移動一個取樣點，再計算兩個半窗口的相異性，直到窗口包含最后一個樣點，即所有的取樣點都參與計算為止。利用歐氏距離的平方計算相異性的公式如下：

式中，SED是歐氏距離的平方；k是樣點的變量數；X是一個半窗口的均值；m和n分別為兩個半窗口；w是窗口的寬度。

以樣點位置為橫坐標，歐氏距離的平方為縱坐標作圖。當曲線出現兩個明顯且穩定的波峰時，根據陡峭的峰值寬度可以有效地判定核桃－谷子間作系統的影響域。本研究中，當窗口寬度為2時，以中心點代表區域范圍，則歐氏距離的平方值與距樹距離的對應關系如表1所示。影響域內距中心林帶不同距離處土壤水分和土壤養分的含量差異較大，即果樹與農作物之間的競爭比影響域外更加激烈，為主要競爭區域。

圖2 移動窗口法原理示意圖Fig.2 Schematic diagram of moving split window techniques

表1 歐氏距離平方值與距樹距離的對應關系Table 1 Correspondence between SED and the distance from the tree

3 結果與分析

3.1 土壤水分影響域

采用移動窗口法對土壤水分的影響域進行分析，窗口寬度分別設為2、4、6、8。窗口寬度為4、6、8時，曲線波動較大，波峰位置不明顯，無法判定影響域范圍和影響寬度。當窗口寬度為2時，有兩個明顯且穩定的波峰出現（圖3），說明在窗口寬度為2時能夠辨析土壤水分的影響域。0～40 cm土層，在SED5和SED8有兩個明顯的波峰，峰寬分別是SED4－SED6和SED6－SED9，通過半波峰法判定影響域范圍是距樹南側1.25m到距樹北側1.25m，影響寬度為2.5 m。0～100 cm土層，在SED4和SED8出現兩個波峰，由此判定影響域范圍是距樹南側1.75 m到距樹北側1.25 m，影響寬度為3 m。40～100 cm土層，在SED6和SED9有兩個明顯的波峰，但兩個半窗口之間的相異系數較小，無法準確地判斷影響域范圍。以上結果表明，綜合考慮各個土壤層次，土壤水分的影響域范圍是距樹南側1.75 m到距樹北側1.25m，影響寬度為3m；即核桃－谷子間作系統中，土壤水分的主要競爭區域是距樹南側1.75 m到距樹北側1.25 m。

圖3 窗口寬度為2時核桃－谷子間作系統土壤水分的相異性Fig.3 SED of soil moisture in walnut－millet intercropping systems when window width was2

3.2 土壤養分影響域

采用移動窗口法對土壤養分的不同指標（有機質、全氮、速效磷和速效鉀）進行分析，曲線變化情況見圖4。與土壤水分相同，當窗口寬度為4、6、8時，曲線沒有明顯且穩定的波峰，土壤養分的影響域范圍和影響寬度難以判定。圖4是窗口寬度為2時土壤養分的影響域，隨著窗口向右移動，曲線變化趨于穩定，有兩個明顯且穩定的波峰，能夠準確判定影響域。

核桃－谷子間作系統中土壤有機質在0～40 cm土層的影響域范圍是距樹南側2.25 m到距樹北側2.25 m，影響寬度為4.5 m，波峰分別出現在SED4和SED10。0～100 cm土層在SED4和SED9有兩個明顯的波峰。影響域范圍分別是距樹南側1.75 m到距樹北側2.25 m，影響寬度為4 m。40～100 cm土層，距中心林帶不同距離的土壤有機質含量差異并不顯著，無法準確量化影響域范圍（圖4a）。

核桃－谷子間作系統中土壤全氮在表層土壤0～40 cm相異性較大，SED6和SED8歐氏距離的平方值達到最大，呈現出兩個穩定波峰，其影響域范圍是距樹南側1.25 m到距樹北側1.25m，影響寬度為2.5 m。0～100 cm土層，影響域范圍是距樹南側1.75m到距樹北側2.25 m，影響寬度為4 m。40～100 cm土層，曲線變化趨于穩定，有兩個穩定的波峰，但是與0～40 cm和0～100 cm相比，距中心林帶不同距離處土壤全氮含量相異性較小，波峰位置不明顯，難以判定影響域（圖4b）。在影響域范圍以外，中心林帶南側有小波峰出現，說明此處取樣點的全氮含量有較大的相異性，可能原因是該試驗示范區2013年種植具有固氮功能的大豆，大豆植株附近土壤全氮含量增加，使距樹行不同距離處土壤全氮的含量具有明顯差異。

核桃－谷子間作系統中土壤速效鉀在0～40 cm土層，波峰位置為SED4和SED8，影響域范圍是距樹南側1.75m到距樹北側2.25m，影響寬度為4m。0～100 cm土層，波峰出現在SED6和SED9，影響域范圍是距樹南側1.75m到距樹北側2.25m，影響寬度為4m（圖4c）。

同上所述，核桃－谷子間作系統中土壤速效磷在0～40、0～100 cm土層，均有兩個穩定的波峰。0～40 cm土層，在SED4和SED9出現波峰，影響域范圍是距樹南側2.25m到距樹北側2.25m，影響域寬度為4.5 m。0～100 cm土層，影響域范圍距樹南北兩側各2.25 m，影響域寬度為4.5 m。40～100 cm土層，有兩個穩定的波峰，但是與0～40 cm和0～100 cm相比，總體波動不明顯（圖4d）。

圖4 核桃－谷子間作系統土壤養分的相異性Fig.4 SED of soil nutrients in walnuts－millet intercropping systems

綜合考慮土壤養分的四個不同指標，核桃－谷子間作系統中土壤養分的主要競爭區域為距樹南側2.25m到距樹北側2.25 m。對比土壤水分和土壤養分的影響域，總體而言，除0～40 cm土層的土壤全氮外，在0～40 cm和0～100 cm土層，土壤水分的影響域寬度均小于土壤養分各指標的影響域寬度。土層深度為40～100 cm，兩個半窗口之間歐氏距離的平方明顯小于0～40 cm和0～100 cm土層，難以判定該土層土壤水分和土壤養分的影響域。這主要是因為核桃根系主要集中在0～100 cm土層［17］，谷子根重90%以上集中在0～40 cm土層［18－21］，核桃和谷子的競爭主要集中在0～40 cm土層。土壤養分的不同指標影響域范圍進行對比，全氮的影響寬度最小，可能原因是上季種植具有固氮作用的農作物大豆，增加了土壤全氮的含量，進而在一定程度上緩解了果農間作系統中全氮的競爭。

4 討論

影響域的判定方法主要有多元序列［24］、主成分分析［25］、聚類分析［26］、除趨勢對應分析（DCA）［27］以及移動窗口法［23］等。本文主要采用移動窗口法對影響域進行定量判定。李麗光［28］等利用移動窗口法研究發現岷江上游花椒地和林地邊界土壤水分影響域，王曉等［29］辨析了濕地－干草原生態交錯帶的寬度和邊界，石培禮等［30］發現了岷江冷杉林線的影響域范圍和寬度。雖然眾多研究者利用移動窗口法成功對影響域進行定量判定，但此方法仍然具有一定的局限性，從本文實際研究可知，只有窗口寬度為2時才能夠有效地判定影響域，即適當的窗口寬度才能準確地判斷影響域的范圍和影響寬度。

通過對核桃－谷子間作系統中水肥的影響域研究初步發現，在垂直方向上，40～100 cm土層基本不存在競爭，0～40 cm土層是土壤水分和土壤養分的主要競爭區域。由于谷子根重90%以上存在于0～40 cm土層［17］，所以本研究中將土壤水分和土壤養分的監測分為0～40、40～100 cm兩個土層，要想精確地判定核桃－谷子間作系統在垂直方向上的影響域以及主要競爭區域，需要更加精細地劃分土層進行進一步研究。隨著核桃與谷子生長、降雨等條件的動態變化，間作系統中土壤水分和土壤養分的含量、分布、影響域及主要競爭區域也不斷發生變化，而本文只是研究了谷子開花期（雨季）核桃－谷子間作系統的主要競爭區域，因此劃分多尺度、多時間、多空間、多類型果農間作系統中水肥影響域，并判定主要競爭區域，是果農間作系統高效可持續經營研究中一個重要的方向。

5 結論

在晉西黃土區，對主要競爭區域的水肥進行科學調控，是緩解果農間作系統中種間競爭、提高資源利用效率、增加當地農民收入關鍵技術之一。在研究區，農業生產完全依靠天然降水，水資源緊缺是制約果農間作系統中經濟效益的主要因素。根據本研究可知，土壤水分的主要競爭區域是距樹南側1.75 m到距樹北側1.25 m，建議在此區域改變壟溝進行降雨就地富集，引導非耕地地區徑流向此處富集以及起壟覆膜等措施，最大程度利用降水，提高水資源利用效率和降水生產力。同時，果農間作系統中土壤養分含量對作物的產量也具有較大影響［31］，因此根據土壤養分的主要競爭區域，建議在距樹南側1.75 m到距樹北側2.25m主要施用氮肥和鉀肥；在距樹南北兩側各2.25m施用磷肥；若施用肥料品種為氮磷鉀復合肥，建議在距樹南北兩側各2.25m進行施肥。有目的性的施肥，可以提高作物產量以及肥料利用率，同時防止過度施肥造成資源浪費和農田污染。

［1］許華森，畢華興，王若水，等.晉西黃土區蘋果農作物間作系統經濟效益［J］.中國水土保持科學，2014，12（4）：81-85.

［2］孟平，張勁松.太行山低山丘陵區果-草復合系統生態經濟效應研究［J］.中國生態農業學報，2003，11（2）：111-113.

［3］Ong C K，Corlett J E，Singh R P，et al.Above and below ground interactions in agroforestry systems［J］.Forest Ecology and Management，1991，45（1）：45-57.

［4］趙薈，王晶.復合農林系統種間關系研究進展［J］.林業建設，2011，（6）：34-39.

［5］田陽，云雷，畢華興，等.晉西黃土區果農間作光競爭研究［J］.水土保持研究，2013，20（4）：288-292.

［6］Peng X B，Zhang Y Y，Cai J，etal.Photosynthesis，growth and yield of soybean and maize in a tree-based agroforestry intercropping system on the Loess Plateau［J］.Agroforestry Systems，2009，76（3）：569-577.

［7］彭曉邦，蔡靖，姜在民，等.光能競爭對農林復合生態系統生產力的影響［J］.生態學報，2009，29（1）：545-552.

［8］Mcintyre B D，Riha SJ，Ong CK.Competition forwater in a hedgeintercrop system［J］.Field Crops Research，1997，52（1）：151-160.

［9］云雷，畢華興，任怡，等.黃土區果農復合系統種間水分關系研究［J］.水土保持通報，2008，28（6）：110-114.

［10］廖文超，畢華興，趙云杰，等.晉西蘋果＋大豆間作土壤水分分布及其對大豆生長的影響［J］.中國水土保持科學，2014，12（1）：24-28.

［11］蔡崇法，王峰，丁樹文，等.間作及農林復合系統中植物組分間養分競爭機理分析［J］.水土保持研究，2000，7（3）：219-221，252.

［12］高國治，王明珠，張斌.低丘紅壤南酸棗-花生復合系統物種間水肥光競爭的研究Ⅱ.南酸棗與花生利用光能分析［J］.中國生態農業學報，2004，（2）：92-94.

［13］Moreno G，Obrador JJ，Garcia A.Impact of evergreen oaks on soil fertility and crop production in intercropped dehesas［J］.Agriculture，Ecosystems＆Environment，2007，119（3）：270-280.

［14］云雷，畢華興，馬雯靜，等.晉西黃土區林草復合系統刺槐根系分布特征［J］.干旱區資源與環境，2012，26（2）：151-155.

［15］馬秀玲，陸光明，徐兆齡，等.農林復合系統中林帶和作物的根系分布特征［J］.中國農業大學學報，1997，2（1）：109-116.

［16］張勁松，孟平.果農復合系統中果樹根系空間分布特征［J］.林業科學，2002，38（4）：30-33.

［17］張永清.幾種谷類作物根土系統的研究［D］.太谷：山西農業大學，2005：8-9.

［18］馬長明，翟明普，劉春鵬.單作與間作條件下核桃根系分布特征研究［J］.北京林業大學學報，2009，31（6）：181-186.

［19］云雷，畢華興，馬雯靜，等.晉西黃土區核桃花生復合系統核桃根系空間分布特征［J］.東北林業大學學報，2010，38（7）：67-70.

［20］彭曉邦，張碩新.商洛低山丘陵區核桃-黃芩農林復合系統中植物根系分布特征［J］.西北農業學報，2013，22（9）：171-177.

［21］何春霞，孟平，張勁松，等.華北低丘山區核桃-決明子復合模式的根系分布［J］.林業科學研究，2013，26（6）：715-721.

［22］葛曉改，肖文發，曾立雄，等.不同林齡馬尾松凋落物基質質量與土壤養分的關系［J］.生態學報，2012，32（3）：852-862.

［23］石培禮，李文華.生態交錯帶的定量判定［J］.生態學報，2002，22（4）：586-592.

［24］Stohlgren T J，Bachand R R.Lodgepole pine（Pinus contorta）ecotones in rocky mountain national park，colorado，USA［J］.Ecology，1997，78（2）：632-641.

［25］許華森，云雷，畢華興，等.刺槐＋苜蓿復合系統土壤養分分布特征及邊界影響域——以晉西黃土區為例［J］.中國水土保持科學，2011，9（5）：48-53.

［26］Fortin M.Edge detection algorithms for two-dimensional ecological data［J］.Ecology，1994，75（4）：956-965.

［27］Choesin D，Boerner R E J.Vegetation boundary detection：a comparison of two approaches applied to field data［J］.Plant Ecology，2002，158（1）：85-96.

［28］李麗光，何興元，李秀珍，等.岷江上游花椒地／林地邊界土壤水分影響域的定量判定［J］.應用生態學報，2007，17（11）：2011-2015.

［29］王曉，張克斌，楊曉暉，等.半干旱區濕地-干草原交錯帶邊界判定及其變化［J］.生態學報，2012，32（16）：5121-5127.

［30］石培禮，劉興良.游動分割窗技術在生態交錯帶定量判定中的應用—以四川巴郎山岷江冷杉林線為例［J］.植物生態學報，2002，26（2）：189-194.

［31］高路博，畢華興，許華森，等.晉西幼齡蘋果×大豆間作的土壤中水分，養分空間分布特征及對大豆的影響［J］.中國農學通報，2013，29（24）：36-42.

Distance of edge influence of soil moisture and nutrients in walnut-millet intercropping systems in the Loess Plateau in west Shanxi Province

WANG Dan1，BIHua-xing1，2，3，XU Hua-sen4，CHANG Yi-fang1，CAIZhi-cai1
（1.College of Soil and Water Conservation，Beijing Forestry University，Beijing，100083，China；2.Ji county station，Chinese National Ecosystem Research Network（CNERN），100083，Beijing，China；3.The Collaborative Innovation Center for Promoting Ecological Functions of Forestry and Economic Trees，Beijing，102206，China；4.College of Resource and Environment，China Agricultural University，Beijing 100094，China）

In order to divide the main competition area for intercropping systems and provide a theoretical guidance for the sustainable and efficient management of intercropping systems on the Loess Plateau in west Shanxi Province，the walnut-millet intercropping system was studied.Soil moisture and nutrients in the intercropping systems were monitored.Moving split window technique was used to detect the distance of edge influence of soil moisture and nutrients in the intercropping systems.The results showed that when the window width reached 2，the distance of edge influence of soil moisture and nutrients could be detected by moving split window technique.It is difficult to determine the distance of edge influence in 40～100 cm soil layer because soil moisture and nutrients content at different distances from the tree row were quite similar.Combining soil moisture in 0～40 cm and 0～100 cm soil layers，distance of edge influence of soil moisture was from 1.75m away in the south to1.25m away in the north side of the tree，that was to say influencing width was 3m，which means the main competition area of soil moisture ranged from 1.75m south of the tree to 1.25m north of the tree.Different soil nutrient indicators had different distance of edge influence.Taking into account of the distribution of four soil nutrient indicators（organic matter，total N，available P and available K）in 0～40 cm and 0～100 cm soil layers，distance of edge influence of soil nutrients was 2.25m south away from the tree to 2.25m north away from the tree and the influencing width was4.5m，which means the main competition area was from 2.25m south of thetree to 2.25 m north of the tree.Changing ridge and furrow，covering ridge and fertilizing selectively in the main competition area of soil moisture and nutrients was proposed to ease the competition between walnuts and millet so as to improve efficiency of resource use and increase the income of local farmers.

walnut-millet intercropping systems；soil moisture；soil nutrients；distance of edge influence；the main competition area；the Loess Plateau in west Shanxi Province

S344.2

A

1000-7601（2016）06-0150-06

10.7606／j.issn.1000-7601.2016.06.23

2016-01-30

國家自然科學基金（31470638）

王丹（1990—），女，河南鶴壁人，在讀碩士研究生，研究方向為林業生態工程。E-mail：beilinwangdan＠126.com。

畢華興，教授，博士生導師，研究方向為林業生態工程。E-mail：bhx＠bjfu.edu.cn。