Mapping the organic carbon stocks of surface soils using local spatial interpolator

Abstract

The largest uncertainties are associated with estimating the soil organic carbon (SOC) stock because of natural soil variability and data scarcity. Thus, a local spatial geostatistical hybrid approach, the geographically weighted regression kriging (GWRK), was used in the present study to overcome some of these uncertainties. This study was designed to estimate the SOC stock (kg C m⁻²) for the surface 0 to 15 cm depth using the state of Pennsylvania as the study region. A total of 920 soil profiles were extracted from the National Soil Survey Center database and were divided into calibration (80%) and validation (20%) periods. Some soil parameters that include clay content, bulk density (ρ_b), total nitrogen (TN) content, pH, Ca²⁺, Na⁺, extractable acidity (EXACID), and cation exchange capacity (CEC) were used as covariates for estimating the SOC stock. These covariates exhibited spatial autocorrelation (Moran's Index, I = 0.62 to 0.89). Further, residuals of geographically weighted regression were spatially autocorrelated, and hence support the use of the GWRK approach. Validation results concluded that the performance of the GWRK approach was the best with the lowest values of root mean square error, mean estimation error and mean absolute estimation error. The estimated SOC stock for the surface 0 to 15 cm depth ranged from 1.41 to 3.94 kg m⁻². Results from this study show that the GWRK captures spatial dependent relationships, and addresses spatial non-stationarity issues, hence this approach improves the estimations of SOC stock.