This study addresses the issue of optimal management of aquifers using a mathematical simulation- optimization model which relies on the stability of water quality and quantity, considering salinity. In this research first we developed a hydrological model (SWAT) to estimate recharge rates and its spatiotemporal distribution. Then, groundwater simulation of the basin was simulated and calibrated using MODFLOW 2000 and water quality was simulated and calibrated using MT3DMS. Afterwards, a multi-objective optimization model (MOPSO) and embed simulation models as tools to assess the objective function was carried out in order to produce a simulation-optimization model. Finally, a sustainability index to assess Pareto front's answers and three management scenarios (continuing previous operation, 30% increasing and reduction in previous operation) was developed. The results show that the majority of Pareto optimal answers have more sustainability index than a 30% reduction of operation with the best answer of 0.059. Relatively, the sustainability index of 30% reduction of operation is 0.05.