王红岩, 李强子, 吴利桥, 丁雷龙,3, 张 磊, 杜 鑫

(1.中国科学院 遥感与数字地球研究所, 北京 100101;2.珠江水资源保护科学研究所, 广州 510611; 3.中国地质大学, 北京 100083)



基于RS和GIS的水土流失敏感性评价及动态监测

王红岩1, 李强子1, 吴利桥2, 丁雷龙1,3, 张 磊1, 杜 鑫1

(1.中国科学院 遥感与数字地球研究所, 北京 100101;2.珠江水资源保护科学研究所, 广州 510611; 3.中国地质大学, 北京 100083)

以Landsat-5 TM影像为主要数据源，结合高精度的土地利用数据和DEM高程，基于遥感和GIS技术，通过选取关键影响因子，对红水河干流区1990年，2000年和2010年的水土流失敏感性进行了评价及其时空动态监测分析。结果表明：红水河干流区水土流失敏感性主要以轻度敏感和中度敏感为主，1990—2010年近30 a来水土流失敏感性总体呈现好转趋势；在水土流失不同敏感性等级变化中，以水土流失中度敏感转为轻度敏感面积最大，1990—2000年和2000—2010年转化面积分别为2 736.24 km2和2 462.09 km2，1990—2000年水土流失敏感性向轻度敏感、不敏感转化的地区主要集中在红水河干流区中段区域，但2000—2010年水土流失敏感性等级的变化区域与之明显不同，主要发生在红水河干流区的东段。总体看来，在响应国家水土流失治理措施下，红水河干流区水土流失防治效果显著，生态环境逐渐转好。

水土流失; 敏感性评价; 动态监测; 遥感; GIS

水土流失是当今全球面临的一个严峻的环境与灾害问题，严重的水土流失已造成了生态环境的急剧恶化，受到了越来越多的关注。水土流失导致土壤表层营养成分流失、生产力下降、土地退化、陆地生态环境恶化等问题的发生，另外农业用地的化肥农药残留随着水土流失进入河道水库，是水环境恶化的主要污染源之一[1-3]。水土流失敏感性是指自然状态下发生土壤流失可能性的大小[4]，开展水土流失敏感性评价，有助于识别土壤侵蚀易形成的区域，评价土壤侵蚀对人类活动的敏感程度[5]，同时也为区域生态环境保护、生态安全维护、生态环境规划和管理提供科学依据[6]。近年来随着3S技术的发展，遥感技术在土地植被覆盖信息获取方面具有大面积、多时相、多波段、低成本的特色和信息丰富、动态性强的优势，GIS技术具有强大的多要素相关空间分析的强大特点。鉴于此，国内外学者基于3S技术，在水土流失敏感性评价及土壤侵蚀潜在危险度评价等方面开展了大量研究工作[7-12]。但是在我国红水河流域的研究相对较少，同时对于区域水土流失敏感性动态变化研究较少。

红水河流域资源富庶，自然条件复杂多样，但生态环境很脆弱。红水河流域是西南地区水土流失的重灾区，生产条件恶劣，生存环境恶化，改善生态环境已经刻不容缓。在此背景下，本文基于遥感和GIS技术，选取关键影响因子，开展红水河流域水土流失敏感性评价，探讨了水土流失敏感区的空间分布格局，同时分析了1990—2010年水土流失敏感性动态变化，以期为红水河流域水土流失治理、有效控制该区域土壤侵蚀提出参考建议，为评价区域的水土保持规划制定和生态环境建设提供重要的科学依据和技术支撑。

1 研究区概况

红水河是珠江流域西江水系干流，上游系发源于云南沾益县的南盘江，在贵州省望谟县庶香村与北盘江汇合后始称红水河，自西向东横穿广西中部，至象州县石龙镇三江口与柳江汇合后改称黔江。红水河干流全长659 km，流域面积3.3万km2。河流流经高原、低山和丘陵，其上游为三叠纪的砂页岩，中下游为石炭纪、二叠纪的石灰岩。红水河流域是典型的岩溶地貌区，岩溶地貌面积占山地总面积的65.1%。属亚热带季风气候区，年均气温为16.9～21.0℃，年降水量达1 244.0～1 590.7 mm。土壤类型主要包括赤红壤、红壤、黄壤、石灰土等。红水河流域矿产资源丰富，河川径流量大，河流落差大，蕴藏着丰富的水能资源。本文中研究区定义为红水河干流区，其涵盖干流区两侧第一道山脊线以内的范围，从大藤峡沿干流上朔至天生桥一级库尾段，包括了红水河上游10个梯级电站流域范围，总面积为14 193 km2，跨越广西、贵州、云南三个省区。

2 数据及方法

2.1 遥感数据及其预处理

本文以Landsat-5 TM数据为基本遥感数据源，该数据的空间分辨率为30 m，收集整理了覆盖红水河干流区的1990年，2000年，2010年三个时期的影像。获得红水河干流区数字高程模型DEM数据，其分辨率为30 m。影像的预处理主要包括大气校正和几何校正两个方面。大气校正采用ENVI/IDL软件提供的FLAASH大气校正模块进行，其基本原理是改进的MORTRAN4+模型。几何校正参考1∶50 000比例尺地形图，采用二次多项式进行校正，RMS误差控制在一个像元之内。

2.2 水土流失敏感性评价指标

红水河干流区主要为水蚀分布区，水土流失主要受降雨侵蚀力、土壤可蚀性、土地利用类型、地形、植被覆盖度、水土保持措施等因子的影响[13-14]。由于降雨侵蚀力和土壤可蚀性变化不大，本文主要选取坡度、植被覆盖度和土地利用类型为水土流失敏感性评价及其动态变化分析的指标(表1)。

表1 水土流失敏感性评价因子及其分级

2.2.1 坡度因子提取 根据水土流失敏感性评价分级标准，将研究区坡度划分为<8°,8°～15°,15°～25°,25°～35°,>35°共5个级别。通过获取的红水河干流区DEM栅格数据，在ArcGIS中Spatial Analysis工具条下的Surface Analysis中的坡度工具(Slope)功能来生成以度为单位的坡度图，基于划分标准进行等级划分，获得红水河干流区坡度等级分布图。

2.2.2 植被覆盖度因子提取 归一化植被指数(NDVI)与植被覆盖度存在显著的正相关关系，本文基于红水河干流区TM遥感影像，计算获得NDVI分布图，利用Gutman模型(见公式1)提取植被覆盖度，基于水土流失划分标准进行了等级划分，获得研究区植被覆盖度等级分布图。

(1)

式中：其中FC——植被覆盖度；NDVIsoil——完全是裸土或无植被区域的NDVI值；NDVIveg——完全被植被覆盖的像元的NDVI值。

2.2.3 土地利用类型因子提取 基于获取的红水河干流区1990年，2000年，2010年三个时期的影像数据，主要采用面向对象的SVM自动识别算法进行研究区土地覆被信息提取[15-16]。通过空间数据层建立、分类划分体系、专题信息提取及专题信息后处理等四大流程获得高精度的红水河干流区土地利用类型分布图。根据红水河干流区土地利用类型特点，研究区土地利用类型划分为针叶林、阔叶林、灌木林、人工林、草地、水田、旱地、居住地、工矿地、河流、水库/坑塘和裸露地12类。

2.3 水土流失敏感性综合评价

根据红水河流域的实际情况及各个因子对水土流失影响的大小，参考前人研究结果，确定红水河干流区水土流失敏感性评价各因子的权重，其中，坡度因子权重0.3，植被覆盖度权重为0.4，土地利用类型权重为0.3[17]。利用ArcGIS空间分析模块功能，采用多因子加权求和模型(见公式2),将各评价因子与其相应的权重叠加分析和代数运算，再采用自然分界法(Natural break)将运算结果划分5级(表2)，得到红水河干流区1990年水土流失敏感性分级图，其中，自然分界法是利用统计学的Jenk最优化法得出的分界点，能够使各级的内部方差之和最小[18]。根据1990年研究区自然分界法的划分阈值，对2000年和2010年的各因子加权求和结果进行划分，得到研究区2000，2010年水土流失敏感性分级图。

(2)

式中：S——水土流失敏感性综合评价结果；Wi——第i个评价指标因子的权重；Ai——第i个指标因子的数量化值；n——评价指标因子的个数(本文n=3)。

表2 水土流失敏感性综合评价分级

3 结果与分析

3.1 红水河干流区水土流失敏感性空间分布

基于获取的红水河干流区的TM遥感数据及DEM数据，依据水土流失敏感性分级遥感评价方法，获得1990年，2000年，2010年红水河干流区水土流失敏感性等级分布图，再分别对研究区水土流失不同敏感性等级的面积进行统计分析，分析结果见表3。从表3中可以看出，红水河干流区水土流失敏感性等级主要以轻度和中度为主，分别占了总面积的92.66%，94.83%，93.61%。在1990年代，水土流失中度敏感区所占比重最大，但在2000和2010年代，水土流失轻度敏感区所占比重最大，表明1990年代水土流失发生可能性要比2000年代、2010年代大。

从表3中可以看出，红水河干流区1990年，2000年，和2010年三个年代水土流失敏感性空间分布存在较大变化。在红水河干流区西段区域主要是原始阔叶林、灌木林和高山草地，人为扰动相对较少，植被覆盖度较高，在1990年代，水土流失轻度敏感区分布最多，在2000—2010近二十年来，水土流失敏感性等级变化不大，主要以轻度和中度水土流失为主。红水河干流区中段区域2000年代和1990年代水土流失敏感性存在明显变化，主要由原来的中度敏感区向水土流失轻度敏感区的转化。红水河干流区东段区域2010年代和2000年代水土流失敏感性等级存在明显变化，主要由原来的水土流失中度敏感区向轻度敏感和不敏感的转化。这主要是因为国家对水土流失治理工作的重视，加强了水土流失的治理，以及红水河流域的开发利用与生态环境治理保护相结合的开发战略，在红水河中、下游大力发展针叶、针阔混交林，建立速生丰产用林，在下游平原发展混农林业等措施减小了水土流失发生的可能性。

从表3中可以看出，近三十年来红水河干流区水土流失敏感性有所好转，中度和高度敏感区面积减少，水土流失不敏感和轻度敏感区面积增加。其中1990—2000年代水土流失敏感性等级变化明显，中度敏感区面积减少了1 388.12 km2，在2000—2010年代中度敏感区面积也有所减少，减少575.57 km2。同时，高度敏感区面积也存在较大幅度的减少，说明红水河干流区生态环境逐渐好转，生态质量得到较明显的改善。

表3 1990，2000，2010年代红水河干流区水土流失敏感性类型面积统计

3.2 红水河干流区水土流失敏感性的动态监测

对水土流失敏感性分级图做了进一步统计分析，表4和表5分别为1990—2000年代和2000—2010年代不同水土流失敏感等级的转移矩阵。从表4中可以得出，1990—2000年水土流失敏感性向轻度敏感、中度敏感和不敏感转化明显。1990—2000年，有2 736.24 km2的水土流失中度敏感区转为轻度敏感，有276.74 km2的水土流失轻度敏感区转为不敏感，并且主要分布在红水河干流区的中段以及东部的大藤峡水电站区段，这些区域加强了水土流失治理工作，石山荒山绿化，加大了植树造林的面积，控制了水土流失，改善了生态环境。同时2000年代的水土流失高度敏感区的增加主要是中度敏感区的转移，增加的面积达247.40 km2，主要分布在红水河干流区西段区域，说明这些区域植被大面积遭破坏，植被覆盖率下降，水土流失发生可能性大大增加。

从表5中可以得出，2000—2010年水土流失敏感性等级变化主要还是向轻度、不敏感和中度的转化，与1990—2000年变化一致，但变化分布区域有所不同。2000—2010年水土流失敏感性等级变化中转为不敏感和轻度敏感区域主要分布在红水河干流区东段。其中：408.95 km2的水土流失轻度敏感区转化为不敏感，但同时有256.17 km2不敏感区转为水土流失轻度敏感区；有2 462.09 km2的中度敏感区转变为轻度敏感，同时水土流失高度敏感区中有60.42 km2转为轻度敏感，有262.22 km2转为水土流失中度敏感区。总体来说，2000—2010年红水河干流区水土流失得到控制，水土流失敏感性有所减轻，生态环境逐渐变好。

表4 红水河干流区1990-2000年水土流失敏感性等级转移矩阵 km2

表5 红水河干流区2000-2010年水土流失敏感性等级转移矩阵 km2

4 结论与讨论

本文基于遥感和GIS技术，分析了红水河干流区1990—2010年近30 a的水土流失敏感性分布及其时空动态变化，研究结果为掌握该区域水土流失发生可能性及治理措施的实施提供一定的参考价值。研究结果表明：

1) 通过对红水河干流区1990年，2000年，2010年水土流失不同敏感性等级面积统计分析发现，红水河干流区水土流失敏感性主要以轻度和中度为主，分别占了总面积的92.66%，94.83%，93.61%。

2) 通过水土流失不同敏感性等级面积转移矩阵分析得出，1990—2010年来，红水河干流区水土流失敏感性向轻度、中度和不敏感转化明显，水土流失持续好转。式中：水土流失中度敏感区转为轻度敏感面积最大，1990—2000年和2000—2010年转化面积分别为2 736.24 km2和2 462.09 km2。但1990—2000年和2000—2010年水土流失敏感性向轻度敏感、不敏感转化的地区分布不同，在1990—2000年，变化主要集中在红水河干流区中段区域，而2000—2010年水土流失的转化主要发生在红水河干流区的东段。

3) 虽然近30 a来红水河干流区水土流失敏感性有所改善，但在有些区域水土流失中度敏感转为高度敏感，对于这些区域需要加大水土流失治理力度，采取植树造林等措施进一步有效的遏制研究区的水土流失发生的可能性。

4) 水土流失的影响因子是多方面的，本研究从地形坡度、植被覆盖度及土地利用类型等3个因子考虑，参考前人研究结果定义了水土流失敏感性评价中三个因子的权重，研究结果有一定的片面性。在今后的研究中，将从自然演化和人为干预两个方面，选取较为全面的影响因子进行评价，探索红水河干流区水土流失敏感性影响因子的权重计算与评价方法。

[1] 徐涵秋.水土流失区生态变化的遥感评估[J].农业工程学报,2013,29(7):91-97.

[2] 葛德祥,庞治国,王义成,等.二滩库区蓄水前后水土流失变化遥感监测与分析[J].中国水利水电科学研究院学报,2010,8(1):44-51.

[3] 贺缠生,傅伯杰,陈利顶.非点源污染的管理及控制[J].环境科学,1998,19(5):87-91.

[4] 胡良军,李锐,杨勤科.基于GIS的区域水土流失评价研究[J].土壤学报,2001,38(2):167-175.

[5] 陈建军,张树文,李洪星,等.吉林省土壤侵蚀敏感性评价[J].水土保持通报,2005,25(3):49-53.

[6] 张玉娟,刘丹丹,王延亮.基于RS和GIS东北农业主产区水土流失敏感性评价：以黑龙江省宾县为例[J].测绘与空间地理信息,2013,36(8):37-39.

[7] 张朝琼,郜红娟,张帮云.基于GIS的仁怀市生态敏感性评价[J].水土保持研究,2013,20(4):179-182.

[8] 张伟,王家卓,任希岩,等.基于GIS的山地城市生态敏感性分析研究[J].水土保持研究,2013,20(3):44-47.

[9] 朱志玲,吴咏梅,张敏.基于GIS的宁夏生态环境敏感性综合评价[J].水土保持研究,2012,19(4):101-105.

[10] 魏兴萍.基于RS和GIS的重庆南川区水土流失变化研究[J].水土保持研究,2009,16(5):60-65.

[11] 李婷,于青秀,张世熔.基于RS和GIS的涪江流域上游地区土壤侵蚀定量估算[J].四川农业大学学报,2011,29(1):84-88.

[12] 周璟,张旭东,何丹,等.基于GIS与RUSLE的武陵山区小流域土壤侵蚀评价研究[J].长江流域资源与环境,2011,20(4):468-474.

[13] 崔保山,杨志峰.湿地生态系统健康评价指标体系Ⅱ.方法与案例[J].生态学报,2002,22(8):1231-1239.

[14] 吴建国,常学向.荒漠生态系统健康评价的探索[J].中国沙漠,2005,25(4):604-611.

[15] Zhang L, Huang X. Object-oriented subspace analysis for airborne hyperspectral remote sensing imagery [J]. Neurocomputing,2010,73(4):927-936.

[16] Wu B, Xiong Z, Chen Y, et al. Classification of quickbird image with maximal mutual information feature selection and support vector machine [J]. Procedia Earth and Planetary Science, 2009,1(1):1165-1172.

[17] 闫利会,周忠发,王媛媛.丹霞地区水土流失遥感评价及空间分异：以贵州省赤水市为例[J].湖北农业科学,2012,51(20):4491-4495.

[18] 王万忠,焦菊英.中国的土壤侵蚀因子定量评价研究[J].水土保持通报,1996,16(5):1-20.

Dynamic Monitoring and Sensitivity Evaluation of Soil Erosion Based on RS and GIS

WANG Hongyan1, LI Qiangzi1, WU Liqiao2, DING Leilong1,3, ZHANG Lei1, DU Xin1

(1.InstituteofRemoteSensingandDigitalEarth,ChineseAcademyofSciences,Beijing100101,China; 2.ScientificInstituteofPearlRiverWaterResourcesProtection,Guangzhou510611,China; 3.ChinaUniversityofGeosciences,Beijing100083,China)

With the support of remote sensing and GIS, the sensitivity of soil erosion in Red River district was studied based on the Landsat-5 TM image and combining with high-precision land use data and DEM. Based on these, the dynamic monitoring of soil erosion sensitivity in 1990—2010 was analyzed. The results showed that the soil erosion sensitivity in Red River district was mainly at mild and moderate degree, the sensitivity of soil erosion had become better from 1990 to 2010. The largest change area of the soil erosion sensitivity was from moderate level to mild degree, the transformed areas of 1990—2000 and 2000—2010 were 2 736.24 km2and 2 462.09 km2, respectively. The soil erosion sensitivity conversion which was the other erosion sensitivity converted to mild, and micro soil erosion sensitivity from 1990 to 2000 mainly concentrated in the middle of the Red River district, however, the change area from 2000 to 2010 was different from that from 1990—2000, which was mainly in the east of the Red River district. Overall, in response to the national soil erosion control measures, the performed water and soil conservation measures in Red River district were very effective and the eco-environment gradually improved.

soil erosion; sensitivity evaluation; dynamic monitoring; RS; GIS

2014-04-28

2014-05-29

红水河综合利用规划环境影响后评价研究(1261320610403)

王红岩(1986—),女,山东省聊城市人,博士,助理研究员,主要从事生态遥感、荒漠化遥感监测与评价研究。E-mail:wanghy@radi.ac.cn

吴利桥(1970—),男,湖北通城人，高级工程师，主要从事水资源保护、环境规划、环境影响评价工作。E-mail:wulq2004@126.com

S157

1005-3409(2015)02-0064-05