本研究使用的NDVI数据部分下载于国家青藏高原科学数据中心(https://data.tpdc.ac.cn/zh-hans/), 空间分辨率为30 m。数据集包含1990、1995、2000、2005、2010、2015、2019年共7年84期月度NDVI影像。为了有效减少大气中云、气溶胶、云阴影、视角及太阳高度角的影响, 数据基于最大值法(Max Value Composition, MVC)合成。数据集已经过辐射校正、大气校正等预处理, 后期使用ArcMap工具对数据集进行裁剪、镶嵌等处理, 通过软件计算将流域植被覆盖较好的4— 10月NDVI数据进行平均, 得到黑河流域年度NDVI数据集用于研究分析。

本研究使用的气温数据下载于中国气象数据网(http://data.cma.cn/data), 从该网站下载黑河流域10个气象台站的日值气温数据, 其中上游台站5个(阿柔站、垭口站、大沙龙站、景阳岭站、司大隆站), 中游台站3个(大满超级站、花寨子站、张掖湿地站), 下游台站2个(混合林站、荒漠站)。对数据进行整理, 挑选所需数据, 本研究使用的大于0 ℃的年积温是对站点全年大于0 ℃的日气温进行累加的结果。

2.2.1 植被覆盖度计算

植被覆盖度(Fractional Vegetation Cover, FVC)是指植被的叶、茎、枝在地面的垂直投影面积占单位面积的百分比^[33], 与NDVI之间存在极强的线性关系^[34]。像元二分模型计算植被覆盖度的计算公式为

f_vc= NDVI-NDVIsoilNDVIveg+NDVIsoil。 (1)

式中: f_vc为植被覆盖度; NDVI为像元的归一化植被指数值; NDVI_soil为无植被覆盖像元的NDVI值; NDVI_veg为完全由植被覆盖的像元的NDVI值。

像元二分模型的主要思想为: 假设像元光谱由纯植物光谱与纯裸土光谱线性混合而成, 将植被覆盖度定义为任一像元的NDVI在纯植物像元植被指数(NDVI_veg)与裸土植被指数(NDVI_soil)之间的线性插值。因此, 计算植被覆盖度的关键是NDVI_soil和NDVI_veg的取值^{[35, 36, 37]}。结合黑河流域植被覆盖的实际情况和前人研究经验^{[38, 39, 40]}, 经过反复试验对比, 本研究选定土地利用单元内累积频率为0.5%的值为NDVI_soil, 土地利用单元内累积频率为99.5%的值为NDVI_veg, 代入式(1), 通过 ArcGIS 10.7运算即得到研究区的历年植被覆盖度。

为了详细比较不同时间和空间位置上的植被覆盖情况, 本研究根据FVC值对研究区内的植被覆盖情况进行了分级。依据黑河流域的植被结构特点并结合前人的相关研究^{[40, 41]}, 经多方案试验并参考野外考查结果, 将植被覆盖情况分为5级(表1)。

表1

Tab.1

表1(Tab.1)

表1 植被覆盖度等级划分 Tab.1 Classification of vegetation coverage 植被覆盖度 等级划分 无植被覆盖度 0.00< FVC≤ 0.10 低植被覆盖度 0.10< FVC≤ 0.20 中植被覆盖度 0.20< FVC≤ 0.35 中高植被覆盖度 0.35< FVC≤ 0.50 高植被覆盖度 0.50< FVC

表1 植被覆盖度等级划分 Tab.1 Classification of vegetation coverage

2.2.2 植被覆盖度变化趋势计算

植被覆盖度变化趋势的估算通常使用一元线性回归分析方法^[42], 该方法不仅可以满足在若干时间节点的长时段内模拟变化趋势的要求, 还可以为认识空间演变规律提供指示信息。本文采用该方法反映植被覆盖度的变化规律, 计算公式为

θslope=n∑ni=1(ifvci)-∑ni=1i∑ni=1fvcin∑ni=1i2-(∑ni=1i)2。 (2)

式中: n为研究时段的年数; fvci为第 i年的植被覆盖度; θslope是植被覆盖度在研究年段期间的变化斜率, 其正负可以指示植被覆盖度的变化趋势。当 θslope> 0时, 代表植被覆盖情况改善, 呈增加趋势; 当 θslope=0时, 代表植被覆盖度没有变化, 保持稳定; 当 θslope< 0时, 代表植被覆盖情况退化, 呈减少趋势。此外, θslope值越大代表植被覆盖度增加的速度越快, 即植被覆盖情况改善越迅速, 反之则退化越迅速。因此, 通过 θslope值的正负及大小, 可以判断一个地区地表植被的活动情况。

依据自然地理常识及前人研究^[43], 选取黑河流域每年植被生长状况较好的生长季(4— 10月)的植被覆盖数据进行研究。根据研究年份生长季的植被覆盖度数据制成研究区不同时期的植被覆盖度等级空间分布图(图2)。黑河流域的植被覆盖度表现出明显的空间差异, 流域内以无植被覆盖区为主, 高、中高、中和低植被覆盖区面积相对较小。在空间分布上, 黑河流域的植被覆盖度呈现由南向北递减的趋势, 高、中高植被覆盖主要分布在上游山区, 中游平原地区主要为中高、低植被覆盖, 下游除额济纳旗和部分地区为中、低植被覆盖外, 其他区域主要为低植被覆盖或无植被覆盖, 地貌景观主要为荒漠。植被覆盖情况出现区域差异是黑河流域水资源空间分布格局造成的: 上游祁连山区降水量大, 适宜的温度促进了冰川消融, 水资源丰富; 下游额济纳旗地区是荒漠草原区, 降水量小, 气温高, 蒸发量大, 水资源匮乏。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 黑河流域不同时期植被覆盖度等级空间分布Fig.2 Classification of vegetation coverage in Heihe River Basin during different periods

根据年植被覆盖度数据进行面积统计, 制成研究区不同时期各等级植被覆盖度面积统计表(表2)。由表2可知, 1990— 2019年, 黑河流域有植被覆盖区域的面积呈微弱上升的趋势; 1990— 2000年间, 无植被覆盖区面积基本维持稳定; 2000年之后, 无植被覆盖区面积逐渐减少, 无植被覆盖区可能向其他的植被覆盖类型转化。1990— 2000年, 低植被、中植被覆盖区面积略有减少, 从2005年开始, 低植被覆盖区面积逐年上升, 而中植被覆盖区面积整体呈波动上升趋势, 指示生态环境在逐渐恢复, 恢复区主要分布于黑河流域下游, 可能得益于2000年前后开始实施的黑河分水计划。中高植被覆盖区面积自1990年开始一直呈上升趋势, 大部分由无植被覆盖区转换而来, 高植被覆盖区面积30 a间变化不大, 无明显变化趋势。

表2

Tab.2

表2(Tab.2)

表2 黑河流域不同时期各等级植被覆盖度面积统计表 Tab.2 Area statistics of each vegetation coverage in Heihe River Basin during different periods 年份 各等级植被覆盖度面积/万km2 无植被 覆盖 低植被 覆盖 中植被 覆盖 中高植 被覆盖 高植被 覆盖 1990 10.67 1.77 0.72 0.54 0.60 1995 10.87 1.53 0.73 0.60 0.57 2000 10.64 1.58 0.71 0.76 0.60 2005 10.32 1.59 0.93 0.84 0.62 2010 10.12 1.70 0.99 0.92 0.58 2015 10.15 1.77 0.74 1.13 0.50 2019 10.07 1.78 0.71 1.14 0.61

表2 黑河流域不同时期各等级植被覆盖度面积统计表 Tab.2 Area statistics of each vegetation coverage in Heihe River Basin during different periods

将1990— 2019年植被覆盖度影像数据代入式(2), 在ArcMap 10.7 软件中进行制图, 得到黑河流域近30 a以来的植被覆盖变化程度空间分布(图3)。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 1990— 2019年黑河流域植被覆盖变化程度空间分布Fig.3 Distribution of vegetation coverage changes in Heihe River Basin from 1990 to 2019

结果显示: 1990— 2019年流域大部分区域的植被覆盖度维持不变或增加, 仅少数地区的植被覆盖度有所减少(表3)。植被覆盖度增加区域的面积远大于退化区域。植被覆盖度极显著增加的区域主要分布在祁连山东南侧、中游城市周边及靠近黑河两侧的区域; 植被覆盖度显著减少或极显著减少区域的面积较小, 主要分布于黑河下游和中游的部分地区, 上游也有零星出现, 主要对应区域为中游张掖市山丹县、甘州区、临泽县和高台县, 以及下游部分荒漠地区。

表3

Tab.3

表3(Tab.3)

表3 1990— 2019年黑河流域植被覆盖变化程度统计表 Tab.3 Statistics of vegetation coverage changes in Heihe River Basin from 1990 to 2019 植被覆盖变化程度 面积/万km2 比例/% 极显著退化 0.07 0.51 显著退化 0.03 0.23 基本不变 11.10 77.61 显著增加 2.26 15.81 极显著增加 0.83 5.84

表3 1990— 2019年黑河流域植被覆盖变化程度统计表 Tab.3 Statistics of vegetation coverage changes in Heihe River Basin from 1990 to 2019

大于0 ℃的年积温是影响植被生长的重要因素, 会极大地影响植被生长的适宜性。本文对1990— 2019年黑河流域大于0 ℃年积温的变化进行了趋势分析(图4), 结果显示流域平均大于0 ℃年积温在下游最大, 约4 650 ℃; 中游次之, 约3 000 ℃; 上游最小, 约1 340 ℃。1990— 1995年, 流域大于0 ℃年积温有一个较大幅度的升高, 自2000年至今, 年积温呈波动上升趋势, 2005年和2015年的年积温最低, 中游增加幅度最大, 下游次之, 上游最小。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 黑河流域大于0 ℃年积温变化趋势图Fig.4 Trend chart of the 0 ℃ annual accumulated temperature change in Heihe River Basin

气温是控制地球表面一切生物物理化学过程的主要因子, 植被生长需要适宜的温度范围, 温度过高或过低都会阻碍植被生长。此外, 气温的升高或降低会导致流域内其他自然要素的变化, 进而影响植被的生长。温度升高会导致冰川消融、积雪融化, 使流域径流和降水量增大, 改变流域水资源供给, 进而影响植被的生长变化。温度过高会提高蒸散发, 增加植被对水资源的需求, 使土壤水分减少, 加剧土地沙漠化, 导致植被覆盖退化。相反, 当气温低于植被生长的需要则会使植被生长受阻^[44]。

研究区现有的水资源分布格局为从南部祁连山区到北部内蒙古额济纳旗逐渐减少, 这是由研究区降水、地形及海拔梯度共同决定的^{[45, 46, 47]}。水资源分布格局在一定程度上影响了该地区植被的分布区域, 使地面植被覆盖从南向北形成了由温带山地森林草原到荒漠植被的分布规律(图2)。

在全球气温持续升高的大背景下, 我国平均气温的增加速度明显高于全球均温增速, 其中以我国北部和西部地区尤为突出^[48]。气温变化会引起降水和径流的改变, 从而影响整个流域的地面植被分布。黑河的径流量主要依靠祁连山区的冰川融水补充, 近30 a来, 黑河流域的气温在波动中趋于上升(图4), 下游大于0 ℃ 年积温增长为103.95 ℃, 中游年积温增长最多, 为150.10 ℃, 上游增长相对较少, 但也达到了55.04 ℃。温度升高会导致上游地区降水增多, 冰川积雪融化加快, 黑河径流量增加, 进而对地下水的侧向补给加强, 地下水水位抬升, 缓解河流两侧植被的需水压力, 使流域上游及两侧的植被覆盖度明显增加。中游气温上升提高了植被需水量, 导致到达下游的水量减少, 加之下游分布了广阔的荒漠地区, 气温上升加剧了蒸散发, 因此下游植被退化最为严重(图3)。