To better understand the spatial dynamics of non-point source (NPS) phosphorus loading with soil property at watershed scale, integrated modeling and soil chemistry is crucial to ensure that the indicator is functioning properly and expressing the spatial interaction at two depths. Developments in distributed modeling have greatly enriched the availability of geospatial data analysis and assess the NPS pollution loading response to soil property over larger area. The 1.5 km-grid soil sampling at two depths was analyzed with eight parameters, which provided detailed spatial and vertical soil data under four main types of landuses. The impacts of landuse conversion and agricultural practice on soil property were firstly identified. Except for the slightly bigger total of potassium (TK) and cadmium (Cr), the other six parameters had larger content in 20-40 cm surface than the top 20 cm surface. The Soil and Water Assessment Tool was employed to simulate the loading of NPS phosphorus. Overlaying with the landuse distribution, it was found that the NPS phosphorus mainly comes from the subbasins dominated with upland and paddy rice. The linear correlations of eight soil parameters at two depths with NPS phosphorus loading in the subbasins of upland and paddy rice were compared, respectively. The correlations of available phosphorus (AP), total phosphorus (TP), total nitrogen (TN) and TK varied in two depths, and also can assess the loading. The soil with lower soil organic carbon (SOC) presented a significant higher risk for NPS phosphorus loading, especially in agricultural area. The Principal Component Analysis showed that the TP and zinc (Zn) in top soil and copper (Cu) and Cr in subsurface can work as indicators. The analysis suggested that the application of soil property indicators is useful for assessing NPS phosphorus loss, which is promising for water safety in agricultural area.
Copyright © 2012 Elsevier B.V. All rights reserved.