接壩山區陡坡地段闊葉林改造對水文生態環境的影響

塔 莉, 梁文俊

(1.河北環境工程學院 生態學系, 河北 秦皇島 066102; 2.山西農業大學, 山西 太谷 030800)

接壩地區地處灤河上游，對下游地區的樹木生長和水源涵養影響較重，是渾善達克沙地與北京地區之間的一道天然綠色屏障，阻擋著渾善達克沙地向北京地區進軍的風沙，擔負著護衛京津生態安全這一重任，當今社會對生態環境質量的要求日益增高，人們對綠色生活的更加向往，但是目前清新的空氣、清潔的水源和宜人的氣候等生態產品明顯短缺，而植樹造林是我們提高環境質量、保護地球的一種最好、最有效、最直接的方法，森林能夠通過枯落物層與土壤層的生態功能來調節大氣水分循環，從而起到涵養水源、保持水土、凈化水質等作用[1-2]。接壩地區始終致力于森林的培育，植被覆蓋率高，但是坡度較高的陽面陡坡地段仍是植被重建和治理水土流失的重點、難點。目前國內對該地區的研究主要集中在林分結構、生物多樣性、森林健康評價、針葉純林水文效應等方面，而對陡坡地段植被水文生態研究較少。本文以陡坡地段較為常見的兩種植被類型與引針入闊以后的林分作為研究對象，探究其枯落物層與土壤層水文生態效應，旨在為該地區造林與改善水文生態環境提供一定的科學依據[3-4]。

1 研究區概況

研究區設在木蘭圍場國有林場管理局，位于河北省圍場滿族蒙古族自治縣境內，區域地理坐標為北緯41°35′—42°40′，東經116°32′—117°14′，地處渾善達克沙地南緣，與塞罕壩均屬灤河上游地區，是陰山、大興安嶺、燕山余脈的匯接地段，同時也是連接壩上高原與冀北山地的結合部，其地質構造主要趨向于內蒙古臺背斜區，海拔相差較大，從750～1 998 m。這一地區屬半干旱向半濕潤過度、寒溫帶向中溫帶過度、大陸性季風型山地氣候，無霜期67～128 d，年平均氣溫-1.4～4.7 ℃,極端最高氣溫38.9 ℃，極端最低氣溫-42.9 ℃，年均降水量380～560 mm，土壤類型主要是黃棕壤與黑棕壤。研究區內植被主要喬木樹種有華北落葉松(Larixprincipis-rupprechtiiMayr.)、油松(Pinustabulaeformis)、云杉(Piceaasperata)、黑樺(Betuladahurica)、山楊(Populusdavidiana)、山杏(Armeniacasibirica)、蒙古櫟(QuercusMongolica)、白樺(Betulaplatyphylla)等；林下灌木主要有榛子(Corylusheterophylla)、胡枝子(Lespedezabicolor)、沙棘(Hippophaerhamnoides)、繡線菊(Spiraeasalicifolia)等；草本植物主要烏蘇里苔草(Carexussuriensis)、歪頭菜(Viciaunijuga)、寬葉苔草(Carexsiderosticta)、小花油點草(Tricyrtismacropoda)、野艾蒿(Artemisialavandulaefolia)等。

2 標準地設置與研究方法

試驗地設置在木蘭圍場國有林場管理局北溝林場哈叭氣營林區，于2017年6—9月選擇繡線菊灌叢、山杏林與油松—山杏混交林三種植被類型作為研究對象，并在每個植被內設置3塊大小為30 m×30 m的標準地并調查林分的基本情況，其中植被覆蓋度通過實地測定植被的投影面積并結合該林場的遙感地圖來確定(表1)。

表1 研究區3種植被類型標準地概況

2.1 枯落物層水文測定

在標準樣地4個角和中心位置各取面積為50 cm×50 cm的小樣方5個，按照分解程度(未分解、半分解)分層取樣并測定厚度、鮮重。將枯落物帶回實驗室，隨后在設定75 ℃的烘箱進行烘干至恒重，統計各樣地內枯落物的蓄積量。利用室內浸泡法對枯落物的持水能力進行測定，在尼龍網里放置烘干的枯落物并浸在水中，分別測定0.25，0.5，1，2，4，6，8，10和24 h時枯落物樣品重量的變化，從而對枯落物持水量、持水過程進行研究[5-7]。為體現枯落物的真實攔蓄量一般都用有效攔蓄量表示：

W=(0.85Rm-R0)M

式中：W——有效攔蓄量(t/hm2)；Rm——最大持水率(%)；R0——自然含水率(%)；M——枯落物儲量(t/hm2)。

2.2 土壤層水文測定

利用環刀法和雙環法對土壤層的物理性質進行測定，在林分標準地中選取適宜地點進行土壤剖面的挖取，將土壤分為3層，分別為0—10 cm，10—20 cm和20—30 cm，并對土壤容重、非毛管孔隙度、毛管孔隙度等持水能力指標進行測定[8-10]。土壤持水量計算公式為：

W=10 000P·h

式中：W——土壤持水量(t/hm2)；P——土壤孔隙度(%)；h——土壤厚度(m)。

3 結果與分析

3.1 枯落物蓄積量

由表2可知，研究區3種植被類型林分枯落物厚度變化范圍為11.5～14.8 mm，總蓄積量變動范圍為6.24～9.23 t/hm2，其中油松—山杏混交林均為最大，繡線菊灌叢均為最小，未分解層蓄積量均低于半分解層，這種現象是由于枯落物蓄積量影響因素較多，決定枯落蓄積量則是枯枝落葉進入量與分解量，進入量受到林分樹種組成、郁閉度、樣地的水平及垂直結構、本身的厚度和林下植被生長情況等因素影響，分解量則受到樹種組成、溫度、濕度、光照等因素影響，綜上可知引針入闊有利于林分枯落物的形成。

表2 研究區3種植被類型枯落物層厚度與蓄積量

3.2 枯落物水文生態效應

3.2.1 枯落物最大持水量 由表3可知，研究區3種植被類型中油松—山杏混交林最大持水量最大，為26.30 t/hm2，是最大持水量最小繡線菊灌叢(16.61 t/hm2)的1.6倍，這是由于油松—山杏混交林形成的枯落物蓄積量明顯高于繡線菊灌叢；3種植被類型最大持水率變化范圍在265.78%～281.27%之間，依次為：油松—山杏混交林(281.27%)>山杏林(268.13%)>繡線菊灌叢(265.78%)。油松—山杏混交林最大持水量和最大持水率均高于陡坡地區原有的兩種林分，而最大持水量反映了林分攔蓄能力的潛力，表明引入針葉樹種有助于提高該地區最大持水潛力的提高。

表3 研究區3種植被類型最大持水量與最大持水率

3.2.2 枯落物有效攔蓄量 由表4可知，研究區3種植被類型枯落物未分解層有效攔蓄量變化范圍為6.32～7.78 t/hm2，差距不是特別明顯，半分解層有效攔蓄量變化范圍為7.32～13.54 t/hm2，枯落物的兩個層次大小排序一致，均表現為：油松—山杏混交林>山杏林>繡線菊灌叢，在測定林分攔蓄能力時一般采用有效攔蓄量來代表林分真實攔蓄能力，而采用最大持水量代表林分攔蓄能力會導致結果偏高，綜合以上3種植被類型枯落物有效攔蓄量變化規律，可知油松—山杏混交林真實攔蓄能力最強，即在山杏林引入油松后提高了林分枯落物層有效攔蓄能力，枯落物層有效攔蓄量受到林分樹種組成、枯落物蓄積量、郁閉度、溫度、濕度與林下植被生長情況等因素影響，引針入闊要遵循去掉較弱的闊葉樹植入優良針葉樹種的原則，在一定程度上提高了林分整體質量，并且闊葉引入針葉樹種后提高了林分郁閉度，有效增加枯落物層蓄積量對林分有效攔蓄量有一定的促進作用。

3.2.3 枯落物持水過程 由圖1可知，研究區3種植被類型枯落物未分解層、半分解層持水過程變化趨勢較一致，在枯落物持水初期，枯落物較干燥，其吸水速度較快，林分枯落物蓄積量越多，在一定時間內持水的水量越大，試驗中枯落物持水量在浸泡前0.5 h內迅速增加，持水量隨著時間推移逐漸變緩，3種植被類型在浸泡10 h時基本處于飽和狀態，室內浸泡法模擬的持水過程與室外降雨過程中枯落物表現趨勢基本一致，即降雨初期，枯落物由于較干燥攔蓄降雨的功能較強，隨后由于枯落物濕度增加，一定時間內持水量減少，最后枯落物蓄水飽和將不再吸水。3種植被類型枯落物未分解層、半分解層表現的持水規律均為：油松—山杏混交林>山杏林>繡線菊灌叢。

表4 研究區3種植被類型枯落物有效攔蓄量

圖1 研究區3種植被類型枯落物持水量變化過程

對3種植被類型枯落物各層持水量與浸泡時間的關系進行擬合分析，兩者存在較好的函數關系(表5)，具體表達式為：

Q=aln(t)+b

式中：Q，t，a，b——枯落物持水量(g/kg)、浸泡時間(h)、方程系數和常數項。

表5 研究區3種植被類型枯落物持水量、持水率與浸泡時間關系

3.2.4 枯落物吸水速率 由圖2可知，3種植被類型枯落物兩個層次吸水速率變化趨勢較為一致，均與浸水時間呈“倒J”形曲線：枯落物浸泡初期吸水速率非常高，后迅速下降，達到4 h后下降速度明顯減緩，在10 h后沒有較大變化。3種植被類型枯落物未分解層、半分解層在浸泡初期，油松—山杏混交林吸水速率最大，高于其他兩種林分，繡線菊灌叢遠遠低于喬木林分，3種植被類型枯落物未分解層、半分解層吸水速率與浸泡時間的關系進行擬合，兩者具有較好函數關系(表5)，表達式為：

V=ktn(R2>0.99)

式中：V，t，k，n——枯落物吸水速率〔g/(kg·h)〕、浸泡時間(h)、方程系數和指數，其中n是由持水率與時間進行函數擬合得出。

圖2 研究區3種植被類型枯落物吸水速率變化過程

3.3 土壤層物理性質及水文生態效應

3.3.1 土壤物理性質 由表6可知，研究區3種植被類型土壤容重大小排序為：山杏林(1.15 g/cm3)>油松—山杏混交林(1.12 g/cm3)>繡線菊灌叢(1.10 g/cm3)，土壤垂直面上隨著深度的增加土壤容重呈現增大的趨勢。

非毛管孔隙度決定土壤有效持水量的大小，3種植被類型大小排序為：油松—山杏混交林(6.90%)>山杏林(6.09%)>繡線菊灌叢(5.56%)，說明油松—山杏混交林有效持水量最高，這是由于引針入闊時要提前進行土地的整理并且采取減弱干擾喬木與灌木等干擾因素，增加了針葉樹種生長空間，有利于非毛管孔隙度的形成，因此引針入闊提高了土壤層的有效攔蓄量。

表6 研究區3種植被類型土壤物理性質及持水量

3.3.2 土壤層孔隙度與蓄水能力 由表6可知，不同植被土壤結構有一定的差異，土壤孔隙一定程度上影響了土壤通透性與土壤持水量。非毛管孔隙度決定了林分土壤層有效涵養能力的強弱，有效持水量依次為：油松—山杏混交林(69.00 t/hm2)>山杏林(60.87 t/hm2)>繡線菊灌叢(55.60 t/hm2)；飽和持水量依次為：油松—山杏混交林(458.70 t/hm2)>山杏林(433.07 t/hm2)>繡線菊灌叢(394.67 t/hm2)。林分土壤蓄水能力是判定林分水文生態效應的重要指標之一，而決定林分土壤蓄水能力的主要因子是土壤結構，綜合以上飽和持水量和有效持水量表現可知，油松—山杏混交林的土壤持水能力最強，即油松—山杏混交林土壤層水文生態效應最好，說明引針入闊提高了土壤層蓄水能力，改善了水文生態環境。

3.3.3 土壤層的滲透性 由表7可知，3種植被類型土壤層初滲速率差距比較大，油松—山杏混交林初滲速率最大為29.78 mm/min，繡線菊灌叢初滲速率最小為22.38 mm/min，前者為后者的1.33倍，隨著時間的推移，入滲速率逐漸減慢，直至趨于穩滲，3種植被類型穩滲速率變化范圍為0.96～1.29 mm/min，依次為：油松—山杏混交林(1.29 mm/min)>山杏林(1.06 mm/min)>繡線菊灌叢(0.96 mm/min)，林分土壤滲透性越好，其水源涵養能力越強，反映出林分土壤層水文生態效應越好，說明引針入闊有效提高了水源涵養能力，改善了水文生態環境。對3種植被類型入滲速率與時間進行函數擬合(圖3)，兩者具有較好的冪函數關系：

f=at-b(R2>0.97)

式中：f——入滲速率(mm/min)；a，b——常數；t——入滲時間(min)。

表7 研究區3種植被類型土壤滲透的速率及滲透數學模型

圖3 研究區3種植被類型土壤入滲過程

4 結 論

(1) 研究區枯落物蓄積量依次為：油松—山杏混交林(9.23 t/hm2)>山杏林(8.55 t/hm2)>繡線菊灌叢(6.24 t/hm2)；有效攔蓄量依次為：油松—山杏混交林(21.32 t/hm2)>山杏林(18.73 t/hm2)>繡線菊灌叢(13.64 t/hm2)，有效攔蓄量決定了林分枯落物層真實攔蓄能力，但是受到的影響因子較多，與蓄積量、自然含水量和降雨特征等均有關；枯落物持水量、枯落物持水速率與浸水時間存在較好的函數關系，與陳波等[11]研究結果相一致。

(2) 研究區3種植被類型土壤容重隨著土壤深度的增大而呈現增大的趨勢，非毛管孔隙度則與之相反，這與涂志華等[12]研究結果一致；3種植被類型中油松—山杏混交林土壤層飽和持水量與有效持水量均為最大，主要取決于油松—山杏混交林土壤的孔隙特性，土壤結構決定林分土壤蓄水能力，土壤蓄水能力反映了林分土壤層水文生態環境是否良好，研究中油松—山杏混交林的土壤持水能力最強，即油松—山杏混交林土壤層水文生態效應最好，說明闊葉林改造提高了土壤層蓄水能力，改善了水文生態環境；3種植被類型中油松—山杏混交林初滲速率最大為29.78 mm/min，繡線菊灌叢初滲速率最小為22.38 mm/min，林分土壤滲透性越好，其水源涵養能力越強，反映出林分土壤層水文生態效應越好，這說明闊葉林改造有效地提高了其水源涵養能力，改善了水文生態環境；3種植被類型土壤入滲速率與入滲時間之間存在較好的冪函數關系(R>0.97)，這與韓雪成等[13]的研究結果一致。

(3) 綜合枯落物層與土壤層持水能力可知，3種不同植被類型中油松—山杏混交林水源涵養功能最強，而繡線菊灌叢水源涵養能力最差，說明在陡坡地段闊葉林改造在一定程度上增強了林分水源涵養能力，改善了林分水文生態環境，因此從增強涵養水源的角度出發，建議在該地區灌木林與闊葉林中引入適當密度油松，從而改善該地區水文生態環境。

(4) 影響枯落物有效攔蓄量的因子較多，包括蓄積量、自然含水量和降雨特征等。本研究未將主要因素進行明確，需要進一步的進行試驗來確定。但是本研究只引入了油松，未將該地區的主要樹種華北落葉松與云杉進行引入，需要進一步引入探究。