Digital soil mapping of soil carbon at the farm scale: A spatial downscaling approach in consideration of measured and uncertain data

Abstract In this paper a spatial downscaling method is explored for generating appropriate farm scale digital soil maps. The digital soil map product to be downscaled is an Australian national extent soil carbon map (100 m grid resolution). Taking into account the associated prediction uncertainties of this map, we used a simulation approach based on Gaussian random fields to generate plausible mapping realisations that were in turn downscaled to 10 m resolution for a farm in North-western NSW, Australia. We were able to derive both a downscaled map of soil carbon and associated prediction variance with this approach. Building further upon this development, we then incorporated a bias correction step into the spatial downscaling procedure which permits the inclusion of field observations as a way to moderate the downscaling results to better reflect actual conditions on the ground. Based on an independent validation dataset, it was found that incorporating field observations increase the concordance correlation coefficient to 0.8 from 0.2. This relatively lower correlation achieved using spatial downscaling alone was due to the national scale mapping for the study area being positively biased in the area of interest. It was found that downscaling that incorporates observational data was marginally better if not comparable to using a point-based digital soil mapping approach. The advantage of spatial downscaling is that it can be implemented in situations of data scarcity. This will be ideal for on farm soil monitoring in situations where detailed soil mapping is initially not available. For example, soil carbon auditing schemes requiring prior soil information for implementation of design-based soil sampling could potentially be universally applied with such a spatial downscaling approach.

[1]  Budiman Minasny,et al.  Farm-Scale Soil Carbon Auditing , 2016 .

[2]  Alfred E. Hartemink,et al.  Total soil organic carbon and carbon sequestration potential in Nigeria , 2016 .

[3]  Ruiliang Pu,et al.  Downscaling Thermal Infrared Radiance for Subpixel Land Surface Temperature Retrieval , 2008, Sensors.

[4]  L. Montanarella,et al.  A map of the topsoil organic carbon content of Europe generated by a generalized additive model , 2015 .

[5]  Laura Poggio,et al.  National scale 3D modelling of soil organic carbon stocks with uncertainty propagation — An example from Scotland , 2014 .

[6]  A. Gobiet,et al.  Multi-variable error correction of regional climate models , 2013, Climatic Change.

[7]  R. Lark,et al.  Geostatistics for Environmental Scientists , 2001 .

[8]  Z. Libohova,et al.  Equal-area spline functions applied to a legacy soil database to create weighted-means maps of soil organic carbon at a continental scale , 2012 .

[9]  Edzer J. Pebesma,et al.  Real-time automatic interpolation of ambient gamma dose rates from the Dutch radioactivity monitoring network , 2009, Comput. Geosci..

[10]  Alain Dassargues,et al.  Conceptual model uncertainty in groundwater modeling: Combining generalized likelihood uncertainty estimation and Bayesian model averaging , 2008 .

[11]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[12]  M. Richardson,et al.  The Radiometric Map of Australia , 2009 .

[13]  B. Minasny,et al.  Digital Soil Map of the World , 2009, Science.

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

[15]  R. Tibshirani,et al.  Generalized additive models for medical research , 1986, Statistical methods in medical research.

[16]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[17]  Budiman Minasny,et al.  A general method for downscaling earth resource information , 2012, Comput. Geosci..

[18]  N. Scott,et al.  Monitoring land-use change effects on soil carbon in New Zealand: quantifying baseline soil carbon stocks. , 2002, Environmental pollution.

[19]  Alfred E. Hartemink,et al.  Digital Mapping of Soil Organic Carbon Contents and Stocks in Denmark , 2014, PloS one.

[20]  M. Déqué,et al.  Frequency of precipitation and temperature extremes over France in an anthropogenic scenario: Model results and statistical correction according to observed values , 2007 .

[21]  Edzer J. Pebesma,et al.  Applied Spatial Data Analysis with R - Second Edition , 2008, Use R!.

[22]  A. Wood,et al.  Simulation of Stationary Gaussian Processes in [0, 1] d , 1994 .

[23]  Dominique Arrouays,et al.  GlobalSoilMap France: High-resolution spatial modelling the soils of France up to two meter depth. , 2016, The Science of the total environment.

[24]  Max Kuhn,et al.  caret: Classification and Regression Training , 2015 .

[25]  Steve Weston,et al.  Foreach Parallel Adaptor for the 'parallel' Package , 2015 .

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

[27]  John M. Antle,et al.  Spatial heterogeneity, contract design, and the efficiency of carbon sequestration policies for agriculture , 2003 .

[28]  Laura Poggio,et al.  Downscaling and correction of regional climate models outputs with a hybrid geostatistical approach , 2015 .

[29]  Dominique Arrouays,et al.  GlobalSoilMap : Basis of the global spatial soil information system , 2014 .

[30]  Budiman Minasny,et al.  Digital soil mapping: A brief history and some lessons , 2016 .

[31]  Roger Bivand,et al.  Bindings for the Geospatial Data Abstraction Library , 2015 .

[32]  Rattan Lal,et al.  The knowns, known unknowns and unknowns of sequestration of soil organic carbon , 2013 .

[33]  S. McNeill,et al.  Development and calibration of a soil carbon inventory model for New Zealand , 2014 .

[34]  Edzer J. Pebesma,et al.  Multivariable geostatistics in S: the gstat package , 2004, Comput. Geosci..

[35]  Nicolai Meinshausen,et al.  Quantile Regression Forests , 2006, J. Mach. Learn. Res..

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

[37]  Mike Grundy,et al.  Soil and landscape grid of Australia. , 2015 .

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

[39]  Budiman Minasny,et al.  Using model averaging to combine soil property rasters from legacy soil maps and from point data , 2014 .

[40]  J. Chilès,et al.  Geostatistics: Modeling Spatial Uncertainty , 1999 .

[41]  X. R. L N K C Z V E A Y H G U T X R L N K C Z V E A Y H [General method]. , 2000, Diabetes & metabolism.

[42]  Budiman Minasny,et al.  Utilizing portable X-ray fluorescence spectrometry for in-field investigation of pedogenesis , 2016 .