Contaminant Plume Research Articles

Effectiveness and impact factors of passive convergence-permeable reactive barrier (PC-PRB): Insights from tracer simulation study

Passive convergence-permeable reactive barrier (PC-PRB) represents a green and sustainable technology for in-situ remediation of contaminated groundwater. A laboratory-scale PC-PRB tracer simulation system was established to quantify its contaminant plume capture performance using image analysis method. Results indicate that PC-PRB captures the plume 65% wider than C-PRB, which means that fewer PRB sizes and materials volume would be necessary to treat an equivalent contaminated plume. This improvement is due to a significant drawdown within the PC-PRB's passive well, known as the passive hydraulic decompression-convergent flow effect. We further evaluated the effects of water pipe length, hydraulic gradient, and media particle size on PC-PRB's plume capture performance. Results indicate that an increased water pipe length enhances the PC-PRB's plume capture capacity due to greater well drawdown. PC-PRB not only captures the plume but also acts as a hydraulic barrier. The retardation effect of PC-PRB on plume migration increases with water pipe length. Conversely, both hydraulic gradient and media particle size impact the plume capture capacity of PC-PRB by modifying groundwater flow velocity and pollutant dispersion. An increase in either hydraulic gradient or media particle size decreases the plume capture performance of PC-PRB. Therefore, PC-PRB technology may be more effective in contaminated sites characterized by low hydraulic gradients and permeability. Overall, PC-PRB demonstrates significant effectiveness in enhancing plume capture performance, which can notably reduce remediation costs and environmental footprint, broadening its application scope.

The role of mobile-immobile sorbents on flow and colloid-facilitated contaminant transport through porous media: Two-dimensional modeling

A mathematical model based on conceptual equilibrium has been devised to study the role of mobile-immobile sorbents on colloid-facilitated contaminant transport through a two-dimensional, saturated, and homogeneous groundwater system with transient flow conditions. The primary drawback of the kinetic model arises from its inability to handle the large number of parameters involved in multi-dimensional regions at bigger scales when several transport processes and transient flow conditions are included. In order to prevent model complexity, the current study is the first attempt to investigate the equilibrium approach with transient flow conditions in a two-dimensional domain. Conventionally, colloids enhance the contaminant transport by adsorbing contaminants through multiple sorption reactions. However, it was noticed that even in the presence of colloids (Cc=100mg/l), BTCs get attenuated due to an increase in k3, which indicates a stronger affinity of pollutants with the static solid matrix. Therefore, the present study suggests that the mere existence of colloids may not guarantee enhanced contaminant movement. A contaminant plume with 50 mg/l of colloids travels 2 m ahead of the contaminant without colloids, and this distance increases to 4.7 m for 100 mg/l colloids. Therefore, colloids can either be favorable in the case of contaminated aquifer rehabilitation or adverse in the event of groundwater contamination spread.

Quantification of contaminant mass discharge and uncertainties: Method and challenges in application at contaminated sites

Contaminant mass discharge (CMD) estimation involves combining multilevel concentration and flow measurements to quantify the contaminant mass passing through a control plane downgradient of a point source. However, geological heterogeneities and limited data introduce uncertainties that complicate CMD estimation and risk assessment. Although CMD is increasingly used in groundwater management, methods for quantifying and handling these uncertainties are still needed. This study develops and tests a CMD estimation method based on Bayesian geostatistics to quantify CMD uncertainties using data from a control plane perpendicular to the contaminant plume.By combining geostatistical conditional simulations of the spatial concentration distribution with the flow, an ensemble of CMD realizations is generated, from which a cumulative distribution function is derived. A key element of this approach is the use of a macrodispersive transport model to simulate the spatial concentration trend. This ensures that the estimated concentration reflects the expected physical behavior of the contaminant plume while also allowing the integration of site-specific conceptual information.The method is applicable to plumes with dissolved contaminants, such as chlorinated solvents, petroleum hydrocarbons, Per- and polyfluoroalkyl substances (PFAS) and pesticides. Site-specific conceptual understanding is used to inform the prior probability density functions of the structural model parameters and to define acceptable simulated concentration limits. We applied the method at three sites contaminated with chlorinated ethenes, demonstrating its robustness across varying information levels and data availability.Our results shows that strong site-specific conceptual knowledge and high sampling density constrain the CMD uncertainty (CV = 21 %) and results in estimated model parameters and a spatial concentration distribution that agrees well with the conceptual model. For a site with less data and limited conceptual knowledge, CMD and concentration distribution estimates are still feasible, though with higher uncertainty (CV = 41 %). Extending the method to account for multiple source zones and complex plume migration improved parameter identification and reduced the 95 % CMD confidence interval by 11 % ([4950–8750] to [5090–8480] g yr−1), while also providing a spatial concentration distribution in better agreement with the plume conceptualization.This study highlights the importance of integrating site-specific conceptual knowledge in CMD estimation, particularly for less-sampled sites. The method can furthermore assist in identifying remediation targets, evaluating remedial effectiveness, and optimizing sampling strategies.

Spatial thinking skills used by hydrogeology practitioners and students while completing a hydrogeology task

AbstractA typical classroom exercise in hydrogeology is to develop a conceptual model of a contaminated site, identify groundwater flow direction(s), and predict the location and mass of a contaminant plume. This requires knowledge of key hydrogeological concepts and is highly visuospatial in nature. Among multiple discrete spatial thinking skills identified by cognitive science, the combination of visual penetrative ability and working in multiple frames of reference were identified to significantly predict performance on a hydrogeology task and showed that together with hydrogeology knowledge, these spatial thinking skills account for 49% of the variability on task performance. Seventy-two hydrogeology practitioners and students with varying levels of expertise were administered multiple spatial thinking tests and an assessment of hydrogeology knowledge before completing a hydrogeology task that was developed for the study. Using spatial thinking and knowledge test scores as predictor variables, a hierarchical regression analysis was conducted with performance on the hydrogeology task as the outcome variable. The resulting model predicts that at low levels of hydrogeology knowledge, the identified spatial thinking skills account for more than a 25% difference on the hydrogeology task. This study provides empirical evidence that visual penetrative ability and working in multiple frames of reference are important skills in hydrogeology; thus, instructors are encouraged to recognize that underdeveloped spatial thinking skills could present hurdles for students and that targeted spatial thinking training may yield positive results for both weak and strong spatial thinkers.

Unveiling heavy metal(loid) contamination and migration at an abandoned smelting site: Integrated geophysical and hydrological analyse

The heavy metal contamination at abandoned smelting sites is receiving increasing attention. This study focuses on an abandoned lead–zinc smelting site and uses several methods (such as soil sampling, geophysical exploration, and groundwater monitoring) to systematically investigate the contamination characteristics, distribution patterns, and migration mechanisms of heavy metal(loid) in soil and groundwater. A geological structure model of the study area was established through drilling data and geophysical techniques (Ground penetrating radar (GPR) and Electrical resistivity tomography (ERT)), revealing the relationships among soil resistivity, electromagnetic wave reflection characteristics, and geological structure. This study revealed that heavy metal contamination in soil was concentrated mainly northwest of the site and was closely related to historical smelting activities. The presence of heavy metal(loid)s in groundwater had a complex correlation with soil physicochemical properties and environmental elements, and its distribution was influenced by both hydrogeological conditions and soil properties. This study also investigated the effects of the soil water content (SWC) and soil hydraulic conductivity (Ks) on heavy metal(loid)s migration and reported that high SWC and Ks contribute to the dissolution and migration of heavy metals. In addition, a 3D visualization model of heavy metal contamination plumes was established via Voxler software. This model intuitively demonstrated the migration and diffusion of heavy metal(loid)s in the underground environment, offering a novel perspective on their subsurface behavior. The research results provide important information for understanding the migration of heavy metal(loid)s in soil-groundwater systems.

Time‐lapse inversion of self‐potential data through particle filtering

AbstractIn environmental sciences, comprehending the movement of subsurface contaminants is crucial for formulating effective remediation measures. The self‐potential (SP) method has become a common tool for delineating landfill contamination plumes. Contaminant diffusion or migration represents dynamic processes, with corresponding SP responses evolving over time. However, conventional SP interpretation approaches have predominantly relied on static single‐frame inversion, overlooking the temporal correlation in time‐series SP data and resulting in cumulative errors. To tackle this challenge, we introduce a novel method for time‐lapse inversion of SP data leveraging particle filtering. This approach recursively refines the priori state model through posteriori observations to achieve precise estimations of dynamic models. Specifically, a spherical polarization model is deployed to establish the state equations of underground contaminant diffusion and transport models, whereas the observation model is derived through forward modeling. The proposed method is validated using two synthetic examples and one lab‐measured dataset. The findings demonstrate the efficacy of the time‐lapse inversion algorithm in precisely estimating dynamic models, outperforming static single‐frame inversion based on the particle swarm optimization algorithm. The posteriori distribution of particles approximates a bell‐shaped distribution, with the true state closely positioned near the peak probability. Therefore, we affirm that conducting time‐lapse inversion of time‐series SP data through particle filtering is an effective and dependable approach for accurately estimating dynamic model states.

Characterizing Flow Dynamics in a Two-dimensional Flow Cell: Developing Methodologies for Enhanced Subsurface Contaminant Remediation Research

This study used a two-dimensional (2-D) flow cell apparatus (25 cm × 25 cm × 1 cm) to characterize water flow and contaminant transport in porous media. The primary objective was to enhance the understanding of flow dynamics and contaminant behavior in controlled 2-D environments, which is crucial for optimizing the application of remediation technologies to treat contaminants in saturated and unsaturated soil conditions. The 2-D cell was designed with transparent walls to simulate subsurface conditions using medium grained sand, allowing for precise control and observation of fluid movement under saturated and quasi-saturated (trapped air) conditions. Visualization techniques and flow measurement tools were employed to capture detailed data on flow patterns, velocities, and dispersion characteristics. Analysis of breakthrough curves for salt tracer tests using electrical conductivity measurements, alongside varying flow rates, flow fractions, and port configurations, provided significant insights into fluid behavior and contaminant transport. The study also explored the effects of varying salt concentrations (and solution densities) on flow dynamics, contributing to a better understanding of the role of contaminant concentration in contaminant plume development and remediation processes. Visual analysis using dyed contaminant complemented the quantitative data, offering real-time observations of flow patterns and contaminant distribution. A custom-built apparatus was used to control the hydraulic head of multiple outlet ports to sample water from different depths. The experiments demonstrated that gravity plays a crucial role in solute (salt) breakthrough, with lower ports showing earlier breakthrough. Flow fields formed vertically, from bottom-to-top, and horizontally, from inlet-to-outlet, with lower portions advancing sooner. However, multi-port configurations, especially involving five ports, disrupted the sequential breakthrough order, possibly due to convection in the clear well that formed the outlet boundary of the flow cell. These findings will support research students in conducting similar experiments with flow and transport that varies with depth, with future insights for professionals in developing and optimizing remediation strategies related to subsurface contaminant flow. Special thanks are extended to Liam Price, Julia-Barnes James, and Brooke Belfall who were helpful during the research and the experimentation process.

Contaminant mass discharge estimation of a sulfonamide plume by use of hydraulic profiling tool (HPT) and fluorescence techniques

The contaminant mass discharge is a relevant metric to evaluate the risk that a groundwater plume poses to water resources. However, this assessment is often vitiated by a high uncertainty inherent to the assessment method and often limited number of measurement points to carry out the assessment. Direct-Push techniques in combination with profiling tools and dedicated sampling can be an interesting alternative to increase the measurement point density and hence reduce the mass discharge uncertainty. The main objective of our study was to assess if DP logging and sampling could be employed to get a reasonable estimate of contaminant mass discharge in a large sulfonamide contaminant plume (> 1500 m wide), compared to a more traditional approach based on monitoring wells. To do so, an Hydraulic Profiling Tool (HPT) logging with a dedicated site calibration was used to estimate the hydraulic conductivity field. The sulfonamide concentrations were inferred from the compound fluorescence properties measured by laboratory spectrofluorometry (λEx / λEm = 255/340 nm) and a dedicated log-log linear regression model. Our results show that HPT-derived hydraulic conductivity values are in good agreement with the monitoring well results, and within the order of magnitude reported in similar studies or indirect geophysical techniques. Fluorescence appears as a powerful proxy for the sulfonamide concentration levels. Ultimately, the contaminant mass discharge estimate from HPT and fluorescence techniques lies within a factor 2 from the estimate by monitoring wells, with 549 [274–668] and 776 [695–879] kg/yr respectively. Overall, this study highlights that DP logging tools combined with indirect methods (correlation with fluorescence) could provide a relevant contaminant mass discharge estimate for some optically active substances, given that a proper calibration phase is carried out.

Aqueous benzene, toluene, ethylbenzene and xylene (BTEX) removal using core-shell structure activated carbon ball as a permeable reactive barrier material

Permeable reactive barriers (PRBs) are passive and sustainable treatment systems for remediating the diffusion of contaminant plumes in groundwater. Several conventional reactive materials such as activated carbon (AC) have long been used as reactive media for PRBs. AC, which is known for its high adsorption capability and cost-effectiveness, is commonly used to remove multiple pollutants from groundwater. Unfortunately, among the reactive materials, AC can fill in the barrier and pose practical problems, such as a pressure drop, solid losses during handling, and safe disposal of filled sorbents, because of its low particle strength. In this study, AC balls were prepared using zeolite as the core and powdered AC, quartz, and calcite as the shell. AC ball with excellent mechanical strength and high permeability properties in the form of a core–shell layer is a good alternative to conventional reactive materials. The adsorption characteristics of benzene, toluene, ethylbenzene, and p-xylene (BTEX) from solutions using AC balls were investigated. The adsorption equilibrium is in the order of X > E > T > B. The adsorption capacity for B (4.90 mg/g), T (3.37 mg/g), E (1.16 mg/g) and X (0.46 mg/g) were observed at solution pH = 7. To validate the proposed models, batch experiments indicated that the pseudo-2nd-order (11.61 and 10.54 mg/g) and Langmuir models (8.81 and 5.27 mg/g) were the most suitable for describing the kinetics and equilibrium of benzene and toluene, respectively. Regeneration experiments were performed using chemical extraction (methanol) and microwave (MW) heating. When the AC ball was reused six times, the removal efficiency for benzene was maintained above 81 % using methanol, whereas that for benzene was increased to 91 % using MW. A series of experiments revealed that AC balls have considerable reusability. The AC ball can be effectively used as a medium for the removal of benzene in a continuous flow system such as column and sandbox. Based on these results, AC balls are a potential reactive medium for field-scale PRB practical remediation applications.

Coupled hydrogeophysical inversion through ensemble smoother with multiple data assimilation and convolutional neural network for contaminant plume reconstruction

In the field of groundwater, accurate delineation of contaminant plumes is critical for designing effective remediation strategies. Typically, this identification poses a challenge as it involves solving an inverse problem with limited concentration data available. To improve the understanding of contaminant behavior within aquifers, hydrogeophysics emerges as a powerful tool by enabling the combination of non-invasive geophysical techniques (i.e., electrical resistivity tomography—ERT) and hydrological variables. This paper investigates the potential of the Ensemble Smoother with Multiple Data Assimilation method to address the inverse problem at hand by simultaneously assimilating observed ERT data and scattered concentration values from monitoring wells. A novelty aspect is the integration of a Convolutional Neural Network (CNN) to replace and expedite the expensive geophysical forward model. The proposed approach is applied to a synthetic case study, simulating a tracer test in an unconfined aquifer. Five scenarios are compared, allowing to explore the effects of combining multiple data sources and their abundance. The outcomes highlight the efficacy of the proposed approach in estimating the spatial distribution of a concentration plume. Notably, the scenario integrating apparent resistivity with concentration values emerges as the most promising, as long as there are enough concentration data. This underlines the importance of adopting a comprehensive approach to tracer plume mapping by leveraging different types of information. Additionally, a comparison was conducted between the inverse procedure solved using the full geophysical forward model and the CNN model, showcasing comparable performance in terms of results, but with a significant acceleration in computational time.

The Result of Electrical Resistivity Investigation of Solid Waste Landfill Conducted at Nsukka Municipality Using Electrical Resistivity Measurments

Solid waste landfill management has been a significant issue for Nigerian urban areas and other developing countries across the globe.Similar to most other cities, Nsukka also generates waste on a daily basis, much of which is dumped in poorly designed and positioned dumping sites. The majority of the disposal sites are found on roadsides, at marketplaces, on farms, and in residential neighborhoods, among other places. The road infrastructure and groundwater are under danger, and the beauty of the impacted communities are not spared. Undoubtedly, the unchecked citation of boreholes as the source of potable water in the majority of our rural and urban communities—given that the government doesn't seem to be providing water to the people—has become a significant challenge. An investigation using electrical resistivity method was conducted around a solid waste dumpsite at Nsukka in Nsukka L.G.A of Enugu State, Nigeria with an aim to investigate the level of groundwater contamination and the objectives to determine the subsurface geoelectric layers, depth to water table, lithology delineation and map out the contamination zones. The scope of this study provides an overview of some of the approaches used to assess the aquifer vulnerability and aquifer potential using Vertical Electrical Sounding (Schlumberger array) and 2D resistivity imaging (Wenner array) in different locations around Nsukka municipality dumpsite. Both methods were used for this study in order to provide a geophysical database for exploration of the study area’s groundwater resources and also they are less expensive and less time consuming. VES has proved to be effective in solving groundwater problems in most places in Nigeria (Ezeh and Ugwu, 2010; Ugwu and Ezeh, 2012; Nzemeka et al. [1,2]. Electrical Sounding (VES) and 2D resistivity imaging were carried out with a digital read out resistivity meter (ABEM SAS 1000) to acquire data in the area and were interpreted using the Schlumberger automatic INTERPEX analysis software and the RES2DINV software respectively, which generates model curves using initial layer parameters and display the variations of electrical resistivities respectively. A total of eight (8) sounding and six (6) 2D resistivity imagings were carried out in the area. A contaminant leachate plume was delineated in 2D resistivity sections as low resistivity zones while the VES shows the depth of aquifer. In 2D pseudosections where bluish colours with low resistivities (less than 20.80\(\Omega\)m) with the depth ranging from 1.28m to 17.1m in the Line 1 and 2 are seen as contaminated zones. The rest of the lines are not contaminated because of their high resistivities (greater than 20.80\(\Omega\)m). The result of the electrical resistivity survey also showed 4 - 5 layers geo-electric sections and an AA and AK type sounding curves. The VES result shows that VES 1A, 1B, 2A and 2B which are carried out on line 1 & 2 of the wenner lines showed signs of contamination with low resistivity values less than 20.80\(\Omega\)m complementing the wenner results. The contamination has not yet got to where the aquifer is located on the lines. Since the depth to the aquifer ranges from 30.26m to 155.43m while maximum depth of contamination is 17.1m. It is believed that the leachate has not percolated down to the aquiferous zones as such aquifers are presumed to be free. As such, it is recommended that boreholes around the study area should not be less than 30m deep to avoid exploiting polluted water.

Insights into In Situ Benthic Caging Tests for Ecotoxicity Assessments Targeting Discharging Groundwater Contaminant Plumes

While in situ toxicity testing with caged organisms has been used to assess surface water and sediment contamination, no successful application to benthic organisms exposed to highly contaminated groundwater plumes discharging to surface waters has been reported. The objective of this study was to demonstrate and evaluate this application using four sets of tests performed at three previously reported contaminated groundwater sites, which include one river site affected by volatile organic contaminant plumes, and two sites, one pond and one small urban stream, impacted by landfill plumes. The study examined multiple cage designs and orientations and two test organisms: an amphipod (Hyalella azteca) and midge larvae (Chironomus riparius; only one study). Cages were deployed for between 5 and 28 days and assessed for organism survival and growth. At all sites and for some deployment conditions, cages exposed to high contaminant concentrations in the plume footprint had greater mortality compared to those exposed to lower or background concentrations. Organism growth was less clear as a metric of toxicity. Vertically oriented cages typically showed high mortality to plume contaminants, but some were also affected by other non-target groundwater conditions (e.g., low dissolved oxygen, other contaminant sources), while horizontally oriented cages were rarely responsive to either groundwater influence. A hybrid cage design showed much promise in its single study. Useful observations on the test organisms and on potentially problematic site conditions were also made. The informed use of in situ toxicity cages could be an additional beneficial tool for groundwater contaminated site assessments.

Dense Contaminants Mixing Into the Saltwater Wedge in Coastal Aquifers: Laboratory and Numerical Investigations

AbstractThe saltwater‐freshwater mixing zones in coastal aquifers can host complex physical exchange and biogeochemical transformations. The land‐sourced dense contaminant plumes could be transferred into the mixing zone of the saltwater wedge due to the density effect prior to discharge to the sea. However, the mixing process between dense contaminants and the saltwater wedge has not received much attention, largely due to the lack of physical evidence. This study used laboratory experiments and numerical simulations to investigate the transport and discharge behaviors of variable‐density contaminant plumes in tidally influenced unconfined coastal aquifers. Results demonstrate that the highly dense contaminants mix with the underlying saltwater and finally merge with the saltwater wedge. This process significantly extends the contaminant discharge durations, thereby reducing the peak value of contaminant efflux. The dense contaminants are elongated along the landward margin of the saltwater wedge, leading to a larger spreading area (Ms) than that of constant‐density contaminants. The sensitivity analysis indicates that the high density of contaminants acts as a trigger to induce the mixing of them and wedges. The higher hydraulic conductivity, lower dispersivities and reduced inland freshwater flux significantly increase the residence times (Rt) and discharge duration (Dt) by enhancing the mixing of dense contaminants with seawater. In contrast, both Rt and Dt values are not only non‐monotonic functions of tidal amplitudes but also less sensitive to tidal effects. Compared with the non‐tidal condition, however, the addition of tides significantly increases both Rt and Ms values of dense contaminant plumes. The results presented herein provide valuable insights into the mechanisms of dense contaminants mixing into saltwater wedges, which could guide practitioners in designing effective strategies to protect coastal environments from land‐sourced contaminants.

Wykorzystanie metody radaru penetrującego grunt w badaniu gleby zanieczyszczonej węglowodorami w obszarze Navodari - Rumunia

Ground penetrating radar (GPR) is a very useful geophysical method for use in hydrogeologic and near-surface mapping studies. It can be used to study contaminants in groundwater, subsurface faulting, and underground cavities (natural or man-made), all of which pose potentially dangerous geological hazards. The GPR technique is similar in principle to seismic reflection and sonar techniques. The propagation of the radar signal depends on the frequency-dependent electrical properties of the ground.Electrical conductivity of the soil or rock materials along the propagation paths introduces significant absorptive losses which limit the depth of penetration into the earth formations and is primarily dependent upon the moisture content and mineralization present.Reflected signals are amplified, and transformed to the audio-frequency range, recorded, processed, and displayed. From the recorded display, subsurface features such as soil/ soil, soil/rock, and unsaturated/saturated interfaces can be identified. In addition, the presence of floating hydrocarbons on the water table, the geometry of contaminant plumes, and the location of buried cables, pipes, drums, and tanks can be detected. The GPR data are presented as a two-dimensional depth profile along a scanned traverse line in which the vertical axis is two-way travel time measured in nanoseconds. The location of hydrocarbon contamination in the ground using the GPR method is based mainly on information taken from reflected signals. In the cases investigated in Romania contaminated sites (Navodari area), such signals were very rarely recorded. A long time after spillage, contamination takes the form of plumes with different size and distribution, which depends on the geological and hydraulic properties of the ground. The survey discussed in this paper was carried out using the GPR system-Noggin with two antennas (250 and 500mHz) Data collected were processed using software(EKKO_Project™ GPR Data Analysis) to produce 2D radargram in time scale. The presence of contaminant plumes as well as the water table are observed in the GPR sections at depths approximately of 0.5 to 1.5 m. In the GPR section, the oil contaminated layer exhibits discontinuous, subparallel, and chaotic high amplitude reflection patterns. Promising results were also obtained in the GPR survey where three obvious reflection patterns representing the top sand-silt layer, oil-contaminated zone and, the underlying thick soft clay were detected in all 2D radargrams of the GPR traverse lines.

3D GeoRemediation: A Digital Hydrogeophysical–Chemical Clone and Virtual Hydraulic Barrier with Groundwater Circulation Wells (GCWs) for Groundwater Remediation

Identification of contamination sources and delineation of plumes in the geological environment stand as pivotal elements in reconstructing the conceptual site model (CSM) and devising remediation strategies tailored to specific physicochemical traits. This study endeavors to showcase the capabilities of a 3D digital interface, seamlessly integrating multi-source data, to elucidate site-specific contamination dynamics and steer the implementation of remediation strategies harmoniously aligned with the ethos of remediation geology. In a site historically marred by chlorinated solvent contamination, the digitization of stratigraphic, piezometric, chemical, and membrane interface probe (MIP) data underpins geomodeling endeavors and yields a meticulously crafted, data-driven CSM. The hydrogeochemical and hydrogeophysical data were interpolated to build a volumetric, digital 3D model illustrating data-driven elements. The comprehensive 3D clone adeptly delineates secondary contamination sources and renders visible the contamination plume within a georeferenced framework, mirroring the nuanced interplay of stratigraphic nuances and groundwater path. A data-centric approach to modeling facilitates the design of the first hydraulic virtual barrier leveraging groundwater circulation well (GCW) technology, its geometry finely attuned to intercept the contamination plume originating from source dissolution and aligning with preferential groundwater flow trajectories. Conventional hydrochemical monitoring and multilevel sampling substantiate the discernible reduction in chlorinated solvent concentrations across various depths within the aquifer horizon, affirming the efficacy of GCWs in their virtual barrier configuration. The findings highlight the effectiveness and limited groundwater consumption of the virtual barrier compared to the on-site pump-and-stock system. This research underscores the potency of a multi-faceted evidence-driven puzzle in conceptualizing contamination mechanisms within the geological milieu, thereby fostering the application of cutting-edge, effective, and sustainable remediation strategies.

Determination of contamination sources and geochemical reactions in groundwater of a mine area using Cu, Zn, and S–O isotopes

To determine contamination sources and pathways, the use of multiple isotopes, including metal isotopes, can increase the reliability of environmental forensic techniques. This study differentiated contamination sources in groundwater of a mine area and elucidated geochemical processes using Cu, Zn, S–O, and O–H isotopes. Sulfate reduction and sulfide precipitation were elucidated using concentrations of dissolved sulfides, δ34SSO4, δ18OSO4, and δ66Zn. The overlying contaminated soil was possibly responsible for the contamination of groundwater at <5 mbgl, which was suggested by low δ65Cu values (0.419–1.120‰) reflecting those of soil (0.279–1.115‰). The existence of dissolved Cu as Cu(I) may prevent the increase in δ65Cu during leaching of contaminated soil in the sulfate-reducing environment. In contrast, the groundwater at >5 mbgl seemed to be highly affected by the contamination plume from the adit water, which was suggested by high SO42− concentrations (407–447 mg L−1) and δ65Cu (0.252–2.275‰) and δ66Zn (−0.105‰–0.362‰) values at a multilevel sampler approaching those of the adit seepages. Additionally, the O–H isotopic ratios were distinguished between <5 mbgl and >5 mbgl. Using δ65Cu and δ66Zn to support the determination of groundwater contamination sources may be encouraged, particularly where the isotopic signatures are distinct for each source.

2D model of groundwater flow and total dissolved HCH transport through the Gállego alluvial aquifer downstream the Sardas landfill (Huesca, Spain)

The organic pollutants disposed at the Sardas landfill in Sabiñánigo (Huesca, northeastern Spain) by the INQUINOSA lindane factory have reached the Gállego alluvial aquifer and could affect the Sabiñánigo reservoir. The daily oscillations of the reservoir water level produce a tidal effect on the piezometric heads of the aquifer. These oscillations are transmitted in a damped way with a time lag, thus attesting that the silting sediments of the reservoir and the natural silts of the Gállego alluvial are interposed between the reservoir water and the layer of sands and gravels. A 2D finite element groundwater flow and total dissolved hexachlorocyclohexane (HCH) transport model through the Gállego alluvial aquifer is presented here. The flow model was constructed to: (1) Quantify the tidal effect, produced by the daily fluctuations of the reservoir water level on the aquifer; (2) Estimate the hydrodynamic parameters of the layer of sands and gravels; and 3) Estimate the vertical hydraulic conductivity of the silting sediments and silts; and (4) Quantify aquifer/reservoir interactions. The flow model reproduces the dynamics of the tidal effect and attests that groundwater velocity and flow direction changes daily in response to the oscillations of the reservoir level. Model results reproduce the measured well hydrographs and the Darcy velocity derived from tracer tests and confirm the validity of the conceptual model. The transport model of total dissolved HCH simulates the time evolution of the contaminant plume. The computed concentrations of total dissolved HCH and the contaminant mass outflux are very sensitive to changes in the source terms and the distribution coefficient, Kd of HCH. The best fit to the measured HCH plumes in September 2010 and December 2020 is obtained with a Kd ranging from 1 to 3 L/kg. The computed flux of dissolved HCH leaving the Sardas site in 2020 towards the Sabiñánigo reservoir ranges from 0.6 kg/year for Kd = 3 L/kg to 3.1 kg/year for Kd = 1 L/kg. The findings of this study will be most useful for planning and designing remedial and containment actions at the Sardas site and other similar lindane-affected sites.

Unravelling the influence of heterogeneities and abstraction on groundwater flow and solute transport in a fractured carbonate aquifer, Sicily, Italy

The migration of a groundwater contaminant plume consisting of light nonaqueous phase liquids (LNAPLs) along the Ionian coastline of Sicily, Italy, has been reported to follow a trajectory that is inconsistent with the regional hydraulic gradient. The influence of several faults affecting the fractured carbonate bedrock aquifer, and groundwater abstraction from a well, were hypothesized to be responsible for the anomalous trajectory of the contaminant plume. A conceptual hydrogeological model was developed for the study area that incorporated structural information derived from geophysical surveys and the mapping of fractures in bedrock outcrops. This conceptual model was incorporated into numerical groundwater flow and solute transport models to simulate the groundwater transport of the light nonaqueous phase liquids. Four model scenarios representing different levels of complexity were tested to assess the relative influence of geological heterogeneity and groundwater abstraction on the migration of the contaminant plume. Results show that underground major discontinuity systems, invoking the presence of the faults in the model domain accounted for the observed migration of the contaminant plume, act as conduits for groundwater flow. Conversely, groundwater abstraction from a well was found to result in relatively minor, localized impacts to the migration of the contamination plume. This study demonstrates the importance of incorporating geological heterogeneity into groundwater modelling and environmental risk assessments associated with the storage of LNAPLs.

Laboratory-scale characterization of slow-release permanganate gel (SRP-G) for the in-situ treatment of chlorinated-solvent groundwater plumes

Significant challenges remain for the remediation of chlorinated-solvent plumes in groundwater, such as trichloroethene (TCE) and tetrachloroethene (PCE). A novel slow-release permanganate gel (SRP-G) technique may show promise for the in-situ treatment (remediation) of chlorinated contaminant plumes in groundwater. A series of laboratory experiments were conducted to characterize the primary physical factors that influence SRP-G gelation processes to optimize SRP-G performance for plume treatment. Specifically, experiments were conducted to quantify gel zeta potential, particle size distribution, and viscosity to determine SRP-G gelation characteristics and processes. These experiments tested various concentrations of two SRP-G amendment solutions (NaMnO4 and KMnO4) prepared with 30-wt.% and 50-wt.% colloidal silica to determine such influences on zeta potential, particle size distribution, and viscosity. The results of this study show that SRP-G solutions with low zeta potential and relatively high pH favor more rapid SRP-G gelation. The concomitant interaction of the predominantly negatively charged colloidal silica particles and the positively charged dissociated cations (Na+ and K+) in the SRP-G solution had the effect of stabilizing charge imbalance via attraction of particles and thereby inducing a greater influence on the gelation process. Gel particle size distribution and changes in viscosity had a significant influence on SRP-G solution gelation. The addition of permanganate (NaMnO4 or KMnO4) increased the average particle size distribution and the viscosity of the SRP-G solution and decreased the overall gelation time. SRP-G amendments (NaMnO4 or KMnO4) prepared with 50-wt.% colloidal silica showed more effective gelation (and reduced gelation time) compared to SRP-G amendments prepared with 30-wt.% colloidal silica. Under the conditions of these experiments, it was determined that both the 7-wt.% NaMnO4 solution and 90 mg/L KMnO4 solution using 50-wt.% colloidal silica would be the optimal injection SRP-G solution concentrations for this in-situ treatment technique.

Quantitative assessment of the impact factors upon remediation performance of the PRB technique in conjunction with pumping/recharge

Installation of pumping/recharge wells adjacent to contaminant plume and permeable reactive barrier (PRB) can effectively improve PRB's performance by artificially altering groundwater flow field attracting more contaminants to pass through PRB within a shorten time. However, the impact factors affecting the performance of PRB in conjunction with pumping/recharge (abbreviated as PRB-PR technique) have not been quantitatively assessed. In this study, numerical approach was adopted and multiple scenarios were implemented based on the developed MODFLOW and MT3MDS models to quantitatively evaluate the influences of impact factors on PRB-PR's remediation performance over a hypothetical contaminated site. Simulation results indicated that: (1) PRB should be installed nearby the downstream boundary of contaminant plume; (2) PRB's contaminant removal efficiency (η), pumping/recharge rate (Q), and aquifer hydraulic conductivity (K) affect PRB-PR's performance significantly, and a linear relationship exists between contaminant removal rate (R) and the product of K and Q no matter η changes (a series of fitting equations representing the linear relationship between R and KQ under various η were yielded based on data fitting); and (3) PRB's hydraulic conductivity, pumping/recharge wells' locations, and pumping/recharge wells' quantity do not affect PRB-PR's performance seriously. The equations indicating PRB-PR's remediation performance and key impact factors can be utilized to calculate the minimum values of Q and η that can meet the requirements of remediation objectives of real contaminated sites, contributing to supervising the installation and operation of the PRB-PR technique accordingly.