区域森林水源涵养功能评价模型的开发与应用_余恩旭.pdf

第 卷 第 期 年 月中 国 水 土 保 持 科 学 区域森林水源涵养功能评价模型的开发与应用余恩旭，张明芳，徐亚莉，孙鹏森，蒙作主（电子科技大学资源与环境学院，成都；中国林业科学研究院森林生态环境与保护研究所，北京）摘要：区域尺度上毁林、虫害、火灾以及造林等多种类型森林变化伴随森林植被恢复生长，使森林及其水源涵养功能恢复具有时空累积性。而传统森林水源涵养功能评价方法难以直接刻画森林时空累积变化及其对水源涵养功能的影响，且存在模型计算量大、效率低等问题。针对上述难题，集成自主开发的森林恢复评价模型（模型）与经典的综合蓄水法能力，结合 二次开发平台，设计与实现区域森林水源涵养功能评价模型（）用于不同森林恢复情境下森林景观和森林水源涵养功能的时空动态评价和预测。该模型通过输入降水、植被类型、土壤、森林恢复措施等数据，生成森林变化、水源涵养量变化、恢复措施布局等空间分布图以及相应的统计数据，可以快速、准确定量评估区域尺度森林变化及其恢复过程中水源涵养功能的时空累积变化，辨识影响森林恢复的关键因子。区域森林水源涵养功能评价模型应用于杂谷脑流域，结果表明：）年流域森林总水源涵养量为.；）年间在自然恢复和人工恢复情景下流域森林土壤层、枯枝落叶层、林冠层及总水源涵养量均呈增加趋势；）人工恢复效果较自然恢复更明显，预计到 年人工恢复措施下森林水源涵养总量较自然恢复措施下约高.。本模型的开发及成功应用能够为林业部门和水资源管理部门快速准确评估森林恢复工程对水源涵养功能的提升效果提供高效的评价工具，同时为我国森林资源适应性管理提供科学依据，为水土保持工程的评价和管理提供科学支撑。关键词：森林恢复；水源涵养；二次开发；综合蓄水能力法；水土保持中图分类号：文献标志码：文章编号：（）：收稿日期：修回日期：项目名称：四川省科技厅杰出青年科技人才项目“碳中和背景下川西区域森 林 水 碳 生 态 功 能 的 协 同 与 权 衡”（）；国家自然科学基金“非稳态条件下西南亚高山流域森林恢复的水文响应定量”（）第一作者简介：余恩旭（），男，硕士研究生。主要研究方向：森林生态系统。：通信作者简介：张明芳（），女，博士。主要研究方向：生态水文信息采集与分析，流域生态恢复与规划管理等。：，（，；，）：，中国水土保持科学 年 ，（），（）（），（），（），.，.，.，.，.，：；水源涵养是重要的森林生态系统服务功能之一。但长期的人类活动干扰破坏森林生态系统平衡，造成森林水源涵养功能严重退化，进而引发水土流失等生态环境问题。采取封育、抚育、疏伐等森林恢复措施对受损或退化森林生态系统进行重建和修复是提升森林水源涵养功能、遏制水土流失的有效手段。快速、准确定量评估和预测不同恢复措施对森林植被结构和水源涵养功能的影响有助于筛选有效的森林恢复措施，优化森林恢复模式，为区域未来森林生态系统保护以及水土保持工作提供科学支撑。森林水源涵养功能效益的评价方法主要分为实验法和模型法。实验法多为小尺度上的评价，对多尺度水源涵养功能动态过程的刻画相对有限。模型法较直观地反映森林变化对水源涵养功能的影响，但现有成熟的模型模拟数据需求量大、过程复杂、耗时长，难以实现快速地森林变化的水文效益评价，且水文模型缺少对森林植被生长变化的刻画和预测，无法直接评价森林时空累积变化及其对水源涵养功能的影响，难以满足森林资源管理经营的实际需求。是美国 公司的旗舰产品，为面向矩阵进行运算的交互式数据语言，则为由 开发的遥感影像处理软件。第 期余恩旭等：区域森林水源涵养功能评价模型的开发与应用 为进行 二次开发的功能组件包，其支持 开发语言。和 二次开发和混合编程技术被广泛运用于森林资源动态监测、森林水源涵养功能评价以及空间信息管理与决策等方面。针对水源涵养功能评价，李姗姗基于 模型，运用水量平衡法，采用 和 技术，开发区域水源涵养量估算系统，实现大尺度长时间序列的区域水源涵养量估算；但是目前缺少能够同时模拟和预测森林变化动态以及水源涵养等生态功能的时空变化的综合性快速评价模型。针对上述问题，笔者基于森林等量采伐面积（，）概念，创新性地开发森林等量恢复面积（，）预测模型以用于森林恢复动态定量。该模型不仅能够快速定量时空累积的森林恢复水平，同时考虑森林植被的动态演替和水文功能恢复过程，能够反映生态功能的时空累积变化。笔者通过集成 模型和传统的综合蓄水能力法等水源涵养功能评价方法，基于 二次开发平台，研究、设计与实现区域森林水源涵养功能评价模型（，），并应用于杂谷脑流域。该模型可为林业部门和水资源管理部门快速准确地评价和预测森林生长动态以及森林水源涵养量变化提供技术支撑。模型设计.评价流程区域森林变化与水源涵养功能评价步骤涉及评价数据管理和处理、森林变化模拟与预测、森林水源涵养功能评价与预测。输入数据包括空间数据（、植被 森林类型、降水栅格、林相图、土壤数据等）与非空间数据（森林恢复措施、系数等）；森林变化模拟预测采用基于森林演替过程的 模型；森林水源涵养量采用综合蓄水能力法计算（含林冠层截留量、土壤层蓄水量和枯枝落叶层持水量）；最后可基于森林变化情景模拟预测区域森林水源涵养量变化（图）。图 评价流程图 .模型设计.数据管理和数据处理研究前需要对空间数据和非空间数据进行数据处理和数据管理，以满足模型各部分数据的需求。数据管理模块主要管理文档、空间数据和数据库。数据处理模块主要包括矢量栅格数据处理、遥感影像处理、字段处理、格式转换以及预处理。其中预处理核心功能为日最大降雨查找、土地类型和龄组划分、林地与非林地划分。.森林变化模拟与预测森林变化模拟与预测包含森林区域筛选、森林变化模拟（区域变化和随机变化）和森林变化情景预测。运用森林区域筛选结合地形数据提取特定区域用于森林恢复措施布局。森林变化情景预测通过输入土地覆被数据、优势树种生长曲线、恢复措施和 系数，输出森林恢复措施布局和森林景观变化分布图。森林变化预测的核心为采用 模型刻画森林变化过程。模型的雏形为 模型，模型通过 恢复系数表征不同类型森林干扰（林火、病虫害、采伐等）后森林变化的水文功能变化状中国水土保持科学 年况，用于计算区域尺度各类干扰后森林的时空累积变化。森林等量采伐面积的减少意味着森林等量恢复面积的增加。据此，笔者改进 模型，基于森林植被的水文功能恢复过程构建出森林等量恢复面积（）模型用于表征造林或自然更新后森林的时空累积恢复情况。模型通过建立森林植被生长（年龄或树高等）与水文功能恢复速率的关系来确定 系数，并由此定量不同森林恢复措施作用后流域森林的时空变化，系数等于，表明森林水文功能完全未恢复；系数等于 ，表明森林水文功能已完全恢复。图 为 计算示意图，计算公式如下所示：；（）（）。（）式中：为累积森林等量恢复面积，；为 系数；为恢复林分面积，；为森林等量恢复面积，；为恢复林分编号；为恢复时间长度。基于时间序列计算出各年份不同恢复措施下（自然恢复和人工恢复）的 系数，得到流域 系数曲线。：，；：；：，图 系数计算示意图 .森林水源涵养功能评价与预测该模块通过输入土地覆被、降水、土壤、水源涵养功能参数和 水文功能系数等数据实现对区域森林水源涵养功能现状的评价和预测，采用综合蓄水能力法定量评价森林水源涵养功能现状，计算公式见（）（）。森林水源涵养功能预测是基于初始年份天然针叶林（流域顶级森林群落类型为恢复终极目标）的林冠层、枯枝落叶层、土壤层以及总水源涵养量，结合 模型和文献参数生成各类型森林植被在不同森林恢复措施下的 水文功能系数曲线，该系数定义为等效水文功能恢复系数（，），构建森林水源涵养功能预测模型，实现未来各层次水源涵养量的预测，计算公式见（）（）：；（）；（）；（）；（）；（）；（）；（）。（）式中：、分别为森林生态系统总水源涵养量、林冠层截留量、枯枝落叶层持水量和土壤层蓄水量，；为林冠截留率，；为年日最大降雨量，；为枯落物最大持水量，；为土壤非毛管孔隙度，；为土壤厚度，；为植被面积，为植被类型；、和 分别为第 年 类森林植被 措施组合的土壤层蓄水量、枯枝落叶层最大持水量、林冠层截留量和总水源涵养量，；为第 年 类森林植被 措施组合的 系数；、和 分别为初始年份天然针叶林的林冠层、枯枝落叶层、土壤层以及总水源涵养量，。模型实现笔者基于 和 开发和实现 模型，模型框架和主界面如图 所示，模型数据层主要存储和管理模型的输入数据和相关参数等。业务逻辑层主要提供数据访问、数据处理、模型计算与模拟的功能。表现层是基于、.、和.进行开发，借助 进行 设计与美化，实现软件菜单设计、窗口设计、工作目录、工具箱、基本地图显示和操作、专题制 第 期余恩旭等：区域森林水源涵养功能评价模型的开发与应用图 模型框架 图、图层渲染、统计、图表绘制和时空分析等功能。模型应用笔者以岷江杂谷脑流域数据为实例，应用 模型进行流域森林恢复措施空间布局，评价和预测森林恢复对杂谷脑森林水源涵养功能的影响。本研究涉及到的数据包括林相图（川西林业局 提供）、土地覆被（，空间分辨率，：）、降水（：，插值，空间分辨率）、土壤（，万，：）和地形（，空间分辨率 ，：）数据，恢复措施及水源涵养参数见表。表 不同森林植被类型对应的恢复措施及水源涵养参数 森林类型恢复措施优势树种林冠截留率 枯枝落叶最大持水量 （）非毛管孔隙度 天然针叶林 自然恢复 冷杉 .人工针叶林 疏伐、抚育、密度调整，云杉、岷江柏 .天然针阔混交林 疏伐、修枝 冷杉、桦木 .人工针阔混交林 抚育（割灌除草辅以人工补植）（）云杉、桦木 .天 然 阔 叶 林 自然恢复 桦木.灌木林 清林 高山柳 .注：数据来源于参考文献收集；数据来源于野外实测和实验计算。：，中国水土保持科学 年.森林恢复模拟与预测研究区森林变化以森林恢复为主。依次使用森林变化模拟与预测模块中的森林区域筛选、森林变化模拟功能生成杂谷脑流域森林恢复布局图（图）。该区域森林恢复模式分自然恢复和人工恢复 种模式。自然恢复下，每种森林类型均采用自然演替进行恢复；人工恢复措施下，每种森林类型按照表、图 措施进行恢复。，严重退化灌丛地生态重建与功能提升技术；，箭竹 桦木 岷江冷杉混交林结构调整与功能提升技术；，箭竹 岷江冷杉 粗枝云杉林生态调控与功能提升技术；，岷江柏人工林恢复与功能提高；，自然恢复；，云杉人工林结构调整与水源涵养功能提升技术。：；：；：；：；：；：图 杂谷脑流域森林恢复措施布局图 基于区域样地实验数据，采用“空间代时间”的方法构建不同森林植被类型 恢复措施组合对应的等效水文功能恢复系数（，），运用模型的 模块计算得到不同森林恢复措施下流域等效水文功能恢复曲线（图）。恢复措施实施 年内，自然恢复和人工恢复措施下的森林等效水文功能恢复系数总体均呈上升趋势，且分别达到.和.。结合流域森林恢复布局图和森林演替过程，运用森林变化预测模块，模拟 年起未来森林景观变化情况。经自然恢复 年后，幼龄灌木林、幼龄天然针阔混交林变化较明显，幼龄灌木林全部变为中龄灌木林，幼龄天然针阔混交林全部变为近熟天然针阔混交林；经人工恢复措施实施后，灌木林则全部演替为人工针叶林，其余类型变化与自然恢复下基本一致（图）。图 自然恢复和人工恢复 系数曲线 （）.森林恢复对水源涵养功能的影响评价与预测应用森林水源涵养量评价模块，计算杂谷脑流域 年林冠层截留量、土壤层蓄水量、枯枝落叶层持水量和总水源涵养量分别为.、.、.和.。进一步基于 系数对未来 年流域森林水源涵养量进行预测，结果表明自然恢复和人工恢复下，流域森林林冠层、枯枝落物层、土壤层和总水源涵养量均呈上升趋势。年自然恢复和人工恢复情景下流域森林总水源涵养量相较于 年分别提升.和.。总体而言，人工恢复措施下各层次水源涵养能力较自然恢复高.左右（表）。讨论相比李珊珊开发的区域水源涵养评价系统，模型在水源涵养功能预测、退化森林恢复措施布局等方面具有较强优势。研究中使用的降水、地形、林相图等数据均来自官方渠道，模型构建使用的水源涵养参数来自文章报道和野外实测，使用的方法是应用广泛的经典水源涵养计算方法，据此认为模型和结果可靠；然而，由于水文参数