甘蔗作为中国重要的糖料作物，据统计，中国的甘蔗种植面积在全球排名第三位，约占6.05%。而中国70%的甘蔗种植在坡耕地上^[1]。广西壮族自治区(以下简称“广西”)是中国最重要的甘蔗生产基地，2020年，甘蔗种植面积和产量分别占全国约65%，70%^[2]。广西甘蔗种植区多为山地丘陵地貌，大部分为坡耕地且种植甘蔗，小部分坡度大的坡顶为林地且多种植桉树，其余土地利用类型较少。调查表明，甘蔗生长周期长，主要种植于土层较薄的赤红壤和红壤的旱作坡耕地上，岩性抗侵蚀性弱，由于对蔗糖需求增加，甘蔗种植扩张到坡度较大的山坡地，甘蔗生长前期植被覆盖度较小，加上降雨集中且多暴雨，而甘蔗为剑形叶片，大部分降雨穿透植株直接击打土壤表面，且顺坡种植等不合理的耕作方式，易引发严重的土壤侵蚀。但目前有关蔗区土壤侵蚀特征，侵蚀速率等方面数据不够明晰，有必要进一步探究，为蔗区土壤侵蚀防治提供依据。

目前，研究土壤侵蚀的方法有人工模拟降雨法、径流小区观测法、核素示踪法和侵蚀针法等。付兴涛等^[3]采用人工模拟降雨法，研究了南方红壤丘陵区作物覆盖坡耕地上不同雨强下坡长对其土壤侵蚀的影响。Hartanto等^[4]在印度尼西亚加里曼丹中部的森林集水区设置径流小区，发现影响土壤侵蚀的重要生态因子包括冠层盖度、幼树密度、凋落物深度和木本碎屑，但降雨对径流大小和土壤侵蚀量影响更大。Rabesiranana等^[5]运用放射性同位素¹³⁷Cs和²¹⁰Pb估算马达加斯加土壤侵蚀速率为18 t/(hm²·a)，比当地农业气候条件下的可持续土壤流失率高出两倍。赵珩钪等^[6]利用侵蚀针法，探究了黄土丘陵沟壑区小流域退耕坡地不同植被恢复条件下多种不同植物群落的土壤侵蚀特征。但这些方法需要投入大量资金、时间和人力，局限在较小范围，很难大规模推广使用。为此，许多评价侵蚀、水土流失量的模型与方程被提出来。经验模型有：通用土壤流失方程(universal soil loss equation，USLE)、修正的通用土壤流失方程、中国土壤流失方程等，物理模型有欧洲土壤侵蚀模型、水蚀预报项目、荷兰的土壤侵蚀模型、水土评价模型等^[7]。但是，物理模型所需数据较多，参数计算较复杂，而USLE/RUSLE模型机制较简单且数据易获取，在国内外广泛应用。Sinha等^[8]利用USLE估算印度阿杜拉河和马哈伦吉河耕地不同作物的土壤流失，发现土壤流失量已超过容许限度，作物因素即人类活动对土壤流失量更为重要。Lu等^[9]利用USLE，RS和GIS相结合的方法，探讨巴西亚马逊流域土壤侵蚀风险与土地利用、土地覆盖之间的关系，森林向其他土地覆盖转变是土壤侵蚀风险增加的重要因素。朱青等^[10]采用RUSLE模型和随机森林法定量分析了赣江上游流域土壤侵蚀及其影响因子。但目前有关甘蔗种植区土壤侵蚀的估算主要利用田间监测法和径流小区法。谢如林等^[11]利用田间监测法，研究了广西旱坡地甘蔗种植不同施肥模式下地表径流引起的土壤侵蚀及氮、磷流失特点。杨任翔等^[12]基于径流小区观测，探究南方高强度、高频次降雨下甘蔗种植模式对赤红壤坡面土壤侵蚀的影响。而利用GIS技术和土壤侵蚀模型等定量估算、分析蔗区土壤侵蚀特征、土壤侵蚀速率等鲜有报道。因此，本文将GIS技术与RUSLE模型相结合，选取广西典型赤红壤集约化蔗区那辣小流域，定量估算该流域的平均土壤侵蚀速率，确定土壤侵蚀等级和不同子流域的土壤侵蚀状况，探讨不同坡度、土地利用类型与土壤侵蚀强度的分布特征，以期为广西集约化蔗区小流域防治土壤侵蚀提供科学参考。

1 研究区概况

那辣小流域位于广西崇左市扶绥县境内，位于107°39′29″—104°40′17″E，22°20′36″—22°20′50″N，流域总面积为1.29 km² 〔划分为3个子流域，子流域1(SW₁)，子流域2(SW₂)和子流域3(SW₃)，其面积分别为0.22 km²，0.62 km²，0.45 km²〕，流域出口排水进入连接左江流域的客兰水库。流域为山地丘陵地貌，坡度0°~35°(表 1)，属南亚热带湿热气候区，年平均气温20~25 ℃，雨量充沛，全年降雨主要集中在5—9月，年降雨量在1 100 mm左右，土壤类型为赤红壤，第四纪红土母质，pH值4.7~5.2。流域坡耕地种植甘蔗，占土地利用总面积的82.85%，林地种植桉树，占土地利用总面积的10.99%，道路和沟渠分别占5.40%，0.77%(附图1，见封3)。该小流域2—3月种植甘蔗并在12月至次年3月收割，一般连续种植3 a再翻耕新种。

2 数据来源与研究方法 2.1 数据来源

本研究主要使用以下数据：①流域2020年的月降雨量和年降雨量，来源于流域内监测站收集的降雨数据。②流域土壤类型和属性数据来源于1∶100万全国土壤类型数据和中国土壤数据库(http://vdb3.soil.csdb.cn/)，流域实测的土壤理化性质；③数字高程模型DEM，约为3 m空间分辨率，来源于流域测绘所得数据；④流域2020年土地利用图，由无人机航拍影像图进行监督分类获得；⑤2020年流域植被覆盖度数据获取来源于使用无人机在甘蔗4个生长期拍摄流域的四景可见光正射影像图。

2.2 RUSLE模型各因子的确定

本文采用RUSLE模型^[13]来估算那辣小流域的土壤侵蚀速率。其公式为：

$ A=R \times K \times L \times S \times C \times P $

(1)

式中：A为土壤侵蚀速率〔t/(hm²·a)〕; R为降雨侵蚀力因子〔MJ·mm/(hm²·h·a)〕; K为土壤可蚀性因子〔t·hm²·h/(hm²·MJ·mm)〕; LS为地形因子(L为坡长因子，S为坡度因子)，无量纲；C为植被覆盖与作物管理因子，无量纲，范围0~1；P为水土保持措因子，无量纲，范围0~1。

(1) 降雨侵蚀力因子R。R因子代表降雨对土壤侵蚀的影响，本研究采用Wischmeier^[14]提出的简化算法，计算公式为：

$ R=\sum\limits_{i=1}^{12} 1.735 \times 10^{\left(1.5 \lg \left(\frac{p_{i} 2}{p}\right)-0.8188\right]} $

(2)

式中：R为降雨侵蚀力因子，国际单位为〔MJ·mm/(hm²·h·a)〕; P_i为月降水量(mm); P为年降水量(mm)。本研究采用研究区2020年的年降雨量和月降雨量数据代入公式计算出流域内年降雨侵蚀力图。

(2) 土壤可蚀性因子K。土壤可蚀性反映了土壤对侵蚀的敏感程度，这里采用Williams^[15]提出的侵蚀—生产力评价模型EPIC中土壤机械组成和有机碳含量的计算公式：

$ \begin{aligned} K=&\left\{0.2+0.3 \exp \left[-0.0256 \mathrm{SAN}\left(1-\frac{\mathrm{SIL}}{100}\right)\right]\right\} \times\left(\frac{\mathrm{SIL}}{\mathrm{CLA}+\mathrm{SIL}}\right)^{0.3} \times \\ &\left(1-\frac{0.25 C}{C+\exp (-2.95 C+3.72)}\right) \times\left(1-\frac{0.7 \mathrm{SN}_{1}}{\mathrm{SN}_{1}+\exp \left(-5.51+22.9 \mathrm{SN}_{1}\right)}\right) \end{aligned} $

(3)

式中：K为土壤可蚀性因子，为美制单位，乘以系数0.131 7转为国际单位〔t·hm²·h/(hm²·MJ·mm)〕; SAN为砂粒(0.050~2.000 mm)质量分数(%); SIL为粉粒(0.002~0.050 mm)质量分数(%); CLA为黏粒(＜0.002 mm)质量分数(%); C为土壤有机碳质量分数(%)；SN₁=1-SAN/100；通过1∶100万全国土壤类型数据，根据研究区所在扶绥县的土壤类型通过ArcGIS软件把流域的土壤类型裁剪出来，并根据中国土壤数据库找到研究区所在扶绥县的土种并结合野外调查和流域实测土壤理化性质，以确定研究区的土壤类型和土壤的砂粒、粉粒、黏粒和有机碳质量分数，代入公式(3)求出土壤可蚀性因子K值并赋值到土壤类型图中。

(3) 坡度坡长因子LS。坡度、坡长是影响土壤侵蚀的地形因子，本文的坡长因子L计算采用Desmet等提出的经典公式，坡度因子缓坡采用McCool的坡度公式，陡坡用刘宝元的坡度公式^[16]：

$ S= \begin{cases}10.8 \sin \theta+0.03 & \left(\theta<5^{\circ}\right) \\ 16.8 \sin \theta-0.50 & \left(5^{\circ} \leqslant \theta<10^{\circ}\right) \\ 21.9 \sin \theta-0.96 & \left(\theta \geqslant 10^{\circ}\right)\end{cases} $

(4)

$ L=\left(\frac{\lambda}{22.1}\right)^{n}, \quad n= \begin{cases}0.2 & \left(\theta \leqslant 1^{\circ}\right) \\ 0.3 & \left(1^{\circ}<\theta \leqslant 3^{\circ}\right) \\ 0.4 & \left(3^{\circ}<\theta \leqslant 5^{\circ}\right) \\ 0.5 & \left(\theta>5^{\circ}\right)\end{cases} $

(5)

式中：θ为坡度(°)；n为坡长指数；量纲为1；λ为坡长(m)。先将分辨率为3 m的DEM数据在ArcGIS中提取坡长和坡度，再代入公式(4)，(5)分别进行计算，得到坡长因子L和坡度因子S，相乘得到地形因子LS。

(4) 植被覆盖与管理因子C。C为植被覆盖与管理因子，在C因子的计算过程中需要以无人机在甘蔗4个生长期拍摄流域的四景可见光正射影像图为基础，借助ArcGIS和ENVI软件，根据过绿减过红指数ExG-ExR利用像元二分法模型算出流域的植被覆盖度^[17]，再根据蔡崇法等^[18]建立的覆盖度与C值的关系来计算C因子。

$ \mathrm{ExG}-\mathrm{ExR}=3 \times G-2.4 \times R-B $

(6)

式中：ExG-ExR为过绿减过红植被指数; R，G，B分别为红、绿、蓝3个波段的像元值，把流域的正射影像图利用ENVI软件的波段计算器计算。

$ \begin{aligned} f_{c}=&\left\{\left(b_{1}<b_{\min }\right) \times 0+\left(b_{1}>b_{\max }\right) \times 1+\right.\\ &\left(b_{\min }<b_{1}<b_{\max }\right) \times \\ & {\left.\left[\left(b_{1}-b_{\min }\right) /\left(b_{\max }-b_{\min }\right)\right]\right\} } \end{aligned} $

(7)

式中：f_c为植被覆盖度(%)；b₁是植被指数的值，本研究用过绿减过红植被指数的值；b_min，b_max分别是植被指数的最小值和最大值，利用ENVI软件，根据过绿减过红植被指数的统计直方图，确定纯植被和纯裸地的值，选取植被指数频率累计5%为最小值(b_min)，植被指数频率累计95%为最大值(b_max)。

$ \left\{\begin{array}{cl} C=1 & \left(f_{c}=0\right) \\ C=0.6508-0.3436 \lg f_{c} & \left(0<f_{c}<78.3 \%\right) \\ C=0 & \left(f_{c} \geqslant 78.3 \%\right) \end{array}\right. $

(8)

式中：C为研究区的植被覆盖与管理因子。

(5) 水土保持措施因子P。水土保持措施因子P，本文通过调查那辣小流域水土保持现状并根据相关研究报道^[19]，对每一种土地利用类型赋予P值见表 2。

3 结果与分析 3.1 土壤侵蚀空间分布特征

根据土壤侵蚀模型中各因子的运算关系，在ArcGIS 10中将获取的各土壤侵蚀影响因子栅格图均转化为同一坐标下像元大小为0.5 m×0.5 m后进行叠加计算，得到各像元土壤侵蚀量，并根据水利部制定的《土壤侵蚀分类分级标准》(SL190-2007)进行等级划分^[20]，获得那辣小流域土壤侵蚀强度分级图。由封3附图2和表 3可知，2020年，那辣小流域平均土壤侵蚀速率为22.97 t/(hm²·a)，属轻度侵蚀水平，是水利部规定的南方红壤丘陵区土壤允许流失量5 t/(hm²·a)的4.6倍。轻度侵蚀等级及以下的面积占比达78.09%，中度侵蚀等级及以上的侵蚀面积占21.91%，说明从总体上看流域土壤侵蚀程度较轻。流域内中度侵蚀等级及以上的侵蚀量占侵蚀总量的81.39%，侵蚀量为轻度侵蚀等级及以下的4.37倍，表明流域局部区域土壤侵蚀比较严重，因此应该重视防治那辣小流域的土壤侵蚀。

由表 4可知，3个子流域的平均土壤侵蚀速率由大到小为SW₃＞SW₂＞SW₁，子流域SW₃的平均土壤侵蚀速率为25.31 (t/hm²·a)，是SW₁的1.67倍，SW₂与SW₃相差不大。然而子流域SW₂在3个子流域中的侵蚀面积和侵蚀量占比最大，分别为47.90%，50.25%，土壤侵蚀量分别是子流域SW₁，SW₃的4.42，1.31倍。

3.2 土壤侵蚀与坡度的关系

坡度是引发土壤侵蚀的众多影响因素之一，因此探究不同坡度下土壤侵蚀的变化规律和特征是防治那辣小流域水土流失的关键举措。根据水利部颁布的关于土壤侵蚀坡度等级划分标准(SL190-2007)，将流域坡度按照0°~5°，5°~8°，8°~15°，15°~25°，25°~35°进行重分类，将分类后的坡度等级图与土壤侵蚀速率图进行空间叠加分析，得出不同坡度级别上的土壤侵蚀分布状况。从表 5中可知，随着坡度的增加，平均土壤侵蚀速率不断增大，其值在5.64~112.64 t/(hm²·a)之间。在那辣小流域中，8°~15°坡度级别的侵蚀面积和侵蚀量最大，其值分别为49.91 hm²，1 420.39 t/a，侵蚀强度属于中度侵蚀；25°~35°的侵蚀面积和侵蚀量最小，其值分别为0.37 hm²，41.58 t/a，侵蚀强度属于极强烈侵蚀。不同侵蚀强度的土壤侵蚀量在各坡度级别所占比例差异较大，随着坡度的增加，微度侵蚀在各坡度级别所占侵蚀量比例逐渐减小，而其他侵蚀强度，以15°坡度级别为中心向上和向下土壤侵蚀量在各坡度级别所占比例呈减小趋势。流域内土壤侵蚀主要发生在5°~25°的坡度范围内，占侵蚀总量的92.36%，是预防和治理土壤侵蚀的重点区域，应采取合理有效的防治措施。由图 1可知，不同侵蚀强度在各个坡度级别的侵蚀面积所占比例差异较大。总体上，微度侵蚀在各个坡度级别上所占的面积比例总体呈减少趋势，而其他侵蚀强度所占的面积比例则呈现出先增加后减少趋势。说明坡度分级标准可以明显区分不同土壤侵蚀等级与坡度的变化关系。

3.3 土壤侵蚀与土地利用的关系

将流域土壤侵蚀强度图与土地利用类型图在ArcGIS中进行叠加分析，得到不同土地利用类型的土壤侵蚀强度分布情况。由表 6可知，坡耕地(甘蔗)、林地(桉树)、道路3种土地利用类型的土壤侵蚀强度面积比例主要处于微度和轻度侵蚀，而沟渠处于微度侵蚀。

从表 7可知，在不同土地利用类型中，侵蚀面积最大的是坡耕地(甘蔗)，占82.85%，林地(桉树)、沟渠和道路的侵蚀面积共占17.15%，流域平均侵蚀速率由大到小为：林地(桉树)＞道路＞坡耕地(甘蔗)＞沟渠，林地(桉树)平均土壤侵蚀速率为53.59 t/(hm²·a)，分别是坡耕地(甘蔗)、道路、沟渠的2.84，2.12，27.91倍，而侵蚀量最大的是坡耕地(甘蔗)，其值为2 013.05 t/a，分别是林地(桉树)、道路、沟渠的2.65，11.46，1 059.50倍。

4 讨论与结论 4.1 讨论

本文利用RUSLE模型计算2020年广西典型赤红壤集约化蔗区那辣小流域平均土壤侵蚀速率为22.97 t/(hm²·a)，属轻度侵蚀水平，是水利部规定的南方红壤丘陵区土壤允许流失量5 t/(hm²·a)的4.6倍。轻度侵蚀等级及以下的面积占比达78.09%，从总体上看流域土壤侵蚀程度较轻，中度侵蚀等级及以上的侵蚀量占侵蚀总量的81.39%，侵蚀量为轻度侵蚀等级及以下的4.37倍，说明局部区域土壤侵蚀较严重。根据Li等^[21]实测那辣小流域甘蔗坡面沟蚀量表明，2017—2019年甘蔗不同种植模式下年均坡面沟蚀速率为18.9~32.0 t/(hm²·a)，流域甘蔗坡面侵蚀主要是沟蚀，而RUSLE模型主要用于估算坡面细沟或细沟间侵蚀形成的平均土壤流失量，本文用RUSLE模型估算的坡耕地(甘蔗)土壤侵蚀速率18.86 〔t/(hm²·a)〕与实测的坡面沟蚀速率相近，表明利用RUSLE模型估算流域土壤侵蚀总体结果较可靠。子流域SW₃的平均土壤侵蚀速率是SW₁的1.67倍，SW₂与SW₃相差不大，可能原因是子流域SW₃的坡度比SW₁，SW₂大，表明子流域坡度的差异影响土壤侵蚀速率。Abdalla等^[22]在南非甘蔗种植区的试验研究表明，坡度与土壤侵蚀速率呈正相关关系，与本研究结果相符。子流域SW₂，SW₃的林地(桉树)面积大于SW₁，流域林地(桉树)土壤侵蚀速率较大。有研究表明，桉树人工林产生较高坡面径流量和泥沙流失量主要是因为造林抚育措施及大量级降水^[23]。此外，本研究发现那辣小流域5°~25°的坡度范围是土壤侵蚀防治的重点区域。

流域平均土壤侵蚀速率由大到小为：林地(桉树)＞道路＞坡耕地(甘蔗)＞沟渠，林地(桉树)平均土壤侵蚀速率为53.59 t/(hm²·a)，分别是坡耕地(甘蔗)、道路、沟渠的2.84，2.12，27.91倍。王永平等^[24]利用¹³⁷Cs示踪法研究了滇池东大河小流域不同土地利用类型平均土壤侵蚀速率为耕地〔15.39 t/(hm²·a)〕＞草地〔13.09 t/(hm²·a)〕＞林地〔0.82 t/(hm²·a)〕，其林地种植马尾松、云南松和萌芽灌木等乔灌混交林，冠层植被丰富，植被覆盖起到较大抑制水土流失的作用。而本研究的林地土壤侵蚀速率与之相比大许多，主要是因为流域山坡上林地为纯桉树人工林，坡度大，冠层植被稀疏，地表覆盖度较低，种植、砍伐等人为活动易引发严重的土壤侵蚀。王会利等^[25]采用径流场观测发现，炼山后种植桉树第一年土壤侵蚀量达350 t/hm²。有研究表明，南方红壤区林下土壤侵蚀主要原因是人为干扰频繁、林分结构单一、林下植被匮乏^[26]。流域道路的土壤侵蚀速率较大，但其侵蚀量和侵蚀面积占流域侵蚀总量和总面积的5%左右，流域道路几乎没有水土保持措施，Li等^[27]研究表明，2017年那辣小流域农田从小块地变成大块地之后，道路密度、交通强度增加，受大型卡车等运输车辆产生的车辙形成细沟，导致道路侵蚀量显著增加。流域坡耕地(甘蔗)的侵蚀量最大，其值为2 013.05 t/a，分别是林地(桉树)、道路、沟渠的2.65，11.46，1 059.50倍。可能原因是坡耕地面积占流域土地利用总面积的82.85%，坡耕地种植甘蔗，坡度较大，降雨集中且多暴雨，种植方式大部分是顺坡种植，加上翻耕除草，新植甘蔗地表较为裸露，植被覆盖不足以抵挡降雨产生径流对表层土壤的冲刷，导致土壤侵蚀比较严重。有研究表明，在种植甘蔗模式下，每年的泥沙流失量为27~30 t/hm²左右^[11]，甘蔗坡面径流量和侵蚀量为：裸坡＞顺坡蔗地＞横坡蔗地，顺坡蔗地坡面径流量和侵蚀量是横坡蔗地的1.9，2.3倍^[12]，甘蔗地表覆盖是影响地表径流的最重要因素，甘蔗地表凋落物不会改变地表粗糙度，在最低凋落物量(4 t/hm²)下水土流失量减少，入渗率增加^[28]。甘蔗坡面土壤侵蚀和养分流失主要在萌芽期和苗期^[29]，要增加地表覆盖等措施，有针对性地对甘蔗生长前期，特别是大的降雨期坡面土壤侵蚀进行防治。由于各子流域的海拔、坡度、土地利用类型面积和种植管理方式等的差异，引发各子流域的土壤侵蚀速率也不同，可以根据流域的情况采取针对性的措施来防治土壤侵蚀。例如，流域甘蔗可以采取横坡种植，新植宿根按合适的比例种植，合理规划土地利用类型等，尽量减轻流域的土壤侵蚀。

RUSLE经验模型结构相对简单，各因子易于理解，所需数据较少，可以对所研究地区的数据变化进行模型各因子数据输入计算。但在利用RUSLE模型因子计算过程中，由于模型参数的数据精度不一，导致在ArcGIS的栅格计算器计算过程中难免产生误差，而且修正的通用土壤流失方程是一个经验模型，其应用存在较大的局限性，在未来的研究中考虑改进模型的一些参数或使用其他更优化的模型，使得评估结果更具有科学性和可靠性。

4.2 结论

(1) 2020年那辣小流域平均土壤侵蚀速率为22.97 t/(hm²·a)，表明流域土壤侵蚀较严重，应重视流域的水土流失治理工作。

(2) 流域内土壤侵蚀主要发生在5°~25°的坡度范围内，是预防和治理水土流失的重点区域，应采取合理有效的防治措施。

(3) 研究区不同土地利用类型中，平均侵蚀速率最大的是林地(桉树)，为53.59 t/(hm²·a)，分别是坡耕地(甘蔗)、道路、沟渠的2.84，2.12，27.91倍，应合理利用并规划管理好流域的坡耕地(甘蔗)、林地(桉树)，道路等土地利用类型，采取具有针对性的、有效的水土保持措施，以改善流域的土壤侵蚀状况。