Excessive exploitation of groundwater resources can increase the concentration of pollutants in addition to the progressive drawdown of groundwater table. In this research, to achieve aquifer quantitative and qualitative (QQ) sustainable development, an optimal scenario for withdrawing from operation wells is proposed. At the frst step, the aquifer QQ simulation was carried out with the GMS model. The developed code in MATLAB2018b in the second step provides the link between the simulation and the NSGA-II optimization tools. In the third step, a multi-objective coupled optimization-simulation model based on GMS and NSGA-II developed. Finally, optimal scenario was chosen based on applying the multiple criteria decision-making (MCDM) and Berda Aggregation Method (BAM). The results show that reducing the current withdrawal rate to 51.55% can establish the QQ stability of the aquifer. This decrease in groundwater abstraction has led to a 4.6 m increase in groundwater level (GWL) over 3 years (average 19 cm per month). The spatial and temporal distribution of nitrate concentration after applying the optimal discharge of wells shows the nitrate concentration in central and eastern parts of the aquifer has greatly reduced. Developed sustainable management model can be used to provide a real operation planning of wells to improvement of the QQ status of groundwater in each unconfned aquifer.
