《Frontiers in Soil Science》：Mapping soil vulnerability based on the soil erodibility factor (K) in the El Harrar Watershed (Eastern Pre-Rif, Morocco)
编辑推荐：
水蚀是土壤退化的主要原因，威胁着土壤可持续性和农业生产力，尤其是在环境脆弱地区。其强度受自然因素和人为因素的相互作用影响，包括土壤性质、地形、植被覆盖、气候变率和土地利用方式。本研究旨在通过整合土壤理化特征、基于GIS的制图、统计分析和野外观察，利用土壤可蚀性
水蚀是土壤退化的主要原因，威胁着土壤可持续性和农业生产力，尤其是在环境脆弱地区。其强度受自然因素和人为因素的相互作用影响，包括土壤性质、地形、植被覆盖、气候变率和土地利用方式。本研究旨在通过整合土壤理化特征、基于GIS的制图、统计分析和野外观察，利用土壤可蚀性因子（K）评估摩洛哥El Harrar流域（东前里夫）土壤对水蚀脆弱性的空间变异性。研究人员将先前建立的K因子值按照土壤类型对应关系分配给已制图的土壤单元，以表征整个流域土壤可蚀性的空间分布。结果表明，流域45.15%的区域具有低土壤可蚀性（K < 0.25），主要对应于钙镁质土壤；而23.78%的区域具有极高可蚀性（K > 0.45），主要与铁硅铝质土壤相关。探索性相关分析显示，土壤可蚀性与有机质（R = ?0.70）和黏粒含量（R = ?0.43）呈负相关关系，表明这些土壤属性含量越高，土壤结构稳定性越强，可蚀性越低。这些发现强调了在评估土壤水蚀脆弱性空间变异性时考虑土壤理化特征的重要性，并为识别摩洛哥东前里夫地区土壤保持、可持续土地管理和国土规划的优先区域提供了有用基础。
**摩洛哥东前里夫El Harrar流域基于土壤可蚀性因子（K）的土壤脆弱性制图研究解读**
**一、研究背景与问题**
土壤是农业、生态系统功能和粮食安全不可或缺的自然资源，然而水蚀导致的土壤退化已成为地中海国家（包括摩洛哥）面临的重大环境问题之一。据统计，摩洛哥约40%的土地受到水蚀威胁；西班牙已有约20%的土地退化，另有1%严重退化。摩洛哥尤其受水蚀影响，特别是在里夫山脉地区，约60%进入水库的泥沙来自仅占国土面积6%的区域。根据国家流域管理计划，里夫地区年均土壤流失量超过20 t ha?1 yr?1，前里夫地区为10–20 t ha?1 yr?1。水蚀的演变受气候、地形、岩性、土壤、植被和人为因素的综合影响，其中坡度、降雨侵蚀力、土壤质量和植被覆盖是关键自然控制因子，而不合理的农业实践、过度放牧和土地利用变化可能显著加速土壤退化。
土壤可蚀性表示土壤因降雨和地表径流而分离和搬运的敏感性，是理解水蚀空间变化最重要方面之一。可蚀性因子（K）由土壤质地、有机质含量、土壤结构和渗透性决定。以往研究已表明，这些理化性质间的相互作用主要贡献于土壤侵蚀，但其效应随气候条件、母质、地形、土壤类型和土地利用而异。将K因子与地理信息系统（GIS）集成已成为描述土壤可蚀性空间变异和识别土壤保持优先区域的热门方法。然而，摩洛哥以往研究多集中于基于RUSLE、PAP/CAR和EPM等经验与空间模型估算土壤流失量，对里夫和前里夫地区侵蚀过程提供了良好理解，但普遍侧重于土壤流失速率和整体侵蚀严重程度，而非由土壤理化性质控制的土壤可蚀性内在空间变异性。此外，摩洛哥北部很少有流域研究同时评估土壤单元、K的空间分布、土壤性质与土壤的关系，以及这些差异如何影响土壤和地面侵蚀证据。为此，本研究旨在通过整合土壤理化数据、基于GIS的空间分析、统计关系和野外观察，解决El Harrar流域土壤可蚀性空间变异性问题。研究目标包括：（i）将已知土壤可蚀性因子（K）值赋予El Harrar流域各土壤单元并制图；（ii）将K与主要土壤理化组分关联；（iii）讨论这些发现对摩洛哥东前里夫地区土壤保持、可持续土地管理和国土规划的影响。
**二、研究内容与主要技术方法**
研究人员在摩洛哥东前里夫丘陵地带的El Harrar流域（面积89 km2）开展研究，该流域地质构造以黏土和砂岩为主，土壤类型以钙镁质土壤（45.15%）和变性土（23.78%）为主，自然植被覆盖极为有限，农业用地高强度利用且集中于坡地。研究采用的技术方法主要包括：（1）土壤制图与GIS空间分析：从1:50,000 Taza土壤图数字化土壤单元，在ArcGIS Pro环境中利用GIS技术生成土壤可蚀性空间分布图；K值数据来源于Tribak（2000）等已发表的Oued Lahdar流域土壤分析数据库，按土壤类型对应关系赋予各土壤单元，并采用Manrique标准划分可蚀性等级。（2）统计相关分析：利用Python和MATLAB计算K值与黏粒含量、粉粒加极细砂含量、有机质等土壤变量间的Pearson相关系数，基于7个制图土壤单元的代表值进行土壤单元层次的相关性探索。（3）野外观察验证：通过实地勘察和Google Earth Pro高分辨率影像验证土壤单元对应关系，并记录沟蚀、细沟侵蚀、河岸侵蚀和坡地失稳等侵蚀现象。
**三、研究结果**
**3.1 El Harrar流域土壤特征**
研究识别出七个主要土壤类别：钙镁质土壤占流域面积45.15%，铁硅铝质土壤占23.78%，石质土占11.46%，变性土占9.20%，疏松土占4.18%，弱发育崩积土占3.18%，弱发育冲积土占3.06%。这种异质空间分布反映了流域的地质和地貌特征，表明土壤理化性质多样。野外剖面观察显示，钙镁质土壤具有较深厚的黏土至粉砂黏土层次，有机质含量相对较高；铁硅铝质土壤呈红色且无碳酸盐剖面，有机质含量较低；变性土显示深厚的富黏粒层次；石质土和疏松土通常浅薄且弱发育；弱发育冲积土和崩积土受近期沉积过程影响呈现不同质地。
**3.2 土壤可蚀性因子（K）的空间变异性**
El Harrar流域K值范围为0.23至0.46，变化高度显著。钙镁质土壤占流域面积45.15%，具有相对较高的有机质（2.01%）和黏粒（38.5%），可蚀性最低（K = 0.23）。中等可蚀性土壤（崩积土和石质土）约占14.64%。高脆弱性土壤（K = 0.35–0.45）约占16.44%，包括变性土、冲积土和疏松土，主要分布于溪流和河岸沿线及流域南部。极高脆弱性土壤（K > 0.45）占23.78%，主要由Jbel En nehir山坡上的铁硅铝质土壤组成，其有机质含量为0.95%，黏粒含量为28.4%。
**3.3 统计相关分析**
基于7个土壤单元的代表值，相关分析显示：有机质（OM）与K值负相关最强（R = ?0.70），表明有机质在土壤结构稳定性和抗颗粒分离方面起重要作用；黏粒含量与K值负相关程度略低（R = ?0.43），但黏粒对土壤结构稳定性和抗蚀性有积极贡献；粉粒加极细砂含量的负相关较弱（R = ?0.19）。相关热图进一步展示了K与所研究土壤性质间的负相关关系，其中有机质和黏粒的关联较强，粉粒加极细砂、土壤结构和渗透性的关联较弱。
**四、讨论与结论总结**
讨论部分指出，El Harrar流域K值的差异反映了土壤理化性质与地貌条件的混合作用。钙镁质土壤K值较低因其有机质（2.01%）和黏粒（38.5%）含量高，有助于维持土壤强度并降低颗粒分离倾向；铁硅铝质土壤K值最高，主要因其有机质（0.95%）和黏粒（28.4%）含量低，更易受水蚀。野外观察证实了沟蚀、细沟侵蚀、河岸侵蚀和坡地退化存在于高和极高侵蚀敏感区域，与K因子空间分布一致。本研究K值范围（0.23–0.46）与摩洛哥北部其他流域报道值相符，但比Amzaz流域（0.10–0.82）和Asra流域（0.11–0.82）的报道范围更窄，差异可能归因于母质、土壤组成、地形和土地覆被的不同。统计结果表明有机质和黏粒是控制土壤可蚀性的重要因素，但需谨慎解释，因分析基于仅7个制图土壤单元的代表值，而非更具代表性的样本级数据集；观察到的系数表明土壤单元水平趋势，不应视为因果关系，未来需补充野外采样和实验室工作以提供更精确的统计证据。
研究结论表明：El Harrar流域土壤可蚀性高度受土壤理化性质空间变异控制。制图的K值范围为0.23至0.46，主要土壤类型间差异显著。钙镁质土壤因有机质（2.01%）和黏粒（38.5%）含量相对较高而可蚀性最低，铁硅铝质土壤因有机质（0.95%）和黏粒（28.4%）含量较低而可蚀性最高。探索性统计分析显示，有机质（R = ?0.70）和黏粒（R = ?0.43）与K值负相关，可能是影响土壤稳定性和抗蚀性的重要参数。本研究除提供土壤可蚀性空间分布外，还为将土壤理化信息纳入东前里夫地区侵蚀风险评估和土地管理策略提供了科学依据。可蚀性制图是识别需要土壤保持和可持续土地管理区域的有用工具。由于评估基于已发表的土壤单元数据，未来分析应通过样本级测量和野外验证改进空间数据。基于这些发现，保护措施应强调有机改良以增加土壤有机质、植被覆盖保护和恢复，以及实施农业实践（如在耕作坡地建设和维护梯田）。土壤可蚀性图还应纳入国土规划和环境决策，以确保高度敏感的 Mediterranean 地区土壤可持续性和自然资源维持。
