Integrated application of multiple indicators and geographic information system-based approaches for comprehensive assessment of environmental impacts of toxic metals-contaminated agricultural soils and vegetables

Conventional monitoring and mapping approaches are laborious, expensive, and time-consuming because they need a large number of data and consequently extensive sampling and experimental operations. Therefore, due to the growing concern about the potential of contamination of soils and agricultural products with heavy metals (HMs), a field experiment was conducted on 77 farm lands in an area of 2300 ha in the southeast of Shiraz (Iran) to investigate the source of metal contamination in the soils and vegetables and to model spatial distribution of HMs (iron, Fe; manganese, Mn; copper, Cu; zinc, Zn; cadmium, Cd; nickel, Ni, and lead, Pb) over the region using geographic information system (GIS) and geostatistical (Ordinary Kriging, OK) approaches and compare the results with deterministic approaches (Inverse Distance Weighting, IDW with different weighting power). Furthermore, some ecological and health risks indices including Pollution index (PI), Nemerow integrated pollution index (NIPI), pollution load index (PLI), degree of contamination (Cdeg), modified contamination degree (mCd), PIaverage and PIvector for soil quality, multi-element contamination (MEC), the probability of toxicity (MERMQ), the potential ecological index (RI), total hazard index (THI) and total carcinogenic risk index (TCR) based on ingestion, inhalation, and dermal exposure pathways for adults and children respectively for analyzing the noncarcinogenic and carcinogenic risks were calculated. Experimental semivariogram of the mentioned HMs were calculated and theoretical models (i.e., exponential, spherical, Gaussian, and linear models) were fitted in order to model their spatial structures and to investigate the most representative models. Moreover, principal component analysis (PCA) and cluster analysis (CA) were used to identify sources of HMs in the soils. Results showed that IDW method was more efficient than the OK approach to estimate the properties and HMs contents in the soils and plants. The estimated daily intake of metals (DIM) values of Pb and Ni exceeded their safe limits. In addition, Cd was the main element responsible for ecological risk. The PIave and PIvector indices showed that soil quality in the study area is not suitable. According to mCd values, the soils classified as ultra-high contaminated for Cu and Cd, extremely high for Zn and Pb, very high, high, and very low degree of contamination for Ni, Mn, and Fe, respectively. 36, 60, and 4 % of the sampling sites h