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.

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.

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.

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.

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate. (c) 2005 Elsevier Ltd. 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.