王雪珊, 沈庆松,, 高凤杰, 张兴义, 张少良, 王 力

(1.东北农业大学 资源与环境学院, 哈尔滨 150030; 2.中国科学院 东北地理与农业生态研究所,哈尔滨 150081; 3.西北农林科技大学 水土保持研究所, 陕西 杨凌 712100)

土壤速效磷(Available phosphorus,AP)是衡量土壤供磷能力和磷素流失风险的重要指标[1-2]，该指标能够反映农田土壤肥力特征，研究其空间分布规律对于区域施肥管理和农业面源污染防控等具有重要意义[3]。受生物地球化学循环过程、自然环境和人类活动等因子影响[4-5]，农田土壤AP在空间分布上存在高度空间异质性。研究证实AP预测的精度直接影响农田土壤磷素的管理水平[6]，因此利用有限信息，改进模拟方法提高AP空间分布预测模型的精度十分迫切。

空间插值方法是利用采样点数据通过数学模型对研究区内未知点进行模拟预测，并能够获得连续的属性值空间分布特征[7]。回归克里格法[8]将多元线性回归模型与克里格方法相结合，同时引入多个环境因子参与建模，并考虑土壤属性空间分布过程中的结构性因子与随机性因子，其空间插值结果更加精确[9-10]。Li[10]和Zhang等[11]提取地形和土地利用等变量利用回归克里格法预测土壤有机质的空间分布特征，研究结果表明回归克里格方法要优于单变量空间插值方法。但是空间非平稳性的存在导致回归克里格方法不能有效地捕捉土壤属性空间变异的局部特征[12-13]。地理加权回归模型利用局部加权最小二乘法进行逐点参数估计和空间预测，能够在引入环境变量的同时基于空间非平稳性假设揭示更多的局部信息，该模型已逐渐成为土壤属性空间预测的有效手段[14-15]。……