Spatial prediction of soil organic matter in the presence of different external trends with REML-EBLUP

Abstract The objective of our study was to compare the performance of the empirical best linear unbiased predictor (E-BLUP) with residual maximum likelihood (REML) with that of regression kriging (RK) for predicting soil organic matter (SOM) with the presence of different external drifts. The study was conducted on a 933 km 2 area in Pinggu district of Beijing. Terrain attributes (elevation, slope and topographic wetness index) calculated from DEM were used as external drift variable. The root mean squared errors (RMSE) and the mean squared deviation ratio (MSDR) were used to assess the accuracy of prediction and the goodness of fit of the theoretical estimate of error respectively. RK resulted in both the most and least accurate predictions. REML-EBLUP provided more correct residual variogram models than RK for each trend model. Our results have shown that when the value of adjusted R 2 is greater than 0.45, there is litter difference in the ability to increase the accuracy between REML-EBLUP and RK; and when the value is less than 0.45, the performance of REML-EBLUP is significantly better than RK. It also suggested that topographical data can further improve the accuracy of the spatial predictions of SOM by using RK and REML-EBLUP.

[1]  P. Goovaerts.  Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall , 2000 .

[2]  F. J. Pierce,et al.  Soil Carbon Maps , 2003 .

[3]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.

[4]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[5]  J. Qi,et al.  Total carbon mapping in glacial till soils using near-infrared spectroscopy, Landsat imagery and topographical information , 2007 .

[6]  Sabine Grunwald,et al.  Incorporation of spectral data into multivariate geostatistical models to map soil phosphorus variability in a Florida wetland , 2007 .

[7]  A. McBratney,et al.  Further results on prediction of soil properties from terrain attributes: heterotopic cokriging and regression-kriging , 1995 .

[8]  Pierre Goovaerts,et al.  Fine-resolution mapping of soil organic carbon based on multivariate secondary data , 2006 .

[9]  Alex B. McBratney,et al.  A comparison of prediction methods for the creation of field-extent soil property maps , 2001 .

[10]  N. Cressie,et al.  Statistics for Spatial Data. , 1992 .

[11]  R. M. Lark,et al.  Geostatistical mapping of geomorphic variables in the presence of trend , 2006 .

[12]  W. P. Miller,et al.  Soil variations along a climatic gradient in an Andean agro-ecosystem , 2006 .

[13]  Yusuke Takata,et al.  Spatial prediction of soil organic matter in northern Kazakhstan based on topographic and vegetation information , 2007 .

[14]  Biao Huang,et al.  Temporal and spatial variability of soil organic matter and total nitrogen in an agricultural ecosystem as affected by farming practices , 2007 .

[15]  R. M. Lark,et al.  Reply to “Standardized vs. customary ordinary cokriging…” by A. Papritz , 2008 .

[16]  Jan Seibert,et al.  Topographical influences on soil properties in boreal forests , 2007 .

[17]  R. Lark,et al.  On spatial prediction of soil properties in the presence of a spatial trend: the empirical best linear unbiased predictor (E‐BLUP) with REML , 2006 .

[18]  R. M. Lark,et al.  A comparison of some robust estimators of the variogram for use in soil survey , 2000 .

[19]  Budiman Minasny,et al.  On digital soil mapping , 2003 .

[20]  B. Minasny,et al.  The Matérn function as a general model for soil variograms , 2005 .

[21]  B. Minasny,et al.  Spatial prediction of soil properties using EBLUP with the Matérn covariance function , 2007 .

[22]  Baoguo Li,et al.  Spatial and temporal patterns of soil organic matter in the urban–rural transition zone of Beijing , 2007 .

[23]  Margaret A. Oliver,et al.  The spatial prediction of soil mineral N and potentially available N using elevation , 2005 .

[24]  H. D. Patterson,et al.  Recovery of inter-block information when block sizes are unequal , 1971 .

[25]  Frank Canters,et al.  A multiple regression approach to assess the spatial distribution of Soil Organic Carbon (SOC) at the regional scale (Flanders, Belgium) , 2008 .

[26]  K. Beven,et al.  A physically based, variable contributing area model of basin hydrology , 1979 .

[27]  Ingrid Kögel-Knabner,et al.  Soil organic matter fractions as early indicators for carbon stock changes under different land-use? , 2005 .

[28]  Margaret A. Oliver,et al.  Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood , 2007 .

[29]  R. Reese.  Geostatistics for Environmental Scientists , 2001 .

[30]  Gerard B. M. Heuvelink,et al.  About regression-kriging: From equations to case studies , 2007, Comput. Geosci..

[31]  B. Diekkrüger,et al.  Geostatistical co-regionalization of soil hydraulic properties in a micro-scale catchment using terrain attributes , 2006 .