Optimization of Groundwater Use in the Goksu Delta at Silifke, Turkey

Creative Commons License

Gördü F., Yurtal R. , Motz L. H.

First International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management,, Fas, Morocco, 23 - 25 April 2001, vol.1, no.72, pp.1-10

  • Publication Type: Conference Paper / Full Text
  • Volume: 1
  • City: Fas
  • Country: Morocco
  • Page Numbers: pp.1-10


An optimization model was developed to manage the supplemental use of groundwater in the coastal aquifer subject to saltwater intrusion in the Goksu Delta at Silifke, Turkey. The response of the aquifer system was linked to the optimization model using the response matrix method. A calibrated groundwater simulation model using the SUTRA (Saturated-Unsaturated Transport) code was run to generate aquifer response coefficients at specific well locations. It was assumed that pumping occurs from two wells, and a linear optimization model was constructed under steady-state conditions to maximize the total pumping rates from these two wells subject to water demands and chloride concentration and drawdown limitations. The model was solved for five levels of chloride concentration limits to develop an optimal pumping strategy for the different chloride concentration levels. To verify the optimization model, hydraulic heads and chloride concentrations estimated by the optimization model were compared with those computed by the simulation model. The GAMS (General Algebraic Modeling System) code was used to execute the optimization model. A trade-off curve was constructed at the conclusion of the study that shows the change in maximum pumping rate with respect to the different chloride concentration levels. After plotting the results of both the simulation and optimization models, the predicted hydraulic heads from the optimization model match those generated by the simulation model very closely. However, due to nonlinear effects, the correlation between the chloride concentrations predicted by the optimization model and those calculated by the simulation model is not as good.