Incorporating Auxiliary Data of Different Spatial Scales for Spatial Prediction of Soil Nitrogen Using Robust Residual Cokriging (RRCoK)

Auxiliary data has usually been incorporated into geostatistics for high-accuracy spatial prediction. Due to the different spatial scales, category and point auxiliary data have rarely been incorporated into prediction models together. Moreover, traditionally used geostatistical models are usually sensitive to outliers. This study first quantified the land-use type (LUT) effect on soil total nitrogen (TN) in Hanchuan County, China. Next, the relationship between soil TN and the auxiliary soil organic matter (SOM) was explored. Then, robust residual cokriging (RRCoK) with LUTs was proposed for the spatial prediction of soil TN. Finally, its spatial prediction accuracy was compared with that of ordinary kriging (OK), robust cokriging (RCoK), and robust residual kriging (RRK). Results show that: (i) both LUT and SOM are closely related to soil TN; (ii) by incorporating SOM, the relative improvement accuracy of RCoK over OK was 29.41%; (iii) by incorporating LUTs, the relative improvement accuracy of RRK over OK was 33.33%; (iv) RRCoK obtained the highest spatial prediction accuracy (RI = 43.14%). It is concluded that the recommended method, RRCoK, can effectively incorporate category and point auxiliary data together for the high-accuracy spatial prediction of soil properties.