Spatially apportioning the source-oriented ecological risks of soil heavy metals using robust spatial receptor model with land-use data and robust residual kriging

Previous ecological risk assessments were mainly concentration-oriented rather than source-oriented. Moreover, land use is usually related to source emissions but was rarely used to improve the source apportionment accuracy. In this study, the land-use effects of heavy metals (HMs) in surface (0–20 cm) and subsurface (20–40 cm) soils were first explored using ANOVA in a suburb of Changzhou City, China; next, based on robust absolute principal component scores-robust geographically weighted regression (RAPCS/RGWR), this study proposed RAPCS/RGWR with land-use type (RAPCS/RGWR-LUT) and compared its source apportionment accuracy with those of basic RAPCS/RGWR and commonly-used absolute principal component scores/multiple linear regression (APCS/MLR); then, the source-oriented ecological risks were apportioned based on RAPCS/RGWR-LUT and Hakanson potential ecological risk index method; finally, this study proposed robust residual kriging with land-use type (RRK) for spatially predicting the source-oriented ecological risks, and compared its spatial prediction accuracy with those of robust ordinary kriging (ROK) and traditionally-used ordinary kriging (OK). Results showed that: (i) by incorporating land-use effects, RAPCS/RGWR-LUT obtained higher source apportionment accuracy than RAPCS/RGWR and APCS/MLR; (ii) the two most important external input sources of the ecological risks were 'atmospheric deposition' (PERIₛᵤᵣfₐcₑ = 47.11 and PERIₛᵤbₛᵤᵣfₐcₑ = 35.27) and 'agronomic measure' (PERIₛᵤᵣfₐcₑ = 28.93 and PERIₛᵤbₛᵤᵣfₐcₑ = 20.37); (iii) the biggest ecological risk factor was soil Cd (ERₛᵤᵣfₐcₑ = 57.14 and ERₛᵤbₛᵤᵣfₐcₑ = 47.62), which was mainly contributed by 'atmospheric deposition' (ERₛᵤᵣfₐcₑ=33.14 and ERₛᵤbₛᵤᵣfₐcₑ=25.71); (iv) RRK obtained higher spatial prediction accuracy than ROK and OK; (v) the high-risk areas derived from 'atmospheric deposition' were mainly located in the southwest of the study area, and the high-risk areas derived from 'agronomic measure' were scattered in the agricultural land in the north and south of the study area. The above information provided effective spatial decision support for reducing the source-oriented input of the ecological risks of soil HMs in a large-scale area.