Beach water table fluctuations due to spring-neap tides: moving boundary effects

Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection-diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring-neap tides (a bichromatic tidal system with the fundamental frequencies wt and wt) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency omega(1) - omega(2), i.e., the spring-neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations. (C) 2000 Elsevier Science Ltd. All rights reserved.

Spring-neap tide-induced beach water table fluctuations in a sloping coastal aquifer

Predictions of water table fluctuations in coastal aquifers are needed for numerous coastal and water resources engineering problems. Most previous investigations have been based on the Boussinesq equation for the case of a vertical beach. In this note an analytical solution based on shallow water expansion for the spring- neap tide- induced water table fluctuations in a coastal aquifer is presented. Unlike most previous investigations, multitidal signals are considered with a sloping coastal aquifer. The new solution is verified by comparing with field observations from Ardeer, Scotland. On the basis of the analytical approximation the influences of higher- order components on water table elevation are examined first. Then, a parametric study has been performed to investigate the effects of the amplitude ratio (lambda), frequency ratio (omega), and phases (delta(1) and delta(2)) on the tide- induced water table fluctuations in a sloping sandy beach.

Numerical modelling of tide-induced beach water table fluctuations

Field studies have shown that the elevation of the beach groundwater table varies with the tide and such variations affect significantly beach erosion or accretion. In this paper, we present a BEM (Boundary Element Method) model for simulating the tidal fluctuation of the beach groundwater table. The model solves the two-dimensional flow equation subject to free and moving boundary conditions, including the seepage dynamics at the beach face. The simulated seepage faces were found to agree with the predictions of a simple model (Turner, 1993). The advantage of the present model is, however, that it can be used with little modification to simulate more complicated cases, e.g., surface recharge from rainfall and drainage in the aquifer may be included (the latter is related to beach dewatering technique). The model also simulated well the field data of Nielsen (1990). In particular, the model replicated three distinct features of local water table fluctuations: steep rising phase versus flat falling phase, amplitude attenuation and phase lagging.

Beach water table fluctuations due to wave run-up: Capillarity effects

High-frequency beach water table fluctuations due to wave run-up and rundown have been observed in the field [Waddell, 1976]. Such fluctuations affect the infiltration/exfiltration process across the beach face and the interstitial oxygenation process in the beach ecosystem. Accurate representation of high-frequency water table fluctuations is of importance in the modeling of (1) the interaction between seawater and groundwater, more important, the effects on swash sediment transport and (2) the biological activities in the beach ecosystem. Capillarity effects provide a mechanism for high-frequency water table fluctuations. Previous modeling approaches adopted the assumption of saturated flow only and failed to predict the propagation of high-frequency fluctuations in the aquifer. In this paper we develop a modified kinematic boundary condition (kbc) for the water table which incorporates capillarity effects. The application of this kbc in a boundary element model enables the simulation of high-frequency water table fluctuations due to wave run-up. Numerical tests were carried out for a rectangular domain with small-amplitude oscillations; the behavior of water table responses was found to be similar to that predicted by an analytical solution based on the one-dimensional Boussinesq equation. The model was also applied to simulate the water table response to wave run-up on a doping beach. The results showed similar features of water table fluctuations observed in the field. In particular, these fluctuations are standing wave-like with the amplitude becoming increasingly damped inland. We conclude that the modified kbc presented here is a reasonable approximation of capillarity effects on beach water table fluctuations. However, further model validation is necessary before the model can confidently be used to simulate high-frequency water table fluctuations due to wave run-up.

Tide-induced water table fluctuations in coastal aquifers bounded by rhythmic shorelines

Previous studies on tidal water table dynamics in unconfined coastal aquifers have focused on the inland propagation of oceanic tides in the cross-shore direction based on the assumption of a straight coastline. Here, two-dimensional analytical solutions are derived to study the effects of rhythmic coastlines on tidal water table fluctuations. The computational results demonstrate that the alongshore variations of the coastline can affect the water table behavior significantly, especially in areas near the centers of the headland and embayment. With the coastline shape effects ignored, traditional analytical solutions may lead to large errors in predicting coastal water table fluctuations or in estimating the aquifer's properties based on these signals. The conditions under which the coastline shape needs to be considered are derived from the new analytical solution.

A new analytical solution for water table fluctuations in coastal aquifers with sloping beaches

The Boussinesq equation appears as the zeroth-order term in the shallow water flow expansion of the non-linear equation describing the flow of fluid in an unconfined aquifer. One-dimensional models based on the Boussinesq equation have been used to analyse tide-induced water table fluctuations in coastal aquifers. Previous analytical solutions for a sloping beach are based on the perturbation parameter, epsilon(N) = alphaepsilon cot beta (in which beta is the beach slope, alpha is the amplitude parameter and epsilon is the shallow water parameter) and are limited to tan(-1) (alphaepsilon) much less than beta less than or equal to pi/2. In this paper, a new higher-order solution to the non-linear boundary value problem is derived. The results demonstrate the significant influence of the higher-order components and beach slope on the water table fluctuations. The relative difference between the linear solution and the present solution increases as 6 and a increase, and reaches 7% of the linear solution. (C) 2003 Elsevier Ltd. All rights reserved.

Swash-aquifer interaction in the vicinity of the water table exit point on a sandy beach

The coupling of sandy beach aquifers with the swash zone in the vicinity of the water table exit point is investigated through simultaneous measurements of the instantaneous shoreline (swash front) location, pore pressures and the water table exit point. The field observations reveal new insights into swash-aquifer coupling not previously gleaned from measurements of pore pressure only. In particular, for the case where the exit point is seaward of the observation point, the pore pressure response is correlated with the distance between the exit point and the shoreline in that when the distance is large the rate of pressure drop is fast and when the distance is small the rate decreases. The observations expose limitations in a simple model describing exit point dynamics which is based only on the force balance on a particle of water at the sand surface and neglects subsurface pressures. A new modified form of the model is shown to significantly improve the model-data comparison through a parameterization of the effects of capillarity into the aquifer storage coefficient. The model enables sufficiently accurate predictions of the exit point to determine when the swash uprush propagates over a saturated or a partially saturated sand surface, potentially an important factor in the morphological evolution of the beach face. Observations of the shoreward propagation of the swash-induced pore pressure waves ahead of the runup limit shows that the magnitude of the pressure fluctuation decays exponentially and that there is a linear increase in time lags, behavior similar to that of tidally induced water table waves. The location of the exit point and the intermittency of wave runup events is also shown to be significant in terms of the shore-normal energy distribution. Seaward of the mean exit point location, peak energies are small because of the saturated sand surface within the seepage face acting as a "rigid lid'' and limiting pressure fluctuations. Landward of the mean exit point the peak energies grow before decreasing landward of the maximum shoreline position.

Water table waves in an unconfined aquifer: Experiments and modeling

[1] Comprehensive measurements are presented of the piezometric head in an unconfined aquifer during steady, simple harmonic oscillations driven by a hydrostatic clear water reservoir through a vertical interface. The results are analyzed and used to test existing hydrostatic and nonhydrostatic, small-amplitude theories along with capillary fringe effects. As expected, the amplitude of the water table wave decays exponentially. However, the decay rates and phase lags indicate the influence of both vertical flow and capillary effects. The capillary effects are reconciled with observations of water table oscillations in a sand column with the same sand. The effects of vertical flows and the corresponding nonhydrostatic pressure are reasonably well described by small-amplitude theory for water table waves in finite depth aquifers. That includes the oscillation amplitudes being greater at the bottom than at the top and the phase lead of the bottom compared with the top. The main problems with respect to interpreting the measurements through existing theory relate to the complicated boundary condition at the interface between the driving head reservoir and the aquifer. That is, the small-amplitude, finite depth expansion solution, which matches a hydrostatic boundary condition between the bottom and the mean driving head level, is unrealistic with respect to the pressure variation above this level. Hence it cannot describe the finer details of the multiple mode behavior close to the driving head boundary. The mean water table height initially increases with distance from the forcing boundary but then decreases again, and its asymptotic value is considerably smaller than that previously predicted for finite depth aquifers without capillary effects. Just as the mean water table over-height is smaller than predicted by capillarity-free shallow aquifer models, so is the amplitude of the second harmonic. In fact, there is no indication of extra second harmonics ( in addition to that contained in the driving head) being generated at the interface or in the interior.

Response of coastal groundwater table to offshore storms

Large groundwater table fluctuations were observed in a coastal aquifer during an offshore storm. The storm induced significant changes of the mean shoreline elevation, characterized by a pulse-like oscillation. This pulse propagated in the aquifer, resulting in the water table fluctuations. A general analytical solution is derived to quantify this new mechanism of water table fluctuation. The solution is applied to field observations and is found to be able to predict reasonably well the observed storm-induced water table fluctuations. Based on the analytical solution, the damping characteristics and phase shift of the oscillation as it propagates inland are examined.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

Reply to comment by A. G. J. Hilberts and P. A. Troch on "Influence of capillarity on a simple harmonic oscillating water table: Sand column experiments and modeling"

We thank Hilberts and Troch [2006] for their comment on our paper [Cartwright et al, 2005]. Before proceeding with our specific replies to the comments we would first like to clarify the definitions and meanings of equations (1)-(3) as presented by Hilberts and Troch [2006]. First, equation (1) is the fundamental definition of the (complex) effective porosity as derived by Nielsen and Perrochet [2000]. Equations (2) and (3), however, represent the linear frequency response function of the water table in the sand column responding to simple harmonic forcing. This function, which was validated by Nielsen and Perrochet [2000], provides an alternative method for estimating the complex effective porosity from the experimental sand column data in the absence of direct measurements of h_(tot) (which are required if equation (1) is to be used).

Groundwater waves in a coastal aquifer: A new governing equation including vertical effects and capillarity

Groundwater waves, that is, water table fluctuations, are a natural phenomenon in coastal aquifers. They represent an important part of the interaction between the ocean and aquifer and affect the mass exchange between them. This paper presents a new groundwater wave equation. Because it includes the effects of vertical flows and capillarity, the new equation is applicable to both intermediate-depth aquifers and high-frequency waves. Compared with the wave equation derived by Nielsen ed al. [1997], the present equation provides a closer representation of groundwater waves. In particular, it predicts high-frequency water table fluctuations as observed in the field. A validation of the new equation has been carried out by comparing the analytical solutions to it with predictions from direct simulations using the numerical model SUTRA. The effects of various physical parameters and their relative importance are also discussed.

BeachWin: modelling groundwater effects on swash sediment transport and beach profile changes

Field and laboratory observations have shown that a relatively low beach groundwater table enhances beach accretion. These observations have led to the beach dewatering technique (artificially lowering the beach water table) for combating beach erosion. Here we present a process-based numerical model that simulates the interacting wave motion on the beach. coastal groundwater flow, swash sediment transport and beach profile changes. Results of model simulations demonstrate that the model replicates accretionary effects of a low beach water table on beach profile changes and has the potential to become a tool for assessing the effectiveness of beach dewatering systems. (C) 2002 Elsevier Science Ltd. All rights reserved.

Analytical solution for tidal propagation in a coupled semi-confined/phreatic coastal aquifer

An analytical solution is derived for tidal fluctuations in a coupled coastal aquifer system consisting of a semi-confined aquifer, a thin semi-permeable layer and a phreatic aquifer. Based on the solution, we study the interactions (via leakage) between the confined and unconfined aquifers in response to tides. The results show that, under certain conditions, leakage from the confined aquifer can affect considerably the tidal water table fluctuation in the phreatic aquifer and vice versa. Ignoring these effects could lead to errors in estimating aquifer properties based on tidal signals. (C) 2002 Elsevier Science Ltd. All rights reserved.