生态修复前后祁连山地区植被覆盖变化
2023-02-03吴晶晶杜达石朱许丽车曦晨
吴晶晶,焦 亮,*,张 华,杜达石,朱许丽,车曦晨
1 西北师范大学 地理与环境科学学院, 兰州 730070 2 西北师范大学 甘肃省绿洲资源环境与可持续发展重点实验室, 兰州 730070
植被覆盖动态变化可用于监测生态环境的优劣,实时监测植被覆盖的动态变化可为生态环境的保护提供决策方向和建议[1—2]。相关研究表明植被覆盖的变化是受气候因素和人类活动共同影响的[3],气候因素在大尺度方向影响植被的生长状况,而人类活动则是在小尺度短时间内影响植被覆盖度的变化,其中生态修复工程对植被覆盖恢复具有重大影响。退耕还林工程、天然资源保护工程使得西南地区的森林覆盖面积增大和林木再生,改善了西南地区的森林生态系统[4—5];“三北”防护林和天然林保护工程促进了黄河流域地区的植被增长[6—7];生态补水工程促使河道两侧地下水埋深大幅抬升,遏制了海水倒灌对湿地生态系统的破坏趋势[8];封地育林、退耕还林等工程使城镇生态系统得到显著改善[9—12];石漠化治理计划、草原生态保护计划、森林生态系统补偿基金、水土保持重点建设工程等13项生态修复工程促使了旱区植被变绿、改善了当地的生态系统的恢复力[13]。同时,生态修复工程对山地森林系统也具有重大影响,青藏高原植被覆盖总面积的61.16%是由人类活动引起的[14];天然林保护、退耕还林还草等工程使秦岭地区大范围植被生态系统得到抚育,林地、草地和水域面积大幅度增加[15];适度的生态工程及农业生产活动对黄河流域的植被恢复起促进作用[16]。探究人类活动对生态环境的影响,可以更好地调控其对生态环境的干扰程度,为政府的决策提供有力依据。
祁连山作为中国西部一条重要的生态屏障[17],具有涵养水源、保护森林、调节气候、保持生物多样性、维持生态系统平衡等多重重要的生态功能[18]。目前对祁连山生态环境变化的研究主要包括气候变化、冰川冻土变化、生物多样性变化以及祁连山土壤植被变化等方面。同时,祁连山为保护生态环境也实施了退耕还林还草、天然林保护、冰川环境保护、山水林田湖草沙冰等一系列生态修复工程[19—23],但是,目前生态工程实施后对生态系统的保护和恢复起到的作用还没有一个系统的评估。另外,影响植被覆盖变化的主要因素为气候因素和人类活动[24—25],但气候变化和人类活动对于植被覆盖变化的贡献率还存在一定的不确定性。
鉴于此,本文以祁连山地区为研究对象,利用谷歌地球引擎(Google Earth Engine)云平台的Landsat NDVI数据集,通过JavaScript编写程序代码研究祁连山生态修复前后的植被覆盖度的变化。同时,为了定量化评估气候因素和人类活动的贡献率,采用多元回归残差方法对祁连山地区生态修复前后植被覆盖度的变化进行分析,进而探究生态修复工程对祁连山地区的作用,并为以后的祁连山生态修复方案提供科学依据和发展方向。
1 研究区概况
祁连山位于甘肃省和青海省的交界处(图1),连接黄土、蒙新、青藏三大高原,地理位置位于93.4°—103.4°E,35.8°—40.0°N[26]。祁连山地势西高东低,山脉走向表现为西北—东南方向,平均海拔在4000—5000m。祁连山地区为典型的大陆性高寒半湿润山地气候,冬季寒冷干燥,夏季温凉湿润,年平均气温在4℃以下,全年日照时数1744h,年降水量维持在在400mm左右,且主要集中在5—9月。祁连山水系呈辐射格状分布,且水系多以冰川融水为主要补给[27]。祁连山地具有明显的垂直地带性,主要的植被有青海云杉(Piceacrassifolia)、祁连圆柏(Sabinaprzewalskii)、油松(Pinustabuliformis)、金露梅(Potentillafruticosa)、银露梅(Potentillaglabra)、箭叶锦鸡儿(Tetraogallushimalayensis)、杜鹃(Rhododendronsimsii)等[28—29]。
图1 研究区位置图Fig.1 Location of the study area
2 数据与方法
2.1 数据来源
1986—2020年Landsat NDVI数据集来源谷歌地球引擎云平台[30—31](https://earthengine.google.com/)。该数据源为Landsat TM/OLI影像,空间分辨率为30×30m。通过JavaScript编写程序代码在Google Earth Engine云平台处理数据,包括辐射定标、几何校正、掩膜处理、去云处理、NDVI计算、NDVI年最大合成、像元二分模型植被覆盖度反演等。根据美国地质调查局官方介绍,不同Landsat卫星传感器之间是经过校正的,所有一级Landsat卫星数据之间具有一致性,可以用于时间序列分析[32]。
本文所用到的气象数据来源于祁连山及周边的22个气象站,包括逐年平均气温、降水量、相对湿度及干旱指数。逐年平均气温、降水量、相对湿度数据从中国气象数据共享服务网(http://data.cma.cn/)获得。干旱指数(SPEI)数据来自荷兰皇家气象研究所(http://climexp.knmi.nl/,于2020年10月17日访问),干旱指数是基于气象数据(降水、温度、相对湿度、太阳辐射、水气压和风速)计算得出的。
2.2 研究方法
时间序列划分:本文依据祁连山地区发生的重大生态环境安全保护事件的时间点,即1999年甘肃率先开展退耕还林还草工程和2017年受中央督导组指示开展的生态修复工程,将研究时段1986—2020年划分为三个阶段,即生态修复前(1999年以前)、生态修复中(2000—2016年)和生态修复后(2017—2020年)。
像元二分模型:本文利用像元二分模型合成的最大NDVI近似估算植被覆盖度[33—34]计算公式为:
(1)
式中,FVC为像元植被覆盖度大小,NDVI为像元归一化植被指数值,NDVIveg为纯植被像元的NDVI值,NDVIsoil为纯裸土像元的NDVI值。
根据祁连山地区植被覆盖的实际情况,采用等间距分级法[35],将植被覆盖度分为5个等级,即低盖度(FVC<20%)、中低盖度(20%≤FVC<40%)、中盖度(40%≤FVC<60%)、中高盖度(60%≤FVC<80%)、高盖度(FVC≥80%)。
多元回归残差分析:利用多元回归残差方法分析气候因素和人类活动对植被覆盖度的影响,即贡献率[36—40]。该方法主要有以下3个步骤:①以年FVC为因变量,气温、降水量、相对湿度和干旱指数为自变量建立二元线性回归方程,计算模型中的各项参数;②基于1986—1999年气候因素数据以及回归模型的参数计算得到FVC的预测值(FVCCC),用来表示气候因素对植被NDVI的影响;③计算FVC观测值与FVCCC之间的差值,即FVC残差(FVCHA),用来表示人类活动对NDVI的影响。具体计算公式如下:
FVCCC=a×T+b×P+c×S+d×Q+e
(2)
式中,a、b、c、d、e是常数,T为气温,P为降水量,S为相对湿度,Q为干旱指数,FVCcc为预测值,FVCobs为实际观测值,FVCHA为残差值。
趋势分析法:使用Thril-Sen斜率表示气象因素和人类活动对植被覆盖度的正反向作用[41—44],计算公式如下:
(3)
式中,median表示中位数函数;xi和xj为序列数据,表示不同年份的贡献率;tj和ti为时间序列数据;当β> 0时表示影响因子为正向促进作用;当β< 0时表示影响因子为反向抑制作用。
偏相关系数:本文利用偏相关系数研究植被覆盖度与平均气温、年降水量、相对湿度、干旱指数的关系[45]。当在分析变量x1和x2之间的相关程度时,控制了变量x3的线性作用后,x1和x2之间的一阶偏相关系数定义为:
(4)
式中,r12(3)表示固定自变量3后自变量1与因变量2之间的偏相关系数,其中r12、r13、r23分别为自变量1、自变量2、因变量3之间的相关系数。rxy>0表示两个变量之间呈现正相关,即相关程度强;rxy<0则表示两个变量之间呈现负相关关系,即相关程度弱。
3 结果与分析
3.1 祁连山地区植被覆盖时空变化特征
3.1.1植被覆盖度的年际变化
通过对祁连山植被覆盖度年际变化的分析可知,祁连山地区1986—2020年的FVC总体呈现为波动增加的趋势,增长倾向为0.022/10a。生态修复前后FVC平均值有显著的差异,其中生态修复前祁连山地区FVC表现为减少趋势,减少趋势为0.008/10a;生态修复中祁连山地区FVC为增加趋势,增速为0.030/10a;生态修复后祁连山植被快速增长,增加速率为0.137/10a(图2)。
图2 祁连山地区平均植被覆盖度年际变化图Fig.2 Annual variation of mean vegetion coverage in Qilian Mountains 不同的小写字母表示不同时期之间差异显著(P<0.05)
祁连山地区实行生态修复工程前,1999年高FVC面积为1.35×105km2,仅占0.70%;中高FVC面积为5.11×106km2,占26.46%;中FVC面积为2.64×106km2,占总面积的13.71%;中低FVC面积为3.45×106km2,占总面积的17.87%;低FVC面积为7.97×106km2,占总面积的41.27%(图3)。在实行退耕还林还草生态修复工程过程中,2016年高FVC面积达到了5.34×106km2,占27.65%;中高FVC面积为1.94×106km2,占10.06%;中FVC面积为1.87×106km2,占总面积的9.68%;中低FVC面积为3.24×106km2,占总面积的16.80%;低FVC面积则减少到了6.92×106km2,占35.85%(图3)。实行生态修复工程后,2020年高FVC面积为5.33×106km2,占27.61%;中高FVC面积为1.90×106km2,占9.86%;中FVC面积为2.02×106km2,占总面积的10.49%;中低FVC面积为3.33×106km2,占总面积的17.29%;低FVC的面积为6.70×106km2,占总面积的34.75%(图3)。
从图3中可以看出,生态修复后祁连山地区高等级植被覆盖面积增加,低等级植被覆盖面积减少。其中,高等级FVC(中高、高FVC)面积增加了1.99×106km2,占祁连山总面积的10.3%,低FVC面积减少了1.26×106km2,面积占比减少了6.52%,表明祁连山地区植被向好的方向转化。
图3 祁连山植被覆盖等级分布Fig.3 Grade distribution of vegetation coverage in Qilian Mountains
3.1.2植被覆盖度的空间变化
祁连山地区中高覆盖度分布在东南地区,西北方向地区以低覆盖度为主。自实施生态修复工程后,各个等级的FVC相互转化,且多为低等级FVC向高等级FVC方向转移(图4)。其中在生态修复过程中,2016年祁连山的中高覆盖度面积增长迅速,中低覆盖度向中高覆盖度方向转化(图4);生态修复后,2020年低植被覆盖度面积减少了1.39×104km2,中高FVC向高FVC转移面积达4.64×104km2,转移速率为1.36×102km2/10a(图4)。对比实施生态修复工程前后,总体程度上祁连山地区的植被状态恢复良好,低等级FVC面积减少,高等级FVC面积增加。
图4 祁连山植被覆盖等级转移图Fig.4 The transfer map of vegetation cover grade in Qilian Mountains
3.2 植被覆盖变化驱动力分析
祁连山地区植被覆盖度受气候因素与人类活动共同作用,因此,利用偏相关系数探究植被覆盖度与气候因素的相关关系,其中祁连山地区的FVC变化与降水量、气温、干旱指数表现为显著正相关(P<0.01),与相对湿度为负相关。在研究时段中干旱指数、气温对植被生长影响最大,也就是说植被覆盖度受当地干旱程度制约最大,平均气温对植被覆盖的影响程度表现为增强的趋势,而干旱指数、降水量则为减弱的趋势。在不同研究时段内,影响植被覆盖度变化的主导因子各不相同,其中,生态修复前祁连山植被覆盖度受降水量和干旱指数影响最大,在生态修复过程中以及实施生态修复后,平均气温制约着植物的生长发育(图5)。
图5 气候因素与植被覆盖度(FVC)间的相关系数Fig.5 Correlation Coefficient between climatic factors and vegetation coverage (FVC)*表示在0.05水平上显著相关;**表示在0.01水平上显著相关
为了定量化评估气候因素和人类活动对祁连山植被覆盖度的贡献率,采用多元回归残差分析将气候因素的影响与人类活动的作用剥离,从而得到以下结果:在研究时段内,气候因素对植被覆盖的贡献率呈减少趋势,减少的幅度为1.52/10a;人类活动对植被覆盖的贡献率则为增加趋势,增幅为1.52/10a(图6)。
气候因素和人类活动对植被覆盖的影响分为正向促进作用和反向抑制作用。本文利用Thril-Sen斜率表示气候因素和人类活动对植被覆盖度的正反向作用。研究时段的34年间,其中气候因素有15年对植被覆盖表现为正向促进作用,人类活动则有19年对植被覆盖变化表现为促进作用。生态修复过程中,气候因素对祁连山植被覆盖变化的影响为反向抑制作用,Thril-Sen斜率β为-0.0023,人类活动则为正向促进作用,Thril-Sen斜率β为0.0023;生态修复后,气候因素对祁连山植被覆盖变化的影响为反向抑制作用,斜率β为-0.015,人类活动则为正向促进作用,Thril-Sen斜率β为0.015。这表明人类活动对植被覆盖度的影响作用逐渐增加。同时,实行生态修复工程后,人类活动对植被覆盖的正向促进作用大于反向抑制作用(图7)。
图6 气候因素和人类活动对植被覆盖度的贡献率Fig.6 Contribution rate of climate factors and human activities to vegetation coverage
图7 气候因素及人类活动对植被覆盖度的贡献率斜率Fig.7 Slope of contribution rate of climate factors and human activities to vegetation coverage
4 讨论
4.1 祁连山植被覆盖度时空变化
祁连山植被覆盖度在生态修复前后有明显变化。在时间尺度上,1986—2020年间祁连山地区植被覆盖呈现增加趋势,增加倾向为0.022/10a。祁连山作为中国西部重要的生态屏障,随着生态文明理念的普及,国家和地方意识到生态环境状况决定着该地的经济发展和民生状况,相关部门采取了一系列的生态保护措施,包括退耕还林工程、天然林保护工程与多项水利补偿工程等[46—48],这些工程的实施一定程度上促进了区域植被覆盖度的增加。甘肃省作为1999年退耕还林政策的首批试点,其生态保护工作的实施使得人类砍伐森林、践踏草丛的现象减少,区域植被生长逐渐趋于稳定,这可能是造成植被覆盖变化在2001年前波动较大,而2001年后变化趋势较为稳定的原因。2011年实施的天保工程人工造林项目为祁连山地区植被覆盖的增加奠定了良好基础,截止到2016年祁连山地区的植被覆盖度达到了38.00%。另外,实施生态修复工程后祁连山地区植被快速增长,2020年祁连山植被覆盖度高达46.11%,此结果得益于退耕还林还草、天保工程人工造林项目等对祁连山植被恢复奠定的坚实基础,2017年实施的山水林田湖草项目工程也更推动了祁连山地区的植被恢复。
在空间尺度上,祁连山东南地区多为中高覆盖度,低覆盖度多分布在西北方向地区。这是由于祁连山西北方向地区海拔较高,随着海拔的逐渐升高,气温变化幅度较大,制约了植被生长发育。东南地区多为低海拔、中海拔,适宜的温度以及丰沛的降水促进了植物的生长,又因退耕还林和天保工程项目的实施,减少了人类活动对植被覆盖的反向抑制作用,进而使得祁连山东南地区植被增加。祁连山地区植被覆盖结构变化表现为低植被覆盖度面积逐渐减少,中低覆盖度向中高覆盖度方向转化,低植被覆盖面积减少了1.26×106km2,面积占比减少了6.52%。其中,刚察县、祁连县和天峻县的中低覆盖度向中高覆盖度面积转移最多,其原因为地处祁连山中段腹地,环绕青海湖,且境内有多条河流通过,优越的地形、气候条件为植被生长提供了良好的生长环境,同时县内人民积极响应国家号召,完成森林生态效益补偿以及农牧民的搬迁工作,减少对生态环境的干扰,该项目工程的实施促使祁连山区森林资源进一步得到了有效保护,为保障祁连山可持续发展提供了强有力的生态支撑。
4.2 祁连山植被覆盖变化驱动作用
1986—2020年间,祁连山地区干旱指数、降水和气温均呈现增加趋势,而相对湿度则为减少趋势,表明区域气候逐渐呈现暖干化趋势,由相关性分析结果发现,研究时段内植被整体增长趋势与气候的变化趋势有密切关系。
本文研究结果指出影响祁连山植被变化的主要气候因素为气温和降水,该结论与大量研究结果相一致[49—53],但本文发现祁连山地区整体植被覆盖变化与干旱指数表现为显著正相关(P<0.01),与相对湿度为负相关,说明干旱也是制约祁连山地区植被生长的重要因素。依据本文划分的时间段,祁连山植被覆盖度同干旱指数、降水、气温的相关性方向以及程度均有不同,其中气温对植被覆盖的影响程度表现为增加的趋势,干旱指数、降水量则为减少趋势,造成这种现象的原因为祁连山地区冰川分布广泛,冰川融水流出形成祁连山水系,水资源丰富,区域植被对降水的依赖相对较小,而在全球变暖的大背景下,气温表现为逐渐增加的趋势,同时,植被生理作用对气温的干扰更为敏感,过高或者过低的温度会限制植被的呼吸作用、光合作用,从而影响植被的正常生长[54—55]。
本文研究结果指出人类活动对祁连山地区植被恢复具有重要意义。相关研究也指出人类活动对青藏高原植被恢复和黄河流域植被覆盖增加起到促进作用[56—58],文中残差分析的结果显示人类活动对区域植被覆盖情况主要产生正面促进作用,且产生正面影响的趋势正在逐渐增大,说明在实施相关生态修复工程的作用下,研究区植被覆盖有了明显改善,特别是实施的退耕还林还草、封山育林、天保工程、森林生态效益补偿等生态修复工程,一方面为植被的生长发育筑牢了良好的基础,另一方面也为植被的自然恢复提供有力的生态支撑。生态修复工程实施前,林区植被受到人类活动干扰较大,主要表现在人工伐木伐林以及周边村庄牲畜践踏,导致灌木林受到严重破坏,部分土地因植被破坏出现沙漠化趋势,而在实施生态修复工程后,森林资源、野生动植物得到严格的保护使得林业用地面积不断增加,表现出林木生长良好,植被恢复迅速等生态环境逐渐好转的现象,这也充分肯定了生态修复工程对于植被覆盖变化的正向促进作用。
研究时段内气候因素对植被覆盖的贡献率为93.40%,人类活动贡献率为6.60%,说明气候因素是影响祁连山区植被变化的主导因素,同时,受祁连山地区特殊的地形影响,气温、降水、干旱指数、相对湿度对植被覆盖的影响具有明显区域特征,但本文研究指出人类活动对植被覆盖变化的影响程度逐渐增加,且增加倾向为1.52/10a,说明人类活动在小尺度短时间范围内改变植被覆盖状况。
经过34年的生态恢复,祁连山地区植被覆盖面积增加,生态环境恢复良好。在探究气候因素对植被覆盖影响作用的同时,也应该肯定正确的人类活动对植被覆盖度的贡献。生态修复工程的实施,改变了区域植被覆盖低的现象,促进了区域植被覆盖的增长,对祁连山地区的植被覆盖增加产生了积极影响。因此,研究该区植被覆盖特征并量化环境驱动因子的影响,能够针对不同环境背景条件下合理优化植被格局,对于维系该区植被建设成果,实现区域生态保护同资源利用协调发展具有重要意义。
5 结论
本文基于Landsat NDVI数据集利用Google Earth Engine云平台研究了祁连山生态修复前后的植被覆盖的时空动态变化,并采用多元回归残差分析计算气候因素及人类活动对植被覆盖变化的贡献率,主要结论如下:
(1)1986—2020年祁连山地区植被覆盖表现为波动增加的趋势,平均上升幅度为0.022/10a,多年平均植被覆盖度为37.36%,这也表明经过34年的生态恢复,祁连山的植被覆盖度状况明显改善。
(2)祁连山地区植被覆盖空间上表现为东南地区以高植被覆盖度为主,西北地区多为低植被覆盖。实行生态修复工程后,祁连山植被覆盖多为低等级覆盖向高等级覆盖方向转移,其中,中低植被覆盖向中高植被覆盖转化面积最大。
(3)祁连山植被覆盖的变化受气候因素和人类活动共同作用。其中,多元回归残差分析结果显示气候因素对于植被覆盖的作用程度表现为下降趋势,而人类活动的影响则为增加趋势,其中生态修复工程对植被覆盖变化具有重要作用。