Laboratory-scale experimental and modelling investigations of 222Rn profiles in chemically heterogeneous LNAPL contaminated vadose zones

Numerical simulation and experimental study of three-phase distribution characteristics of leaked light non-aqueous phase liquid from buried pipelines in soils containing groundwater and gas

Leakage accidents of buried pipelines have become increasingly common due to the prolonged service of some pipelines which have been in use for more than 150 years. Therefore, there is an urgent need for accurate prediction of pollution scope to aid in the development of emergency remediation strategies. This study investigated the distribution of a light non-aqueous phase liquid in soils containing gas and water through numerical simulations and laboratory experiments. Firstly, a three-dimensional porous medium model was established using ANSYS FLUENT, and for the first time, the distribution of gas and groundwater in soil environments was simulated in the model. Subsequently, the distribution of the three phases of diesel, gas, and water in soil was studied with different leakage velocities and it was found that the leakage velocity played a significant role in the distribution. The areas of diesel in soils at 60 min were 0.112 m2, 0.194 m2, 0.217 m2, and 0.252 m2, with corresponding volumes of 0.028 m3, 0.070 m3, 0.086 m3, and 0.106 m3, respectively, for leakage velocities of 1.3 m/s, 3.4 m/s, 4.6 m/s, and 4.9 m/s. Calculation formulas for distribution areas and volumes were also developed to aid in future prevention and control strategies under different leakage velocities. The study also compared the distribution areas and volumes of diesel in soils with and without groundwater, and it was found that distribution scopes were larger in soils containing groundwater due to capillary force. In order to validate the accuracy of the numerical simulation, laboratory experiments were conducted to study the diffusion of oil, gas, and water under different leakage velocities. The results showed good agreement between the experiments and the simulations. The research findings are of great significance for preventing soil pollution and provide a theoretical basis for developing scientifically sound soil remediation strategies.

Radon deficit technique applied to the study of the ageing of a spilled LNAPL in a shallow aquifer

A recent diesel spill (dated January 2019 ± 1 month) in a refilling station is investigated by the Radon deficit technique. The primary focus was on quantifying the LNAPL pore saturation as a function of duration of ageing, and on proposing a predictive model for on-site natural attenuation. A biennial monitoring of the local fluctuating shallow aquifer has involved the saturated zone nine times, and the vadose zone only once. Rn background generally measured in external and upstream wells is elaborated further due to the site characteristics, using drilling logs and phreatic oscillations. Notably, this study marks the first application of the Rn deficit method to produce a detailed Rn background mapping throughout the soil depth. Simultaneously, tests are performed on LNAPL surnatant samples to study diesel ageing. In particular, they are focused on temporal variations of LNAPL viscosity (from an initial 3.90 cP to 8.99 cP, measured at 25 °C, after 34 months), and Rn partition coefficient between the pollutant and water (from 47.7 to 80.2, measured at 25 °C, after 14 months). Rn diffusion is also measured in different fluids (0.092 cm2 s-1, 1.14 × 10-5 cm2 s-1, and 2.53 × 10-6 cm2 s-1 at 25 °C for air, water and LNAPL, respectively) directly. All parameters and equations utilized during this study are introduced, discussing their influence on Radon deficit technique from a theoretical point of view. Experimental findings are used to mitigate the effect of LNAPL ageing and of phreatic oscillations on determination of LNAPL saturation index (S.I.LNAPL). Finally, S.I.LNAPL dataset is discussed and elaborated to show the pollutant attenuation across subsurface over time, induced by natural processes primarily. The proposed predictive model for on-site natural attenuation suggests a half-removal time of one year and six months. The significance of such models lies in their capability to assess site-specific reactions to pollutants, thereby enhancing the effectiveness of remediation efforts over time. These experimental findings may offer a novel approach to application of Rn deficit technique and to environmental remediation of persistent organic compounds.

Vertical Light Non-Aqueous Phase Liquid (LNAPL) distribution by Rn prospecting in monitoring wells

In the frame of a collaboration between the Italian National Research Council (CNR) and Mares s.r.l., a study, about the possibility of determining radon vertical distribution at different soil depths in order to trace light non-aqueous phase liquid (LNAPL) contaminations, was developed. The radon deficit technique, based on the preferential solubility of soil gas radon into non-polar fluids, such as refined hydrocarbons, has been investigated by various theoretical and applied research so far. According to international scientific literature, radon deficit can be used both for geochemical prospection of the spatial irregular NAPL dispersion and for monitoring of remediation activities. Even though it is well known that this type of pollutants can be distributed along the vertical soil profile-firstly due to their density in comparison to water density, and secondly due to fluctuations of shallow aquifers, soil pore size, aging of contamination, and so on-the vertical localization of the plume still represents a scientific challenge. In this article, a method to determine the radon vertical profile is tested and applied to assess the potential use of the radon deficit technique in the vertical detection of pollutant presence for the first time in a fuelling station. Two LNAPL-contaminated sites were selected for a pilot test. Experimental findings seem to support the use of vertical radon geochemical prospection to delimit the depth range of a LNAPL pollution directly. Systematic data collection and modeling may lead to a 3D reconstruction of the dispersion of contaminant in different soil levels.

Understanding of Co-boiling between Organic Contaminants and Water during Thermal Remediation: Effects of Nonequilibrium Heat and Mass Transport

In situ thermal desorption (ISTD) provides an efficient solution to remediation of soil and groundwater contaminated with nonaqueous phase liquids (NAPLs). Establishing a relationship between the subsurface temperature rise and NAPL removal is significant to reduce energy consumption of ISTD. However, the co-boiling phenomenon between NAPL and water poses a great challenge in developing this relationship due to the nonequilibrium heat and mass transport effects. We performed a systematic experimental investigation into the local temperature rise patterns at different distances from a NAPL pool and under different degrees of superheat by selecting four representative NAPLs (i.e., trichloroethylene, tetrachlorethylene, n-hexane, and n-octane) according to their density and boiling point relative to water. The patterns of temperature rise indicated that the underground temperature field can be divided into three zones: the zone of local thermal equilibrium, the nonequilibrium zone affected by co-boiling, and the zone unaffected by co-boiling. We developed a pattern-recognition-based approach, which considers the effects of local heat and mass transport to establish a qualitative correlation between the temperature rise and NAPL removal. Our results give deeper insights into the understanding of subsurface temperatures in ISTD practice, which can serve as the guideline for more accurate and sustainable remediation.

Using groundwater monitoring wells for rapid application of soil gas radon deficit technique to evaluate residual LNAPL

The application of the 222Radon (Rn) deficit technique using subsurface soil gas probes for the identification and quantification of light non-aqueous phase liquids (LNAPL) has provided positive outcomes in recent years. This study presents an alternative method for applying this technique in the headspace of groundwater monitoring wells. The developed protocol, designed for groundwater monitoring wells with a portion of their screen in the vadose zone, is based on the use of portable equipment that allows rapid measurement of the Rn soil gas activity in the vadose zone close to the water table (i.e., smear zone) where LNAPL is typically expected. The paper first describes the step-by-step procedure to be followed for the application of this method. Then, a preliminary assessment of the potential of the method was carried out at two Italian sites characterized by accidental gasoline and diesel spills into the subsurface from underground storage tanks. Although the number of tests conducted does not allow for definitive conclusions, the results obtained suggest that, from a qualitative point of view, Rn monitoring in the headspace of monitoring wells is a promising, fast, and minimally invasive screening method that could also potentially reduce the costs associated with field data acquisition. This method proves to be suitable for detecting the presence of LNAPL in both the mobile and residual phases with results consistent with the other lines of evidence available at the sites, such as groundwater and soil gas monitoring. Future efforts should be directed toward evaluating the accuracy of this method for a quantitative assessment of residual LNAPL saturations.

Influence of advection on the soil gas radon deficit technique for the quantification of LNAPL

The Radon (Rn) deficit technique is a rapid, low-cost, and non-invasive method to identify and quantify light non-aqueous phase liquids (LNAPL) in the soil. LNAPL saturation is typically estimated from Rn deficit using Rn partition coefficients, assuming equilibrium conditions. This work examines the applicability of this method in the presence of local advective fluxes that can be generated by groundwater fluctuations or biodegradation processes in the source zone. To this end, a one-dimensional analytical model was developed to simulate the steady-state diffusive-advective transport of soil gas Rn in the presence of LNAPL. The analytical solution was first validated against an existing numerical model adapted to include advection. Then a series of simulations to study the effect of advection on Rn profiles were carried out. It was found that in high-permeability soils (such as sandy soils), advective phenomena can significantly affect Rn deficit curves in the subsurface compared with those expected, assuming either equilibrium conditions or a diffusion-dominated transport. Namely, in the presence of pressure gradients generated by groundwater fluctuations, applying the traditional Rn deficit technique (assuming equilibrium conditions) can lead to an underestimation of LNAPL saturation. Furthermore, in the presence of methanogenesis processes (e.g., in the case of a fresh LNAPL of petroleum hydrocarbons), local advective fluxes can be expected above the source zone. In such cases, Rn concentrations above the source zone can be higher than those above background areas without advective phenomena, resulting in Rn deficits higher than 1 (i.e., Rn excess), and thus leading to a wrong interpretation regarding the presence of LNAPL in the subsurface if advection is not considered. Overall, the results obtained suggest that advection should be considered in the presence of pressure gradients in the subsurface to ensure an effective application of the soil gas Rn-deficit technique for quantitative estimation of LNAPL saturation.

Groundwater and soil contamination by LNAPL: State of the art and future challenges

Contamination by Light Non-Aqueous Phase Liquids (LNAPL) represents a challenge due to the difficulties encountered in its underground assessment and recovery. The major risks arising from subsoil LNAPL accumulation face human health and environment, gaining a social relevance also in the frame of a continuously changing climate. This paper reports on a literature review about the underground contamination by LNAPL, with the aims of providing a categorization of the aspects involved in this topic, analyzing the current state of the art, underlying potential lacks and future perspectives. The review was focused on papers published in the 2012-2022 time-interval, in journals indexed in Scopus and WoS databases, by querying "LNAPL" within article title, abstract and/or key words. 245 papers were collected and classified according to three "key approaches" -namely laboratory activity, field based-data studies and mathematical simulations- and subordinate "key themes", so to allow summarizing and commenting the main aspects based on the application setting, content and scope. Results show that there is a wide experience on plume dynamics and evolution, detection and monitoring through direct and indirect surveys, oil recovery and natural attenuation processes. Few cues of innovations were found regarding both the use of new materials and/or specific field configuration for remediation, and the application of new techniques for plume detection. Some limitations were found in the common oversimplification of the polluted media in laboratory or mathematical models, where the contamination is set within homogeneous porous environments, and in the low number of studies focused on rock masses, where the discontinuous hydraulic behavior complicates the address and modeling of the issue. This paper represents a reference for a quick update on the addressed topic, along with a starting point to develop new ideas and cues for the advance in one of the greatest environmental banes of the current century.

Migration and redistribution of LNAPL in inclined stratified soil media

Leakage of light non-aqueous phase liquid (LNAPL) into soil can cause serious environmental issues. In this study, a two-dimensional device with adjustable dip angles was designed to investigate the migration and redistribution of LNAPL in natural inclined stratified soil media by the light transmission visualization (LTV) technology. The captured experimental images were processed to obtain the diesel distribution based on gray value which could represent the LNAPL saturation distribution. LNAPL may not be able to penetrate through the fine-coarse interface due to the capillary barrier effects. In this case, the vertical and horizontal migration distances (V and H), contaminated area (S), as well as deviation angle (γ) of centroid increased with the dip angle. Increasing the leakage amount to more than 30 mL would result in LNAPL breakthrough at the 10°-inclined interface, leading to much larger V, H, S, and γ than those at 10 mL, while 20-mL LNAPL failed to break through. In the latter case, a lower leakage rate than 10 mL/min would cause larger H and γ but similar V or S in the long term. This study could enrich the understanding of LNAPL contamination in vadose zone, providing reference for the prediction and treatment in realistic inclined contaminated sites.

1D_RnDPM: A freely available 222Rn production, diffusion, and partition model to evaluate confounding factors in the radon-deficit technique

The radon-deficit technique is a powerful tool to detect and delineate sub-surface accumulations of organic contaminants. Field measurements of ²²²Rn in soil air, however, are affected by several confounding factors that can lead to the misinterpretation of results. Among the most influential are: vertical and lateral changes of lithology, fluctuating contaminant saturations with depth, varying water saturation ratios along the soil profile and atmospheric (and, therefore, soil) thermal oscillations. To evaluate and minimize the effect of these confounding factors on the interpretation of the results of the Rn deficit technique, a Matlab® based multi-layer model of ²²²Rn production-partition-diffusion in unsaturated porous media (1D_RnDPM: One-Dimensional ²²²Rn Diffusion and Partition Model) has been developed and is freely available as Supplementary Material in this work. A laboratory protocol has also been proposed to obtain site-specific input parameters for the model, i.e., ²²²Rn equilibrium concentration (as determined by the accumulation chamber method), soil bulk density and soil solid-phase density. The model predictions have been contrasted with field information obtained from successive sampling campaigns in which ²²²Rn in soil air was measured at a site where the vadose zone, consisting of an anthropogenic backfill underlain by a silt layer, is affected by a complex mixture of benzene, phenol, (poly) chlorobenzenes, (poly) chlorophenols and hexachlorocyclohexane isomers, among other compounds. The model has successfully predicted the vertical profile of ²²²Rn concentrations in soil air, including the effect of the oscillations of the water table and of ground-level temperature. The results also underline that ²²²Rn measurements in subsoil air are representative only of local conditions around the sampling point, an expected result given that ²²²Rn maximum effective diffusion length is very limited. As a consequence, the influence of a highly fluctuating water table at the site goes undetected at the sampling depths used in the field campaigns. MAIN FINDINGS: The combination of a numerical model and a laboratory protocol allows to predict the activity of ²²²Rn along the soil profile and to assess the influence of site-specific confounding factors.

Modeling of soil gas radon as an in situ partitioning tracer for quantifying LNAPL contamination

In the last decades radon (Rn) has been widely proposed as a naturally occurring tracer for non-aqueous phase liquids (NAPL) in the soil. This work examines the feasibility of using soil gas data collected at some distance from the source zone for the application of the Rn deficit technique for the identification and quantification of NAPL contamination. To this end, we used a steady-state 1-D analytical solution that is based on a 3-layer model that allows to simulate the transport and distribution of Rn in the source zone, capillary fringe and overlying unsaturated soil. The analytical solution was first validated against a more detailed numerical model available in the literature. Then, a series of simulations were carried out to evaluate the vertical concentration profiles of Rn in soil gas above the source zone and in background location not impacted by NAPL. Simulation results showed that the parameters that most influence the migration and distribution of Rn in the subsurface are the distance of the soil gas probe from the source zone and, to a lower extent, the type of contamination (e.g. diesel or gasoline) and the soil type. On the basis of these results, we developed some easy-to-use nomographs to estimate the residual NAPL phase based on the observed radon deficit in soil gas and on the probe to source distance and soil and NAPL characteristics. According to the obtained results, the radon deficit technique results a feasible method for a qualitative identification of residual NAPL when radon in soil gas is measured at distances lower than 2 m from the contaminated zone. However, for an accurate quantitative estimation of the NAPL phase content, soil gas probes should be preferably located at distances lower than 1 m from the source zone.

Effect of NAPL mixture and alteration on 222Rn partitioning coefficients: Implications for NAPL subsurface contamination quantification

Soils and groundwater are often contaminated by complex organic mixtures also called Non Aqueous Phase Liquids (NAPLs). Several techniques such as drilling, monitoring of soil gas or injection of tracers are traditionally used to quantify NAPLs in aquifers but are complex to perform. The use of natural soil gas such as ²²²Rn could be an easy and cheap alternative. This method requires the knowledge of the radon NAPL-water partitioning coefficients (K_n-w.). Once spilled on soil, NAPL will undergo degradation (evaporation, effects of sun light among others) and this degradation could impact the K_n-w. This study aims at investigating the partitioning coefficients of complex NAPLs such as commercial diesel fuel and gasoline in relation to degradation such as evaporation and UV-degradation. For that purpose, batch experiments and GCMS investigations were carried out. The results show different K_n-w for the commercial diesel fuel (60.7 ± 6.1) and gasoline (37.4 ± 5.6). The results also show different K_n-w behaviors in relation with degradation. Degraded diesel fuel display opposite K_n-w values (74.8 ± 7.5 and 25.1 ± 2.5 for UV degraded and evaporated diesel fuel, respectively), compared to fresh one. Degraded gasoline shows no significant variations of the K_n-w compared to fresh one. The molecular investigation reveals the removal of the most volatile fraction for the evaporation treatment, whereas UV-degradation do not have pronounced effects on the chromatogram pattern. For the gasoline molecular investigation, no difference is observed between the treatments excepted a very slight removal of the lightest compounds under evaporation. These results show that NAPL degradation have effects on the K_n-w for diesel fuel and no significant effects for gasoline, at least with these degradation paths. This K_n-w variation will have in fine effects on ²²²Rn activity interpretation and NAPL subsurface quantification.

Field performance of the radon-deficit technique to detect and delineate a complex DNAPL accumulation in a multi-layer soil profile

The performance of the radon (²²²Rn)-deficit technique has been evaluated at a site in which a complex DNAPL mixture (mostly hexachlorocyclohexanes and chlorobenzenes) has contaminated all four layers (from top to bottom: anthropic backfill, silt, gravel and marl) of the soil profile. Soil gas samples were collected at two depths (0.8 m and 1.7 m) in seven field campaigns and a total of 186 ²²²Rn measurements were performed with a pulse ionization detector. A statistical assessment of the influence of field parameters on the results revealed that sampling depth and atmospheric pressure did not significantly affect the measurements, while the location of the sampling point and ground-level atmospheric temperature did. In order to remove the bias introduced by varying field temperatures and hence to be able to jointly interpret ²²²Rn measurements from different campaigns, ²²²Rn concentrations were rescaled by dividing each individual datum by the mean ²²²Rn concentration of its corresponding field campaign. Rescaled ²²²Rn maps showed a high spatial correlation between ²²²Rn minima and maximum contaminant concentrations in the top two layers of the soil profile, successfully delineating the surface trace of DNAPL accumulation in the anthropic backfill and silt layers. However, no correlation could be established between ²²²Rn concentrations in superficial soil gas and contaminant concentration in the deeper two layers of the soil profile. These results indicate that the ²²²Rn-deficit technique is unable to describe the vertical variation of contamination processes with depth but can be an effective tool for the preliminary characterization of sites in which the distance between the inlet point of the sampling probe and the contaminant accumulation falls within the effective diffusion length of ²²²Rn in the affected soil profile.

A numerical model to optimize LNAPL remediation by multi-phase extraction

Light non-aqueous phase liquid (LNAPL) contaminated sites pose a risk to human health and the natural environment. Multi-phase extraction (MPE) is one of the most widely used technologies to remediate these sites. Thus, it is important to optimize MPE systems to improve their effectiveness and cost-efficiency. In this study, we developed a numerical model to optimize LNAPL mass removal by MPE, in which the aquifer domain was simplified as a cylinder with a single MPE extraction well located at the center. A dual-pump extraction system was applied to the model, which involved vacuum enhanced recovery to remove volatilized gaseous phase contaminants and a submerged pump to remove NAPL and contaminants in groundwater. After the model was validated with field data, the results showed that the contaminant extraction rate varied with the LNAPL thickness and submerged pump position. For benzene selected as the contaminant of concern, greater LNAPL extraction rates were achieved when the initial LNAPL thickness was large (>1.5 m) or in aquifers of high permeability (>2.45 × 10^-10 m²). Importantly, it was discovered that in highly permeable coarse sand and gravel, the submerged pump ought to be placed within the LNAPL layer, whereas the pump should be placed below the water-NAPL interface in fine to medium sand aquifers. It was also found that an optimal liquid pumping rates exist, beyond which contaminant mass removal rates do not increase. Furthermore, it was found that in aquifers contaminated with thin LNAPL layers, mass transfer modelling that assumes equilibrium between the phases may greatly overestimate the accumulated mass of contaminants removed and, therefore, non-equilibrium modelling should be adopted. Finally, a cost analysis was carried out to compare the costs of remediating a contaminated site with MPE and by an alternative chemical oxidation approach. The MPE technology was found to be more cost effective when the initial thickness of LNAPL was relatively thin. In summary, the numerical model developed in this study is a useful tool for optimizing MPE system design.

Applicability and limitations of the radon-deficit technique for the preliminary assessment of sites contaminated with complex mixtures of organic chemicals: A blind field-test

A blind field test with 136 independent measurements of radon (²²²Rn) in soil air retrieved from a depth of 0.8 m in a decommissioned lindane (γ-hexachlorocyclohexane) production plant was undertaken to evaluate the performance of the ²²²Rn-deficit technique as a screening methodology for the location and delineation of subsurface accumulations of complex mixtures of organic contaminants. Maps of ²²²Rn iso-concentrations were drawn and interpreted before direct analytical information regarding concentrations of hexachlorocyclohexanes, chlorobenzenes and BTEX compounds in soil, groundwater and soil air were disclosed to the authors. The location and extension of pollution hot spots inferred from the ²²²Rn campaigns agrees remarkably well with the analytical data obtained from the intrusive sampling campaigns and with the location of contaminant source zones (chemical reactor and waste-storage area) and geological sinks of those contaminants (paleochannel). Two main limitations to the applicability of the ²²²Rn-deficit technique were identified and assessed: The statistically significant variation of ²²²Rn concentrations with diurnal changes of ground-level air temperature and the maximum depth of investigation in the absence of significant advective and co-advective transport of radon. If the influence of those two factors is accounted for and/or minimized (by averaging replicated measurements during the workday and in different days), the ²²²Rn-deficit technique has the potential to be an efficient technique which delivers information in quasi-real time, with a much higher spatial density than that of intrusive techniques, at a much faster rate and at a significantly lower cost. MAIN FINDINGS: The ²²²Rn-deficit technique is an effective tool for real-time site characterization only limited by diffusion length of radon and diurnal temperature variations.