The long-term application of phosphate fertilizers in agricultural production leads to a large accumulation of phosphorus in the soil. When it exceeds a certain limit, phosphorus will migrate to surrounding water bodies through surface runoff and other mechanisms, potentially causing environmental risks such as eutrophication of water bodies and increasing the risk of water source pollution. This study takes Shiyan City, the water resources area of the mid-route of the South-to-North Diversion Project (MSDP), as the study area. Based on 701 sampling points of topsoil, geostatistics and geodetectors were used to explore the spatial heterogeneity and influencing factors of available phosphorus (AP) in the topsoil of the area. The results show that the effective phosphorus content in the topsoil of the study area ranges from 0.30 to 146.00 mg/kg, with an average value of 14.28 mg/kg, showing strong variability characteristics. Geostatistical analysis shows that among all theoretical models, the exponential model has the best fitting effect, with a lump gold effect of 0.447 and a range of 82,000 m. The soil available phosphorus content shows an increasing trend from the Central Valley lowlands to the surrounding mountainous hills. Among them, elevation is the main controlling factor for the spatial variation of available phosphorus in the topsoil, followed by soil types, planting systems, annual precipitation, and organic matter. The non-linear enhancement or dual-factor enhancement among various environmental factors reveals the diversity and complexity of spatial heterogeneity affecting available phosphorus content in cultivated soil. This study could provide scientific references for maintaining ecological security in the water source area of the MSDP, improving the precise management of AP, and enhancing cultivated land quality.
Read full abstract