ORCA
Online Research @ CardiffÂ

# Mathematical modelling of the Germasogeia aquifer

 Benham, Graham, Frolkovic, Peter, Ivanov, Tihomir, Mondal, Raka, Mondal, Sourav, Rottschafer, Vivi and Kaouri, Katerina 2017. Mathematical modelling of the Germasogeia aquifer. Presented at: 125th European Study Group with Industry, Limassol, Cyprus, 5-9 December 2016.

 Preview
PDF - Accepted Post-Print Version

## Abstract

Two challenges related to improving the management of the Germasogeia aquifer were presented to the Study Group by the Cyprus Water Development Department (WDD), the public organisation responsible for managing the water resources in Cyprus. The first challenge was how to optimally recharge the aquifer in order to compensate for the extraction of drinking and irrigation water whilst preventing sea water intrusion. In order to address this challenge we developed model for the water in the aquifer. Note that by exploiting the long, thin nature of the aquifer we only develop two-dimensional models in this work. We first develop a simple model based on Darcy flows for porous media which gives the water table height for given dam seepage rate, recharge and extraction rates; we neglect seawater intrusion. We then use the steady version of this model to develop an optimized recharge strategy with which we can identify minimal recharge required for a desired extracted water volume such that the minimum prescribed water table is respected. We explore 4 different scenarios and we find that in certain cases there can be a considerable reduction in the amount of recharged water compared to the current empirical strategy the Water Development Department is employing, where water is recharged and extracted in equal proportions. To incorporate the effects of seawater intrusion, which can be very damaging to the water quality, we next develop transient twodimensional models of saturated-unsaturated groundwater flow and solve them numerically using the open source software SUTRASuite and the commercial package ANSYS FLUENT; the position of the water table and the seawaterfreshwater interface are determined for various extraction/recharge strategies. Data from the WDD are used in some of the simulations. The second important challenge we were asked to look at was to predict the transport of pollutants in the aquifer in the case of an accidental leakage. An advection-diffusion equation for the contaminant concentration is introduced and simulations are undertaken using the commercial package COMSOL. The concentration profiles of the contaminant are studied and we find that the effect of contamination varies depending on where the contamination site is; the closer the contamination site is to the dam, the larger the extent of contamination will be.

Item Type: Conference or Workshop Item (Paper) Unpublished Mathematics 3 July 2018 03 Jul 2018 11:16 http://orca.cardiff.ac.uk/id/eprint/111561

### Actions (repository staff only)

 Edit Item