本文所采用的植被NDVI数据为MODIS MOD13Q1产品,来源于美国国家航空航天局(NASA),空间分辨率为250 m×250 m。该数据已经过水、云、气溶胶等的去除处理,首先对其进行镶嵌、格式和投影转换,然后采用最大值合成法合成年NDVI数据; 气象数据来源于中国气象数据网(http:∥data.cma.cn/),选取2001—2017年黑河上游及周边14个气象站点,利用Kring方法插值为空间分辨率均为1 km影像; 高程数据来源于地理空间数据云(http:∥www.gscloud.cn/)中的SRTM 90 m DEM产品,坡度、坡向数据则基于DEM数据使用ArcGIS软件生成; 中国植被图集和土壤类型空间分布数据来源于中国科学院资源环境科学数据中心(http:∥www.resdc.cn/)。

由图2可知,2001—2017年NDVI的年际变化呈波动增加趋势,年增长速率为0.001 7。NDVI多年均值为0.46,最大值出现在2010年为0.49,2001年为最小值0.42。多年NDVI变化在均值上下波动,向植被覆盖增加的方向变化。黑河上游地区植被低覆盖面积占比由2001年23.83%降低到2017年的14.66%; 中低覆盖多年面积平均占比在26.55%左右; 中等覆盖面积占比由2001年的18.21%增长到2017年的22.48%,增长相对较为明显; 中高植被覆盖区面积占比基本均大于30%,是研究区植被覆盖的主要类型; 高植被覆盖面积占比较小,处于3%左右。

2001年和2017年黑河上游植被均呈现东南高、西北低的空间分布特征(图3)。NDVI的高值区域主要集中在黑河上游东南部的八宝河流域,以高寒草甸、高寒草原、灌丛为主要的植被类型。NDVI的低值区域分布在西北部,植被类型多为半荒漠化草原。整体来看,黑河上游植被NDVI整体呈现改善的趋势,改善区域主要位于研究区中部和北部。

图2 黑河上游2001-2017年NDVI年际变化

图3 2001年和2017年黑河上游NDVI空间分布格局

2001—2017年黑河上游流域植被NDVI年均值变化率空间分布见图4,根据其变化率将植被变化分为严重退化、轻度退化、基本不变、轻度改善和高度改善5种情况,各类面积统计见表3。由图4和表3可知,近17 a黑河上游植被变化在空间上整体呈比较稳定的状态,覆盖度逐渐增加。其中37.16%的面积处于基本不变的状态,轻度改善和高度改善的面积之和远大于轻度退化、严重退化区,表明黑河上游植被覆盖改善趋势明显。

图4 植被NDVI变化率空间分布

表3 2001-2017年黑河上游NDVI变化统计

通过因子探测计算出来q值的大小可以说明该自然因子对植被生长解释力的强弱。探测结果见图5,各自然因子对植被NDVI的影响力排序为:日照时数(0.44)>高程(0.42)>降水(0.39)>相对湿度(0.38)>气温(0.22)>平均风速(0.21)>土壤类型(0.16)>植被类型(0.15)>地表温度(0.10)>坡度(0.07)>坡向(0.01)。日照时数、高程、相对湿度、降水是主要影响因子,解释力均超过35%,表明黑河上游植被发育与高程密切相关,同时水热条件也是重要影响因素,这与已有研究结论一致^[29]。气温、平均风速作为次级影响因子解释力均在20%以上; 土壤类型、植被类型、地表温度解释力超过10%,坡度、坡向的单个因素解释力很小。

图5 NDVI空间分布影响解释力指标

生态探测可以反映各因子对植被NDVI影响的统计差异性显著,探测结果见表4。从表中可以看出:日照时数与高程对植被NDVI空间分布的影响无显著性差异,其他因子对植被NDVI空间分布的影响具有显著性差异; 降水与相对湿度、高程对植被空间分布无显著性差异,与其他因子对植被空间分布存在显著差异; 气温与平均风速、土壤类型对植被空间分布无显著性差异,与其他因子对植被NDVI空间分布具有显著性差异; 植被类型与平均风速、地表温度、土壤类型对植被空间分布无显著影响,与其他因子存在显著的差异。

表4 生态探测结果

交互探测结果表明(表5),双因子交互作用对黑河上游植被NDVI空间分布的影响力大于单因子作用,因子之间的交互作用均呈现非线性增强或双因子协同增强作用。由表5可知,高程与除坡向外各因子之间交互效应增强了二者对植被NDVI空间分布的影响,日照时数与除坡度、坡向外各因子之间的交互作用增强了二者对植被NDVI空间分布的影响,即存在显著双因子增强关系,这进一步表明了高程与日照时数的主导影响作用。降水与地表温度、坡度和坡向交互作用显著增强了降水对植被NDVI的影响。坡向、坡度与各因子之间交互作用显著增强了坡度坡向对植被的影响,这说明坡度坡向虽然对植被的直接影响较小,但是与其他因子之间的叠加效应大大增强了单一因子对植被的影响。

表5 自然因子交互作用

利用地理探测器风险因子探测分析得到适合植被生长的各自然因子的范围(表6),不同自然因子的植被NDVI均值差异明显。统计检验表明,气象因子对植被NDVI均值的影响存在一个最适宜范围。随着气温、相对湿度、降水、平均风速、日照时数、地表温度的增加,分别在2.12～3.08℃,53%～55%,457.80～510.12 mm,2.30～2.38 m/s,2 571.65～2 726.74 h,7.14～8.03℃范围内最适合植被生长,NDVI均值分别达到0.572,0.617,0.661,0.620,0.650,0.576。在黑河上游区域,热量成为限制高海拔地区的主要因子,日照时数和地表温度的最佳范围为植被生长提供一个良好的热量条件。

表6 自然因子适应类型或范围

高程因子作为影响植被空间分布的第二大因子,在研究区内随高程的变化,植被的类型和生长具有明显的地带性,坡度坡向的不同会给植被生长带来明显的差异。统计检验表明,随着高程的增加NDVI均值呈现先增后减的趋势,结合已有研究^[25],高程为3 400 m时植被指数最高,而后植被指数减小。本研究得出高程在3 008～3 471 m范围内适合植被生长,NDVI均值达到0.601。黑河上游源区位于祁连山,降水随高程增加而增加,与研究结果保持一致^[30-31],随着降水量的增加,NDVI均值不断增大。水热条件的良好耦合关系对促进植被生长提供基础条件。坡向通过影响日照时数和接受的太阳辐射量从而影响坡面温度; 在北方地区背向也会影响不同坡向的水分含量,对坡面植被的生长造成一定的影响。坡度通过地表陡缓程度对植被生长产生影响。风险探测结果表示,坡度在<2°,坡向为北坡时最适合植被生长,NDVI均值分别达到0.509,0.506。坡度坡向虽在本研究中单一因子的影响力较小,但与其他因子的交互作用大大增强了单一因子的影响力。

植被、土壤类型的不同会引起NDVI均值的波动变化。黑河上游主要植被类型有针叶林、灌丛、草原、草甸等,其中针叶林NDVI均值最大达到0.640。植被生长的基本条件是水分的充裕情况,植被的空间分布情况往往受限于土壤含水量,而土壤水分又受到地形和太阳辐射的影响。土壤类型主要分为半淋溶土、钙层土、干旱土、漠土、初育土、半水成土、水成土、高山土、冰川雪被9类,其中半水成土最适合植被生长,为黑河上游植被提供生长需要的水分,NDVI均值达到0.732。

黑河上游地区处于生态环境脆弱区和气候变化的敏感区,是探究影响植被变化因素的理想区域。面对生态环境问题日益恶化,如土地荒漠化、草场退化和水资源短缺等,实施了一系列生态系统修复工程,植被覆被发生了显著变化。本文结果表明:2001—2017年黑河上游植被NDVI增加范围广、趋势明显,植被改善的区域占总面积的49.25%,空间上呈东南高西北低的分布特征。自生态保护方案实施以来,黑河上游地区植被覆盖得到较好的改善^[23-24]。

本研究中,日照时数、高程和降水是影响黑河上游NDVI变化主要自然因子,并且影响力表现为日照时数>高程>降水。由此说明:黑河上游地形复杂多变,光热条件可显著影响植被的生长^[32],日照时数的增加可以使植被最大有效化地进行光合作用^[33]; 黑河上游植被分布垂直地带性较为明显,从高到低依次是高山灌丛草甸、青海云杉林、针阔混交林、山地荒漠和山地草原。张立峰^[34]、李国亮^[32]等研究了黑河流域植被变化以及对地形因子的响应,海拔高度对植被的生长有明显的影响; 黑河上游水分主要来自山区降雨、冰川积雪融水,其中降雨占比超过90%,降水的空间分布在很大程度上决定了水分的空间分布情况,植被的生长受水分条件的制约,体现了植被随水而生的特点^[35]。水资源充沛的条件有利于植被的生长,以草本为主的类型则受到降水量的作用更加显著。

(1)黑河上游植被时空变化显著。2001—2017年植被NDVI年际变化呈波动增加趋势,植被由低覆盖向高覆盖转化,中高覆盖的面积占比最高; 空间分布呈东南高西北低,植被改善区域主要分布在中部和西北部地区。

(2)各因子对植被NDVI的影响程度:日照时数>高程>降水>相对湿度>气温>平均风速>土壤类型>植被类型>地表温度>坡度>坡向。其中,日照时数和高程成为影响植被变化的主要自然因子且解释力超过40%; 其次是降水和相对湿度指数解释力均在30%以上。

(3)植被NDVI空间分布情况对于自然因子的响应具有明显的差异,除日照时数与高程,降水与相对湿度、高程,气温与平均风速、土壤类型,植被类型与平均风速、地表温度、土壤类型外,其余各因子之间对植被NDVI变化影响均具有显著的差异。

(4)黑河上游地区植被生长受到自然因子之间的交互影响,因子之间均呈现相互增强和非线性增强的关系。

(5)风险因子探测分析发现促进植被生长的各主要因子最适宜范围:半水成土,针叶林植被类型,气温为2.12～3.08℃,降水为457～510 mm,日照时数为2 571.65～2 726.74 h,平均风速为2.30～2.38 m/s,相对湿度为53%～55%,地表温度为7.14～8.03℃范围内,高程为3 008～3 471 m,坡度<2°,坡向为北坡适合黑河上游植被生长。该成果为高海拔地区植被生长变化的驱动因素分析提供了参考和借鉴。