A coupled model for simulating surface water and groundwater interactions in coastal wetlands

Coastal wetlands are characterized by strong, dynamic interactions between surface water and groundwater. This paper presents a coupled model that simulates interacting surface water and groundwater flow and solute transport processes in these wetlands. The coupled model is based on two existing (sub) models for surface water and groundwater, respectively: ELCIRC (a three-dimensional (3-D) finite-volume/finite-difference model for simulating shallow water flow and solute transport in rivers, estuaries and coastal seas) and SUTRA (a 3-D finite-element/finite-difference model for simulating variably saturated, variable-density fluid flow and solute transport in porous media). Both submodels, using compatible unstructured meshes, are coupled spatially at the common interface between the surface water and groundwater bodies. The surface water level and solute concentrations computed by the ELCIRC model are used to determine the boundary conditions of the SUTRA-based groundwater model at the interface. In turn, the groundwater model provides water and solute fluxes as inputs for the continuity equations of surface water flow and solute transport to account for the mass exchange across the interface. Additionally, flux from the seepage face was routed instantaneously to the nearest surface water cell according to the local sediment surface slope. With an external coupling approach, these two submodels run in parallel using time steps of different sizes. The time step (Delta t(g)) for the groundwater model is set to be larger than that (Delta t(s)) used by the surface water model for computational efficiency: Delta t(g) = M x Delta t(s) where M is an integer greater than 1. Data exchange takes place between the two submodels through a common database at synchronized times (e. g. end of each Delta t(g)). The coupled model was validated against two previously reported experiments on surface water and groundwater interactions in coastal lagoons. The results suggest that the model represents well the interacting surface water and groundwater flow and solute transport processes in the lagoons. Copyright (C) 2011 John Wiley & Sons