Contaminant plume classification system based on mass discharge

Quantification of contaminant mass discharge from point sources in aquitard/aquifer systems based on vertical concentration profiles and 3D modeling

Point sources with contaminants, such as chlorinated solvents, per- and polyfluoroalkyl substances (PFAS), or pesticides, are often located in low-permeability aquitards, where they can act as long-term sources and threaten underlying groundwater resources. We demonstrate the use of a 3D numerical model integrating comprehensive hydrogeological and contamination data to determine the contaminant mass discharge (CMD) from an aquitard into the underlying aquifer. A mature point source with a dissolved chlorinated solvent in a clayey till is used as an example. The quantitative determination is facilitated by model calibration to high-resolution vertical concentration profiles obtained by direct-push sampling techniques in the aquifer downgradient of the contaminant source zone. The concentration profiles showed a plume sinking with distance from the source characteristic for such aquitard/aquifer settings. The sinking is caused by the interplay between infiltrating water and horizontal groundwater flow. The application of 3D solute transport modeling on high-resolution profiles allowed for determining the infiltration rate, the hydraulic conductivity in the aquitard, and, ultimately, the CMD. Different source zone conceptualizations demonstrate the potential effects of fractures and sorption in source zones in aquitards on CMD development. Fractures in the aquitard had a minor influence on the current CMD determined with the presented approach. Still, fractures with hydraulic apertures larger than 10 μm were crucial for the temporal development of the CMD and plume. A thorough characterization of the source zone conditions combined with high-resolution concentration profiles and detailed modeling is valuable for shedding light on the probable future development of groundwater contamination arising from sources in aquitard/aquifer settings and evaluating remedial actions.

A novel concept for estimating the contaminant mass discharge of chlorinated ethenes emanating from clay till sites

Interest in using contaminant mass discharge (CMD) for risk assessment of contaminated sites has increased over the years, as it accounts for the contaminant mass that is moving and posing a risk to water resources and receptors. The most common investigation of CMD involves a transect of multilevel wells; however, this is an expensive undertaking, and it is difficult to place it in the right position in a plume. Additionally, infrastructure at the site needs to be considered. To derive an initial CMD estimate at a contaminated site and to allow for the prioritization of further investigations and remedial actions, the ProfileFlux method has been developed. It is targeted at former industrial sites with a source zone in a low conductivity layer with primarily vertical flow overlying an aquifer with primarily horizontal groundwater flow. The ProfileFlux method was developed for mature chlorinated solvent plumes, typically originating from more than 30 to 50-year-old spills, as the usage of chlorinated solvents is primarily historical. Thus, it is assumed that the contaminant had time to distribute in the low conductivity layer by mainly diffusive processes. Today the contamination is continuously released to the underlying aquifer, where advection and dispersive (other than diffusive) processes are of higher importance. The approach combines high-resolution, depth-discrete vertical concentration profiles and a simple 2D flow and transport model to estimate CMD by comparing measured and simulated concentration profiles. The study presented herein includes a global sensitivity analysis, in order to identify crucial field parameters, and of particular importance in this regard are source length, groundwater flux and infiltration. The ProfileFlux method was tested at a well-examined industrial site primarily contaminated with trichloroethylene, thereby allowing a comparison between CMD from the ProfileFlux method and the traditional transect method. CMD was estimated at 117-170 g/year, when using the ProfileFlux method, against 143 g/year with the transect method, thus validating ProfileFlux method's ability to estimate CMD. In addition, applying the method identified weak points in the conceptual site model. The method will be incorporated into a user-friendly online tool directed at environmental consultants and decision-makers working on the risk assessment and prioritization of contaminated sites with the specific hydrogeological conditions of an aquifer with an overlying low permeability layer.

Experimental study of thermally enhanced recovery of high-viscosity DNAPL in saturated porous media under non-isothermal conditions

Thermal enhancement is known to be an efficient way to decrease the residual saturation of some common dense non-aqueous phase liquids (DNAPLs) after pumping. However, the effect of transient heat transfer during the recovery of a high-viscosity contaminant, such as coal tar, in highly permeable porous media is still unknown. A 2D tank experimental setup allowing monitoring of temperature and saturation fields during DNAPL pumping has been developed. Experiments were run under isothermal and non-isothermal conditions, at low and high flow rates. We investigated the presence of viscous fingering and how that influences the shape of the cone of depression, as well as the residual saturation. The saturation fields show that less viscous fingering occurs in pre-heated cases and that heating increases the recovery efficiency. Increasing the temperature increases the critical velocity and the viscosity ratio and helps to stabilize the interface between the non-wetting and wetting phase. Observations were first made on an oil and ethanol fluid pair because its properties were known, before extending the experiments to a coal tar and water fluid pair. Residual oil saturation after pumping was decreased by 6-16% in all pre-heated conditions. Pumping at low flow rate in these conditions leaves the smallest oil residual saturation (20%) after pumping. A low flow rate increases the recovery efficiency by reducing viscous fingering and by spreading the generated heat to a larger part of the tank. Finally, results on coal tar pumping show that the high thermal conductivity of water helps in keeping the temperature high during pumping. The residual coal tar saturation was reduced from 40% at 20 °C to 28% when pre-heating the tank. Operating at a low flow rate and with a uniform temperature is the key to recovering the highest amount of a viscous DNAPL such as coal tar from the soil and satisfying cleanup goals when using thermally enhanced pumping.

Actual and Predictive Transport Modeling of Fluoride Contamination of the Sfax-Agareb Coastal Aquifer in the Mediterranean Basin

The Tunisian coast is suffering from several active or abandoned polluted sites, suspected to have released high concentrations of various contaminants infiltrating the environment and probably causing groundwater degradation. Within this scope, this study comes to assess and model the Sfax-Agareb aquifer contamination by fluoride released through phosphogypsum leachate percolation. For that, a spatial-temporal monitoring of fluoride contents was achieved for the period between October 2013 and October 2014. Observed data show that hazardous water contamination is occurring especially close to the phosphogypsum leachate collection basins. At this level, groundwater fluoride concentrations may reach up to 29 mg/L. Flow and transport modeling to evaluate fluoride contamination plume expansion in 2030 was achieved using MODFLOW and MT3DMS software packages based on a homogeneous and isotropic aquifer conceptual model. Flow and transport model calibrations were assessed by varying hydraulic conductivity, effective porosity, and dispersivity and then validated through observed data for two reference dates (October 2013 and October 2014). Based on the Tunisian NT 106-002 liquid discharge norms, fluoride contamination front was set at 3 mg/L. Numerical simulation shows that, in 2014, plume contamination by fluoride in the saturated zone extended over 250 m from the contamination source. In 2030, the spatial extent of this contamination will extend over a distance of 900 m towards the sea, following the aquifer flow direction. At the control piezometer, fluoride concentrations will increase from 29 mg/L in 2014 to 86 mg/L in 2030. This study, using numerical modeling, gives new insights for short- and medium-term prediction of eventual fluoride concentrations in the saturated zone of the Sfax-Agareb aquifer.

Source strength functions from long-term monitoring data and spatially distributed mass discharge measurements

Source strength functions (SSF), defined as contaminant mass discharge or flux-averaged concentration from dense nonaqueous phase liquid (DNAPL) source zones as a function of time, provide a quantitative model of DNAPL source-zone behavior. Such information is useful for calibration of screening-level models to assist with site management decisions. We investigate the use of historic data collected during long-term monitoring (LTM) activities at a site in Rhode Island to predict the SSF based on temporal mass discharge measurements at a fixed location, as well as SSF estimation using mass discharge measurements at a fixed time from three spatially distributed control planes. Mass discharge based on LTM data decreased from ~300 g/day in 1996 to ~70 g/day in 2012 at a control plane downgradient of the suspected DNAPL source zone, and indicates an overall decline of ~80% in 16 years. These measurements were compared to current mass discharge measurements across three spatially distributed control planes. Results indicate that mass discharge increased in the downgradient direction, and was ~6 g/day, ~37 g/day, and ~400 g/day at near, intermediate, and far distances from the suspected source zone, respectively. This behavior was expected given the decreasing trend observed in the LTM data at a fixed location. These two data sets were compared using travel time as a means to plot the data sets on a common axis. The similarity between the two data sets gives greater confidence to the use of this combined data set for site-specific SSF estimation relative to either the sole use of LTM or spatially distributed data sets.

Assessing the chemical contamination dynamics in a mixed land use stream system

Traditionally, the monitoring of streams for chemical and ecological status has been limited to surface water concentrations, where the dominant focus has been on general water quality and the risk for eutrophication. Mixed land use stream systems, comprising urban areas and agricultural production, are challenging to assess with multiple chemical stressors impacting stream corridors. New approaches are urgently needed for identifying relevant sources, pathways and potential impacts for implementation of suitable source management and remedial measures. We developed a method for risk assessing chemical stressors in these systems and applied the approach to a 16-km groundwater-fed stream corridor (Grindsted, Denmark). Three methods were combined: (i) in-stream contaminant mass discharge for source quantification, (ii) Toxic Units and (iii) environmental standards. An evaluation of the chemical quality of all three stream compartments - stream water, hyporheic zone, streambed sediment - made it possible to link chemical stressors to their respective sources and obtain new knowledge about source composition and origin. Moreover, toxic unit estimation and comparison to environmental standards revealed the stream water quality was substantially impaired by both geogenic and diffuse anthropogenic sources of metals along the entire corridor, while the streambed was less impacted. Quantification of the contaminant mass discharge originating from a former pharmaceutical factory revealed that several 100 kgs of chlorinated ethenes and pharmaceutical compounds discharge into the stream every year. The strongly reduced redox conditions in the plume result in high concentrations of dissolved iron and additionally release arsenic, generating the complex contaminant mixture found in the narrow discharge zone. The fingerprint of the plume was observed in the stream several km downgradient, while nutrients, inorganics and pesticides played a minor role for the stream health. The results emphasize that future investigations should include multiple compounds and stream compartments, and highlight the need for holistic approaches when risk assessing these dynamic systems.

Screening-level estimates of mass discharge uncertainty from point measurement methods

The uncertainty of mass discharge measurements associated with point-scale measurement techniques was investigated by deriving analytical solutions for the mass discharge coefficient of variation for two simplified, conceptual models. In the first case, a depth-averaged domain was assumed, consisting of one-dimensional groundwater flow perpendicular to a one-dimensional control plane of uniformly spaced sampling points. The contaminant flux along the control plane was assumed to be normally distributed. The second case consisted of one-dimensional groundwater flow perpendicular to a two-dimensional control plane of uniformly spaced sampling points. The contaminant flux in this case was assumed to be distributed according to a bivariate normal distribution. The center point for the flux distributions in both cases was allowed to vary in the domain of the control plane as a uniform random variable. Simplified equations for the uncertainty were investigated to facilitate screening-level evaluations of uncertainty as a function of sampling network design. Results were used to express uncertainty as a function of the length of the control plane and number of wells, or alternatively as a function of the sample spacing. Uncertainty was also expressed as a function of a new dimensionless parameter, Ω, defined as the ratio of the maximum local flux to the product of mass discharge and sample density. Expressing uncertainty as a function of Ω provided a convenient means to demonstrate the relationship between uncertainty, the magnitude of a local hot spot, magnitude of mass discharge, distribution of the contaminant across the control plane, and the sampling density.

Membrane interface probe protocol for contaminants in low-permeability zones

Accurate characterization of contaminant mass in zones of low hydraulic conductivity (low k) is essential for site management because this difficult-to-treat mass can be a long-term secondary source. This study developed a protocol for the membrane interface probe (MIP) as a low-cost, rapid data-acquisition tool for qualitatively evaluating the location and relative distribution of mass in low-k zones. MIP operating parameters were varied systematically at high and low concentration locations at a contaminated site to evaluate the impact of the parameters on data quality relative to a detailed adjacent profile of soil concentrations. Evaluation of the relative location of maximum concentrations and the shape of the MIP vs. soil profiles led to a standard operating procedure (SOP) for the MIP to delineate contamination in low-k zones. This includes recommendations for: (1) preferred detector (ECD for low concentration zones, PID or ECD for higher concentration zones); (2) combining downlogged and uplogged data to reduce carryover; and (3) higher carrier gas flow rate in high concentration zones. Linear regression indicated scatter in all MIP-to-soil comparisons, including R(2) values using the SOP of 0.32 in the low concentration boring and 0.49 in the high concentration boring. In contrast, a control dataset with soil-to-soil correlations from borings 1-m apart exhibited an R(2) of ≥ 0.88, highlighting the uncertainty in predicting soil concentrations using MIP data. This study demonstrates that the MIP provides lower-precision contaminant distribution and heterogeneity data compared to more intensive high-resolution characterization methods. This is consistent with its use as a complementary screening tool.

The uncertainty of mass discharge measurements using pumping methods under simplified conditions

Mass discharge measurements at contaminated sites have been used to assist with site management decisions, and can be divided into two broad categories: point-scale measurement techniques and pumping methods. Pumping methods can be sub-divided based on the pumping procedures used into sequential, concurrent, and tandem circulating well categories. Recent work has investigated the uncertainty of point measurement methods, and to a lesser extent, pumping methods. However, the focus of this study was a direct comparison of uncertainty between the various pumping method approaches that have been used, as well as a comparison of uncertainty between pumping and point measurement methods. Mass discharge measurement error was investigated using a Monte Carlo modeling analysis as a function of the contaminant plume position and width, and as a function of the pumping conditions used in the different pumping tests. Results indicated that for the conditions investigated, uncertainty in mass discharge estimates based on pumping methods was 1.3 to 16 times less than point measurement method uncertainty, and that a sequential pumping approach resulted in 5 to 12 times less uncertainty than the concurrent pumping or tandem circulating well approaches. Uncertainty was also investigated as a function of the plume width relative to well spacing. For a given well spacing, uncertainty decreased for all methods as the plume width increased, and comparable levels of uncertainty between point measurement and pumping methods were obtained when three wells were distributed across the plume. A hybrid pumping technique in which alternate wells were pumped concurrently in two separate campaigns yielded similar uncertainty to the sequential pumping approach. This suggests that the hybrid approach can be used to capitalize on the advantages of sequential pumping yet minimize the overall test duration.

Field demonstration and evaluation of the Passive Flux Meter on a CAH groundwater plume

This study comprises the first application of the Passive Flux Meter (PFM) for the measurement of chlorinated aliphatic hydrocarbon (CAH) mass fluxes and Darcy water fluxes in groundwater at a European field site. The PFM was originally developed and applied to measurements near source zones. The focus of the PFM is extended from near source to plume zones. For this purpose, 48 PFMs of 1.4 m length were constructed and installed in eight different monitoring wells in the source and plume zone of a CAH-contaminated field site located in France. The PFMs were retrieved, sampled, and analyzed after 3 to 11 weeks of exposure time, depending on the expected contaminant flux. PFM evaluation criteria include analytical, technical, and practical aspects as well as conditions and applicability. PFM flux data were compared with so-called traditional soil and groundwater concentration data obtained using active sampling methods. The PFMs deliver reasonable results for source as well as plume zones. The limiting factor in the PFM applicability is the exposure time together with the groundwater flux. Measured groundwater velocities at the field site range from 2 to 41 cm/day. Measured contaminant flux data raise up to 13 g/m(2)/day for perchloroethylene in the plume zone. Calculated PFM flux averaged concentration data and traditional concentration data were of similar magnitude for most wells. However, both datasets need to be compared with reservation because of the different sampling nature and time. Two important issues are the PFM tracer loss during installation/extraction and the deviation of the groundwater flow field when passing the monitoring well and PFM. The demonstration of the PFM at a CAH-contaminated field site in Europe confirmed the efficiency of the flux measurement technique for source as well as plume zones. The PFM can be applied without concerns in monitoring wells with European standards. The acquired flux data are of great value for the purpose of site characterization and mass discharge modeling, and can be used in combination with traditional soil and groundwater sampling methods.

The impact of low-temperature seasonal aquifer thermal energy storage (SATES) systems on chlorinated solvent contaminated groundwater: modeling of spreading and degradation

Groundwater systems are increasingly used for seasonal aquifer thermal energy storage (SATES) for periodic heating and cooling of buildings. Its use is hampered in contaminated aquifers because of the potential environmental risks associated with the spreading of contaminated groundwater, but positive side effects, such as enhanced contaminant remediation, might also occur. A first reactive transport study is presented to assess the effect of SATES on the fate of chlorinated solvents by means of scenario modeling, with emphasis on the effects of transient SATES pumping and applicable kinetic degradation regime. Temperature effects on physical, chemical, and biological reactions were excluded as calculations and initial simulations showed that the small temperature range commonly involved (ΔT<15 °C) only caused minor effects. The results show that a significant decrease of the contaminant mass and (eventually) plume volume occurs when degradation is described as sediment-limited with a constant rate in space and time, provided that dense non-aqueous phase liquid (DNAPL) is absent. However, in the presence of DNAPL dissolution, particularly when the dissolved contaminant reaches SATES wells, a considerably larger contaminant plume is created, depending on the balance between DNAPL dissolution and mass removal by degradation. Under conditions where degradation is contaminant-limited and degradation rates depend on contaminant concentrations in the aquifer, a SATES system does not result in enhanced remediation of a contaminant plume. Although field data are lacking and existing regulatory constraints do not yet permit the application of SATES in contaminated aquifers, reactive transport modeling provides a means of assessing the risks of SATES application in contaminated aquifers. The results from this study are considered to be a first step in identifying the subsurface conditions under which SATES can be applied in a safe or even beneficial manner.

Mass discharge in a tracer plume: evaluation of the Theissen Polygon Method

A tracer plume was created within a thin aquifer by injection for 299 d of two adjacent "sub-plumes" to represent one type of plume heterogeneity encountered in practice. The plume was monitored by snapshot sampling of transects of fully screened wells. The mass injection rate and total mass injected were known. Using all wells in each transect (0.77 m well spacing, 1.4 points/m(2) sampling density), the Theissen Polygon Method (TPM) yielded apparently accurate mass discharge (M(d) ) estimates at three transects for 12 snapshots. When applied to hypothetical sparser transects using subsets of the wells with average spacing and sampling density from 1.55 to 5.39 m and 0.70 to 0.20 points/m(2) , respectively, the TPM accuracy depended on well spacing and location of the wells in the hypothesized transect with respect to the sub-plumes. Potential error was relatively low when the well spacing was less than the widths of the sub-plumes (>0.35 points/m(2) ). Potential error increased for well spacing similar to or greater than the sub-plume widths, or when less than 1% of the plume area was sampled. For low density sampling of laterally heterogeneous plumes, small changes in groundwater flow direction can lead to wide fluctuations in M(d) estimates by the TPM. However, sampling conducted when flow is known or likely to be in a preferred direction can potentially allow more useful comparisons of M(d) over multiyear time frames, such as required for performance evaluation of natural attenuation or engineered remediation systems.

Groundwater remediation: the next 30 years

Groundwater remediation technologies are designed, installed, and operated based on the conceptual models of contaminant hydrogeology that are accepted at that time. However, conceptual models of remediation can change as new research, new technologies, and new performance data become available. Over the past few years, results from multiple-site remediation performance studies have shown that achieving drinking water standards (i.e., Maximum Contaminant Levels, MCLs) at contaminated groundwater sites is very difficult. Recent groundwater research has shown that the process of matrix diffusion is one key constraint. New developments, such as mass discharge, orders of magnitude (OoMs), and SMART objectives are now being discussed more frequently by the groundwater remediation community. In this paper, the authors provide their perspectives on the existing "reach MCLs" approach that has historically guided groundwater remediation projects, and advocate a new approach built around the concepts of OoMs and mass discharge.

Relative contribution of DNAPL dissolution and matrix diffusion to the long-term persistence of chlorinated solvent source zones

The relative contribution of dense non-aqueous phase liquid (DNAPL) dissolution versus matrix diffusion processes to the longevity of chlorinated source zones was investigated. Matrix diffusion is being increasingly recognized as an important non-DNAPL component of source behavior over time, and understanding the persistence of contaminants that have diffused into lower permeability units can impact remedial decision-making. In this study, a hypothetical DNAPL source zone architecture consisting of several different sized pools and fingers originally developed by Anderson et al. (1992) was adapted to include defined low permeability layers. A coupled dissolution-diffusion model was developed to allow diffusion into these layers while in contact with DNAPL, followed by diffusion out of these same layers after complete DNAPL dissolution. This exercise was performed for releases of equivalent masses (675 kg) of three different compounds, including chlorinated solvents with solubilities ranging from low (tetrachloroethene (PCE)), moderate (trichloroethene (TCE)) to high (dichloromethane (DCM)). The results of this simple modeling exercise demonstrate that matrix diffusion can be a critical component of source zone longevity and may represent a longer-term contributor to source longevity (i.e., longer time maintaining concentrations above MCLs) than DNAPL dissolution alone at many sites. For the hypothetical TCE release, the simulation indicated that dissolution of DNAPL would take approximately 38 years, while the back diffusion from low permeability zones could maintain the source for an additional 83 years. This effect was even more dramatic for the higher solubility DCM (97% of longevity due to matrix diffusion), while the lower solubility PCE showed a more equal contribution from DNAPL dissolution vs. matrix diffusion. Several methods were used to describe the resulting source attenuation curves, including a first-order decay model which showed that half-life of mass discharge from the matrix-diffusion dominated phase is in the range of 13 to 29 years for TCE. Because the mass discharge rate shifts significantly over time once DNAPL dissolution is complete, a Power-Law model was shown to be useful, especially at later stages when matrix diffusion dominates. An assessment of mass distribution showed that while relatively small percentages of the initial source mass diffused into the low permeability compartment, this mass was sufficient to sustain concentrations above drinking water standards for decades. These data show that relatively typical conditions (e.g., 50-year-old release, moderate to high solubility contaminant) are consistent with late stage sources, where mass in low permeability matrices serves as the primary source, and fit the conceptual model that mass in low permeability zones is important when evaluating source longevity.