Drying front in a sloping aquifer: Nonlinear effects

[1] The profiles for the water table height h(x, t) in a shallow sloping aquifer are reexamined with a solution of the nonlinear Boussinesq equation. We demonstrate that the previous anomaly first reported by Brutsaert [1994] that the point at which the water table h first becomes zero at x = L at time t = t(c) remains fixed at this point for all times t > t(c) is actually a result of the linearization of the Boussinesq equation and not, as previously suggested [Brutsaert, 1994; Verhoest and Troch, 2000], a result of the Dupuit assumption. Rather, by examination of the nonlinear Boussinesq equation the drying front, i.e., the point x(f) at which h is zero for times t greater than or equal to t(c), actually recedes downslope as physically expected. This points out that the linear Boussinesq equation should be used carefully when a zero depth is obtained as the concept of an average'' depth loses meaning at that time.

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.

Experimental observations of watertable waves in an unconfined aquifer with a sloping boundary

Observations of horizontal and vertical variations in piezometric head in a homogeneous, laboratory aquifer are presented and discussed. The observed fluctuations are induced by a simple harmonic oscillation in the clear water reservoir acting across a sloping boundary. The data qualitatively supports existing theories in that higher harmonics are generated in the active forcing zone and that a significant increase in the inland, asymptotic watertable over height (relative to that found for the vertical boundary case) is observed. The observed overheight is shown to be accurately reproduced by existing small-amplitude perturbation theory. Detailed measurements in the vicinity of the sloping boundary reveal that the signal of generated higher harmonics is strongest near the sand surface and that vertical flows are significant in this region. The aquifer is of finite-depth and is influenced by capillary effects, the experimental data therefore exposes limitations of theories which are based on the assumption of a shallow aquifer free of capillary effects. The dispersive properties of the measured pressure wave in the aquifer are comparable to those found from field observations and likewise do not agree with those predicted by the capillary free, shallow aquifer theory. Although some improvement is obtained, discrepancies between the data and theory persist even when a finite-depth aquifer and capillary effects are considered in the theoretical model. Further sand column experiments eliminate a truncated capillary fringe as a possible contributor to these discrepancies. However, the neglect of horizontal flows in the fringe may have caused the discrepancies. (C) 2004 Elsevier Ltd. All rights reserved.

Tidal influence on behaviour of a coastal aquifer adjacent to a low-relief estuary

The tidal influence on groundwater hydrodynamics, salt-water intrusion and submarine groundwater discharge from coastal/estuarine aquifers is poorly quantified for systems with a mildly sloping beach, in contrast to the case where a vertical beach face is assumed. We investigated the effect of beach slope for a coastal aquifer adjacent to a low-relief estuary, where industrial waste was emplaced over the aquifer. The waste was suspected to discharge leachate towards the estuary. Field observations at various locations showed that tidally induced groundwater head fluctuations were skewed temporally. Frequency analysis suggested that the fluctuation amplitudes decreased exponentially and the phase-tags increased Linearly for the primary tidal signals as they propagated inland. Salinisation zones were observed in the bottom part of the estuary and near the beach surface. Flow and transport processes in a cross-section perpendicular to the estuary were simulated using SEAWAT-2000, which is capable of depicting density-dependent flow and multi-species transport. The simulations showed that the modelled water table fluctuations were in good agreement with the monitored data. Further simulations were conducted to gain insight into the effects of beach slope. In particular the limiting case of a vertical beach face was considered. The simulations showed that density difference and tidal forcing drive a more complex hydrodynamic pattern for the mildly sloping beach than the vertical beach, as well as a profound asymmetry in tidally induced water table fluctuations and enhanced salt-water intrusion. The simulation results also indicated that contaminant transport from the aquifer to the estuary was affected by the tide, where for the mildly sloping beach, the tide tended to intensify the vertical mass exchange in the vicinity of the shorelines, (c) 2005 Elsevier B.V. All rights reserved.

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.

Two-dimensional approximation for tide-induced watertable fluctuations in a sloping sandy beach

The prediction of watertable fluctuations in a coastal aquifer is important for coastal management. However, most previous approaches have based on the one-dimensional Boussinesq equation, neglecting variations in the coastline and beach slope. In this paper, a closed-form analytical solution for a two-dimensional unconfined coastal aquifer bounded by a rhythmic coastline is derived. In the new model, the effect of beach slope is also included, a feature that has not been considered in previous two-dimensional approximations. Three small parameters, the shallow water parameter (epsilon), the amplitude parameter (a) and coastline parameter (beta) are used in the perturbation approximation. The numerical results demonstrate the significant influence of both the coastline shape and beach slopes on tide-driven coastal groundwater fluctuations. (c) 2004 Elsevier Ltd. All rights reserved.

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.

Assessment of groundwater recharge and water fluxes of the Guarani Aquifer System, Brazil

The groundwater recharge and water fluxes of the Guarani Aquifer System in the state of Sao Paulo in Brazil were assessed through a numeric model. The study area (6,748 km(2)) comprises Jacar,-Gua double dagger A(0) and Jacar,-Pepira River watersheds, tributaries of the Tiet River in the central region of the state. GIS based tools were used in the storage, processing and analysis of data. Main hydrologic phenomena were selected, leading to a groundwater conceptual model, taking into account the significant outcrops occurring in the study area. Six recharge zones were related to the geologic formation and structures of the semi-confined and phreatic aquifer. The model was calibrated against the baseflows and static water levels of the wells. The results emphasize the strong interaction of groundwater flows between watersheds and the groundwater inflow into the rivers. It has been concluded that lateral groundwater exchanges between basins, the deep discharges to the regional system, and well exploitation were not significant aquifer outflows when compared to the aquifer recharge. The results have shown that the inflows from the river into the aquifer are significant and have the utmost importance since the aquifer is potentially more vulnerable in these places.

Tidal effects on sea water intrusion in unconfined aquifers

A variable-density groundwater model is used to analyse the effects of tidal fluctuations on sea-water intrusion in an unconfined aquifer. It is shown that the tidal activity forces the sea-water to intrude further inland and it also creates a thicker interface than would occur without tidal effects. Moreover, the configuration of the interface is radically changed when the tidal fluctuations are included. This is because of the dramatic changes in the flow pattern and velocity of the groundwater near the shoreline. For aquifer depths much larger than tidal amplitudes, the tidal fluctuation does not have much effect on how far the sea-water intrudes into the aquifer; nevertheless, a significant change in the configuration of concentration contours because of the effect of tidal fluctuations is observed. This change is more noticeable at the top of the aquifer, near the water table, than at the bottom of the aquifer, and is caused by the infiltration of salt water into the top of the aquifer at higher tidal levels. A flatter beach slope, therefore, intensifies this phenomenon. The interface configurations do not change noticeably over the course of a tidal cycle. Neglecting tidal fluctuation effects results in an inaccurate evaluation of the water table elevation at the land end of the aquifer, although no distinguishable difference is seen between the water tables near the shoreline. Where the landward boundary condition is a constant head, the effects of tidal fluctuations on sea-water intrusion are more pronounced than for cases where the landward boundary condition is a specified flux. Also it is shown that the effects of tidal fluctuations are more significant for a sloping beach than for a vertical shoreline and the salt water intrudes further inland for the sloping case. (C) 1999 Elsevier Science B.V. All rights reserved.

Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

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.

Tidal fluctuations in a leaky confined aquifer: Dynamic effects of an overlying phreatic aquifer

Tidal fluctuations in a leaky confined coastal aquifer are damped significantly due to leakage into an overlying phreatic aquifer. Jiao and Tang [1999] presented an analytical solution to a simple model describing this phenomenon. Their solution assumes that the tidal fluctuations in the overlying phreatic aquifer are negligible (i.e,, a static phreatic aquifer), Here we examine dynamic effects of the overlying aquifer based on a new approximate analytical solution. The numerical results indicate that the dynamic effects can be significant for a relatively large leakage and a high transmissivity of the phreatic aquifer.

Sudden drawdown and drainage of a horizontal aquifer

Drainage of a saturated horizontal aquifer following a sudden drawdown is reanalyzed using the Boussinesq equation. The effect of the finite length of the aquifer is considered in detail. An analytical approximation based on a superposition principle yields a very good estimate of the outflow when compared to accurate numerical solutions. An illustration of the new analytical approach to analyze basin-scale field data is used to demonstrate possible field applications of the new solution.