The fate of DNAPL contaminants in non-consolidated subsurface systems - Discussion on the relevance of effective source zone geometries for plume propagation

Mapping and monitoring dense non-aqueous phase liquid source zone by fused surface and cross-borehole electrical resistivity tomography

Effective characterization of dense non-aqueous phase liquid (DNAPL) source zones is crucial for remediating polluted sites. DNAPL often reside as residuals or pools within high-permeability lenses and above impermeable layers due to soil heterogeneity, gravity, and capillary barriers. Given the high cost of drilling, electrical resistivity tomography (ERT) techniques-including surface ERT and cross-borehole ERT, are commonly used for DNAPL source zone mapping and monitoring. However, the low spatial resolution of ERT increases uncertainty in source zone investigations. This study proposes a method for improving DNAPL mapping and monitoring by fusing surface and cross-borehole ERT data. Sandbox experiments were conducted to simulate a heterogeneous DNAPL source zone, employing both ERT methods for static mapping and dynamic monitoring. Reflective light imaging (RLM) was used to visualize DNAPL migration and provide saturation data, allowing for the quantification of ERT's effectiveness in characterizing DNAPL distribution. The results indicate that individual ERT methods face significant challenges in DNAPL source zone mapping due to background interference. Surface ERT alone tends to underestimate the extent of deeper DNAPL source zones. However, fusing surface and cross-borehole ERT results in a complementary enhancement of vertical spatial resolution, thereby improving the characterization of DNAPL source zones. The fusion of static and time-lapse ERT data substantially enhances DNAPL source zone mapping and monitoring capabilities. By calculating the ratio of the ERT-monitored area to the actual area using resistivity change contours (5 %, 10 %, 15 %), it was found that fusing surface and cross-borehole ERT data improved monitoring resolution by 50.48 % compared to surface ERT alone and by 22.95 % compared to cross-borehole ERT. Principal component analysis (PCA) was effective in fusing time-lapse data, while the weighted average method (WAM) outperformed PCA for static resistivity data fusion.

The effect of Tween 80 on monochlorobenzene migration in bentonite

Several studies conducted at industrial sites have documented the infiltration of dense non-aqueous phase liquids (DNAPLs) into clay layers, a phenomenon potentially influenced by the coexistence of chemicals like surfactants in some common pollutants. Bentonite (Ben), monochlorobenzene (MCB), and Tween 80 (T80) were selected as reference components to investigate the influences of nonionic surfactants on DNAPLs migration in clays. Results showed that T80 promotes MCB dissolution and encourages MCB adsorption on Ben. This process reduces the hydrophilicity of Ben, resulting in water loss and shrinkage, which creates cracks and facilitates the migration of MCB within the clay. Tw80 notably enhances MCB solubility, as indicated by a molar solubilization ratio of 7.80. The MCB adsorption on Ben (QMCB) displays a linear increase with raising the T80 adsorption on Ben (QT80), especially when QT80 are below the thresholds, e.g., 408.4 mg/g at pH 3 and 339.3 mg/g at pH 7; however, QMCB is decreased with increasing adsorbed T80 further. The average fracture ratio, crack length, and crack width of cracked samples in the cracking experiments were 0.794%, 11.29 mm, and 0.209 mm, respectively. The findings here contribute to understanding the role of surfactants in VOC transport in contaminated sites.

A data-driven approach for simplifying the estimation of time for contaminant plumes to reach their maximum extent

Globally there exist a very large number of contaminated or possibly contaminated sites where a basic preliminary assessment has not been completed. This is largely, among others, due to limited simple methods/models available for estimating key site quantities such as the maximum plume length, further denoted as Lmax and the corresponding time T=TLmax, at which the plume reaches its maximum extent L=Lmax. An approach to easily obtain an estimate of TLmax in particular is presented in this work. Limited availability of high-quality field data, particularly of TLmax, necessitates the use of synthetic data, which constrains the overall model development works. Taking BIOSCREEN-AT (transient 3D model) as a base model, this work proposes second-order polynomial models, with only two parameters, for estimating Lmax and TLmax. This reformulation of the well established solution significantly reduces data requirement and workload for initial site assessment purposes. A global sensitivity analysis (Morris, 1991), using a large number of random synthetic data, identifies the first-order decay rate constants in the plume λEFF and at the source γ as dominantly most influential for TLmax. For Lmax, the first-order decay rate constant λEFF and groundwater velocity v are the two important parameters. The sensitivity analysis also identifies that these parameters non-linearly impact TLmax or Lmax. With this information, the proposed polynomial models (each for Lmax and TLmax) were trained to obtain model coefficients, using a large amount of synthetic data. For verification, the developed models were tested using four datasets comprising over 100 sample sets against the results obtained from BIOSCREEN-AT and the developed BIOSCREEN-AT-based steady-state model. Additionally, the developed models were evaluated against two well documented field sites. The proposed models largely simplify estimation, particularly, of TLmax, for which only very limited field or literature information is available.

Utilizing machine learning algorithm for finely three-dimensional delineation of soil-groundwater contamination in a typical industrial park, North China: Importance of multisource auxiliary data

Intensive industrial activities cause soil contamination with wide variations and even perturb groundwater safety. Precision delineation of soil contamination is the foundation and precondition for soil quality assurance in the practical environmental management process. However, spatial non-stationarity phenomenon of soil contamination and heterogeneous sampling are two key issues that affect the accuracy of contamination delineation model. Taking a typical industrial park in North China as the research object, we constructed a random forest (RF) model for finely characterizing the distribution of soil contaminants using sparse-biased drilling data. Results showed that the R2 values of arsenic and 1,2-dichloroethane predicted by RF (0.8896 and 0.8973) were greatly higher than those of inverse distance weighted model (0.2848 and 0.2908), indicating that RF was more adaptable to actual non-stationarity sites. The back propagation neural network algorithm was utilized to establish a three-dimensional visualization of the contamination parcel of subsoil-groundwater system. Multiple sources of environmental data, including hydrogeological conditions, geochemical characteristics and anthropogenic industrial activities were integrated into the model to optimize the prediction accuracy. The feature importance analysis revealed that soil particle size was dominant for the migration of arsenic, while the migration of 1,2-dichloroethane highly depended on vertical permeability coefficients of the soil. Contaminants migrated downwards with soil water under gravity-driven conditions and penetrated through the subsoil to reach the saturated aquifer, forming a contamination plume with groundwater flow. Our findings afford a new idea for spatial analysis of soil-groundwater contamination at industrial sites, which will provide valuable technical support for maintaining sustainable industry.

Research on the detection of leakage points in vertical barrier walls using a combined method of ERT and tracer methods

In order to improve the detection accuracy of vertical barrier leakage, three contamination leakage working conditions, including point leakage only, point and vertical leakage, and horizontal leakage, were simulated by small-scale soil tank tests. The dynamic evolution of soil resistivity over time was monitored by electrical resistivity tomography (ERT). The accuracy of the ERT detection results was validated through thermal tracer method, chloride tracer method, and soil true resistivity experiments. The results indicate that the resistivity profiles at different times can more accurately reflect information on the location of leakage points, the extent of contamination plumes, and the migration pathways of pollutants under different working conditions. The extent of anomalous areas in resistivity profiles is a crucial factor in representing the geometric shape of pollution leakage. However, the preferential seepage or lateral migration of contaminant in the soil significantly reduces the detection accuracy of ERT for identifying leakage points. The thermal tracer method and the chloride tracer method can produce better complementary interpretations of ERT monitoring results. The measurement points near the leakage point exhibit faster temperature response rates, which can serve as a characteristic for identifying the location of leakage points. Compared with the thermal tracer method, the chloride tracer method can monitor the migration of contaminants over a larger range. Therefore, the proposed combined diagnostic detection method in this paper presents a feasible solution with promising engineering applications in leakage detection for vertical barrier barriers.

Hydrodynamically-enhanced transfer of dense non-aqueous phase liquids into an aqueous reservoir

The use of surfactants represents a viable strategy to boost the removal yield of Dense Non-Aqueous Phase Liquids (DNAPLs) from groundwater and to shorten the operational timing of the remediation process. Surfactants, in general, help in reducing the interfacial tension at the DNAPL/water interface and enhance the solubility of the pollutant in the water phase through the formation of dispersed systems, such as micelles and emulsions. In this paper, we show that a suitable choice of a surfactant, in this case belonging to the bio-degradable class of ethoxylated alcohols, allows for the formation of hydrodynamic interfacial instabilities that further enhances the dissolution rate of the organic pollutant into the water phase. In a stratified configuration (denser organic phase at the bottom and lighter water phase on top), the instabilities appear as upward-pointing fingers that originate from the inversion of the local density at the interface. This inversion stems from the synergetic coupling of two effects promoted by the ethoxylated surfactant: i) the enhanced co-solubility of the DNAPL into the water (and viceversa), and (ii) the differential diffusion of the DNAPL and the surfactant in the aqueous phase. By dissolving into the DNAPL, the surfactant also reduces locally the surface tension at the liquid-liquid interface, thereby inducing transversal Marangoni flows. In our work, we carefully evaluated the effects of the concentration of different surfactants (two different ethoxylated alcohols, sodium dodecylsulphate, cetyltrimethyl ammonium bromide, N-tetradecyl-N, N-dimethylamine oxide and bis(2-ethylhexyl) sulfosuccinate sodium salt) on the onset of the instabilities in 3 different DNAPLs/water stratifications, namely chloroform, trichloroethylene and tetrachloroethylene, with a special emphasis on the trichloroethylene/water system. By means of a theoretical model and nonlinear simulations, supported by surface tension, density and diffusivity measurements, we could provide a solid explanation to the observed phenomena and we found that the type of the dispersed system, the solubility of the DNAPL into the water phase, the solubility of the surfactant in the organic phase, as well as the relative diffusion and density of the surfactant and the DNAPL in the aqueous phase, are all key parameters for the onset of the instabilities. These results can be exploited in the most common remediation techniques.

An approach for quantification of the heterogeneity of DNAPL source zone geometries

Many studies have investigated the migration and entrapment processes of source zones from dense non-aqueous phase liquid (DNAPL) contamination under different conditions. However, the characterization of occupying area by source zone (or source shape) in water-saturated aquifers is still rudimentarily considered. In this study, we demonstrated this issue (1) by providing a brief review of existing approaches for source shape consideration, (2) by proposing an approach with simple shape parameters based on the non-uniformity of source widths, and (3) by providing exemplary applications of our proposed approach on shapes already published in previous research works. Our literature review suggested that the source zone in mathematical approaches is generally characterized as simple geometrical shapes (arbitrary lines or rectangles) or system-defined parameters that contrast to complex and discontinuous shapes observed in the real world. But the characterization of such complex shapes is still not possible with acceptable efforts. Therefore, we proposed an approach to parameterize the source shape by considering the variation of width and midpoints over the depth of the entire source zone and formulate four parameters based on population statistics (mean, standard deviation). To illustrate the suitability of our approach, we applied it to the results of lab experiments, and by analyzing these complex shapes, we highlighted the potential for improving the characterization techniques of non-uniformity of the source zones.

Prediction of Dichloroethene Concentration in the Groundwater of a Contaminated Site Using XGBoost and LSTM

Chlorinated aliphatic hydrocarbons (CAHs) are widely used in agriculture and industries and have become one of the most common groundwater contaminations. With the excellent performance of the deep learning method in predicting, LSTM and XGBoost were used to forecast dichloroethene (DCE) concentrations in a pesticide-contaminated site undergoing natural attenuation. The input variables included BTEX, vinyl chloride (VC), and five water quality indicators. In this study, the predictive performances of long short-term memory (LSTM) and extreme gradient boosting (XGBoost) were compared, and the influences of variables on models’ performances were evaluated. The results indicated XGBoost was more likely to capture DCE variation and was robust in high values, while the LSTM model presented better accuracy for all wells. The well with higher DCE concentrations would lower the model’s accuracy, and its influence was more evident in XGBoost than LSTM. The explanation of the SHapley Additive exPlanations (SHAP) value of each variable indicated high consistency with the rules of biodegradation in the real environment. LSTM and XGBoost could predict DCE concentrations through only using water quality variables, and LSTM performed better than XGBoost.

A feasibility study for the treatment of 1,2-dichloroethane-contaminated groundwater using reedbed system and assessment of its natural attenuation

A reedbed system planted with Phragmites australis was implemented to treat chlorinated hydrocarbon-contaminated groundwater in an industrial plant area. Reedbed commissioning was conducted from July 2016 to November 2016 to treat contaminated groundwater via a pump-and-treat mechanism. Combination of horizontal and vertical reedbed systems was applied to treat 1,2-dichloroethane (1,2 DCA) under four parallel installations. The 2-acre horizontal and vertical reedbed systems were designed to treat approximately 305 m³/day of pumped groundwater. Initial concentration of 1,2 DCA was observed at 0.362 mg/L to 4320 mg/L, and the reedbed system successfully reduced the concentration up to 67.9%. The average outlet concentration was measured to be 2.08 mg/L, which was lower than the site-specific target level of 156 mg/L. Natural attenuation analysis was conducted using first-order decay kinetics, showing an average natural attenuation rate of 0.00372/year. Natural attenuation of 1,2 DCA was observed in shallow monitoring wells, which was indicated by the reduction trend of 1,2 DCA concentration, thereby confirming that the reedbed system worked well to remove 1.2 DCA from contaminated groundwater at the shallow profile.

Aquifer flushing using a SDS/1-butanol based in-situ microemulsion: Performance and mechanism for the remediation of nitrobenzene contamination

In-situ microemulsion flushing is an effective remediation technology for the removal of dense non-aqueous phase liquids (DNAPLs) from aquifers. Nitrobenzene (NB) is a typical DNAPL pollutant that is responsible for the serious contamination of many groundwater systems, while its removal using the flushing method has rarely been studied. In this study, bench scale, 1-D column and 2-D tank experiments were conducted to establish an efficient salt-free sodium dodecyl sulfate (SDS)/1-butanol based in-situ microemulsion flushing system for NB contaminated aquifers. Results showed that the NB/SDS/1-butanol/water microemulsion increased dissolved NB concentrations by more than 15-fold compared to the SDS-only solution. The formulation also presented good solubilization capacity at low temperature (5 ℃) and with clay media. NB was effectively removed from the aquifer by solubilization and mobilization via the formation of the microemulsion with the injected SDS/1-butanol solution. The flushing system also reduced the tailing phenomenon in later remediation stages, and exhibited weak reagent adsorption onto aquifer media. Furthermore, the vertical DNAPL migration to deeper aquifer was effectively controlled. Therefore, the constructed in-situ microemulsion flushing system is a highly efficient treatment method for NB contaminated aquifers, with this study providing valuable reference information on the optimal reagent parameters and the remediation mechanism.

Upscaled modeling of complex DNAPL dissolution

A straightforward, upscaled DNAPL mass dissolution model is developed using relatively simple input consisting of characteristic dimensions and saturations of a DNAPL accumulation. Multiple accumulations are aggregated into a single source zone volume. Physically, the dissolution process is a combination of flow through the mass (advective component) and flow around the mass (dispersive component). The contribution of each component is based on initial characteristic length scales and the average initial saturation. Changes over time with the depletion of mass are captured with a changing relative permeability and a power law relationship for the fraction of initial mass remaining. The utility of the upscaled process model is demonstrated with data from three studies: numerical simulation of multiple pools, two-dimensional test cell experiments with mixed architecture and with heterogeneous soil, and a controlled field study of multicomponent DNAPL release and depletion. Use of the model successfully reproduced the observed multistage mass discharge in each study and illuminated the governing processes. The power law exponent was relatively constant for the various conditions and relative permeability changes were integral to the success. The numerical and experimental studies were run to complete mass depletion which the upscaled model matched. The input parameters are minimal and are found in typical DNAPL source zone characterization data.

Towards predicting DNAPL source zone formation to improve plume assessment: Using robust laboratory and numerical experiments to evaluate the relevance of retention curve characteristics

We conducted multiple laboratory trials in a robust and repeatable experimental layout to study dense non-aqueous phase liquid (DNAPL) source zone formation. We extended an image processing and analysis framework to derive DNAPL saturation distributions from reflective optical imaging data, with volume balance deviations < 5.07%. We used a multiphase flow model to simulate source zone formation in a Monte Carlo approach, where the parameter space was defined by the variation of retention curve parameters. Integral and geometric measures were used to characterize the source zones and implemented into a multi-criteria objective function. The latter showed good agreement between observation data and simulation results for effective DNAPL saturation values > 0.04, especially for early stages of DNAPL migration. The common hypothesis that parameters defining the DNAPL-water retention curves are constant over time was not confirmed. Once DNAPL pooling started, the optimal fit in the parameter space was significantly different compared to the earlier DNAPL migration stages. We suspect more complex processes (e.g., capillary hysteresis, adsorption) to become relevant during pool formation. Our results reveal deficits in the grayscale-DNAPL saturation relationship definition and laboratory estimation of DNAPL-water retention curve parameters to overcome current limitations to describe DNAPL source zone formation.

Assessing the impact of water infiltration on LNAPL mobilization in sand column using simplified image analysis method

Laboratory-scale column experiments were carried out to assess the influence of water infiltration on pooled light non-aqueous phase liquid (LNAPL) redistribution in porous media. A simplified image analysis method (SIAM) was used to evaluate the saturation distributions of the LNAPL and water in the entire domain under dynamic conditions. The experiments were conducted for high/low LNAPL volumes LNAPL volumes differentiated as low and high volumes. High resolution SIAM images of the soil column during LNAPL migration and water infiltration events were captured and analyzed. Results indicated that the capillary fringe is about 6-7 cm which was consistent with the capillary height derived from empirical equations. Moreover, SIAM provided an estimate of the field capacity (30%) of the sand. Once the LNAPL infiltration stage was started, the LNAPL was observed to rapidly migrate through the vadose zone. For the case of large LNAPL volume, the LNAPL penetrated further into capillary fringe zone. Analysis of SIAM images showed that the LNAPL redistribution was observed to vary significantly with the rate of infiltration. For higher water infiltration intensity, the injected water exerted a larger hydrodynamic force on the entrapped LNAPL forcing it move further downward into the capillary zone and the saturated zone. Overall, this study demonstrated that the SIAM technique is an accurate and cost-effective tool for the visualization of the time-dependent NAPL/water movement in laboratory-scale experiments and dynamic changes in fluid saturation in porous media.

Influence of recharge rates on steady-state plume lengths

A large number of potentially contaminated sites reported worldwide require cost- and time-effective assessment of the extent of contamination and the threats posed to the water resources. A significant risk assessment metric for these sites can be the determination of the maximum (i.e., steady-state) contaminant plume length (L_max). Analytical approaches in the literature provide an option for such an assessment, but they include a certain degree of uncertainty. Often, the causes of such uncertainties are the simplifications in the analytical models, e.g., not considering the influence of hydrogeological stresses such as recharge, which impact the plume development significantly. This may lead to an over- or underestimation of L_max. This work includes the influence of the recharge for the effective estimation of L_max. For that, several two-dimensional (2D) numerical simulations have been performed by considering different aquifer thicknesses (1 m- 4 m) and recharge rates (ranging from 0 to 3.6 mm/day). From the numerical results of this work, it has been deduced that 1) the application of the recharge shortens L_max, and the recharge entering the aquifer top causes the plume to tilt, 2) the reduction percentage in L_max depends on the recharge rate applied and the aquifer thickness, and 3) the reduction percentage varies in a non-linear manner with respect to the recharge rate for a fixed aquifer thickness. Based on these results, a hybrid analytical-empirical solution has been developed for the estimation of L_max with the inclusion of the recharge rate. The proposed hybrid analytical-empirical solution superimposes an empirically obtained correction factor onto an analytical solution. Although extensive confirmation steps of the developed model are required for including the effect of the recharge on aquifer hydraulics, the proposed expression improves the estimation of the L_max significantly. The hybrid analytical-empirical solution has also been confirmed with a selection of limited field contamination sites data. The hybrid model result (L_hyb) provides a significant improvement in the estimation, i.e., an order of magnitude lower mean relative error compared to the analytical model.

Layered Double Hydroxides with Intercalated Permanganate and Peroxydisulphate Anions for Oxidative Removal of Chlorinated Organic Solvents Contaminated Water

The contamination by chlorinated organic solvents is a worldwide problem as they can deeply penetrate aquifers, accumulating in the sub-surface as lenses of highly hazardous pollutants. In recent years, so called in situ oxidation processes have been developed to remediate chlorinated organic solvents from groundwater and soil by injecting solutions of oxidising agents such as permanganate or peroxydisulphate. We here present modified layered double hydroxides (LDHs) with intercalated oxidising agents that might serve as new reactants for these remediation strategies. LDHs might serve as support and stabiliser materials for selected oxidising agents during injection, as the uncontrolled reaction and consumption might be inhibited, and guarantee that the selected oxidants persist in the subsurface after injection. In this study, LDHs with hydrotalcite- and hydrocalumite-like structures intercalated with permanganate and peroxydisulphate anions were synthesised and their efficiency was tested in batch experiments using trichloroethene or 1,1,2-trichloroethane as the target contaminants. All samples were characterised using powder X-ray diffraction, thermal analysis coupled with mass spectrometry to directly analyse evolving gases, and Fourier-transform infrared spectroscopy. Additionally, particle size distribution measurements were carried out on the synthesised materials. Results of the batch experiments confirmed the hypothesis that oxidising agents keep their properties after intercalation. Permanganate intercalated LDHs proved to be most efficient at degrading trichloroethene while peroxydisulphate intercalated Ca,Al-LDHs were the most promising studied reactants degrading 1,1,2-trichloroethane. The detection of dichloroethene as well as the transformation of the studied reactants into new LDH phases confirmed the successful degradation of the target contaminant by oxidation processes generated from the intercalated oxidising agent.

'Emulsion locks' for the containment of hydrocarbons during surfactant flushing

Reversible double water in oil in water (W/O/W) emulsions were developed to contain subsurface hydrocarbon spills during their remediation using surfactant flushing. Double emulsions were prepared by emulsifying CaCl₂ solutions in canola oil, and subsequently by emulsifying the W/O emulsions in aqueous sodium alginate solutions. The formation of double emulsions was confirmed with confocal and optical microscopy. The double emulsions reversed and gelled when mixed with the surfactants sodium dodecyl sulfate (SDS) and cocamidopropyl betaine (CPB). Gels can act as 'emulsion locks' to prevent spreading of the hydrocarbon plume from the areas treated with surfactant flushing, as shown in sand column tests. Shear rheology was used to quantify the viscoelastic moduli increase (gelation) upon mixing the double emulsion with SDS and CPB. SDS was more effective than CPB in gelling the double emulsions. CPB and SDS could adsorb at the interface between water and model hydrocarbons (toluene and motor oil), lowering the interfacial tension and rigidifying the interface (as shown with a Langmuir trough). Bottle tests and optical microscopy showed that SDS and CPB produced W/O and O/W emulsions, with either toluene or motor oil and water. The emulsification of motor oil and toluene in water with SDS and CPB facilitated their flow through sand columns and their recovery. Toluene recovery from sand columns was quantitated using Gas-Chromatography Mass-Spectroscopy (GC-MS). The data show that SDS and CPB can be used both for surfactant flushing and to trigger the gelation of 'emulsion locks'. Ethanol also gelled the emulsions at 100 mL/L.

Quantification of Uncertainties from Image Processing and Analysis in Laboratory-Scale DNAPL Release Studies Evaluated by Reflective Optical Imaging

Subsurface DNAPL (dense non-aqueous phase liquid) contamination from (un-) intentional spilling typically leads to severe environmental hazards. A large number of studies have demonstrated the relevance of DNAPL source zone geometry for the determination of contaminant plume propagation in groundwater. Optical imaging represents a promising non-invasive method for identifying DNAPL saturation without disturbing multiphase flow dynamics. However, workflow and image analysis methodologies have not been sufficiently developed or described for general application to related experimental efforts. For example, the choice of dye(s) used for phase colorization affects image processing and can bias final estimations of DNAPL saturations. In this study, we perform a series of DNAPL migration and entrapment studies in transparent tanks that are filled with three different types of porous media. Different dyes are used and raw images are acquired. Subsequently, these are used to evaluate a suite of image processing and analysis approaches, which are organized into a workflow. Our approach allows for us to identify key image processing and analysis steps that introduce the most error. Applicable dye configurations led to uncertainties of up to 41% depending on the selection of processing steps. Based on these findings, it was possible to delineate a flexible framework for image processing and analysis that has the potential for transfer and application in other tank experiment setups.