黄土高原沟道侵蚀分析研究.docxVIP

? ? 陕西省黄土高原沟道侵蚀分析研究 ? ? 左飞航，张少伟，权红花 （1. 自然资源部陕西基础地理信息中心，陕西 西安 710054） 黄土高原的土壤侵蚀现象是一个重大的环境问题[1]。陕西省黄土高原地区的地形异常破碎，植被覆盖度较差，土壤侵蚀极为严重[2]。陕西省北部是我国最典型的黄土高原区域，研究该区域的土壤侵蚀特征，对于整个黄土高原来说极具代表性。另一方面，总书记提出黄河流域生态保护和高质量发展的战略部署，使得该区域的研究更具时代意义。本文利用侵蚀沟道数据、地理国情监测数据[3]、数字高程模型（DEM）和遥感影像数据等，首先评价了研究区的土壤侵蚀状况；然后以侵蚀沟道为单元，对陕西省黄土高原土壤侵蚀强度和分布特征进行了分析研究。 1 研究区概况与数据来源 1.1 研究区概况 陕西省分为北部黄土高原区、中部关中平原区和南部秦巴山区，其中黄土高原区在陕西省境内共包括6个地市、46个县区，地貌类型为高原沟壑区和丘陵沟壑区，总面积约为10.6 万km2。陕西省黄土高原地区水土流失严重，属于重点防治区域。20年来我国大范围推广退耕还林还草工程，陕西省黄土高原的植被覆盖不断增加，已经极大改善了区域土壤侵蚀情况，水土流失也呈现出了新的特征，因此有必要利用新的数据和方法进行研究和总结。 1.2 数据来源 近年来，我国开展了水利普查[4]、地理国情监测[5]等自然资源调查工作，积累了丰富的自然资源数据。本文利用全国第一次水利普查数据获取陕西省黄土高原侵蚀沟道数据；利用基础性地理国情监测成果获取土地利用数据，提取耕地数据；利用10 m分辨率的精细化DEM计算坡度；利用中分辨率遥感影像Landsat8 OLI 影像估算植被覆盖度；进而研究侵蚀沟道的土壤侵蚀强度特征。 2 研究方法 2.1 侵蚀沟道数据获取方法 沟道数据的提取以DEM和高分辨率遥感影像为数据源，首先利用空间水文分析方法进行初步提取，再采用人工后期处理的方式，提取侵蚀沟道范围面和沟道线[6]。空间水文分析的具体步骤包括在DEM 数据的基础上进行洼地填充、水流方向提取、汇流累积量计算、参考阈值分区域选取、栅格河网生成、出水口数据生成、原始沟道面数据生成、栅格数据矢量化等[7]。侵蚀沟道范围面是指侵蚀沟的汇水区域。提取沟道线层时，按照由高到低的方向提取矢量线，并计算高差。 2.2 土壤侵蚀评价方法 本文采用SL 190-2007《土壤侵蚀分类分级标准》，主要因素包括坡度、土地利用类型和植被覆盖度。土壤侵蚀强度分为轻度、中度、强烈、极强烈、剧烈5个等级，具体如表1所示。 表1 土壤侵蚀分级指标 1）坡度计算方法。利用DEM 计算坡度，并按照表1划分坡度等级。 2）土地利用类型提取方法。地理国情监测中的地表覆盖包括种植土地（01）、林草覆盖（03）、房屋建筑区（05）、铁路与道路（06）、构筑物（07）、人工堆掘地（08）、荒漠与裸露地表（09）、水域（10）等类别，括号中表示类别代码。本文将种植土地（01）中的水田（0110）和旱地（0120）作为耕地数据，其他类作为非耕地数据。 3）植被覆盖度计算方法。利用Landsat8 OLI影像计算归一化植被指数（NDVI），首先选取2015—2018 年植物生长期（7 月）的影像，再利用近红外波段和红光波段计算NDVI 值。 采用二分法估算植被覆盖度，将植被指数换算为植被覆盖度，即 式中， NDVImin、 NDVImax分别为最小、最大NDVI 值。 在影像上选取20个无植被覆盖度的裸土区和沙漠区，再选取20 个高植被覆盖度的植被区作为样本区域，每个区域不少于4 个像元；分别计算无植被覆盖度和高植被覆盖度样本区域的NDVI 平均值；取无植被覆盖度和高植被覆盖度样本区域的NDVI 算术平均值作为NDVImin、NDVImax。 2.3 统计分析方法 1）频数分析。频数是指某一特征出现的次数。本文对侵蚀沟道的面积和沟道线长度进行分级，并统计了每个级别对应沟道的频数。 2）空间（纬度）分布。为了研究侵蚀特征随纬度的变化情况，本文将研究区域按纬度划分为若干纬度分带，统计了每个分带上的侵蚀沟道面积、沟道线长度、沟道纵比、侵蚀强度分级等指标。 3）相关性分析。本文中各项指标的统计值均随纬度的变化而呈现出规律性的变化，且多为线性变化，因此采用决定系数R2表示拟合的相关性。 2.4 沟道侵蚀强度分析方法 1）重度以上侵蚀比。将重度、强烈、极强烈侵蚀定义为重度以上侵蚀，将沟道内重度以上侵蚀面积与沟道面积的比例定义为重度以上侵蚀比。 式中，S1为重度以上侵蚀面积；S 为侵蚀沟道面积。 2）沟道侵蚀强度。以沟道为单元，以侵蚀强度等级的面积为权重，计算各等级的加权平均值，作为沟道侵蚀强度。 式中， Di为侵蚀强度对应的等级，无侵蚀（0）、轻度侵蚀（1）、