土壤侵蚀是土壤与粮食安全面临的主要威胁之一。为阻控土壤侵蚀，需进行土壤侵蚀影响因子的定量化研究。土壤可蚀性 (soil erodibility) 能反映土壤对侵蚀外营力剥蚀和搬运的敏感程度，是组成土壤侵蚀模型的重要因子之一，常使用指标K表示。然而，由于实地测量K值样点密度稀疏，亟需发展K值空间尺度扩展方法。如何利用有限的样点数据提升K值空间分布图的精度，是目前区域尺度K值研究的瓶颈。
南京土壤所梁音研究员课题组围绕土壤可蚀性进行了长期研究。之前基于第二次全国土壤普查数据利用1：50万土壤类型图完成的全国水蚀区土壤可蚀性K值计算及宏观分布，然而通过土壤分类单元对K因子赋值存在精度不高的问题。对于区域尺度K值空间由点及面的扩展方法仍有待提升。
此次研究基于成土因素空间分布与土壤类型具有一致性的思路，利用随机森林回归模型与63项成土因素指标进行K值拟合，结果发现，坡度、高程、最大降雨量、温度季节性、夏冬季节地表温度以及NDVI指标是预测K值的重要环境变量。通过5折交叉验证法对随机森林调参的结果显示，冗余变量对模型精度的影响较小，同时优化了模型的决策树个数 (ntree) 与属性特征个数 (mtry)。最终采用优化的随机森林模型预测得到250米空间分辨率的太湖流域片土壤可蚀性K值空间分布图 (图1c)，并在90%置信区间下得到K值的预测上、下限。该方法显著提升了K值图的空间分辨率，比目前公认的全国土壤可蚀性因子图的平均图斑面积细化了100倍。这一工作通过模型学习方法有效提升了我国土壤可蚀性因子空间尺度扩展精度，并形成利用成土因素预测K值的理论依据，制图结果可服务于水土保持相关部门的管理与决策工作等。
该研究成果发表在水土保持领域专业期刊 (Q1) International Soil and Water Conservation Research上，南京土壤研究所博士后田芷源为论文第一作者，研究员梁音为论文通讯作者。该研究得到了江苏省资助招收博士后研究人员、太湖流域片国家级重点防治区土壤侵蚀因子率定更新项目等资助。
文章链接
图1 土壤可蚀性K值在太湖流域片 (a) 样点分布，(b) 克里格插值，以及 (c) 随机森林模型预测图
本文版权归原作者所有，文章内容不代表平台观点或立场。如有关于文章内容、版权或其他问题请与我方联系，我方将在核实情况后对相关内容做删除或保留处理！