Baseline estimates of soil organic carbon by proximal sensing : Comparing design-based, model-assisted and model-based inference

For baselining and to assess changes in soil organic carbon (C) we need efficient soil sampling designs and methods for measuring C stocks. Conventional analytical methods are time-consuming, expensive and impractical, particularly for measuring at depth. Here we demonstrate the use of proximal soil sensors for estimating the total soil organic C stocks and their accuracies in the 0-10 cm, 0-30 cm and 0-100 cm layers, and for mapping the stocks in each of the three depth layers across 2837 ha of grazing land. Sampling locations were selected by probability sampling, which allowed design-based, model-assisted and model-based estimation of the total organic C stock in the study area. We show that spectroscopic and gamma attenuation sensors can produce accurate measures of soil organic C and bulk density at the sampling locations, in this case every 5 cm to a depth of 1 m. Interpolated data from a mobile multisensor platform were used as covariates in Cubist to map soil organic C. The Cubist map was subsequently used as a covariate in the model-assisted and model-based estimation of the total organic C stock. The design-based, model-assisted and model-based estimates of the total organic C stocks in the study area were similar. However, the variances of the model-assisted and model-based estimates were smaller compared to those of the design-based method. The model-based method produced the smallest variances for all three depth layers. Maps helped to assess variability in the C stock of the study area. The contribution of the spectroscopic model prediction error to our uncertainty about the total soil organic C stocks was relatively small. We found that in soil under unimproved pastures, remnant vegetation and forests there is good rationale for measuring soil organic C beyond the commonly recommended depth of 0-30 cm.

[1]  D. Brus Balanced sampling: A versatile sampling approach for statistical soil surveys , 2015 .

[2]  R. Storn,et al.  Differential Evolution: A Practical Approach to Global Optimization (Natural Computing Series) , 2005 .

[3]  Interpolation and optimal linear prediction , 1989 .

[4]  D. J. Brus,et al.  Sampling for Natural Resource Monitoring , 2006 .

[5]  Budiman Minasny,et al.  A conditioned Latin hypercube method for sampling in the presence of ancillary information , 2006, Comput. Geosci..

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

[7]  Yves Tillé,et al.  Doubly balanced spatial sampling with spreading and restitution of auxiliary totals , 2013 .

[8]  D. J. Brus,et al.  Using regression models in design‐based estimation of spatial means of soil properties , 2000 .

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

[10]  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 .

[11]  R. V. Viscarra Rossel,et al.  Sensing of soil bulk density for more accurate carbon accounting , 2016 .

[12]  C. Braak,et al.  Model-free estimation from spatial samples: A reappraisal of classical sampling theory , 1990 .

[13]  L. A. Stone,et al.  Computer Aided Design of Experiments , 1969 .

[14]  D. Brus,et al.  Design‐based versus model‐based estimates of spatial means: Theory and application in environmental soil science , 1993 .

[15]  Budiman Minasny,et al.  Optimizing stratification and allocation for design-based estimation of spatial means using predictions with error , 2015 .

[16]  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 .

[17]  A. Thomson,et al.  Evaluation of Three Field-Based Methods for Quantifying Soil Carbon , 2013, PloS one.

[18]  R. B. Jackson,et al.  THE VERTICAL DISTRIBUTION OF SOIL ORGANIC CARBON AND ITS RELATION TO CLIMATE AND VEGETATION , 2000 .

[19]  Achim Dobermann,et al.  Sampling optimization based on secondary information and its utilization in soil carbon mapping , 2006 .

[20]  A. Savitzky,et al.  Smoothing and Differentiation of Data by Simplified Least Squares Procedures. , 1964 .

[21]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[22]  R. V. Rossel,et al.  Visible and near infrared spectroscopy in soil science , 2010 .

[23]  Marc Voltz,et al.  A comparison of kriging, cubic splines and classification for predicting soil properties from sample information , 1990 .

[24]  D. J. Brus,et al.  Random sampling or geostatistical modelling? Choosing between design-based and model-based sampling strategies for soil (with discussion) , 1997 .

[25]  Richard Webster,et al.  Estimating temporal change in soil monitoring: I. Statistical theory , 1995 .

[26]  D. Brus,et al.  A comparison of kriging, co-kriging and kriging combined with regression for spatial interpolation of horizon depth with censored observations , 1995 .

[27]  M. Kramer,et al.  Emerging land use practices rapidly increase soil organic matter , 2015, Nature Communications.

[28]  Neil McKenzie,et al.  Proximal Soil Sensing: An Effective Approach for Soil Measurements in Space and Time , 2011 .

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

[30]  George E. Rayment,et al.  Soil Chemical Methods – Australasia , 2011 .

[31]  J. J. de Gruijter,et al.  An R package for spatial coverage sampling and random sampling from compact geographical strata by k-means , 2010, Comput. Geosci..

[32]  Richard Webster,et al.  Predicting soil properties from the Australian soil visible–near infrared spectroscopic database , 2012 .

[33]  Alex B. McBratney,et al.  Multivariate calibration of hyperspectral γ‐ray energy spectra for proximal soil sensing , 2007 .

[34]  Carl-Erik Särndal,et al.  Model Assisted Survey Sampling , 1997 .

[35]  R. Lal,et al.  The Depth Distribution of Soil Organic Carbon in Relation to Land Use and Management and the Potential of Carbon Sequestration in Subsoil Horizons , 2005 .

[36]  J. MacQueen Some methods for classification and analysis of multivariate observations , 1967 .

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

[38]  G. Marsily Quantitative Hydrogeology: Groundwater Hydrology for Engineers , 1986 .

[39]  R. Isbell Australian Soil Classification , 1996 .