Fate of arsenic, phosphate and ammonium plumes in a coastal aquifer affected by saltwater intrusion.

A severe groundwater contamination with extensive plumes of arsenic, phosphate and ammonium was found in a coastal aquifer beneath a former fertilizer production plant. The implementation of an active groundwater remediation strategy, based on a comprehensive pump and treat scheme, now prevents the migration of the dissolved contaminants into the marine environment. However, due to the site's proximity to the coastline, a seawater wedge was induced by the pumping scheme. Additionally the groundwater flow and salinity patterns were also strongly affected by leakage from the site's sewer system and from a seawater-fed cooling canal. The objective of this study was to elucidate the fate of arsenic and its co-contaminants over the site's history under the complex, coupled hydrodynamic and geochemical conditions that prevail at the site. A detailed geochemical characterisation of samples from sediment cores and hydrochemical data provided valuable high-resolution information. The obtained data were used to develop various conceptual models and to constrain the development and calibration of a reactive transport model. The reactive transport simulations were performed for a sub-domain (two-dimensional transect) of an earlier developed three-dimensional flow and variable density solute transport model. The results suggest that in the upper sub-oxic zone the influx of oxygenated water promoted As attenuation via co-precipitation with Al and Fe oxides and copper hydroxides. In contrast, in the deeper aquifer zone, iron reduction, associated with the release of adsorbed As and the dissolution of As bearing phases, provided and still provides to date a persistent source for groundwater pollution. The presented monitoring and modelling approach could be broadly applied to coastal polluted sites by complex contaminant mixture containing As.

Performance of different assessment methods to evaluate contaminant sources and fate in a coastal aquifer.

The present study deals with the application of different monitoring techniques and numerical models to characterize coastal aquifers affected by multiple sources of contamination. Specifically, equivalent freshwater heads in 243 monitoring wells were used to reconstruct the piezometric map of the studied aquifer; flow meter tests were carried out to infer vertical groundwater fluxes at selected wells; deuterium and oxygen isotopes were used to identify the groundwater origin, and tritium was analyzed to estimate the residence time; compound-specific isotope analyses and microbial analyses were employed to track different sources of contamination and their degradation; numerical modelling was used to estimate and verify groundwater flow direction and magnitude throughout the aquifer. The comparison of the information level for each technique allowed determining which of the applied approaches showed the best results to locate the possible sources and better understanding of the fate of the contaminants. This study reports a detailed site characterization process and outcomes for a coastal industrial site, where a comprehensive conceptual model of pollution and seawater intrusion has been built using different assessment methods. Information and results from this study encourages combining different methods for the design and implementation of the monitoring activities in real-life coastal contaminated sites in order to develop an appropriate strategy for control and remediation of the contamination.

Reactive modelling of 1,2-DCA and DOC near the shoreline.

1,2-Dichloroethane (1,2-DCA) was found to be the most abundant compound among chlorinated hydrocarbons detected in a petrochemical plant in southern Italy. This site is located near the coastline, and it is set above an unconfined coastal aquifer, where seawater intrusion is present. The presence of organic and inorganic contaminants at this site has required the implementation of remediation strategies, consisting of pumping wells (hydraulic barrier) and a horizontal flow barrier. The purpose of this work was to assess the influence of salt water intrusion on the degradation rate of 1,2-DCA. This was done on a three-dimensional domain relative to a limited portion of a well characterized field site, accounting for density-dependent flow and reactive transport modelling of 1,2-DCA and Dissolved Organic Carbon (DOC). The modelling procedure was performed employing SEAWAT-4.0 and PHT3D, to reproduce the complex three-dimensional flow and transport domain. In order to determine the fate of 1,2-DCA, detailed field investigations provided intensive depth profile information. Different, kinetically controlled degradation rates were simulated to explain the observed, selective degradation of pollutants in groundwater. Calibration of the model was accomplished by comparison with the two different sets of measurements obtained from the MLS devices and from pumping wells. With the calibrated model, it was possible to distinguish between dispersive non-reactive processes and bacterially mediated reactions. In the non-reactive model, 1,2-DCA sorption was simulated using linear sorption coefficient determined with field data and 1,2-DCA degradation was simulated using a first order decay coefficient using literature data as initial guess. Finally, on the reactive transport model, where a two-step approach with partial equilibrium approach was implemented, the effects of neglecting the cation exchange capacity, omitting density-dependent flow, and refining the vertical discretization of the model were investigated. Comparison of results from various scenarios shows that geochemical changes in inorganic constituents can be used to improve the site's conceptual model, and establishes that natural degradation processes can be suitable for 1,2-DCA as a remediation option.

Modelling the fate of styrene in a mixed petroleum hydrocarbon plume.

Severe petroleum hydrocarbon contamination (styrene and the BTEX compounds: benzene, toluene, ethylbenzene and the isomers of xylene) from leaking sewers was detected in a Quaternary aquifer below a chemical plant in the Padana Plain, Italy. From 1994, active pump and treat remediation has been employed. The site is bordered by canals which, in combination with variable pumping rates and groundwater flow directions, control groundwater levels. In this study we sought to determine the fate of styrene at the site within a mixed styrene/BTEX plume where the hydraulic boundaries induced strong seasonal variations in flows. In order to determine the fate of styrene, detailed field investigations provided intensive depth profile information. This information was then incorporated into a staged flow and reactive transport modelling. Three sets of measurements were obtained from sampling multilevel samplers (MLSs) under different hydraulic conditions at the site. These included measurements of BTEX, styrene, all major ions, pH and redox potential. A three-dimensional transient flow model was developed and calibrated to simulate an unconfined sandy aquifer with a variable flow field. Subsequently a reactive, multi-component transport model was employed to simulate the fate of dissolved BTEX and styrene along a selected flow line at the site. Each petroleum hydrocarbon compound was transported as independent species. Different, kinetically controlled degradation rates and a toxicity effect were simulated to explain the observed, selective degradation of pollutants in groundwater. Calibration of the model was accomplished by comparison with the three different sets of measurements obtained from the MLS devices. The results from various scenarios show that the detailed simulation of geochemical changes can be very useful to improve the site's conceptual model.