Characterization of subsoil heterogeneity, estimation of grain size distribution and hydraulic conductivity at the Krauthausen test site using Cone Penetration Test

A Cone Penetration Test (CPT) survey with a high spatial resolution was performed in order to investigate the stratigraphy as well as the spatial variability of various soil properties of the Krauthausen test site. Analyses of the CPT measurements showed the subsurface to be dominated by a planar layered structure. Variogram analysis of the various CPT parameters disclosed that within each layer the soil properties have an anisotropic spatial correlation structure. A correlation analysis of the measured CPT data and co-located grain size distributions from soil samples was performed. Since the correlation coefficients were greater equal to 0.7, a reliable empirical relationship between the data sets could be developed. Based on this empirical relationship grain size distributions were estimated at CPT locations. The statistical processing of estimated and measured grain size distributions with respect to their spatial correlation structure disclosed good agreement between the data sets. The estimated grain size distributions from CPT data were used to estimate the hydraulic conductivity in the aquifer. The results provide detailed information of the spatial heterogeneity of the hydraulic conductivity at Krauthausen test site. The validation of these results, using a prior investigation of hydraulic conductivity statistics, suggests the CPT a fast and inexpensive tool for the estimation of three dimensional hydraulic conductivity fields with sufficient accuracy.     

Stochastic evaluation of mass discharge from pointlike concentration measurements

The contaminant mass discharge crossing a control plane is an important metric in the assessment of natural attenuation at contaminated sites. For risk-assessment purposes, the mass discharge must be estimated together with a level of uncertainty. We present a conditional Monte Carlo approach that allows estimating the statistical distribution of mass discharge. The approach is based on conditioning multiple realizations of the hydraulic conductivity field on all data available. We jointly determine a first-order decay coefficient in each realization, leading to conditional statistical distributions of all estimated parameters and the total mass discharge. The resulting statistical distribution of contaminant mass discharges can be used in the assessment of risks at the contaminated site. The method is applied to data of hypothetical test cases, which gives the opportunity to compare estimation results to the true field. As concentration data, we account for pointlike measurements obtained in multi-level sampling wells. The obtained empirical distribution of mass discharge crossing the multi-level sampling fence could be well fitted by a log-normal distribution.     

Combined use of field and laboratory testing to predict preferred flow paths in an heterogeneous aquifer

Elevated nitrate concentrations within a municipal water supply aquifer led to pilot testing of a field-scale, in situ denitrification technology based on carbon substrate injections. In advance of the pilot test, detailed characterization of the site was undertaken. The aquifer consisted of complex, discontinuous and interstratified silt, sand and gravel units, similar to other well studied aquifers of glaciofluvial origin, 15-40 m deep. Laboratory and field tests, including a conservative tracer test, a pumping test, a borehole flowmeter test, grain-size analysis of drill cuttings and core material, and permeameter testing performed on core samples, were performed on the most productive depth range (27-40 m), and the results were compared. The velocity profiles derived from the tracer tests served as the basis for comparison with other methods. The spatial variation in K, based on grain-size analysis, using the Hazen method, were poorly correlated with the breakthrough data. Trends in relative hydraulic conductivity (K/K(avg)) from permeameter testing compared somewhat better. However, the trends in transient drawdown with depth, measured in multilevel sampling points, corresponded particularly well with those of solute mass flux. Estimates of absolute K, based on standard pumping test analysis of the multilevel drawdown data, were inversely correlated with the tracer test data. The inverse nature of the correlation was attributed to assumptions in the transient drawdown packages that were inconsistent with the variable diffusivities encountered at the scale of the measurements. Collectively, the data showed that despite a relatively low variability in K within the aquifer under study (within a factor of 3), water and solute mass fluxes were concentrated in discrete intervals that could be targeted for later bioremediation.     

Estimating contaminant mass discharge: a field comparison of the multilevel point measurement and the integral pumping investigation approaches and their uncertainties

In this field study, two approaches to assess contaminant mass discharge were compared: the sampling of multilevel wells (MLS) and the integral groundwater investigation (or integral pumping test, IPT) that makes use of the concentration-time series obtained from pumping wells. The MLS approached used concentrations, hydraulic conductivity and gradient rather than direct chemical flux measurements, while the IPT made use of a simplified analytical inversion. The two approaches were applied at a control plane located approximately 40m downgradient of a gasoline source at Canadian Forces Base Borden, Ontario, Canada. The methods yielded similar estimates of the mass discharging across the control plane. The sources of uncertainties in the mass discharge in each approach were evaluated, including the uncertainties inherent in the underlying assumptions and procedures. The maximum uncertainty of the MLS method was about 67%, and about 28% for the IPT method in this specific field situation. For the MLS method, the largest relative uncertainty (62%) was attributed to the limited sampling density (0.63 points/m(2)), through a novel comparison with a denser sampling grid nearby. A five-fold increase of the sampling grid density would have been required to reduce the overall relative uncertainty for the MLS method to about the same level as that for the IPT method. Uncertainty in the complete coverage of the control plane provided the largest relative uncertainty (37%) in the IPT method. While MLS or IPT methods to assess contaminant mass discharge are attractive assessment tools, the large relative uncertainty in either method found for this reasonable well monitored and simple aquifer suggests that results in more complex plumes in more heterogeneous aquifers should be viewed with caution.     

A pragmatic approach for estimation of source-zone emissions at LNAPL contaminated sites

When considering natural attenuation as a remediation strategy at a site contaminated by a light non-aqueous phase liquid (LNAPL), it is important to consider the emission of contaminants from the source zone. A quantification of source-zone emissions is essential both for comparison with down-gradient mass fluxes to provide an estimate of fractional mass flux reduction, as well as for estimating the source lifetime. Because the spatial distribution of LNAPL at a field site is strongly dependent on both the spill circumstances and the heterogeneity of the geologic materials, which can be problematic for in-situ determination, alternative methods for estimating source-zone emissions are needed. In this work, a three-dimensional multiphase flow and transport modelling approach is used to investigate the relationship between the lateral extent of an LNAPL body and the emission of contaminants to groundwater at a contaminated site. For simulations involving an LNAPL release in an aquifer comprised of heterogeneous porosity and permeability distributions that were generated geostatistically, it is shown that a simple linear relationship exists between the lateral extent of the LNAPL body in the capillary fringe and the emission to the aqueous phase. The parameters describing the relationship are found to be linear functions of the groundwater flow velocity and the vertical infiltration rate. This site-specific relationship provides a simple method to estimate contaminant emissions to groundwater at LNAPL contaminated sites.     

Optimal reduction of the ozone monitoring network over France

Ozone is a harmful air pollutant at ground level, and its concentrations are measured with routine monitoring networks. Due to the heterogeneous nature of ozone fields, the spatial distribution of the ozone concentration measurements is very important. Therefore, the evaluation of distributed monitoring networks is of both theoretical and practical interests. In this study, we assess the efficiency of the ozone monitoring network over France (BDQA) by investigating a network reduction problem. We examine how well a subset of the BDQA network can represent the full network. The performance of a subnetwork is taken to be the root mean square error (rmse) of the hourly ozone mean concentration estimations over the whole network given the observations from that subnetwork. Spatial interpolations are conducted for the ozone estimation taking into account the spatial correlations. Several interpolation methods, namely ordinary kriging, simple kriging, kriging about the means, and consistent kriging about the means, are compared for a reliable estimation. Exponential models are employed for the spatial correlations. It is found that the statistical information about the means improves significantly the kriging results, and that it is necessary to consider the correlation model to be hourly-varying and daily stationary. The network reduction problem is solved using a simulated annealing algorithm. Significant improvements can be obtained through these optimizations. For instance, removing optimally half the stations leads to an estimation error of the order of the standard observational error (10 μg m^?3). The resulting optimal subnetworks are dense in urban agglomerations around Paris (Île-de-France) and Nice (Côte d’Azur), where high ozone concentrations and strong heterogeneity are observed. The optimal subnetworks are probably dense near frontiers because beyond these frontiers there is no observation to reduce the uncertainty of the ozone field. For large rural regions, the stations are uniformly distributed. The fractions between urban, suburban and rural stations are rather constant for optimal subnetworks of larger size (beyond 100 stations). By contrast, for smaller subnetworks, the urban stations dominate.     

Small-scale variability of metals in soil and composite sampling

Soil pollution data is also strongly scattering at small scale. Sampling of composite samples, therefore, is recommended for pollution assessment. Different statistical methods are available to provide information about the accuracy of the sampling process. Autocorrelation and variogram analysis can be applied to investigate spatial relationships. Analysis of variance is a useful method for homogeneity testing. The main source of the total measurement uncertainty is the uncertainty arising from sampling. The sample mass required for analysis can also be estimated using an analysis of variance. The number of increments to be taken for a composite sample can be estimated by means of simple statistical formulae. Analytical results of composite samples obtained from different fusion procedures of increments can be compared by means of multiple mean comparison. The applicability of statistical methods and their advantages are demonstrated for a case study investigating metals in soil at a very small spatial scale. The paper describes important statistical tools for the quantitative assessment of the sampling process. Detailed results clearly depend on the purpose of sampling, the spatial scale of the object under investigation and the specific case study, and have to be determined for each particular case.     

New methodology to investigate potential contaminant mass fluxes at the stream-aquifer interface by combining integral pumping tests and streambed temperatures

The spatial pattern and magnitude of mass fluxes at the stream-aquifer interface have important implications for the fate and transport of contaminants in river basins. Integral pumping tests were performed to quantify average concentrations of chlorinated benzenes in an unconfined aquifer partially penetrated by a stream. Four pumping wells were operated simultaneously for a time period of 5 days and sampled for contaminant concentrations. Streambed temperatures were mapped at multiple depths along a 60m long stream reach to identify the spatial patterns of groundwater discharge and to quantify water fluxes at the stream-aquifer interface. The combined interpretation of the results showed average potential contaminant mass fluxes from the aquifer to the stream of 272microgm(-2)d(-1) MCB and 71microgm(-2)d(-1) DCB, respectively. This methodology combines a large-scale assessment of aquifer contamination with a high-resolution survey of groundwater discharge zones to estimate contaminant mass fluxes between aquifer and stream.     

Estimating spatially-variable first-order rate constants in groundwater reactive transport systems

Numerical reactive transport models are often used as tools to assess aquifers contaminated with reactive groundwater solutes as well as investigating mitigation scenarios. The ability to accurately simulate the fate and transport of solutes, however, is often impeded by a lack of information regarding the parameters that define chemical reactions. In this study, we employ a steady-state Ensemble Kalman Filter (EnKF), a data assimilation algorithm, to provide improved estimates of a spatially-variable first-order rate constant λ through assimilation of solute concentration measurement data into reactive transport simulation results. The methodology is applied in a steady-state, synthetic aquifer system in which a contaminant is leached to the saturated zone and undergoes first-order decay. Multiple sources of uncertainty are investigated, including hydraulic conductivity of the aquifer and the statistical parameters that define the spatial structure of the parameter field. For the latter scenario, an iterative method is employed to identify the statistical mean of λ of the reference system. Results from all simulations show that the filter scheme is successful in conditioning the λ ensemble to the reference λ field. Sensitivity analyses demonstrate that the estimation of the λ values is dependent on the number of concentration measurements assimilated, the locations from which the measurement data are collected, the error assigned to the measurement values, and the correlation length of the λ fields.     

Towards optimal sampling schedules for integral pumping tests

Ethanol use as a gasoline additive is increasing, as are the chances of groundwater contamination caused by gasoline releases involving ethanol. To evaluate the impact of ethanol on dissolved hydrocarbon plumes, a field test was performed in which three gasoline residual sources with different ethanol fractions (E0: no ethanol, E10: 10% ethanol and E95: 95% ethanol) were emplaced below the water table. Using the numerical model BIONAPL/3D, the mass discharge rates of benzene, toluene, ethylbenzene, xylenes, trimethylbenzenes and naphthalene were simulated and results compared to those obtained from sampling transects of multilevel samplers. It was shown that ethanol dissolved rapidly and migrated downgradient as a short slug. Mass discharge of the hydrocarbons from the E0 and E10 sources suggested similar first-order hydrocarbon decay rates, indicating that ethanol from E10 had no impact on hydrocarbon degradation. In contrast, the estimated hydrocarbon decay rates were significantly lower when the source was E95. For the E0 and E10 cases, the aquifer did not have enough oxygen to support complete mineralization of the hydrocarbon compounds to the extent suggested by the field-based mass discharge. Introducing a heterogeneous distribution of hydraulic conductivity did little to overcome this discrepancy. A better match between the numerical model and the field data was obtained assuming partial degradation of the hydrocarbons to intermediate compounds. Besides depending on the ethanol concentration, the impact of ethanol on hydrocarbon degradation appears to be highly dependent on the availability of electron acceptors.     

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland-Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

The significance of heterogeneity on mass flux from DNAPL source zones: an experimental investigation

Understanding the process of mass transfer from source zones of aquifers contaminated with organic chemicals in the form of dense non-aqueous phase liquids (DNAPL) is of importance in site management and remediation. A series of intermediate-scale tank experiments was conducted to examine the influence of aquifer heterogeneity on DNAPL mass transfer contributing to dissolved mass emission from source zone into groundwater under natural flow before and after remediation. A Tetrachloroethylene (PCE) spill was performed into six source zone models of increasing heterogeneity, and both the spatial distribution of the dissolution behavior and the net effluent mass flux were examined. Experimentally created initial PCE entrapment architecture resulting from the PCE migration was largely influenced by the coarser sand lenses and the PCE occupied between 30 and 60% of the model aquifer depth. The presence of DNAPL had no apparent effect on the bulk hydraulic conductivity of the porous media. Up to 71% of PCE mass in each of the tested source zone was removed during a series of surfactant flushes, with associated induced PCE mobilization responsible for increasing vertical DNAPL distributions. Effluent mass flux due to water dissolution was also found to increase progressively due to the increase in NAPL-water contact area even though the PCE mass was reduced. Doubling of local groundwater flow velocities showed negligible rate-limited effects at the scale of these experiments. Thus, mass transfer behavior was directly controlled by the morphology of DNAPL within each source zone. Effluent mass flux values were normalized by the up-gradient DNAPL distributions. For the suite of aquifer heterogeneities and all remedial stages, normalized flux values fell within a narrow band with mean of 0.39 and showed insensitivity to average source zone saturations.     

Use of tandem circulation wells to measure hydraulic conductivity without groundwater extraction

Conventional methods to measure the hydraulic conductivity of an aquifer on a relatively large scale (10-100 m) require extraction of significant quantities of groundwater. This can be expensive, and otherwise problematic, when investigating a contaminated aquifer. In this study, innovative approaches that make use of tandem circulation wells to measure hydraulic conductivity are proposed. These approaches measure conductivity on a relatively large scale, but do not require extraction of groundwater. Two basic approaches for using circulation wells to measure hydraulic conductivity are presented; one approach is based upon the dipole-flow test method, while the other approach relies on a tracer test to measure the flow of water between two recirculating wells. The approaches are tested in a relatively homogeneous and isotropic artificial aquifer, where the conductivities measured by both approaches are compared to each other and to the previously measured hydraulic conductivity of the aquifer. It was shown that both approaches have the potential to accurately measure horizontal and vertical hydraulic conductivity for a relatively large subsurface volume without the need to pump groundwater to the surface. Future work is recommended to evaluate the ability of these tandem circulation wells to accurately measure hydraulic conductivity when anisotropy and heterogeneity are greater than in the artificial aquifer used for these studies.     

A hybrid method for inverse characterization of subsurface contaminant flux

The methods presented in this work provide a potential tool for characterizing contaminant source zones in terms of mass flux. The problem was conceptualized by considering contaminant transport through a vertical "flux plane" located between a source zone and a downgradient region where contaminant concentrations were measured. The goal was to develop a robust method capable of providing a statement of the magnitude and uncertainty associated with estimated contaminant mass flux values. In order to estimate the magnitude and transverse spatial distribution of mass flux through a plane, the problem was considered in an optimization framework. Two numerical optimization techniques were applied, simulated annealing (SA) and minimum relative entropy (MRE). The capabilities of the flux plane model and the numerical solution techniques were evaluated using data from a numerically generated test problem and a nonreactive tracer experiment performed in a three-dimensional aquifer model. Results demonstrate that SA is more robust and converges more quickly than MRE. However, SA is not capable of providing an estimate of the uncertainty associated with the simulated flux values. In contrast, MRE is not as robust as SA, but once in the neighborhood of the optimal solution, it is quite effective as a tool for inferring mass flux probability density functions, expected flux values, and confidence limits. A hybrid (SA-MRE) solution technique was developed in order to take advantage of the robust solution capabilities of SA and the uncertainty estimation capabilities of MRE. The coupled technique provided probability density functions and confidence intervals that would not have been available from an independent SA algorithm and they were obtained more efficiently than if provided by an independent MRE algorithm.     

Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

Dimension reduction and source identification for multispecies groundwater contamination

Assessment of chemical contamination at large industrial complexes with long and sometimes unknown histories of operation represents a challenging environmental problem. The spatial and temporal complexity of the contaminant may be due to changes in production processes, differences in the chemical transport, and the physical heterogeneity of the soil and aquifer materials. Traditional mapping techniques are of limited value for sites where dozens of chemicals with diverse transport characteristics may be scattered over large spatial areas without documentation of disposal histories. In this context, a site with a long and largely undocumented disposal history of shallow groundwater contamination is examined using principal component analysis (PCA). The dominant chemical groups and chemical "modes" at the site were identified. PCA results indicate that five primary and three transition chemical groups can be identified in the space of the first three eigenvectors of the correlation matrix, which account for 61% of the total variance of the data. These groups represent a significant reduction in the dimension of the original data (116 chemicals). It is shown that each group represents a class of chemicals with similar chemo-dynamic properties and/or environmental response. Finally, the groups are mapped back onto the site map to infer delineation of contaminant source areas for each class of compounds. The approach serves as a preliminary step in subsurface characterization, and a data reduction strategy for source identification, subsurface modeling and remediation planning.     

Hexadecane and pristane degradation potential at the level of the aquifer—evidence from sediment incubations compared to in situ microcosms

Monitored natural attenuation is widely accepted as a sustainable remediation method. However, methods providing proof of proceeding natural attenuation within the water-unsaturated (vadose) zone are still relying on proxies such as measurements of reactive and non-reactive gases, or sediment sampling and subsequent mineralisation assays, under artificial conditions in the laboratory. In particular, at field sites contaminated with hydrophobic compounds, e.g. crude oil spills, an in situ evaluation of natural attenuation is needed, because in situ methods are assumed to provide less bias than investigations applying either proxies for biodegradation or off-site microcosm experiments. In order to compare the current toolbox of methods with the recently developed in situ microcosms, incubations with direct push-sampled sediments from the vadose and the aquifer zones of a site contaminated with crude oil were carried out in conventional microcosms and in situ microcosms. The results demonstrate the applicability of the in situ microcosm approach also outside water-saturated aquifer conditions in the vadose zone. The sediment incubation experiments demonstrated turnover rates in a similar range (vadose, 4.7 mg/kg*day; aquifer, 6.4 mg_hexadecane/kg_soil/day) of hexadecane degradation in the vadose zone and the aquifer, although mediated by slightly different microbial communities according to the analysis of fatty acid patterns and amounts. Additional experiments had the task of evaluating the degradation potential for the branched-chain alkane pristane (2,6,10,14-tetramethylpentadecane). Although this compound is regarded to be hardly degradable in comparison to n-alkanes and is thus frequently used as a reference parameter for indexing the extent of biodegradation of crude oils, it could be shown to be degraded by means of the incubation experiments. Thus, the site had a high inherent potential for natural attenuation of crude oils both in the vadose zone and the aquifer.