黃土高原典型退耕草地植被特征對土壤入滲過程的影響

蔣忙舟

(軌道交通工程信息化國家重點實驗室(中鐵一院)，西安 710043)

土壤入滲性能對于降水、蒸發、下滲、徑流等水文循環過程以及坡面產流產沙過程有著重要的影響。土壤水分入滲過程往往受到植被特征和土壤性狀等多種因素的影響。植被地上部分可通過攔蓄降雨延緩產流，削弱降雨動能，抑制雨滴打擊造成的土壤孔隙堵塞，增加地表糙度、降低流速、增強入滲時間等方式來提高土壤入滲能力，植物根系則通過生長穿插分割土壤造成的土體裂隙來增強土壤入滲能力。植被也可通過改善土壤性狀來影響土壤入滲過程。植物生長過程中根系死亡及枯落物降解過程可降低土壤容重，增加土壤黏結力，提高土壤有機質，促進團聚體的形成，增加土壤孔隙，進而改善土壤結構，增強土壤入滲能力。

黃土高原是我國生態安全屏障，也是世界上水土流失最為嚴重的地區之一。充分利用黃土入滲速率強，實施“攔蓄降雨，就地入滲”是該地區水土流失治理的重要措施。1999年“退耕還林(草)”工程的實施，植被得以顯著恢復，植被特征及其對土壤性狀的影響也發生顯著性變化，導致其對土壤入滲能力的影響也存在很大差異。然而植被恢復影響土壤入滲能力的相關研究多集中于土地利用類型差異，且多以生態系統為研究單元，植被演替過程中土壤入滲能力變化相關研究仍不系統，尤其是演替過程中典型物種對土壤入滲過程的影響相對薄弱?；诖?，本研究選取黃土高原撂荒演替后期白羊草和鐵桿蒿草本植被，系統研究不同種植密度及其導致的植被特征和土壤性狀差異對土壤入滲過程的影響，量化初始產流時間、入滲速率和入滲量與植被特征及土壤性狀的相互關系，構建基于植被特征及土壤性狀的典型草地土壤入滲量模擬方程，以期為揭示植物群落的水土保持功能、評價植物群落促進生態恢復效果提供科學依據。

1 材料與方法

1.1 試驗處理

選取黃土高原撂荒演替后期典型物種白羊草和鐵桿蒿，采用穴播種植的方式，按照5，10，15，20，25，30株/m種植密度進行種植。經過1個完整的生長季，采用網格法測定不同種植密度下白羊草和鐵桿蒿植被蓋度及其土壤生物結皮蓋度(表1)。

白羊草和鐵桿蒿均種植于土槽內(長×寬×高為2 m×0.5 m×0.5 m)，每個種植密度設置2個重復，并設置1個裸地對照，共計26個土槽。種植前對土槽進行填土，土壤容重控制在1.2 g/cm。土壤類型為黃綿土，取自延安市安塞區坡耕地表層0—40 cm。填土前將供試土壤過2 mm篩，去除草根和礫石。為保證土壤水分均勻下滲，土槽底部鋪10 cm細砂，并用紗布隔開土和砂。填土過程中，上部土壤分4層填裝，每層填土時鋪平、壓實，并在填裝下層土壤之前將表土打毛，消除土壤之間的分層現象，填土厚度為40 cm。

填土結束后按照預設種植密度種植白羊草和鐵桿蒿，種植方式為穴播。按種植密度將土槽劃分為不同大小的正方形格子進行種植，在格子中心位置挖取0.2～0.3 cm深度的種穴進行播種。種植結束后，將土槽放置于光照充足的地方，并定期灑水，為種子萌芽提供適宜水熱條件。

1.2 降雨入滲過程

采用人工模擬降雨系統進行降雨入滲試驗。降雨設備采用組合式噴頭，高度18 m，雨滴大小0.4～6.0 mm，降雨均勻度>80%，雨滴終點速度近似天然降雨。降雨強度設置為60 mm/h，降雨歷時60 min。降雨前試驗土槽坡度調整為15°。降雨開始后，記錄初始產流時間，降雨過程中收集徑流泥沙，降雨前10 min每隔2 min收集1次徑流泥沙，后50 min每隔5 min收集1次徑流泥沙。降雨結束后采用袖珍式黏結力儀測定土壤黏結力，每個土槽重復6～9次；靜置3天后測定土壤容重(環刀法)，重復3次，并采集土壤混合樣用以測定土壤有機質(重鉻酸鉀外加熱氧化法)，最后挖取整個土槽植物根系測定根系生物量(水洗法，表1)。土壤理化性見表2。

表1 不同種植密度下植被特征的變化

表2 不同種植密度下土壤理化性狀的變化

1.3 土壤入滲參數指標計算

土壤入滲速率和入滲總量均是表征土壤入滲性能的重要參數，可分別按公式(1)和公式(2)計算。

(1)

(2)

式中：為土壤入滲率(mm/min)；為降雨強度(mm/min)；為土槽坡度(°)；為降雨時間內產生的徑流量(kg)；為土槽面積(m)；為入滲歷時(min)；為累積入滲量(mm)。

Kostiakov土壤入滲模型能夠很好地擬合實測入滲過程，基于上述測定土壤入滲速率及對應時間，可按公式(3)擬合出入滲系數()和衰減系數()，用以評估植被特征對土壤入滲過程的影響。

(3)

式中：為時刻的入滲速率(mm/min)；為第1時間單位末尾時的入滲速率(mm/min)；為入滲歷時(min)；為時刻的累積入滲量(mm)；、為經驗入滲系數。

2 結果與分析

2.1 不同種植密度下初始產流時間的變化

不同種植密度坡面初始產流時間存在顯著差異(圖1)。白羊草草地和鐵桿蒿草地不同種植密度下初始產流時間分別為0.76～5.74，0.87～2.08 s。對于白羊草草地，初始產流時間在最大種植密度(30株/m)時最大，是其他種植密度1.15～7.60倍；鐵桿蒿草地初始產流時間在種植密度25株/m時最大，是其他種植密度的1.04～2.40倍?？傮w而言，隨著種植密度的增加，白羊草草地和鐵桿蒿草地坡面初始產流時間均呈冪函數增加(≥0.84，<0.01，圖1a、圖1b)。由于植被類型的差異，白羊草草地和鐵桿蒿草地其初始產流時間也存在顯著差異。白羊草草地最大種植密度下初始產流時間是最小種植密度的7.6倍，鐵桿蒿草地最大種植密度初始產流時間是最小種植密度的2.4倍。這一結果表明，相對于鐵桿蒿草地，白羊草草地可顯著延長初始產流時間。當種植密度≥25株/m時，白羊草草地初始產流時間隨種植密度依然表現為增加趨勢，而鐵桿蒿草地則趨于平穩，甚至略有降低，這在一定程度上表明白羊草草地在延緩初始產流時間方面仍有潛力。除種植密度為5株/m時白羊草草地初始產流時間略低于鐵桿蒿草地外(12.71%)，其他種植密度下白羊草草地初始產流時間均不同程度高于鐵桿蒿草地，是其1.27～2.86倍?？傮w而言，白羊草草地不同種植密度可平均延緩初始產流時間1.42 s，是鐵桿蒿草地的1.78倍。

圖1 不同種植密度下白羊草草地和鐵桿蒿草地初始產流時間的變化

2.2 不同種植密度下坡面入滲過程

白羊草和鐵桿蒿草地土壤入滲速率均在降雨前10 min迅速減小，隨后減小幅度減緩，降雨30 min后趨于穩定，即土壤穩定入滲速率(圖2)。白羊草和鐵桿蒿草地降雨30 min后土壤穩定入滲速率分別為0.45～0.82，0.34～0.60 mm/min，分別較降雨前10 min土壤平均入滲速率減少7.71%～16.98%和17.25%～25.37%?？傮w而言，土壤入滲速率隨降雨歷時的增加呈冪函數降低關系(0.219≤≤0.89，<0.01或<0.05)。在不同種植密度下，白羊草草地和鐵桿蒿草地土壤穩定入滲速率和土壤平均入滲速率均在種植密度為20株/m時最大，分別為0.82，0.82 mm/min和0.67，0.60 mm/min。其中白羊草草地在種植密度為20株/m時土壤平均入滲速率和土壤穩定入滲速率是其他種植密度的1.01～1.73，1.1～1.80倍；鐵桿蒿草地在種植密度為20株/m時土壤平均入滲速率和土壤穩定入滲速率是其他種植密度的1.24～1.77，1.26～1.76倍；此外，白羊草草地和鐵桿蒿草地土壤穩定入滲速率和土壤平均入滲速率均隨種植密度的增加呈增大趨勢，且白羊草草地土壤穩定入滲速率和土壤平均入滲速率均與種植密度存在冪函數關系(≥0.71，<0.01，圖3)。由于植被類型的差異，白羊草草地和鐵桿蒿草地土壤入滲過程也存在顯著差異。白羊草草地土壤穩定入滲率和土壤平均入滲速率在不同種植密度下均高于鐵桿蒿草地，其中土壤平均入滲速率是鐵桿蒿草地的1.16～1.70倍，土壤穩定入滲速率是鐵桿蒿草地的1.27～1.79倍。

圖2 白羊草和鐵桿蒿草地土壤入滲速率隨降雨歷時的變化

圖3 種植密度與土壤平均入滲速率和穩定入滲速率的關系

不同種植密度間入滲量總量存在顯著差異(圖4)。當植被種植密度為20株/m時，白羊草草地和鐵桿蒿草地入滲總量均最大，分別為12.35，10.07 mm，是其他種植密度的1.04～1.74，1.25～1.77倍?？傮w而言，隨種植密度的增加，白羊草草地和鐵桿蒿草地入滲總量均呈增加趨勢，且可表示為種植密度的冪函數。受植被類型差異的影響，白羊草草地不同種植密度下土壤入滲總量顯著高于鐵桿蒿草地，是鐵桿蒿草地的1.16～1.70倍。

圖4 不同種植密度下土壤入滲總量的變化

2.3 植被特征及土壤性狀對土壤入滲過程的影響

種植密度可顯著影響植被特征及其土壤理化性狀，進而對土壤入滲過程產生影響。對于白羊草草地和鐵桿蒿草地，初始產流時間均隨植被蓋度、根系生物量密度和土壤有機質含量的增加呈冪函數增加(≥0.18，<0.05)(圖5)；平均入滲速率、穩定入滲速率及入滲總量均隨土壤有機質含量增加呈顯著的冪函數增加(0.29≤≤0.74，<0.05)(圖6)。Kostiakov土壤入滲模型可用于非均質土壤中入滲速率的模擬，基于公式(3)，通過輸入入滲時間及其對應的土壤入滲速率，計算不同種植密度下白羊草草地和鐵桿草地入滲系數()和衰減系數()。白羊草草地入滲系數()和衰減系數()變化范圍分別為0.36～0.70 mm/min和0.26～0.31，鐵桿蒿草地則分別為0.33～0.67 mm/min和0.29～0.33(表3)。白羊草草地和鐵桿蒿草地入滲系數()和衰減系數()均與種植密度呈增加趨勢，但白羊草草地入滲系數()與種植密度存在冪函數關系；衰減系數()與生物結皮蓋度呈增加趨勢，但并不顯著(>0.05)。白羊草草地和鐵桿蒿草地入滲系數()均隨土壤容重的增加呈冪函數減小，隨土壤有機質含量的增加呈顯著的冪函數增大?？紤]到本研究中土壤容重變化范圍僅為1.20～1.24，容重對入滲系數()的影響并沒有統計學意義?？傮w而言，植被類型及種植密度對土壤入滲的影響可表征為植被特征和土壤性狀的函數關系。逐步回歸結果表明，土壤入滲量()可表征為根系生物量密度(RMD)和土壤有機質(SOC)的擬合函數(公式(4))，其為0.66，NSE為0.65，擬合效果較好。

表3 土壤累積入滲量隨降雨歷時的擬合結果

圖5 植被蓋度、根系生物量密度及土壤有機質與初始產流時間的函數關系

圖6 土壤有機質與土壤入滲性能參數的函數關系

=2.77×RMDSOC

(=0.663，NSE=0.653)

(4)

3 討 論

3.1 植被特征和土壤性狀差異對入滲過程的影響

不同種植密度導致植被特征和土壤性狀的差異是影響坡面土壤入滲過程的主要原因。對于初始產流時間而言，由于白羊草和鐵桿蒿的覆蓋，降雨后并未立即產流，存在不同程度的滯后。這主要是由于植被覆蓋攔蓄降雨、根系生長增加土壤孔隙度和改變土壤理化性狀所造成的。且隨著種植密度、覆蓋度、根系生物量密度和有機質的增加，其延緩初始產流時間的作用增強，初始產流時間可表示為植被蓋度、根系生物量和有機質的冪函數(0.95≥≥0.18，<0.05)。對于平均入滲速率、穩定入滲速率和入滲量總量而言，隨種植密度的增加整體呈增加趨勢。但當種植密度為20株/m時，白羊草草地和鐵桿蒿草地平均入滲速率、穩定入滲速率和入滲量總量均最大，這主要是由于隨種植密度的增加，植物生長受到種間競爭影響，存在一定的水分和養分限制，從而導致植被生長趨于緩慢，這也在一定程度上說明對于不同種植密度導致的植被特征變化，其對土壤入滲的影響存在閾值或上限，當種植密度過大時，土壤入滲性能趨于穩定。

白羊草草地和鐵桿蒿草地平均入滲速率、穩定入滲速率、入滲量總量及入滲系數()均隨土壤有機質的提高而增大，且存在顯著冪函數關系。這主要是由于有機質可促進土壤團聚體的形成，進而使土壤結構進一步改善，土壤孔隙進一步增加，土壤入滲性能增強。作為Kostiakov土壤入滲方程中表征入滲衰減速率的指標，其隨生物結皮蓋度的增加呈降低趨勢，很可能是由于生物結皮遇水膨脹后糙率增大，可降低流速從而加強土壤入滲，其加強作用與生物結皮組成有關。

3.2 植被類型對土壤入滲性能的影響

由于物種的不同，白羊草草地和鐵桿草地入滲過程也存在顯著差異。白羊草在不同種植密度下初始產流時間均在不同程度高于鐵桿蒿草地，其平均初始產流時間是鐵桿蒿的2.01倍。植物根系生長發育過程中通過穿插和擠壓土壤，增加土壤孔隙，從而提高入滲，延緩坡面產流；而相較于直根系鐵桿蒿而言，須根系白羊草根系大都集中分布于表層30 cm土層中，其根系生物量是鐵桿蒿的1.18倍，因而對坡面產流具有較強的延緩作用。

對于平均入滲速率、穩定入滲速率和總入滲量，白羊草草地土壤入滲能力均不同程度高于鐵桿蒿草地。一方面是由于白羊草根系集中分布于表層土壤導致的土壤孔隙度較高；另一方面，生物結皮尤其是苔蘚結皮可顯著降低流速，從而增強土壤入滲。本研究中白羊草草地生物土壤結皮蓋度是鐵桿蒿草地的1.49倍，與此同時，白羊草草地有機質含量是鐵桿蒿草地的1.72倍，這也表明白羊草草地較高的有機質含量能更強地促進土壤團聚體形成，進而增強土壤入滲性能。此外，白羊草相對于鐵桿蒿具有較大根基，其地表起伏/糙度較大，因而也可降低流速增強土壤入滲。對于Kostiakov模型，白羊草草地入滲系數()是鐵桿蒿草地的1.26倍，衰減系數()較鐵桿蒿草地減少6.25%，這與白羊草草地上述土壤入滲參數變化一致，表明白羊草草地入滲性能總體高于鐵桿蒿草地。

4 結 論

(1)種植密度增大可顯著延緩產流時間，白羊草草地和鐵桿蒿草地的初始產流時間和植被蓋度、根系生物量和有機質均存在冪函數關系，且隨植被蓋度、根系生物量和有機質的增加，初始產流時間增大。

(2)種植密度增大可顯著增強土壤入滲能力，過高的種植密度一定程度加劇植物種內競爭，導致種植密度為20株/m時，土壤入滲能力最強。白羊草樣地和鐵桿蒿草地平均入滲速率、穩定入滲速率、入滲量總量及入滲系數()均隨土壤有機質的提高而呈冪函數增大，衰減系數()隨生物結皮蓋度呈降低趨勢。

(3)植被類型差異可顯著影響土壤入滲能力。白羊草草地有較高的根系生物量、生物結皮蓋度和有機質含量，其初始產流時間、平均入滲速率、穩定入滲速率、入滲量總量及入滲系數()均不同程度高于鐵桿蒿草地，衰減系數()低于鐵桿蒿草地?？傮w而言，土壤入滲量可表征為根系生物量密度和土壤有機質的擬合函數(為0.66，NSE為0.65)。