Using analytic element models to delineate drinking water source protection areas

Since 1999, Ohio EPA hydrogeologists have used two analytic element models (AEMs), the proprietary software GFLOW and U.S. EPA's WhAEM, to delineate protection areas for 535 public water systems. Both models now use the GFLOW2001 solution engine, integrate well with Geographic Information System (GIS) technology, have a user-friendly graphical interface, are capable of simulating a variety of complex hydrogeologic settings, and do not rely upon a model grid. These features simplify the modeling process and enable AEMs to bridge the gap between existing simplistic delineation methods and more complex numerical models. Ohio EPA hydrogeologists demonstrated that WhAEM2000 and GFLOW2000 were capable of producing capture zones similar to more widely accepted models by applying the AEMs to eight sites that had been previously delineated using other methods. After the Ohio EPA delineated protection areas using AEMs, more simplistic delineation methods used by other states (volumetric equation and arbitrary fixed radii) were applied to the same water systems to compare the differences between various methods. GIS software and two-tailed paired t-tests were used to quantify the differences in protection areas and analyze the data. The results of this analysis demonstrate that AEMs typically produce significantly different protection areas than the most simplistic delineation methods, in terms of total area and shape. If the volumetric equation had been used instead of AEMs, Ohio would not have protected 265 km2 of critical upgradient area and would have overprotected 269 km2 of primarily downgradient land. Since an increasing number of land-use restrictions are being tied to drinking water protection areas, this analysis has broad policy implications.     

Positive solution of two-dimensional solute transport in heterogeneous aquifers

The transport of contaminants in aquifers is usually represented by a convection-dispersion equation. There are several well-known problems of oscillation and artificial dispersion that affect the numerical solution of this equation. For example, several studies have shown that standard treatment of the cross-dispersion terms always leads to a negative concentration. It is also well known that the numerical solution of the convective term is affected by spurious oscillations or substantial numerical dispersion. These difficulties are especially significant for solute transport in nonuniform flow in heterogeneous aquifers. For the case of coupled reactive-transport models, even small negative concentration values can become amplified through nonlinear reaction source/sink terms and thus result in physically erroneous and unstable results. This paper includes a brief discussion about how nonpositive concentrations arise from numerical solution of the convection and cross-dispersion terms. We demonstrate the effectiveness of directional splitting with one-dimensional flux limiters for the convection term. Also, a new numerical scheme for the dispersion term that preserves positivity is presented. The results of the proposed convection scheme and the solution given by the new method to compute dispersion are compared with standard numerical methods as used in MT3DMS.     

Magnetic cleansing of weathered/tarry oiled feathers--the role of pre-conditioners

Iron powder has previously been demonstrated to be effective in the removal, via magnetic harvesting, of a wide variety of oil contaminants from feathers and plumage. This study investigates the efficacy of magnetic cleansing for the removal from feathers of tarry contamination that has been allowed to weather. Clusters of feathers from Mallard duck (Anas platyrhnchos) and Little Penguin (Eudyptula minor) were completely immersed in a tarry contaminant and allowed to weather from one to fourteen days. The contaminant was removed using a magnetic cleansing protocol and the removal efficacy assessed gravimetrically. For one, seven and fourteen days of weathering, a final removal (after fourteen treatments) of more than 99% and 97% was achieved for duck feathers and penguin feathers, respectively. Repeating the experiments (for a seven-day weathering period) for both duck and penguin feathers, with the judicious application of a pre-conditioner (olive oil), further improved removal efficacy. A convenient method to screen for improved pre-conditioning agents is suggested.     

Multiscale characterization of a heterogeneous aquifer using an ASR operation

Heterogeneity in the physical properties of an aquifer can significantly affect the viability of aquifer storage and recovery (ASR) by reducing the recoverable proportion of low-salinity water where the ambient ground water is brackish or saline. This study investigated the relationship between knowledge of heterogeneity and predictions of solute transport and recovery efficiency by combining permeability and ASR-based tracer testing with modeling. Multiscale permeability testing of a sandy limestone aquifer at an ASR trial site showed that small-scale core data give lower-bound estimates of aquifer hydraulic conductivity (K), intermediate-scale downhole flowmeter data offer valuable information on variations in K with depth, and large-scale pumping test data provide an integrated measure of the effective K that is useful to constrain ground water models. Chloride breakthrough and thermal profiling data measured during two cycles of ASR showed that the movement of injected water is predominantly within two stratigraphic layers identified from the flowmeter data. The behavior of the injectant was reasonably well simulated with a four-layer numerical model that required minimal calibration. Verification in the second cycle achieved acceptable results given the model's simplicity. Without accounting for the aquifer's layered structure, high precision could be achieved on either piezometer breakthrough or recovered water quality, but not both. This study demonstrates the merit of an integrated approach to characterizing aquifers targeted for ASR.     

Conventional and combined pump-and-treat systems under nonuniform background flow

We investigate the performance of vertical hydraulic barriers in combination with extraction wells for the partial hydraulic isolation of contaminated aquifer areas. The potential advantage of such combinations compared to a conventional pump-and-treat system has already been demonstrated in a previous study. Here we extend the scope of the performance analysis to the impact of uncertainty in the regional flow direction as well as to highly heterogeneous aquifer transmissivity distributions. In addition, two new well-barrier scenarios are proposed and analyzed. The hydraulic efficiency of the scenarios is rated based on the expected (mean) reduction of the pumping rate that is required to achieve downgradient contaminant capture. The uncertain spatial distribution of aquifer transmissivity is considered by means of unconditioned Monte Carlo simulations. The significance of uncertain background flow conditions is incorporated by computing minimized pumping rates for deviations of the regional flow direction up to 30 degrees from a normative base case. The results give an answer on how pumping rates have to be changed for each barrier-well combination in order to achieve robust systems. It is exposed that in comparison to installing exclusively wells, the barrier-supported approach generally yields savings in the (average) pumping rate. The particular efficiency is shown to be highly dependent on the interaction of variance and integral scale of transmissivity distribution, well and barrier position, as well as direction of background flow.     

Macro- and meiofaunal recolonisation of dredged material used for habitat enhancement: temporal patterns in community development

In recent years, dredged material has become regarded as a potential resource and used to create and/or improve intertidal habitats (i.e., beneficial use). This paper presents the results of a sampling programme to investigate the long-term (42 months post-recharge) macro- and meiofaunal recolonisation processes of a beneficial use scheme in south-east England. While univariate indices of community structure indicated that the scheme’s meiofaunal community was never significantly different from that of a nearby reference area, such attributes for macrofauna were continually significantly below those of the reference area, although this was not the case for all reference stations. Multivariate analyses revealed that macro- and meiofaunal community structures were always significantly different from those of the reference communities. We discuss the factors responsible for these observations and propose that assessing recovery of a beneficial use scheme should be undertaken using pre-defined criteria in addition to comparisons with a reference site.     

Atmospheric lead fallout over the last century recorded in Gulf of Lions sediments (Mediterranean Sea)

Six marine sediment cores from the Gulf of Lions continental slope (700-1700 m water depth) were analyzed for stable lead isotopes and (210)Pb geochronology in order to reconstruct lead atmospheric fallout pattern during the last century. The detrital lead contribution is 25 microg g(-1) and the mean sediment anthropogenic inventory is 110+/-7 microg cm(-2), a little bit higher than atmospheric deposition estimate. Anthropogenic lead accumulation in sediments peaked in early 1970s (1973+/-2) in agreement with lead emissions features. For the period 1986-1997, the sediment signal also reflect the decrease of atmospheric lead described by independent atmospheric fallout investigations. The anthropogenic Pb deposition in the late 1990s was similar to the 1950s deposition, attesting thus of the output of European environmental policies.     

The use of the marine biotic index AMBI in the assessment of the ecological status of the Obidos lagoon (Portugal)

This study reports the longitudinal distribution of sediment properties, including inorganic and organic contaminants, and the structure of the benthic community in Obidos lagoon, a coastal system permanently connected to the sea and with negligible freshwater sources. Sediments from the upper to central lagoon consist of fine particles (91%) and from the lower lagoon of sands (94%). Chemical composition is strongly correlated to the percentage of fine particles. Contamination is relatively low in those sediments suggesting the effect of diffuse sources. The increase in organic matter content from down- to upstream areas was associated with the dominance of opportunistic species, while sensitive and indifferent species to organic enrichment were mainly associated to the clean sandy downstream area. The marine biotic index (AMBI) was suitable for the discrimination of the biological and environmental gradients in the Obidos lagoon and was highly related with the gradient of organic matter content in this system.     

A 3D hydrodynamic fate and transport model for herbicides in Sacca di Goro coastal lagoon (Northern Adriatic)

Sacca di Goro is a shallow coastal microtidal lagoon with a surface area of 26 km2, and an average depth of about 1.5m. Fresh water pollutant loads from Po River branches and several drainage canals lead to anthropogenic eutrophication, frequent summer anoxia crises and chemical contamination. Such events not only affect the lagoon ecosystem but also cause serious economic losses, the lagoon being the second largest producer of clams in Italy. The present work aims at using a fate model coupled with COHERENS 3D hydrodynamic model to simulate and to explain the spatial distribution and temporal variations of s-triazines herbicides in the Sacca di Goro lagoon. The simulation results of spatial and temporal dynamic behaviour of atrazine, simazine and terbuthylazine have been compared with experimental data obtained during an annual monitoring programme.     

Are Taxonomic Distinctness measures compliant to other ecological indicators in assessing ecological status?

Assessing the ecological status, a concept implemented in the European Water Framework Directive [EC, 2000. Directive of the European Parliament and of the Council 2000/60/EC establishing a framework for community action in the field of water policy PE-CONS 3639/1/00, p. 72], requires the application of methods capable of distinguishing different levels of ecological quality. The Average Taxonomic Distinctness has been used as tool in this context, and we tested the robustness of Taxonomic Distinctness measures applying it in different scenarios (estuarine eutrophication, organic pollution, and re-colonisation after physical disturbance), analysing simultaneously its compliance to other types of ecological indicators. Results show that, in most of the case studies, only Total Taxonomic Distinctness was relatively satisfactory in discriminating between disturbed situations. Other Taxonomic Distinctness measures have not proved to be more sensitive than other ecological indicators (Shannon-Wiener, Margalef, and Eco-Exergy indices). Therefore, this approach does not seem to be particularly helpful in assessing systems' ecological status with regard to the WFD implementation.