Optimization design of quality monitoring network of Urmia plain using genetic algorithm and vulnerability map

Abstract Contamination and seawater intrusion are esteemed as chief issues of coastal aquifers resulting from unscientific utilization, non-standard monitoring network and improper management. In the present study, an optimum monitoring network with appropriate numbers and standard spatial distribution was designed based on a vulnerability map of the Urmia coastal aquifer. A vulnerability map was extracted using a modified GALDIT-iP model and searching optimum network was done based on genetic algorithm (GA). The maximum value of the correlation between electrical conductivity (EC) and vulnerability index, the minimum number of monitoring wells and the highest value of Nash-Sutcliff were utilized for the simultaneous optimization model. The W-weighting coefficient was considered for economical goals and three targets were defined in a general objective function. The results showed that the W-weighting coefficient has a significant effect to determine optimal solution, and the best weighting was opted considering the most optimal response based on the precise of the monitoring network and vulnerability index. An acceptable optimization and validation process was obtained with the predictions of the validation results. For W = 1, the final value of the objective function was obtained 1.791 using 91 wells with a correlation coefficient of 0.935 and Nash-Sutcliff of 0.979, resulting in the appropriate spatial distribution of the wells and reduction of 18 wells from the existing monitoring wells.


Introduction
Underground water performs a fundamental role to supply domestic, agricultural and industrial water requirements.On the other hand, the development of human societies besides CONTACT Mahdi Majedi-Asl mehdi.majedi@gmail.com,majedi@maragheh.ac.ir increasing agricultural and industrial activities plays a major role in environmental pollution, especially in the water issue.Meanwhile, underground water resources are highly exposed to pollutant and vulnerability.Therefore, it is more crucial to predict vulnerability precisely and scientifically to achieve an authentic output.Many methods have been developed to evaluate vulnerability of aquifer as process-based methods, statistical methods and overlay/index methods (Tesoriero et al. 1998;Babiker et al. 2005).Data shortage and computational complexity are two main constraints of process-based method (Barbash and Resek 1996).
Insufficient water quality observations, data accuracy and careful selection of spatial variables are three significant limitations of statistical approaches.Although the subjectivity in assigning numerical values to the descriptive entities and relative weights for the different attributes can be counted as overlay-index's drawbacks, but some controlling factors such as rainfall and depth to groundwater which are applied handling the movement of pollutants from the ground surface into the saturated zone can be available over large areas, which makes them suitable for regional scale assessments (Thapinta and Hudak 2003).Some common methodologies of overlay/index method are DRASTIC (Aller et al. 1987), EPIK (Doerfliger et al. 1999), SINTACS (Vrba and Zaporozec 1994), GOD (Foster 1987) and GALDIT (Chachadi and Lobo-Ferreira 2001).As a new method, GALDIT has been utilized by Lobo-Ferreira et al. (2005) in Portugal, Pedreira et al. (2015) in Greece, Kallioras et al. (2011) and Kanani et al. (2017) in India, Kura et al. (2015) for Kapas Island in Malaysia, Luoma et al. (2017) in Finland, Tasnim and Tahsin (2016) in the United States of America, Santha Sophiya and Syed (2013), Gontara et al. (2016) and Trabelsi et al. (2016) in Tunisian coast, Nakhaei et al. (2015), Kardan Moghaddam and Javadi (2016), Bordbar et al. (2018Bordbar et al. ( , 2019) ) and Barzegar et al. (2021) in Iran.Implementing management scenarios requires the utilization of valid data extracted and gathered from observation wells at study zone.Ordinary and heuristic methods are two categories of optimization process.The former one is time-consuming and is unable to respond to large and non-linear hybrid systems.The latter has attached more attention of researchers recently because of more efficient and acceptable responses.While the answer has to be affordable in view of economical and temporal aspects, it gets entirely an estimate of the aquifer condition (Baalousha 2010).As an approach of heuristic types, the genetic algorithm (GA) was developed by Holland (1975).The GA opts the most optimal replies among obtained ones to reach the supreme and the most optimized answer.Randomness of wells distribution besides their large numbers lead to reduction of validity and uneconomical situation of aquifer management (Khader and McKee 2014).Therefore, optimization of monitoring network management via number and distribution is performed (Ayvaz and Elc ¸i 2018).Baalousha (2010) used a combined method of vulnerability mapping and a geostatistical data to design the most optimum monitoring network of water quality.
The preference for monitoring of high pollution potential regions was determined by vulnerability mapping and geostatistical approach was utilized to the spatial distribution of monitored regions in different sites.The DRASTIC was applied to a proper vulnerability map of the regions.The Kriging and the GA were used by Chadalavada and Datta (2008) to optimize the quality monitoring network.Dhar and Patil (2012) applied the NSGA-II algorithm to the design of the network and ordinary Kriging approach to display the temporal and spatial density for optimization of the underground quality monitoring network.The NSGA-II opted as a superior and effective approach to design a quality monitoring network.Fisher (2013) optimized the water level monitoring networks in the Eastern Snak River Plain Aquifer using a Kriging-based GA method.The reduction of 10 to 80 wells along with maintaining the spatial pattern of water level of the aquifer was the output of the research.Nakhaei et al. (2014) used the GA to optimize the rate of pumping from the coastal aquifer of Urmia plain.The results showed that using optimized pumping rates than initial ones raises the level of underground water at most regions of the aquifer.Khoramdel et al. (2015) conducted research using the GA with a combination of Kriging and Particle Swarm Optimization (PSO) methods to optimize the monitoring network of Astaneh-Kooch Isfahan plain.The optimization results proved the reduction of observation well numbers from 57 to 42 without data deficit for analyzing.The GA was utilized by Sadatineghad et al. (2019) to re-design of a monitoring network of Koohdasht aquifer.The results were stranded for the best distribution of 28 wells.Ayvaz and Elci (2018) performed research to optimize an aquifer using the GA and optimization approach for two incompatible goals.At the first stage, the Nash-Sutcliff coefficient was used to maximize the conformity of all existing well data and new-designed network data.At the second stage, the minimizing of the wells number was done with regards to cost restrictions.Utterly, these two goals were combined in an objective function.Reduction of wells number without increasing in error and transporting mapping points from less pollution to regions with pollutant density were the research's outputs.Mahmoodpour et al. ( 2021) modified the vulnerability map of the Talesh plain aquifer using DRASTIC and the GA methods to optimize the quality of the monitoring network of underground water.The validation's optimization showed reliable and acceptable results.Nadjla et al. (2021) applied modified GALDIT method for mapping the groundwater vulnerability to saline intrusion in Algeria.Bordbar et al. (2022) used hybrid machine learning models for improving the GALDIT model for coastal aquifer vulnerability mapping.The results showed all machine learning models performed well in assessing the coastal aquifer vulnerability with respect to the GALDIT model.The best result was obtained by the logistic model tree model.Yang et al. (2022) used GALDIT and TOPSIS for seasonal seawater intrusion mapping.They claimed that their method can help planning and management decisions to prevent groundwater pollution in the studied region.The innovation of this research compared to previous studies is the modification and improvement of the vulnerability map based on GALDIT indices with utilization of two parameters namely hydraulic gradient and pumping rate.Another distinguishing feature of this research is the redesign of the monitoring network with the aim of reducing the number of wells based on genetic algorithm along with weighting the objective function, that is, the number of wells.

Case study
Located at the West Azarbaijan province, Urmia plain has an area of 51762 km 2 and has been expanded between the geographical coordinates of 44 14 0 to 47 53 0 east longitude and 35 40 0 to 38 30 0 north latitude.The research zone is one of the Urmia lake basins covering an area of 2166.2 km 2 (4.18% of Urmia lake basin area).The shares of plain and mountainous areas are 748.96and 1249.2 km 2 , respectively.The highest and lowest height of the research zone are 2000 and 1270 m, respectively.Figure 1 displays the location of the study area.A physical and hydrogeological description of the aquifer, including a geological cross section is illustrated in Figure 2.
The agricultural section has the highest water consumption of underground water in Urmia plain.The amount of water withdrawal from underground resources was 367.86 Mm 3 during 2014-2015.The share of the well and spring were 357.32 and 10.54 Mm 3 , respectively.The negative value of water balance illustrates the lower inlet water volume than outlet volume.Qualitative data of underground water included measured EC provided with Water Organization of West Azarbaijan.The wells of number 109 were used to compile observed data during 2014-2015.The spatial distribution of measured EC using IDW method has been presented in Figure 1.

The method of modified GALDIT vulnerability assessment
Six intrinsic physical characteristics of the media are included in GALDIT index to provide vulnerability map.The GALDIT is abbreviated from the following six hydrogeological factors: groundwater occurrence (G), groundwater system hydraulic conductivity (A), water table head (L), distance from shoreline (D), impact of present status of saltwater intrusion (I) and saturated media depth (T) (Chachadi and Lobo-Ferreira 2001).This model is utilized for the coastal and Mediterranean zones.The GALDIT approach has been modified by adding a hydraulic gradient map and pumping rate (GALDIT-iP) to provide a precise vulnerability map with maximum conformity with existing conditions to achieve appropriate management options for maintenance and utilization of the Urmia plain aquifer.
The GALDIT method seems to be the only example of a large-scale profiling approach to assess the vulnerability of coastal aquifers to saltwater intrusion.Hence, creating a more flexible method is the aim of the present study.The large population of coastal areas, the increasing demand for fresh water and the prominence of underground water as the main source of water supply in coastal areas have been the reasons for choosing the pumping rate for this research.Obviously, there is a direct relationship between the water level and the amount of pumping.The reason for choosing hydraulic gradient is its decisive role in coastal aquifers as the most important controlling variable of sea water, which provides the main source of groundwater discharge to the sea.On the other hand, the limited data of the water level collected from piezometers and observation wells in some areas and the distribution of wells in a wide range and the accessibility of their statistics compared to the water level, causes more accuracy and confidence in a region.Sometimes, the presence of defects in the piezometer and the network of observation wells causes a deficiency in the measurement of the water level.On the other hand, the large diameter of the pumping wells and the small storage in the well result in misleading water level measurements.Also, in karst systems, their special channel and fissure system cannot represent the real water level.Therefore, in this research, the pumping rate has been introduced to solve the aforementioned problems for measuring the water level.
The map of vulnerability index of GALDIT-iP was obtained in ArcGIS software using Equation (1), overlapping of raster layers gotten from parameters from Figure 3 with regards to their order and weighting coefficients.

Groundwater occurrence (G)
Seawater intrusion depends on aquifer type.Unconfined, confined and leakage are tree types of aquifers.The highest score belongs to confined one of value 10.The study area is unconfined type based on the Water Organization balance report (Figure 3a).

Aquifer hydraulic conductivity (A)
Hydraulic conductivity is the capability of water transportation through soil or bedrock of aquifer that is created due to the effective porosity of soil particles (Castany 1982).There is a straight relationship between hydraulic conductivity and the penetrability of seawater into the aquifer.The value of hydraulic conductivity was calculated by dividing the transmissivity coefficient by the saturated thickness of the aquifer (Figure 3b).

The difference between underground level and sea level (L)
Level difference plays a significant role to provide required pressure for the movement of seawater into the aquifer.Using Ghyben-Herzberg relationship, for every meter of freshwater withdrawal above the average sea level, the saltwater level rises by 40 m (Figure 3c).

Distance from the shoreline (D)
It is called the air distance between the desired location and the line perpendicular to the shoreline.If the aquifer hydrogeological conditions are susceptible to transfer, it will happen the highest intrusion of seawater into the aquifer when the aquifer distance is short and close to the shoreline (Figure 3d).
2.1.6.Impact of existing status of seawater intrusion in the area (I) While the amount of bicarbonate in groundwater is much higher than in seawater, the most abundant ion in seawater and the least ion in groundwater is chloride.Equation ( 2) has been developed by Revelle (1941) as an index for detecting seawater intrusion into the coastal aquifer and is used as a ranking criterion of the existence of the seawater intrusion effect: Figure 3e illustrates the magnitude of the seawater intrusion map.

Thickness of the mapped aquifer (T)
Aquifer thickness has a remarkable role to estimate the amount of seawater into the coastal zone that is called the thickness of the saturation part of the aquifer.The magnitude of seawater intrusion is directly related to the thickness of the aquifer.The thickness of the aquifer is defined as the distance between the water table up to the aquifer impermeable layer (Figure 3f).Table 1 shows the ranking of this parameter.

Hydraulic gradient (i)
Distance from the shoreline and the level difference of underground water with sea level define the hydraulic gradient.The hydraulic slope of the sea in coastal aquifers provides the main source of groundwater discharge to the sea with overuse of groundwater and the most significant controlling parameter is seawater (Ferguson and Gleeson 2012;Werner et al. 2012;Ketabchi and Ataie-Ashtiani, 2015;Comte et al. 2016;Holding and Allen 2015).Figure 3g displays the map of the hydraulic gradient of the study area.

Pumping rate (P)
The pumping rate can be replaced by the water level in the area of high population and more fresh water consumption which is extracted from underground water resources (Javadi et al. 2012).The pumping wells provide worthy information with the annual pumping rate known due to defects in the network of observation wells in some areas (Figure 3h).

Genetic algorithm
Developed by Holland (1975) and introduced by Goldberg (1989), GA is an evolutionary adaptive method and based on the genetic process of living organisms that is widely used to solve search and optimization problems.The GA is best suited to solve combinatorial optimization problems that cannot be solved using more conventional operational research methods.Hence, it can be applied to large, complex problems that are nonlinear with multiple local optima.In GA, the process starts with a initial population involving set of solutions.Each solution, called chromosome, is encoded in a string.Each chromosome composed of genes which take a quantity of values.
First, the population is randomly generated.A new population is then produced at each iteration of generation by employing genetic operators such as mutation, crossover, reproduction and selection.Each solution set is assessed based on a fitness function, and this procedure is stopped when the fitness value reaches a threshold which can be determined by the user.The main aim of the optimization procedure is to minimize or maximize the fitness.Figure 4 illustrates the flowchart of the GA (Savic et al. 1999).

Modeling process
A database of EC parameters was provided from observation wells in the study area.
Then, the GA was applied to design an optimized monitoring network and determine an efficient wells number and their location.

Database
The number of 109 data were obtained from monitoring wells in the existing study area.
The IDW interpolation approach was utilized to extract the raster layer of the study zone.
The IDW had the least error and the maximum conformity and correlation than other interpolation methods.The obtained raster layer put EC data into the 741 regular cells with the same distance with mesh size 1000 Â 1000 m in the study area.

Optimization model
The defined optimized model follows three objective functions.The first goal illustrates the conformity of the selected monitoring wells with the vulnerability of the modified GALDIT which is obtained from the maximum correlation of vulnerability map and calculated values.The second goal stands for minimizing the number of monitoring wells due to their uneconomical situation.The third target is minimizing the difference between calculated values of EC and measured ones in all aquifer zones using the Nash-Sutcliff index to achieve a precise monitoring network.The equations of the optimization model were defined as follows: where x i is the vulnerability index of point ith, x is the mean value of the vulnerability index, y i is the qualitative parameter at point ith, y is the mean value of the qualitative parameter, m and n are number of selected wells via optimization algorithm and the total number of wells at the study zone, respectively; y s is predicted value of the qualitative parameter and y s is the mean value of the qualitative parameter.The Z n , with the best value þ1, is an objective optimization function that changes from À1 to þ1 based on three functions Z 1 , Z 2 and Z 3 to maximize the correlation between vulnerability index and qualitative parameter of underground water at the first function.Despite the different distinct nature of each three functions, their ideal values equal one; hence the w-weighting coefficient was used in Z n function for explaining the importance of Z 2 compared to two other functions, that is, Z 1 and Z 3 .The different values of the weighting coefficient lead to optimized answers that are calculated using sensitivity analysis of w.

Simulation
The restrictions and the values of objective functions for the optimization problem are determined using simulation.For the first stage, the values of EC are calculated using the layer of the IDW interpolation method all around the study area.Then, the vulnerability index of the calculated points is utilized to compute the value of correlation (Z 1 ).The Nash-Sutcliff coefficient (Z 3 ) is computed using interpolated values of EC and measured values of potential database points.Utterly, the values of objective function Z 2 are calculated based on the known number of monitoring wells network.

Optimization process
The crossover probability (P cross ) is utilized for the obtained answers for each step at the GA.High values were assigned to P cross , that is, between 0.5 and 1, to eliminate repeatable answers.To diversify the number of obtained answers, mutation probability (P mut ) from 0.05 to 0.25 was considered (Fisher, 2013;Ayvaz and Kentel 2015).These phases continue to meet stop conditions.In the present study, a trial and error process was performed using MATLAB software, and the values of population, crossover probability, mutation probability and the maximum values of repetitions were determined as 50, 0.95, 0.05 and 1000, respectively.

Analysis of the vulnerability map
Data processing was performed in ArcGIS software and the layers of classified vulnerability effective parameters were extracted (Figure 3).Table 1 was utilized for weighting and combination of layers for each parameter; finally, a vulnerability map was provided using Eq. ( 1) and overlap of the weighted layers (Figure 5).The intended index was illustratively categorized which has been presented in Table 2.The computed values explain the vulnerability of the aquifer relative seawater intrusion so that larger numbers indicate greater aquifer vulnerability.The vulnerability index of the whole plain was estimated between 3 and 8 and was classified in five groups: 3-3.75, 3.75-4.5, 4.5-5.25, 5.25-6 and 6-8.The counterpart vulnerability level for the aforementioned classification was very low, low, average, high and very high, respectively.According to Figure 4, while the northern, central part and coastal margin of Urmia plain have high to very high vulnerability index, in the southern zones and far from the coastal margin, the severity of the vulnerability has decreased.

Analysis of the GA output
The analysis of the GA is performed based on different W-weighting values (Table 3).According to Table 3, while the objective function was computed as 1.927 for W ¼ 0.01 with correlation coefficient 0.952 and Nash-Sutcliff coefficient 0.978, the number of monitoring wells was determined as 234.It has 125 wells more than the exiting monitoring network.The value of the objective function was computed as 1.887 for W ¼ 0.2, with the maximum values of Nash-Sutcliff function of 0.9843 and well numbers of 178 (69 wells more than the existing monitoring network).Using 91 wells with a correlation coefficient of 0.935 and Nash-Sutcliff of 0.979, the value of the objective function was calculated as 1.791 for W ¼ 1.0.The number of wells decreased from 109 to 91.There is no significant difference between well number at the monitoring network for a tenfold increase in weight (W ¼ 1 and W ¼ 10).But a remarkable decrease is recognizable for a  tenfold rise from W ¼ 0.1 to W ¼ 1.0.Although the objective function Z n changes from 1.791 to À99.617 (a distinguishing difference) for W ¼ 1 to W ¼ 1000, the difference in the number of selected wells is not large.The lowest well number relates to W ¼ 1000, but this does not meet the best answer because the values of functions Z 1 and Z 3 determine various locations for wells that optimize the final objective function.Also, for similar values of W, the less vales of the wells number, has a more negative effect on the Z 2 and Z 3 .Hence, the optimization model during the optimization process tends to keep the wells number constant and search for better well positions to modify the values of Z 1 and Z 3 .There is no more difference between well number for W ¼ 1 and W ¼ 10, but according to the final objective function and various Z 1 and Z 3 functions for these two weights, different locations are obtained for opted wells with a very impressive ultimate result.Figure 5 presents the conformity of the vulnerability map and interpolated map of EC values along with proposed locations of wells for W ¼ 1 and W ¼ 10.The spatial distribution of wells reveals that the location of proposed wells has been located in regions with critical values in terms of vulnerability and EC.A comparison of the spatial distribution of monitoring wells according to the raster map of the presented EC values and the optimized situation in Figure 6 shows that in areas where aquifer salinity is higher and conditions are more critical, the number of selected wells is higher than the current situation and on the contrary in areas that are less at risk, the number of selected wells is less than the current situation.It can be concluded that although the wells number is less in W ¼ 10 and is more economical, according to the comparison of selected locations, the accuracy of the selected positions is higher for W ¼ 1.0.According to the optimization model, it can be claimed that the ultimate monitoring network is opted by examining the appropriate weight, optimal wells number and correct spatial distribution of wells.Less wells number for high weighting is economical but the accuracy of qualitative results declines due to the reduction of correlation function (Z 1 ) and Nash-Sutcliff function (Z 3 ).

Results verification
The verification of the results can be assessed for the accuracy of the monitoring network for various time series.In the present study, the values of EC in 2016 were used to verify the monitoring network extracted from 2015. Figure 7 illustrates the verification of the results using the correlation coefficient.The values of the correlation coefficients for W ¼ 0.1, 0.5, 0.6, 0.7, 0.8, 0.9 and 100 during verification period are higher than the optimization period, indicating results precisely.Verification of results in terms of the Nash-Sutcliff coefficient for 2015 and 2016 is shown in Figure 8. Moriasi et al. (2007) presented the following classification of Nash-Sutcliff coefficient: 0.75-1 (very good), 0.65 À 0.75 (good), 0.5-0.65 (acceptable), <0.5 (unacceptable).According to Figure 8, all weights are included in the very good category.
The values of root mean square error (RMSE), bias percent (PBIAS) and standard deviation (S) are calculated using the following equations: where n is the number of selected wells, y i is the value of the qualitative parameter at the point ith, y s is the predicted value of the qualitative parameter, y s is the mean value of the qualitative parameter.The results of Eqs. ( 7) to ( 9) have been presented in Table 4. Singh et al. (2004) named an answer as a solution with less error if the value of the RMSE is less than half of S. For all the data in Table 4, both verification and optimization periods for all weights have above-mentioned condition.Santhi et al. (2001) and Van Liew et al. (2003) claimed that the values of r(Z 1 ) more than 0.5 are acceptable.Thus, all weights in Table 4 are allowable precise.PBIAS index shows the average deviation of the predicted values from the observed values in terms of percentage.The closer the value of PBIAS to zero indicates less difference and more accurate answer.Moriasi et al. (2007) classified PBIAS values less than 10% as very good estimations.
Table 5 shows the changes of mean EC values for various weights during optimization and the existing monitoring network.According to Table 5, there is an adverse relationship between different weights and wells number so that the decrease in the number of monitoring wells corresponds to the increase in the EC average.In other words, the optimization model firstly opts for high values of EC and the maximum correlation coefficient with vulnerability index; then the points of low values of EC and correlation coefficient are chosen.Therefore, it is better to choose a weight with high average values of EC.

Conclusion
Inappropriate spatial distribution of monitoring wells is one of the fundamental issues for monitoring networks.Incorrect results are obtained due to incomplete information in critical points and/or excess information in non-critical points because of monitoring wells distribution.Also, the extra number of monitoring wells is not economically viable.The optimization process is known as an appropriate solution for designing a monitoring network.The mission of the present study is designing an optimized monitoring network at the study area.
The time-series 2015 and 2016 were used for optimization and verification processes, respectively.GALDIT-IP index was developed and utilized to obtain vulnerability map of the plain.Increasing correleation coefficient from 0.42 to 0.66 demonstrated the admissible correction of the new monitoring network.The economical objective function was applied at the main objective function to obtain optimal answers based on the direct relationship between wells number and the cost of monitoring network design.The number of wells was reduced by 17% in the presented study.Jason (2013), Mahmoudpour et al. (2021) and Sadatinejad et al. (2019) performed an analysis using the GA method to optimize well numbers.Their results led to reduction of well numbers besides acceptable precision.The application of three targets into the one objective function using w-weighting coefficients illustrated the relative importance of goals.The results of optimization evaluation in different weights showed the importance of w in determining the optimal solution.The value of w should be opted based on the most optimized answer with the most appropriate relationship among all objective functions.RMSE, PBIAS and standard deviation were utilized to determine the optimal solution.Utterly, a comparison between mean values of computed and measured EC was used to identify more pollutant areas.From this study, the following conclusions can be reached: Based on the different weights of W, the results of the genetic algorithm were analyzed.Choosing the optimal solution highly depends on the weighting factor W. Therefore, the value of W should be chosen according to the most balanced solution.The most balanced and optimal solution means an acceptable relationship between the cost and the spatial distribution of wells in the region.Different weights lead to the creation of different optimal quality monitoring networks.In order to design the spatial distribution of monitoring wells, it is better to choose important and critical areas such as points with poorer quality, because these data better reflect the quality of underground water in more polluted areas.The design of the underground water quality monitoring network should be periodically optimized because the successive evaluations of the monitoring network play an important role in determining the long-term evaluation of the underground water quality and the factors affecting it.This can be of great help in planning and applying methods to improve the quality of underground water.

Disclosure Statement
No potential conflict of interest was reported by the author(s).

Figure 1 .
Figure 1.General location of the study area with geological map and spatial distribution of EC values.

Figure 2 .
Figure 2. Geological section map of Urmia plain aquifer.

Figure 3 .
Figure 3. (a) Type of aquifer; (b) Hydraulic conductivity map; (c) Level difference map of underground water with sea level; (d) Distance map from the shoreline; (e) Developing map of seawater intrusion; (f) Aquifer thickness map; (g) Hydraulic gradient map, (h): Pumping rate map.

Figure 4 .
Figure 4. Overview process of the GA.

Figure 5 .
Figure 5.The vulnerability map of coastal aquifer of Urmia plain using modified GALDIT-iP.

Figure 6 .
Figure 6.Interpolation map of EC matching for (a) w ¼ 1 and (b) w ¼ 10 and vulnerability map matching for (c) w ¼ 1 and (d) w ¼ 10.

Figure 7 .
Figure 7. Validation of correlation coefficient of optimal wells for different w values.

Figure 8 .
Figure 8. Validation of Nash-Sutcliffe coefficient model of optimal wells for different w values.

Table 2 .
Classification and percentages of vulnerable areas obtained from GALDIT-iP index.

Table 3 .
Results obtained for different w values.

Table 4 .
Validation outcomes for different weights w values.

Table 5 .
The change of mean EC values in monitored groundwater quality (unit: mSiemens/cm).