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.

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.

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.

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.

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.

Measurement and modelling of controlled beach groundwater levels under wave action

New laboratory scale experimental data are presented on the forcing of beach groundwater levels by wave run-up. The experimental setup simulates a coastal barrier dividing the ocean from a relatively constant back beach water level, conditions approximating a closed off lagoon system or beach aquifer. The data are critically compared to an advanced numerical model for simulating wave and beach groundwater interaction in the coastal zone, and provide the first experimental verification of such a model. Overall model-data comparisons are good, but some systematic discrepancies are apparent, and reasons for these are discussed.

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.

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.

Collapse behavior of soil in a Brazilian region affected by a rising water table

Collapsible soils are usually nonsaturated, low density, and metastable-structured soils that are known to exhibit a volume reduction following an episode of moisture increase or suction reduction. This paper describes the collapsible behavior of clayey sand based on controlled soil suction tests carried out on undisturbed samples from the city of Pereira Barreto, in the State of Sao Paulo, Brazil. Foundation settlements due to soil collapse are common in this region and occurred during the filling of the reservoir of the Tres Irmaos Dam, which induced the elevation of the groundwater table in different parts of Pereira Barreto. This paper shows that collapse strains depend on the stress and soil suction acting in the sample and that saturation is not necessary for a collapse to occur. The influence of soil suction, gradual wetting, and the wetting and drying cycle on the collapsible behavior of the soil is also shown and discussed.

Groundwater fluctuations and footslope ferricrete soils in the humid tropical zone of southern Cameroon

This paper discusses the relationship between the differentiation of ferruginous accumulations and the variable water saturation of footslope soil patterns. An analysis of the slope morphology of a typical hill in the forest zone of southern Cameroon and a seasonal survey of the levels of groundwaters, springs and rivers were considered in relation to the petrology of different soil patterns. The study site is a tabular hillock whose slopes present a progressive development from steep to gentle slopes. The variable residence time of water within the soil, creating an alternation of reducing and oxidizing conditions, affects oil chemistry, structure and lateral extension of the soil patterns. The ferruginous soil patterns, being formed on the footslopes, gradually increase in extent with decreasing slope angle and the relative rise of the groundwater level. The steep footslopes, where groundwater has a shorter residence time, show a soft mottled clay pattern, restricted to the bottom part of the slope. The moderate footslopes exhibit a deep permanent and a temporary perched groundwater table. The latter, with its regular capillary fringe, contributes to more reducing conditions within isolated domains in the soil patterns, and thus to the alternation with oxidizing conditions, generating a continuous hard soil pattern (massive carapace). The more gently dipping footslopes exhibit groundwater levels near the surface and also a significant amplitude of groundwater fluctuation. Iron, previously accumulated in moderate footslope patterns, is reduced, remobilized, and leached. The soil patterns formed develop into a variegated carapace, more extended along the slope, containing less iron, but nevertheless more hardened, due to the important fluctuations of the groundwater table. These patterns are limited to the zone of groundwater fluctuation and deteriorate as the water fluctuation zone recedes. Copyright (c) 2005 John Wiley & Sons, Ltd.

A mathematical view of water table fluctuations in a shallow aquifer in Brazil

Detailed monitoring of the groundwater table can provide important data about both short- and long-term aquifer processes, including information useful for estimating recharge and facilitating groundwater modeling and remediation efforts. In this paper, we presents results of 4 years (2002 to 2005) of monitoring groundwater water levels in the Rio Claro Aquifer using observation wells drilled at the Rio Claro campus of São Paulo State University in Brazil. The data were used to follow natural periodic fluctuations in the water table, specifically those resulting from earth tides and seasonal recharge cycles. Statistical analyses included methods of time-series analysis using Fourier analysis, cross-correlation, and R/S analysis. Relationships could be established between rainfall and well recovery, as well as the persistence and degree of autocorrelation of the water table variations. We further used numerical solutions of the Richards equation to obtain estimates of the recharge rate and seasonable groundwater fluctuations. Seasonable soil moisture transit times through the vadose zone obtained with the numerical solution were very close to those obtained with the cross-correlation analysis. We also employed a little-used deep drainage boundary condition to obtain estimates of seasonable water table fluctuations, which were found to be consistent with observed transient groundwater levels during the period of study.

The Extended Water Footprint and illegal groundwater use in the Upper Guadiana Basin (Spain)

Tablas de Daimiel National Park is located in the Upper Guadiana Basin and represents one of the largest and most important wetlands in Europe. The long term ecological integrity of this wetland is inherently associated with the maintenance of a shallow groundwater table, namely the Western Mancha aquifer (WMA) or Aquifer 23. The intensive use of groundwater, mainly for irrigation, has led over the last decades to deep socio‐economic changes. Such intensive use has also lowered the water table of Aquifer 23, drastically reducing the flooded area of the wetland and threatening its ecological integrity. A number of plans and measures have been developed and implemented since the declaration of overexploitation of Aquifer 23 in the year 1987. The most recent one is the Special Plan for the Upper Guadiana (SPUG), approved in 2008. This Plan is the main measure to comply with achieving the objective of good quantitative and qualitative status required under the Water Framework Directive (2000). This paper offers a new type of integrated analysis which allows assessing under a common lens the physical, economic and social dimensions of groundwater use in the area. The first objective is to calculate the groundwater footprint of agricultural production in the Upper Guadiana basin and its evolution during 2000‐2008. For this purpose, we have applied the Extended Water Footprint (EWF) methodology ‐a novel approach based on the classical Water Footprint (WF) approach‐ that includes an assessment of the water productivity from an economic and social perspective. Compared to the classical WF, the EWF allows for a more complete overview of the sector, providing new insights for policy decisions (e.g. to define options and possibilities on water re‐allocation in order to achieve both better ecosystem conservation and social equity). The second objective is to use the EWF to compare the existing authorized and non‐authorized or illegal use of water. This allows us to discuss current initiatives by public authorities in relation to the existing frame of water rights

Chemical changes in argisols under irrigated grape production in the central São Francisco river valley, Brazil

This study compares the chemical composition of the solution and exchange complex of soil in a 3-year-old irrigated vineyard (Vitis vinifera L., Red Globe cultivar) with that of adjacent clearing in the native hyperxerophyllic 'caatinga' vegetation. The soils are classified as Plinthic Eutrophic Red-Yellow Argisol; according to Soil Taxonomy they are isohyperthermic Plinthustalfs. Detailed physiographic characterization revealed an impermeable gravel and cobble covering the crystalline rocks; the relief of this layer was more undulating than the level surface. Significant higher concentrations of extractable Na, K, Mg and Ca were observed within the vineyard. Lower soil acidity, higher Ca/Mg ratios, as well as lower sodium adsorption and Na/K ratios reflected additions of dolomitic lime, superphosphate and K-bearing fertilizers. As the water of the São Francisco River is of good quality for irrigation (C1S1), the increases in Na were primarily attributed to capillary rise from the saline groundwater table. None of the soil in the study area was found to be sodic. About 62% of the vineyard had an Ap horizon with salinity levels above 1.5dSm-1 (considered detrimental for grape production); according to average values for this horizon, a potential 13% reduction in grape production was predicted. Differences in chemical composition in function of distance to the collector canals were observed in the clearing, but not in the vineyard. The influence of differences in the elevations of the surface and impermeable layers, as well as pediment thickness, was generally weaker under irrigation. Under irrigation, soil moisture was greater in points of convergent surface waterflow; the effect of surface curvature on chemical properties, though less consistent, was also stronger in the vineyard.