蓄水坑灌坑深對土壤氧氣分布的影響

雷 震，郭向紅，雷 濤，路明杰，孫西歡,2，馬娟娟，張少文

(1.太原理工大學水利科學與工程學院，太原 030024；2.晉中學院，山西 晉中 030600)

0 引 言

良好的土壤氧氣狀況是保證果樹正常生長發育的必要條件。適宜的土壤氧氣條件能夠改善果樹根系呼吸作用[1]，促進果樹生長發育，進而實現果樹增產提質[2-4]。傳統農業灌溉模式按照濕潤土壤方式不同可分為溝灌、畦灌、漫灌、淹灌，水從地表面進人田間并借重力和毛細管作用浸潤土壤，是目前應用最為廣泛的灌溉模式，但灌溉水對空氣的驅替作用會引起土壤根區發生間歇性低氧脅迫問題，不利于果樹生長和產量形成[5]。在改善土壤氧氣狀況方面，國內外學者主要圍繞增氧灌溉開展了大量研究工作。已有研究表明，增氧灌溉能改善作物根區的生長環境[6]，促進根系生長[7]，提高作物的光合作用效率[8]，增加作物生長速率，提高作物產量及品質[9, 10]。這些研究工作主要集中在土壤氧氣狀況的定性研究方面，而定量方面研究報道較少[11]。通過建立土壤氧氣分布模型，能為準確描述土壤氧氣分布特征提供量化工具，但現有研究報道較少，有待進一步深入。蓄水坑灌法[12, 13]是一種有助于改善土壤氧氣狀況的新型中深層立體灌溉方法，蓄水坑壁能夠增加土壤與大氣間的接觸面積，進而改善土壤通氣條件，可以有效避免傳統灌溉方法的低氧脅迫問題。坑深[14]是蓄水坑灌法的重要技術參數之一，不同坑深會對土壤氧氣產生多大的影響還不是很清楚。開展不同蓄水坑深條件下土壤氧氣分布特征研究，對于完善蓄水坑灌技術具有重要的理論和實踐意義。近20 a，國內外學者主要圍繞蓄水坑灌條件下土壤水氮遷移分布特征[15]、果樹生理響應[16]、土壤微生物活性[17]等方面進行了大量研究工作，然而在土壤氧氣方面的研究報道較少。不同坑深條件下土壤氧氣分布特征還尚不明確，哪種坑深條件更有利于改善土壤通氣條件也有待進一步深入研究。

本文擬通過蓄水坑灌條件下田間試驗，探明不同坑深條件下土壤氧氣一維垂向、徑向及二維空間分布特征，構建蓄水坑灌土壤氧氣分布量化模型。研究成果可進一步豐富蓄水坑灌理論，為合理確定蓄水坑深提供理論依據，為合理調控果園根區土壤氣體環境提供參考。

1 材料與方法

1.1 試驗條件

試驗于2019年4-10月在位于山西省晉中市太谷縣的山西農科院果樹研究所蘋果園中進行。試驗區的地理坐標為東經112°32′，北緯37°23′，海拔781.9 m。所在區域屬于典型的大陸性半干旱氣候，多年平均氣溫9.8 ℃，無霜期175 d，平均降雨量460 mm。受季風影響，年內降雨嚴重不均，6-9月間降雨量大，占60%以上，12月至次年2月降雨很少，僅占2%左右。試驗區土壤為粉(沙)壤土，質地比較均一。土壤干容重為1.47 g/cm3，飽和含水率為49.21%，田間持水率為30%(體積含水率)，灌溉水源為地下水。試驗用果樹品種為“長富二號”，均選擇7 a生的矮砧蘋果樹為試材。

1.2 試驗設計

本試驗選取園區長勢優良，生長形態接近的6棵蘋果樹進行試驗研究。試驗以蓄水坑坑深為控制因子，設坑深40 cm(T1)和60 cm(T2)2個處理，重復3次，每棵樹為一個試驗小區。蘋果樹種植方向為南北向，株行距分別為2 m和4 m。采用蓄水坑灌法進行灌溉，布置方式為每株果樹周邊均勻布置4個直徑為30 cm的圓柱體蓄水坑，坑壁距樹干中心位置60 cm。坑底采用土工布墊底作隔水處理。各處理灌水上、下限均為90%和60%田間持水率。

1.3 測定指標及方法

土壤氧飽和度：采用13.05.01型土壤原位氧氣測定儀(上海賽弗生物科技有限公司)進行測定。在距離樹干中心30、50、100、120 cm徑向位置處布置土壤氧氣測點。在各個徑向位置處，沿垂向每20 cm一層布置測點，測定深度為120 cm，共6層。測點布置見圖1。

圖1 測點布置(單位：cm)Fig.1 The layout of measuring points

1.4 數據處理

選用數據為9月27日3個重復的平均值。采用Microsoft Office Excel 2016軟件進行原始數據處理，采用1stopt 1.5軟件進行數學建模，采用SPSS 25進行統計學分析，采用AutoCAD 2019、Origin 2019 軟件進行繪圖，圖中數據以平均值±標準差形式呈現。

2 結果與分析

2.1 不同坑深條件下土壤氧氣一維垂向分布特征

圖2為蓄水坑灌條件下距樹干不同徑向位置(r為30、50、100、120 cm)土壤氧氣一維垂向分布特征。由圖2可知，在距樹干同一徑向位置處，不同坑深處理后土壤氧飽和度隨土層深度增加呈現指數型下降趨勢。以40 cm坑深處理為例，進行土壤氧氣一維垂向分布特征分析。地表處土壤氧飽和度為1，則在各個徑向位置處，土層20、40、60、80、100、120 cm深度位置處土壤氧飽和度分別比地表處下降了1.9%～2.8%、3.1%～3.6%、3.9%～4.3%、4.4%～4.9%、4.8%～5.7%、5.2%～5.9%。這表明，0～40 cm深度范圍內土壤氧飽和度下降速率較40～120 cm范圍的快，并且40～120 cm深度范圍的土壤氧飽和度下降速率隨深度增加呈下降趨勢。這是由于土壤氧氣主要來源于大氣，距離地表越近時，大氣進入土體后的擴散路徑越短，土壤越容易與大氣進行氣體交換，土壤氧氣供應越充足，氧飽和度越高。由統計分析還可知，在r為30、120 cm徑向位置處，2個處理間垂向數據存在顯著性差異(P<0.05); 在r為50、100 cm徑向位置處，2個處理間垂向數據存在極顯著性差異(P<0.01)。坑深60 cm的處理，各層(20、40、60、80、100、120 cm)土壤氧飽和度分別比40 cm處理高0.10%～1.33%、0.30%～1.29%、0.80%～1.20%、1.00%～1.32%、1.00%～1.88%、1.20%～1.98%。由此說明，在任意徑向位置處，60 cm坑深處理下各層土壤氧飽和度均略微高于40 cm坑深處理。這是由于60 cm坑深處理條件下土壤臨空面積大于40 cm坑深處理，相對40 cm坑深處理，經60 cm坑深處理后的土壤坑壁面積(土壤-大氣交換界面)能夠增加33.3%，氣體交換更容易發生，土壤氧氣補給會更充足。

圖2 不同坑深條件下土壤氧氣一維垂向分布特征Fig.2 One-dimensional vertical distribution characteristics of soil oxygen under different pit depths

為了進一步描述蓄水坑灌不同坑深條件下土壤氧氣一維垂向分布特征，建立了土壤氧飽和度一維垂向分布模型OS(z)：

OS=aeb z

(1)

式中：OS為土壤氧飽和度；z為土層深度，cm；a、b為擬合參數。

模型參數見表1。由表1可知，不同坑深的土壤氧飽和度一維垂向分布模型OS(z)的決定系數R2為0.869 4～0.983 8，均值為0.923 9。說明該模型能夠較好地模擬蓄水坑灌不同坑深條件下土壤氧飽和度一維垂向分布特征。在OS(z)模型中，參數b能夠反映土壤氧飽和度沿垂向變化速率的快慢程度，其絕對值越大，土壤氧飽和度變化速率越快。由表1可知，在距離樹干30、50、100、120 cm特征位置處，40 cm坑深處理的參數b分別是坑深60 cm處理的2.70、1.12、1.12、2.49倍。這說明40 cm坑深處理的土壤氧飽和度沿垂向變化速率要高于60 cm坑深處理，這與圖2分析結果一致。

表1 不同坑深條件下土壤氧氣垂向分布擬合參數Tab.1 Fitting parameters of vertical distribution of soil oxygen under different pit depth conditions

2.2 不同坑深條件下土壤氧氣一維徑向分布特征

圖3為蓄水坑灌距地表不同土層深度處(20、40、60、80、100、120 cm)土壤氧氣一維徑向分布特征。以坑深60 cm處理為例，進行土壤氧氣一維徑向分布特征分析。在各土壤深度位置(20、40、60、80、100、120 cm)處，距離樹干50、100、120 cm位置處的土壤氧飽和度相對30 cm徑向位置處分別變化了0%～1.93%、0%～1.29%、0%～0.80%、0%～0.4%、0%～0.3%、0%～0.2%。蓄水坑壁處(徑向距離60、90 cm)土壤氧飽和度非常接近大氣。由此說明，在同一土層深度處，不同土壤氧飽和度隨徑向距離增加呈先增加后減小的趨勢，在坑壁處達到最大，但變化并不太明顯。這可能與果樹根系呼吸強度、土壤-大氣接觸面積及氣象條件有關。距離樹干越近處，根系分布越密集，根系呼吸強度越大，土壤氧氣含量越低。蓄水坑增加了土壤-大氣接觸面積，坑中氣體更容易向周圍土壤中擴散，坑壁周邊土壤氧氣能得到及時的補給和更新，使近坑壁位置處土壤氧飽和度較高。蓄水坑壁形成的臨空面與大氣直接接觸，受氣象條件的影響，坑壁周邊水分蒸散量較大，造成土壤孔隙度增加，附近土壤氧氣更易得到補給，因而近坑壁處土壤氧飽和度比距離坑壁遠處的高。由統計分析可知，在20、40 cm垂向位置處，2處理的數據差異不顯著；在60、80、100、120 cm垂向位置處，2處理間數據差異極顯著(P<0.01)。60 cm處理的各徑向位置(30、50、100、120 cm)土壤氧飽和度分別比40 cm處理大0.20%～1.94%、1.20%～1.60%、0.84%～1.20%、0.10%～1.80%。由此說明，在任意垂向位置處，60 cm坑深處理下各徑向位置土壤氧飽和度均高于40 cm坑深處理。這是因為60 cm坑深處理較40 cm坑深處理具有更大的氣體交換界面，土壤-大氣間氣體交換會更頻繁，擴散進入到土壤中氧氣量相對就更多。

圖3 不同坑深條件下土壤氧氣一維徑向分布特征Fig.3 One-dimensional radial distribution of soil oxygen under different pit depth conditions

為了進一步描述蓄水坑灌不同坑深條件下土壤氧氣一維徑向分布特征，以蓄水坑中心點為坐標原點，建立了土壤氧飽和度一維徑向分布模型OS(x)：

OS=ced x2

(2)

式中：OS為土壤氧飽和度；x為距離蓄水坑中心點的水平距離，cm；c、d為擬合參數。

模型參數見表2。由表2可知，模型OS(x)決定系數R2都在0.8以上，擬合精度較高。分析模型參數還可知，模型指數部分趨近于常數1，氧飽和度主要取決于參數c。由此表明，徑向距離對土壤氧飽和度的影響較小，這也與前文圖中反映的規律一致。

表2 不同坑深條件下土壤氧氣徑向分布擬合參數Tab.2 Fitting parameters of soil oxygen radial distribution under different pit depth conditions

2.3 不同坑深條件下土壤氧氣二維空間分布特征

圖4為蓄水坑灌不同坑深條件下土壤氧氣二維空間分布特征。由圖4可以看出，不同坑深處理下，土壤氧飽和度隨垂向深度增加呈現下降的趨勢，隨徑向距離增加呈現先增加后減小的趨勢。采用統計學方法對坑兩側數據進行分析可知，坑兩側氧飽和度相似，不存在顯著性差異。且由圖4可以看出，影響土壤氧飽和度的主要因素是土層深度，在較淺土層處土壤氧飽和度隨深度的變化幅度較大。不同坑深處理下，坑兩側的土壤氧飽和度基本上是以通過坑中心點的縱軸為基準對稱分布的。這表明，蓄水坑灌條件下，土壤氧飽和度主要與距蓄水坑坑壁的距離及蓄水坑深有關，而與距樹干中心點的距離關系不大。由圖4還可知，40 cm坑深處理土壤氧飽和度在0.98以上的區域主要集中在垂向0～60 cm和坑壁周邊20 cm范圍內；60 cm坑深處理土壤氧飽和度在0.98以上的區域主要集中在垂向0～80 cm和蓄水坑坑壁周邊20 cm范圍內。由此說明，在相同氧飽和度條件下，60 cm坑深處理相對40 cm坑深處理氧氣分布范圍更大。60 cm蓄水坑相比40 cm蓄水坑，更有助于提高大氣與土壤的氣體交換能力，提高土壤氧飽和度，改善果園根區的土壤含氧狀況。

圖4 不同坑深條件下土壤氧氣二維空間分布特征Fig.4 Two-dimensional spatial distribution of soil oxygen under different pit depth conditions

為了進一步對土壤氧飽和度二維分布特征進行定量描述，對式(1)和式(2)采用相乘形式進行耦合，構建了土壤氧飽和度二維空間分布模型OS(x,z)：

OS=ke(m x2+n z)

(3)

式中：OS為土壤氧飽和度；x為距蓄水坑中心點的水平距離，cm；z為土壤深度，cm；k、m、n均為擬合參數。

擬合參數見表3。由表3可知，在坑深40、60 cm處理下，模型OS(x,z)的決定系數R2分別為0.955 7和0.808 6，表明該模型能夠較好地描述蓄水坑灌條件下土壤氧氣二維空間分布特征。

表3 不同坑深條件下土壤氧氣二維空間分布擬合參數Tab.3 Two-dimensional spatial distribution fitting parameters of soil oxygen under different pit depth conditions

3 結 論

(1)蓄水坑灌不同坑深處理下土壤氧飽和度隨土層深度增加均呈現指數型下降趨勢。在同一徑向不同深度位置處，60 cm坑深處理的土壤氧飽和度能比40 cm坑深處理下的數值高0.1%～1.98%，處理間數據差異顯著。建立了不同坑深處理下土壤氧飽和度一維垂向指數型分布模型OS(z)，模型決定系數R2為0.869 4～0.983 8，模擬精度較高。

(2)蓄水坑灌不同坑深處理下土壤氧飽和度均隨徑向距離增加呈現先升高后下降趨勢。在同一深度不同徑向位置處，60 cm坑深處理氧飽和度能比40 cm坑深處理高0.1%～1.94%。在60、80、100、120 cm垂向位置處，2處理間數據差異極顯著。建立了不同坑深處理下土壤氧飽和度一維徑向指數型分布模型OS(x)，擬合精度較高。

(3)在同一垂向和徑向位置處，60 cm坑深處理土壤氧飽和度均高于40 cm坑深處理。 40 cm坑深和60 cm坑深處理對土壤氧氣的影響范圍分別為垂向0～40 cm坑壁20 cm范圍、垂向0～80 cm坑壁20 cm范圍，且60 cm坑深處理對土壤氧飽和度的影響范圍更廣。蓄水坑灌條件下，土壤氧飽和度主要與距蓄水坑坑壁的距離及蓄水坑深有關。建立了不同坑深(40、60 cm)處理下土壤氧飽和度二維空間指數型分布模型OS(x，z)，模型決定系數R2為0.955 7和0.808 6，模擬精度較高。