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.

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.

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.

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.

A composite DOP approach to excluding bottom reflectance in mapping water parameters of shallow coastal zones from TM imagery

In this paper we proposed a composite depth of penetration (DOP) approach to excluding bottom reflectance in mapping water quality parameters from Landsat thematic mapper (TM) data in the shallow coastal zone of Moreton Bay, Queensland, Australia. Three DOPs were calculated from TM1, TM2 and TM3, in conjunction with bathymetric data, at an accuracy ranging from +/-5% to +/-23%. These depths were used to segment the image into four DOP zones. Sixteen in situ water samples were collected concurrently with the recording of the satellite image. These samples were used to establish regression models for total suspended sediment (TSS) concentration and Secchi depth with respect to a particular DOP zone. Containing identical bands and their transformations for both parameters, the models are linear for TSS concentration, logarithmic for Secchi depth. Based on these models, TSS concentration and Secchi depth were mapped from the satellite image in respective DOP zones. Their mapped patterns are consistent with the in situ observed ones. Spatially, overestimation and underestimation of the parameters are restricted to localised areas but related to the absolute value of the parameters. The mapping was accomplished more accurately using multiple DOP zones than using a single zone in shallower areas. The composite DOP approach enables the mapping to be extended to areas as shallow as <3 m. (C) 2004 Elsevier Inc. All rights reserved.

Sugarcane water use from shallow water tables: implications for improving irrigation water use efficiency

The resource potential of shallow water tables for cropping systems has been investigated using the Australian sugar industry as a case study. Literature concerning shallow water table contributions to sugarcane crops has been summarised, and an assessment of required irrigation for water tables to depths of 2 m investigated using the SWIMv2.1 soil water balance model for three different soils. The study was undertaken because water availability is a major limitation for sugarcane and other crop production systems in Australia and knowledge on how best to incorporate upflow from water tables in irrigation scheduling is limited. Our results showed that for the three soils studied (representing a range of permeabilities as defined by near-saturated hydraulic conductivities), no irrigation would be required for static water tables within 1 m of the soil surface. Irrigation requirements when static water tables exceeded 1 m depth were dependent on the soil type and rooting characteristics (root depth and density). Our results also show that the near-saturated hydraulic conductivities are a better indicator of the ability of water tables below 1 m to supply sufficient upflow as opposed to soil textural classifications. We conclude that there is potential for reductions in irrigation and hence improvements in irrigation water use efficiency in areas where shallow water tables are a low salinity risk: either fresh, or the local hydrology results in net recharge. (C) 2003 Elsevier B.V. All rights reserved.

New approximation for free surface flow of groundwater: capillarity correction

An existing capillarity correction for free surface groundwater flow as modelled by the Boussinesq equation is re-investigated. Existing solutions, based on the shallow flow expansion, have considered only the zeroth-order approximation. Here, a second-order capillarity correction to tide-induced watertable fluctuations in a coastal aquifer adjacent to a sloping beach is derived. A new definition of the capillarity correction is proposed for small capillary fringes, and a simplified solution is derived. Comparisons of the two models show that the simplified model can be used in most cases. The significant effects of higher-order capillarity corrections on tidal fluctuations in a sloping beach are also demonstrated. (c) 2004 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.

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.

A trace element study of siderite-jasper banded iron formation in the 3.45 Ga Warrawoona Group, Pilbara Craton - Formation from hydrothermal fluids and shallow seawater

Shale-normalised rare earth element and yttrium (REE + Y) patterns for siderite-jasper couples in a banded iron formation of the 3.45 Ga Panorama Formation, Warrawoona Group, eastern Pilbara Craton, display distinct positive Y and Eu anomalies and weak positive La and Gd anomalies, combined with depleted light REE relative to middle and heavy REE. Ambient seawater and hydrothermal fluids are identified as major sources of REE + Y for the BIF. In the case of siderites, strong correlations between incompatible trace elements and trace element ratios diagnostic of seawater indicate variable input from a terrigenous source (e.g. volcanic ash). We propose a volcanic caldera setting as a likely depositional environment where jasper and siderite precipitated as alternating bands in response to episodic changes in ambient water chemistry. The episodicity was either driven by fluctuations in the intensity of hydrothermal activity or changes in magma chamber activity, which in turn controlled relative sea level. In this context, precipitation of jasper probably reflects background conditions during which seawater was saturated in silica due to evaporative conditions, while siderites were deposited most likely during intermittent periods of enhanced volcanic activity when seawater was more acidic due to the release of exhalative phases (e.g. CO2). © 2005 Elsevier B.V. All rights reserved.

A modern soft-bottom, shallow-water crinoid fauna (Echinodermata) from the Great Barrier Reef, Australia

A recent preliminary survey revealed that 12 species of unstalked crinoids occur on a gentle sandy slope (12-18 m depth) at Lizard Island, Great Barrier Reef, Australia; five of which are also found on coral reefs. The other seven appear to constitute a unique assemblage restricted to unconsolidated substrates, where most cling to algae or hide beneath rubble or sponges. Members of this assemblage exhibit all of the basic feeding postures found among reef-dwelling species. However, Comatula rotalaria, which lacks anchoring cirri and bears uniquely differentiated short and long arms, exhibits a posture different from other living crinoids. Quantitative transects reveal apparent depth-related differences in species composition: C. rotalaria dominated the 12 transects in 12-13 m (84% of 82 specimens), while Comatella nigra, Comatula cf. purpurea, Amphimetra cf. tessellata and Zygometra microdiscus accounted for 96% of 54 specimens observed along 12 transects in 16-17 m.

Quantifying the effects of bycatch reduction devices in Queensland's (Australia) shallow water eastern king prawn (Penaeus plebejus) trawl fishery

This study presents results from an experimental 10-day research charter that was designed to quantify the effects of (a) a turtle excluder device (TED), (b) a radial escape section bycatch reduction device (BRD) and (c) both devices together, on bycatch and prawn catch rates in the Queensland shallow water eastern king prawn (Penaeus plebejus) trawl fishery. The bycatch was comprised of 250 taxa, mainly gurnards, whiting, lizard fish, flathead, dragonets, portunid crabs, turretfish and flounders. The observed mean catch rates of bycatch and marketable eastern king prawns from the standard trawl net (i.e., net with no TED or BRD) used during the charter were 11.06 kg/hectare (ha(-1)) (S.E. 0.90) swept by the trawl gear and 0.94 kg ha(-1), respectively. For the range of depths sampled (20.1-90.7 m), bycatch rates declined significantly at a rate of 0.14 kg ha-1 for every 1 m increase in depth, while prawn catch rates were unaffected. When both the TED and radial escape section BRD were used together, the bycatch rate declined by 24% compared to a standard net, but at a 20% reduction in marketable prawn catch rate. The largest reductions were achieved for stout whiting Sillago robusta (57% reduction) and yellowtail scad Trachurus novaezelandiae (32% reduction). Multidimensional scaling and analysis of similarities revealed that bycatch assemblages differed significantly between depths and latitude, but not between the different combinations of bycatch reduction devices. Despite the lowered prawn catch rates, the reduced bycatch rates are promising, particularly for S. robusta, which is targeted in another fishery. Prawn trawl operators are not permitted to retain S. robusta and the devices examined herein offer the potential to significantly reduce the incidental fishing mortality that this species experiences. (c) 2006 Elsevier B.V. All rights reserved.

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.