A Two-Stage Approach to the Estimation of High-Resolution Soil Organic Carbon Storage with Good Extension Capability

Soil organic carbon storage (SOCS) estimation is a crucial branch of the atmospheric–vegetation–soil carbon cycle study under the background of global climate change. SOCS research has increased worldwide. The objective of this study is to develop a two-stage approach with good extension capability to estimate SOCS. In the first stage, an artificial neural network (ANN) model is adopted to estimate SOCS based on 255 soil samples with five soil layers (20 cm increments to 100 cm) in Luoding, Guangdong Province, China. This method is compared with three common methods: The soil type method (STM), ordinary kriging (OK), and radial basis function (RBF) interpolation. In the second stage, a linear model is introduced to capture the regional differences and further improve the estimation accuracy of the Luoding-based ANN model when extending it to Xinxing, Guangdong Province. This is done after assessing the generalizability of the above four methods with 120 soil samples from Xinxing. The results for the first stage show that the ANN model has much better estimation accuracy than STM, OK, and RBF, with the average root mean square error (RMSE) of the five soil layers decreasing by 0.62–0.90 kg·m −2 , R 2 increasing from 0.54 to 0.65, and the mean absolute error decreasing from 0.32 to 0.42. Moreover, the spatial distribution maps produced by the ANN model are more accurate than those of other methods for describing the overall and local SOCS in detail. The results of the second stage indicate that STM, OK, and RBF have poor generalizability ( R 2 < 0.1), and the R 2 value obtained with ANN method is also 43–56% lower for the five soil layers compared with the estimation accuracy achieved in Luoding. However, the R 2 of the linear models built with the 20% soil samples from Xinxing are 0.23–0.29 higher for the five soil layers. Thus, the ANN model is an effective method for accurately estimating SOCS on a regional scale with a small number of field samples. The linear model could easily extend the ANN model to outside areas where the ANN model was originally developed with a better level of accuracy.