杨韶洋,刘霞,姚孝友,吴迪,张荣华,齐斐,张洪达

(1.山东省土壤侵蚀与生态修复重点实验室,山东农业大学林学院水土保持系,271018,山东泰安;2.美国密西西比州立大学森林资源系,39759,美国密西西比州斯塔克维尔;3.南京林业大学林学院,210000,南京;4.水利部淮河水利委员会水土保持监测中心站,233001,安徽蚌埠)

降雨侵蚀力(rainfall erosivity)是评价降雨引起土壤侵蚀潜在能力的一个动力指标[1],研究区域降雨侵蚀力空间插值方法,对提高预测精度,了解区域降雨侵蚀力空间分布规律具有重要意义。空间插值方法主要有确定性插值方法(反距离权重插值、样条曲线插值)和地统计插值方法(克里格插值),其中,地统计插值方法既充分考虑了样本点的方向、位置和距离,又能够对数据中存在的趋势和异向性进行处理,选出最优、最适合的模型进行拟合[2]。因此,近年来国内外诸多学者采用地统计插值方法中的普通克里格法(Ordinary Kriging)对不同区域降雨侵蚀力的空间分布特征进行了研究[3-8],并取得了一定的成果;但在土壤养分、温度空间预测的相关研究探索中发现,由多元线性回归和普通克里格法结合而成的回归克里格法(Regression Kriging)比普通克里格法更能提高空间预测精度[9-15],而目前应用回归克里格法对降雨侵蚀力空间分布特征的研究鲜见报道。

笔者借助沂蒙山区88个雨量站点1980—2010年日降水数据和降雨侵蚀力日雨量公式,计算站点年均降雨侵蚀力,采用回归克里格法和普通克里格法对沂蒙山……