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.

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.

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.

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.

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.

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.

Modelling the fate of chromated copper arsenate in a sandy soil

A pulse of chromated copper arsenate (CCA, a timber preservative) was applied in irrigation water to an undisturbed field soil in a laboratory column. Concentrations of various elements in the leachate from the column were measured during the experiment. Also, the remnants within the soil were measured at the end of the experiment. The geochemical modelling package, PHREEQC-2, was used to simulate the experimental data. Processes included in the CCA transport modelling were advection, dispersion, non-specific adsorption (cation exchange) and specific adsorption by clay minerals and organic matter, as well as other possible chemical reactions such as precipitation/dissolution. The modelling effort highlighted the possible complexities in CCA transport and reaction experiments. For example, the uneven dosing of CCA as well as incomplete knowledge of the soil properties resulted in simulations that gave only partial, although reasonable, agreement with the experimental data. Both the experimental data and simulations show that As and Cu are strongly adsorbed and therefore, will mostly remain at the top of the soil profile, with a small proportion appearing in leachate. On the other hand, Cr is more mobile and thus it is present in the soil column leachate. Further simulations show that both the quantity of CCA added to the soil and the pH of the irrigation water will influence CCA transport. Simulations suggest that application of larger doses of CCA to the soil will result in higher leachate concentrations, especially for Cu and As. Irrigation water with a lower pH will dramatically increase leaching of Cu. These results indicate that acidic rainfall or significant accidental spillage of CCA will increase the risk of groundwater pollution.

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.

Steepness expansion for free surface flows in coastal aquifers

Free surface flow of groundwater in aquifers has been studied since the early 1960s. Previous investigations have been based on the Boussinesq equation, derived from the non-linear kinematic boundary condition. In fact, the Boussinesq equation is the zeroth-order equation in the shallow-water expansion. A key assumption in this expansion is that the mean thickness of the aquifer is small compared with a reference length, normally taken to be the linear decay length. In this study, we re-examine the expansion scheme for free surface groundwater flows, and propose a new expansion wherein the shallow-water assumption is replaced by a steepness assumption. A comparison with experimental data shows that the new model provides a better prediction of water table levels than the conventional shallow-water expansion. The applicable ranges of the two expansions are exhibited. (c) 2004 Elsevier B.V. All rights reserved.

Application of a novel phytotoxicity assay for the detection of herbicides in Hervey Bay and the Great Sandy Straits

A novel phytotoxicity assay was incorporated into an environmental assessment of Hervey Bay and the Great Sandy Straits, to investigate the role of run-off associated herbicides in the deteriorated health of intertidal seagrass meadows. Dose response curves of common herbicides were performed and their toxicity equivalents elucidated to assist in analysis. The results of the assay were reproducible and corresponded strongly with results of chemical analyses. The incorporation of the assay into the assessment of surface waters added an important aspect to the study by allowing investigation of the toxicity of cumulative herbicide concentrations and yielding biologically relevant data. The highest herbicide concentration detected during the study was equivalent to 0.23 mu g 1(-1) diuron; a concentration known to inhibit photosynthetic efficiency of the assay biomaterial by approximately 3%. (c) 2004 Elsevier Ltd. All rights reserved.

Tidal impacts on riparian salinities near estuaries

Groundwater-dependent riparian biota is known to be sensitive to changes in soil and groundwater salinity in estuarine systems. The groundwater flow and salinity behaviour in a phreatic aquifer adjoining a partially penetrating, tidal. estuary is investigated through two-dimensional numerical experiments for a lateral cross-section, which explore the influence of factors, such as aquifer and soil materials, tidal amptitudes, and regional groundwater hydraulic gradients. The density contrast between estuarine water and the fresh groundwater drives saltwater penetration of the aquifer even in the case of a marked groundwater hydraulic gradient towards the estuary. We show that tidal fluctuations in estuaries can significantly affect the groundwater salinity distribution in adjacent density-stratified phreatic aquifers. This has consequences for the expected distribution of salinity-sensitive biota in the hyporheic zone as well as vegetation and fauna dependent on water in the riparian soil and aquifer. The shape of the dense saltwater wedge propagating into the adjacent groundwater system is also modified by the estuarine tidal signal, although this effect appears to have only minor influence on the maximum distance penetrated into the aquifer (i.e., location of the 'toe' of the wedge). Tide-induced changes to riparian groundwater salinity are advection-driven, as evidenced by the modified time-averaged groundwater flow dynamics. (c) 2006 Elsevier B.V. All rights reserved.