Local variations in the model parameters can play an important explanatory role in the spatial modeling of soil organic carbon (SOC) stock. Linear regression models assume parameters to be spatially invariant and are unable to account for the spatially varying relationships in the variables. A recently developed approach, geographically weighted regression kriging (GWRK), was used in this study to examine the relationships between environmental variables and SOC stock for the state of Pennsylvania, USA. The specific objectives were to (i) estimate the SOC stock (kgCm−2) to 1.0-m depth, and (ii) compare the GWRK results with those obtained from regression kriging (RK). Data for 878 georeferenced soil profiles, extracted from National Soil Survey Center database, were divided into calibration (n=702) and validation (n=176) datasets. Environmental variables including temperature, precipitation, elevation, slope, geology, land use, and normalized difference vegetation index were explored and included as independent variables to establish the model for estimating the SOC stock. Results using Pennsylvania as a case study conclude that GWRK was the least biased and more accurate compared to RK for estimating the SOC stock based on the lowest root mean square error (2.61 vs. 4.61kgm−2), and high R2 (0.36 vs. 0.23) values. Higher stock was consistent with higher precipitation and cooler temperature of the region. Total SOC stock ranged from 1.12 to 1.18Pg for the soils of Pennsylvania. Forests store the highest SOC stock (64% of the total), followed by croplands (22%), wetlands (2.3%), and shrubs (2%). Results show that GWRK enhances the precision for estimating the SOC stock compared to the RK since the former takes into account the spatial non-stationarity coupled with spatial autocorrelation of the residuals.
Read full abstract