摘要:长江上游森林林冠截留量与林分郁闭度呈正相关,当亚高山冷杉林的林分郁闭度为 0.7 时(5~7 月),平均截留率为 24%,当郁闭度在 0.3 时(5~7 月),平均截留率降为 9.5%;从枯落物持水量来看,箭竹冷杉林最大 (6.0mm),藓类冷杉林最小 (2.8mm),主要因为前者有较多的落叶伴生树种和灌木。岷江冷杉原始林的土壤最大持水量、枯落物最大持水量、苔藓层最大持水量比皆伐后形成的其他森林类型要大 2.3~17.2 倍,从而具有更好的水源涵养功能。植被对径流的影响初步结论是,森林大流域的年径流量常常大于少林或无林流域的径流量;不同采伐强度径流量比较是皆伐迹地>择伐迹地>原始森林。与全国其他森林类型的蒸发散研究比较显示,长江上游森林的相对蒸散率较低,为 30%~40%,这主要是由于海拔较高,降水量大所致。
关键词:林冠截留;枯落物持水量;径流量;相对蒸散率
中图分类号:S715.2;
文献标识码:A
文章编号:1000-3037(2001)05-0451-06
长江上游(宜昌以上)占长江流域面积的 50% 以上,控制流域面积约 100万km2,由于其地理位置特殊、地势地貌陡峻、河流比降巨大,降雨丰沛且强度大,长江上游分布的森林植被成为长江流域生态环境保护的绿色生态屏障和社会经济可持续发展的重要基础。尽管长江洪灾发生有气候方面的客观原因,但是人为因素造成的中上游地区森林过度采伐、乱砍滥伐、毁林开荒引起水土流失加剧和江湖、水库淤积,降低了调蓄洪能力,也是造成洪涝灾害的重要原因。长江上游天然林由于遭到长期过量采伐,面积锐减,质量下降,仅川西林区天然林面积在 1985 年至 1995 年的 10 年间减少了 35.1%,而人工更新的森林面积有限,结构不合理,致使长江上游森林的生态功能尤其是水源涵养功能不断削弱,水土流失逐年加剧,江河含沙量急增,严重影响到长江上游社会经济发展,也是造成长江中下游水患的重要原因之一。
因此,研究该区森林植被水文功能的格局与过程,从根本上揭示森林植被对该区水文过程与循环的调节作用刻不容缓。揭示森林水文生态功能及其随自然地理、植被条件和人为干扰施加的环境影响下的变化机制不仅有助于了解森林生态系统中水分的运转过程与机制以及对系统结构、生物地球化学循环、能量代谢和生产力的影响,而且能正确评价和认识森林的作用,为森林合理经营利用,保护自然环境和水资源以及维持人类与自然的协调统一所需生存环境的稳定提供持续发展的科学理论。本文主要从该区森林对大气降水的截留、对地表径流的影响以及森林蒸发散 3 个方面阐述该区森林生态系统水文生态功能和水量平衡的格局与过程。
1 研究点概况
自 1960 年始,中国林业科学研究院和四川省林业科学研究院在川西高山米亚罗林区(北纬 31°40′,东经 102°48′)进行了定位与半定位观测研究。该区地形属高山深谷型,山峰高度多在 4000m 以上,河谷多在 3000m 以下,比高 2000m 左右,坡度多在 35°以上,有的高达 60°~70°。该区土壤属山地棕色灰化土和褐色森林土,孔隙大,渗透性好。地处低纬度、高海拔,受西风带大气环流的影响,年降雨变异为 700~900mm,平均为 857mm,年平均相对湿度 70%。年平均气温为 6.5℃,7 月最高温 31.2°C,1 月最低温 -16.1℃。
由于气候分异,森林植被垂直分布明显。海拔 2500m 以下为落叶阔叶林,2500~3000m 为针阔混交林带,3000~3700m 为亚高山暗针叶林带,海拔 3700m 为高山岷江冷杉、红杉针叶疏林,3900~4400m 为高山灌丛草甸。
2 结果分析
2.1 森林对大气降水的截持作用
2.1.1 林冠截留与降水量的关系
林冠截留是森林水文研究的热点之一。多年来国内外许多学者对此作了大量研究,现已清楚地划分为降雨过程中湿润林冠上的雨水蒸发和雨后林冠蓄水蒸发和树干蓄水蒸发。目前,除获取大量有关林冠截留量和截留率的实测数据外,许多学者根据影响林冠截留的各种因子和林冠截留量的数量关系推导出降水截留的半经验和经验理论模型[1~6]。由于建立的模型有各自的侧重点,导致可比性太差;另外,由于模型要求许多因子的测定,容易产生误差,会降低模型的精度。这些存在的问题主要是缺乏对林冠截流过程的本质认识。刘世荣、温远光[7]等人综合大量的研究结果,认为林冠截留量与降水量存在极紧密的正相关关系,但是不同类型的森林,两者相关的线形是不同的,有的表现为极紧密的直线相关,有的则以幂函数关系较佳。方差分析结果表明,不同森林类型之间林冠截留与降水的关系无显著差异。
冷杉林林冠截留量与降雨量的关系是,林冠截留量随降雨量的增加而增加,截留率随降雨的增加而减小,一年中林冠截留的幅度变化在 20%~70%。林冠最大截留量出现在 2~3 月,最小出现在 5~6 月;截留量与林分郁闭度的关系是,当林分郁闭度在 0.7 时 (5~7 月),平均截留率为 24%,当郁闭度在 0.3 时 (5~7 月),平均截留率为 9.5%(图 1)。
图 1 川西米亚罗亚高山冷杉林林冠降水截留的季节动态
Fig..Seasonal variation of subalpine forest interception in west Sichuan province
2.1.2 枯枝落叶层的持水作用
森林枯枝落叶层不仅能吸滞通过植被而降落在地面之水,而且还能增加地表层的粗糙度,从而减缓水流速度,使降水通过地被物缓慢渗透到土壤中,变地表径流为地下径流,减少土壤表层的流失,因此枯枝落叶层对森林的水源涵养和水土保持具有重要的作用。从表 1 可见,不同森林植被类型之间枯落物持水能力不同,岷江冷杉原始林的枯落物最大持水量、土壤最大持水量、苔藓层最大持水量比皆伐后形成的其他森林类型要大得多,尤其是枯落物、苔藓层最大持水量为其他森林类型的 2.3~17.2 倍和 3.6~6.1 倍,从灌丛到次生林或人工林,最后演替顶极群落,枯枝落叶层、苔藓层、土壤层最大持水量逐渐增大,并朝顶极群落水源涵养能力方向恢复,显示出岷江冷杉原始林具有良好的水文生态功能。
表 1 各森林类型的蓄水能力
Table.1 Water-retaining capacity of different forest types
| 林分类型 | 50cm 土壤最大 持水量 (t/hm2) | 枯落物的最大 持水量 (t/hm2) | 苔藓层最大持 水量 (t/hm2) |
|---|
| 泯江冷杉原始林 | 3666.0 | 39.5 | 44.5 |
| 泯江冷杉次生林 | 3357.5 | 5.8 | 8.7 |
| 粗枝云杉林 | 3504.1 | 2.3 | 7.8 |
| 桦树林 | 3490.4 | 16.5 | 12.3 |
| 灌木林 | 3464.5 | 4.9 | 7.1 |
枯枝落叶层干重的差异,反映了不同林型下凋落物的积累。其中箭竹、杜鹃、高山栎—冷杉林之所以大于藓类—冷杉林,主要因为前者有较多的落叶伴生树种和灌木,加上不同林型的气温、水分状况、土壤酸度的盐基含量等条件的差异,导致枯枝落叶层累积量和苔藓层的变化(表 2)。
表 2 不同林型枯落物厚度、重量及持水量
Table.2 Litter depth, weight and water-holding capacity
| 林型 | 枯分枝落叶层干重 (t/hm2) | 枯枝落叶层厚 (cm) | 持水量 (mm) |
|---|
| 未解物 | 半分解物 | 合计 | 枯枝落叶层 | 苔藓层 | 合计 |
|---|
| 藓类-冷杉林 | 7.8 | 2.1 | 9.9 | 4.5 | 4.0 | 8.5 | 2.8 |
| 箭竹-冷杉林 | 9.5 | 10.1 | 19.6 | 4.5 | 5.7 | 12.2 | 6.0 |
| 杜鹃-冷杉林 | 8.9 | 5.7 | 14.6 | 4.0 | 4.0 | 8.0 | 4.4 |
| 高山栎-冷杉林 | 10.4 | 3.0 | 13.4 | 2.0 | 1.0 | 3.0 | 3.9 |
根据刘世荣等人[7]对全国 79 个点枯落物持水量与现存量的相关分析表明,两者存在极强的直线正相关关系,即:
M=-0.472346+0.331645W (r=0.9266**n=79)
式中,M 为持水量 (mm);W 为枯落物现存量 (t/hm2)。
2.1.3 不同林分类型的土壤理化性质及蓄水能力评价
土壤理化性质与森林状态密切相关,一方面土壤不断供给森林所需的养分,另一方面森林也不断地改造着土壤。森林结构变化,往往引起土壤理化性质的显著变化,而变化后土壤,其理化性质恢复需要较为漫长的时间。当岷江冷杉原始林皆伐后,迹地表土层平均都比原森林覆盖下的土壤减薄 5~10cm;迹地土壤容重自 20cm 以下就大于 1.0,比伐前林地土壤提高 20cm;迹地土壤 pH 值一般比林下土壤提高 0.3~0.8;从有机质的绝对积累数字看,迹地比林地少 49%;迹地土壤的含水率比林下土壤低 36%;迹地土壤的石砾含量较林下土壤要高 3%~5% 以上。可见,迹地土壤的宜林特性是变劣了。但是,随着灌草入侵和森林群落的形成,土壤理化性质逐渐得到改善和恢复(表 3)。
表 3 岷江冷杉原始林采伐后形成的不同林分类型的土壤理化性质
Table.3 Soil physical and chemical properties of different forest types developed after cutting
| 林分类型 | 泯江冷杉原始林 | 泯江冷杉次生林 | 粗枝云杉林 | 桦树林 | 灌木林 |
|---|
| 土层厚度 (cm) | 120 | 110 | 120 | 95 | 85 |
| 腐殖质厚度 (cm) | 5 | 3 | 2 | 8 | 2 |
| 黑土厚度 (cm) | 20 | 8 | 6 | 7 | 4 |
| 容重 (g/cm2) | 0.714 | 0.789 | 0.765 | 0.820 | 0.82 |
| 总孔隙 (%) | 73.1 | 70.3 | 71.2 | 69.0 | 69.3 |
| 非毛管孔隙 (%) | 20.1 | 17.1 | 11.4 | 25.3 | 19.9 |
| 毛管孔隙 (%) | 53.0 | 53.2 | 59.8 | 43.7 | 49.4 |
生物量测定结果表明,岷江冷杉原始林最大,其次是皆伐后经 30~40 年形成的岷江冷杉次生林、粗枝云杉人工林、桦木次生林,最小是皆伐后经 5 年形成的灌木次生林。岷江冷杉原始林的苔藓层生物量高达 7.426t/hm2,是其他类型的 3~7 倍,苔藓层是最为重要的水源涵养功能层,苔藓持水量是其自身重的 6 倍,其土壤的腐殖质层厚、容重轻、总孔隙和非毛管孔隙度大,表现出较高的土壤肥力,良好土壤结构和通气状况,有利于林木生长和涵养水源。
2.2 森林植被覆盖率对地表径流的影响
森林对流域径流量调节受森林植被的类型、年龄、覆盖率的影响,以及气候、地形、地质和土壤的影响。一般的研究结论是,在温暖湿润地区,森林采伐后年径流量增加;在半湿润、半干旱和高温干燥的热带地区,森林采伐后年径流量减少。不同植被类型对径流的影响目前尚属研究热点,并无定论。然而,由于森林采伐及不同采伐强度导致植被覆盖率变化,从而影响地表径流确有比较明显的结论,但现实观测中,径流的产生必须和降雨联系在一起。表 4 列出了不同采伐强度和降雨强度影响下地表径流的动态。
表 4 川西亚高山天然林区降雨强度与径流强度关系表
Table.4 Rainfall intensity and runoff rate of different cutting
types in subalpine forest of west Sichuan province
| 单位:雨强=mm/min;径流强度=mm/m2/min |
| 雨强 | 520 | 369 | 300 | 228 | 191 | 158 | 153 |
|---|
| 皆伐迹地 | 3.6 | 2.16 | 2.19 | 1.92 | 1.95 | 1.35 | 0.87 | | 择伐迹地 | - | 1.18 | 0.55 | 0.44 | 0.42 | 0.36 | 0.36 | | 原始森林 | - | 0.75 | 0.54 | 0.42 | 0.42 | 0.38 | 0.28 |
|
通过表 4 的分析,可以得到以下基本结论,即在径流产生的过程中,径流强度随雨强的增大而增大。由于采伐方式或强度影响了植被的覆盖率,从而使平均径流系数皆伐迹地>择伐迹地>原始森林。这说明,由于森林的采伐而增加了地表径流量。采伐强度越大,径流量增加越大。分别对择伐迹地和皆伐迹地的后期调查表明,伐后植被恢复很快,地表径流量和径流系数逐渐减少。皆伐迹地的平均径流系数为 0.193%,翌年下降为 0.152%;择伐迹地为 0.090%,翌年为 0.076%。由于长江上游地处高山深谷切割很深,气候湿润寒冷,森林的总蒸发量较小,降雨强度不大,所以,大面积的森林流域的年径流量常常大于少林或无林流域的径流量。
对同一流域,不同时期由于森林覆盖率变化不大,因而流域年径流量的变化主要由于当年降雨的影响所致,森林覆盖率的影响并没有完全反映出来(表 5)。但是,不同流域之间对比,则结果非常明显。两流域年降水量变化不大,引起年径流量差别的主要原因是植被。即较高的植被覆盖率具有较高的年径流量。这似乎和同一流域中由于不同采伐措施引起植被变化而改变径流量的结论相互矛盾。而事实上,在表 5 中同一流域不同年代对比的结论也和同一年代不同流域对比的结论相反。森林植被对径流的影响究竟是正效应还是副效应?怎样的对比研究是相对科学的?这恐怕只有等更多的试验验证了。目前的结论可归纳为:①森林大流域的年径流量常常大于少林或无林流域的径流量;②森林采伐改变径流过程观测结果为,岷江流域森林覆盖率下降 10%~15%,地表径流与地下径流的比值由 7:2 变为 8:2,导致洪水量增加,而枯水量减少。
表 5 长江上游森林植被率的变化对流域年径流量的影响(单个流域对比)
Table.5 Comparison of annual valley runoffs under different coverage in the upper reaches of the Yangtze River
流域 名称 | 水文 站名 | 东经 | 北纬 | 面积 (km2) | 年代 (20 世纪) | 森林植被 率 (%) | 年降水 量 (mm) | 年径流 量 (mm) | 森林增建径流 的变化 (mm) |
|---|
| 文锦河斜江 | 跃子站 | 103°29' | 30°46' | 354 | 60 | 87.47 | 1254 | 1168 | |
| 70 | 82.80 | 1454 | 1396 | +228 |
| 80 | 82.81 | 1305 | 1214 | +46 |
| 大邑 | 103°31' | 30°35' | 264 | 60 | 56.89 | 1123 | 783 | - |
| 70 | 54.86 | 1141 | 794 | +11 |
注:文献来源于杨雨行,曾平江.我国森林植被变化对河川径流泥沙的影响.中国林学会.全国森林与水环境学术讨论会论文。宜昌:1993.
2.3 不同类型森林生态系统蒸发散特征比较
对森林蒸散实测研究工作始于 20 世纪 60 年代。以往的研究,森林蒸散量是作为生态系统水量平衡的残余项来进行计算,直接测定少。表 6 列出了自 60 年代开始在我国各地利用能量平衡或水量平衡方法直接测定的森林蒸散量比较,其中包括长江上游四川米亚罗的测定数值,便于和全国其他地方的测定进行比较。
表 6 我国不同类型森林生态系统蒸发散的比较
Table.6 Comparison of evapotranspiration among different types of forest ecosystem