基于GIS和RS的荒漠区土壤侵蚀研究

武文一，刘 瑛，梁继业

（1．中国电力工程顾问集团中南电力设计院，武汉430071；2．三峡大学 水利与环境学院，湖北 宜昌443002；3．塔里木大学，新疆 阿拉尔843300）

水土流失和土地荒漠化已经成为当今研究的热点问题［1－11］。我国是土壤侵蚀最为严重的国家之一，其范围遍及全国各地。土壤侵蚀作为一种最活跃、最敏感的生态致灾因子［12］，关系着山区、丘陵区及风沙区水土资源的开发、利用和保护，关系着江河、湖泊的利用和整治，涉及整个区域的生态环境、经济持续发展和社会稳定［13］。研究区地处荒漠和绿洲的生态过渡带、生态环境脆弱，区域的土壤侵蚀状况与区域的稳定、发展密切相关。对研究区土壤侵蚀的研究，将为区域的水土流失治理和可持续发展规划提供参考依据。

1 研究区概况

研究区属高寒内陆沙漠型气候，由于受地理位置和地形的制约，区内终年干旱、少雨、多风、寒冷。《青海省地面气象资料30 a整编（1961－1990）》资料显示：区域年平均气温1．4℃，年平均降水量为83．5 mm，年平均蒸发量为2 176．6 mm。年平均风速达2．1 m／s，多为西北风，年平均日照时数为3 244．8 h，无霜期为108 d；平均极端低温为－34．2℃，冻土深度可达1．72 m。研究区植被为温带荒漠区、温带半灌木、灌木荒漠地带性植被，植被覆盖率低，主要种类有红砂、柽柳、骆驼刺、白刺、芨芨草以及碱篷等耐旱、耐寒、耐风沙的种类。区域土壤基质主要为粗细沙、沙粒和角砾。区域土壤侵蚀类型以风力侵蚀为主。该区域内人烟稀少，地形以戈壁滩、沙丘、高山为主。区内矿产资源丰富，工业以冶金、煤炭、电力化工为主。本研究选取人类活动较为剧烈的煤矿区及其周边区域作为研究对象。

青海鱼卡矿区位于青海省海西蒙古族藏族自治洲鱼卡乡。鱼卡乡地域东至羊肠子沟与德令哈毗邻，西起南八仙与冷湖行委相连，南靠察尔汗盐湖与格尔木相接，北侧翻过祁连山与甘肃省相连，地处东经94°49′38″－95°04′37″，北纬38°00′29″－38°09′40″。

2 数据来源及处理

解译使用的信息源主要为美国陆地卫星Landsat 5遥感影像T M，空间分辨率30 m，数据的获取时间为2005年8月28日。选取此时间段的遥感数据，主要是考虑到该时间段具有植被发育好、地表信息丰富的特点，有利于对各生态环境因子的判读。研究区土壤侵蚀类型、强度及分布的判定主要是在野外实地调查、定点记录的基础上，采用GIS建模运算来获取。其中，风蚀的判定是在野外定点调查的基础上，参考中国植被图集（1∶1000 000）、中国1∶100万土壤图、地形图（1∶50 000）、风速、降水等气象资料，对遥感影像进行目视解译来获取。水蚀的判别主要是在风蚀划分的基础上，结合土地利用现状图并参考水文、气象、土壤、植被和DEM（1∶50 000）等图件资料，利用Arc GIS 9．2的空间分析模块，对相关图件进行插值和叠加分析等操作来计算RUSLE方程的各个因子。最后利用修正后的通用土壤流失方程RUSLE来计算不同区域的土壤侵蚀模数。将所得结果参照土壤水蚀侵蚀强度分级标准［14］，进行重分类得到水蚀不同侵蚀等级的分布图。由于现行的风蚀研究指标分级标准主要针对沙地制定，而本研究区内的地貌类型主要是戈壁，因此结合研究区的实际情况制定了该区的风蚀研究指标体系（表1）。

表1 研究区风蚀强度分级

首先对T M影像进行必要的辐射校正和增强处理，然后结合DEM、土壤图、植被图及野外调查的结果，采用监督分类和目视解译相结合的方法，结合地形地貌、土壤、植被、气象资料等数据，对土壤侵蚀和土地利用／覆盖等进行解译，得到研究区土地利用和土壤侵蚀数据。利用Arc GIS 9．2的空间分析模块，将土壤数据、土地利用数据和土壤侵蚀数据进行空间叠加分析，进而对研究区土壤类型、土地利用类型与土壤侵蚀的耦合关系进行研究。

3 结果与分析

3．1 研究区土壤侵蚀现状分析

由表2中可以看出，从侵蚀形态看，研究区兼有风力侵蚀和水力侵蚀两种侵蚀形态，以风力侵蚀为主。风蚀总面积达68 422．04 h m2，占研究区面积的90%以上；水蚀总面积仅为4 780．82 h m2，主要分布在河流、居民点、工矿附近和山前倾斜洪冲积平原区。从侵蚀程度来看，研究区60%以上的区域侵蚀程度达到中度以上，表明该区土壤侵蚀程度较为严重。值得注意的是，虽然该区环境条件恶劣，植被盖度很低，但仍然有近30%的区域为轻度侵蚀，这主要是由于土壤表层石砾含量较高，而表面又多为大的砾石所覆盖，使得这些区域的土壤侵蚀模数相对较低；但从另外一个角度来看，由于土壤中细粒物质含量较低，一旦这些物质被侵蚀掉，土地生产力就会降至极低的水平甚至接近于零，整个生态系统就会处于崩溃的边缘。因此该区域土壤侵蚀程度的变化对整个研究区乃至其下游的绿洲都有着十分重要的影响。根据青海省生态功能区划，研究区被划分为沙漠化控制区［15］，因此，该区的土壤侵蚀必须给予足够重视，尽量避免人类活动对地表植被的破坏。

3．2 土壤类型与土壤侵蚀的耦合关系分析

由表2可以看出，土壤组成为灰棕漠土、高山漠土和洪积盐土的地区风蚀强烈。这些土壤的粘结性差，土层薄，多为砂砾质，土壤贫瘠，生物积累弱，植被盖度低，极易发生风蚀。研究区大部分位于鱼卡盆地内，盆地略呈三角形，由西北向东南方向延伸，该区盛行的西北风常年对鱼卡盆地进行吹蚀，加剧了该区的风蚀。另外，鱼卡盆地土壤母质多为第三纪和第四纪松散的堆积物，土壤类型多为洪积盐土，该土类表面盐分聚性强，上层积累大量易溶盐，导致该区植被稀疏。上述原因综合作用下，使得上述地类的风蚀较为严重。

表2 土壤类型与土壤侵蚀的关系h m2

表3 不同土地利用类型与土壤侵蚀的关系 h m2

土壤组成为灰棕漠土、高山漠土、高山荒漠化草原土的区域水蚀强烈。一方面，这些土壤的母质多为坡积物、湖积物、洪冲积物，土壤肥力低，植被稀疏，土质多为粗砾碎屑物质或砂砾物质；另一方面，研究区的主要水源为冰川融水，其自东向西横穿研究区。两者共同作用，加剧了分布在河流附近的灰棕漠土、高山漠土和高山荒漠化草原土的土壤侵蚀程度。特别是广泛分布于柴达木山的高山荒漠化草原土，其土壤侵蚀最为严重。高山荒漠化草原土主要分布在柴达木山上，该山体的干旱剥蚀明显，山麓遍布碎石，岩石裸露，且多沙化。长期水力侵蚀作用使得使该区的土壤侵蚀严重。

3．3 土地利用与土壤侵蚀的耦合关系分析

由表3可以看出，石漠、沙地、戈壁和砾漠风蚀强烈，其中极强度和剧烈侵蚀类型主要分布在沙地和戈壁两种土地利用类型上。这些区域风蚀严重的主要原因是该区植被稀疏、地表多被沙物质覆盖，为风蚀提供了丰富的沙源。同时这些土地利用类型的植被稀疏，地表粗糙度低，起沙风速低。因此，这些土地利用类型更易发生风蚀［16］。

沼泽、荒草地、盐碱地和沟道等地类的水蚀强烈，其中极强度侵蚀类型主要位于沼泽和沟道两种土地利用类型上。沼泽的侵蚀以水蚀为主，这主要是由于该区水资源条件相对于其它地类要好，水蚀强度要大于其它地类。由于沟道多为雨水冲刷形成，该区雨水为短历时降雨，短时间内会形成季节性河流，进而对沟道进行冲刷，使得沟道的侵蚀强度很大。

4 结 论

通过上述分析可知，研究区土壤侵蚀严重，兼有水蚀和风蚀两种侵蚀形态，侵蚀类型以风蚀为主。从土壤类型来看，灰棕漠土、高山漠土侵蚀强烈，兼有风力侵蚀和水力侵蚀。洪积盐土风蚀强烈，高山荒漠化草原土水蚀较为强烈；从土地利用类型来看，极强度风蚀和剧烈风蚀主要分布在沙地和戈壁两种土地利用类型上，极强度水蚀主要位于沼泽和沟道两种土地利用类型上。

鱼卡盆地是洪积盐土的主要分布区，该区土地利用类型以沙地和戈壁为主。独特的地理条件和地形特点使得鱼卡盆地风蚀严重。在今后的区域生态整治中，对人为活动剧烈且侵蚀最严重的盆地应特别保护，限制人为活动的范围和程度。特别是在人类活动（如矿产资源的开发）过程中，应充分考虑到该区的生态环境特点，采取有效措施防治水土流失，实现区域环境和经济的可持续发展。

［1］ 王思远，刘纪远，张增祥，等．不同土地利用背景下土壤侵蚀空间分布规律研究［J］．水土保持学报，2001，15（3）：48－51．

［2］ 范丽丽，沈珍瑶，刘瑞民．基于GIS的大宁河流域土壤侵蚀评价及其空间特征研究［J］．北京师范大学学报：自然科学版，2007，43（5）：563－566．

［3］ 邹亚荣，张增祥，周全斌．农牧交错带土地利用的土壤侵蚀状况分析［J］．水土保持通报，2004，24（5）：35－38．

［4］ 赵文武，傅伯杰，吕一河，等．多尺度土地利用与土壤侵蚀［J］．地理科学进展，2006，25（1）：24－33．

［5］ Russell S H，William W D．Landscape Erosion and Evolution Modeling［M］．New Yor k：Kluwer Academic／Plenu m Publishers，2001．

［6］ Stanley S W，Pierre C．U．S．Soil erosion rates－myth and reality［J］．Science，2000，289：248－250．

［7］ 周平，蒙吉军．鄂尔多斯市1988－2000年土壤水力侵蚀与土地利用时空变化关系［J］．自然资源学报，2009，29（10）：1706－1716．

［8］ 张鲁，周跃，张丽彤．国内外土地利用与土壤侵蚀关系的研究现状与展望［J］．水土保持研究，2008，15（3）：43－48．

［9］ 仲晓雷，郭成久，范昊明，等．辽宁省不同土地利用背景下土壤侵蚀空间分布规律［J］．水土保持研究，2008，15（3）：217－219．

［10］ 郝慧梅，任志远．区域LUCC的土壤侵蚀响应研究：以榆林市为例［J］．干旱区研究，2008，25（4）：583－591．

［11］ 邸利，张仁陟，张富，等．基于RS与GIS的定西市安定区土地利用变化与土壤侵蚀研究［J］．干旱区资源与环境，2011，25（2）：40－45．

［12］ 何兴元，胡志斌，李月辉．GIS支持下岷江上游土壤侵蚀动态研究［J］．应用生态学报，2005，16（12）：2271－2278．

［13］ 张洪江．土壤侵蚀原理［M］．北京：中国林业出版社，2000．

［14］ 中华人民共和国水利部．土壤侵蚀分类分级标准SL190－96（B）［S］．

［15］ 青海省环境科学研究设计院．青海省生态功能区划报告［R］．2003．

［16］ 刘明义，宋志超，张黎辉，等．沙地风蚀动力因子分析［J］．中国水土保持，2000（7）：28－30．