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.

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.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

CXTANNEAL: an improved program for estimating solute transport parameters

CXTANNEAL is a program for analysing contaminant transport in soils. The code, written in Fortran 77, is a modified version of CXTFIT, a commonly used package for estimating solute transport parameters in soils. The improvement of the present code is that it includes simulated annealing as the optimization technique for curve fitting. Tests with hypothetical data show that CXTANNEAL performs better than the original code in searching for optimal parameter estimates. To reduce the computational time, a parallel version of CXTANNEAL (CXTANNEAL_P) was also developed. (C) 1999 Elsevier Science Ltd. All rights reserved.

Turbulence in Small Sub-Tropical Estuary with Semi-Diurnal Tides

In natural estuaries, contaminant transport is driven by the turbulent momentum mixing. The predictions of scalar dispersion can rarely be predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted at high frequency and continuously for up to 50 hours per investigation in a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such small estuarine systems with shallow water depths (less than 0.5 m at low tides), and a thorough post-processing technique was applied. The estuarine flow is always a fluctuating process. The bulk flow parameters fluctuated with periods comparable to tidal cycles and other large-scale processes. But turbulence properties depended upon the instantaneous local flow properties. They were little affected by the flow history, but their structure and temporal variability were influenced by a variety of mechanisms. This resulted in behaviour which deviated from that for equilibrium turbulent boundary layer induced by velocity shear only. A striking feature of the data sets is the large fluctuations in all turbulence characteristics during the tidal cycle. This feature was rarely documented, but an important difference between the data sets used in this study from earlier reported measurements is that the present data were collected continuously at high frequency during relatively long periods. The findings bring new lights in the fluctuating nature of momentum exchange coefficients and integral time and length scales. These turbulent properties should not be assumed constant.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

Wave-induced seepage flux into anisotropic seabeds

Considerable effort has been devoted to quantifying the wave-induced soil response in a porous seabed in the last few decades. Most previous investigations have focused on the analysis of pore pressure and effective stresses within isotropic sediments, despite strong evidence of anisotropic soil behaviour reported in the literature. Furthermore, the seepage flux, which is important in the context of contaminant transport, has not been examined. In this paper, we focus on water wave-driven seepage in anisotropic marine sediments of finite thickness. The numerical results predict that the effects of hydraulic anisotropy and anisotropic soil behaviour on the wave-driven seepage in marine sediment are significant. Copyright (C) 2001 John Wiley & Sons, Ltd.

Tidal influence on behaviour of a coastal aquifer adjacent to a low-relief estuary

The tidal influence on groundwater hydrodynamics, salt-water intrusion and submarine groundwater discharge from coastal/estuarine aquifers is poorly quantified for systems with a mildly sloping beach, in contrast to the case where a vertical beach face is assumed. We investigated the effect of beach slope for a coastal aquifer adjacent to a low-relief estuary, where industrial waste was emplaced over the aquifer. The waste was suspected to discharge leachate towards the estuary. Field observations at various locations showed that tidally induced groundwater head fluctuations were skewed temporally. Frequency analysis suggested that the fluctuation amplitudes decreased exponentially and the phase-tags increased Linearly for the primary tidal signals as they propagated inland. Salinisation zones were observed in the bottom part of the estuary and near the beach surface. Flow and transport processes in a cross-section perpendicular to the estuary were simulated using SEAWAT-2000, which is capable of depicting density-dependent flow and multi-species transport. The simulations showed that the modelled water table fluctuations were in good agreement with the monitored data. Further simulations were conducted to gain insight into the effects of beach slope. In particular the limiting case of a vertical beach face was considered. The simulations showed that density difference and tidal forcing drive a more complex hydrodynamic pattern for the mildly sloping beach than the vertical beach, as well as a profound asymmetry in tidally induced water table fluctuations and enhanced salt-water intrusion. The simulation results also indicated that contaminant transport from the aquifer to the estuary was affected by the tide, where for the mildly sloping beach, the tide tended to intensify the vertical mass exchange in the vicinity of the shorelines, (c) 2005 Elsevier B.V. All rights reserved.

Modelling the fate of chromated copper arsenate in a sandy soil

A pulse of chromated copper arsenate (CCA, a timber preservative) was applied in irrigation water to an undisturbed field soil in a laboratory column. Concentrations of various elements in the leachate from the column were measured during the experiment. Also, the remnants within the soil were measured at the end of the experiment. The geochemical modelling package, PHREEQC-2, was used to simulate the experimental data. Processes included in the CCA transport modelling were advection, dispersion, non-specific adsorption (cation exchange) and specific adsorption by clay minerals and organic matter, as well as other possible chemical reactions such as precipitation/dissolution. The modelling effort highlighted the possible complexities in CCA transport and reaction experiments. For example, the uneven dosing of CCA as well as incomplete knowledge of the soil properties resulted in simulations that gave only partial, although reasonable, agreement with the experimental data. Both the experimental data and simulations show that As and Cu are strongly adsorbed and therefore, will mostly remain at the top of the soil profile, with a small proportion appearing in leachate. On the other hand, Cr is more mobile and thus it is present in the soil column leachate. Further simulations show that both the quantity of CCA added to the soil and the pH of the irrigation water will influence CCA transport. Simulations suggest that application of larger doses of CCA to the soil will result in higher leachate concentrations, especially for Cu and As. Irrigation water with a lower pH will dramatically increase leaching of Cu. These results indicate that acidic rainfall or significant accidental spillage of CCA will increase the risk of groundwater pollution.

Useful representations of series solutions for transport problems

Simple techniques are presented for rearrangement of an infinite series in a systematic way such that the convergence of the resulting expression is accelerated. These procedures also allow calculation of required boundary derivatives. Several examples of conduction and diffusion-reaction problems illustrate the methods.

Evaluating Flexible Transport Solutions

Flexible transport services (FTS) have been of increasing interest in developed countries as a bridge between the use of personal car travel and fixed route transit services. This paper reports on findings from a recent study in Queensland Australia, which identified lessons from an international review and implications for Australia. Potential strategic directions, including a vision, mission, key result areas, strategies, and identified means of measuring performance are described. Evaluation criteria for assessing flexible transport proposals were developed, and approaches to identifying and assessing needs and demands outlined. The use of emerging technologies is also a key element of successful flexible transport services.