Effects of aquifer geometry on seawater intrusion in annulus segment island aquifers

Seawater intrusion in island aquifers was considered analytically, specifically for annulus segment aquifers (ASAs), i.e., aquifers that (in plan) have the shape of an annulus segment. Based on the Ghijben–Herzberg and hillslope-storage Boussinesq equations, analytical solutions were derived for steady-state seawater intrusion in ASAs, with a focus on the freshwater–seawater interface and its corresponding watertable elevation. Predictions of the analytical solutions compared well with experimental data, and so they were employed to investigate the effects of aquifer geometry on seawater intrusion in island aquifers. Three different ASA geometries were compared: convergent (smaller side is facing the lagoon, larger side is the internal no-flow boundary and flow converges towards the lagoon), rectangular and divergent (smaller side is the internal no-flow boundary, larger side is facing the sea and flow diverges towards the sea). Depending on the aquifer geometry, seawater intrusion was found to vary greatly, such that the assumption of a rectangular aquifer to model an ASA can lead to poor estimates of seawater intrusion. Other factors being equal, compared with rectangular aquifers, seawater intrusion is more extensive, and watertable elevation is lower in divergent aquifers, with the opposite tendency in convergent aquifers. Sensitivity analysis further indicated that the effects of aquifer geometry on seawater intrusion and watertable elevation vary with aquifer width and distance from the circle center to the inner arc (the lagoon boundary for convergent aquifers or the internal no-flow boundary for divergent aquifers). A larger aquifer width and distance from the circle center to the inner arc weaken the effects of aquifer geometry, and hence differences in predictions for the three geometries become less pronounced.