Spatial non‐stationary effects of explanatory variables on soil bulk density in the critical zone of the Chinese Loess Plateau

Shaozhen Liu,Yunqiang Wang,Zixuan Zhou,Pingping Zhang,Yali Zhao,Lan Xu,Lijun Qi,Hui Sun,Jingxiong Zhou,Zhisheng An

DOI: https://doi.org/10.1111/ejss.13247

IF: 4.178

2022-05-16

European Journal of Soil Science

Abstract:Bulk density (BD) is one of the most important soil physical properties in hydrological processes of Earth's critical zone. The heterogeneity and related controlling factors of surface BD have received extensive research attention. However, it is unclear whether the effects of drivers on BD are homogeneous or vary spatially. This study quantified topsoil BD with a high‐density sampling strategy (10 m × 10 m; n = 1134) across the ShuangChaGou watershed in the critical zone of the Chinese Loess Plateau (CLP‐CZ). Stepwise multiple linear regression (SMLR), regression kriging (RK), random forest (RF), and geographically weighted regression (GWR) were established as modelling approaches to (i) compare the performance of these four models in explaining BD variation, (ii) assess the main factors driving BD, and (iii) evaluate the spatial non‐stationary relationship between BD and its explanatory variables. GWR explained the BD variation to the same extent as the RK model (35%), outperforming both the SMLR (15%) and RF models (19%). The variation in topsoil BD depended mainly on: soil organic carbon (SOC) content, normalized difference vegetation index (NDVI), plan curvature, and slope length. The spatial distribution of the local R2, regression coefficient magnitude, significance, and sign calculated by GWR revealed that the effects of variables accounting for variation in BD varied spatially. The explanatory power of these variables in terms of BD variation was strongest in northern and northeastern areas of the watershed and weakest in southern areas of the watershed. Both significant negative and positive correlations between BD and silt content, NDVI, and plan curvature were found. This study enriched traditional research linked to the complicated relationships between BD and controlling factors by incorporating spatial non‐stationarity, which is useful for improving the accuracy of hydrological modeling in the CLP‐CZ and other similar regions around the world. This article is protected by copyright. All rights reserved.

soil science