Mapping soil organic carbon content over New South Wales, Australia using local regression kriging

Abstract Various digital soil mapping techniques ranging from simple linear models to complex machine learning techniques have been employed for soil organic carbon (SOC) mapping. When SOC mapping over a large region is required, the usual approach has to employ a model calibrated for the whole area. An alternative is to use a series of locally calibrated models to map smaller areas that collectively cover the large region of interest. The accuracy of the SOC products generated by these two approaches can potentially vary. However, performance of whole-area calibrated models versus locally calibrated models in mapping SOC of large extents has seldom been explored in detail, particularly with respect to the type of model being employed. Our study aims to fill this gap by evaluating the SOC prediction performance of three common models, multiple linear regression (MLR), Regression tree model; Cubist and Support Vector Regression (SVR) that are calibrated locally and for the whole study area. This study was carried out using eight identified local areas in New South Wales (NSW), Australia and across the whole state entirely. Every model was calibrated separately for each local area and for the entire state. The local and whole-area models were validated using the same test data set over 50 realizations. In particular, local prediction accuracy of whole-area calibrated models was compared to that of locally calibrated models. The models were tested separately for the standard soil depth layers including 0–5 cm, 5–15 cm, 15–30 cm, 30–60 cm; 60–100 cm. The results show that SVR models have a superior performance out of three tested models for all standardized depth layers. In general the local models outperform the whole-area models for all three tested models with respect to the accuracy of predictions. All models displayed area specific performances proving the importance of inclusion of prevailing local conditions in SOC modelling and mapping. Therefore, we introduce a moving window approach where a hybrid series of locally calibrated models and a whole-area calibrated model can be used against using one calibrated model for the modelling very large mapping extents. Moving window approach provides more accurate results having the lowest error compared to the whole-area model. Also it provides the least biased predictions. Therefore, this novel approach provides a promising way of increasing the efficiency and accuracy of digital soil mapping.

[1]  Cristiano Ballabio,et al.  Spatial prediction of soil properties in temperate mountain regions using support vector regression , 2009 .

[2]  Patrick Bogaert,et al.  Regional‐scale characterization of the geomorphic control of the spatial distribution of soil organic carbon in cropland , 2014 .

[3]  Budiman Minasny,et al.  Using Additional Criteria for Measuring the Quality of Predictions and Their Uncertainties in a Digital Soil Mapping Framework , 2011 .

[4]  Fan-Rui Meng,et al.  Using artificial neural network models to produce soil organic carbon content distribution maps across landscapes , 2010 .

[5]  Sandra L. Berry,et al.  Continental‐scale measurement of the soil organic carbon pool with climatic, edaphic, and biotic controls , 2006 .

[6]  M. B. Beck,et al.  Modelling Change in Environmental Systems , 1995 .

[7]  Ian D. Moore,et al.  Terrain attributes: estimation methods and scale effects , 1993 .

[8]  Baron-Hay Court,et al.  Use of airborne gamma radiometric data for soil mapping , 1996 .

[9]  Rattan Lal,et al.  Predicting the spatial variation of the soil organic carbon pool at a regional scale. , 2010 .

[10]  Nicholas M Holden,et al.  Soil organic carbon across scales , 2015, Global change biology.

[11]  Stéphane Canu,et al.  Advanced Spatial Data Analysis and Modelling with Support Vector Machines , 2000 .

[12]  B. Henderson,et al.  Australia-wide predictions of soil properties using decision trees , 2005 .

[13]  Rosa Francaviglia,et al.  Estimation of soil organic matter by geostatistical methods: Use of auxiliary information in agricultural and environmental assessment , 2014 .

[14]  Philippe Lagacherie,et al.  Digital Soil Mapping: A State of the Art , 2008 .

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

[16]  Bernhard Schölkopf,et al.  A tutorial on support vector regression , 2004, Stat. Comput..

[17]  Robert Tibshirani,et al.  The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd Edition , 2001, Springer Series in Statistics.

[18]  Budiman Minasny,et al.  High resolution 3D mapping of soil organic carbon in a heterogeneous agricultural landscape , 2014 .

[19]  Damien J. Field,et al.  Digital soil mapping of soil organic carbon stocks under different land use and land cover types in montane ecosystems, Eastern Himalayas , 2014 .

[20]  Budiman Minasny,et al.  Digital mapping of soil carbon , 2013 .

[21]  H. Jenny Factors of Soil Formation: A System of Quantitative Pedology , 2011 .

[22]  Budiman Minasny,et al.  The neuro-m method for fitting neural network parametric pedotransfer functions , 2002 .

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

[24]  Gerard B. M. Heuvelink,et al.  The GlobalSoilMap project specifications , 2014 .

[25]  Ian H. Witten,et al.  Induction of model trees for predicting continuous classes , 1996 .

[26]  J. Gallant,et al.  A multiresolution index of valley bottom flatness for mapping depositional areas , 2003 .

[27]  Budiman Minasny,et al.  Digital mapping of a soil drainage index for irrigated enterprise suitability in Tasmania, Australia , 2014 .

[28]  G. Heuvelink,et al.  Mapping Soil Properties of Africa at 250 m Resolution: Random Forests Significantly Improve Current Predictions , 2015, PloS one.

[29]  Cristiano Ballabio,et al.  Topsoil Organic Carbon Map of Europe , 2014 .

[30]  R. M. Lark,et al.  Sampling and analytical plus subsampling variance components for five soil indicators observed at regional scale , 2009 .

[31]  Alex B. McBratney,et al.  Modelling soil attribute depth functions with equal-area quadratic smoothing splines , 1999 .

[32]  Karin Viergever,et al.  Using knowledge discovery with data mining from the Australian Soil Resource Information System database to inform soil carbon mapping in Australia , 2009 .

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

[34]  Budiman Minasny,et al.  Predicting and mapping the soil available water capacity of Australian wheatbelt , 2014 .

[35]  Budiman Minasny,et al.  Mapping continuous depth functions of soil carbon storage and available water capacity , 2009 .

[36]  G. Heuvelink,et al.  SoilGrids1km — Global Soil Information Based on Automated Mapping , 2014, PloS one.

[37]  G. Heuvelink,et al.  A generic framework for spatial prediction of soil variables based on regression-kriging , 2004 .

[38]  Budiman Minasny,et al.  VESPER 1.5 - spatial prediction software for precision agriculture. , 2002 .

[39]  Gangcai Liu,et al.  Spatial prediction of soil organic matter content integrating artificial neural network and ordinary kriging in Tibetan Plateau , 2014 .

[40]  Vincent Chaplot,et al.  Estimating carbon stocks at a regional level using soil information and easily accessible auxiliary variables. , 2010 .

[41]  Geoff Holmes,et al.  Generating Rule Sets from Model Trees , 1999, Australian Joint Conference on Artificial Intelligence.

[42]  Simon E. Cook,et al.  Use of airborne gamma radiometric data for soil mapping , 1996 .

[43]  Gerard B. M. Heuvelink,et al.  Small scale digital soil mapping in Southeastern Kenya , 2008 .

[44]  David G. Rossiter,et al.  Soil organic carbon stocks in the Limpopo National Park, Mozambique: Amount, spatial distribution and uncertainty , 2014 .

[45]  Budiman Minasny,et al.  Mapping and comparing the distribution of soil carbon under cropping and grazing management practices in Narrabri, north-west New South Wales , 2010 .

[46]  D. W. Reeves The role of soil organic matter in maintaining soil quality in continuous cropping systems , 1997 .

[47]  M. Wiesmeier,et al.  Digital mapping of soil organic matter stocks using Random Forest modeling in a semi-arid steppe ecosystem , 2011, Plant and Soil.

[48]  L. Lin,et al.  A concordance correlation coefficient to evaluate reproducibility. , 1989, Biometrics.

[49]  Laura Poggio,et al.  Regional scale mapping of soil properties and their uncertainty with a large number of satellite-derived covariates , 2013 .

[50]  Budiman Minasny,et al.  Analysis and prediction of soil properties using local regression-kriging , 2012 .

[51]  Budiman Minasny,et al.  Quantitative models for pedogenesis — A review , 2008 .

[52]  R. Webster,et al.  Baseline map of organic carbon in Australian soil to support national carbon accounting and monitoring under climate change , 2014, Global Change Biology.

[53]  Budiman Minasny,et al.  Constructing a soil class map of Denmark based on the FAO legend using digital techniques , 2014 .

[54]  Markus Neteler,et al.  Advanced Spatial Data Analysis With GRASS GIS 7 , 2013 .

[55]  C. Nickerson A note on a concordance correlation coefficient to evaluate reproducibility , 1997 .

[56]  Brian K. Slater,et al.  Soil Series Mapping By Knowledge Discovery from an Ohio County Soil Map , 2013 .

[57]  Inakwu O. A. Odeh,et al.  Catchment scale mapping of measureable soil organic carbon fractions , 2014 .

[58]  Maohua Wang,et al.  Support vector machines regression and modeling of greenhouse environment , 2009 .