Spatial multi-scale variability of soil nutrients in relation to environmental factors in a typical agricultural region, eastern China.

Identifying the sources of spatial variability and deficiency risk of soil nutrients is a crucial issue for soil and agriculture management. A total of 1247 topsoil samples (0-20 cm) were collected at the nodes of a 2×2 km grid in Rizhao City and the contents of soil organic carbon (OC), total nitrogen (TN), and total phosphorus (TP) were determined. Factorial kriging analysis (FKA), stepwise multiple regression, and indicator kriging (IK) were appled to investigate the scale dependent correlations among soil nutrients, identify the sources of spatial variability at each spatial scale, and delineate the potential risk of soil nutrient deficiency. Linear model of co-regionalization (LMC) fitting indicated that the presence of multi-scale variation was comprised of nugget effect, an exponential structure with a range of 12 km (local scale), and a spherical structure with a range of 84 km (regional scale). The short-range variation of OC and TN was mainly dominated by land use types, and TP was controlled by terrain. At long-range scale, spatial variation of OC, TN, and TP was dominated by parent material. Indicator kriging maps depicted the probability of soil nutrient deficiency compared with the background values in eastern Shandong province. The high deficiency risk area of all nutrient integration was mainly located in eastern and northwestern parts.