Modelling the effects of tides and storm surges on coastal aquifers using a coupled surface-subsurface approach

Coastal aquifers are complex hydrologic systems because many physical processes interact: (i) variably saturated flow, (ii) spatial-temporal fluid density variations, (iii) tidal fluctuations, (iv) storm surges overtopping dykes, and (v) surface runoff of storm water. The HydroGeoSphere model is used to numerically simulate coastal flow dynamics, assuming a fully coupled surface-subsurface approach, accounting for all processes listed above. The diffusive wave approximation of the St. Venant equation is used to describe surface flow. Surface flow and salt transport are fully coupled with subsurficial variably saturated, variable-density flow and salt transport through mathematical terms that represent exchange of fluid mass and solute mass, respectively. Tides and storm surges induce a time-variant head that is applied to nodes of the surface domain. The approach is applied to real cases of tide and storm surge events. Tide simulation results confirm the existence of a recirculating zone, forming beneath the upper part of the intertidal zone. By monitoring the exchange fluid flux rates through the beach, it was found that the major inflow to the aquifer takes place at the upper part of the intertidal zone, which explains the formation of the recirculating zone. The recirculating zone is forming particularly during rising tide. Results from a storm surge simulation show that plume fingers develop below the flooded land surface. Natural remediation by seaward flowing freshwater is relatively slow, such that reducing the salt concentration in the aquifer down to drinking water standards takes up to 10 years.