吉林省遼河流域不同植被類型土壤水源涵養能力分析

于忠亮，付世萃，王梓默，潘艷艷，苑景淇，楊 帆，張大偉

(1.吉林省林業科學研究院，吉林 長春 130117； 2.吉林省林業調查規劃院，吉林 長春 130022)

在水資源日益短缺和水質遭受嚴重破壞的背景下，水源涵養功能得到了越來越多的關注，水源涵養林既可以涵養水源，又可以保育土壤[1]。土壤作為植物生長發育的基礎物質，在整個森林生態系統中扮演著重要角色，其質量是衡量退化生態系統功能恢復與維持的關鍵指標之一[2-3]。森林水文運動過程主要包括林冠層、枯枝落葉層和土壤層三部分，林冠層的截留能力與植被郁閉度有關，枯枝落葉層的截留能力受其樹種組成和葉片儲存厚度制約，土壤層的截留能力與土壤孔隙大小有關[4-5]。研究發現地上持水量僅占林分水源涵養能力的15%左右，土壤層占林分水源涵養的85%左右，是涵養水源的主體[6]。土壤具有減少地表徑流、控制水土流失的作用，對植被修復和生態保護有重要意義。

遼河是我國七大河流之一，隨著經濟迅猛發展，水資源的過度開發利用和水體污染的日益加重導致河流兩岸的植被及土壤受到極大的干擾，對原有土壤水源涵養功能造成嚴重影響。眾多學者對不同林型[7]、不同林分密度[8]、近林齡不同種之間[9]等的土壤水源涵養能力進行了研究，然而關于不同植被類型的土壤水源涵養能力卻少有論述。本研究選取吉林省遼河流域內6種植被類型的樣地，研究不同植被類型土壤在0～60 cm深度的結構特征、持水能力和滲透能力，分析不同植被類型土壤的水源涵養能力，為水源涵養林建設提供理論支持。

1 材料與方法

1.1 研究區域概況

遼河流經河北、內蒙古、吉林、遼寧四省(自治區)，吉林省內的遼河流域地處東經122°05′～125°35′、北緯42°37′～44°41′，位于吉林省西南部，流經四平市、雙遼市、公主嶺市、遼源市、梨樹縣、伊通滿族自治縣、東遼縣、龍山區等。屬于典型的溫帶大陸性季風氣候區，年內溫差較大，四季分明，年均氣溫6.15 ℃，年均降水量600 mm，年日照時長2 150～3 000 h。流域兩側以次生林、針闊混交林為主，主要樹種有山楊(Populusdavidiana)、樟子松(Pinussylvestrisvar.mongolica)、水曲柳(Fraxinusmandschurica)、黃刺玫(Rosaxanthina)、紫穗槐(Amorphafruticosa)、蒙古櫟(Quercusmongolica)，主要草本有豬毛蒿(Artemisiascoparia)、馬齒筧(Portulacaoleracea)、鴨跖草(Commelinacommunis)等。

1.2 樣地設置與采樣

2019年7月15—26日對布設在吉林省遼河流域的51個調查點進行了調查，有效調查點位46個(圖1)。依據海拔、坡度坡向、土壤類型、植被生長狀況等選取具有代表性的樣地，包括闊葉林樣地5個、針葉林樣地5個、針闊混交林樣地5個、灌木林樣地8個、草地樣地10個、農田樣地5個，其中闊葉林、針葉林和針闊混交林樣地大小為30 m×30 m、灌木林和農田樣地大小為10 m×10 m、草地樣地大小為5 m×5 m。統計各樣地海拔、坡度坡向、土壤類型、植被組成等相關指標，樣地概況見表1。同時采用100 cm3環刀，按照0～20、20～40、40～60 cm土層深度分層采集各樣地土壤樣本，共采集土樣114個，在實驗室進行土壤指標測定。

圖1 吉林省遼河流域調查點分布

表1 各樣地類型概況

1.3 土壤指標測定

土壤結構特征指標包括土壤容重、土壤總孔隙度、土壤毛管孔隙度和土壤非毛管孔隙度，計算公式分別為

D=(m3-m)/v

(1)

Ptot=(m1-m3)/(m3-m)×D×100%

(2)

Pcap=(m2-m3)/(m3-m)×D×100%

(3)

Pnon=Ptot-Pcap

(4)

上四式中：D為土壤容重，g/cm3；m、m1、m2、m3分別為環刀質量、浸泡12 h后環刀及土壤的質量、靜置5 h后環刀及土壤的質量、烘干后環刀及土壤的質量，g；v為環刀體積，cm3；Ptot、Pcap、Pnon分別為土壤總孔隙度、土壤毛管孔隙度、土壤非毛管孔隙度，%。

土壤持水能力指標包括土壤含水量、最大持水量、毛管持水量，計算公式分別為

Cw=(m0-m3)/(m3-m)

(5)

Cmax=(m1-m3)/(m3-m)×103

(6)

Ccap=(m2-m3)/(m3-m)×103

(7)

上三式中：Cw為土壤含水量，%；m0為環刀及土壤的質量，g；Cmax、Ccap分別為最大持水量、毛管持水量，g/kg；其他符號意義同上。

土壤滲透能力指標包括土壤滲透速率和土壤滲透系數。土壤滲透速率依據《森林土壤濾率性的測定》(LY/T 1218—1999)計算；根據達西定律，分別計算初始溫度時土壤滲透系數(K0)和10 ℃時土壤滲透系數(K10)。

1.4 數據處理及水源涵養能力計算

采用Excel 2020和SPSS 21.0軟件對外業數據進行整理，對相關數據進行單因素方差分析(one-way ANOVA)和Pearson相關分析。采用坐標綜合評定法評定不同植被類型土壤的結構特征、土壤持水能力和土壤滲透能力三項指標，分別用p1、p2、p3表示，三項指標值相加得到不同林型土壤的水源涵養能力指標(p)。坐標綜合評定法基于多維空間理論，將評定指標視為多向量決定的空間點，比較各點與最佳點間的距離，可以對量綱不同的指標進行綜合評價，指標值越小代表功能越優。計算公式為

dij=aij/mj

(8)

(9)

式中：dij為原始數據的相對值；aij為原始數據；mj為每個指標中的最優值；pi為第i個數據到標準點的距離。

2 結果與分析

2.1 不同植被類型土壤的結構特征

土壤容重和孔隙度是反映土體松緊虛實的重要指標，不僅會對土壤通氣性、透水性、植物根系產生影響，也是土壤水源涵養能力的集中體現。由表2可知，在0～20、20～40、40～60 cm土層，土壤容重最大的是針闊混交林，土壤容重分別為1.37、1.41、1.45 g/cm3，土壤容重最小的是農田，土壤容重分別為1.20、1.26、1.18 g/cm3。闊葉林、針葉林、針闊混交林、灌木林和草地的土壤容重隨著土層深度的增加而增大，而農田并未表現出此規律。在0～20、20～40、40～60 cm土層，土壤總孔隙度、毛管孔隙度最大的均為闊葉林，最小的均為農田；在0～20、20～40、40～60 cm土層，非毛管孔隙度最大的為灌木林，最小的為針闊混交林。

2.2 不同植被類型土壤的持水能力

由表3可知，在0～20、20～40 cm土層，含水量最大的均是闊葉林，含水量最小的分別是草地、農田。在40～60 cm土層，含水量最大的是農田，最小的是針葉林。在0～20、20～40 cm土層，最大持水量最大的是闊葉林，在40～60 cm土層，最大持水量最大的是草地，3個土層最大持水量最小的均是農田。在0～20、20～40、40～60 cm土層，毛管持水量最大的是闊葉林，最小的是農田；在0～20、40～60 cm土層，非毛管持水量最大的是灌木林，在20～40 cm土層，非毛管持水量最大的是農田，在0～20 cm土層，非毛管持水量最小的是針闊混交林，在20～40、40～60 cm土層，非毛管持水量最小是針葉林。除農田、灌木林外，其余樣地土壤含水量隨土層深度的增加而減小。闊葉林和針闊混交林在土層>40 cm時的最大持水量下降，針葉林、草地和農田在0～60 cm土層的最大持水量相對穩定，灌木林在土層>20 cm時的最大持水量下降顯著。闊葉林、針葉林和針闊混交林的最大持水量和毛管持水量隨土層深度的增加而減小。

表3 不同植被類型土壤的持水能力

2.3 不同植被類型土壤的滲透能力

由表4可知，土壤滲透能力表現為針闊混交林>針葉林>灌木林>闊葉林>草地>農田，滲透速率分別為11.5、8.75、8.10、4.13、3.25、1.25 mm/min，初始溫度時的土壤滲透系數分別為6.08、4.01、3.49、1.99、1.56、0.91，10 ℃溫度時土壤滲透系數分別為0.72、0.53、0.48、0.38、0.30、0.15。針闊混交林的土壤滲透速率和滲透系數明顯高于其他植被類型。

表4 不同植被類型土壤的滲透能力

2.4 不同植被類型土壤的水源涵養能力

基于坐標綜合判定法評定不同植被類型的土壤結構特征、持水能力、滲透能力和水源涵養能力(表5)。不同植被類型土壤的結構特征指標值表現為闊葉林>針闊混交林>針葉林>灌木林>農田>草地。不同植被類型的土壤持水能力表現為闊葉林>針闊混交林>針葉林>灌木林>草地>農田，滲透能力表現為針闊混交林>針葉林>灌木林>闊葉林>草地>農田，三項指標值相加得到的水源涵養能力表現為闊葉林>針闊混交林>針葉林>灌木林>草地>農田，闊葉林的土壤水源涵養能力最優，農田的土壤水源涵養能力最差。

表5 不同植被類型土壤的結構特征、持水能力、滲透能力和水源涵養能力

3 討 論

除農田外其他植被類型土壤容重均表現為隨著土層深度增加而增加，孔隙度隨著土層深度的增加而減小。土壤結構特征主要是由外界環境因素和植被類型決定的，本研究的所有樣地均位于吉林省遼河流域內，外界因素對土壤結構特征的影響很少，可以認為土壤結構特征的差異是植被類型不同造成的。闊葉林、針葉林、針闊混交林、灌木林、草地的土壤上層腐殖質含量隨著枯落物腐爛分解而增加，表層土壤形成土壤團粒體結構，導致土壤孔隙度增大，土質疏松，土壤容重變小。由于農田種植的通常為玉米和大豆，二者均為淺根系植物，枯落物腐爛分解較少，對土壤結構特征的影響較小，因此并未呈現上述規律。這個結果與劉西剛等[10]、時鐘瑜等[11]提出的森林土壤結構特征隨土層變化的規律一致。

土壤水分是植物生長發育所必須，在促進水循環方面起重要作用，不同植被類型土壤的持水能力存在明顯差異。試驗結果表明，隨著土層深度的增加，土壤含水量和最大持水量逐漸減少，在同一土層中，闊葉林的土壤持水能力最強，農田的土壤持水能力最弱。林地的枯枝落葉中含有大量的木質素、酚類物質化合物，不易被分解吸收，使土壤生物可降解性低[12]，可以有效減少地表徑流，促進更多的降水滲入土壤。然而，農田長期受耕種，土壤層的原有結構遭到破壞，表層土壤呈現出退化、板結、孔隙度降低、持水能力變差。流域內實施退耕還林可以有效改善土壤持水能力。

土壤滲透能力是評價水源涵養能力的重要指標，與土壤結構、孔隙度、土壤濕度和溫度有關[13]。試驗結果表明不同植被類型的土壤滲透能力差異較大，灌木林、針葉林和針闊混交林的滲透能力較好，農田的滲透能力較差。由于不同植被類型土壤的地上凋落物和地下根系生長發育不同，因此造成地表土壤結構和孔隙度的差異，進一步引起土壤滲透能力的差異。劉宇等[14]發現土壤的滲透速度與土壤孔隙度呈正相關關系，本研究也呈現類似的關系，但未達到一致，筆者認為評價土壤的滲透能力需要綜合考慮土壤濕度、溫度和0～60 cm土壤層的結構特征。

4 結 論

(1)除農田外其他樣地類型呈現土壤容重隨土層深度的增加而增大，土壤總孔隙度和毛管孔隙度隨土層深度的增加而減小的規律，同一土層中土壤總孔隙度和毛管孔隙度大小基本表現為闊葉林>針闊混交林>針葉林>灌木林>草地>農田。

(2)不同植被類型土壤的持水能力表現為闊葉林>針闊混交林>針葉林>灌木林>草地>農田，土壤滲透能力表現為針闊混交林>針葉林>灌木林>闊葉林>草地>農田。

(3)由土壤結構特征、持水能力、滲透能力綜合得到土壤水源涵養能力，不同植被類型的土壤水源涵養能力表現為闊葉林>針闊混交林>針葉林>灌木林>草地>農田，闊葉林的土壤水源涵養能力最優，農田的土壤水源涵養能力最差。因此，流域內實施退耕還林措施可以改善水源涵養能力。

(4)本研究僅對該流域內不同植被類型的土壤水源涵能力進行評價，為了全面評價區域水源涵養能力，應該對植被林冠層、枯落物層及林地產流產沙進行長期監測。