Monitoring of vegetation coverage change in hilly and gully regions on the Loess Plateau

安塞区(36°30′45″~37°19′03″N，108°05′44″~109°26′18″E)位于陕西省北部，延安市北部，地处西北内陆黄土高原腹地，属典型的黄土高原梁峁丘陵沟壑区。土地总面积2 950 km²(其中90%为坡地^[14])，地势西北高、东南低，海拔为997~1 731 m。地貌类型以黄土梁涧、黄土峁状丘陵、黄土梁峁丘陵和沟谷阶地为主。气候属中温带大陆性半干旱季风气候，年均气温8.8 ℃，年均降水量505.3 mm，降水多集中在7−9月，年均蒸发量1 000.0 mm，远大于降水量，干旱频发。土壤以黄绵土、绵沙土等黄土性土和灰褐土为主，缺乏有机质和黏粒物质，胶结能力差，抗蚀能力低。因此，研究区水土流失严重，水土流失面积比例约97%，属于强度水土流失区。处于森林草原植被带向暖温带落叶阔叶林带过渡地带，植被及森林资源分布南北差异明显。天然次生林集中分布在南川四乡，人工林由小叶杨Populus simonii、山杏Armeniaca sibirica、刺槐Armeniaca sibirica及白榆Ulmus pumila等树种构成^[15]，多分布在北部。

以Landsat卫星影像为数据源，由中国科学院地理空间数据云平台下载(http://www.gscloud.cn/)2000年的Landsat TM影像和2017年的Landsat OLI影像，遥感影像空间分辨率为30 m，轨道号分别为127/34和127/35。2期影像成像时间均为树木枝叶繁茂的6月，质量较好，没有山区积雪、云和雾等影响，而且成像时间相近，避免不同植被生长状况造成的误差，有利于开展植被覆盖度的遥感分析比较。为了开展植被覆盖度地带性特征分析，本研究还采用了地面分辨率为30 m的ASTER GDEM V1数字高程产品。

在ENVI 5.1软件支持下对影像进行预处理。对影像进行辐射定标和大气校正，消除因传感器、大气条件、太阳位置和角度等引起的测量值与真实光谱反射率之间的差异。研究区沟壑纵横，复杂地形下遥感影像受传感器方位与太阳高度、方位等影响，阴坡接收到的照度较弱，亮度值较低，阳坡接收到的照度较强，亮度值较高。地形校正是复杂地形区遥感影像预处理的重要步骤。本研究应用了SCS+C地形校正^[16]，消除地形起伏引起的影像辐射亮度值的变化，使影像更好地反映地物的光谱特征。最后采用相同投影的矢量边界数据裁剪获得研究区的遥感影像。

像元二分模型是植被覆盖度估算中常见的模型之一，其本质是简单的线性像元分解模型^[17]。一般情况下，像元包含的遥感信息与大气、土壤和植被类型等因素相关。像元二分模型设置了S_veg和S_soil 2个参数，S_veg反映了像元内的植被信息，包括植被结构和种类等，S_soil反映了包括土壤湿度、亮度、颜色和类型等与土壤相关的像元信息，从而削弱了大气、植被和土壤等因素对应影像信息的影像，保留了植被覆盖度信息。像元二分模型简单而且结果可靠，采用该模型估算植被覆盖度。归一化植被指数(I_NDV)与植被覆盖度之间存在较好的线性关系，能反映植被生长状态及植被覆盖等信息，构建的植被覆盖度计算公式^[18]为：

式(1)中：F_C为植被覆盖度(%)；I_NDV为归一化植被指数；I_NDVmin为无植被覆盖像元的I_NDV值；I_NDVmax为完全被植被所覆盖的像元的I_NDV值。根据I_NDV统计直方图，以1%确定阈值。

计算出来的植被覆盖度影像，将大于1的赋值为1，小于0的赋值为0，然后对植被覆盖度进行分级。根据已有研究^[19]，将植被覆盖度分为5级：低植被覆盖度(F_C＜10%)、较低植被覆盖度(10%≤F_C＜30%)、中等植被覆盖度(30%≤F_C＜50%)、较高植被覆盖度(50%≤F_C＜70%)和高植被覆盖度(F_C≥70%)。

对2000和2017年植被覆盖度采用差值法获取植被覆盖度变化值(D)，设置阈值进行标准化分级，将全区分为7个植被覆盖度变化等级：D≤−30%为极显著减少，−30%＜D≤−10%为显著减少，−10%＜D≤−5%为较显著减少，−5%＜D≤5%为无显著变化，5%＜D≤10%为较显著增加，10%＜D≤30%为显著增加，D＞30%为极显著增加。

地形是植被分布的重要影响因素。利用ArcGIS软件从DEM数据中提取海拔、坡度和坡向数据，结合研究区实际情况进行分级，与植被覆盖度叠加，分析植被覆盖度在不同地形环境中的分布和变化差异。安塞区海拔为997~1 731 m，将海拔分为6级。结合地形对植物长势的影响，借鉴水土流失调查中采用8°作为缓坡和斜坡界线的方法^[20-21]，将坡度分为6级；依据地球公转以及太阳入射角，北半球西南坡最干燥，东北坡最潮湿，按湿润程度由高到低依次为东北向、北向、西北向、东向、西向、东南向、南向、西南向，考虑坡向对植被长势的影响，将其分为5级^[22-23]。地形要素分类情况见表1。

2000和2017年安塞区平均植被覆盖度分别为24.98%和53.34%，呈明显上升趋势(表2)。利用2018年8月研究区内40个样点测量植被覆盖度数据进行精度验证，相关系数为0.729 0，说明估算值与实际值具有一定的相关性。2000年安塞区植被覆盖度整体低，其中较低植被覆盖度占57.57%，低植被覆盖度占18.18%，中等植被覆盖度占16.34%。较高植被覆盖度和高植被覆盖度均很少，面积分别占4.28%和3.63%。由2000年安塞区植被覆盖度分级图(图1)可见：安塞区植被及森林资源分布南北差异显著，南部存在高植被覆盖度的天然次生林，北部植被覆盖度较低。

Figure 1.  Vegetation coverage classification in 2000 and 2017

2017年安塞区植被覆盖度获得较大改善，植被覆盖度总体上较2000年大幅度提升。低和较低植被覆盖度面积减少，其中低植被覆盖度面积仅占2.95%。中等、较高和高植被覆盖度面积分别占27.65%、25.41%和27.47%。由2017年植被覆盖度分级图(图1)可见：安塞区植被覆盖度呈现由南向北递减趋势，北部主要是低和较低植被覆盖度区域，南部是高植被覆盖度的天然次生林地带，河流沿岸和城镇周边为较高植被覆盖度区。

2000−2017年安塞区植被覆盖状况明显好转，区域内植被覆盖度呈现不同程度的增加(图2)。植被覆盖度显著和极显著增加面积占77.38%，其中极显著增加比例较大，占44.70%。植被覆盖度极显著增加区域分布在安塞区的河流沿岸，植被覆盖度显著增加区域主要集中在安塞区的西南部和北部。安塞区有177.83 km²植被覆盖度较显著减少，植被覆盖度无显著变化区域集中在南部的天然次生林地带，这些区域植被覆盖状况良好，改善空间不大。

Figure 2.  Vegetation coverage change classification map in 2000 and 2017

人类活动是植被覆盖变化的影响因素。1999年以来黄土高原大规模实施的退耕还林(草)、“三北”防护林、天然林保护等一系列重点林业工程是植被覆盖度增加的重要因素之一。根据《中国林业统计年鉴》获得安塞区2001−2017年的造林数据(表3)。17 a间，安塞区在退耕还林、天然林保护等林业生态工程的带动下，完成人工造林1 148.95 km²，飞播造林74.18 km²，封山育林99.53 km²，修复退化林13.34 km²。安塞区始终将人工造林、飞播造林和封山育林三者结合，增加森林植被覆盖度。2001−2003年造林面积最多，但2004−2006年，由于干旱缺水，大量沙棘Hippophae rhamnoides人工林死亡导致植被覆盖度明显下降^[24]。封山育林按照自然规律，利用植被自我恢复功能培育森林资源，对迅速恢复和扩大植被意义重大。此外，自2016年起安塞区采用更新改造技术开展退化林修复，以提升森林质量，实现森林可持续经营。

根据不同海拔植被覆盖度的统计(表4)，2000年安塞区植被覆盖度随海拔的增加大致呈降低趋势，2017年植被覆盖度随海拔升高先增加后减少，在海拔1 100~1 200 m内植被覆盖度最高，较高和高植被覆盖主要集中在海拔1 100~1 300 m区域。

比较各海拔等级的植被覆盖度，2000年研究区整体植被覆盖度低，各海拔等级内，较低植被覆盖度所占面积比例最大，为49.84%~60.32%；高植被覆盖面积比例最少，为1.88%~5.44%。2017年海拔＜1 300 m区域，面积比例最大的是高植被覆盖度；海拔＞1 300 m区域，中等植被覆盖度面积比例最大；海拔＞1 500 m时，30.10%为较低植被覆盖度。2000和2017年，各海拔等级间的植被覆盖度面积随海拔增加呈先增加后减少的趋势，拐点海拔区间为1 200~1 400 m(表4)。

根据不同坡度植被覆盖度的统计(表5)，2000年安塞区植被覆盖度随坡度增加逐渐增加，增幅很小；2017年植被覆盖度随坡度增加先增加后减少，坡度为25°~35°时，植被覆盖度最大。

比较各坡度等级的植被覆盖度，2000年较低植被覆盖度面积比例最大，为52.90%~63.67%；较高植被覆盖度面积较少，高植被覆盖度面积仅占1.97%~3.78%。2017年，坡度＜15°和＞35°区域中，中等植被覆盖度面积比例最大，在27%左右；坡度15°~25°、25°~35°，高植被覆盖度面积比例最大，分别为29.04%和30.21%。各植被覆盖度面积随海拔增加呈先增加后减少的趋势，15°~25°坡度的面积最大，但2017年低植被覆盖度在8°~15°坡度时，面积最大(表5)。

根据不同坡向植被覆盖度的统计(表6)，安塞区植被覆盖度由高到低依次为阴坡、半阴坡、半阳坡、平地和阳坡。分析各坡向等级内的植被覆盖度，2000年各坡向内，较低植被覆盖度面积比例为50%以上，但2017年，平地、阴坡和半阴坡中以高植被覆盖度为主要类型，面积比例分别为28.84%、34.8%和30.78%；半阳坡、阳坡以中等植被覆盖度为主要类型，面积比例分别为32.63%和29.54%。说明阴坡、半阴坡植被覆盖改善优于阴坡、半阴坡。2000年，阴坡植被长势好，面积大，但较低植被覆盖度在阳坡面积为477.43 km²，大于阴坡的390.71 km²。2017年，低和高植被覆盖度在阴坡面积最大，较低、中等和高植被覆盖度在阳坡面积最大。

地形控制了太阳辐射和降水的空间再分配过程，一定程度上制约了植被的生长和空间分布，尤其是在地形复杂的山区。黄土高原地域广阔、地形复杂，由于地理本底(如水热条件)和社会经济因素(如土地利用类型)的影响，不同地形类型区域的植被格局不同。

海拔通过影响温度、湿度等资源的再分配，以及人类活动频率影响植被覆盖率。安塞区海拔＜1 100 m的区域，地势平坦、水热条件好，人工造林以及自然植被恢复的成活率高，所以植被覆盖度很快得到提高。随海拔升高，气温降低，海拔1 400 m以上主要为丘陵沟壑分布区，沟壑纵横，植被稀疏，侵蚀剧烈，土壤保水能力差，有机质少，给人工造林和植被自然恢复增加了难度，植被覆盖度增幅受限。

坡度影响地表的物质流动和能量循环，一般坡度越大的地区，积温越少，投影面积相同条件下降水量越小，土层持水性能也越差。研究表明：15°~35°坡度内，土壤有机质含量最高，蓄水条件好，是植被覆盖的稳定区域；而坡度＜5°时，人类活动在平坦土地上较为频繁，对植被覆盖有一定影响^[25]。安塞区坚持实施坡度＞25°以上的坡耕地先退耕还林，所以15°~35°坡度范围内，植被恢复好，以高植被覆盖度为主；坡度＜15°区域人类活动频繁，耕地面积大，使得植被覆盖度较低，以中等植被覆盖度为主；坡度＞35°区域，植被生长受水热条件限制，所以中等植被覆盖度面积比例最大。

坡向影响坡面接受的太阳辐射以及地表与风向的夹角，从而影响光、热、水、土等因子的分配^[26]。安塞区处于干旱半干旱区域，阴坡日照少、温度低，蒸发量较少，土壤保持相对较高的湿度，所以各坡向植被覆盖度由高到低依次为阴坡、半阴坡、半阳坡和阳坡。2017年平地、阴坡和半阴坡以高植被覆盖度为主，半阳坡、阳坡的中等植被覆盖度面积比例最大，符合此规律。近年来，安塞区重视重点区域的荒坡治理，采取承包造林的方式，阳坡栽植侧柏Platycladus orientalis，阴坡栽植油松Pinus tabuliformis的思路，选用营养钵苗木造林，局部区域治理成效显著。安塞区内阴坡、半阴坡和半阳坡面积相近，但2017年较低、中等和高植被覆盖度在阳坡面积最大，说明近年来，人工造林使得阳坡的植被在较低、中等和高植被覆盖度上的改善程度大于其他坡向^[27]。

安塞区是黄土高原生态环境治理的重点区域，植被覆盖度从2000年的24.98%增长到2017年的53.34%，河流沿岸植被显著增加，整体植被呈改善趋势。地形是影响植被覆盖度时空动态变化的重要自然因子。安塞区内，坡度25°以上区域治理重点是保持水土，应根据山地肥力引种沙棘、苜蓿Medicago、紫穗槐Amorpha fruticosa等根系发达的植物。海拔1 300~1 400 m、坡度15°~25°的区域占整个研究区面积比例的35%以上，这些区域应采取不同林业工程的实施模式，适地适树、因地制宜，才能加快实现黄土高原大绿化进程。