Comparison of four numerical models for simulating seepage from salt marsh sediments

A boundary integral equation method (BIEM) model and three differently formulated finite element method (FEM) models were implemented to explore the spatial and temporal patterns in marsh pore water seepage that each generated. The BIEM model is based on the Laplace equation coupled to a dynamic free-surface condition that assumes that, as the water-table changes, the aquifer instantaneously loses or gains an amount of water equal to the change in head times the specific yield. The FEM models all implement a simplified Richards equation that allows gradual desaturation or resaturation and thus flow in both the saturated and unsaturated zones of the aquifer. Two of the FEM models are based on the governing equation for the USGS model SUTRA and thus take into account fluid and aquifer compressibility. One of these was modified to take into account the effect of tidal loading on the total stress, which is assumed to be constant in the derivation of the original version of SUTRA. The third FEM model assumes that neither the fluid or aquifer matrix is compressible so that changes in storage are due solely to changes in saturation. The unmodified SUTRA model generated instantaneous boundary fluxes that were up to two orders of magnitude greater, and spatially more uniform, than those of the other models. The FEM model without compressibility generated spatial and temporal patterns of the boundary fluxes very similar to those produced by the BIEM model. The SUTRA model with the tidal stress modification gave fluxes similar in magnitude to the BIEM and no compressibility models but with distinctly different distributions in space and time. These results indicate that accurate simulation of seepage from marsh soils is highly sensitive to aquifer compressibility and to proper formulation of the effect of tidal loading on the total stress in the aquifer. They also suggest that accurate simulation may require total stress correction not only for tidal loading but for changes in the water table as well. Finally, to aid the development of methods for the measurement of compressibility, we present a schematic, pore-scale model to illustrate the factors that may govern the compressibility of marsh soils.