Effects of pollution caused by seawater intrusion into groundwater in coastal aquifers cannot be ignored. Identification of areas exposed to this pollution by preparing vulnerability maps is one way of preventing aquifer pollution. In its primary section, the present study compared three different index ranking methods of DRASTIC, GALDIT, and SINTACS to select an optimal model for determining vulnerability of the Gharesoo-Gorgan Rood coastal aquifer. Initial results led to selection of the GALDIT model for vulnerability assessment of the selected coastal aquifer. Since this type of models use a rating system, the model must be modified and optimized in various regions to show the vulnerable areas more accurately. In the next step, and for the first time, the ratings in this index were modified using the Wilcoxon nonparametric statistical method and its weights were optimized employing particle swarm optimization (PSO) and single-parameter sensitivity analysis (SPSA) methods. Finally, in order to select the best hybrid model, the total dissolved solids (TDS) parameter was used to determine correlation coefficients. Results indicated that the GALDT model modified by the Wilcoxon-PSO method has the strongest correlation (0.77) with the TDS parameter. Moreover, the correlations of the Wilcoxon-GALDIT and Wilcoxon-SPSA models were 0.66 and 0.73, respectively. Final results of the Wilcoxon-PSO model revealed that the northwestern and western areas of the study region needed considerable protection against pollution. In general, we can conclude that by combining statistical, mathematical, and metaheuristic methods, we can obtain more accurate results for preparing vulnerability maps.
Keywords: DRASTIC; GALDIT; GIS; Gharesoo-Gorgan Rood; PSO; SINTACS; Wilcoxon.