Modeling and Predictive Mapping of Soil Organic Carbon Density in a Small-Scale Area Using Geographically Weighted Regression Kriging Approach

Different natural environmental variables affect the spatial distribution of soil organic carbon (SOC), which has strong spatial heterogeneity and non-stationarity. Additionally, the soil organic carbon density (SOCD) has strong spatial varying relationships with the environmental factors, and the residuals should keep independent. This is one hard and challenging study in digital soil mapping. This study was designed to explore the different impacts of natural environmental factors and construct spatial prediction models of SOC in the junction region (with an area of 2130.37 km(2)) between Enshi City and Yidu City, Hubei Province, China. Multiple spatial interpolation models, such as stepwise linear regression (STR), geographically weighted regression (GWR), regression kriging (RK), and geographically weighted regression kriging (GWRK), were built using different natural environmental variables (e.g., terrain, environmental, and human factors) as auxiliary variables. The goodness of fit (R-2), root mean square error, and improving accuracy were used to evaluate the predicted results of the spatial interpolation models. Results from Pearson correlation coefficient analysis and STR showed that SOCD was strongly correlated with elevation, topographic position index (TPI), topographic wetness index (TWI), slope, and normalized difference vegetation index (NDVI). GWRK had the highest simulation accuracy, followed by RK, whereas STR was the weakest. A theoretical scientific basis is, therefore, provided for exploring the relationship between SOCD and the corresponding environmental variables as well as for modeling and estimating the regional soil carbon pool.