收稿日期:2004-07-12;修订日期:2005-01-18基金项目:北京市自然科学基金重点项目“官厅库区土壤水分承载量及植被恢复机理研究”(6011003)。作者简介:谭炳香(1966-),女,副研究员,研究方向为遥感在林业中的应用,现主要从事遥感和GIS技术的森林资源监测与评估,森林参数遥感估测等方面的研究。

基于遥感数据的流域土壤侵蚀强度快速估测方法谭炳香1,李增元1,王彦辉2,于澎涛2,柳立兵3(1.中国林业科学研究院资源信息研究所,北京　100091;2.中国林业科学研究院森林生态环境与保护研究所,北京　100091;3.北京市延庆县水利局水土保持工作站,北京　102100)

摘要:以北京延庆县境内的妫水河流域为例,提出了一种基于遥感数据的土壤侵蚀强度快速估测方法。首先,利用遥感数据和植被指数模型提取流域内土地利用类型信息和植被覆盖度信息;其次,利用数字高程模型数据生成坡度图;然后,结合土壤侵蚀强度分级指标,将坡度图与土地类型图、植被覆盖度图空间叠加,判断和计算侵蚀强度等级,结果获得了流域土壤侵蚀强度等级图;最后,计算了流域的年平均侵蚀模数。结果表明,妫水河流域的土壤侵蚀以微度和轻度为主,所占面积比例为74.88%,极度和剧烈侵蚀很少,不到总面积的2%。整个流域的年侵蚀模数估计为1746.1t/km2·a。

关　键　词:遥感;植被覆盖;土壤侵蚀;流域中图分类号:TP79 文献标识码:A 文章编号:1004-0323(2005)02-0215-06

1　引　言植被具有截持降雨、减缓径流、保土固土等生态功能,尤其是对控制土壤流失起着决定性的作用,植被盖度(特别是地面覆盖度)的大小直接影响着土壤流失的强弱。已有的大量对比观测和试验发现,植被防治土壤侵蚀的能力随植被覆盖度的增加而增强。在其它条件一定时,侵蚀量与植被覆盖度成反比关系〔1〕。因此,进行植被覆盖调查和监测,是预测水土流失和制定水土保持规划的最基础的工作。但是,常规的植被覆盖地面调查费时、费力、成本高,而且不利于大范围的植被监测调查。近些年来,遥感技术的快速发展和广泛应用为植被覆盖信息的提取提供了可靠的技术途径。相对于传统的地面调查,遥感方法具有耗资少、周期短、覆盖广等特点,遥感信息已成为无实测资料的广大地区植被覆盖估算的一种极其重要的信息源〔2～5〕。目前,有关植被覆盖的遥感研究主要集中在植被类型分类制图方面,而对植被覆盖度的应用研究相对不足。应用遥感技术能够及时、快速地获取流域内的大范围区域植被覆盖和土壤侵蚀强度的空间分布状况,为流域水土保持规划提供科学依据,加强水源涵养型植被建设,改善水源区生态环境,增强植被的水源涵养与土壤保护功能。因此,本文以北京市延庆县的妫水河流域为例,试图提出一种基于遥感数据的土壤侵蚀强度快速估测方法,包括流域植被覆盖遥感监测、流域土壤侵蚀强度图绘制、流域年平均侵蚀模数计算等环节。

2　试验区及其数据获取2.1　试验区概况官厅水库是北京市仅次于密云水库的第二大水源,是首都水资源保护和水环境治理的重点地区。妫水河是官厅水库的一级流域和主要水源之一,主要分布在北京市延庆县境内,代表着华北半干旱区的水源地,具有相对丰富的前期调查研究基础,开展华北半干旱区水源涵养型植被研究较为理想。妫水河流域是官厅水库主要水源之一,位于北京延庆县境内。延庆县属延怀盆地的一部分,北、东、南三面环山,一般海拔600～800m,西面是官厅水库。平原面积5.4万hm2,占总面积的27.2%,山区面积14.5万hm2,占总面积的72.8%,是华北平原到坝

第20卷　第2期2005年4月遥　感　技　术　与　应　用REMOTESENSINGTECHNOLOGYANDAPPLICATION Vol.20　No.2Apr.2005上高原的过渡地带。县内山脉大致走向北东与东西向,由中部北起佛爷顶,经九里梁至凤坨梁,形成一自然分水岭,分水岭以西为山前平原区,以东为山后区。平原区的高程均在510m以下,由东北至西南倾斜,妫水河从东向西流入官厅水库,总流域面积1073.64km2,地形上构成一个典型的开启式盆地。除平原谷地耕作区外,延庆县大部分地区植被属落叶阔叶林向针阔混交林过渡的植被带,除落叶阔叶林外,还有暖性针叶林,次生的落叶灌丛和草本群落,其中以次生的落叶灌丛占最大面积。其分布特点为:盆地边缘有喜温耐旱的灌草植物,南部有少量喜温的果、柿树种。浅山丘陵大部分阳坡植被稀疏,覆盖度小,35°以上陡坡岩石裸露,少灌无草,水土流失严重,沟壑纵横,人为破坏严重,植被残破,为稀疏植被。2.2　遥感数据获取获取遥感图像时,时相的选择至关重要,因为最佳时相的选择是提高侵蚀调查与监测质量的基础和保障〔5〕。因此,我们选用了2001年5月19日的Landsat5TM多波段遥感图像(空间分辨率30m)和全色波段图像(15m)。5月中下旬的北京延庆地区乔灌木树种已展叶并开始旺盛生长,而农作物尚未覆盖地表,坡耕地在此季节较易从图像上判读出来,故土地类型的分类效果会最佳。以最新版的1∶5万的地形图为基准,采用人机交互方式选择地面控制点,利用二次多项式和最近邻内插法对图像进行几何精度纠正,经重新选点检验,将误差控制在一个像元内。利用妫水河流域的边界线对遥感影像进行掩膜,得妫水河流域的TM图像。将TM多波段图像与全色波段图像融合,获得15m空间分辨率的多波段图像。利用计算机自动分类和目视判读相结合的方法从遥感图像中提取土地利用类型信息,利用植被指数模型求算植被覆盖度。2.3　流域数字高程模型从延庆县水利局获得了数字化的妫水河流域1∶1万等高线数据,利用遥感图像处理软件中的地形分析工具生成了妫水河流域的数字高程模型(DEM),空间分辨率采样为与TM图像一样大小的15m;再由DEM数据生成坡度图。3　技术路线与方法3.1　技术路线以辐射纠正和几何纠正的妫水河流域遥感图像为信息源,提取土地利用类型信息和植被覆盖度信息;利用遥感影像处理软件由1∶1万数字高程模型(DEM)数据生成坡度图。最后,结合土壤侵蚀强度分级指标,将坡度等级图与土地类型图和植被覆盖度等级图分别进行空间叠加,最后生成土壤侵蚀强度等级图,统计各侵蚀强度等级的面积,并计算年平均侵蚀模数(如图1所示)。

3.2　植被覆盖度的估算植被覆盖度是指植被冠层的垂直投影面积与土地面积之比,它是描述森林植被覆盖特征的一个重要基础数据,广泛应用于森林水文、森林生态等方面的研究,也是许多生态模型中的基本参数。准确、快速地获取地表植被覆盖度及其变化信息,对揭示地表特征时空变化规律、探讨变化的驱动因子、分析评价区域生态环境等都具有重要的现实意义。植被覆盖度的监测方法分为地面实测和遥感监测两类。由于植被覆盖度具有显著的时空分异特性,大范围内森林植被覆盖的高密度地面动态实测的方法是不现实的,因而遥感监测已成为区域植被覆盖动态监测的主要手段〔6,7〕。遥感估算植被覆盖度的方

法大致可归纳为两种:经验模型法和植被指数转换法。其中,经验模型法主要是通过建立实测植被覆盖度数据与植被指数之间的经验模型来求取大面积植被覆盖。由于经验模型法基于特定区域的实测数据,在局地使用效果可能较好,但在大面积推广应用时则受到限制。植被指数转换法则是在对光谱信号进行分析的基础上,通过建立植被指数和植被覆盖度的转换关系来直接估算植被覆盖度。由于植被指数转换法不依赖于植被覆盖度的实测数据,因而实际工作中使用较多。目前比较广泛使用的是Gutman

216 遥　感　技　术　与　应　用 第20卷提出的一种根据归一化差值植被指数(NDVI)来计算植被覆盖度的方法。其基本原理是根据健康植被在红光波段有一个吸收峰,在近红外波段有一个高反射峰,而土壤在可见光到近红外波段呈近似线性变化的差异,首先计算出植被指数NDVI,然后再假设遥感器所测得的反射辐射R完全由植被总反射辐射和土壤总反射辐射构成〔7〕,即:R=fRv+(1-f)Rs(1)其中:Rv、Rs分别为植被、土壤的总反射辐射;f为植被覆盖度。则:f1=(R-Rs)/(Rv-Rs)或f2=(d-ds)/(dv-ds)(2)其中:d为植被与土壤混合光谱反射率;dv、ds分别为纯植被和纯土壤宽波段发射率。基于上述原理,利用植被指数可以计算植被覆盖度,即:　f3=(NDVI-NDVIs)/(NDVIv-NDVIs)(3)其中:NDVI为所求地块或像元的归一化植被指数;NDVIv、NDVIs分别为纯植被和纯土壤的植被指数。该式适合用于植被覆盖处的叶面积指数较均一的场合。也可得:f4=(NDVI-NDVImin)/(NDVImax-NDVImin)(4)其中:NDVI为所求像元的归一化植被指数;NDVImin、NDVImax分别为研究区内NDVI的最小、最大值。一般说来,归一化植被指数NDVI比较真实地表现了影像数据上植被的分布。但由于NDVI本身无法完全分离植被和土壤的影响,不管是裸土还是部分植被覆盖的像元,其NDVI值都可能是一个较小的整数。所以在实际应用中,往往根据经验确定NDVImin、NDVImax,或从计算得到的影像像元的NDVI中提出NDVImin、NDVImax

,从而计算植被覆

盖度。本文采取后一种方法获得NDVImin、NDVImax,然后求算植被覆盖度。3.3　土壤侵蚀强度的估测北京市水利局根据延庆县的地形、地貌、土地利用和植被状况制定了土壤侵蚀强度分级指标(见表1),将土壤侵蚀程度划分为6级:(1)微度侵蚀:无侵蚀或者土壤流失不明显。植被覆盖度>75%的成片林、灌、草地和坡度<5°的平地,包括山间、山前平地、河流阶地以及水体、裸地等,但不包括坡度3°～5°的耕地,只包括<3°的耕地。(2)轻度侵蚀:水土流失比较明显。坡度在3°～8°的坡耕地;或植被覆盖度为60%～75%、坡度为5°～25°的坡地;或覆盖度在45%～60%、坡度5°～15°的坡地;或覆盖度为30%～45%、坡度在5°～8°的坡地。(3)中度侵蚀:土壤侵蚀十分明显。坡度8°～15°的坡耕地;或者坡度5°～15°且植被覆盖度<30%的坡地;或者坡度8°～25°且植被覆盖度为30%～45%的坡地;或者坡度15°～35°且植被覆盖度为45%～60%的坡地;或坡度>25°且植被覆盖度为60%～75%的坡地。(4)强度侵蚀:土壤侵蚀强烈,主要分布于植被覆盖度比较少的山地。植被覆盖度为45%～60%且坡度>35°的坡地;或者植被覆盖度为30%～45%且坡度为25°～35°的坡地;或者植被覆盖度<30%,但坡度在15°～25°之间的坡地;坡度为15°～25°之间的坡耕地,皆为强度侵蚀。(5)极强度侵蚀:土壤侵蚀十分强烈,主要分布于山地高坡。植被覆盖度30%～45%且坡度>35°的坡地;或者坡度25°～35°但植被覆盖度<30%的表1　延庆县土壤侵蚀强度分级指标植被覆盖度坡度(°)0～33～55～88～1515～2525～35>35