基于时间序列的新疆某流域植被覆盖度动态变化研究

摘要 以新疆某流域为研究对象，基于MODIS NDVI数据，结合DEM引入地形位指数，分析研究区2011—2020年植被覆盖度时空变化特征。结果表明，研究区植被覆盖度变化在这10年间趋于改善，且以轻微改善为主;各年度植被盖度均值由0.429 3增长至0.529 2;就地形因素而言，植被盖度随地形复杂程度的增加呈现先增后减的变化趋势，当地势复杂程度高于某一阈值时，植被盖度出现下降。

关键词 植被覆盖度;NDVI;地形位指数;时空变化特征

中图分类号 Q 948  文献标识码 A

文章编号 0517-6611（2022）16-0050-03

doi：10.3969/j.issn.0517-6611.2022.16.014

开放科学（资源服务）标识码（OSID）：

Study on Dynamic Change of Vegetation Coverage in a River Basin of Xinjiang Based on Time Series

LI Xiao-ran1，2，3，4， REN Xuan1，2，3，4， TANG Meng-ying1，2，3，4 et al

（1.Xinjiang Academy of Environmental Protection Science， Urumqi， Xinjiang 830011;2.Xinjiang Key Laboratory for Environmental Pollution Monitoring and Risk Warning， Urumqi， Xinjiang 830011;3.Xinjiang Engineering Technology Research Center for Cleaner Production，Urumqi， Xinjiang 830011;4.Junggar Desertoasis Ecotone Station for Scientific Observation and Research of National Environmental Protection， Urumqi， Xinjiang 830011）

Abstract Taking a river basin in Xinjiang as the research object， based on MODIS NDVI data， combined with the topographic index introduced by DEM， the temporal and spatial variation characteristics of vegetation coverage in the study area from 2011 to 2020 were analyzed.The results showed that the change of vegetation coverage of the study area tends to improve in the past 10 years， and it was mainly realized as a slight improvement.The average of vegetation coverage increased from 0.429 3 to 0.529 2.In terms of topographical factors，the vegetation coverage increased first and then decreased with the increase of terrain complexity. When the terrain complexity was higher than a certain threshold， the vegetation coverage decreased.

Key words Vegetation coverage;NDVI;Terrain index;Temporal and spatial variation characteristics

植被是土壤和气候条件共同作用的产物，同时又能通过光合作用反作用于气候，对调节气候有重要意义[1-2]。陆地生态质量的变化必然会反映在植被类型、数量及质量的变化上，而植被覆盖度指数是目前既能表征植被的茂密程度，又代表植物进行光合作用面积较为准确的计量手段，是衡量生态系统质量变化的重要因子之一[3-6]。

近年来国内外学者对植被覆盖度的相关内容开展了大量的研究，如马琳雅等[7]对甘南草地类型基于MODIS数据对2001—2011年植被覆盖度进行动态变化分析，结果表明，甘南草地植被覆盖度水平升高;白庆坤等[8]从植被覆盖度的角度分析了盖度变化对蒙古高原干旱的响应;焦世文等[9-10]研究了NDVI的时空变化特征及其与环境因子的关系。

笔者基于GIS与遥感手段，以MODIS NDVI数据为数据源，分析了2011—2020年研究区天然植被覆盖度的时空变化特征，并引入地形位指数，综合考虑海拔与坡度对植被覆盖度的影响程度，以期为研究区生态质量评估提供科学依据。

1 资料与方法

1.1 研究区概况

研究区位于新疆维吾尔自治区最西端的天山山脉腹地，北部以科古琴为界，西以阿拉爱格尔山口与哈萨克斯坦共和国接壤，南部以天山山脉为界，东以那拉提山为分水岭、与开都河流域相隔。研究区为我国北方防沙带生态保护和修复重点工程布局的关键区域之一，是我国35个生物多样性保护优先区之一，同时也是我国西北区域最重要的生态安全屏障，对维护国家重要生态功能区具有重要意义。

1.2 数据来源及处理 该研究所使用的NDVI数据来自2015—2020年新疆生态状况变化遥感调查评估项目和生态环境部卫星环境应用中心，DEM数据来自地理空间数据云，其空间分辨率为30 m，生态系统类型数据来自中国科学院新疆生态与地理研究所。

1.3 研究方法

植被覆盖度的反演目前有效方法为像元二分模型[11]，该模型是通过遥感卫星传感器所获取的信息可以分为绿色植物所提供的信息和裸土所提供的信息2个部分，植被覆盖度计算公式如下：

FVC=（NDVI-NDVI s）/（NDVI-NDVI v）   （1）

式中，FVC为植被覆盖度，NDVI s和NDVI v分别为纯裸土和纯植被情况下NDVI像元值。

该研究以流域范围去除农田、城镇、湿地及永久冰川积雪外的自然生态系统植被覆盖度为研究区，以2011—2020年为研究时段，为削弱单年份气温、降水等因素对植物生长的影响，以2014—2016年各年度植物生长季植被覆盖度平均值作为2015年植被盖度水平值，以2018—2020年各年度植物生长季平均值作为2020年植被盖度水平值。根据前人研究基础，将植被覆盖度分为低覆盖度（0～0.2）、较低覆盖度（>0.2～0.4）、中覆盖度（>0.4～0.6）、较高覆盖度（>0.6～0.8）、高覆盖度（>0.8～1.0）5个等级[12]。分析2015—2020年植被覆盖度变化，将这5年植被覆盖度变化情况共分为9级，分级标准见表1。

同时，引入地形位指数[13]因子，从高程、坡度等地形因素综合分析植被覆盖度高低及变化特征。详见公式（2）。

T=lg（EE+1）×（SS+1）

（2）

式中，T为空间某个栅格单元的地形位指数，E、E分别为某个栅格单元的高程和研究区内的平均高程，S、S分别为某个栅格单元的坡度和研究区内的平均坡度。该指数越小，表明某区域高程与坡度均较小;数值越大，表明某区域高程与坡度均较大。

2 结果与分析

2.1 植被覆盖度年度特征

借助ArgGIS 10.2软件，计算得到研究区内2011—2020年各年度植被覆盖度最大值及均值。从图1可以看出，研究区植被覆盖度这10年间呈略微改善的趋势，10年间植被盖度最大值的最低值出现在2014年，为0.824 2;植被盖度最大值的最高值出现在2015年，为0.863 8。其中2013—2017年植被盖度最大值出现转折，呈现“降—增—降—增”的变化趋势，在2017年后趋于稳定。研究区内植被覆盖度均值在0.429 3～0.529 2，2014年均值较2013年变化率为-14.92%，2015年均值较2014年变化率为19.05%，其他各年度较前一年变化率均在正负8%以内，变化不大。

分年度统计研究区各等级植被盖度面积占比情况，结果如表2所示。从表2可以看出，研究区植被盖度以较高覆盖度和中覆盖度为主，且较高覆盖度除2014年外，其他各年份占比均接近50%;高覆盖度占比普遍较低，研究区植被高覆盖度区域分布较少，高覆盖度区域各年度面积占比均在5%以下。

2.2 植被覆盖度变化分析

根据研究区各年度植被覆盖度（图1），计算得到2015—2020年植被覆盖度变化情况（表3）。由表3可知，研究区内植被盖度以基本稳定为主并呈改善趋势，改善程度以轻微改善为主，占评价区总面积的24.66%;其次为中度改善，占研究区总面积的8.70%。轻微及中度改善区域主要分布在依连哈比尔尕山及哈尔克他乌山前山带，流域内平原区零散分布，显著改善区域集中分布在南部荒漠区域内。

退化区域面积较小，占研究区总面积的5.45%，且以轻度退化和中度退化为主。轻度与中度退化相交分布，主要集中在研究区西北部博罗科努山南脉前山带及中部平原区的农田绿洲外围。轻度和中度退化区域生态系统类型以荒漠和草地为主。研究区内存在零星极重度退化区域，主要分布在高山区，极重度退化区域生态系统类型以裸土、裸岩为主，初步考虑为稀疏草地退化为裸土、裸岩等生态系统。

植被覆盖度基本稳定区为退化与改善程度均在10%以内的区域，其面积占研究区总面积的56.95%，该区域分布较广，研究区前山带及谷地内平原区均有分布。

2.3 植被覆盖度变化的地形分异特征

地形位指数综合考虑坡度和海拔因子对植被覆盖度变化的影响，研究区地形位指数在0.095 98～1.344 43，基于ArcGIS空间分析功能，将研究区地形位指数通过自然间断点分级法分为5级，其中Ⅰ级（0.095 981～0.296 711）、Ⅱ级（>0.296 711～0.482 754）、Ⅲ级（>0.482 754～0.668 797）、Ⅳ级（>0.668 797～0.848 048）、Ⅴ级（>0.848 048～1.344 426）。不同地形梯度内植被覆盖度变化情况如表4～5所示。

研究区内不同地形梯度中，Ⅰ 梯度仅占评价区总面积的7.78%，该梯度内植被覆盖度变化以基本稳定为主;Ⅱ、Ⅲ、Ⅳ梯度分别占评价区总面积的22.06%、26.06%和26.92%，且植被覆盖度变化均以基本稳定为主;V梯度内植被覆盖度变化以改善为主，改善区域面积占研究区总面积的9.72%。基本稳定和改善区域在Ⅰ～Ⅳ梯度时随着地形位指数增加，分布面积逐渐增加，在V梯度时面积占比出现下降。退化区域在Ⅱ梯度占比最大、基本稳定区域在Ⅲ梯度占比最大、改善区域在Ⅳ梯度占比最大。

3 结论与讨论

该研究基于NDVI数据，采用像元二分模型模拟研究区植被覆盖度，为准确计算天然的植被覆盖度，将农田、城镇、湿地、冰川永久积雪等生态系统类型作为非评价区，对2011—2020年天然植被覆盖度的时空变化特征进行了研究。结果表明，2011—2020年研究区内盖度平均值由0.429 3增长到0.529 2，增长率为23.27%，呈现以较高覆盖度和中覆盖度为主的分布趋势，天然植被盖度逐步向较高覆盖度等级转化。研究时段内，植被覆盖度变化以基本稳定为主，部分区域发生不同程度的改善，改善程度主要以轻微改善和中度改善为主;盖度退化的以轻度退化为主。随着地形位指数的增加，植被覆盖度基本稳定和改善的面积占比均呈现先增后减的变化趋势，说明随着地势复杂程度的增加，植被盖度变化先出现好转，当地势复杂程度高于某一阈值时，植被盖度出现下降。