土壤可蚀性因子(K)是衡量土壤在降雨作用下被侵蚀难易程度的重要参数,其值越大,土壤越易受侵蚀。在处理好HWSD数据后(HWSD数据简单处理方法),借鉴如下图公示,计算K因子。

在HWSD数据集中,K命名可能与公示上不太相同,如下所示:

T_SILT: Real (表层粉沙粒含量); (单位:% wt.)
T_CLAY: Real (顶层粘土含量); (单位:% wt.)
T_OC: Real (顶层有机碳含量); (单位:% weight)

因此,在ArcGIS 栅格计算器中,公示如下所示。其中,SN1需要提前使用栅格计算器计算好。

0.1317 * (0.2+0.3*Exp(-0.0256*"Tsand.tif"*(1-"TSILT.tif"/100)) )* Power(("TSILT.tif"/("Tclay.tif"+"TSILT.tif")),0.3)*(1-((0.25*"Toc.tif")/("Toc.tif"+ Exp(3.72-2.95*"Toc.tif") )))*(1-((0.7*"SN1.tif")/("SN1.tif"+ Exp(-5.51+22.9*"SN1.tif") )))

此外,除了EPIC方法计算K因子,还可以使用其他数据和其他方法计算K因子,如下所示是USLE中K因子计算方法。该计算方法的结果可以直接在青藏高原数据中心等下载别人分享的成果。

我对比了东北某地区 EPIC计算结果和USLE计算结果的均值,EPIC-K均值为0.0358,青藏高原数据中心下载的数据K因子均值为0.0354(去空值前),均值差别不大。

Logo

DAMO开发者矩阵,由阿里巴巴达摩院和中国互联网协会联合发起,致力于探讨最前沿的技术趋势与应用成果,搭建高质量的交流与分享平台,推动技术创新与产业应用链接,围绕“人工智能与新型计算”构建开放共享的开发者生态。

更多推荐