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.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Analytical solution for tidal propagation in a coupled semi-confined/phreatic coastal aquifer

An analytical solution is derived for tidal fluctuations in a coupled coastal aquifer system consisting of a semi-confined aquifer, a thin semi-permeable layer and a phreatic aquifer. Based on the solution, we study the interactions (via leakage) between the confined and unconfined aquifers in response to tides. The results show that, under certain conditions, leakage from the confined aquifer can affect considerably the tidal water table fluctuation in the phreatic aquifer and vice versa. Ignoring these effects could lead to errors in estimating aquifer properties based on tidal signals. (C) 2002 Elsevier Science Ltd. All rights reserved.

Two-dimensional approximation for tidal dynamics in coastal aquifers: Capillarity correction

Most previous investigations on tide-induced watertable fluctuations in coastal aquifers have been based on one-dimensional models that describe the processes in the cross-shore direction alone, assuming negligible along-shore variability. A recent study proposed a two-dimensional approximation for tide-induced watertable fluctuations that took into account coastline variations. Here, we further develop this approximation in two ways, by extending the approximation to second order and by taking into account capillary effects. Our results demonstrate that both effects can markedly influence watertable fluctuations. In particular, with the first-order approximation, the local damping rate of the tidal signal could be subject to sizable errors.

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.

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.

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.

Hydraulic tracer studies in a pilot scale subsurface flow constructed wetland

Current design procedures for Subsurface Flow (SSF) Wetlands are based on the simplifying assumptions of plug flow and first order decay of pollutants. These design procedures do yield functional wetlands but result in over-design and inadequate descriptions of the pollutant removal mechanisms which occur within them. Even though these deficiencies are often noted, few authors have attempted to improve modelling of either flow or pollutant removal in such systems. Consequently the Oxley Creek Wetland, a pilot scale SSF wetland designed to enable rigorous monitoring, has recently been constructed in Brisbane, Australia. Tracer studies have been carried out in order to determine the hydraulics of this wetland prior to commissioning it with sealed sewage. The tracer studies will continue during the wetland's commissioning and operational phases. These studies will improve our understanding of the hydraulics of newly built SSF wetlands and the changes brought on by operational factors such as biological films and wetland plant root structures. Results to date indicate that the flow through the gravel beds is not uniform and cannot be adequately modelled by a single parameter, plug flow with dispersion, model. We have developed a multiparameter model, incorporating four plug flow reactors, which provides a better approximation of our experimental data. With further development this model will allow improvements to current SSF wetland design procedures and operational strategies, and will underpin investigations into the pollutant removal mechanisms at the Oxley Creek Wetland. (C) 1997 IAWQ. Published by Elsevier Science Ltd.

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.

Critical review of acid in situ leach uranium mining: 2. Soviet Block and Asia

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 sulphuric acid chemistry and no restoration of groundwater following mining. ISL mines in the former Soviet Union generally used acid reagents and were operated without due consideration given to environmental protection. At many former mine sites, the extent of groundwater contamination is significant because of high salinity, heavy metal and radionuclide concentrations compared with pre-mining and changes in the hydrogeological regime caused by mining. After the political collapse of the Soviet Union by the early 1990s, most uranium mines were shut down or ordered to be phased out by government policy. Programmes of restoration are now being undertaken but are proving technically difficult and hampered by a lack of adequate financial resources. The history and problems of acid ISL sites in countries of the former Soviet Union and Asia are presented in this study.

Can the nitrogenous composition of xylem sap be used to assess salinity stress in Casuarina glauca?

It is predicted that dryland salinity will affect up to 17 Mha of the Australian landscape by 2050, and therefore, monitoring the health of tree plantings and remnant native vegetation in saline areas is increasingly important. Casuarina glauca Sieber ex Spreng. has considerable salinity tolerance and is commonly planted in areas with a shallow, saline water table. To evaluate the potential of using the nitrogenous composition of xylem sap to assess salinity stress in C. glauca, the responses of trees grown with various soil salinities in a greenhouse were compared with those of trees growing in field plots with different water table depths and groundwater salinities. In the greenhouse, increasing soil salinity led to increased allocation of nitrogen (N) to proline and arginine in both stem and root xylem sap, with coincident decreases in citrulline and asparagine. Although the field plots were ranked as increasingly saline-based on ground water salinity and depth-only the allocation of N to citrulline differed significantly between the field plots. Within each plot, temporal variation in the composition of the xylem sap was related to rainfall, rainfall infiltration and soil salinity. Periods of low rainfall and infiltration and higher soil salinity corresponded with increased allocation of N to proline and arginine in the xylem sap. The allocation of N to citrulline and asparagine increased following rainfall events where rain was calculated to have infiltrated sufficiently to decrease soil salinity. The relationship between nitrogenous composition of the xylem sap of C. glauca and soil salinity indicates that the analysis of xylem sap is an effective method for assessing changes in salinity stress in trees at a particular site over time. However, the composition of the xylem sap proved less useful as a comparative index of salinity stress in trees growing at different sites.

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.

Confined-unconfined flow in a horizontal confined aquifer

Mixed confined and unconfined groundwater flow occurs in a bounded initially dry aquifer when the hydraulic head at the side boundary suddenly rises above the elevation of the aquifer's top boundary. The flow problem as modelled by the Boussinesq equation is non-trivial because of the involvement of two moving boundaries. The transformed equation (based on a similarity transformation) can, however, be dealt with more easily. Here, we present an approximate analytical solution for this flow problem. The approximate solution is compared with an 'exact' numerical solution and found to be a very accurate description for describing the mixed confined and unconfined flow in the confined aquifer. (C) 2002 Elsevier Science B.V. All rights reserved.

Effect of soil composition and dissolved organic matter on pesticide sorption

The effect of the solid and dissolved organic matter fractions, mineral composition and ionic strength of the soil solution on the sorption behaviour of pesticides were studied. A number of soils, chosen so as to have different clay mineral and organic carbon content, were used to study the sorption of the pesticides atrazine (6-chloro-N-2-ethyl-N-4-isopropyl-1,3,5-triazine-2,4-diamine), 2,4-D ((2,4-dichlorophenoxy) acetic acid), isoproturon (3-(4-isopropylphenyl)1,1-dimethylurea) and paraquat (1,1'-dimethyl-4,4'-bipyridinium) in the presence of low and high levels of dissolved organic carbon and different background electrolytes. The sorption behaviour of atrazine, isoproturon and paraquat was dominated by the solid state soil components and the presence of dissolved organic matter had little effect. The sorption of 2,4-D was slightly affected by the soluble organic matter in the soil. However, this effect may be due to competition for adsorption sites between the pesticide and the soluble organic matter rather than due to a positive interaction between the pesticide and the soluble fraction of soil organic matter. It is concluded that the major factor governing the sorption of these pesticides is the solid state organic fraction with the clay mineral content also making a significant contribution. The dissolved organic carbon fraction of the total organic carbon in the soil and the ionic strength of the soil solution appear to have little or no effect on the sorption/transport characteristics of these pesticides over the range of concentrations studied. (C) 2002 Elsevier Science B.V. All rights reserved.

Influence of hysteresis on tidal capillary fringe dynamics in a well-sorted sand

Nielsen and Perrochet [Adv. Water Resour. 23 (2000) 503] presented experimental data for cyclic water movement in the vadose zone above an oscillating watertable. The response of the watertable to cyclic forcing was characterised by the ratios of the forcing head to watertable amplitudes and their associated phase lag. They found that their non-hysteretic Richards' equation model failed to represent the observed behaviour of these parameters. This paper explores the effect on the simulated capillary fringe dynamics (in terms of these parameters) of including varying degrees of hysteresis in the moisture retention curve used in a numerical model of their experiment. It is clear that hysteresis can indeed account for observed discrepancies between simulation and experiment and that the effect of hysteresis varies with the frequency of oscillation. The use of a single-valued mean retention curve, as advocated by some authors, fails to provide a match between the simulated and observed behaviour of the Nielsen and Perrochet parameters, but is shown to be adequate for predicting time-averaged soil moisture profiles. (C) 2003 Elsevier Ltd. All rights reserved.