Effects of aquitard windows on groundwater fluctuations within a coastal leaky aquifer system: An analytical and experimental study

Abstract

Low-permeability aquitards are important sedimentary units that protect the underlying confined aquifers from contamination by seawater and anthropogenic sources in coastal leaky aquifer systems. Therefore, the integrity of the aquitard is of crucial importance. However, the aquitard integrity may be vulnerable to horizontal heterogeneities in the form of aquitard windows. In this study, an innovative analytical model was developed to describe groundwater flow in a coastal leaky aquifer system with horizontal heterogeneities in hydraulic properties of either the aquitard or the confined aquifer. Consequently, the effects of aquitard windows, which manifest as horizontally local heterogeneities in hydraulic conductivity and/or specific storage, on the tidal wave propagation can be considered in the analysis. Analytical solutions were derived for tidal-head or constant-flux inland boundary conditions. The derivation process avoided previous methodologies that involved the assembly of large-size matrices. The newly developed analytical model can cover multiple existing analytical models that are associated with the assumptions of homogeneity in hydraulic properties of the aquifer and the aquitard, infinite horizontal extent of the coastal leaky aquifer system, impermeability of the aquitard, or ignorance of aquitard storativity. Analytical investigations revealed that high-permeability aquitard windows acting as preferential-flow conduits tended to dampen the tide-induced head propagation within the underlying confined aquifer. Moreover, an aquitard with a higher storativity (but the same permeability) and of considerable size can result in an apparent propagation bias. The newly developed analytical model in this study was further corroborated through sandbox experiments. Both piezometer measurements and parameter estimates indicated the non-negligible effects of high-permeability aquitard windows on the tidal wave propagation. The sandbox investigation also highlighted the importance of geostatistical inverse modeling in characterizing local heterogeneities within the confined aquifer and in detecting the location and size of aquitard windows, when measurement data of tide-induced aquifer head fluctuations were provided.