Study on the Factors Affecting the Humus Horizon Thickness in the Black Soil Region of Liaoning Province, China

Understanding the spatial variability and driving mechanisms of humus horizon thickness (HHT) degradation is crucial for effective soil degradation prevention in black soil regions. The study compared ordinary kriging interpolation (OK), inverse distance weighted interpolation (IDW), and regression kriging interpolation (RK) using mean error (ME), mean absolute error (MAE), root mean square error (RMSE), and relative RMSE to select the most accurate model. Environmental variables were then integrated to predict HHT characteristics. Results indicate that: (1) RK was superior to OK and IDW in characterizing HHT with the smallest ME (11.45), RMSE (14.98), MAE (11.45), and RRMSE (0.44). (2) The average annual temperature (0.29), precipitation (0.27), and digital elevation model (DEM) (0.21) were the primary factors influencing the spatial variability of HHT. (3) The HHT exhibited notable variability, with an increasing trend from the southeast towards the central and northern directions, being the thinnest in the southeast. It was thicker in the northeast and southwest regions, thicker but less dense along the southern Bohai coast, thicker yet sporadically distributed in the northwest (especially Chaoyang and Fuxin), and thick with aggregated distribution over a smaller area in the northeastern direction (e.g., Tieling). These findings provide a scientific basis for accurate soil management in Liaoning Province.