岷江雜谷腦流域典型天然林和人工林林地水文效應(yīng)研究

鄧力濠 張明芳,2 師嘉祺 蔣志魏 張 晟 余恩旭 楊 晨

（1. 電子科技大學(xué)資源與環(huán)境學(xué)院，四川 成都 611731；2. 電子科技大學(xué)信息地學(xué)中心，四川 成都 611731）

雜谷腦流域地處青藏高原東緣橫斷山系北段高山峽谷地帶的岷江上游區(qū)域。該地區(qū)的天然林資源豐富，在涵養(yǎng)水源、調(diào)節(jié)河川徑流、改善區(qū)域氣候、防止土壤侵蝕、減少地質(zhì)和洪澇災(zāi)害等方面起著重要的作用[1]，是長江上游和成都平原重要的生態(tài)屏障地區(qū)，也是長江上游地區(qū)和成都平原的重要水源涵養(yǎng)區(qū)[2]。由于20世紀(jì)中葉以來受大規(guī)模的采伐以及地質(zhì)災(zāi)害的影響，雜谷腦流域森林生態(tài)系統(tǒng)開始退化，引發(fā)了水土流失、河道淤積等一系列的生態(tài)環(huán)境問題。隨著天然林保護(hù)工程和退耕還林工程的啟動，該區(qū)域森林生態(tài)系統(tǒng)得到較好的保護(hù)，進(jìn)入森林恢復(fù)階段。開展區(qū)域森林恢復(fù)工程的生態(tài)水文效應(yīng)評價，篩選水源涵養(yǎng)林恢復(fù)與重建的最佳模式，對于該區(qū)域森林生態(tài)保護(hù)以及長江上游森林生態(tài)功能恢復(fù)與提升具有重要意義。盡管在岷江上游開展人工林或天然林地被物和土壤持水性能的研究由來已久[3-5]，但對不同恢復(fù)類型的林地生態(tài)水文效應(yīng)缺乏系統(tǒng)的比較研究。本研究通過野外測定和室內(nèi)實(shí)驗(yàn)分析研究了岷江上游雜谷腦流域米亞羅林區(qū)4種天然林（天然落葉闊葉林、天然常綠/落葉闊葉林、天然針闊混交林、天然針葉林）和3種人工林（按恢復(fù)類型分低郁閉人工針葉林、中高郁閉人工針葉林和人工天然針闊混交林）的林下地被物和土壤層的生態(tài)水文特性。研究結(jié)果對該區(qū)域的人工營林、天然林保護(hù)等相關(guān)森林管理決策的制定具有系統(tǒng)的指導(dǎo)意義，可為岷江上游森林生態(tài)系統(tǒng)服務(wù)功能的恢復(fù)與提升提供參考。

1 研究區(qū)概況

研究區(qū)域位于岷江雜谷腦上游流域米亞羅林區(qū)（31°24′～31°55′N，102°35′～103°40′E）。雜谷腦流域是四川盆地與青藏高原的過渡地帶，具有典型的高山峽谷地貌，上游流域面積4 632 km2，海拔為1 789～5 632 m，平均海拔3 814 m。受高原地形影響，夏季濕潤多雨，冬季寒冷干燥，屬冬寒夏涼的高山氣候。1月份和7月份月平均氣溫分別為8 ℃和12.6 ℃，≥10 ℃的年積溫為1 200～1 400 ℃。年降水量在620～1 470 mm，年蒸發(fā)量為1 000～1 900 mm。土壤特征表現(xiàn)為土層淺薄，土壤貧瘠，供養(yǎng)不平衡，土壤結(jié)構(gòu)松散或板結(jié)，且石礫含量高[6]。土層厚度多在30～50 cm，薄的僅10～20 cm，有白色基巖裸露。土壤類型的分布表現(xiàn)出明顯的垂直地帶性，主要土壤類型為山地暗棕壤、棕色針葉林土和高山草甸土[7]。

米亞羅林區(qū)植被亦呈現(xiàn)明顯的垂直地帶特征，自河谷至山頂依次分布的植被有闊葉林、針闊混交林、針葉林、高山灌叢、高山草甸。喬木優(yōu)勢種有岷江冷杉（Abies faxoniana）、川西云杉（Picea likiangensisvar.balfouriana）、川滇高山櫟（Quercus aquifolioides）、紅樺（Betula albosinensis）、白樺（Betula platyphylla）、山楊（Populus davidiana）、榆樹（Ulmus pumila）。灌木優(yōu)勢種有花楸（Sorbus hupehensis）、絨毛杜鵑（Rhododendron pachytrichum）、峨眉薔薇（Rosa omeiensis）。草本種類主要有窄葉鮮卑花（Sibiraea angustata）、金露梅（Dosiphora fruticosa）、川西翠雀花（Delphinium caeruleum）等。原始森林分布于海拔2 400～4 200 m，以亞高山暗針葉林為主。該地區(qū)天然林在20世紀(jì)60年代間遭到大規(guī)模采伐，而后采伐量逐年減少，直至1998年全面禁伐封育，天然林生境正處于逐漸恢復(fù)的過程中。早期采伐后部分跡地上陸續(xù)進(jìn)行了以粗枝云杉（Picea asperata）為主的人工更新，未進(jìn)行人工更新的采伐跡地最初為草本覆蓋，2～3 a后形成灌叢，10～30 a后則逐漸自然演替為川滇高山櫟、山楊等占優(yōu)勢的天然次生闊葉林。

2 研究方法

2.1 樣地設(shè)置

通過實(shí)地勘察并結(jié)合遙感解譯影像，本研究 以 海 拔 在2 660～3 300 m的7種 典 型 森 林 類型為對象，分別為天然次生常綠/落葉闊葉林（NEDB）、天然落葉闊葉林（NDB）、天然針闊混交林（NCB）、天然針葉林（NC）、低郁閉人工針葉林（LCCP）、中高郁閉人工針葉林（HCCP）和人工天然針闊混交林（DCP）。在各植被類型中設(shè)置20 m×20 m的一級標(biāo)準(zhǔn)樣地（圖1），再將其劃分16塊5 m×5 m的二級樣地，采用GPS定位，調(diào)查樣地內(nèi)喬木層、灌木層和草本層植物種類、高度、郁閉度等因子（表1）。

圖 1 樣地空間分布Fig. 1 Spatial distribution of sample plots

表 1 不同植被類型樣地基本特征Table 1 Basic characteristics of different vegetation types

野外調(diào)查表明各調(diào)查樣地由于森林類型和恢復(fù)演替階段不同以及人類活動的影響，其林內(nèi)植被組成和林下植被狀況有明顯差異。其中人工天然針闊混交林為采伐跡地經(jīng)過人工更新云杉林后，由于立地條件或人為因素（栽植技術(shù)不當(dāng)、未撫育間伐等）導(dǎo)致林地郁閉度降低、質(zhì)量和功能下降、植被組成退化形成的人工林。低郁閉度人工針葉林是指郁閉度低于20%的人工針葉林；中高郁閉度人工針葉林指郁閉度高于20%且未出現(xiàn)退化跡象的人工針葉林。樣地內(nèi)樹種組成比例見圖2。

圖 2 不同植被類型樹種組成比例Fig. 2 Composition proportion of tree species of different vegetation types

2.2 地被物及土壤樣本的采集

在各類型林地20 m×20 m的一級標(biāo)準(zhǔn)樣地中，按對角線法設(shè)置4個50 cm×50 cm的枯落物樣方，用密封袋收集地表苔蘚和枯落物，現(xiàn)場稱質(zhì)量（M1，單位：g），標(biāo)號帶回?？紤]到該區(qū)域林木根系集中分布于0～30 cm土層，30 cm以下土層石礫含量較高，在各林分對角線坡中位置挖土壤剖面，用60 mm3的環(huán)刀分別在0～10、10～20、20～30 cm取3層原狀土壤樣本，每層取3個重復(fù)，裝入鋁盒標(biāo)號帶回分析。

2.3 地被物水文特征測定

將野外采集的7組原狀地被物樣本，分別裝入3 mm孔徑篩網(wǎng)，篩去土壤砂礫，再放入孔徑1 mm的篩網(wǎng)，在清水中浸泡24 h，然后稱質(zhì)量，得到枯落物濕質(zhì)量（M2，單位：g）；最后在65 ℃條件下烘干24 h測定干質(zhì)量（M0，單位：g）。自然含水率、最大持水率、單位面積蓄積量和有效攔蓄量計(jì)算方法分別見式（1）～（4）。

式中：W為有效攔截量（t/hm2）；Rm為最大持水率（%）；R0為平均自然含水率（%）；M為枯落物蓄積量（t/hm2）；A為采樣面積（m2）。

2.4 土壤層水文特征測定

土壤容重、孔隙度、持水率和持水量使用環(huán)刀法1次取樣連續(xù)測定，將裝有原狀土壤的環(huán)刀在水中浸12 h稱質(zhì)量（m1，單位：g）；然后放于干砂上2 h，此時環(huán)刀中土壤的非毛管水已全部流出，稱質(zhì)量（m2，單位：g）；最后對環(huán)刀中土壤再次取樣，放入鋁盒中烘干，將環(huán)刀中的濕土質(zhì)量轉(zhuǎn)換成烘干土質(zhì)量（m0，單位：g）。該方法參考《森林土壤定位研究方法》[8]，計(jì)算方法見（5）～（12）。

式中：W1為單位面積土壤飽和持水量（t/hm2）；W2為單位面積土壤非毛管最大持水量（t/hm2）；D為土壤容重（g/cm3）；C1為土壤最大持水率（%）；C2為土壤毛管持水率（%）；K為土壤總孔隙度（%）；K1為毛管孔隙度（%）；K2為非毛管孔隙度孔（%）；V為環(huán)刀容積（cm3）；m為環(huán)刀重（g）。

2.5 數(shù)據(jù)處理

由于單因素方差分析不適用于本研究問題，本研究使用SPSS 22.0軟件對不同植被類型地被物和土壤的各項(xiàng)生態(tài)水文特征進(jìn)行Kruskal-Wallis檢驗(yàn)。Kruskal-Wallis檢驗(yàn)是非參數(shù)檢驗(yàn)方法，用于多個連續(xù)型獨(dú)立樣本的比較，基本原理是根據(jù)秩和來構(gòu)造檢驗(yàn)統(tǒng)計(jì)量[9]。

3 結(jié)果與分析

3.1 地被物蓄積量

由于生境和立地條件不同，研究區(qū)不同植被類型林地的地被物厚度及儲量存在顯著性差異（圖3）。總體來看地被物厚度和儲量成正比，表現(xiàn)為低郁閉度人工針葉林＜天然落葉/常綠闊葉林＜天然針闊混交林＜中高郁閉度人工針葉林＜人工天然針闊混交林＜天然落葉闊葉林＜天然針葉林，天然林地被物儲量整體較人工林高。在4種天然林中，天然針葉林由于苔蘚層和枯落物腐蝕質(zhì)層厚且針葉難于分解，其林下地被物儲量最大，為16.74 t/hm2；而天然常綠/落葉闊葉林是采伐后形成的天然次生林，林齡短、郁閉度低，林地地被物蓄積量最小，為8.88 t/hm2。3種人工林地被物儲量差異顯著（表2），人工天然針闊混交林由于有大量闊葉落葉和枯枝堆積，林地地被物蓄積量最大，為10.97 t/hm2；中高郁閉度人工針葉林由于林分密度大，枯落物多且分解速度較慢，林下大量枯枝落葉積累，地被物蓄積量與人工天然針闊混交林相當(dāng)；而低密度人工林林下主要以稀疏的草本和灌叢落葉為主，林分密度低，枯枝落葉量少，地被物蓄積量最小，僅為3.27 t/hm2。

圖 3 不同植被類型林下地被物厚度與蓄積量Fig. 3 Thickness and accumulation of understory ground cover in different vegetation types

3.2 地被物持水能力

由表2可知，人工林和天然林的平均自然含水率分別在100.49%～140.64%和126.34%～176.5%之間，最大持水率變動范圍分別為229.05%～279.72和251.15%～347.23%。不同植被類型之間的自然含水率和最大持水率并無顯著性差異，整體看天然林的自然含水率和最大持水率較人工林高。不同植被類型地被物最大持水量表現(xiàn)出顯著性差異，表現(xiàn)為低郁閉度人工針葉林＜人工天然針闊混交林＜天然落葉/常綠闊葉林＜天然落葉闊葉林＜中高郁閉度人工針葉林＜天然針闊混交林＜天然針葉林，其中天然針葉林最大39.67 t/hm2，表明天然針葉林的地被物具有頂級的水文生態(tài)功能，而低郁閉度人工針葉林最小僅7.72 t/hm2，這與其較低的蓄積量有關(guān)。不同植被類型之間有效攔蓄量差異顯著（表3），表現(xiàn)為低郁閉度人工針葉林＜天然落葉/常綠闊葉林＜人工天然針闊混交林＜天然落葉闊葉林＜天然針葉林＜中高郁閉度人工針葉林＜天然針闊混交林，其中天然針闊混交林最大11.2 t/hm2，而低郁閉度人工針葉林最小僅2.58 t/hm2。對比自然含水率占最大持水率的比例和有效攔蓄量占最大持水量的比例發(fā)現(xiàn)，4種天然林間無顯著性差異，而3種人工林間差異顯著。

3.3 土壤持水能力

由表4可知，不同植被類型土壤水文特征變動規(guī)律為，隨土壤深度的增加，土壤容重顯著增大，最大持水率、毛管持水率、總孔隙度、非毛管孔隙度和持水量顯著下降。人工林與天然林之間容重、飽和持水量和非毛管持水量存在顯著性差異（表5），總體上天然林土壤水源涵養(yǎng)能力更強(qiáng)。人工林之間土壤各項(xiàng)水文特征均未表現(xiàn)出顯著性差異，說明3種人工林的土壤持水能力相近。累計(jì)土壤0～20 cm飽和持水量和非毛管持水量，飽和持水量大小排序?yàn)椋禾烊会樔~林＞天然針闊混交林＞人工天然針闊混交林＞中高郁閉度人工針葉林＞低郁閉度人工針葉林＞天然常綠/落葉闊葉林＞天然落葉闊葉林；非毛管持水量大小排序?yàn)椋禾烊会樔~林＞天然常綠/落葉闊葉林＞天然針闊混交林＞人工天然針闊混交林＞低郁閉度人工針葉林＞中高郁閉度人工針葉林＞天然落葉闊葉林。在3種人工林中，人工天然針闊混交林0～20 cm土壤飽和持水量和非毛管持水量都最大。

表 3 不同植被類型地被物生態(tài)水文指標(biāo)差異顯著性檢驗(yàn)Table 3 Significance test of ecological and hydrological indexes of different vegetation types

表 4 不同植被類型土壤生態(tài)水文指標(biāo)Table 4 Soil ecological hydrological index of different vegetation types

4 結(jié)論與討論

地被物的持水能力是反映森林生態(tài)系統(tǒng)水源涵養(yǎng)功能的重要指標(biāo)之一。最大持水量代表地被物對水的最大涵蓄能力，而有效攔蓄量則反映地被物對降雨攔蓄的真實(shí)指標(biāo)[10]。影響地被物有效攔蓄量的因子主要有地被物蓄積量、最大持水率和平均自然含水率。自然含水率取決于光照、溫度、水分、坡度、坡向等立地條件，而最大持水率取決于苔蘚和枯枝落葉的質(zhì)和量。米亞羅林區(qū)海拔較高，年平均溫低于10 ℃，濕度較高且蒸散較低[11]，導(dǎo)致地被物自然含水率較高。各個類型林地地被物蓄積量差異明顯，故林地地被物最大持水量存在顯著性差異，其中天然針葉林地被物蓄積量最大，其最大持水量最高；低密度人工林林地地被物蓄積量最小，其最大持水量最低。天然林的有效攔蓄量平均值高于人工林，自然含水率和最大持水率約為人工林的1.26倍、1.24倍，表明天然林林地被物水源涵養(yǎng)能力更強(qiáng)。對比3種人工林發(fā)現(xiàn)，低郁閉度人工針葉林除自然含水率外，其余地被物水文特征值均低于其他2種人工林，其差異主要來源于地被物的蓄積量；人工天然針闊混交林林地林分結(jié)構(gòu)退化明顯，郁閉度低，地被物表層接收的太陽輻射強(qiáng)，蒸發(fā)量高，同時地被物中樺木（Betulaceaesp.）、漆樹（Toxicodendronsp.）等闊葉落葉比重高，表面積大，易于蒸發(fā)，因此人工天然針闊混交林地被物自然含水率明顯低于其他林地；中高郁閉度人工針葉林的林分結(jié)構(gòu)較其他類型人工林更加完整，郁閉度高于50%，其地被物自然含水率、最大持水率和有效攔蓄量更高。

森林土壤的最大可能蓄水量與土壤總孔隙度密切相關(guān)[12]，而土壤非毛管孔隙是土壤重力水移動的主要通道。因此，土壤的持水能力與土壤的非毛管孔隙更相關(guān)[13]。人工林與天然林之間容重、總孔隙度、非毛管孔隙度、飽和持水量和非毛管持水量存在顯著性差異。天然林尤其是天然針葉林和天然針闊混交林以冷杉老齡林為主，林下根系發(fā)達(dá)，土壤土質(zhì)疏松，孔隙多，有機(jī)質(zhì)豐富，持水性能好；而人工林以中林齡云杉為主，樹種單一，林下植被少，植被根系欠發(fā)達(dá)，營養(yǎng)物質(zhì)循環(huán)功能差，土壤相對貧瘠，且受營林作業(yè)影響土壤板結(jié)嚴(yán)重，土壤持水性能差。天然落葉闊葉林分布于低海拔地區(qū)，土壤石礫含量高，且受道路施工等人類活動影響較多，土壤流失較明顯，因此天然落葉闊葉林相對其他類型森林土壤持水能力最低。在3種人工林中，人工天然針闊混交林0～20 cm土壤飽和持水量和非毛管持水量最大，說明人工天然針闊混交林淺層土壤的水源涵養(yǎng)能力反而較中高/低郁閉人工針葉純林強(qiáng)。3種人工林0～20 cm土壤非毛管持水量均小于天然常綠/落葉闊葉林，由此可見，在人工造林恢復(fù)模式下，盡管人工林地被物蓄積量大，其持水量恢復(fù)快于自然恢復(fù)，但土壤攔蓄降水能力不如天然次生林[14]。

岷江雜谷腦流域植被恢復(fù)正處于重要階段，通過分析不同植被類型地被物和土壤的持水能力差異，能夠?yàn)殚L江上游地區(qū)森林水源涵養(yǎng)效益評估和森林恢復(fù)措施的制定提供科學(xué)依據(jù)。本研究表明天然林地被物和土壤的絕對水源涵養(yǎng)能力明顯高于人工林。3種人工林中，中高郁閉度人工針葉林地被物水源涵養(yǎng)能力最強(qiáng)，人工天然針闊混交林淺層土壤的水源涵養(yǎng)能力更高。從林分的樹種組成來看，中高郁閉度人工針葉林中云杉與冷杉占比為57.4%和42.6%；人工天然針闊混交林中云杉、冷杉和白樺占比為81.0%、9.5%和9.5%。鑒于林地水源涵養(yǎng)的主體為土壤層，因此在人工造林時采用針葉樹種與闊葉樹種混交種植方式，將有利于改善林分結(jié)構(gòu)，增加人工林下植物的多度和豐富度，最終提高人工林土壤攔蓄降水的能力[15]?？傊卺航s谷腦上游地區(qū)在保護(hù)天然林的同時，合理的人工更新針闊混交林將有利于促進(jìn)現(xiàn)有群落盡快向頂級群落恢復(fù)，提升該區(qū)域森林涵養(yǎng)水源、調(diào)節(jié)河川徑流、改善區(qū)域氣候等生態(tài)功能。