Comparison of Sharp Interface to Variable Density Models in Pumping Optimisation of Coastal Aquifers

A number of models have been developed to simulate seawater intrusion in coastal aquifers, which differ in the accuracy level and computational demands, based on the approximation level of the application. In this paper, four seawater intrusion models are employed to calculate the optimal pumping rates in a coastal aquifer management problem. The first model considers both fluid flow and solute transport processes and assumes a variable-density transition zone between saltwater and freshwater. The implementation of the model in simulation-optimisation routines is impractical, due to the computational time required for the simulation. The second model neglects the dispersion mechanism and assumes a sharp interface between saltwater and freshwater. The sharp interface model is significantly faster than the variable density model, however, it may introduce errors in the estimation of the seawater intrusion extent. The remaining two models are modifications of the second model, which intent to correct the inaccuracies of the simplified sharp interface approximation. All four models are utilised to simulate an unconfined coastal aquifer with multiple pumping wells and an optimisation method is used to calculate the maximum allowed pumping rates. The optimisation results are then analysed, in order to examine if the three sharp interface models could provide feasible solutions in the area of the variable density optimum, which is considered as a benchmark solution.