1、第 卷 第 期 年 月中 国 水 土 保 持 科 学 甘肃省植被时空变化及其未来发展潜力靳 峰，戈文艳，秦 伟，韩剑桥，马 涛，郑 安（兰州大学 资源与环境学院，兰州；甘肃省水利厅，兰州；西北农林科技大学，陕西杨凌；中国科学院 水利部 水土保持研究所，陕西杨凌；中国水利水电科学研究院 水利部水土保持生态工程技术研究中心，北京；甘肃省水土保持科学研究所，兰州）摘要：为探究甘肃省生态治理成效及未来植被提升空间，基于全省 年生长季（月）数据及 年土地利用 覆被数据，采用线性回归方程及基于滑动窗口的相似栖息地潜力模型分析全省植被覆盖的时空变化特征及未来植被恢复潜力。结果表明：）年间，甘肃省生长季 以.
2、速率显著上升，年后有所减缓，年以后再次大幅提升，近 年全省植被覆盖整体持续改善。）年甘肃省生长季 总体由东南向西北递减，年均变化速率介于.；区域的植被覆盖有所改善，以陇东和陇南为主的 省域植被覆盖显著快速提升。）陇南与河西分别为全省同期 的高值和低值区，且后者的植被覆盖增加最弱；陇中、陇东的植被覆盖变化相对缓慢。）省内植被恢复指数成效显著，当前植被恢复潜力实现度较高，剩余植被恢复潜力度较低。其中，.的裸地具备转化为绿地的潜力，草地和林地则分别有.和.的面积具备进一步提升的潜力。总之，近 年来甘肃省植被覆盖整体明显改善，但其空间分异特征尚未改变；剔除非自然植被后，全省整体植被覆盖度仅有 左右的提
3、升空间。关键词：植被覆盖；时空变化；基于滑动窗口的相似栖息地潜力模型；生态治理；植被恢复潜力中图分类号：；文献标志码：文章编号：（）：收稿日期：修回日期：项目名称：国家自然科学基金“黄土高原小流域植被格局对侵蚀产沙影响的定量表征与耦合关系”（）；甘肃省水利厅水土保持科研项目“甘肃省水土保持率目标值复核研提与市县分解”（）第一作者简介：靳峰（），男，硕士。主要研究方向：水土保持管理。：通信作者简介：秦伟（），男，博士，教授级高工。主要研究方向：土壤侵蚀与水土保持。：，（，；，；，；，；，；，）：，第 期靳峰等：甘肃省植被时空变化及其未来发展潜力 （），（）（）（）（），.，），.，），）（）（）
4、（），.，.，：；植被作为陆地生态系统的重要组成部分，在保持土壤、调节气候、固碳释氧等方面发挥着重要作用。退耕还林（草）等生态工程的实施显著提升中国西部地区的植被覆盖状况。研究植被变化特征不仅是开展生态与地表过程研究的基础，也可为生态工程效益评估与治理潜力预测等提供必要依据。人类活动和气候变化是植被覆盖变化的主要驱动因素。随着气候变化和人类活动的不断增强，我国北方尤其是干旱半干旱区的植被覆盖发生显著变化 近年来温带大陆性气候区内整体植被趋向良好；大面积实施的退耕还林（草）工程改善黄土高原脆弱的生态环境，促进植被的显著增长，从而减少流域内的输沙量；然而，也有学者指出，干旱半干旱地区的植被虽在整体
5、上呈增加趋势，但在局部区域上却呈现减少的变化趋势。在生境条件差异明显的甘肃省内，植被变化的空间异质性特征显著。甘肃省作为典型气候变化敏感区和生态脆弱区，其水土流失和沙化土壤分布面积大、分布广，是水土保持和生态修复的重点区域。甘肃省大规模实施的国家水土保持重点工程及退耕还林还草等生态治理，使林草植被得到有效恢复，国土绿化面积逐步提升。然而，现有研究多关注省内典型地区的植被演变及影响因素，全区植被恢复的总体状况及不同地区的差异有待探析。此外，对全省植被恢复潜力进行预测，有助于人工调整植被格局以提升甘肃省水土保持功能和效益。鉴于上述背景，笔者采用 年增强型植被指数（，）分析甘肃省近 年来植被覆盖的时
6、空变化特征，并利用基于滑动窗口的相似栖息地潜力模型（中国水土保持科学 年 ，）评估全省植被恢复潜力，以期更加科学地推进甘肃省植被恢复和水土保持等生态建设，为促进我国西北生态屏障安全提供依据。研究区概况甘肃省位于我国西北部（图），全区年均降水介于 ，主要集中在夏季（月），年均气温 ，均呈东南向西北递减的分异格局。境内植被资源丰富且类型多样，包括常绿阔叶林和落叶阔叶林混交带、落叶阔叶林带、森林草原带、草原带、荒漠草原带、荒漠带 个植被地带。受自然区位影响，全省气象和地质灾害等较为频发，是我国生态最为脆弱的地区之一。图 甘肃省地理位置及生态功能分布 数据与方法.数据来源与处理）植被覆盖分析采用来自美
7、国国家航天航空局发布的空间分辨率 、时间分辨率 的（）产品数据集进行。笔者共下载 年的 期影像，并利用 软件进行处理以获得真实 值（）；）地形数据采用中国科学院地理空间数据云平台提供的 数据，分辨率 ，在 软件中提取坡度、坡向；）土壤数据采用中国科学院土壤科学研究所绘制的全国 万土壤类型矢量图，在 软件中按研究区域裁切后转为栅格图；）土地利用数据采用 全球地表覆盖数据 数据库提供的 年土地利用数据。由于人为因素存在较大不确定性，笔者仅研究自然恢复下植被覆盖潜力的改善。为此，以甘肃省 年土地利用为基准，仅保留林草地和裸地进行植被恢复潜力评价，共有 的土地参与后续运算。上述数据参照按 数据分辨率统
8、一重采样为 。.研究方法.最大值合成法生长季植被覆盖率可反映 第 期靳峰等：甘肃省植被时空变化及其未来发展潜力植被生长全过程。为此，通过最大值合成法（，）获得 年甘肃省生长季（月）各月的 影像，并以当年生长季 的均值代表当年植被覆盖的 值：（）。（）式中：为最终采用的第 月的 影像（，）；为当月两景 影像（，）。，。（）式中：为第 年的 影像（，）；，为第 年第 月的 影像。.趋势分析 采用线性回归方程（）对甘肃省 年生长季 数据进行趋势分析。由于全国第 期生态工程于 年结束，第 期生态工程于 年开始，笔者以 和 年为节点，分阶段分析甘肃省植被覆盖的时间变化特征，其显著性水平设为.：。（）式中
9、：为 与时间变量拟合的线性回归方程斜率，和 分别表示 序列随时间递减和递增；为时间变量，取 的整数；为研究时间，；为第 年的生长季月平均。.基于滑动窗口的相似栖息地潜力模型植被恢复潜力（，）是指一定尺度地表单元内的植被覆盖现状与其对应的自然最佳植被状态间的差距。基于滑动窗口的相似栖息地潜力模型（，）是在每个局部窗口建立植被恢复潜力测量模型，能减少空间异质性及缺失量等因素影响，从而提高植被未来预测精度。因此，选择 分片评价小区域的植被恢复潜力指数（，）：，（，）（，）（，）（，）。（）式中：，（，）为第 行、列网格的 值；，为对应网格的环境变量，本研究的环境变量主要包括地形和土壤；（，）（，）为
10、每个网格在 个环境变量中最大的 值；（，）为每个网格的现状 值。为避免数据不确定性，计算 时使用 最大分布的分位值代替最大值，并将 年 值进行综合，构建 最大值，以消除气候的年际变化影响。实际 值与 的比值为植被恢复潜力实现 程 度（，），反映区域植被覆盖恢复潜力：，（，），（，），。（）式中：，为第 年评价区域内第 行、列网格单元的 值；为第 年评价区域内第 行、列网格单元的 值；，为研究时段内评价区域中第 行、列网格单元的最大 值。同时，将其 与实际 值的差值定义为剩余植被恢复潜力度（，），反映当前区域植被恢复可提升程度值的大小：，（，），（，），。（）式中：，为第 年评价区域内第 行、列
11、网格单元的；为第 年评价区域内第 行、列网格单元的 值；，为研究时段内评价区域中第 行、列网格单元的最大 值。结果与分析.时空变化 年甘肃省平均生长季 值为.，随时间呈波动性显著上升趋势（.），其增速为.（图）。在 年前，以.增速显著上升（.），年增速减缓且不显著（.，.），年至今，值增长更为迅速且显著（.，中国水土保持科学 年.）。总体上，近 年甘肃省整体植被覆盖呈明显改善。图 年甘肃省生长季 年际变化 年甘肃省生长季 具有明显空间分异特征（图 和图）。其中，高值区域主要集中分布在南部的草甸草原（）、针叶林（）及落叶与常绿阔叶林（）生态区；图 年甘肃省生长季 空间分异特征 低值区域则广泛分布
12、于西部的山地荒漠（）、高寒荒漠草原（、）和荒漠生态区（、）。整体来讲，陇南较河西地区具有更为优越的生态恢复自然条件。甘肃省 年植被的变化趋势介于.（图），由东南西北逐渐递减，且不同地区间的植被覆盖优劣态势尚未改变（图）。其中，约 区域生长季 变化速率.，表现出显著的快速改善状态（.），主要分布于陇东的黄土高原农业与草原生态区（）、陇南的秦巴山地落叶与常绿阔叶林生态区（）；另有 区域生长季 变化速率为.，集中在河西的内蒙古高原西部、中部荒漠生态区（和）以及塔里木盆地荒漠生态区（）等地；西部荒漠生态区北部（）与草甸草原生态区（）南部则存在少部分植被不显著退化地区。.植被恢复潜力空间特征甘肃省植被恢
13、复潜力指数介于 .，由西北向东南递增，表现为与植被覆盖现状基本一致的空间分布格局（图）。其中，河西地区尤其最北部的内蒙古山地荒漠生态区（）植被恢复潜力较低，说明已达到植被覆盖“天花板”。中部（和）和东南部（和）大部分地区的植 第 期靳峰等：甘肃省植被时空变化及其未来发展潜力图 甘肃省植被恢复潜力空间特征 被恢复潜力均为正值，表现出不同程度的植被恢复潜力。其峰值区集中在中部的内蒙古高原中部草原化荒漠生态区（），相应地区植被改善空间最大。甘肃省约 区域植被恢复潜力实现程度值在.以上（图）。其中，全省的地区植被潜力实现度.，广泛分布在省内东南部的林地、草地与西北荒漠地区。而植被恢复潜力实现程度值.的
14、区域主要分布在中部荒漠草原生态区（）与西部山地荒漠生态区（），说明上述区域植被覆盖度有进一步提升的可能。全省未来植被恢复空间较小（图），地区的剩余植被恢复潜力度值 .。区域剩余植被潜力实现程度值为，说明该区域的植被覆盖已达“天花板”；的区域剩余植被恢复潜力度值在.之间，且主要集中在西部荒漠生态区（如、等），表明该区植被覆盖已接近“天花板”状态；其余 的地区仍具有一定的植被恢复潜力（.），广泛分布在祁连山北部、省内中部内蒙古高原中部（）以及西部北山荒漠区（）。小部分区域（）的潜力实现程度值仅达到.左右，相应区域未来的植被恢复工程建设工作待进一步加强。分析不同地类植被改善程度可知，裸地可改善面积最
15、大（在 中占.），其次为草地（在 中占.），林地的提升空间最小（在 中占.）（图）。.的裸地将具有转变为绿地的潜力，主要集中在祁连山、柴达木盆地、中部草原荒漠生态区的西部及塔里木盆地和内蒙古高原的交界处，相应剩余植被恢复潜力度值主要集中在.以下。林地植被覆盖发生变化的面积较小，仅占其总面积的.，零星分布于省内的东南部黄土高原（）与秦巴山地（）生态区。林地恢复潜力实现程度值也较低（.）。.的草地植被覆盖度具备进一步提高的潜力，但剩余植被恢复潜力度值也较低（.），主要分布在内蒙古高原中部陇中荒漠草原生态区内。甘肃省植被恢复潜力约有的提升空间。全省 可由.上升至.，林草地 可由.上升至.。全省林草地
16、面积占比可由 提升到，表明甘肃省部分地区的植被仍具有一定的恢复空间。中国水土保持科学 年 讨论 年际变化表明甘肃省近 年植被覆盖持续改善，年植被盖度较前期增速降缓，年后再次增速加快。其植被覆盖分阶段的变化趋势主要与生态工程和降水有关，这与基于 的研究 类似。结合我国生态治理进程的综合分析可知，西部各省 年开始实施的首轮退耕还林（草）工程成为甘肃省 年之前植被覆盖改善的主要驱动；年属于降水偏少年份，加之第一期生态工程任务趋于结束，植被覆盖增加有所减缓；年以后，将西北干旱地区作为重点区域的全国第二轮退耕还林（草）工程启动，而同期甘肃省降雨也开始恢复，进一步促进全省植被覆盖的快速改善。这成为 年后全省植被覆盖再次显著增加的主要原因。尽管全省不同地区生长季的 均有不同程度的增加，但其植被覆盖由东南向西北递减总体分布格局及不同地区间的植被覆盖优劣态势并未改变。这主要与全省的气候条件和经济活动紧密相关。陇南地区气候和自然环境相对优越，利于植被生长恢复，始终是全省同期 值最高的地区；中部和东部是省内农业耕作集中区，多山地和草甸，具有植被恢复条件，但植被覆盖变化相对缓慢，主要受退耕、造林、封禁等生态治
