Development of a mathematical and computer model to assess water quality impacts of hazardous compounds from minor gasoline spills in the inter-tidal zone of the Miami River, Florida

The objective of this study was to develop a model to predict transport and fate of gasoline components of environmental concern in the Miami River by mathematically simulating the movement of dissolved benzene, toluene, xylene (BTX), and methyl-tertiary-butyl ether (MTBE) occurring from minor gasoline spills in the inter-tidal zone of the river. Computer codes were based on mathematical algorithms that acknowledge the role of advective and dispersive physical phenomena along the river and prevailing phase transformations of BTX and MTBE. Phase transformations included volatilization and settling. ^ The model used a finite-difference scheme of steady-state conditions, with a set of numerical equations that was solved by two numerical methods: Gauss-Seidel and Jacobi iterations. A numerical validation process was conducted by comparing the results from both methods with analytical and numerical reference solutions. Since similar trends were achieved after the numerical validation process, it was concluded that the computer codes algorithmically were correct. The Gauss-Seidel iteration yielded at a faster convergence rate than the Jacobi iteration. Hence, the mathematical code was selected to further develop the computer program and software. The model was then analyzed for its sensitivity. It was found that the model was very sensitive to wind speed but not to sediment settling velocity. ^ A computer software was developed with the model code embedded. The software was provided with two major user-friendly visualized forms, one to interface with the database files and the other to execute and present the graphical and tabulated results. For all predicted concentrations of BTX and MTBE, the maximum concentrations were over an order of magnitude lower than current drinking water standards. It should be pointed out, however, that smaller concentrations than the latter reported standards and values, although not harmful to humans, may be very harmful to organisms of the trophic levels of the Miami River ecosystem and associated waters. This computer model can be used for the rapid assessment and management of the effects of minor gasoline spills on inter-tidal riverine water quality. ^

Remote sensing of tidal networks and their relation to vegetation

The study of the morphology of tidal networks and their relation to salt marsh vegetation is currently an active area of research, and a number of theories have been developed which require validation using extensive observations. Conventional methods of measuring networks and associated vegetation can be cumbersome and subjective. Recent advances in remote sensing techniques mean that these can now often reduce measurement effort whilst at the same time increasing measurement scale. The status of remote sensing of tidal networks and their relation to vegetation is reviewed. The measurement of network planforms and their associated variables is possible to sufficient resolution using digital aerial photography and airborne scanning laser altimetry (LiDAR), with LiDAR also being able to measure channel depths. A multi-level knowledge-based technique is described to extract networks from LiDAR in a semi-automated fashion. This allows objective and detailed geomorphological information on networks to be obtained over large areas of the inter-tidal zone. It is illustrated using LIDAR data of the River Ems, Germany, the Venice lagoon, and Carnforth Marsh, Morecambe Bay, UK. Examples of geomorphological variables of networks extracted from LiDAR data are given. Associated marsh vegetation can be classified into its component species using airborne hyperspectral and satellite multispectral data. Other potential applications of remote sensing for network studies include determining spatial relationships between networks and vegetation, measuring marsh platform vegetation roughness, in-channel velocities and sediment processes, studying salt pans, and for marsh restoration schemes.

The archaeogenetics of european ancestry

The archaeogenetics of Europe remains deeply controversial. Advances in ancient deoxyribonucleic acid (DNA) analysis have suggested gene flow between Neanderthals and modern humans, who arrived in Europe <50 000 years ago, but have so far failed to support evolution of Neanderthals from a population of Homo heidelbergensis represented by remains in northern Spain. The extent to which European Mesolithic forager populations versus Neolithic pioneers from the Near East contributed to the extant gene pool of Europeans also continues to be contested. Whilst analyses of extant mitochondrial lineages have emphasised late Palaeolithic and Mesolithic expansions, ancient DNA (aDNA) results suggest significant Neolithic dispersals from the southern ‘refugial’ zone into the northern ‘bio-tidal’ zone. However, whether these had a primarily Near Eastern or North Mediterranean source remains a matter for debate. Meanwhile, aDNA has also begun to highlight an important role for later dispersals, especially during the late Neolithic, in shaping the European gene pool.

Data assimilation for a coastal area morphodynamic model: Morecambe Bay

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.

Uses of airborne scanning laser altimetry for river flood modelling and coastal geomorphology

Two ongoing projects at ESSC that involve the development of new techniques for extracting information from airborne LiDAR data and combining this information with environmental models will be discussed. The first project in conjunction with Bristol University is aiming to improve 2-D river flood flow models by using remote sensing to provide distributed data for model calibration and validation. Airborne LiDAR can provide such models with a dense and accurate floodplain topography together with vegetation heights for parameterisation of model friction. The vegetation height data can be used to specify a friction factor at each node of a model’s finite element mesh. A LiDAR range image segmenter has been developed which converts a LiDAR image into separate raster maps of surface topography and vegetation height for use in the model. Satellite and airborne SAR data have been used to measure flood extent remotely in order to validate the modelled flood extent. Methods have also been developed for improving the models by decomposing the model’s finite element mesh to reflect floodplain features such as hedges and trees having different frictional properties to their surroundings. Originally developed for rural floodplains, the segmenter is currently being extended to provide DEMs and friction parameter maps for urban floods, by fusing the LiDAR data with digital map data. The second project is concerned with the extraction of tidal channel networks from LiDAR. These networks are important features of the inter-tidal zone, and play a key role in tidal propagation and in the evolution of salt-marshes and tidal flats. The study of their morphology is currently an active area of research, and a number of theories related to networks have been developed which require validation using dense and extensive observations of network forms and cross-sections. The conventional method of measuring networks is cumbersome and subjective, involving manual digitisation of aerial photographs in conjunction with field measurement of channel depths and widths for selected parts of the network. A semi-automatic technique has been developed to extract networks from LiDAR data of the inter-tidal zone. A multi-level knowledge-based approach has been implemented, whereby low level algorithms first extract channel fragments based mainly on image properties then a high level processing stage improves the network using domain knowledge. The approach adopted at low level uses multi-scale edge detection to detect channel edges, then associates adjacent anti-parallel edges together to form channels. The higher level processing includes a channel repair mechanism.

Progress on “Changing coastlines: data assimilation for morphodynamic prediction and predictability”

The task of assessing the likelihood and extent of coastal flooding is hampered by the lack of detailed information on near-shore bathymetry. This is required as an input for coastal inundation models, and in some cases the variability in the bathymetry can impact the prediction of those areas likely to be affected by flooding in a storm. The constant monitoring and data collection that would be required to characterise the near-shore bathymetry over large coastal areas is impractical, leaving the option of running morphodynamic models to predict the likely bathymetry at any given time. However, if the models are inaccurate the errors may be significant if incorrect bathymetry is used to predict possible flood risks. This project is assessing the use of data assimilation techniques to improve the predictions from a simple model, by rigorously incorporating observations of the bathymetry into the model, to bring the model closer to the actual situation. Currently we are concentrating on Morecambe Bay as a primary study site, as it has a highly dynamic inter-tidal zone, with changes in the course of channels in this zone impacting the likely locations of flooding from storms. We are working with SAR images, LiDAR, and swath bathymetry to give us the observations over a 2.5 year period running from May 2003 – November 2005. We have a LiDAR image of the entire inter-tidal zone for November 2005 to use as validation data. We have implemented a 3D-Var data assimilation scheme, to investigate the improvements in performance of the data assimilation compared to the previous scheme which was based on the optimal interpolation method. We are currently evaluating these different data assimilation techniques, using 22 SAR data observations. We will also include the LiDAR data and swath bathymetry to improve the observational coverage, and investigate the impact of different types of observation on the predictive ability of the model. We are also assessing the ability of the data assimilation scheme to recover the correct bathymetry after storm events, which can dramatically change the bathymetry in a short period of time.

Desiccation tolerance in Pleurostima purpurea (Velloziaceae)

Previous works suggested that Pleurostima purpurea (Velloziaceae-Barbacenioideae) shows a remarkable capacity to endure desiccation of its vegetative tissues. P. purpurea occurs in monocotyledons mats on soil islands in the Pao de Acucar (Sugar Loaf) one of the most recognizable rock outcrops of the world, in Rio de Janeiro, southeastern Brazil. Mats of P. purpurea occur in cliffs by the sea some meters above the tidal zone. Although living in rock outcrops almost devoid of any soil cover, P. purpurea seems to occur preferably on less exposed rock faces and slightly shady sites. Usually, less extreme adaptations to drought would be expected in plants with the habitat preference of P. purpurea. Relying on this observation, we argue if a combination of different strategies of dealing with low water availability can be found in P. purpurea as on other desiccation tolerant angiosperms. This study aims to examine the occurrence of desiccation tolerant behavior in P. purpurea together with the expression of drought avoidance mechanisms during dehydration progression. For this, it was analyzed the gas exchanges, leaf pigments and relative leaf water content during desiccation and rehydration of cultivated mature individuals. P. purpurea behaved like typical drought avoiders under moderated drought condition with stomatal closure occurring around a relative leaf water content up to 90%. During this process, it was observed a delay in the leaf relative water content (RWC(leaf)) decrease comparing to the plant-soil relative water content (RWC(plant-soil)). As soil dehydration worsened, gas exchanges restrictions progressed until a lack of activity which characterizes anabiosis. The loss of chlorophyll occurs before the end of total dehydration, characterizing the presence of poikilochlorophylly. The chlorophyll degradation follows the RWC(leaf) decrease, which achieved the minimum average value of 17% without incurring in leaf abscission. The chlorophyll re-synthesis seems to start well after the full rehydration of the leaf. During all of this process, carotenoid content remained stable. These results are coherent with a combination of drought avoidance and desiccation tolerance in P. purpurea which seems to be coherent with the amplitude of water availability in the rock outcrop habitat where it occurs, suggesting that the periods of water availability are sufficiently long for the success of the costly desiccation tolerant behavior but too short to make a typical drought avoider species win the competition for exploring the rock outcrop substrate where P. purpurea occurs.

Occurrence and behavior of butyltins in intertidal and shallow subtidal surface sediments of an estuarine beach under different sampling conditions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Occurrence and behavior of butyltins in intertidal and shallow subtidal surface sediments of an estuarine beach under different sampling conditions

Contamination by butyltin compounds (BTs) has been reported in estuarine environments worldwide, with serious impacts on the biota of these areas. Considering that BTs can be degraded by varying environmental conditions such as incident light and salinity, the short-term variations in such factors may lead to inaccurate estimates of BTs concentrations in nature. Therefore, the present study aimed to evaluate the possibility that measurements of BTs in estuarine sediments are influenced by different sampling conditions, including period of the day (day or night), tidal zone (intertidal or subtidal), and tides (high or low). The study area is located on the Brazilian southeastern coast, Sao Vicente Estuary, at Pescadores Beach, where BT contamination was previously detected. Three replicate samples of surface sediment were collected randomly in each combination of period of the day, tidal zone, and tide condition, from three subareas along the beach, totaling 72 samples. BTs were analyzed by GC-PFPD using a tin filter and a VF-5 column, by means of a validated method. The concentrations of tributyltin (TBT), dibutyltin (DBT), and monobutyltin (MBT) ranged from undetectable to 161 ng Sn g(-1) (d.w.). In most samples (71%), only MBT was quantifiable, whereas TBTs were measured in only 14, suggesting either an old contamination or rapid degradation processes. DBT was found in 27 samples, but could be quantified in only one. MBT concentrations did not differ significantly with time of day, zones, or tide conditions. DBT and TBT could not be compared under all these environmental conditions, because only a few samples were above the quantification limit. Pooled samples of TBT did not reveal any difference between day and night. These results indicated that, in assessing contamination by butyltin compounds, surface-sediment samples can be collected in any environmental conditions. However, the wide variation of BTs concentrations in the study area, i.e., over a very small geographic scale, illustrates the need for representative hierarchical and composite sampling designs that are compatible with the multiscalar temporal and spatial variability common to most marine systems. The use of such sampling designs will be necessary for future attempts to quantitatively evaluate and monitor the occurrence and impact of these compounds in nature

Iodine speciation in atmospheric aerosols in the marine boundary layer

Iodine chemistry plays an important role in the tropospheric ozone depletion and the new particle formation in the Marine Boundary Layer (MBL). The sources, reaction pathways, and the sinks of iodine are investigated using lab experiments and field observations. The aims of this work are, firstly, to develop analytical methods for iodine measurements of marine aerosol samples especially for iodine speciation in the soluble iodine; secondly, to apply the analytical methods in field collected aerosol samples, and to estimate the characteristics of aerosol iodine in the MBL. Inductively Coupled Plasma – Mass Spectrometry (ICP-MS) was the technique used for iodine measurements. Offline methods using water extraction and Tetra-methyl-ammonium-hydroxide (TMAH) extraction were applied to measure total soluble iodine (TSI) and total insoluble iodine (TII) in the marine aerosol samples. External standard calibration and isotope dilution analysis (IDA) were both conducted for iodine quantification and the limits of detection (LODs) were both 0.1 μg L-1 for TSI and TII measurements. Online couplings of Ion Chromatography (IC)-ICP-MS and Gel electrophoresis (GE)-ICP-MS were both developed for soluble iodine speciation. Anion exchange columns were adopted for IC-ICP-MS systems. Iodide, iodate, and unknown signal(s) were observed in these methods. Iodide and iodate were separated successfully and the LODs were 0.1 and 0.5 μg L-1, respectively. Unknown signals were soluble organic iodine species (SOI) and quantified by the calibration curve of iodide, but not clearly identified and quantified yet. These analytical methods were all applied to the iodine measurements of marine aerosol samples from the worldwide filed campaigns. The TSI and TII concentrations (medians) in PM2.5 were found to be 240.87 pmol m-3 and 105.37 pmol m-3 at Mace Head, west coast of Ireland, as well as 119.10 pmol m-3 and 97.88 pmol m-3 in the cruise campaign over the North Atlantic Ocean, during June – July 2006. Inorganic iodine, namely iodide and iodate, was the minor iodine fraction in both campaigns, accounting for 7.3% (median) and 5.8% (median) in PM2.5 iodine at Mace Head and over the North Atlantic Ocean, respectively. Iodide concentrations were higher than iodate in most of the samples. In the contrast, more than 90% of TSI was SOI and the SOI concentration was correlated significantly with the iodide concentration. The correlation coefficients (R2) were both higher than 0.5 at Mace Head and in the first leg of the cruise. Size fractionated aerosol samples collected by 5 stage Berner impactor cascade sampler showed similar proportions of inorganic and organic iodine. Significant correlations were obtained in the particle size ranges of 0.25 – 0.71 μm and 0.71 – 2.0 μm between SOI and iodide, and better correlations were found in sunny days. TSI and iodide existed mainly in fine particle size range (< 2.0 μm) and iodate resided in coarse range (2.0 – 10 μm). Aerosol iodine was suggested to be related to the primary iodine release in the tidal zone. Natural meteorological conditions such as solar radiation, raining etc were observed to have influence on the aerosol iodine. During the ship campaign over the North Atlantic Ocean (January – February 2007), the TSI concentrations (medians) ranged 35.14 – 60.63 pmol m-3 among the 5 stages. Likewise, SOI was found to be the most abundant iodine fraction in TSI with a median of 98.6%. Significant correlation also presented between SOI and iodide in the size range of 2.0 – 5.9 μm. Higher iodate concentration was again found in the higher particle size range, similar to that at Mace Head. Airmass transport from the biogenic bloom region and the Antarctic ice front sector was observed to play an important role in aerosol iodine enhancement. The TSI concentrations observed along the 30,000 km long cruise round trip from East Asia to Antarctica during November 2005 – March 2006 were much lower than in the other campaigns, with a median of 6.51 pmol m-3. Approximately 70% of the TSI was SOI on average. The abundances of inorganic iodine including iodine and iodide were less than 30% of TSI. The median value of iodide was 1.49 pmol m-3, which was more than four fold higher than that of iodate (median, 0.28 pmol m-3). Spatial variation indicated highest aerosol iodine appearing in the tropical area. Iodine level was considerably lower in coastal Antarctica with the TSI median of 3.22 pmol m-3. However, airmass transport from the ice front sector was correlated with the enhance TSI level, suggesting the unrevealed source of iodine in the polar region. In addition, significant correlation between SOI and iodide was also shown in this campaign. A global distribution in aerosol was shown in the field campaigns in this work. SOI was verified globally ubiquitous due to the presence in the different sampling locations and its high proportion in TSI in the marine aerosols. The correlations between SOI and iodide were obtained not only in different locations but also in different seasons, implying the possible mechanism of iodide production through SOI decomposition. Nevertheless, future studies are needed for improving the current understanding of iodine chemistry in the MBL (e.g. SOI identification and quantification as well as the update modeling involving organic matters).

Corrosión de armaduras en el hormigón armado en ambiente marino aéreo

La presente tesis doctoral surge ante la necesidad de completar el estudio sobre la corrosión en el hormigón armado expuesto al ambiente marino iniciado por D. Miguel Ángel Bermúdez Odriozola en 2007 con la tesis doctoral con título “Corrosión de las armaduras del hormigón armado en ambiente marino: zona de carrera de mareas y zona sumergida”. Una vez realizado el estudio bibliográfico, se ha planteado un estudio experimental en el que se analizan cuatro muelles de distintas edades y localizados en distintos lugares de la geografía española. Uno de ellos contiene adiciones en el hormigón. Una vez inspeccionados, se extraen testigos de cada uno de ellos con el fin de realizar ensayos de caracterización y se analizan sus resultados en paralelo con los coeficientes de difusión. Uno se los objetivos ha sido comprobar el mejor método de ensayo que permita controlar la durabilidad del hormigón en ambiente IIIa. La tesis que se presenta ha demostrado la efectividad del ensayo de penetración de agua para discriminar los hormigones que presentarán un buen comportamiento frente a la difusión de cloruros y por tanto para predecir la vida útil de la estructura. Asimismo, otro objetivo ha sido valorar el uso de adiciones en el hormigón y cómo afectan a la mejora del comportamiento del hormigón expuesto al ambiente marino. Se han obtenido valores del coeficiente de eficacia de cenizas volantes y humo de sílice frente a la difusión de cloruros, muy superiores a los actualmente considerados en términos de resistencia. El ambiente marino aéreo engloba estructuras situadas en primera línea de costa así como aquellas que se sitúan hasta 5km de distancia de la línea del mar. La investigación realizada ha demostrado la necesidad de establecer una nueva subdivisión en este ambiente introduciendo la zona Spray, siguiendo criterios internacionales, que permita diferenciar la zona más desfavorable del mismo que comprende las estructuras ubicadas sobre el nivel del mar tal como los tableros de pantalanes portuarios estudiados. Se han analizado en profundidad los valores del contenido de cloruros en superficie de hormigón en los diferentes tipos de ambiente marino recogidos por la Instrucción EHE- 08 contrastando los requisitos normativos. Como conclusión se han hecho nuevas propuestas normativas que afectan al ambiente sumergido y al ambiente marino aéreo (zona spray). Asimismo se ha realizado un análisis en profundidad de la concentración crítica de cloruros recogida por la Instrucción EHE-08, a partir de los datos obtenidos en estructuras reales. Se ha evaluado el ensayo de resistividad para estimar el valor del coeficiente de difusión de cloruros. Este ensayo ha demostrado su efectividad para discriminar hormigones con adiciones y con muy baja penetración de agua, si bien no resulta sensible para detectar hormigones con elevada permeabilidad. Finalmente en la tesis se ha obtenido un modelo de difusión de cloruros (a/c-D) en ambiente marino aéreo a fin de poder estimar la difusión de cloruros en el hormigón a partir de un requisito normativo como es la relación a/c. This Thesis is focussed to complete the previous study on the corrosion of reinforced concrete exposed to marine environment, carried out by Dr. Miguel Ángel Bermúdez Odrizola in 2007. That Thesis title was “Corrosion of concrete reinforcement located in marine environment: tidal zone and submerged zone”. After the literature review presented in the State of the Art, an experimental programme was carried out. Four docks of different ages and locations have been studied, one of them was a fly ash concrete. The structures were inspected on site and several cores extracted in order to develop characterization tests and obtain diffusion profiles at the laboratory. One of the objectives has been to check the best laboratory test method to control concrete durability when exposed to atmospheric marine environment. This Thesis has proven the effectiveness of the water penetration test to differentiate concretes well behaved against chloride diffusion and also to predict the structure service life. Another objective has been to study the effect of mineral additions in concrete and how they improve the durability of concrete in this environment. Efficiency factors have been obtained for fly ashes and silica fume, much higher than those consider for strength. According to the Spanish Concrete Code EHE 08, the atmospheric marine environment includes structures located at the cost, from the sea line up to 5km inland. The research conducted has shown the need to consider a spray zone which is the most unfavourable area of the atmospheric environment, affecting structures over the sea level as the decks of the wharves here studied. Surface chloride contents of structures in the different marine environments have been deeply analysed, checking the actual standard limits. As a conclusion, new normative values have been proposed for submerged and atmospheric zone (spray zone). A complete analysis of the chloride threshold concentration from the data obtained in real structures has also been performed. It has been evaluated the resistivity test as a tool to estimate the chloride diffusion coefficient. This test has proven effectiveness to discriminate concrete with mineral additions and low water penetration values, while it is not sensitive to detect concretes with high permeability. Finally, a chloride diffusion model (w/c) has been obtained in atmospheric marine environment in order to estimate the chloride diffusion of the concrete knowing a normative requirement as w/c.

Age determination of Holozene marine deposites from the western Mauritania (Table 1)

In an area regarded to be very favourable for the study of Holocene sea level changes one or several eustatic (?) oscillations of sea have been found using sedimentological and ecological methods. After a maximum of +3 m during the Nouakchottian stage (= Middle Flandrian or Late Atlantic) about 5500 YBP a drop of sea to -3.5 ± 0.5 m about 4100 YBP is testified by stromatolitic algae indicating the former sea level within the tidal zone with high accuracy. This evidence is supported by the observation of post-Nouakchottian regressive and transgressive geologic sequences, by buried beach deposits and flooded hardgrounds, post-Nouakchottian marine terraces of different height and age, the cutting off of one large and several small bays from the open sea etc. Possibly, one or two smaller oscillations followed between 4000 and 1500 years B. P. The radiocarbon age of the marine shells dated may be partly somewhat too old or too young.

High-Frequency Turbulence and Sediment Flux Measurements in an Upper Estuarine Zone

In natural estuaries, the predictions of scalar dispersion are rarely predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted continuously at high frequency for 50 hours in the upper zone of a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such shallow water depths (less than 0.4 m at low tides), and a thorough post-processing technique was applied. In addition, some experiments were conducted in laboratory under controlled conditions using water and soil samples collected in the estuary to test the relationship between acoustic backscatter strength and suspended sediment load. A striking feature of the field data set was the large fluctuations in all turbulence characteristics during the tidal cycle, including the suspended sediment flux. This feature was rarely documented.

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.