基于CSLE模型的陕北纸坊沟流域土壤侵蚀评价

马亚亚, 王 杰, 张 超, 刘国彬, 庞国伟

(1.西北农林科技大学 林学院, 陕西 杨凌 712100; 2.西北农林科技大学 水土保持研究所, 陕西 杨凌 712100;3.中国科学院 水利部 水土保持研究, 陕西 杨凌 712100; 4.西北大学 城市与环境学院, 陕西 西安710069)

黄土高原脆弱的生态环境以及人类的社会经济活动导致该区域水土流失日益严重。剧烈的土壤侵蚀使得土壤的表土层随径流流走，土壤中的营养物质大量流失，导致土壤肥力下降，土壤利用率降低以及生态环境恶化[1]。据统计，黄土高原的总面积为6.40×105km2，土壤侵蚀面积达3.90×105km2，土壤侵蚀模数大于等于15 000 t/(km2·a)的剧烈水蚀面积3.67×104km2，占全国同类面积的89%[2]。当前，对土壤侵蚀的研究已经引起了农学、土壤学、水文学等多学科的广泛关注[3]。因此，掌握土壤流失量的动态变化，这对优化配置水土保持措施、保护水土资源具有重要意义。土壤侵蚀模型是检测和预报土壤流失的重要工具。美国的Wischmeier等[4]于1965年建立的通用土壤流失方程USLE是最早建立的土壤侵蚀模型；1993年，Renard等[5]通过对通用土壤流失方程的修改而形成的新一代土壤侵蚀预测模型—修正通用土壤流失方程RUSLE。20世纪80年代以来，我国的专家学者以USLE的基本形式为依据，结合我国当地的地形特征，建立了不同的土壤侵蚀经验模型，其中比较有代表性的模型是江忠善[6]的坡面水蚀模型、范瑞瑜[7]建立的黄河中游地区小流域土壤侵蚀预报模型。此外，具有物理成因的土壤侵蚀模型也相继出现，Knisel等[8]建立的CREAMS模型；美国筹建水蚀预报项目WEPP[9]；澳大利亚Misra等[10]提出的GUEST等，但是由于物理模型自身的局限性，其实用性较低。因此，本文以刘宝元利用美国的通用土壤流失方程式(USLE)和修正通用土壤流失方程式(RUSLE)的有关思想，结合我国的地形特征提出的中国土壤流失方程(CSLE)为研究模型，该模型充分考虑了生物措施、工程措施和耕作措施对土壤侵蚀和水土流失过程和结果的影响，使其更适合中国实际。本研究采用中国水土流失方程(CSLE)作为研究模型，坡度坡长因子、降雨侵蚀力因子、土壤可蚀性因子、生物措施因子、耕作措施因子作为土壤侵蚀的5大因子，利用ArcGIS技术代入影响土壤侵蚀的各因子值，得到土壤侵蚀图，从而探讨土壤侵蚀强度与影响因子之间的空间变化的关系，以期为纸坊沟流域不合理的土地利用方式的调整和优化提供科学依据与理论支撑。

1 材料与方法

1.1 研究区概况

纸坊沟流域隶属于陕西省安塞县沿河湾镇，地理位置在东经109°13′46″—109°16′3″，北纬36°42′42″—36°46′28″，是延河支流杏子河下游的一级支沟，流域呈南北向狭长形，面积8.27 km2，海拔1 425.71～1 041.5 m。年平均降水量为549.1 mm，降水年际变化率较大，降雨量多集中在7—9月，且3个月的降水占到全年降水量的61.1%。流域地形破碎，土壤类型以黄绵土为主，其面积占土壤总面积的63.5%；且沟壑密度8.06 km/km2，水土流失严重，多年平均侵蚀模数达14 000 t/km2。流域植被区划属于暖温带森林草原区，受人为活动的影响，天然森林已经绝迹，取而代之的天然植被，多为干旱草本植物。由于流域集水面积小，地面坡度陡，土质疏松，因此在降雨情况下，特别是暴雨后，产流集中，挟沙量大；短历时高强度的暴雨是造成严重水土流失的主要原因[11]。

1.2 研究方法

1.2.1 CSLE模型及其参数计算

(1) CSLE模型[12]。研究区内主要以水蚀为主，对土壤侵蚀量的计算采用的是中国土壤流失方程CSLE(Chinese soil Loss equation)，其计算表达式为：

A=R·K·L·S·B·E·T

(1)

式中：A——年均土壤流失量〔t/(hm2·a)〕；R——研究区降雨和径流因子〔(MJ·mm)/(hm2·h·a)〕；K——土壤可侵蚀因子〔t·h/(MJ·mm)〕；L——坡长因子；S——坡度因子；B——生物措施因子；E——工程措施因子；T——耕作措施因子。

(2) 降雨侵蚀力因子[13]。计算公式为：

(2)

式中：Pi——多年的月均降水量(mm)；P——多年的年降水量(mm)。

采用Wischmeier[13]提出的经验公式，根据安塞站及周边各气象站点观测的气象数据，得到降雨侵蚀力因子R值，然后采用反距离加权内插法生成研究区降雨侵蚀力R值图。

(3) 土壤可蚀性因子[14-15]。计算公式为：

(3)

(4)

式中：SAN——砂粒含量(%)；CLA——粉粒含量(%)；SIL——黏粒含量(%)；C——有机碳含量。

张科利等[16]通过在全国建立6个观测点(其中包括西北黄土高原的陕西安塞)来研究了我国不同水土流失区的土壤可蚀性值问题。结果表明,国外关于土壤可侵蚀性的研究成果不能直接应用于我国的研究中。虽然估算值明显大于实测值，但估算值与实测值之间存在有良好的线形关系。为此提出了针对国外模型的订正模型，其计算表达式为：

K=-0.013 86+0.515 75

(5)

(4) 坡度坡长因子。采用张宏鸣根据CSLE的计算公式编写的AML程序提取LS因子，同时对DEM值中的NoData值也进行了填充[12]，坡度因子的算法是对McCool等[17]和刘宝元等[18]提出的坡度算法的改良，当坡度β≤5°时，采用等McCool 提出的坡度因子算法，当坡度β>5°时，采用的是刘宝元的坡度因子计算方法；坡长因子的计算，采用的是Wischmeier和Smith[4]提出的经验公式。计算公式为：

(6)

L=(λ/22.13)

(7)

α=β/(β+1)

(8)

β=(sinθ/0.089 6)/3.0(sinθ)0.8+0.56

(9)

式中：S——坡度因子；θ——从DEM提取的坡度值(°)；L——坡长因子；λ——DEM提取的坡长值； 22.13——标准小区的坡长；α——坡长因子指数；β——细沟侵蚀和面蚀的比值。

(5) 水保措施因子(BET)。CSLE模型将水保措施因子区分为生物措施因子(B)、工程措施因子(E)以及耕作措施因子(T)。工程措施因子(T)由于缺少基础数据，本研究中暂不考虑其对土壤侵蚀的影响，将其按照常量对待。生物措施因子(B)和耕作措施因子(T)均参考谢红霞等[16]对黄土高原土壤侵蚀的研究，利用研究区土地利用类型图和基于归一化植被指数求得的植被盖度图，对陕北纸坊沟流域的不同土地利用类型和不同植被盖度赋予不同的B值，利用ArcGIS技术得到研究区B值分布图；对耕作措施因子(T)赋值，没有耕作措施的地方赋值为1，然后得到了研究区T值分布图。

1.2.2 分析土壤侵蚀影响因子空间分布特征 土壤固有结构和地貌特征决定了土壤侵蚀力因子和坡度坡长因子，因此在CSLE中没有必要考虑时间序列上的差异。但降雨侵蚀力和水保措施因子是随着时间的推移而变化，它可以作为土壤侵蚀研究动态变化的关键因素[19]。在区域土壤侵蚀模型中，坡度坡长因子是对土壤侵蚀最为敏感的因子，生物措施因子是区域土壤侵蚀中第二敏感因子，它们随着DEM分辨率的增大，土壤侵蚀模数都会变大。降雨侵蚀力因子和土壤可蚀性因子是土壤侵蚀模型中不敏感因子，改变栅格的大小，土壤侵蚀模数不会发生大的变化。但降雨侵蚀力因子随着时间变化，对土壤侵蚀具有显著的影响。

1.2.3 分析土壤侵蚀影响因子空间分布规律 通过前面对中国水土流失方程各因子的计算，确定了小流域土壤侵蚀所必需的各因子的数据，利用ArcGIC 10.0将上述各因子值的图层叠加运算，得到各计算单元的土壤侵蚀量。根据土壤侵蚀的试验结果，参考中国水利部2007年制定的《土壤侵蚀分类分级标准》[20]，将研究区的土壤侵蚀强度等级划分为微度侵蚀、轻度侵蚀、重度侵蚀、强烈侵蚀、极强烈侵蚀和剧烈侵蚀6类，制作出研究区土壤侵蚀强度等级分布图，分析土壤侵蚀强度分级的分布规律。

2 结果与分析

2.1 土壤侵蚀影响因子空间分布特征

2.1.1 降雨侵蚀力因子 降雨侵蚀力因子反映降雨引起土壤分离和搬运的动力大小，即降雨产生土壤侵蚀的潜在能力[21]。降雨侵蚀力因子R的大小主要取决于研究区降雨强度和持续时间，暴雨越大，持续时间越长，侵蚀越容易发生。由图1可知，从整体来看，2005—2016年流域年降雨量的起伏变化波动不大，但是2013年间降雨量最大。11 a期间的平均降雨量为476.6 mm，波动范围在±403 mm之间，每月最大降雨量与年降水量变化趋势大致相同，2016年的降水俩最低，最高值为2013年的降雨量，虽然各个年份之间的降雨量起伏波动，但总体上呈现出下降趋势。大量研究结果表明：降雨侵蚀力的变化趋势与降雨量变化趋势基本一致，说明研究区内各个年份的降雨量逐渐减少，而降雨侵蚀力也就越来越弱。

图1 纸坊沟流域降雨特征

2.1.2 土壤可蚀性因子 土壤可蚀性因子定义为标准小区内单位降雨侵蚀力引起的土壤流失率，是土壤对雨滴击溅或地表径流等侵蚀介质剥蚀和搬运的敏感程度，反映了土壤颗粒被降雨和径流破坏、输移的敏感性，也是土壤抵抗由降雨、径流产生的侵蚀能力的综合体现[22]。纸坊沟小流域内主要土壤类型为黄绵土，占土壤总面积的63.4%左右；其次为红胶土，其面积占该流域总面积的25.1%，其余为零星分布的潮土、新积土、粗骨土和石质土等。本研究以纸坊沟流域土壤类型图为基础，参考李斌兵等[21]对纸坊沟流域土壤侵蚀强度空间分布的研究成果，获得各土壤类型对应的K值单位，在ArcGIC中将K值作为属性字段添加到土壤类型图中，生成30 m栅格大小的K值栅格图。纸坊沟流域内黄绵土、红胶土、潮土、新积土、粗骨土、石质土的K值分别为：0.37，0.42，0.34，0.22，0.23和0.24 t·hm2·h/(hm2·MJ·mm)。由于黄绵土土质疏松，质地均一，所以抗侵蚀性差。

2.1.3 坡长坡度因子 坡长和坡度因子反映了地形地貌对土壤侵蚀的影响，通常为侵蚀动力的加速因子[23]。坡度和坡长因子是土壤侵蚀研究中最重要的地形因子之一，也是决定径流冲刷能力的主要因素。在其他条件相同的情况下，坡度越陡，坡长越长，土壤发生分离、滑落的潜在危险性就越大，土壤侵蚀就越来越严重。

本研究采用流域的5 m分辨率的DEM数据对研究区坡度因子进行提取。首先在ArcGIC中用Spatial Analyst模块中的Surface Analyst功能得到坡度栅格图，然后在ArcMAP中，利用3D Analyst模块中的Reclasify模块将坡度栅格图分成≤5°，5°～8°，8°～15°，15°～25°，25°～35°，≥35°共6个坡度级，最后对于生成的坡度分级结果，打开其属性表，利用字段计算器求得各个坡度级别的面积(表1)。

研究区坡度坡长因子的空间分布特征与海拔分布具有显著的一致性。从坡度分级结果来看，该流域坡度≤5°的面积为0.53 km2，占全流域面积的6.46%；坡度为5°～8°的面积为0.84 km2，占全流域10.21%；坡度为8°～15°的面积为1.28 km2，占全流域15.53%；坡度为15°～25°的面积为1.46 km2，占全流域17.66%；坡度为25°～35°的面积为1.83 km2，占全流域23.34%；坡度为≥35°的面积为2.21 km2，占全流域面积的26.80%。从流域的坡度分级图可以看出研究区内坡度比较大，全流域坡度大于15°的面积占全流域的67.80%，是土壤侵蚀的主要来源区域。

表1 纸坊沟流域坡度分级统计

坡度坡长因子反映了地形地貌对土壤侵蚀的影响，通常为侵蚀动力的加速因子。在其他条件相同的情况下，坡度越来越陡，坡长越来越长，产流形成的流量越大，土壤侵蚀越来越严重。坡度坡长因子可以利用LP-Tool[14]软件进行计算，然后将坡度因子栅格分布图(图2)与坡长因子栅格分布图(图3)进行相乘叠加，得到该流域的LS值分布图(图4)。流域内5 m分辨率的LS因子平均值为8.79，最大值为138。小于5的坡度坡长因子值的面积占流域面积31.41%，坡度坡长因子值位于5，10的占流域面积的30.89%，19.30%的流域面积的坡长因子值在10～15，坡度坡长因子值在26以下约占流域面积的97.46%。说明大于CSLE模型中规定的LS最大值72的地形因子值很少，5 m分辨率计算出的坡度坡长因子值与实际情况吻合度极高。

图3 纸坊沟流域坡长因子值分布

图4 纸坊沟流域LS因子值分布

2.1.4 水保措施因子 生物措施因子(B)定义为在一定地表覆盖和管理措施下土壤流失量与同等条件下适时翻耕、连续休闲对照地上土壤流失量之比[24]。植被覆盖通常被认为是土壤侵蚀动力的抑制因子，其值小于或等于1。植被覆盖纸坊沟流域主要农作物为玉米、土豆、谷子和大豆，据张岩等[21]研究，黄土高原玉米、土豆、谷子、大豆的C因子分别为0.28，0.47，0.53和0.51，由于缺乏具体农作物信息，取其均值0.45作为研究区耕地的B因子值，水域的B因

子赋值为1，林地为0.05，草地为0.09。利用纸坊沟流域的土地利用图层，将B因子值赋予相应的土地利用类型，得到研究区植被覆盖B因子图层(图5)。从图中可以看出生物措施因子B值接近于1的植被稀疏区主要分布在沟口和流域左岸，例如拐沟、大范家沟、大罗锅沟、正沟。

图5 纸坊沟流域B因子值分布

在CSLE模型中，耕作措施因子定义为采取了一定的耕作措施后，土壤流失量相对于顺坡种植时土壤流失量的比例[25]。农业生产中通常利用等高耕作、等高带状种植、梯田等水土保护措施，调整水流形态、斜坡坡度和表面流的汇流方向，减少径流量，降低径流速率，从而减轻土壤侵蚀。T值通常介于0～1之间，不会发生土壤侵蚀的地区赋值为0；1值被赋予无任何水土保持措施的土地类型。土壤耕作措施因子T被认为是CSLE 模型中最难确定的因子，因此参考谢红霞提出的耕作措施因子的计算方法得到纸坊沟流域的耕作措施因子值(表2)，并借助ArcGIS中的栅格计算器将T因子值赋予相应的坡度分级图中，得到水土保持耕作措施T因子值图。在图中，非耕地所占用的面积占总面积的56.29%，坡度为≤5°的耕地占总面积的16.29%，5°～10°坡度的耕地面积占总面积的25.48%，10°～15°的耕地占总面积的1.94%。

表2 不同坡度下耕作措施因子值

2.2 小流域土壤侵蚀强度的空间分布

通过统计各个侵蚀强度的面积及其所占总面积之比的数据(表3)，从中可以看出，侵蚀强度占比例最大的一类为中度侵蚀，接近总面积的50%左右，其次3类比较严重的侵蚀之和占到了总面积的45%以上，微度侵蚀和轻度侵蚀各占总面积的1.29%和3.97%，从总体来看纸坊沟流域土壤侵蚀仍比较严重，须采取必要的治理措施。从图6—7可看出，流域土壤侵蚀分布特征为沟口区易发生极强烈和剧烈侵蚀，土壤流失比较严重；流域的中部以中度侵蚀为主；沟头侵蚀强度为轻度和中度侵蚀。不同坡度带上的土壤侵蚀强度差异较大，流域中部坡度较缓，坡度大多集中在15°以上，土壤侵蚀以轻度和中度为主，并且微度与轻度的土壤侵蚀主要发生在林地或水域，说明该处是植被覆盖和水土保持比较好的地块。而流域沟头坡度较陡，坡度大多集中在15°～35°，且植被覆盖和水土保持措施较差，所以土壤侵蚀明显以中度为主，并且以主沟为界,流域左岸的侵蚀强度明显的高于右岸，说明流域的侵蚀强度呈明显不对称性。不同土地类型上土壤侵蚀强度的差异性，主要是由于各类土地的土壤理化性质和所处谷坡坡度造成的[26]。在各土地利用类型中，侵蚀面积最小的是水域；其次是林地和草地侵蚀面积较小，由此可知，侵蚀强度较大的区域主要分布在耕地。主要由于长期不合理的耕作方式，使地表土质疏松，导致土壤蓄水和透水能力变差，从而加重土壤流失。

图6 土壤侵蚀等级

图7 研究区土地利用状况

侵蚀强度级别面积/hm2占面积的百分比/%微度1 0751.29轻度3 3003.97中度40 60048.78强烈22 47527.00极强烈11 37513.67剧烈4 4005.29

3 讨 论

3.1 土地利用与土壤侵蚀关系的讨论

纸坊沟流域内土地类型以林地、耕地、草地为主。耕地是土壤侵蚀最易发生的地类，坡地开垦和坡地耕作是造成土壤流失的主要原因。耕地种植区实施退耕还林后，耕地面积不断减少，过度开垦状况有所好转，生态环境不断改善，土壤侵蚀的面积与强度整体呈现改善趋势。在林地动态区，林地面积持续增加，但由于退耕还林初期，大量育林育草，对地表扰动大，林木的覆盖度低，土壤侵蚀强度不降反增，但随着林木覆盖度的不断增加，土壤侵蚀强度慢慢降低。在草地动态区，草地主要向林地转移。由于退耕还草工程的实施，使得部分耕地退化为草地，部分草地来自于林地，土壤侵蚀强度逐渐变小。耕地、林地和草地的不断变化，说明土地利用变化与土壤侵蚀强度具有密切联系，表现为耕地的土壤侵蚀强度较强，林地和草地侵蚀强度相对较弱，说明增加林地和草地，减少耕地，能够显著减弱土壤侵蚀。

3.2 小流域土壤侵蚀防治措施

纸坊沟流域位于黄土丘陵沟壑区，地块狭小分散、沟壑纵横、土质松软、植被稀疏，土壤侵蚀严重；不合理的土地利用方式，脆弱的生态环境，水土流失严重制约了地方经济的可持续发展。为了改善当地生态环境，实现区域经济的可持续发展，按照新时期水土保持规划内容，分析纸坊沟流域内的土壤侵蚀强度，以高效合理的利用土地资源，建设人地协调的绿色生态系统为目标，通过合理安排土地资源的利用方式，确保粮食安全；同时根据纸坊沟流域的自然条件实施造林种草，确保生态安全；贯彻执行国家和地方有关生态文明建设的政策和规定，既要绿水青山，也要金山银山，实现社会、经济、生态效益的并行发展，但是要把生态环境保护放在更加突出的位置。

由坡度分级图和土壤侵蚀强度图分析可知,土壤侵蚀与土地的坡度和植被覆盖关系密切。坡度为15°～25°的区域土壤侵蚀状况严重，一方面由于15°～25°地区范围内遭受人为开垦，另一方面由于15°以上地区以草地和稀疏植被为主，水土保持作用低。因此在纸坊沟流域水土保持规划设计中，应加强对15°～25°坡度地区重点防治，布设相应的措施。纸坊沟流域内地形破碎，植被稀疏，中度以上侵蚀面积比例超过50%，稀疏植被和裸地缺少地表覆盖、草地植被退化和耕地垦殖是土壤侵蚀相对严重的主要原因[27]。如表4所示，5°～20°之间，土壤侵蚀以轻度和中度为主，着重梯田工程措施，提高单位土地面积的生态效益和经济效益；对于土壤侵蚀严重的陡坡地实施有计划的封禁保护，抚育措施，减少人为干预，使其自然恢复；没有修建梯田的坡耕地，可在坡地上种植草地、灌木，实行草粮或草灌带状间作，缩小坡面长度，减缓降雨对地面的冲刷，提高水土保持功能。

研究区土壤侵蚀强度与土地利用和植被覆盖表现出很大的相关性，土地资源的合理利用，是水土保持的一项重要措施和根本途径[28]。因此在纸坊沟流域水土保持规划设计中，合理的安排农、林、牧用地，增加植被覆盖度，例如在农牧交错带，实行以综合防治为主，工程与治理相结合；以保护和促进生态自然修复为主，坚持生态保护与治理相结合[29]，为区域经济的发展提供优良的生态环境保证。在坡度平缓的地区建设高标准的良田，利用改垄、梯田等措施改造坡耕地，提高土壤蓄水保土能力；坡中部利用等高耕作和等高带状间作等工程措施，拦蓄地表径流，减少水土流失；坡上部全部退耕，还林还草，营造生态防护林。林草植被建设是水土保持综合治理、改善黄土高原的生态环境的重要措施之一。实施造林种草，应遵循适地适树，宜林则林，宜草则草的根本性原则。

表4 小流域坡面水土保持技术措施

4 结 论

本研究在基于ArcGIS 技术上采用CSLE模型对纸坊沟流域影响土壤侵蚀的相关因子进行建模和可视化处理，估算了单元土壤侵蚀量，结合水利部水力侵蚀强度分级，分析了纸坊沟流域土壤侵蚀强度空间分布特征[30]。

(1) 流域土壤侵蚀分布特征为沟口区易发生极强烈和剧烈侵蚀，土壤流失比较严重；流域的中部以中度侵蚀为主；沟头侵蚀强度为轻度和中度侵蚀。侵蚀强度占比例最大的一类为中度侵蚀，接近总面积的50%左右。

(2) 同坡度带上的土壤侵蚀强度差异较大，流域中部坡度较缓，坡度大多集中在5°～15°，土壤侵蚀以轻度和中度为主，并且微度与轻度的土壤侵蚀主要发生在林地或水域。

(3) 不同类型土地上土壤侵蚀强度具有差异性，侵蚀面积最小的是水域；林地和草地侵蚀面积较小；耕地的侵蚀强度最大。