Regional-scale, fully coupled modelling of stream-aquifer interaction in a tropical catchment

The planning and management of water resources in the Pioneer Valley, north-eastern Australia requires a tool for assessing the impact of groundwater and stream abstractions on water supply reliabilities and environmental flows in Sandy Creek (the main surface water system studied). Consequently, a fully coupled stream-aquifer model has been constructed using the code MODHMS, calibrated to near-stream observations of watertable behaviour and multiple components of gauged stream flow. This model has been tested using other methods of estimation, including stream depletion analysis and radon isotope tracer sampling. The coarseness of spatial discretisation, which is required for practical reasons of computational efficiency, limits the model's capacity to simulate small-scale processes (e.g., near-stream groundwater pumping, bank storage effects), and alternative approaches are required to complement the model's range of applicability. Model predictions of groundwater influx to Sandy Creek are compared with baseflow estimates from three different hydrograph separation techniques, which were found to be unable to reflect the dynamics of Sandy Creek stream-aquifer interactions. The model was also used to infer changes in the water balance of the system caused by historical land use change. This led to constraints on the recharge distribution which can be implemented to improve model calibration performance. (c) 2006 Elsevier 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.

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.

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

Damping of tidal head fluctuations in a leaky confined coastal aquifer is enhanced by leakage into an overlying phreatic aquifer. We show that the phreatic aquifer is, however, resistant to the leakage flow and in particular, a deep phreatic aquifer can reduce the leakage effects significantly. An analytical solution, based on a vertical flow model for the phreatic aquifer, is derived for quantifying the role of this upper free water body in tidal propagation in the lower semi-confined aquifer. (C) 2002 Elsevier Science B.V. 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.

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.

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.

Confined-unconfined flow in a horizontal confined aquifer

Mixed confined and unconfined groundwater flow occurs in a bounded initially dry aquifer when the hydraulic head at the side boundary suddenly rises above the elevation of the aquifer's top boundary. The flow problem as modelled by the Boussinesq equation is non-trivial because of the involvement of two moving boundaries. The transformed equation (based on a similarity transformation) can, however, be dealt with more easily. Here, we present an approximate analytical solution for this flow problem. The approximate solution is compared with an 'exact' numerical solution and found to be a very accurate description for describing the mixed confined and unconfined flow in the confined aquifer. (C) 2002 Elsevier Science B.V. All rights reserved.

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 the salt-freshwater interface in a coastal aquifer to a wave-induced groundwater pulse: field observations and modelling

Field observations on an unconfined coastal aquifer showed that a groundwater pulse, generated by it moderate (significant wave height, H-sig similar to 4.5 m) wave/storm event, induced significant oscillations in the salt-freshwater interface of the order of several metres in the horizontal direction. A dynamic sharp-interface model is developed to quantify the mechanism of these interface oscillations. The model uses the 50% seawater salinity contour as the location of the equivalent sharp-interface. The model was calibrated against the observed groundwater table fluctuations. It predicted reasonably well the interface oscillations with a slight over-prediction of the oscillation magnitude and a steepening of the interface. The neglect of mixing in the salt-freshwater mixing zone by the sharp-interface model is suggested as a possible contributor to the discrepancies between the model predictions and observations. In contrast with the significant wave effects, there was no observable response of the interface to diurnal or semidiurnal tides. (C) 2004 Elsevier Ltd. All rights reserved.