基于不同插值方法的三江平原白浆土磷素空间分布预测及其适用性分析  

Applicability of spatial interpolation methods to predict total phosphorus in the typical irrigated areas of the Sanjiang Plain

作　　者：张迪[1] 姜柏志 刘国辉 张慧[3] 聂凡 孙琦 纪明元 ZHANG Di;JIANG Baizhi;LIU Guohui;ZHANG Hui;NIE Fan;SUN Qi;JI Mingyuan(College of Resources and Environment,Northeast Agricultural University,Harbin 150030,China;Heilongjiang Agricultural Environment and Cultivated Land Protection Station,Harbin 150000,China;College of Public Administration and Law,Northeast Agricultural University,Harbin 150030,China)

机构地区：[1]东北农业大学资源与环境学院,哈尔滨150030 [2]黑龙江省农业环境与耕地保护站,哈尔滨150000 [3]东北农业大学公共管理与法学院,哈尔滨150030

出　　处：《中国生态农业学报（中英文）》2021年第8期1405-1416,共12页Chinese Journal of Eco-Agriculture

基　　金：国家重点研发计划项目(2017YFD0300500)资助。

摘　　要：土壤磷素含量是反映农业生态系统土壤肥力的重要指标。准确预测磷素空间异质性是评价土壤生产力和质量的关键。本研究采用反距离加权法(IDW)、径向基函数法(RBF)、普通克里金法(OK)、全局多项式法(GPI)、局部多项式法(LPI)、地理加权回归(GWR)和地理加权回归克里金法(GWRK)等插值方法,分别预测了三江平原白浆土典型灌区八五三、七里沁以及大兴灌区土壤磷素分布,并运用交叉验证法,通过平均误差(ME)、均方根误差(RMSE)和改进效果(RI)对各种方法精度进行比较,以期确定同一土壤类型不同采样密度土壤中磷素空间异质性最佳插值方法。对比7种插值方法,在空间插值平滑性方面,LPI、GPI、GWR、GWRK表现较好,在插值速度方面,IDW、RBF、LPI、GPI、OK较快,GWR和GWRK方法运算复杂、速度较慢。IDW、RBF等6种方法与OK相比,根据RI判定,GWRK方法提高了磷素空间分布模拟精度,IDW、GPI和LPI方法降低了磷素空间分布模拟精度,RBF方法在提高磷素空间分布模拟精度上表现不一致。采样密度会影响预测结果,对于本文而言,不论采样密度如何,GWRK方法均方根误差(RMSE)均最低,为最佳插值方法,而RBF方法是在采样密度较低时一种可选方法。GWRK法在本文是最佳的插值方法,但其结果会受到辅助变量多少和各变量之间是否存在共线性的影响。In the late 1990s,the"Dryland to Paddy"project was implemented in the Sanjiang Plain.After planting rice in the Albic soil,the barrier soil layer turns into a favorable soil layer,the low-yield soil becomes high-yield soil,and the Albic soil phosphorus pool increases.After flooding,the availability of phosphorus(closed storage phosphorus[O-P]and iron-bound phosphorus[Fe-P])increases with the decrease in soil redox potential(Eh)and the increase in pH,which substantially affects soil phosphorus heterogeneity.Therefore,we urgently need an optimal interpolation method to improve the prediction accuracy of total phosphorus in the Albic soil of typical irrigation areas of the Sanjiang Plain.This will help evaluate the impact of climate change and land use on the soil phosphorus pools and provide a reference for estimating future soil phosphorus pools.This study used the inverse distance weighting(IDW)method,radial basis function(RBF),ordinary Kriging(OK),global polynomial method(GPI),local polynomial method(LPI),geographic weighted regression(GWR),and geographic weighting regression to Kriging(GWRK)to predict the distribution of soil phosphorus in the Bawusan,Qiliqin,and Daxing irrigation areas of the Sanjiang Plain.The cross-validation method was used to obtain the mean error(ME),root mean square error(RMSE),and relative improvement(RI)to compare the accuracies of the various methods to determine the best interpolation method for assessing the spatial heterogeneity of phosphorus in the same soil type with different sampling densities.Based on the assumptions of regression analysis,this study incorporated 24 environmental variables for exploratory regression analysis,including elevation,pH,organic matter,exchangeable sodium,total nitrogen,available phosphorus,available copper,cultivated layer bulk density,and rice yield.According to the regression results,the auxiliary variables that were significantly correlated with phosphorus were selected for least square analysis.Finally,exchangeable sodium,cation exchange capacit

关 键 词：白浆土 磷素 插值方法 交叉验证 空间自相关性 适用性 

分 类 号：S159.2[农业科学—土壤学]