[免费论文]喀斯特山区土地利用变化与土壤侵蚀耦合关系探究__以毕节试验区为例
主,耕地面积不断减少,林地和建设用地面积不断增加;土壤侵蚀以轻度.中度侵蚀为主,整体呈现改善的趋势,尤其是退耕还林后,土壤侵蚀的面积与强度均不断降低;土地利用变化与土壤侵蚀具有密切的联系,耕地是土壤侵蚀发生的主要土地类型,耕地的减少与土壤侵蚀呈正相关关系,林地面积变化与土壤侵蚀具有较好的相关性,但在退耕还林的初期,土壤侵蚀的强度略有增加,草地增加与土壤侵蚀负相关,增加草地有利于土壤抗侵蚀能力的增强.探究喀斯特山区土地利用与土壤侵蚀的耦合关系,有利于有限的土地资源的合理配置和喀斯特地区的生态保护,为区域土地规划及水土保持等宏观决策部门提供理论支撑.关键字:土地利用变化;土壤侵蚀;叠置剖析;转移矩阵;毕节试验区西南喀斯特山区具有山高坡陡.雨热同期.大面积碳酸盐岩发育.成土过程缓慢而复杂的特点,加之人多地少,人地矛盾突出,不合理的土地利用结构和生产生活方式造成了生态环境和土地资源破坏,植被破坏和土壤保肥性能减弱,使西南喀斯特山区土壤侵蚀更加严重[1-2].土地利用与土壤侵蚀具有密切的相互关系,在土壤侵蚀监测方面,应用广泛的美国修改的通用土壤流失方程(RUSLE),将植被.管理举措.水土保持举措等参数作为土壤侵蚀测定的重要指标[3-4];国内的我国水土流失方程(CSLE),也将生物.工程.耕作等土地利用因子作为土壤侵蚀测定的重要指标[5].运用RS和GIS结合探究两者的关系主要有两种方式:一种是计算不同土地利用类型的面积与土壤侵蚀分布面积的相关性,探究土地利用结构与土壤侵蚀的关系,如喻权刚探究黄土丘陵区发现坡耕地占比越大,年土壤侵蚀量越大[6];另一种方式是土地利用现状图与土壤侵蚀图在GIS空间剖析模块下叠置而发现不同地类的土壤侵蚀程度,如邹亚荣等以广东惠东区为例得出土壤侵蚀分布与地类类型分布是相关联的[7].然而,以往的探究多以某一年份探究土地利用类型与土壤侵蚀之间的耦合关系,长时间序列的土地利用类型变化与土壤侵蚀变化之间相互关系的探究还较为薄弱.探究喀斯特山区土地利用类型变化与土壤侵蚀变化的耦合关系,丰富了喀斯特山区土地利用与土壤侵蚀关系探究的理论,有利于有限的土地资源的合理配置和喀斯特地区的生态保护,为区域土地规划及水土保持等宏观决策提供理论支撑.1探究区概况与探究方法1.1探究区概况及数据来源毕节试验区是国务院于1988年6月批准建立的全国唯一用"扶贫开发.生态建设.人口控制"为主题的试验区.位于东经103°36′-106°44′到北纬26°21′-27°47′之间,地处我国贵州省西北部即乌蒙山的中心位置,靠近川.滇.黔三省交汇处,管辖1区7县,总面积为26853km2,占贵州省总面积的15.25%[8].本探究所用的遥感影像来源于地理空间数据云,包含1990年8月和9月Landsat5TM影像5景及1989年9月的Landsat5TM影像1景,2000年5月和6月的Landsat7ETM+影像6景,2010年8月的Landsat5TM影像6景,成像时间尽量选择各年份相同时间,以便于植被覆盖度及土地利用类型解译.除此之外,还收集了毕节试验区DEM数据以便于坡度专题图制作;矢量格式的行政区划范围数据用来裁剪遥感影像图.1.2探究方法1.2.1土地利用信息提取.结合全国土地利用动态遥感监测规程土地利用分类,参考中科院资源环境数据库的划分系统,结合遥感影像纹理色彩及毕节试验区实际土地利用情况,将毕节试验区土地利用类型划分成为耕地.林地.草地.水域.城乡工矿居民.未利用地6个一级地类,水田.旱地.灌木林地等9个二级地类[9].根据遥感影像的纹理.色彩.形状,结合相关专业图件和地理图件及毕节试验区土地利用分类系统建立毕节试验区遥感影像解译标志后根据解译标志采用人机交互解译方法分别对3期的遥感影像进行解译获得三期相应的土地利用类型和结构分布图(图1)及面积(表1).在ENVI中参考地面资料建立感兴趣区对地类解译结果进行精度评价,得到总体分类精度(OverallAccuracy),制图精度(Producer′sAccuracy),用户精度(User′sAccuracy)Kappa系数(KappaCoeffi-cient),分类混淆矩阵(confusionmatrix),错分和漏分误差等,最终得到OverallAccuracy>80%,KappaCoefficient高于最低精度0.7即符合要求.1.2.2土壤侵蚀信息提取(1)植被覆盖度提取.选用李苗苗改进后的像元二分模型[10],结合归一化植被指数(NDVI),对毕节试验区植被覆盖度进行提取.首先,在ENVI软件中分别对毕节试验区1990年.2000年.2010年计算得到的影像像元的NDVI,然后提取NDVIsoil,NDVIveg来计算各年份的植被覆盖度.同时,参照土壤侵蚀分类分级标准(SL190-2007)[11],在ArcGIS中按≤30%,30%~45%,45%~60%,60%~75%,>75%的标准,将植被覆盖度分为低.中低.中.中高.高度植被覆盖5大类,对应属性编码为1,2,3,4,5,最后,得到1990-2010年3个时间点的毕节试验区植被覆盖图.考虑到植被覆盖度的现实意义,本探究以收集到的已进行数据验证的1∶10万土地利用数据为基础做精度验证.验证结果表明:1990年.2000年.2010年植被覆盖度空间分布数据的总体精度为74.7%,61.9%,65.5%,所得植被覆盖度与用于精度验证的植被覆盖数据间有一定的差距.究其原因,一方面对影像阴影处理方式不同对结果产生影响;另一方面,2000年获取ETM+数据与1990年和2010年的TM数据不同,由于云层原因,遥感影像季相上的选择有所差别,这在一定程度上影响到最终结果.不过由于本探究中所用的人机交互解译,只要总体精度大于60%即满足要求.(2)高程及坡度信息提取.DEM是数字高程模型(DigitalTerrainModel)的简称,按≤5°,5°~8°,8°~15°,15°~25°,25°~35°,>35°的标准,将坡度分为平缓坡.中等坡.斜坡.陡坡.急坡.急陡坡等6大类,对应属性编码为1,2,3,4,5,6,在ArcGIS空间剖析模块(SpatialAnalyst)的支持下,对各DEM进行坡度提取及重分类等操作,生成高程及坡度专题图.(3)土壤侵蚀专题图制作.参考时新玲等的水蚀强度分级指标[12]和水利行业土壤侵蚀分类分级标准(SL190-2007)中水力侵蚀强度判别模型而建立毕节试验区土壤侵蚀强度的分级模型,根据该指标采用地类.植被覆盖度.坡度3个因子综合判读法制定土壤侵蚀专题图[13-14],将毕节试验区的土壤侵蚀强度分为微度侵蚀.轻度侵蚀.中度侵蚀.强度侵蚀.极强度侵蚀.剧烈侵蚀等6个等级.分别将提取出各年份的土地利用图.植被覆盖度图.坡度图做叠置剖析,在其属性表中进行逻辑判断,确定各斑块的土壤侵蚀强度等级,之后将同一属性的斑块间进行合并而最终得到探究区土壤侵蚀强度等级图(图2),并对土壤侵蚀面积进行统计(表2)后结合土地利用统计出不同地类中土壤侵蚀强度指数(表3).2结果与剖析2.1土地利用类型变化状况剖析从土地利用类型分布来看,毕节试验区土地利用类型以林地.耕地和草地为主,总面积达到探究区的99%以上,耕地面积达到了33%以上.从土地利用类型变化过程来看,1990-2010年,毕节试验区土地利用变化较为强烈,其中耕地面积持续减少,由1990年的9445.48km2减少为2010年的9119.76km2,年均减少32.57km2,侧面反映了试验区产业结构调整的成效,种植业在产业结构中的比重正在不断减少;林地与城乡工矿居民用地面积持续增加,其中林地面积由1990年的11296.28km2增加为2010年的11521.22km2,年均增加22.49km2,说明生态环境水平不断提升,尤其是退耕还林工程实施后,林地面积增加显着;草地和水域面积先减少后增加,说明人们对生态环境的重视程度正在不断提升;城乡工矿居民用地面积持续增加,由1990年的36.20km2增加为2010年的93.56km2,年均增加22.49km2,20a间增加了近2倍,说明试验区城镇化进程不断加快.2.2土壤侵蚀变化状况剖析1990-2010年,毕节试验区的土壤侵蚀强度以轻度和中度侵蚀为主,发生土壤侵蚀的面积比例分别为70.22%,71.79%,64.49%,年平均侵蚀模数为221.51t/(km2·a),228.83t/(km2·a),190.75t/(km2·a),土壤侵蚀强度总体呈现先恶化后改善的趋势.不同强度等级的水土流失面积来看,轻度侵蚀变化呈现先减少后增加趋势,其他侵蚀等级则相反,均为先增加后减少趋势.对3期土壤侵蚀分布图进行不同强度土壤侵蚀面积转移矩阵换算发现:1990-2000年,土壤侵蚀变化过程较为平缓,土壤侵蚀状况恶化的总面积(即低等级土壤侵蚀强度转变为高等级土壤侵蚀强度,如轻度土壤侵蚀强度转化为中度土壤侵蚀强度)为2932.97km2,土壤侵蚀状况改善的总面积(即高等级土壤侵蚀强度转变为低等级土壤侵蚀强度,如中度土壤侵蚀强度转化为轻度土壤侵蚀强度)为2055.65km2,总体侵蚀强度状况略有恶化;2000-2010年,土壤侵蚀变化过程较为剧烈,土壤侵蚀状况恶化的总面积为6866.89km2,土壤侵蚀状况改善的总面积为10772.25km2,土壤侵蚀强度在此阶段出现减弱趋向,土壤侵蚀状况得到改善.土壤侵蚀强度在2000-2010年有所减轻,得益于2000年后毕节试验区进行退耕还林.封山育林等生态工程举措的开展.2.3毕节试验区土地利用变化与土壤侵蚀耦合关系2.3.1不同地类的土壤侵蚀强度综合指数土壤侵蚀率反映侵蚀程度但无法反映出侵蚀强烈度,为此引入土壤侵蚀强度综合指数(无量纲,为计算方便放大100倍),对不同地类不同侵蚀情况定量剖析,该指标大小反映出不同地类中土壤侵蚀的侵蚀强度[15],由于篇幅影响,计算过程略.由于水域.城乡工矿居民用地.未利用地在土壤侵蚀判读时定微度侵蚀即不发生侵蚀,故仅列出耕地.林地.草地土地动态区与土壤侵蚀强度变化关系(见表3).结果显示:一级地类中耕地侵蚀最显着,土壤侵蚀强度综合指数最大,其次是草地,再次是林地.二级地类土壤侵蚀强度综合指数由大到小是:旱地>中覆盖草地>灌木林>疏林>有林地>水田>高覆盖草地>低覆盖草地>其他林地.旱地不论是侵蚀程度或者是强度,其土壤侵蚀发生率都最大;中覆盖草地的侵蚀强度大于灌木林地,但是灌木林地的土壤侵蚀程度大于中覆盖草地,因此要加强水土流失防治,应该更注重坡改梯.减少坡地耕种而多增加林草地等生态举措的实行,以防林草地衰退,应提升植被覆盖度.2.3.2土地利用变化与土壤侵蚀变化关系剖析利用转移矩阵剖析的方法,能够较好的模拟土地利用类型变化过程中土壤侵蚀的响应特征,揭示土地利用变化及土壤侵蚀变化的关系.首先,将1990-2010年3期土地利用类型分布图进行叠置剖析,分别获得1990-2000年土地利用类型变化图和2000-2010年土地利用类型变化图,同理,获得1990-2000年土壤侵蚀强度变化图和2000-2010年土壤侵蚀强度变化图;然后,将同期的土地利用类型变化图与同期土壤侵蚀变化图进行叠置剖析后统计出具体转移数据(土地利用动态区)及土壤侵蚀转移数据,并根据此数据计算出各土地利用转移方向下土壤侵蚀强度综合指数,从而反映相应土壤侵蚀的增减情况;最后,分别获得1990-2000年的不同土地利用转移过程中土壤侵蚀增减情况(图3),2000-2010年的不同土地利用转移过程中土壤侵蚀增减情况(图4).耕地是土壤侵蚀发生的主要地类,坡地开垦和坡地耕作是引起土壤侵蚀的主要原因.在耕地动态区,耕地面积持续下降,过度开垦状况有所好转,土壤侵蚀强度综合指数显着变小,表明耕地减少与土壤侵蚀呈正相关关系.在其林地动态区,林地持续增加,土壤侵蚀强度综合指数先增后降,林地增加与土壤侵蚀强度负相关,但由于退耕还林的初期,地表扰动大,林木的覆盖度还很低,土壤侵蚀略有增强,但随着林木覆盖度的增加,土壤侵蚀强度不断降低,应加强坡改梯.植树造林.减少坡耕地增加林草地等生态工程举措实施.在草地动态区,草地主要向林地转移,在退耕还草工程下,部分草地来自于耕地,有部分林地退化为草地,面积先减后增,土壤侵蚀强度综合指数先变大再变小,表明草地增加与土壤侵蚀负相关.增加草地有利于土壤抗蚀力的增强,但存在林地向草地退化现象会加重侵蚀程度,应加强草地的保护.3结论(1)毕节试验区主要土地利用类型为林地.耕地.草地.1990-2010年,毕节试验区土地利用变化较为强烈,林地面积持续增加,生态环境水平不断提升,尤其是退耕还林工程实施后,林地面积增加显着;耕地面积持续减少,反映了试验区产业结构调整的成效,种植业在产业结构中的比重正在不断减少;草地和水域面积先减少后增加,呈现了前期破坏后期整治的过程;城乡工矿居民用地面积持续增加,20a间增加了近2倍,尤其是近10a城镇化进程不断加快.(2)土壤侵蚀以轻度侵蚀与中度侵蚀为主,总体呈现先恶化后改善的趋势,尤其是退耕还林工程实施后,土壤侵蚀的面积和强度均不断降低,应该继续实施退耕还林工程,实施天然林保护,强化封山育林,充分发挥喀斯特山区生态环境的自然修复能力.毕节试验区不同地类的侵蚀大小为:旱地>灌木林地>中覆盖度草地>疏林地>有林地>水田>高覆盖度草地>低覆盖度草地>其他林地.(3)耕地是土壤侵蚀发生的主要土地类型,耕地面积持续下降,过度开垦状况有所好转,土壤侵蚀强度指数显着变小,侵蚀程度减弱,表明耕地的减少与土壤侵蚀呈正相关关系;在林地动态区,林地持续增加,土壤侵蚀强度综合指数先增后降,主要原因是退耕还林的初期,地表扰动大,林木的覆盖度还很低,土壤侵蚀略有增强,但随着林木覆盖度的增加,土壤侵蚀强度不断降低;在草地动态区,草地面积先减后增,土壤侵蚀强度综合指数先变大再变小,表明草地增加与土壤侵蚀负相关,增加草地有利于土壤抗侵蚀能力的增强,但存在林地向草地退化现象会加重土壤侵蚀,应加强草地的保护.参考文献:[1]张鲁,周跃,张丽彤.国内外土地利用与土壤侵蚀关系的探究现状与展望[J].水土保持探究,2008,15(3):43-48.[2]凡非得,王克林,熊鹰,等.西南喀斯特区域水土流失敏感性评价及其空间分异特征[J].生态学报,2011,31(21):6353-6362.[3]崔晓临,雷刚,王涛,等.退耕还林还草工程实施对洛河流域土壤侵蚀的影响[J].水土保持探究,2016,23(5):69-73.[4]孙德亮,赵卫权,李威,等.GIS与RUSLE模型的喀斯特地区土壤侵蚀探究[J].水土保持通报,2016,36(3):271-276.[5]俱战省,文安邦,严冬春,等.137Cs,210Pb和CSLE的三峡库区小流域土壤侵蚀评估[J].水土保持学报,2015,29(3):75-80.[6]喻权刚.遥感信息探究黄土丘陵土地利用与水土流失[M].郑州:黄河水利出版社,1998.[7]邹亚荣,张增祥,周全斌,等.GIS的土壤侵蚀与土地利用关系剖析[J].水土保持探究,2002,9(1):67-69.[8]张凤太,苏维词,赵卫权,等.生态足迹模型的喀斯特高原山地生态系统健康评价探究[J].水土保持通报,2011,31(1):256-261.[9]许月卿,蔡运龙,彭建.土地利用变化的土壤侵蚀效应评价[M].北京:科学出版社,2008.[10]李苗苗.植被覆盖度的遥感估算方法探究[D].北京:我国科遥感应用探究所,2003.[11]SL190-2007.土壤侵蚀分类分级标准[S].北京:我国水利水电出版社,2008.[12]时新玲,李智广.矢量和栅格数据结构的土壤侵蚀强度判别方法探究[J].水土保持通报,2002,22(1):34-38.[13]王志刚,涂人猛,孙佳佳,等.地类.坡度与植被的小流域土壤侵蚀量估算方法:以毕节市毛家湾小流域为例[J].长江科院报,2013,30(5):22-26.[14]蒋忠诚,杨德生,曹建华,等.我国水土流失防治与生态安全.西南岩溶卷[M].北京:科学出版社,2010.[15]甘彩红,赵岩洁,李阳兵.三峡库区腹地流域类型与土壤侵蚀的关系[J].水土保持通报,2015,35(6):291-295.5第6期谢颖颖等:喀斯特山区土地利用变化与土壤侵蚀耦合关系探究
原文链接：http://www.jxszl.com/lwqt/mflw/21605.html
