森林生态系统作为地球上最重要的陆地生态系统, 不仅为人类输出木材、药材及其他工业用品, 而且还具有水源涵养、土壤保持、生物多样性保护、气候调节等功能^[1-2]。由于人口的增长、需求的增加、城镇化进程加快、资源开发力度的增大, 导致了森林生态系统的退化, 服务功能的降低, 进而威胁着人类赖以生存的环境^[3-4]。

径流调节主要指养护水资源的举措达到降低水土流失、控制土壤沙化的目的, 对于缓解洪涝灾害, 减少土壤与营养物质流失, 提高系统生产力有着重要影响^[5]。森林生态系统径流调节效应是十分明显, 大多数径流是以壤中流、地下流的形式流出。地表径流调节受多种因素影响, 包括气候条件(降水量、降水强度)、植被特征(类型、结构、林龄、覆盖率), 土壤条件(土壤厚度、孔隙率), 地质、地形特征以及人类活动影响等^[6]。

目前, 我国水资源和水环境问题十分严峻, 水资源短缺, 水质恶化, 洪灾、旱灾频繁, 已给人们生产生活造成了巨大的负面影响。此外, 经济与人口的增长、城市化进程的加快以及全球气候变暖, 将进一步加剧水资源和水环境问题。由于环境问题的加剧, 径流调节在生态系统服务中的重要性得到了广泛的重视。尽管国内外学者针对森林生态系统的径流调节功能进行了大量的研究, 但由于研究区局限和特定, 并未能反映全国尺度的森林生态系统径流调节特征^[7]。本研究通过整合相关文献资料数据, 建立径流系数与植被回归方程, 研究全国尺度的森林径流调节特征, 在探讨森林降低水土流失, 保护水土资源, 以及森林植被恢复和生态重建方面具有重要的科学意义和应用价值。

1 研究数据与方法 1.1 数据收集与筛选

本研究所用数据主要来自于相关文献资料, 包括从国内公开发表的文献和出版的专著(关于径流小区的地表径流、降水、植被数据)。为了充分、有效地利用所收集的数据, 更好地反映各森林生态系统类型地表径流与降水和植被的关系, 本研究根据样本数据获得的方法的统一性对文献资料进行筛选, 从多篇文献资料中筛选满足条件的样本数据。这些数据较大程度地覆盖了我国森林生态系统的分布区域, 代表了我国的各种森林生态系统类型。植被覆盖度数据是基于像元二分模型^[8]通过Modis影像反演得到。像元二分模型是一种简单实用的遥感估算模型, 它假设一个像元的地表由有植被覆盖部分与无植被覆盖部分地表组成, 而遥感传感器观测到的光谱信息也由这2个组分因子线性加权合成, 各因子的权重是各自的面积在像元中所占的比率, 其中植被覆盖度可以看作是植被的权重。全国2000-2010植被覆盖数据(VC)来源于中国科学院遥感与数字地球研究所。生态系统类型图主要基于Landsat数据采用的面向对象的分类技术, 引入非影像光谱信息强化目标的识别能力；生态系统类型数据来源于中国科学院遥感与数字地球研究所。降水与温度数据, 基于普通薄盘和局部薄盘样条函数插值理论, 来源于中国国家计量信息中心/中国气象局(NMIC/CMA)。

1.2 研究方法

中国幅员辽阔, 自然条件复杂。多样化的气候、地形特征, 维持着不同生态系统类型。就森林生态系统而言, 其分布明显反映了地带性的特征, 从北到南, 大致分布为寒温带和温带山地针叶林、温带针阔混交林、温带落叶阔叶林、亚热带针阔混交林、亚热带常绿阔叶林、热带雨林和季雨林。目前, 我国森林生态系统有多种分类方法, 基于本文的研究目的和空间数据的可得性, 本文依据《中华人民共和国植被图(1:100万)》植被分类方法^[9](表 1), 结合全国生态系统类型图^[10-11]整理得到。

利用筛选后的试验数据, 运用SPSS软件, 进行相关分析。基于拟合结果(回归系数R²)选择合适拟合方程, 建立各森林类型地表径流和降水、植被的函数关系, 揭示其径流调节特征。所建立的降水与径流关系应同时满足统计学、生态学的意义(统计学意义为所建立的趋势线能较好地代表所有数据, 关系达到显著性；生态学意义则是在降水范围内, 径流量均应小于降水量)。

径流调节能力WC表示森林生态系统对雨水的调节能力, 计算公式如下：

式中, RC为径流系数(%), 等于地表径流量R与降雨量P的比值乘以100。

2 结果与分析

由于收集试验数据的不完整性, 影响因子的复杂性, 全面了解各因子对径流调节的影响仍有很大的难度。目前, 关于各因子对径流的影响也尚未得到一致的规律性的结论^[7]。所以, 在分析地表径流与影响因子关系时, 应关注主要的影响因素。研究表明, 降水和植被会显著影响地表径流, 其对径流动态的贡献率分别为69.8%和17.3%^[12]。综合考虑, 本文基于数据可获取性和影响因素重要性, 选取森林类型、降水量和植被覆盖作为区域尺度上分析地表径流的关键因子。

2.1 地表径流与降水的关系

森林作为陆地主要生态系统类型, 通过林冠、林下植被和凋落物对降水截留, 以及蒸散作用和森林土壤的贮存特征, 使森林地表径流减少, 起到水文调节作用^[13]。本文基于所收集的数据^[14-60], 统计分析不同森林类型地表径流与降水的关系。由图 1可见, 各森林类型地表径流与降水呈现极显著的正相关性, 其中相关最高的为常绿针叶林, 决定系数R²达到了0.764(表 2)。回归函数多为幂函数和线性函数, 说明随着降水量的增多, 森林地表的径流量也随之加大。降雨明显作用于各森林类型的地表径流, 其影响的解释力为37%-76%。剩余的24%-63%为其他因子对地表径流的影响, 例如土壤厚度、孔隙度、坡度等等。

2.2 径流系数与植被的关系

在样本分析结果中, 除了发现降水量与地表径流显著相关外, 还发现植被覆盖与径流系数存在显著的负相关特征。由图 2可知, 各森林类型植被覆盖与径流系数呈现极显著的负相关性, 其中相关最高的为落叶阔叶林, 决定系数R²达到了0.56(表 2)。回归函数多为幂函数, 说明随着植被覆盖度的增加, 森林的径流系数随之降低。植被覆盖在一定程度上确定着各森林类型的径流系数, 其作用的解释力为32%-56%。

2.3 径流调节能力

虽然我国降雨时间分布差异明显, 但降雨主要分配在夏季(雨热同期), 同时该时段也是植被的主要生长时期, 所以本文选取4-10月的植被覆盖度平均值来代表年植被覆盖度。基于植被覆盖与径流系数的回归函数模拟了全国森林径流调节能力。

统计分析各森林类型的径流系数和径流调节能力, 用均值表示各类型的状况, 对比分析各类型径流调节特征。由图 3可知, 我国森林类型径流系数的范围为0.83%-4.58%, 径流调节能力在95.42%-99.17%之间。就径流调节能力而言, 落叶针叶林流调节能力最高, 为0.83%；其次为落叶阔叶林, 针阔混交林。

3 总结

本文以森林生态系统定位监测数据为基础, 统计分析地表径流与降雨和植被的相关关系, 得出以下主要结论：

(1)各森林生态系统地表径流与降雨呈显著正相关, 随着降水量的加大, 森林径流量主要呈线性和幂函数式的增加, 降水对地表径流的解释能力达到37%-76%。

(2)径流系数与植被覆盖存在显著的负相关特征, 随着植被覆盖度的增加, 森林径流系数主要呈幂函数式的降低, 植被覆盖对径流系数的解释能力为32%-56%。

(3)通过对比分析各森林类型径流调节特征, 揭示出我国森林径流调节能力：落叶叶针叶林>落叶阔叶林>针阔混交林>常绿针叶林>常绿阔叶林。