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.

Critical review of acid in situ leach uranium mining: 1. USA and Australia

The technique of in situ leach (ISL) uranium mining is well established in the USA, as well as being used extensively in Eastern Europe and the former Soviet Union. The method is being proposed and tested on uranium deposits in Australia, with sulfuric acid chemistry and no restoration of groundwater following mining. Test sites in the USA were required to restore groundwater to ascertain the extent of impacts and compare costs to alkaline ISL mines. The problems encountered include expensive and difficult restoration, gypsum precipitation, higher salinity and some heavy metals and radionuclides after restoration. One of the most critical issues is whether natural attenuation is capable of restoring groundwater quality and geochemical conditions in an acid leached aquifer zone. The history of acid ISL sites in the USA and Australia are presented in this study, with a particular focus on the demonstration of restoration of groundwater impacts.

Tidal effects on groundwater dynamics in unconfined aquifers

The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi-steady-slate rise in the mean water-table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed anti this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water-table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright (C) 2001 John Wiley & Sons, Ltd.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

Petrology and geochemistry of the Bandas del Sur Formation, Las Canadas Edifice, Tenerife (Canary Islands)

The Bandas del Sur Formation preserves a Quaternary extra-caldera record of central phonolitic explosive volcanism of the Las Canadas volcano at Tenerife. Volcanic rocks are bimodal in composition, being predominantly phonolitic pyroclastic deposits, several eruptions of which resulted in summit caldera collapse, alkali basaltic lavas erupted from many fissures around the flanks. For the pyroclastic deposits, there is a broad range of pumice glass compositions from phonotephrite to phonolite. The phonolite pyroclastic deposits are also characterized by a diverse, 7-8-phase phenocryst assemblage (alkali feldspar + biotite + sodian diopside + titanomagnetite + ilmenite + nosean-hauyne + titanite + apatite) with alkali feldspar dominant, in contrast to interbedded phonolite lavas that typically have lower phenocryst contents and lack hydrous phases. Petrological and geochemical data are consistent with fractional crystallization (involving the observed phenocryst assemblages) as the dominant process in the development of phonolite magmas. New stratigraphically constrained data indicate that petrological and geochemical differences exist between pyroclastic deposits of the last two explosive cycles of phonolitic volcanism. Cycle 2 (0.85-0.57 Ma) pyroclastic fall deposits commonly show a cryptic compositional zonation indicating that several eruptions tapped chemically, and probably thermally stratified magma systems. Evidence for magma mixing is most widespread in the pyroclastic deposits of Cycle 3 (0.37-0.17 Ma), which includes the presence of reversely and normally zoned phenocrysts, quenched mafic glass blebs in pumice, banded pumice, and bimodal to polymodal phenocryst compositional populations. Syn-eruptive mixing events involved mostly phonolite and tephriphonolite magmas, whereas a pre-eruptive mixing event involving basaltic magma is recorded in several banded pumice-bearing ignimbrites of Cycle 3. The periodic addition and mixing of basaltic magma ultimately may have triggered several eruptions. Recharge and underplating by basaltic magma is interpreted to have elevated sulphur contents (occurring as an exsolved gas phase) in the capping phonolitic magma reservoir. This promoted nosean-hauyne crystallization over nepheline, elevated SO3 contents in apatite, and possibly resulted in large, climatologically important SO2 emissions.

Experimental investigation of contaminant transport in coastal groundwater

Contaminant transport in coastal aquifers is of increasing interest since, with the development of coastal areas, contaminants from surface sources may enter coastal aquifers and pollute the groundwater flow. Coastal groundwater flow is complicated because of the presence of a freshwater-saltwater diffusion zone and the tidal variation of sea level at the seaward end. This paper investigates experimentally the behaviour of contaminant plumes with different densities in an unconfined coastal aquifer. Experiments were performed in a flow tank filled with glass beads as the porous medium. Results show that the dense contaminant has a more diffusive front than the less dense one in the seaward direction towards the coastline. The plume becomes more diffusive when it travels closer to the saltwater interface. On the contrary, the less dense contaminant presents a relatively sharp outline. It tends to migrate in the upper portion of the aquifer and exits in a concentrated manner over a small discharge area at the coastline, not further seaward under the sea. Non-dimensional parameters show that instabilities occur in our experiments for a density difference of 1.2% or larger between the contaminant and the ambient water. The experimental results provide guidance for field monitoring and numerical modelling. (C) 2002 Elsevier Science Ltd. All rights reserved.

Numerical modelling of contaminant transport in coastal aquifers

This paper employs a two-dimensional variable density flow and transport model to investigate the transport of a dense contaminant plume in an unconfined coastal aquifer. Experimental results are also presented to show the contaminant plume in a freshwater-seawater flow system. Both the numerical and experimental results suggest that the neglect of the seawater interface does not noticeably affect the horizontal migration rate of the plume before it reaches the interface. However, the contaminant will travel further seaward and part of the solute mass will exit under the sea if the higher seawater density is not included. If the seawater density is included, the contaminant will travel upwards towards the beach along the freshwater-saltwater interface as shown experimentally. Neglect of seawater density, therefore, will result in an underestimate of solute mass rate exiting around the coastline. (C) 2002 IMACS. Published by Elsevier Science B.V. All rights reserved.

Tide-induced water table fluctuations in coastal aquifers bounded by rhythmic shorelines

Previous studies on tidal water table dynamics in unconfined coastal aquifers have focused on the inland propagation of oceanic tides in the cross-shore direction based on the assumption of a straight coastline. Here, two-dimensional analytical solutions are derived to study the effects of rhythmic coastlines on tidal water table fluctuations. The computational results demonstrate that the alongshore variations of the coastline can affect the water table behavior significantly, especially in areas near the centers of the headland and embayment. With the coastline shape effects ignored, traditional analytical solutions may lead to large errors in predicting coastal water table fluctuations or in estimating the aquifer's properties based on these signals. The conditions under which the coastline shape needs to be considered are derived from the new analytical solution.

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.

Nitrogen ecophysiology of Heron Island, a subtropical coral cay of the Great Barrier Reef, Australia

Coral cays form part of the Australian Great Barrier Reef. Coral cays with high densities of seabirds are areas of extreme nitrogen (N) enrichment with deposition rates of up to 1000 kg N ha(-1) y(-1). The ways in which N sources are utilised by coral cay plants, N is distributed within the cay, and whether or not seabird-derived N moves from cay to surrounding marine environments were investigated. We used N metabolite analysis, N-15 labelling and N-15 natural abundance (delta(15)N) techniques. Deposited guano-derived uric acid is hydrolysed to ammonium (NH4+) and gaseous ammonia (NH3). Ammonium undergoes nitrification, and nitrate (NO3-) and NH4+ were the main forms of soluble N in the soil. Plants from seabird rookeries have a high capacity to take up and assimilate NH4+, are able to metabolise uric acid, but have low rates of NO3- uptake and assimilation. We concluded that NH4+ is the principal source of N for plants growing at seabird rookeries, and that the presence of NH4+ in soil and gaseous NH3 in the atmosphere inhibits assimilation of NO3-, although NO3- is taken up and stored. Seabird guano, Pisonia forest soil and vegetation were similarly enriched in N-15 suggesting that the isotopic enrichment of guano (delta(15)N 9.9parts per thousand) carries through the forest ecosystem. Soil and plants from woodland and beach environments had lower delta(15)N (average 6.5parts per thousand) indicating a lower contribution of bird-derived N to the N nutrition of plants at these sites. The aquifer under the cay receives seabird-derived N leached from the cay and has high concentrations of N-15-enriched NO3- (delta(15)N 7.9parts per thousand). Macroalgae from reefs with and without seabirds had similar delta(15)N values of 2.0-3.9parts per thousand suggesting that reef macroalgae do not utilise N-15-enriched seabird-derived N as a main source of N. At a site beyond the Heron Reef Crest, macroalgae had elevated delta(15)N of 5.2parts per thousand, possibly indicating that there are locations where macroalgae access isotopically enriched aquifer-derived N. Nitrogen relations of Heron Island vegetation are compared with other reef islands and a conceptual model is presented.

Palaeoclimatic implications of the growth history and stable isotope (delta O-18 and delta C-13) geochemistry of a Middle to Late Pleistocene stalagmite from central-western Italy

The age structure and, stable isotope composition of a stalagmite (CC I) from an upland cave in central-western Italy were studied to investigate regional response to global climatic changes. Four growth phases are constrained by 28 thermal ionization and multi-collector inductively coupled plasma mass spectrometry Th-U ages and reveal intermittent deposition through the period between Marine isotope Stage (MIS) 11 and 3 (similar to380 and similar to43 kyr). Most of the growth took place between similar to380 and similar to280 kyr, a period punctuated briefly by a hiatus in deposition through the glacial maximum of MIS 10. Growth was terminated abruptly at 280 kyr just prior to the MIS 8 glacial maximum. With a present-day chamber temperature of 7.5 degreesC, the timing of hiatuses close to these glacial maxima point to freezing conditions at the time. No deposition was recorded through the entirety of MIS 7 and most of MIS 6, whilst two minor growth phases occurred at similar to141-125 and similar to43 kyr. Growth at 141 kyr indicates temperatures >0 degreesC at a time when MIS 6 ice volumes were close to their maximum. High stable carbon isotope (delta(13)C) values (similar to2.8parts per thousand to +3.1parts per thousand) throughout the stalagmite's growth reflect a persistently low input of biogenic CO2, indicating that the steep, barren and alpine-like recharge area of today ha's been in existence for at least the last similar to380 kyr. During MIS 9, the lowest delta(13)C values occur well after maximum interglacial conditions, suggesting a lag in the development of post-glacial soils in this high-altitude karst. The stable oxygen isotope (delta(18)O) trends match the main structural features of the major climate proxy records (SPECMAP, Vostok and Devils Hole), suggesting that the delta(18)O of CC1 has responded to global-scale climate changes, whilst remarkable similarity exists between CC1 delta(18)O and regional sea-surface temperature reconstructions from North Atlantic core ODP980 and southwest Pacific marine core MD97-2120 through the most detailed part of the CC1 record, MIS 9-8. The results suggest that CC1 and other stalagmites from the cave have the potential to capture a long record of regional temperature trends, particularly in regards to the relative severity of Pleistocene glacial stages. (C) 2004 Elsevier B.V. All rights reserved.

Numerical investigation of seawater intrusion at Gooburrum, Bundaberg, Queensland, Australia

Seawater intrusion in coastal agricultural areas due to groundwater abstraction is a major environmental problem along the northeastern coast of Australia. Management options are being explored using numerical modelling, however, questions remain concerning the appropriate level of sophistication in models, choice of seaward boundary conditions, and how to accommodate heterogeneity and data uncertainty. The choice of seaward boundary condition is important since it affects the amount of salt transported into the aquifers and forms the focus of the present study. The impact of this boundary condition is illustrated for the seawater-intrusion problem in the Gooburrum aquifers, which occur within Tertiary sedimentary strata. A two-dimensional variable-density groundwater and solute-transport model was constructed using the computer code 2DFEMFAT (Cheng et al. 1998). The code was tested against an experiment for a steady-state freshwater-saltwater interface and against the Elder (Elder 1967) free-convection problem. Numerical simulations show that the imposition of the commonly-used equivalent hydrostatic freshwater heads, combined with a constant salt concentration at the seaward boundary, results in overestimated seawater intrusion in the lower Gooburrum aquifer. Since the imposition of this boundary condition allows water flow across the boundary, which subsequently takes salt into the aquifer, a careful check is essential to estimate whether too much mass of salt is introduced.

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.

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.