基于谷歌地球引擎的西欧巴尔干地区集水区优先级划分之改良侵蚀潜力法土壤侵蚀敏感性云端制图

《Land》：Cloud-Based Mapping of Soil Erosion Susceptibility for Catchment Prioritisation in the Western Balkans Using a Modified Erosion Potential Method

　　

土壤侵蚀被认定为西欧巴尔干（Western Balkans, WB）最普遍的土地退化形式之一；该区域面积208,052 km2，可蚀性基质、气候差异及长期土地利用压力导致侵蚀格局高度变异。研究人员通过谷歌地球引擎（Google Earth Engine, GEE）实现改良侵蚀潜力法（Erosion Potential Method, EPM），在30 m分辨率下计算，并依据欧盟水文（EU-Hydro）数据库聚合为9524个集水区，以无量纲侵蚀系数（Z）作为核心敏感性指标。模拟显示侵蚀强度呈强南北不对称：国家尺度Z最高为阿尔巴尼亚0.47、北马其顿0.38、黑山0.37，塞尔维亚0.29与波斯尼亚和黑塞哥维那0.28较低；约19.8%的WB面积属高至极高敏感集水区，集中于亚得里亚海向山前及爱琴河流域源头。验证含1379个野外与正射影像核验点，曲线下面积（Area Under Curve, AUC）=0.828。该框架虽不能替代精细野外建模，但为区域筛查、集水区优先级划分及水土保持干预提供透明、可复现、可扩展基础。

研究背景方面，土壤侵蚀是全球广泛且严重的土地退化过程，影响农业生产力、生态系统功能与集水区水文响应。欧洲农业用地平均土壤流失约3.13 t ha^?1 yr^?1，而巴尔干半岛被视作欧洲最具侵蚀性的部分，前南斯拉夫地区半个多世纪以来主要由Gavrilovi?提出的侵蚀潜力法（Erosion Potential Method, EPM）主导制图。该方法可输出侵蚀系数（Z）与年沉积产量，但传统执行依赖专家野外判定，且西欧巴尔干六国地质图、土壤图、土地利用数据在比例尺与分类体系上不兼容，后2000年农业废弃与覆被变化使历史地图失效。云处理平台谷歌地球引擎（Google Earth Engine, GEE）提供一致全球数据集与弹性算力，使跨国境、次大陆尺度统一建模成为可能。研究人员由此开展首项全WB区域尺度、GEE中改良EPM的土壤侵蚀敏感性评估，以服务集水区优先级划分与土壤保护政策；《Land》刊发本研究。

关键技术方法上，研究人员以GEE为平台，输入为Copernicus GLO-30数字高程模型、CHIRPS降水、ERA5_Land气温、SoilGrids与OpenLandMap土壤属性、CORINE 2018土地覆被、Sentinel-2 Level-2A影像；以EU-Hydro的9524个集水区为聚合单元；验证点集为六国田野与正射影像判读得到的1379个点（602侵蚀点、777弱蚀或稳定点），另用Ukrina、Vrbas、Kolubara、Upper Bregalnica常规GIS-EPM图及Kalimanci、Grad?e、Drenova、Potpe?水库测深记录作交叉比对。

研究结果如下。

1. 2.
   Materials and Methods：研究人员以改良EPM在GEE中实现，Z = Y · Xa · (φ + √I)；基岩抗蚀用SoilGrids的K因子与OpenLandMap土壤pH代理Y = (K + (1 ? sPH))/2，土地保护Xa取自CORINE 2018，可见侵蚀φ由Sentinel-2裸土指数（Bare Soil Index, BSI）经S形函数得到，坡度I由GLO-30导出并除15截断[0,2]；气候势由ERA5_Land与CHIRPS算得，年总侵蚀量Wy = T · H · π · √Z³ · f，输沙比Ru由EU-Hydro形态参数算得，Gy = Wy · Ru；验证分模型间栅格比对、水库沉积量级比对、全区点统计三类。
2. 3.
   Results：WB平均Z=0.33；阿尔巴尼亚、北马其顿、黑山最高，科索沃、塞尔维亚、波黑较低。五类分级中极低位40.8%、低位26.7%、中位23.9%、高位6.9%、极高位1.7%；高敏感集水区集中在亚得里亚向山前与Vardar源头。栅格先行与集水区均值Z高度一致（r=0.981，92.2%集水区类别不变）。比常规GIS-EPM：Ukrina 0.27对0.30、Kolubara 0.34对0.35、Vrbas 0.24对0.22、Upper Bregalnica 0.36对0.40，类面积相似度91.9%–93.9%。水库比对：GEE的Gy均高于实测沉积，Potpe?高约7%、Grad?e超一倍，主因为Ru仅含几何输沙、缺连通性项。点验证AUC=0.828（空间分块交叉验证同值），Z=0.4传统界点准确率79%、Cohen’s Kappa=0.557；最优Youden阈值Z*=0.347；频率比（Frequency Ratio, FR）由极低位0.46单调升至极高位80.1，卡方极显著。

Discussion：研究人员指出栅格先行引入轻微聚合偏倚但不改排序；模型适宜区域筛查而非工程精算；Y代理难分坚硬基岩与可蚀土壤，BSI（Bare Soil Index）把裸露碳酸盐岩误读为产沙是碳酸盐区偏高主因；气候与DEM粗化平滑地形雨影；建议未来分离岩性层、加碳酸盐校正、用LiDAR验坡、以实测产沙率标定Ru、做概率化田野验证。

1. 5.
   Conclusions：研究人员给出首项全WB云端改良EPM制图，区域均Z=0.33，约五分之一面积入高及极高类，9524集水区排名为流域级干预可直接用；三层验证支撑其作为筛查工具，AUC=0.828且空间分块不变，传统Z=0.4界值统计上仍合理。结论明确：框架用于区域筛查与跨境协调；局限在Y代理与碳酸盐裸岩误判；后续七项发展为建统一岩性层、碳酸盐因子、φ现场标定、高分辨气候、LiDAR坡检、连通性显式化、概率田野验证；产出服务欧盟绿色议程与土地退化报告，且GEE架构可随新数据再生评估，长期走向近实时降雨与气候情景的动态评估。