Combining spatial autocorrelation with artificial intelligence models to estimate spatial distribution and risks of heavy metal pollution in agricultural soils

Abstract

Information on spatial distribution and potential sources of heavy metals in agricultural lands is very important for human health and food safety. In this study, pollution degree of lead (Pb), cadmium (Cd), and nickel (Ni) in Yüksekova Plain, located on the border in the southeastern part of Turkey, was evaluated by geoaccumulation index (Igeo), modified contamination factor (mCdeg), and Nemerow pollution index (PINemerow) combined with spatial autocorrelation using deep learning algorithms. A total of 304 soil samples were collected from two different depths (0-20 and 20-40cm) in the study area, which covered 17.5 thousand ha land. Covariates were determined for spatial distribution models of Pb, Cd, and Ni by factor analysis (FA). Spatial distribution models for surface soils were developed using pedovariables (silt, sand, clay lime, organic matter, electrical conductivity, pH, Ca, and Na) determined by the FA and Igeo and mCdeg values by artificial neural network (ANN) and adaptive neuro-fuzzy inference system (ANFIS) models. The estimation success of models for different depths was assessed by root mean square error (RMSE), mean absolute percent error (MAPE), and Taylor diagrams. The RMSE and MAPE values showed a strong correlation between heavy metal contents and the covariates. The RMSE values of ANN-Ni0-20, ANN-Ni20-40, ANN-Pb0-20, ANN-Cd0-20, and ANN-Cd20-40 models (0.01240, 0.07257, 0.0039, 0.00045, 0.00044, and 0.04607, respectively) confirmed the success of the models. Likewise, the MAPE values between 0.2 and 8.5% indicated that all models were very good predictors. In addition, the Taylor diagrams showed that the estimation performance of ANFIS and ANN models are compatible. The IgeoNi and IgeoPb values in both models at both depths indicated that strongly to extremely polluted (4-5) areas are quite high in the study area, while the IgeoCd values revealed that unpolluted areas are widespread. The mCdeg index value showed a moderate to high contamination at the first depth, while very high contamination at the second depth in most of the study area. Spatial distribution of PINemerow revealed that moderate pollution (2-3) is common in both soil depths of the study area. The PINemerow of subsurface layer was between 0.91 and 1 (warning limit class) in a small part of the study area. The results showed that vertical mobility of heavy metals is closely related to pedovariables. In addition, the ANN and ANFIS models are capable of exhibiting the heterogeneity in the spatial distribution pattern of high variation in the data. Thus, the locations with extreme contamination have been accurately determined. The pollution indices calculated considering the commonly used international reference values revealed that heavy metal pollution in some part of the study area reached the detrimental levels for human health and food safety. The results suggested that the pollution indices were more successful than simple heavy metal concentrations in interpreting the pollution risk levels. High-resolution spatial information reported in this study can help policy makers and authorities to reduce heavy metal emissions of pollutants or, if possible, to eliminate the pollution.