基于WaTEM/SEDEM模型的沂河流域土壤侵蚀产沙模拟

基于WaTEM/SEDEM模型,结合临沂水文站和角沂水文站的输沙数据对模型进行校正和验证,分析模拟1975—2015年沂河流域侵蚀产沙的时空变化特征,并进一步研究降水、地形位和土地利用变化对流域侵蚀产沙的影响。结果表明：① 沂河流域输沙能力系数K_tc-low和K_tc-high在40 m和150 m组合下效果最优,模型在沂河流域具有较好的适用性。② 1975—2015年,沂河流域主要以侵蚀为主,微度侵蚀所占面积最大,其次是剧烈侵蚀,沉积主要分布在河谷处;流域侵蚀强度呈现先增加后减少的趋势,侵蚀模数由1975年的30.92 t·hm^-2·a^-1增加至1995年的49.32 t·hm^-2·a^-1再下降至2015年的29.60 t·hm^-2·a^-1;各县（区）平均侵蚀模数为沂水县>费县>沂南县>沂源县>蒙阴县>平邑县>兰山区。③ 沂河流域土壤侵蚀产沙强度的变化是降水、地形、土地利用等综合作用的结果。1975—2015年,流域降雨侵蚀力呈现先降低后升高又降低的变化趋势,各县（区）平均降雨侵蚀力为费县>兰山区>沂南县>蒙阴县>平邑县>沂水县>沂源县,降雨侵蚀力时空变化与流域侵蚀产沙强度时空变化并不完全一致;地形位等级空间分布与流域侵蚀产沙强度空间分布基本一致,侵蚀产沙的优势地形位区间是4~6级,即高程75~428 m,坡度5°~39°;耕地和林地的转化是土壤侵蚀强度转化最主要的原因,林地转化为耕地使侵蚀强度面积升高3389.97 hm²·a^-1,耕地转林地则使侵蚀强度面积降低2216.65 hm²·a^-1,草地与其他土地利用类型的转化对流域侵蚀强度影响较小。该研究可为区域土地利用方式调整和水土流失调控提供参考。