基于MODIS-NDVI 的岷江流域植被时空演变及地形分异研究

王鑫, 薛飞阳, 苏子昕, 青玲萱, 杨存建*

1. 四川师范大学西南土地评价与监测教育部重点实验室, 四川 成都 610068；

2. 四川师范大学地理与资源科学学院, 四川 成都 610068

植被作为地球表面的重要覆盖物之一，同时也是生态系统中不可缺少的一部分，能够通过与水体、大气、土壤等物质的相互作用，起到有效调节该地区生态环境的作用[1]。植被的空间信息分析与研究对于生态环境的建设与保护起到了重要的作用。遥感技术为植被信息获取提供了新的方法与手段[2]。目前随着遥感技术的迅速发展，时间序列MODIS影像，以较高的时间分辨率和优异的多光谱特性，其植被数据已广泛用于全球及区域尺度植被变化研究[3]。而植被指数能够较好地反映地表植被覆盖状况和生长情况，比值植被指数(RVI)、归一化植被指数(NDVI)、差值植被指数(DVI)是较为常见的3 种植被指数形式[4]。其中归一化植被指数（NDVI）应用最为广泛，其原理是通过计算植被反射的近红外波段（NIR）与植被吸收的红光波段（RED）之间的差异对植被进行量化。目前，国内外有许多专家学者利用MODIS-NDVI 数据对不同区域的植被演变特征进行了深入研究。如杨城等[5]基于MODIS-NDVI 数据探讨了青藏铁路沿线的植被演变特征及其响应因子，得出植被演变受降水以及人类活动的影响程度较大。王思等[6]研究了广东省的植被覆盖变化特征，得出不同土地利用类型对该省植被覆盖度的贡献各异的结论。易扬等[7]探究了长江中游地区的植被时空变化特征及其气候响应因子。

岷江作为长江的重要分支之一，其流域面积内地形复杂，生态环境脆弱[9]。该区域内植被作为生态环境的重要调节者，前人对其流域内部分地区进行植被演变研究。但缺乏一定的时效性以及全面性，本文基于MODIS-NDVI 时序数据和地形数据，并借助空间分析方法、变异系数、Hurst 指数、线性趋势分析等方法对岷江流域2003 年至2021 年的植被时空演变特征、未来趋势变化及其地形影响因子进行研究。

1 研究区概况、数据来源与研究方法

1.1 研究区概况

岷江发源于四川省岷山南麓，是长江上游的重要支流。介于99°38′-105°51′E，28°16′-33°39′N 之间，干流全长约711 km，总落差达到3560 m，其流域面积135 881km2。岷江流域内地形差异较大，以都江堰－大邑－洪雅－峨眉－线为界，此线以西，多高山峡谷，地形复杂。其中支流杂谷脑河、黑水河沿河狭长河谷地带，类型主要为干温河谷[10]。其东部为丘陵与成都平原，城镇林立，人口众多。流域内森林资源丰富，主要粮食作物为水稻、小麦、玉米、油菜等（见图1）。

图1 研究区位置图Fig. 1 Location of the study area

1.2 数据来源及预处理

MODIS-NDVI 数据是通过遥感云计算平台(GEE)获取的岷江流域的数据集，时间尺度为2003 年至2021 年，分别包括MOD13Q1 数据与MYD13Q1 数据。空间分辨率为250 m，通过使用这两种数据产品，而将原本的时间分辨率为16 d 提升至8 d。首先，将获取得到的MODIS-NDVI 数据使用ArcMap 软件批量裁剪出岷江流域的NDVI 数据；其次，使用最大合成法（MVC）计算得到月NDVI 数据，用于有效剔除遥感影像中的云影等误差数据，再由月NDVI数据合成年均值NDVI 数据，最后，得到岷江流域逐月及逐年NDVI 值。

地形数据采用NASA 发布的ASTER GDEM V2版数据，空间分辨率为30 m。为保证提取结果的准确性，将其重采样至250 m 空间分辨率。将提取得到的地形因子与植被变化趋势进行结合分析，研究地形因子与植被变化的相关关系。

1.3 研究方法

1.3.1 植被覆盖度稳定性

变异系数是一种描述概率分布离散程度的归一化量度，通过变异系数可以很好地反映植被覆盖度的波动程度[11]，定义为原始数据标准差与原始数据平均数的比，其计算式为：

式中：C 为变异系数，S 为标准差，V 为研究区多年植被覆盖度平均值，n 为研究年限长度， Vi为第i 年植被覆盖度。

1.3.2 趋势分析

采用一元线性回归分析方法，逐像元分析研究区内年均NDVI 的变化趋势，其计算式为：

式中：S 为回归方程的斜率，i 为年份，n 为研究总时间跨度， Vi为第i 年的年均NDVI。当S 小于0时,表示NDVI 处于下降趋势，当S 大于0 时，表示NDVI 处于增长趋势。采用F 检验法对年均NDVI 的变化趋势进行显著性检验，根据检验结果将变化趋势分为6 个等级。

1.3.3 Hurst 指数

基于重标极差法(R/S)的Hurst 指数可用于定量描述时间序列长期依赖性[7]。其计算式为：

对于时间序列{NDVI( t )}，t = 1，2，…，n，定义均值序列：

式中，Hurst 指数即H 值可根据log(R/S)n=a+H×log(n)利用最小二乘法拟合得到。根据H 值判断NDVI 序列是随机变化或存在持续性。H 值的取值范围为0～1，若0.5≤H≤1，表明该时间序列具有可持续性，未来趋势变化与过去时间段内的趋势相同；若H=0.5，表明该时间序列具有随机性，未来趋势变化与过去时间段内的趋势无关；若H≤0.5，表明未来趋势变化可能保持基本稳定或者发生逆转，在研究中将此类归为变化趋势不确定[7]。

2 结果与分析

2.1 岷江流域植被的时间演变特征

通过2003～2021 年的岷江流域植被NDVI 平均值（见图2），对岷江流域的植被覆盖随年际变化的特征进行分析。可以得到，2003～2021 年岷江流域的植被NDVI 均值在0.55—0.60 范围内，最低值出现在2012 年的0.55，最高值为2019 年的0.60。根据一元线性方法所拟合的线性函数可知，近年来岷江流域的植被NDVI 均值整体呈波动上升趋势，增长速率为2.4%·10a-1，表明岷江流域植被覆盖度在稳步增加。

图2 岷江流域植被NDVI 年际变化Fig. 2 Interannual changes of vegetation NDVI in Minjiang River basin

2.2 岷江流域植被的空间演变特征

2.2.1 NDVI 空间分布基本特征

依据2003～2021 年岷江流域年均NDVI 数据，计算19a 平均值得到平均NDVI 的空间分布图(见图3)。可以看出，岷江流域植被分化较为明显，NDVI 值较高的区域主要分布在川西高原的河谷、雅安市、乐山市南部、凉山彝族自治州北部。其次主要分布在果洛藏族自治州南部、德阳市、眉山市、内江市、自贡市、资阳市、宜宾市、重庆市、成都市非中心城区域、乐山市北部。而NDVI 值较低的区域主要分布在川西高原的高海拔山地以及城市城镇分布区，根据影像判断该区域主要为积雪以及建筑物。

图3 岷江流域NDVI 空间分布Fig. 3 Spatial distribution of NDVI in Minjiang River basin

2.2.2 NDVI 空间波动变化特征

根据岷江流域2003—2021 年NDVI 结果，逐像元统计变异系数值，依照计算结果将稳定性分为五个程度[12]，见表1。由空间变化波动性图（见图4）可以看出，在该研究时段内岷江流域大部分区域的植被波动性偏低，其中处于较低波动变化的区域占比最高，达到43.83%。主要分布在雅安市、乐山市、宜宾市、自贡市、内江市、眉山市西部、成都市西部。其次是低波动变化区域占比达到43.56%，主要分布在川西高原地区的山谷、成都市东部、资阳市、眉山市东部、德阳市中部。中等波动变化区域占比为7.89%，主要分布在成都西北缘、乐山市南部。较高波动变化以及高波动变化区域分别占比为2.21%和2.39%，主要分布在川西高原的高海拔山地、成都市城镇区域。

表1 植被NDVI 波动分类及比例Tab. 1 Classification and proportion of vegetation NDVI fluctuation

图4 岷江流域植被NDVI 空间变化波动性Fig. 4 Spatial fluctuation of vegetation NDVI in Minjiang River basin

2.2.3 NDVI 的趋势分析

（1）NDVI变化趋势

采用一元线性回归趋势分析逐像元模拟影像每个栅格该研究时段内的变化趋势（见图5），并对分析结果进行显著性检验，把结果划分为6 个变化等级[20]（见图6），分别为极显著增加、显著增加、不显著增加、极显著减少、显著减少、不显著减少（见表2）。可以看出，岷江流域植被NDVI 呈增长趋势的面积占比为84.49%，而岷江流域面积的15.51%呈下降趋势。在处于增长趋势区域中，极显著增加所占比值最大，主要分布在川西高原的河谷地区、雅安市东部、成都市西北缘以及龙泉山、眉山市西部、乐山市、宜宾市、自贡市、内江市。而下降区域主要分布在川西高原的高海拔山地、成都市向外扩张的新城区、眉山市城区，主要由于城市的快速发展以及城镇扩张等因素所造成[13,14,18]。

表2 植被NDVI 变化趋势显著性统计结果Tab. 2 Statistical results of the significance of NDVI trends

图5 岷江流域植被变化趋势Fig. 5 Vegetation trends in Minjiang River basin

图6 岷江流域植被变化趋势显著性检验结果Fig. 6 Significance test results of vegetation change trend in Minjiang River basin

（2）NDVI未来趋势分析

岷江流域在该研究时间内，植被Hurst 指数小于0.5 的区域占比为57.04%,大于0.5 的区域占比为42.96%。根据Hurst 指数与植被变化趋势的耦合图（见图7），可以看出，植被持续性增加的面积占比为36.05%，主要分布在川西高原的部分河谷地区、成都市西北缘、雅安市中部、乐山市西部、内江市东部、自贡市。植被持续性减少的面积占比为6.91%，主要分布在成都市龙门山附近以及扩张的新城区、眉山市城区。未来趋势不确定所占比值达到最高，为57.04%，主要分布在川西高原。

图7 岷江流域植被未来趋势Fig. 7 Future trends of vegetation in Minjiang River basin

2.3 岷江流域植被的地形效应

因岷江流域地形分异明显，为进一步探究植被变化与地形的关系，将所提取的高程、坡向和坡度进行重分类，探究在不同等级的地形因子上植被的变化情况。其中将高程以1 000 m 为间隔划分为＜1 000 m、1 000～2 000 m、2 000～3 000 m、3 000～4 000 m、＞4 000 m。将坡向划分为阴坡、半阴坡、阳坡、半阳坡和平地[15]。将坡度划分为＜5°，5°～15°，15°～25°，25°～35°，≥35°[16]。

2.3.1 高程

在不同的海拔区域内，植被的变化趋势分布也明显不同。如图8 所示，极显著增加主要发生在1 000 m 以下的海拔范围内，占比达到该流域面积的19%。而不显著增加面积占比达到14%，主要发生在4 000 m 以上的海拔范围内。极显著减少也主要发生在1 000 m 以下的范围内，但其在占比较少，仅占总面积的1%。在每个高程等级内，植被增加的占比均大于植被减少的占比，研究区域内的植被增加的趋势明显。

图8 岷江流域各海拔等级的植被覆盖面积变化Fig. 8 Changes of vegetation cover area at different altitudes in Minjiang River basin

2.3.2 坡向

因平地各植被变化趋势总占比不足1%，将平地不纳入研究范围内，这里比较各个坡向的植被变化趋势。如图9 所示，极显著增加的面积占比在阳坡区域内相对较多，植被减少发生在半阴坡的占比相对较多，但各个坡向的植被变化趋势相差不大，因此坡向因子不作为植被变化主要影响因素[19]。

图9 岷江流域各坡向等级的植被覆盖面积变化Fig. 9 Changes in vegetation cover area of different slope grades in Minjiang River basin

2.3.3 坡度

如图10 所示，植被覆盖的变化主要发生在坡度小于35°的范围内。随着坡度的增加，植被的极显著增加变化趋势的占比逐渐减少。植被的减少变化趋势主要分布在5°～25°的坡度范围内。在各个坡度等级上，植被的增加趋势占比均大于减少趋势的占比。

图10 岷江流域各坡度等级的植被覆盖面积变化Fig. 10 Changes of vegetation cover area of different slope grades in Minjiang River basin

3 讨论

自2003～2021 年以来，岷江流域的植被覆盖整体呈现波动增长趋势，这与李婷[17]等研究结果基本一致。岷江流域是四川省退耕还林还草重大生态工程实施的重点区域，并在天保一期以及天保二期工程的推进下，实现了植被覆盖的稳步增长。其中岷江流域植被改善区域主要分布在川西高原的河谷部分以及盆地内部非城镇区域。而成都市向外扩张的新城区、眉山市和德阳市因为受近年来人口增长、城镇扩张等因素的影响出现了明显的植被退化趋势[18]。

岷江流域的地形分异明显，植被的分布以及变化趋势与其地理环境息息相关,而地形因子又是决定热量、水文、土壤等的主要因素。其中海拔和坡度作为影响岷江流域植被分布以及变化的主要地形因子，这一结论与张诗羽等[19]人的研究结果一致。其中海拔处于1 000 m 以下以及坡度小于25°的地区，植被极显著增加趋势最明显，这得益于各种植被保护措施[20]，后期也应该将保持或进一步完善各项生态修复工程。在未来岷江流域的发展过程中，城镇扩张的同时应注意植被的变化趋势以及生态环境的修复。

4 结论

基于岷江流域2003～2021 年MODIS-NDVI 数据集和DEM 数据，分析该研究时段内岷江流域植被NDVI 的时空变化特征，并探究不同地形因子与植被变化的相关性。得出以下结论：

（1）从植被NDVI 空间分布特征来看，研究时段内的NDVI 均值介于0～0.9 之间。其流域的NDVI值上中游偏高，下游偏低，空间上分布呈现出明显的差异性。从时间变化特征上看，岷江流域植被整体呈现波动上升的趋势，增长速率为2.4%/10a，其中最低值出现在2012 年的0.55，最高值为2019 年的0.60。

（2）从变化趋势来看，岷江流域植被覆盖呈增加趋势和减少趋势的面积分别占84.49%和15.51%。在增加趋势中，极显著和显著增加的面积占48.31%，主要分布在川西高原的河谷地区以及四川盆地的非城镇区域。在减少趋势中，极显著和显著减少的面积占2.38%，主要分布在成都市向外扩张的新城区、眉山市和德阳市。

（3）从稳定性来看，波动性较强的区域主要分布在川西高原的高海拔山地以及成都市，而低波动主要出现在川西高原的相对低海拔区域，植被较为稳定。

（4）从未来趋势来看，岷江流域的大面积区域的未来变化趋势不确定，植被变化呈现持续性减少的主要是成都市向外扩张的新城区。而呈现持续性增加的主要分布在四川盆地的非城镇地区以及川西高原的河谷地区。

（5）从地形影响因子来看，海拔和坡度作为影响植被变化的主要因素，其中低海拔（1 000 m 以下）以及低坡度（5°以下）的植被的极显著增加趋势明显，坡向与植被变化的相关性并未显示出明显的规律。