张 婧，毕如田，丁皓希，文伟杰，荆颖蔷

（山西农业大学资源环境学院，山西太谷 030801）

土壤有机质作为表征耕地肥力大小的重要指标，其含量在农作物的整个生长过程中发挥着非常重要的作用[1]。土壤有机质具有较强的时空变异性，其空间规律一直是许多学者研究的热点问题。大量研究表明，土壤有机质的空间变异会受各种因素的影响，如地形、土壤类型、耕作模式、土地利用方式等[2-3]。但是关于多种因素协同对有机质的空间影响目前还少有研究。

Kriging 空间插值法是一种土壤属性空间预测最普遍的方法[4-5]，但是普通Kriging 插值在平坦、地形单一的区域预测结果较好，而在复杂地区插值效果一般。随着科学技术的进步，Kriging 插值引入其他参数作为辅助变量的方法已经得到广泛的应用，而且预测精度较高。例如，徐占军等[6]基于分区Kriging 插值法对煤炭开采沉陷区土壤有机碳含量进行空间预测，结果显示，分区Kriging 预测结果的精度比普通Kriging 预测的更高；MISHRA 等[7]以印第安纳州为研究区，利用普通Kriging 模型和剖面深度分布函数分别预测了表层土壤有机碳的含量，结果表明，剖面深度分布函数预测精度高于普通Kriging 模型。

地形是引起土壤有机质空间变异的主要因子，而地形起伏度是表征地形特征的重要参考指标之一[8]。地形起伏度是表示某一确定区域中最高海拔和最低海拔之差[9]，它是描述区域地形特征的一个宏观性指标，用以表征地面的起伏状况和切割程度，被广泛应用于自然地质灾害等定量评价方面[10-14]，目前，结合地形起伏度进行有机质空间规律的研究还相对较少。……