任媛媛 張 莉 郁耀闖 張彥軍 張歲岐

大豆種植密度對玉米/大豆間作系統(tǒng)產(chǎn)量形成的競爭效應分析

任媛媛1張 莉2郁耀闖1張彥軍1張歲岐3,*

1寶雞文理學院地理與環(huán)境學院 / 陜西省災害監(jiān)測與機理模擬重點實驗室, 陜西寶雞 721013;2濰坊醫(yī)學院藥學院, 山東濰坊 261053;3中國科學院水利部水土保持研究所 / 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 陜西楊凌 712100

旱區(qū)作物間作系統(tǒng)是優(yōu)化作物群體質(zhì)量、提高產(chǎn)量的一種重要種植方式, 從產(chǎn)量構(gòu)成因子角度探討作物間作競爭優(yōu)勢機制還鮮見報道。本研究設置3個大豆密度梯度和4個種植比例(玉米與大豆以2∶0、0∶2、2∶2和2∶4間作), 分析間作系統(tǒng)的作物競爭指數(shù)、產(chǎn)量構(gòu)成和最終產(chǎn)量的差異性變化, 探討間作群體產(chǎn)量增加的作物競爭機制。結(jié)果表明, 不同大豆密度和種植比例組合下的間作系統(tǒng)產(chǎn)量提高14%~23%。玉米的實際產(chǎn)量損失大于0, 大豆的實際產(chǎn)量損失小于0。間作系統(tǒng)中玉米的穗重、穗長、穗粗、穗粒重、軸重和千粒重均顯著高于單作; 除結(jié)莢長度和主莖節(jié)間長度外, 間作系統(tǒng)中大豆的單株粒重、單株莢數(shù)、單株總粒數(shù)、單株有效粒數(shù)、主莖節(jié)數(shù)和百粒重低于單作或與單作間無顯著差異。間作系統(tǒng)中玉米的競爭比率大于1, 大豆的競爭比率小于1, 在3種大豆密度下玉米和大豆的競爭比率分別為2.08、1.84、1.68和0.49、0.56、0.63, 表明隨著大豆密度的增加, 間作中玉米的競爭比率增加的同時大豆的競爭比率降低。玉米的侵占力大于0, 大豆的侵占力小于0。玉米產(chǎn)量與軸重、千粒重、穗重、穗粒重、穗長、行粒數(shù)呈正相關關系, 與禿尖長呈負相關關系。通徑分析表明, 直接作用中, 穗粒重對玉米產(chǎn)量的貢獻最大(2.18); 間接作用中, 軸重、千粒重通過每穗粒重對玉米產(chǎn)量的貢獻較大(1.64和1.58)。綜上所述, 大豆間作玉米有間作優(yōu)勢, 間作優(yōu)勢來源于每穗粒重。

密度; 間作; 通徑分析; 競爭比率

間作是指在一個生長季內(nèi)同一塊土地上同時種植2種或2種以上的作物。相比常規(guī)單作, 間作系統(tǒng)作物的群體結(jié)構(gòu)特性, 可顯著提高光能利用率[1], 提升根系養(yǎng)分吸收能力[2]; 作物互作效應可防治病蟲害[3], 并有效控制雜草[4]; 作物競爭可提高作物產(chǎn)量[5], 水分養(yǎng)分利用效率[6-7], 同時增加經(jīng)濟效益, 發(fā)掘土地資源生產(chǎn)潛力[8], 在現(xiàn)代農(nóng)業(yè)生產(chǎn)中已廣泛應用。其中, 豆禾間作由于豆科作物的固氮作用可有效的向禾本科作物轉(zhuǎn)移一定量氮素, 而禾本科作物可以提高豆科作物的生物固氮效率[9]。因此, 作物間作在農(nóng)業(yè)生態(tài)系統(tǒng)增產(chǎn)增效中具有舉足輕重的地位。

間作系統(tǒng)中由于作物對資源需求在時間和空間上的分離互補, 促進了不同作物對資源的高效利用[10]。同時, 間作系統(tǒng)中, 作物競爭能力差異及其變化特征受不同間作方式的直接影響[11-12]。種植比例和密度是影響間作系統(tǒng)產(chǎn)量的主要影響因素, Undie等[13]研究發(fā)現(xiàn), 玉米和大豆不同間作比例(1∶1、2∶2和1∶2)的系統(tǒng)產(chǎn)量分別高出同期單作64%、66%、63%和43%、57%、65%。玉米與花生以2∶2、4∶4、6∶6和8∶8比例間作, 其中2∶2間作栽培中玉米產(chǎn)量優(yōu)勢明顯[14]。豌豆與小麥、黑麥、黑小麥以80∶20間作的產(chǎn)量高于60∶40, 主要由于前者間作系統(tǒng)中大豆相對于其他作物的競爭能力高于后者[8]。杜進勇等[15]關于小麥和玉米間作的研究發(fā)現(xiàn), 高密度玉米間作系統(tǒng)產(chǎn)量顯著高于低密度, 主要是由于增加密度引起玉米光合勢和凈光合速率比低密度分別提高12.4%和11.2%。不同密度梯度下玉米和豌豆間作栽培研究發(fā)現(xiàn), 中密度間作系統(tǒng)產(chǎn)量分別高于低密度、高密度間作系統(tǒng)8.37%和9.09%, 相對應豌豆對玉米的競爭力較低密度、高密度提高了36.56%、20.17%[16]。因此間作系統(tǒng)中不同作物的空間配置方式影響間作系統(tǒng)中作物的競爭能力, 進而引起產(chǎn)量潛力的發(fā)揮。

間作系統(tǒng)中作物產(chǎn)量受到產(chǎn)量構(gòu)成因素的協(xié)同作用和相互影響[17-18], 研究間作系統(tǒng)產(chǎn)量構(gòu)成因子對作物競爭的響應機制, 及其與最終產(chǎn)量形成的相互關系是明確作物間作調(diào)節(jié)產(chǎn)量構(gòu)成因子、提升作物產(chǎn)量的關鍵內(nèi)容[18]。合理的間作比例和密度搭配可以挖掘產(chǎn)量構(gòu)成要素的競爭增產(chǎn)優(yōu)勢, 利于作物營養(yǎng)物質(zhì)轉(zhuǎn)移積累、提高作物經(jīng)濟產(chǎn)量。目前, 針對不同間作系統(tǒng)對產(chǎn)量[19]、水分利用[6,20]和養(yǎng)分利用[10,21]等研究已有大量報道{Raza, 2020 #44;Gong, 2020 #45}, 但是從產(chǎn)量構(gòu)成因子分析間作優(yōu)勢的響應機制還鮮見報道?；诖? 本試驗以黃土塬區(qū)特殊自然環(huán)境為背景, 以玉米和大豆為研究對象, 探究種植密度和比例搭配對間作系統(tǒng)作物產(chǎn)量構(gòu)成的影響, 揭示間作優(yōu)勢的產(chǎn)量構(gòu)成因素響應機理, 旨在為有效調(diào)控間作系統(tǒng)作物產(chǎn)量優(yōu)勢和旱作農(nóng)業(yè)間作系統(tǒng)高產(chǎn)高效研究提供理論支撐。

1 材料與方法

1.1 試驗地概況

中國科學院水利部水土保持研究所長武黃土高原農(nóng)業(yè)生態(tài)試驗站位于黃土高原中南部陜甘交界處的陜西省長武縣洪家鎮(zhèn)王東村(107°40′30″E, 35°12′30″N), 海拔1200 m, 屬暖溫帶半濕潤半干旱大陸性季風氣候區(qū), 光照充足, 晝夜溫差較大; 年均降雨量584.1 mm, 且多集中于7月至9月份, 年均氣溫9.1℃, 無霜期171 d, 作物種植多為一年一季; 地下水埋深50~80 m, 田間持水量年均值在20%左右, 屬典型旱作農(nóng)業(yè)區(qū); 地貌屬高原溝壑區(qū), 塬面和溝壑兩大地貌單元各占35%和65%; 地帶性土壤為黑壚土, 土體結(jié)構(gòu)均勻疏松, 是黃土高原溝壑區(qū)典型性土壤。試驗作物生育期內(nèi)的降雨量和平均溫度見圖1。

1.2 試驗方法

1.2.1 試驗設計 試驗于2012—2013年4月至10月進行。玉米品種為鄭單958 (M), 種植密度為9株 m-2; 大豆品種為中黃24, 設9 (低密度S1)、21 (中密度S2)、33 (高密度S3)株 m-23個密度梯度; 玉米與3個密度梯度大豆進行間作(MS1、MS2、MS3)種植, 種植比例為2∶0、0∶2、2∶2和2∶4, 共10個處理, 3次重復, 30個小區(qū), 50 cm等行距覆膜種植。小區(qū)面積為6 m × 4 m, 隨機區(qū)組排列。播前施用底肥 N 90 kg hm-2、P2O5150 kg hm-2, 玉米分別于大喇叭口期和吐絲期各追施氮肥67.5 kg hm-2, 其他農(nóng)田管理與當?shù)叵嗤?/p>

1.2.2 樣品的測定 分別在玉米、大豆成熟期, 在每個小區(qū)劃定6 m2人工測產(chǎn)區(qū), 測定玉米和大豆產(chǎn)量。玉米收獲時每個小區(qū)隨機選取代表性植株手工摘取10個果穗, 測定穗重、穗長、禿尖長、穗粗、行數(shù)、行粒數(shù)、穗粒重、軸粗、軸重和千粒重。大豆收獲時每個小區(qū)隨機選取代表性植株10株, 測定單株莢數(shù)、單株總粒數(shù)、單株有效粒數(shù)、單株粒重、節(jié)數(shù)、結(jié)莢長度、主莖節(jié)間長度和百粒重。

1.2.3 產(chǎn)量及競爭指數(shù)計算 土地當量比(land equivalent ratio)是評價間作系統(tǒng)具有產(chǎn)量優(yōu)勢的重要指標[22], 為2種或多種物種間作產(chǎn)量與相應單作產(chǎn)量比值的總和。具體計算如下:

式中,m、s分別表示玉米和大豆在單作中的產(chǎn)量,mi、si分別表示玉米和大豆在間作中的產(chǎn)量,m、s分別表示玉米和大豆在間作中的土地當量比,表示間作系統(tǒng)的土地當量比。＞1, 表明間作系統(tǒng)存在產(chǎn)量優(yōu)勢;＜1, 表明間作系統(tǒng)產(chǎn)量不會增加。

實際產(chǎn)量損失(actual yield loss)表示在間作系統(tǒng)下作物的相對產(chǎn)量損失或增加[23]。具體計算如下:

式中,m、s分別表示單作中玉米和大豆的種植比例,mi、si分別表示間作中玉米和大豆的種植比例,m、s分別表示玉米和大豆在間作系統(tǒng)中的實際產(chǎn)量損失,表示間作系統(tǒng)的實際產(chǎn)量損失。＞0, 表明間作系統(tǒng)具有間作優(yōu)勢,＜0, 表明間作系統(tǒng)沒有產(chǎn)量優(yōu)勢。

侵占力(aggressivity)指間作系統(tǒng)中作物的相對產(chǎn)量增長大于另一作物產(chǎn)量增長的程度大小[24]。作物的侵占力具體計算如下:

式中,m、s分別表示間作系統(tǒng)中玉米和大豆的侵占力。m＝0, 表明這2種作物的競爭力相同;m＞0, 表明玉米的競爭力高于大豆;s＞0, 表明大豆的競爭力高于玉米。

競爭比率(competitive ratio)為間作系統(tǒng)中作物相對于另一作物的競爭能力[25]。具體計算如下:

式中,m和s分別表示玉米和大豆在間作系統(tǒng)中的競爭比率。m＞1, 表明玉米比大豆的競爭能力強;s＞1, 表明大豆比玉米的競爭能力強。

1.3 數(shù)據(jù)分析

采用Microsoft Excel 2007整理數(shù)據(jù), 采用SPSS 22.0統(tǒng)計軟件分析數(shù)據(jù)。單因素方差分析用于比較種植密度或種植比例下產(chǎn)量及其產(chǎn)量構(gòu)成、競爭指數(shù)的差異, Pearson相關系數(shù)用于反映所測指標之間的相關關系, 通徑分析用于反映產(chǎn)量構(gòu)成對于產(chǎn)量貢獻的直接作用和間接作用。利用SigmaPlot 12.5作圖。

2 結(jié)果與分析

2.1 不同間作方式對產(chǎn)量的影響

間作系統(tǒng)的土地當量比為1.14~1.23, 即單作需增加14%~23%的土地面積才能達到與間作相同的產(chǎn)量, 表現(xiàn)出明顯的間作優(yōu)勢(圖2)。同一種植比例下, 大豆密度對玉米和大豆的土地當量比無顯著影響(> 0.05)。同一密度梯度下, 2∶2間作下玉米的土地當量比顯著高于2∶4間作(< 0.05); 2∶2間作下大豆的土地當量比顯著低于2∶4間作(< 0.05)。

間作系統(tǒng)中玉米的實際產(chǎn)量損失大于0, 大豆的實際產(chǎn)量損失小于0 (圖3-a)。種植比例和大豆密度對玉米的實際產(chǎn)量損失無影響; 同一種植比例下, 大豆密度對大豆的實際產(chǎn)量損失無影響, 大豆中密度時, 2∶4間作下大豆的實際產(chǎn)量損失顯著高于2∶2間作(< 0.05), 大豆低高密度時, 間作比例對大豆的實際產(chǎn)量損失無顯著影響(圖3-a)。對間作系統(tǒng)整體而言, 種植比例和大豆密度組合下的間作系統(tǒng)的實際產(chǎn)量損失大于0 (圖3-b)。同一密度下, 種植比例對間作系統(tǒng)的實際產(chǎn)量損失無顯著影響。同一種植比例下, 大豆密度對間作系統(tǒng)的實際產(chǎn)量損失無顯著影響(圖3-b)。

MS1、MS2、MS3分別表示玉米與低密度(9株 m-2)、中密度(21株 m-2)、高密度(33株 m-2)的大豆間作。白色、灰色符號分別表示玉米與大豆以2:2、2:4種植比例的間作方式。

MS1, MS2, and MS3 indicate maize intercropping soybean with low density (9 plants m–2), medium density (21 plants m–2), and high density (33 plants m–2), respectively. White and gray symbols indicate 2:2 and 2:4 row proportion in maize-soybean intercropping systems.

2.2 不同間作方式對玉米和大豆產(chǎn)量構(gòu)成的影響

2.2.1 不同間作方式對玉米產(chǎn)量構(gòu)成的影響 種植比例和大豆密度間作組合下的穗重、穗長、禿尖長、穗粗、行數(shù)、行粒數(shù)、穗粒重、軸粗、軸重和千粒重高于單作或與單作間無顯著差異(表1)。玉米與低密度大豆以2∶4間作的穗重、穗長、穗粗、行數(shù)、行粒數(shù)、穗粒重、軸粗、軸重和千粒重均顯著高于單作。同一種植比例下(2∶2或2∶4), 大豆密度對間作玉米的穗重、禿尖長、穗粗、行數(shù)、穗粒重和千粒重均無顯著影響。大豆低密度時, 種植比例對穗重、穗長、禿尖長、穗粗、行數(shù)、軸粗和千粒重均無顯著影響, 2∶4種植比例下玉米的行粒數(shù)和軸重均顯著高于2∶2。大豆中密度時, 種植比例對穗重、穗長、禿尖長、穗粗、行數(shù)、行粒數(shù)、穗粒重、軸重和千粒重均無顯著影響, 2∶4種植比例下玉米的軸粗顯著高于2∶2。大豆高密度時, 種植比例對穗重、穗長、禿尖長、穗粗、行數(shù)、行粒數(shù)、穗粒重、軸粗、軸重和千粒重均無顯著影響。

2.2.2 不同間作方式對大豆產(chǎn)量構(gòu)成的影響 除2∶2種植比例下低密度的結(jié)莢長度和主莖節(jié)間長度顯著高于單作外, 其他種植方式下間作大豆的單株粒重、單株莢數(shù)、單株總粒數(shù)、單株有效粒數(shù)、結(jié)莢長度、主莖節(jié)間長度、主莖節(jié)數(shù)和百粒重低于單作或與單作間無顯著差異(表2)。大豆的單株粒重、單株莢數(shù)、單株總粒數(shù)和有效粒數(shù)隨密度的增加呈下降的趨勢。大豆的主莖節(jié)間長度隨密度的增加呈增加的趨勢。大豆密度對結(jié)莢長度和百粒重均無顯著影響。

處理同圖2。白色、灰色符號分別表示玉米與大豆以2:2、2:4種植比例的間作方式。

Treatments are the same as those given in Fig. 2. White and gray symbols indicate 2:2 and 2:4 row proportion in maize-soybean intercropping systems.

表1 不同間作方式對玉米產(chǎn)量構(gòu)成的影響

M表示單作玉米。MS1、MS2、MS3分別表示玉米與低密度(9株 m-2)、中密度(21株 m-2)、高密度(33株 m-2)的大豆間作。2:2和2:4分別表示玉米與大豆以2:2、2:4種植比例的間作方式。同列數(shù)據(jù)不同小寫字母表示同一種植比例下不同密度處理間差異達0.05顯著水平。

M indicates sole-cropped maize. MS1, MS2, and MS3 indicate maize intercropped soybean with low density (9 plants m–2), medium density (21 plants m–2), and high density (33 plants m–2), respectively. 2:2 and 2:4 indicate maize intercropped soybean at 2:2, 2:4 row proportion, respectively. Values within a column followed by different lowercase letters mean significant differences among different density treatments in the same row proportions at the 0.05 probability level.

表2 不同間作方式對大豆產(chǎn)量構(gòu)成的影響

S1、S2、S3分別表示低密度(9株 m-2)、中密度(21株 m-2)、高密度(33株 m-2)下的單作大豆, 其他處理同表1。同列數(shù)據(jù)不同小寫字母表示同一種植比例下不同密度處理間的差異達0.05顯著水平。同列數(shù)據(jù)不同大寫字母表示同一密度梯度下不同種植比例間的差異達0.05顯著水平。

S1, S2, and S3 indicate sole-cropped soybean with low density (9 plants m–2), medium density (21 plants m–2), and high density (33 plants m–2), respectively. Other treatments are the same as those given in Table 2. For each column, different letters denote significant difference at< 0.05 among different soybean density in the same row proportion (lowercase letters), and among different row proportions in the same soybean density (capital letters).

2.3 不同間作方式對作物競爭指數(shù)的影響

間作系統(tǒng)中玉米的侵占力大于0, 大豆的侵占力小于0, 種植比例和大豆種植密度對玉米和大豆的侵占力無顯著影響, 2∶2和2∶4種植比例下玉米和大豆的侵占力分別為0.73、0.71和-0.73、-0.71 (圖4-a, b)。間作系統(tǒng)中玉米的競爭比率大于1, 大豆的競爭比率小于1, 與低中高密度大豆間作, 玉米和大豆的競爭比率分別為2.08、1.84、1.68和0.49、0.56、0.63, 呈現(xiàn)隨大豆密度增加玉米的競爭比率下降而大豆的競爭比率增加的趨勢(圖4-c, d)。同一密度梯度下, 種植比例對玉米和大豆的競爭比率無顯著影響; 2∶4種植比例下高密度玉米的競爭比率顯著高于低密度, 2∶2種植比例下, 大豆密度對玉米和大豆的競爭比率無顯著影響。

處理同圖2。不同小寫字母表示同一種植比例下不同密度梯度間的差異顯著, 不同大寫字母表示同一密度梯度下不同種植比例之間的差異顯著(< 0.05)。

Treatments are the same as those given in Fig. 2. Different letters denote significant difference at< 0.05 among different soybean density in the same row proportion (lowercase letters), and among different row proportions in the same soybean density (capital letters).

2.4 不同間作方式對玉米產(chǎn)量的貢獻

由玉米產(chǎn)量與產(chǎn)量構(gòu)成之間的相關分析可以看出, 玉米產(chǎn)量與穗重、穗長、禿尖長、行粒數(shù)、穗粒重、軸重、千粒重均呈顯著相關關系(表3)。具體表現(xiàn)為玉米產(chǎn)量與軸重、千粒重、穗重、穗粒重、穗長、行粒數(shù)呈極顯著正相關關系, 相關系數(shù)分別為0.89、0.82、0.80、0.73、0.71和0.61, 與禿尖長呈極顯著負相關關系, 相關系數(shù)為-0.48 (<0.01)。軸重與穗重、穗長、禿尖長、行數(shù)、行粒數(shù)、穗粒重、軸粗呈顯著相關關系。千粒重與穗重、穗長、行粒數(shù)、穗粒重、軸重呈極顯著正相關關系, 與禿尖長呈負相關關系(<0.01)。穗重與穗長、穗粗、行數(shù)、行粒數(shù)、穗粒重、軸重、千粒重呈極顯著正相關關系, 與禿尖長呈極顯著負相關關系(<0.01)。穗粒重與穗重、穗長、行數(shù)、行粒數(shù)、軸重、千粒重呈極顯著正相關關系, 與禿尖長呈負相關關系。

為進一步確定產(chǎn)量構(gòu)成對玉米產(chǎn)量的影響效應,對其進行通徑分析發(fā)現(xiàn), 穗粒重對玉米產(chǎn)量的直接作用最大(表4)。具體表現(xiàn)為產(chǎn)量構(gòu)成對玉米產(chǎn)量的直接作用(直接通徑系數(shù))大小依次為穗粒重(2.18)>穗重(–1.71)>禿尖長(0.94)>行粒數(shù)(0.51)>軸重(0.24)>穗長(-0.13)>千粒重(0.05)。間接作用中, 軸重通過穗粒重對玉米產(chǎn)量的貢獻最大(1.64), 其次為千粒重, 通過穗粒重對玉米產(chǎn)量的貢獻(1.58)。

3 討論

3.1 不同間作方式對產(chǎn)量的影響

作物在空間和時間上對資源具有互補效應, 這一互補效應能夠提高產(chǎn)量和資源利用效率, 因此間作系統(tǒng)一直是我國重要的間作模式之一[26]。本研究中, 不同大豆密度和種植比例組合的間作方式的土地當量比為1.14~1.23, 表明玉米大豆間作系統(tǒng)具有產(chǎn)量優(yōu)勢。這與前人關于玉米和豌豆[6]、玉米和花生[27]、谷子和花生[28]的研究結(jié)果一致。宋振偉等[29]的研究表明, 增加作物種植密度提高作物的群體產(chǎn)量以彌補單株產(chǎn)量的減少。作物產(chǎn)量的多少由干物質(zhì)積累決定, 而合理的種植密度是提高產(chǎn)量的重要增產(chǎn)措施之一[30]。本研究大豆密度對間作群體的產(chǎn)量無顯著影響, 這可能是由于間作系統(tǒng)中大豆是競爭劣勢作物[31]。這與前人[32]的研究結(jié)果一致, 表明大豆密度對間作群體產(chǎn)量的影響程度較小。

表3 玉米產(chǎn)量與產(chǎn)量構(gòu)成之間的相關性

*、**分別表示在0.05和0.01水平相關性顯著。

*and**mean significant correlation at the 0.05 and 0.01 probability levels, respectively.

表4 間作玉米產(chǎn)量決定因子的通徑分析

間作系統(tǒng)中玉米的競爭比率大于1, 玉米的侵占力大于0, 表明在玉米和大豆共同生長期, 間作系統(tǒng)中大豆處于競爭劣勢, 玉米是競爭優(yōu)勢作物。前人研究認為, 間作系統(tǒng)中不同作物在生長過程中存在強烈的種間競爭關系, 并且玉米的資源競爭力高于大豆[33]、豌豆[12]。玉米的實際產(chǎn)量損失大于0, 表明間作系統(tǒng)中玉米具有產(chǎn)量優(yōu)勢。大豆的實際產(chǎn)量損失小于0, 表明間作系統(tǒng)中大豆沒有產(chǎn)量優(yōu)勢, 這與間作系統(tǒng)中玉米是競爭優(yōu)勢作物, 大豆是競爭劣勢作物相一致。種間競爭的結(jié)果是間作系統(tǒng)的實際產(chǎn)量損失大于0, 表明玉米和大豆間作具有產(chǎn)量優(yōu)勢的原因是玉米產(chǎn)量增加, 大豆產(chǎn)量不變或降低。這與前人關于小麥和豌豆[8]、燕麥和野豌豆[25]、谷子和大豆[34]間作的研究一致, 表明豆禾間作系統(tǒng)中, 禾本科作物的產(chǎn)量增加, 豆科作物的產(chǎn)量減少, 主要由間作系統(tǒng)中作物的競爭能力差異決定[8]。另外, 高位玉米與低位大豆間作時, 玉米處于生態(tài)位優(yōu)勢, 大豆處于生態(tài)位劣勢, 高稈玉米對大豆的遮陰作用降低了作物的光合作用[34], 并隨著生育期過程逐步加強, 進一步影響大豆的產(chǎn)量構(gòu)成乃至籽粒產(chǎn)量。玉米上部冠層葉片降低玉米中間層葉片和大豆冠層的透光率, 從而影響作物對養(yǎng)分的吸收, 間作系統(tǒng)中通過適當去掉玉米上部葉片提高作物對氮、磷、鉀等養(yǎng)分的吸收和作物產(chǎn)量[35]。

3.2 不同間作方式對產(chǎn)量構(gòu)成的影響

間作玉米的穗重、穗長、穗粗、穗粒重、軸重和千粒重等產(chǎn)量構(gòu)成均顯著高于單作; 除結(jié)莢長度和主莖節(jié)間長度外, 間作方式下大豆的單株莢數(shù)、單株總粒數(shù)、單株有效粒數(shù)、主莖節(jié)數(shù)、百粒重等產(chǎn)量構(gòu)成低于單作或與單作間無顯著差異。表明, 與單作相比, 間作方式下玉米的產(chǎn)量構(gòu)成提高, 大豆的產(chǎn)量構(gòu)成降低或沒有變化, 這與呂越等[36]的研究結(jié)果一致。這也進一步證明玉米大豆具有間作優(yōu)勢的主要原因是間作方式提高了玉米產(chǎn)量的構(gòu)成, 從而提高了玉米產(chǎn)量。不同間作方式對產(chǎn)量構(gòu)成的影響不同, 比如不同密度玉米與小麥間作的研究發(fā)現(xiàn), 低密度間作玉米的穗數(shù)顯著低于高密度20%, 千粒重高于高密度5.0%[37]。本研究結(jié)果表明, 玉米與低密度大豆以2∶4間作的穗重、穗長、穗粗等產(chǎn)量構(gòu)成均顯著高于單作。間作方式對作物產(chǎn)量構(gòu)成的影響還受到覆膜、水分等條件的影響, 比如, 間作玉米的穗粒重、百粒重隨灌水量增加呈先升后降的趨勢[38]。覆膜間作中豌豆的單株粒數(shù)和單株粒重較不覆膜間作平均增加43.7%和60.3%, 覆膜間作玉米的百粒重較不覆膜間作增加10.3%[12]。

3.3 不同間作方式對玉米產(chǎn)量的貢獻

玉米產(chǎn)量是由多基因控制的數(shù)量性狀, 是各農(nóng)藝性狀綜合作用的結(jié)果。提高玉米產(chǎn)量需要綜合考慮穗粒重、穗重、穗長等農(nóng)藝性狀, 這些性狀之間不僅相互作用而且相互影響, 對玉米產(chǎn)量的貢獻大小也不同[39]。玉米產(chǎn)量與穗重、穗長、穗粗、軸粗、行粒數(shù)、穗粒重等穗部相關性狀的產(chǎn)量構(gòu)成因素均存在顯著相關性[40], 故產(chǎn)量構(gòu)成在一定程度上影響了玉米的生長和產(chǎn)量的提高。徐磊等[41]關于黑龍江西部干旱區(qū)玉米產(chǎn)量與產(chǎn)量構(gòu)成的研究發(fā)現(xiàn), 玉米產(chǎn)量與穗粒重呈極顯著正相關。安治良[42]以自主選育的13個玉米雜交的研究發(fā)現(xiàn), 夏玉米產(chǎn)量與穗粒重、出籽率呈極顯著正相關。因此, 不同研究材料、栽培措施和生態(tài)環(huán)境可能影響作物產(chǎn)量與產(chǎn)量構(gòu)成因子之間的相關性。對于黃土塬區(qū)的間作玉米, 本研究發(fā)現(xiàn)間作玉米產(chǎn)量與軸重、千粒重、穗重、穗粒重、穗長、行粒數(shù)均呈顯著正相關關系。

產(chǎn)量構(gòu)成因子既通過本身的直接效應也可以通過其他因子的間接效應實現(xiàn)對玉米產(chǎn)量的貢獻, 因此玉米產(chǎn)量的提高是各個產(chǎn)量構(gòu)成因子綜合作用的結(jié)果, 并且因子之間存在相互作用。安治良[42]發(fā)現(xiàn), 穗粒重、出籽率和倒伏率對玉米產(chǎn)量的直接和間接作用最大。玉米的產(chǎn)量構(gòu)成與單株產(chǎn)量通徑分析發(fā)現(xiàn), 對產(chǎn)量貢獻最大的是千粒重[43]。傳統(tǒng)玉米品種對產(chǎn)量貢獻作用最大的是千粒重, 而新型玉米品種對產(chǎn)量作用貢獻最大的是行粒數(shù)[18]。穗粒重與穗重、穗長、穗粗、軸重、千粒重均呈顯著正相關[44], 并且穗粗通過穗粒重對玉米產(chǎn)量的貢獻最大[41]。本研究對間作玉米產(chǎn)量與產(chǎn)量構(gòu)成之間進行通徑分析發(fā)現(xiàn), 直接作用表現(xiàn)為穗粒重是提高間作玉米產(chǎn)量的主要決定因子, 間接作用中軸重、千粒重通過穗粒重對玉米產(chǎn)量的貢獻最大。因此, 探究影響玉米甚至間作系統(tǒng)的產(chǎn)量構(gòu)成因素對玉米產(chǎn)量的貢獻大小,需要綜合考慮研究材料、栽培方式和環(huán)境等因素, 從而為提高當?shù)刈魑锂a(chǎn)量提供理論依據(jù)。

4 結(jié)論

通過不同種植比例和大豆種植密度對玉米和大豆間作產(chǎn)量的影響分析發(fā)現(xiàn), 種植比例和大豆種植密度對間作產(chǎn)量及產(chǎn)量構(gòu)成具有重要影響。玉米和大豆間作具有產(chǎn)量優(yōu)勢, 玉米與高密度大豆間作的土地當量比最高(2∶2和2∶4種植比例分別為1.22和1.23), 并且提高間作群體產(chǎn)量的原因是玉米產(chǎn)量的提高以及玉米對資源的競爭能力高于大豆(侵占力和競爭比率分別為0.72和1.86)。通過通徑分析發(fā)現(xiàn), 穗粒重對玉米產(chǎn)量的貢獻作用最大(直接作用為2.18)。因此, 玉米與高密度大豆間作的土地當量比高于其他間作方式, 能夠最大化的提高作物的產(chǎn)量, 應該適當進行推廣。

[1] Singh P, Jadhav A, Varshney M. Light interception and light use efficiency in sorghum based intercropping system., 2002, 4: 93–96.

[2] Fustec J, Lesuffleur F, Mahieu S, Cliquet J B. Nitrogen rhizodeposition of legumes: a review., 2010, 30: 57–66.

[3] Zhu Y, Chen H, Fan J, Wang Y, Yan L, Chen J, Fan J X, Yang S, Hu L, Leung H. Genetic diversity and disease control in rice., 2000, 406: 718–722.

[4] Banik P, Midya A, Sarkar B K, Ghose S S. Wheat and chickpea intercropping systems in an additive series experiment: advantages and weed smothering., 2006, 24: 325–332.

[5] Li R, Zhang Z, Tang W, Huang Y, Coulter J A, Nan Z. Common vetch cultivars improve yield of oat row intercropping on the Qinghai-Tibetan plateau by optimizing photosynthetic performance., 2020, 117: 1–13.

[6] Mao L, Zhang L, Li W, van der Werf W, Sun J, Spiertz H, Li L. Yield advantage and water saving in maize/pea intercrop., 2012, 138: 11–20.

[7] Chen P, Du Q, Liu X, Li Z, Sajad H, Lu L, Song C, Wang X, Liu W, Feng Y. Effects of reduced nitrogen inputs on crop yield and nitrogen use efficiency in a long-term maize–soybean relay strip intercropping system., 2017, 12: 1–19.

[8] Lithourgidis A S, Vlachostergios D N, Dordas C A, Damalas C A. Dry matter yield, nitrogen content, and competition in pea-cereal intercropping systems., 2011. 34: 287–294.

[9] 張向前. 不同氮水平下玉米間作大豆和花生的效應研究. 南京農(nóng)業(yè)大學博士學位論文, 江蘇南京, 2013.

Zhang X Q. A Study on the Effect of Maize Intercropped with Soybean and Peanut under Different Nitrogen Fertilization Levels. PhD Dissertation of Nanjing Agricultural University, Nanjing, Jiangsu, China, 2013 (in Chinese with English abstract).

[10] Zhang F, Li L. Using competitive and facilitative interactions in intercropping systems enhances crop productivity and nutrient-use efficiency, 2003, 248: 305–312.

[11] 趙建華, 孫建好, 陳亮之, 李偉綺. 玉米行距對大豆/玉米間作作物生長及種間競爭力的影響大豆科學, 2019, 38: 229–235.

Zhao J H, Sun J H, Chen L Z, Li W Q. Growth and interspecific competition of crops as affected by maize row spacing in soybean/maize intercropping system., 2019, 38: 229–235 (in Chinese with English abstract).

[12] 趙建華, 孫建好, 李偉綺. 覆膜對玉米‖豌豆作物生產(chǎn)力及種間互作的影響.干旱地區(qū)農(nóng)業(yè)研究, 2020, 38(2): 164–169.

Zhao J H, Sun J H, Li W Q. Effect of film mulching on productivity and interspecific interaction in maize–pea intercropping system., 2020, 38(2): 164–169 (in Chinese with English abstract).

[13] Undie U, Uwah D, Attoe E. Effect of intercropping and crop arrangement on yield and productivity of late season maize/soybean mixtures in the humid environment of South Southern Nigeria., 2012, 4: 37–50.

[14] Wang R, Sun Z, Zhang L, Yang N, Feng L, Bai W, Zhang D, Wang Q, Evers J B, Liu Y. Border-row proportion determines strength of interspecific interactions and crop yields in maize/peanut strip intercropping., 2020, 253: 107819.

[15] 杜進勇, 柴強, 王一帆, 范虹, 胡發(fā)龍, 殷文, 李登業(yè). 地上地下互作強度對小麥間作玉米光合特性的影響.作物學報, 2019, 45: 1398–1406.

Du J Y, Chai Q, Wang Y F, Fan H, Hu F L, Yin W, Li D Y. Effect of above- and below-ground interaction intensity on photosynthetic characteristics of wheat–maize intercropping., 2019, 45: 1398–1406 (in Chinese with English abstract).

[16] 任旭靈, 滕園園, 王一帆, 殷文, 柴強. 玉米間作豌豆種間競爭互補對少耕密植的響應.中國生態(tài)農(nóng)業(yè)學報, 2019, 27: 860–869.

Ren X L, Teng Y Y, Wang Y F, Yin W, Chai Q. Response of interspecific competition and complementarity of maize/pea intercropping to reduced tillage and high-density planting., 2019, 27: 860–869 (in Chinese with English abstract).

[17] 于海秋, 徐克章, 陳學求, 武志海, 姜子漣, 沈秀瑛. 玉米主要抗旱性狀的配合力及遺傳參數(shù)分析: I．產(chǎn)量性狀玉米科學, 2003, 11(1): 12–18.

Yu H Q, Xu K Z, Chen X Q, Wu Z H, Jiang Z L, Shen X Y. Analysis of combining ability and hereditary parameter of major drought resistance traits in maize: I. yield traits., 2003, 11(1): 12–18 (in Chinese with English abstract).

[18] 魯珊, 肖荷霞, 毛彩云, 陸建章, 岳金生. 不同類型玉米主要農(nóng)藝性狀的相關和通徑分析.中國農(nóng)學通報, 2019, 35(19): 11–14. .

Lu S, Xiao H X, Mao C Y, Lu J Z, Yue J S. Correlation and path analysis of major agronomic traits of maize types., 2019, 35(19): 11–14 (in Chinese with English abstract).

[19] Raza M A, Feng L Y, Werf W V D, Iqbal N. Optimum strip width increases dry matter, nutrient accumulation, and seed yield of intercrops under the relay intercropping system., 2020, 1: 1–14.

[20] Gong X, Dang K, Lyu S, Zhao G, Tian L, Luo Y, Feng B. Interspecific root interactions and water use efficiency of intercropped proso millet and mung bean., 2020, 115: 1–12.

[21] Raza M A, Khalid M H B, Zhang X, Feng L Y, Khan I, Hassan M J, Ahmed M, Ansar M, Chen Y K, Fan Y F. Effect of planting patterns on yield, nutrient accumulation and distribution in maize and soybean under relay intercropping systems., 2019, 9: 4947.

[22] Mead R, Willey R W. The concept of a ‘land equivalent ratio’ and advantages in yields from intercropping., 1980, 16: 217–228.

[23] Banik P, Sasmal T, Ghosal P K, Bagchi D K. Evaluation of mustard (var. Toria) and legume intercropping under 1:1 and 2:1 row-replacement series systems., 2000, 185: 9–14.

[24] Agegnehu G, Ghizaw A, Sinebo W. Yield performance and land-use efficiency of barley and faba bean mixed cropping in Ethiopian highlands, 2006, 25: 202–207.

[25] Dhima K V, Lithourgidis A S, Vasilakoglou I B, Dordas C A. Competition indices of common vetch and cereal intercrops in two seeding ratio., 2007, 100: 249–256.

[26] 張建華, 馬義勇, 王振南, 齊晶. 間作系統(tǒng)中玉米光合作用指標改善的研究.玉米科學, 2006, 14(4): 104–106.

Zhang J H, Ma Y Y, Wang Z N, Qi J. Research on the improvement of photosynthesis indices of maize in the intercropping system., 2006, 14(4): 104–106 (in Chinese with English abstract).

[27] Ghosh P K. Growth, yield, competition and economics of groundnut/cereal fodder intercropping systems in the semi-arid tropics of India., 2004, 88: 227–237.

[28] Feng L, Sun Z, Zheng M, Muchoki M, Zheng J, Yang N, Bai W, Feng C, Zhang Z, Cai Q, Zhang D. Productivity enhancement and water use efficiency of peanut-millet intercropping., 2016, 48: 1459–1466.

[29] 宋振偉, 齊華, 張振平, 錢春榮, 郭金瑞, 鄧艾興, 張衛(wèi)建. 春玉米中單909農(nóng)藝性狀和產(chǎn)量對密植的響應及其在東北不同區(qū)域的差異.作物學報, 2012, 38: 2267–2277.

Song Z W, Qi H, Zhang Z P, Qian C R, Guo J R, Deng A X, Zhang W J. Effects of plant density on agronomic traits and yield in spring maize Zhongdan 909 and their regional differences in Northeast China., 2012, 38: 2267–2277 (in Chinese with English abstract).

[30] 康彩睿, 謝軍紅, 李玲玲, 王嘉男, 郭喜軍, 彭正凱, 王進斌, Fudjoe S, 王林林. 種植密度與施氮量對隴中旱農(nóng)區(qū)玉米產(chǎn)量及光合特性的影響.草業(yè)學報, 2020, 29(5): 141–149.

Kang C R, Xie J H, Li L L, Wang J N, Guo X J, Peng Z K, Wang J B, Fudjoe S, Wang L L. Effects of planting density and nitrogen fertilizer rate on maize yield and photosynthetic characteristics in arid areas of central Gansu, China., 2020, 29(5): 141–149 (in Chinese with English abstract).

[31] 任媛媛, 王志梁, 王小林, 張歲岐. 黃土塬區(qū)玉米大豆不同間作方式對產(chǎn)量和經(jīng)濟收益的影響及其機制.生態(tài)學報, 2015, 12: 4168–4177.

Ren Y Y, Wang Z L, Wang X L, Zhang S Q. The effect and mechanism of intercropping pattern on yield and economic benefit on the Loess Plateau., 2015, 12: 4168–4177 (in Chinese with English abstract).

[32] 朱元剛, 高鳳菊, 曹鵬鵬, 王樂政. 種植密度對玉米–大豆間作群體產(chǎn)量和經(jīng)濟產(chǎn)值的影響.應用生態(tài)學報, 2015, 26: 1751–1758.

Zhu Y G, Gao F J, Cao P P, Wang L Z. Effect of plant density on population yield and economic output value in maize–soybean intercropping., 2015, 26: 1751–1758 (in Chinese with English abstract).

[33] 呂越, 吳普特, 陳小莉, 王玉寶, 趙西寧. 地上部與地下部作用對玉米/大豆間作優(yōu)勢的影響.農(nóng)業(yè)機械學報, 2014, 45(1): 129–136.

Lyu Y, Wu P T, Chen X L, Wang Y B, Zhao X N. Effect of above- and below-ground interaction on maize/soybean intercropping advantage., 2014, 45(1): 129–136 (in Chinese with English abstract).

[34] 李智, 王宏富, 王鈺云, 楊凈, 魚冰星, 黃珊珊. 谷子大豆間作對作物光合特性及產(chǎn)量的影響.中國農(nóng)業(yè)科技導報, 2020, 22(6): 168–175.

Li Z, Wang H F, Wang Y Y, Yang J, Yu B X, Huang S S. Impact of millet and soybean intercropping on their photosynthetic characteristics and yield., 2020, 22(6): 168–175 (in Chinese with English abstract).

[35] Raza M A, Feng L Y, van der Werf W, Iqbal N, Khan, I., Hassan M J, Ansar M, Chen Y K, Xi Z J, Shi J Y, Ahmed M, Yang F, Yang W. Optimum leaf defoliation: a new agronomic approach for increasing nutrient uptake and land equivalent ratio of maize soybean relay intercropping system., 2019, 244: 1–11.

[36] 呂越, 吳普特, 陳小莉, 王玉寶, 趙西寧. 玉米/大豆間作系統(tǒng)的作物資源競爭.應用生態(tài)學報, 2014, 25: 139–146.

Lyu Y, Wu P T, Chen X L, Wang Y B, Zhao X N. Resource competition in maize/soybean intercropping system., 2014, 25: 139–146 (in Chinese with English abstract).

[37] 王一帆, 秦亞洲, 馮福學, 趙財, 于愛忠, 劉暢, 柴強. 根間作用與密度協(xié)同作用對小麥間作玉米產(chǎn)量及產(chǎn)量構(gòu)成的影響.作物學報, 2017, 43: 754–762.

Wang Y F, Qin Y Z, Feng F X, Zhao C, Yu A Z, Liu C, Chai Q. Synergistic effect of root interaction and density on yield and yield components of wheat/maize intercropping system., 2017, 43: 754–762 (in Chinese with English abstract).

[38] 張作為, 史海濱, 李禎, 李仙岳, 閆建文, 李介均. 不同生育時期非充分灌溉對間作作物產(chǎn)量構(gòu)成因子及收獲指數(shù)的影響干旱地區(qū)農(nóng)業(yè)研究, 2016, 34(4): 31–37.

Zhang Z W, Shi H B, Li Z, Li X Y, Yan J W, Li J J. The influence of deficit irrigation during different growth periods on the yield components and harvest index of intercropped crops., 2016, 34(4): 31–37 (in Chinese with English abstract).

[39] 楊坤, 段正鳳, 李福星, 陳芝能. 玉米穗部性狀及產(chǎn)量間的通徑分析.農(nóng)技服務, 2016, 33(1): 69–71.

Yang K, Duan Z F, Li F X, Chen Z N. Path analysis between ear traits and yield of maize., 2016, 33(1): 69–71 (in Chinese with English abstract).

[40] 湯華, 黃益勤, 嚴建兵, 劉宗華, 湯繼華, 鄭用璉, 李建生. 玉米優(yōu)良雜交種豫玉22產(chǎn)量性狀的遺傳分析.作物學報, 2004, 30: 922–926.

Tang H, Huang Y Q, Yan J B, Liu Z H, Tang J H, Zheng Y L, Li J S. Genetic analysis of yield traits with elite maize hybrid—Yuyu 22., 2004, 30: 922–926 (in Chinese with English abstract).

[41] 徐磊, 譚福忠, 師臣, 周長軍, 鄭巍, 齊國超, 陳剛, 于海峰. 黑龍江省西部干旱區(qū)玉米產(chǎn)量與產(chǎn)量構(gòu)成因素的相關性分析. 黑龍江農(nóng)業(yè)科學, 2020, (7): 1–6.

Xu L, Tan F Z, Shi C, Zhou C J, Zheng W, Qi G C, Chen G, Yu H F. Correlation analysis of maize yield and yield components in arid area of western Heilongjiang province., 2020, (7): 1–6 (in Chinese with English abstract).

[42] 安治良. 夏玉米雜交種農(nóng)藝性狀與產(chǎn)量的相關與通徑分析安徽農(nóng)學通報, 2018, 24(17): 34–35.

An Z L. Genetic correlation and path analysis of yield and agronomic characteristics of summer maize hybrids., 2018, 24(17): 34–35 (in Chinese with English abstract).

[43] 梁曉玲, 阿布來提, 馮國俊, 李進, 李銘東, 李維鼎, 高慧慧. 玉米雜交種的產(chǎn)量比較及主要農(nóng)藝性狀的相關和通徑分析.玉米科學, 2001, 9(1): 16–20.

Liang X L, ABU L T, Feng G J, Li J, Li M D, Li W D, Gao H H. Yield performance of maize hybrids and analysis of correlation between yield and agronomic characteristics., 2001, 9(1): 16–20 (in Chinese with English abstract).

[44] 譚巍巍, 石云素, 宋燕春, 楊德光, 王天宇, 黎裕, 王陽, 李永祥, 劉成, 劉志齋. 不同環(huán)境下多個玉米穗部性狀的QTL分析.中國農(nóng)業(yè)科學, 2011, 44: 233–244.

Tan W W, Shi Y S, Song Y Q, Yang D G, Wang T Y, Li Y, Wang Y, Li Y X, Liu C, Liu Z Z. QTL analysis of ear traits in maize across multiple environments., 2011, 44: 233–244 (in Chinese with English abstract).

Competitive effect of soybean density on yield formation in maize/soybean intercropping systems

REN Yuan-Yuan1, ZHANG Li2, YU Yao-Chuang1, ZHANG Yan-Jun1, and ZHANG Sui-Qi3,*

1College of Geography and Environment Engineering, Baoji University of Arts and Sciences / Shaanxi Key Laboratory of Disaster Monitoring and Mechanism Simulation, Baoji 721013, Shannxi, China;2School of Pharmacy, Weifang Medical University, Weifang 261053, Shandong, China;3State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau / Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, Shaanxi, China

Crop intercropping system in arid areas is an important planting method for optimizing crop population quality and improving crop yield. There are few reports on the mechanism of crop intercropping competitive advantage in terms of yield components. In this study, to explore the mechanism of crop competition for increasing the yield of intercropping systems, three soybean densities and four row proportions (maize and soybean intercropping with 2:0, 0:2, 2:2, and 2:4) were used to investigate the changes of competition index, yield components, and yield of intercropping system. The results showed that the yield of intercropping system with different row proportion and soybean density was increased by 14%–23%. The actual yield loss of maize was greater than 0, the actual yield loss of soybean was less than 0. The ear weight, ear length, ear diameter, grain weight per ear, cob weight, and 1000-grain weight of maize in intercropping systems were significantly higher than those in monoculture. Except pod length and internode length of main stem, the grain weight per plant, pod number per plants, grain number per plant, effective grain number per plants, nodes on main stem, and 100-grain weight of soybean in intercropping systems were lower than that of monoculture. And there was no significant difference between intercropping and monoculture. The competition ratio of maize was greater than one, and the competition ratio of soybean was less than one. The competition ratio of maize and soybean was 2.08, 1.84, 1.68 and 0.49, 0.56, 0.63, respectively, with three soybean densities, indicating that with the increase of soybean density competition ratio of maize increased while competition ratio of soybean decreased in intercrops. The aggressivity of maize was positive value, and that of soybean was less than zero in intercropping system. Maize yield was positively correlated with cob weight, 1000-grain weight, ear weight, grain weight per ear, ear length, and row kernel number negatively correlated with barren tip length. Path analysis revealed that in the direct effect, grain weight per ear contributed the most to maize yield (2.18); in indirect effect, cob weight and 1000-grain weight contributed more to maize yield 1000-grain weight per ear (1.64 and 1.58). In conclusion, maize intercropped with soybean had intercropping advantages that derived from grain weight per ear.

density; intercropping; path analysis; competition ratio

10.3724/SP.J.1006.2021.04226

本研究由陜西省科技廳自然科學基礎研究項目(2019JQ-895), 國家自然科學基金青年科學基金項目(41901025), 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室開放基金(A314021402-1712)和寶雞文理學院博士科研啟動費項目(ZK2017042)資助。

This study was supported by the Shaanxi National Science Foundation (2019JQ-895), the Youth Foundation of National Natural Science Foundation of China (41901025), the Open Foundation of State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau (A314021402-1712), and the PhD Research Startup Foundation of Baoji University of Arts and Sciences (ZK2017042).

張歲岐, E-mail: sqzhang@ms.iswc.ac.cn

E-mail: renyuanyuan0410@163.com

2020-10-12;

2021-01-13;

2021-02-20.

URL: https://kns.cnki.net/kcms/detail/11.1809.S.20210220.1144.002.html