文献参数: 王娟, 史燕东, 赵辉, 等.西南高山峡谷区坡耕地水土流失图斑识别及宜治理性评价[J].水土保持通报,2026,46(1):172-180. Citation:Wang Juan, Shi Yandong, Zhao Hui, et al. Soil erosion mapping patche indentification and conservation suitability evaluation of sloping farmland in high mountain and canyon region of southwest China [J]. Bulletin of Soil and Water Conservation,2026,46(1):172-180.
土壤侵蚀具有类型复杂,分布广泛,侵蚀强度大等特点
[1-2]。它作为世界上最严重的生态环境灾害之一
[3],已成为全球可持续发展进程中面临的突出挑战。土壤侵蚀不仅直接导致土地退化,耕地资源损失,还会引发泥沙淤积,水质恶化等连锁生态问题,治理难度大,生态恢复周期长。在农业方面,其危害表现为作物产量的短期下降,以及土壤肥力、结构和生产力的长期衰退
[4-5],严重制约农业可持续发展和区域经济稳定。近年来,中国依托水土流失动态监测网络,逐步实现对水土流失状况的系统化、常态化掌握。监测数据显示,6°以上坡耕地是中国高强度水力侵蚀分布最为集中的地类,强烈及以上水蚀面积比例超过60%,这类高风险坡耕地主要集中于地形起伏大、降雨条件特殊的云、贵、川、鄂、渝5省(市),水土流失敏感性高,治理需求迫切,成为当前中国水土流失实现“减量,降级”目标的关键区域。
西南高山峡谷区地处中国大横断地带,是中国生态安全屏障体系的重要组成部分,同时也是生态脆弱性突出的典型区域
[6]。该区以高山峡谷地貌为主,受特殊地形构造、地质条件与气候因素共同影响,面临土壤侵蚀和石漠化等土地退化问题的严峻挑战
[6]。其中,川西南高山峡谷保土减灾区作为坡耕地土壤侵蚀治理的重点区域,由于特殊的自然地理条件和水土流失客观实际,存在大量不宜治理或治理潜力较低的区域。目前,该区域的土壤侵蚀研究相对较少,在区域尺度主要采用中国土壤侵蚀模型(Chinese soil loss equation, CSLE)开展
[7-8],导致坡耕地水土流失空间位置定位难度大,且存在侵蚀强度分布边界不清问题,难以满足坡耕地水土流失治理任务的分解和落实,对坡耕地水土保持精准治理的工作带来了一定困难。在微观尺度,研究多集中于区内典型小流域或坡面尺度的侵蚀特征与影响因素分析
[9-11],虽有助于机理认识,却难以支撑区域的坡耕地水土流失监测与治理规划工作。因此,如何准确识别并空间定位坡耕地水土流失图斑,实现精准治理,对保障国家生态安全、粮食安全及推进生态文明建设具有重要意义
[12-13]。
在水土保持领域,部分学者提出“土壤侵蚀图斑”“水保斑”“划分设计图斑”等
[14-15]基础空间单元的新概念,但目前仍处于初步试验阶段,存在与水土保持综合治理工作契合度不高,针对性不强,划分方法不完善,实际应用效果欠佳等问题。自2014年以来,水利部推动水土保持“图斑化”精细管理,并于2021年进一步提出了“水土流失落地图斑”的新概念,但目前还未有相对完善的“水土流失图斑”识别方法,在实际应用中仍然面临着内部要素复杂,空间位置准确率低,高强度侵蚀区域未能准确识别或难以在区域尺度应用等问题。为此,本研究以川西南高山峡谷保土减灾区美姑县为例,基于水土流失遥感监测结果,构建坡耕地水土流失图斑识别方法,并利用坡耕地水土流失图斑开展宜治理坡耕地水土流失图斑评价,以期为区域坡耕地水土流失治理的任务分解、规划制定与工程布局提供科学数据支撑和决策依据。
1 研究区概况
美姑县(102°53′—103°21′E, 28°02′—28°54′N)地处横断山脉东北缘与四川盆地西南边缘的交汇地带,位于四川省凉山彝族自治州东北部,属于典型的高山峡谷向山原过渡区域(
图1)。在地貌类型上,该县以中山山原窄谷地貌为主,整体地势呈现北高南低、东高西低的特征,山体起伏显著,河谷深切,地形破碎度高
[16]。县域内水系极为发育,溪河密布,共有大小河流159条,主要分属金沙江水系和岷江水系,河流密度大,水文网络错综复杂
[17]。据2023年水土流失监测统计,美姑县全域水土流失类型均为水力侵蚀,坡耕地水土流失面积占坡耕地面积的88.96%,水土流失比例高,在西南高山峡谷—川西南高山峡谷保土减灾区这一生态功能区划中,美姑县坡耕地水土流失强度与比例均居首位,具有代表性。
2 数据与方法
2.1 研究数据
(1) 土地利用。土地利用解译和水土保持措施解译采用研究区2023年3—5月2 m分辨率的高分一号卫星(GF-1)遥感影像(https:∥data.cresda.cn/#/home)获取。
(2) 土壤侵蚀因子。降雨侵蚀力
[18]、土壤可蚀性
[19]、坡度因子和坡长因子
[20]、植物覆盖与生物措施
[21-22]、耕作措施数据均来源于2023年水土流失动态监测成果数据库,空间数据分辨率均为10 m。
2.2 研究方法
2.2.1 水土流失遥感监测
采用中国土壤流失方程 (CSLE)计算土壤侵蚀模数
[23],按照《土壤侵蚀分类分级标准(SL 190—2007)》获取5个强度等级的水土流失面积。
2.2.2 坡耕地水土流失图斑识别方法
利用缓冲区空间聚类的分析方法,针对耕地的水土流失空间分布特征,通过设置缓冲参数,将空间离散分布的土壤侵蚀栅格聚类形成集中连片的坡耕地水土流失图斑。具体方法为: ①将土壤侵蚀强度栅格数据转化为二值化栅格图层,将微度等级赋值为0,轻度及以上等级均赋值为1,得到水土流失栅格图层A。 ②将土地利用矢量数据转为土地利用栅格数据,保留耕地、园地、林地和草地这4类1级土地利用类型;土地利用类型具体可细分为水田、水浇地、旱地、果园、茶园、其他园地、有林地、灌木林地、其他林地、天然牧草地、人工牧草地、其他草地;将水土流失栅格二值化图层A和土地利用栅格数据通过空间相乘运算,获取每个水土流失栅格的土地利用类型信息,得到附有土地利用信息的水土流失栅格图层B。 ③将水土流失栅格图层B转为矢量图斑,提取耕地水土流失图斑并导出为独立图层;随后开展缓冲区分析,以每个图斑中心为圆心,分别设置10,20和30 m的正向缓冲半径(依据本研究采用的10 m分辨率空间数据,以10 m为步长设定),统计研究区不同缓冲范围内耕地图斑内的水土流失面积所占比例。结果显示,随着缓冲半径增大,水土流失面积所占比例逐步下降,分别为99.58%,92.34%和88.57%。 ④根据步骤3结果,当缓冲半径为10 m时,所含非水土流失面积极少,与原始耕地水土流失栅格在空间上高度一致,表明该尺度下边界吻合性最佳。因此,选定10 m作为缓冲半径,并进一步设置-10 m的反向缓冲距离,以剔除冗余边界、还原真实地块形态,最终获取更精确的耕地水土流失图斑结果。 ⑤为实现图斑边界与实际地理边界相吻合,采用高精度土地利用解译矢量数据(基于2 m分辨率遥感影像人工目视解译生成,已综合考虑地物边界、地形要素及道路等人工设施)对初步聚合结果进行切割与边界校正。具体校正步骤包括:将水土流失图斑与全地类土地利用矢量数据进行相交处理,即利用不同土地利用类型的解译矢量边界对连片聚合图斑进行分割;保留相交后地类属性为耕地的水土流失图斑,去除非耕地图斑,确保其空间位置与土地利用解译数据一致;最后,依据《年度水土流失动态监测指南》要求,剔除面积小于400 km²的图斑,最终完成图斑优化。 ⑥分别计算耕地水土流失图斑内的水土流失面积及所占比例、中度及以上高强度侵蚀面积及所占比例、图斑内平均坡度信息,得到不同坡度等级的水土流失图斑结果。
2.2.3 坡耕地水土流失图斑识别精度验证方法
(1) 面积精度。对比坡耕地水土流失总面积、各侵蚀强度等级水土流失面积与水土流失遥感监测面积的相对误差。相对误差计算公式如下:
式中:Q为面积相对误差(%); Sa为基于水土流失图斑计算的水土流失面积(km²); Sb为基于遥感监测统计的水土流失面积(km²)。
(2) 野外复核验证。依据《年度全国水土流失动态监测技术指南》,选择该县坡耕地水土流失图斑总数的0.5%共105个水土流失图斑,进行实地复核是否为坡耕地水土流失图斑。结果得到,判断结果确为坡耕地水土流失图斑的准确率为91.43%,其空间位置精度高于85%,可满足后续宜治理评价的相关分析研究。
2.2.4 宜治理评价分析方法
综合考虑坡耕地水土流失图斑治理的自然适宜度和区位适宜度2个因素层,涉及6个因子层,具体详见
表1。其中,自然适宜度基于水土流失图斑空间叠加坡度栅格数据和土壤侵蚀模数栅格数据获取。欧式距离法是一种常用距离定义方法,用于衡量多维空间中两点间的绝对距离
[24]。因此,区位适宜度采用欧氏距离法进行评价坡耕地水土流失图斑到道路、水体和城镇用地的距离。
指标权重反映了各指标对适宜性评价的影响程度,本文采用层次分析法(AHP)与专家打分法相结合的方式,确定不同层次指标的权重及量化标准
[12]。具体过程包括:首先,基于评价指标体系构建层次结构模型;其次,邀请专家使用1~9标度法对同一层次内各指标的相对重要性进行两两比较,构建判断矩阵。本轮共邀请10位专家参与咨询,专家专业背景涵盖水土保持、地理信息系统、遥感监测与土地资源管理等领域,均具备高级职称或10 a以上相关研究与实践经验,以确保判断的科学性与权威性。最后,借助Yaahp软件计算各判断矩阵的权重值,并进行一致性检验,所有判断矩阵的一致性比率(CR)均小于0.1,通过检验,表明专家判断逻辑一致,权重结果可信。其中,适宜性影响因子中的地面坡度、土壤侵蚀强度的分类标准参考《土壤侵蚀分类分级标准(SL 190—2007)》,与水源、道路等距离的分类标准参考《水土保持技术规范》及相关研究,并结合美姑县实际情况和城镇用地分布特征,综合上述AHP分析结果,得到坡耕地水土流失图斑治理适宜性评价因子权重与分级标准(
表2)。
3 结果与分析
3.1 坡耕地水土流失分布特征
通过土地利用遥感解译,美姑县耕地面积603.18 km²,其中水平梯田面积为57.62 km²,所占比例为9.55%。空间叠加坡度栅格数据,统计不同坡度等级下的耕地面积(不包括水平梯田面积),结果得到耕地主要集中在15°~25°和6°~15°区域,分别占美姑县不同坡度等级耕地总面积的42.45%和34.29%(
图2)。
通过遥感监测,美姑县2023年耕地水土流失面积为485.24 km²,以轻度侵蚀和中度侵蚀为主,分别占耕地水土流失总面积的33.87%和28.02%;强烈及以上侵蚀占耕地水土流失总面积的38.11%,高强度侵蚀问题相对突出。总体上,耕地水土流失呈现出随坡度等级升高,水土流失占比逐渐增大趋势,在15°以上坡耕地上的水土流失比例均高于95%,水土流失发生率高。从侵蚀强度看,轻度侵蚀集中在6°~15°,占该县坡耕地轻度侵蚀面积的57.14%;中度侵蚀主要集中在6°~15°和15°~25°,分别占该县坡耕地中度侵蚀面积的36.75%和45.77%;强烈侵蚀集中在15°~25°,占该县坡耕地强烈侵蚀面积的68.21%;极强烈和剧烈侵蚀主要集中在15°~25°和≥25°,15°~25°和≥25°极强烈侵蚀占该县坡耕地极强烈侵蚀面积的40.69%和47.78%,15°~25°和≥25°剧烈侵蚀占该县坡耕地剧烈侵蚀面积的50.05%和49.42%。
3.2 坡耕地水土流失图斑识别
为评估坡耕地水土流失图斑识别的精度,对识别结果与遥感监测数据进行系统对比分析。统计结果表明(
表3),坡耕地水土流失图斑总面积与遥感监测结果相差2.04 km²,相对误差仅为0.42%,显示出较高的一致性。在各侵蚀强度等级上,坡耕地水土流失图斑面积与遥感监测结果的相对误差均控制在5%以内。其中轻度、中度、强烈及以上侵蚀等级的误差分别为0.32%,-0.79%和1.38%。总体上,各侵蚀强度等级的水土流失图斑分级准确率均超过95%,说明坡耕地水土流失图斑识别结果具有较好的可靠性与适用性。
从空间分布来看(
图3),美姑县共识别出坡耕地水土流失图斑约2.1万个,水土流失主要集中分布于县域中部以南区域。强烈及以上高强度侵蚀图斑在空间上呈现集聚特征,主要分布于中部地区。统计显示,高强度侵蚀面积占比超过50%的水土流失图斑数量占总图斑数量的57.53%,进一步印证了美姑县坡耕地水土流失中以高强度侵蚀为主导的现状,也凸显该区域在治理过程中应重点关注高强度侵蚀图斑的识别与防控。
3.3 宜治理坡耕地评价分析
基于综合适宜性评价结果,高适宜治理、中适宜治理、低适宜治理和不适宜治理4类坡耕地水土流失图斑的空间分布如
图4a所示。从数量结构来看,高适宜、中适宜、低适宜和不适宜治理的图斑个数分别为2 817,13 695,4 239和253个,其对应的面积比例分别为19.72%,66.36%,13.50%和0.42%。中适宜治理图斑在数量与面积上均占主导地位,表明该区域治理潜力总体较大,但治理条件存在显著梯度差异。在治理适宜性较高的类别中,高适宜与中适宜治理图斑除受道路可达性和水源获取条件影响外,其差异主要体现在坡度与侵蚀强度两方面。其中,高适宜治理坡耕地水土流失图斑平均坡度小于20°,侵蚀强度以轻度为主,高强度侵蚀比例低于20%,立地条件良好,治理难度较小,治理后易于实现水土流失的“减量,降级”目标。中适宜治理坡耕地水土流失图斑平均坡度集中在15°~28°,侵蚀强度以中度及以上为主,其中超过60%的图斑中高强度侵蚀占比超过50%,表明该类区域侵蚀问题较为突出,治理需求迫切。低适宜与不适宜治理图斑受交通条件限制显著,主要集中在高速公路、国道及其他主干道路难以覆盖的偏远区域,基础设施支撑能力弱,治理难度大、经济性较差,近期治理可行性较低。总体上,以上评价结果初步明确了美姑县坡耕地水土流失治理的优先顺序与空间布局,可为分类施策、分级治理提供科学依据,有助于提升水土保持工程实施的精准性与实效性。针对高适宜与中适宜治理区域,本研究进一步依据坡度等级提出差异化治理措施建议(
图4b)。在2°~6°缓坡区,建议推进高标准农田建设,提升土地生产能力与抗蚀性能;6°~15°区域宜修筑水平梯田,有效截短坡长,减缓径流;15°~25°陡坡耕地可采用坡式梯田或植物篱措施,结合结构性与生物性治理;25°以上极陡坡耕地则应以退耕还林还草为主,或实施保土耕作制度,控制面源侵蚀,恢复植被覆盖。
4 讨 论
(1) 坡耕地宜治理评价合理性分析。在西南高山峡谷区,地形破碎、立地条件异质性强,目前全国尚缺乏统一的山区复杂地形下坡耕地治理适宜性评价指标体系。为此,本文以坡耕地水土流失图斑为基本评价单元,遵循“分坡度,先易后难,分步治理”的原则,从自然适宜度与区位适宜度两个维度构建评价框架,共涵盖6个核心因子。在自然适宜度方面,采用平均坡度反映工程治理的难易程度,以高强度侵蚀占比表征治理过程中实现侵蚀“降级”的挑战性;在区位适宜度方面,以农村道路等一般道路条件评价治理可达性,将高速公路、国道等主干道路视为治理的不利因素(因其出口稀疏,直接接入性差),同时结合河流、沟渠、坑塘水面等水域要素评估灌溉水源的获取便利性。所选指标数据源可靠,易于获取,操作性强,能够支撑区域尺度下的治理优先性初步判别。然而,本研究构建的宜治理评价方法及其后续治理措施配置仍存在一定局限性。在评价层面,当前指标体系尚未纳入土壤质地、有机质含量、耕层厚度及降水量等关键自然因子
[25-26],导致评价结果在反映坡耕地真实地力条件与水文响应方面存在不足,可能影响治理适宜性的精确判别。在治理措施设计方面,当前主要依据坡度等级推荐相应工程类型,如缓坡区推行高标准农田建设,陡坡区建议坡式梯田或退耕还林,该分类方式仍较为粗放,未进一步结合土层厚度、石砾含量、水源保障等立地条件进行差异化设计,导致措施配置的针对性和实操性有所不足。受数据获取限制与模型复杂度制约,未对上述因子开展系统分析,因此在自然立地条件方面的评价仍存在盲区,当前结果更适用于宏观层面的治理分区与初步规划。未来研究应在现有框架基础上,逐步引入土壤属性、气候水文特征及微地形信息,构建多因子、多层级的精细化评价体系;在治理措施设计上,应结合区域土地质量本底,开展治理措施精细化、差异化设计。如在土层薄弱,石砾含量高或水源匮乏区域,优化梯田结构,增强植被恢复的生态适应性,从而提升治理措施的科学性与实际应用效益。
(2) 坡耕地土壤侵蚀防治途径。坡耕地水土流失治理是土地整治与农业可持续发展的重要环节,其核心在于通过系统性工程与生物措施改善农业生产条件
[27],减少水土流失面积,提升土地质量与农业抗灾能力
[28]。美姑县内紫色土广泛分布,有机质含量低,表层土浅薄,抗侵蚀能力弱
[29],区域内大部分坡耕地土壤侵蚀速率远大于土壤恢复能力。即大于《土壤侵蚀分类分级标准(SL 190—2007)》提出的西南紫色土区域土壤允许流失量标准。长期高强度侵蚀不仅造成耕层土壤流失,更会加剧土地质量退化进程
[5]。作为西南高山峡谷区耕地资源的重要组成部分,坡耕地的水土流失与养分损失直接影响区域农业产能与农村经济稳定
[30],也削弱了农业生态系统的韧性与减灾能力。因此,推进坡耕地水土流失治理,必须采取以“保护—改良—合理利用”为主线的综合治理路径,结合水土保持措施体系,实现坡面水土资源的可持续管理。
本研究基于坡耕地水土流失图斑的分布规律与属性特征,提出了针对不同坡度等级的治理措施建议,如缓坡区推进高标准农田建设,中陡坡区布设梯田工程,陡坡区实施退耕还林还草等。需要说明的是,坡耕地治理措施布设成果是基于坡耕地水土流失的分布规律和水土流失图斑属性,制定的可布设措施建议,成果仅作为参考案例并不代表实际规划安排。由于在实际推进中,坡耕地治理策略需综合考虑地方政策导向,行业规划衔接,多部门协同机制、资金投入规模等现实因素。同时还需响应群众在配水、通路、产业结构调整等方面的实际需求,综合实施工程措施、林草生物措施、优化种植模式和坡面水系工程措施等措施体系,建立坡耕地水土流失治理模式体系。此外,应积极探索多元化投入机制,引导社会力量与民间资本参与产业开发型治理项目,加强涉农资金整合与项目协同
[31],提升治理工作的持续性和实效性,从而实现坡耕地水土流失治理从单点治理向系统治理、从政府主导向多元共治的转变。
5 结 论
(1) 美姑县坡耕地水土流失问题较为突出,坡耕地水土流失面积占坡耕地面积的88.96%,是西南高山峡谷川西南高山峡谷保土减灾区中坡耕地流失严重区域的典型代表。全县15°以上坡耕地上的水土流失占比均高于95%,强烈及以上高强度侵蚀占耕地水土流失总面积的38.11%,反映出该区域高强度侵蚀问题相对突出,治理紧迫性高的特点。
(2) 坡耕地水土流失图斑识别共计获取2.1万个坡耕地水土流失图斑,流失面积与遥感监测的面积相对误差为0.42%,各侵蚀强度等级的识别面积相对误差均控制在5%以内,表明图斑识别精度较高、结果可靠。进一步分析发现,近60%坡耕地水土流失图斑中高强度侵蚀比例高于50%,说明高强度侵蚀不仅范围广,而且空间集中性明显,治理难度较大。
(3) 从自然适宜度与区位适宜度两个维度综合评价坡耕地水土流失图斑的治理可行性,结果显示:高适宜治理、中适宜治理、低适宜治理和不适宜治理的坡耕地水土流失图斑面积比例分别为19.72%,66.36%,13.50%和0.42%。结果总体掌握了宜治理坡耕地水土流失图斑的空间分布情况,研究结果可为西南高山峡谷区推进坡耕地分级治理、差异化布设水土保持措施提供科学依据和规划支撑。
国家重点研发计划项目“西南高山峡谷区土壤侵蚀预报优化与智慧监测体系”(2022YFF1302902)