Comparison of interpolation methods for the estimation of groundwater contamination in Andimeshk-Shush Plain, Southwest of Iran

Selection of appropriate interpolation methods for the conversion of discrete samples into continuous maps is a controversial issue in the environmental researches. The main objective of this study was to analyze the suitability of three interpolation methods for the discrimination of groundwater with respect to the water quality index (WQI). The groundwater quality data consisted of 17 variables associated with 65 wells located in Andimeshk-Shush Plain. Three spatial interpolation methods including ordinary kriging (OK), empirical Bayesian kriging (EBK), and inverse distance weighting (IDW) were utilized for modeling the groundwater contamination. In addition, different cross-validation indicators were applied to assess the performance of different interpolation methods. The results showed that the performance differed slightly among different methods, although the best performed interpolation method in this study was the empirical Bayesian kriging. Among the interpolation methods, IDW with weighting power of 4 estimated the most contaminated area, while OK estimated the lowest contaminated area. The weighting power of IDW had a significant influence on the estimation, meaning that the estimated contaminated area was increased when a greater weighting power was selected. The subtraction results indicated that there are slightly spatial differences among the contamination assessment results. Results of both standard deviation (SD) and coefficient of variation (CV) also showed that uncertainty was highest in the southern part of the study area, where the distribution of wells were more intensive than that of the northern part.