Field-scale multi-phase LNAPL remediation: Validating a new computational framework against sequential field pilot trials

Advective and diffusive gas phase transport in vadose zones: Importance for defining vapour risks and natural source zone depletion of petroleum hydrocarbons

Quantifying the interlinked behaviour of the soil microbiome, fluid flow, multi-component transport and partitioning, and biodegradation is key to characterising vapour risks and natural source zone depletion (NSZD) of light non-aqueous phase liquid (LNAPL) petroleum hydrocarbons. Critical to vapour transport and NSZD is transport of gases through the vadose zone (oxygen from the atmosphere, volatile organic compounds (VOCs), methane and carbon dioxide from the zone of LNAPL biodegradation). Volatilisation of VOCs from LNAPL, aerobic biodegradation, methanogenesis and heat production all generate gas pressure changes that may lead to enhanced gas fluxes apart from diffusion. Despite the importance of the gaseous phase dynamics in the vadose zone processes, the relative pressure changes and consequent scales of advective (buoyancy and pressure driven) / diffusive transport is less studied. We use a validated multi-phase multi-component non-isothermal modelling framework to differentiate gas transport mechanisms. We simulate a multicomponent unweathered gasoline LNAPL with high VOC content to maximise the potential for pressure changes due to volatilisation and to enable the joint effects of methanogenesis and shallower aerobic biodegradation of vapours to be assessed, along with heat production. Considering a uniform fine sand profile with LNAPL resident in the water table capillary zone, results suggest that biodegradation plays the key role in gas phase formation and consequent pressure build-up. Results suggest that advection is the main transport mechanism over a thin zone inside the LNAPL/capillary region, where the effective gaseous diffusion is very low. In the bulk of the vadose zone above the LNAPL region, the pressure change is minimal, and gaseous diffusion is dominant. Even for high biodegradation rate cases, pressure build-up due to heat generation (inducing buoyancy effects) is smaller than the contribution of gas formation due to biodegradation. The findings are critical to support broader assumptions of diffusive transport being dominant in vapour transport and NSZD assessments.

On quantifying global carbon emission from oil contaminated lands over centuries

Petroleum releases into the subsurface contribute to global soil carbon emissions. Quantifying releases and changes in releases of carbon from soils over the lifetime of a spill is complex. Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) embodies all key mechanisms for transformation to carbon gases and their release from soils including partitioning, transport and degradation of petroleum components. Quantification of the interconnected behaviours of the soil microbiome, fluid flow, multi-component transport, partitioning, and biodegradation is crucial for understanding NSZD. Volatilization from LNAPL, aerobic biodegradation, methanogenesis, and heat production all lead to release of greenhouse gases to the atmosphere. To estimate carbon emissions, using a validated computational platform, we modelled the long term NSZD of four petroleum hydrocarbon types; crude oil, diesel, jet fuel and gasoline, to span the major products used globally. For two soil types, we estimated 150 years of carbon emissions from annual minor and 25 mostly major petroleum hydrocarbon land release incidents since 1950 - with an estimated released mass of ~9 million tonnes across the circumstances considered. Up to 2100 the mass of carbon emitted to the atmosphere is estimated to range from 4 to 6 Teragrams, with nearly 60 % currently released. Nomographs generated help predict the fate of LNAPL plumes and carbon emissions due to NSZD, which is crucially important to management of soil and groundwater contamination. The method provides a basis to include additionally identified and future petroleum releases. It is noted that the petroleum mixture composition, degradation rates, volatilization, and subsurface characteristics all can influence carbon emission estimations.

LNAPL migration processes based on time-lapse electrical resistivity tomography

Contamination from light non-aqueous phase liquids (LNAPLs) and their derivatives, arising from exploration, production, and transportation, has become a prevalent pollution source. This poses direct threats to human health. However, conventional investigative methods face limitations when applied to studying the extent and migration process of LNAPL contamination, as well as the redistribution of LNAPL during groundwater level fluctuations. Conventional methods lack the ability to rapidly, efficiently, and in real-time acquire information about contaminated areas. Therefore, this study utilizes time-lapse electrical resistivity tomography to investigate the migration mechanism of LNAPL under unsaturated conditions, constant groundwater levels, and groundwater level reductions. A relationship between resistivity and water and oil contents was established and used for inverse calculation of LNAPL content via resistivity inversion. Time-lapse electrical resistivity tomography revealed LNAPL migration in a "concave" shape across three conditions. Groundwater presence notably slowed migration, hindering downward movement and leading to a floating oil band. A robust mathematical model was established to derive the relationship between resistivity and water and oil contents. Finally, LNAPL distribution under unsaturated conditions was inversely obtained from resistivity data, showing highest content at the top leak point, obstructed area, and bottom of soil column. Consequently, time-lapse electrical resistivity tomography demonstrates a notable capacity to characterize the LNAPL migration process. This technique constitutes an effective geophysical method for monitoring and describing the characteristics of LNAPL migration. Its significance lies in enhancing our understanding of remediation for LNAPL-induced groundwater and land contamination.

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.

The role of predicted chemotactic and hydrocarbon degrading taxa in natural source zone depletion at a legacy petroleum hydrocarbon site

Petroleum hydrocarbon contamination is a global problem which can cause long-term environmental damage and impacts water security. Natural source zone depletion (NSZD) is the natural degradation of such contaminants. Chemotaxis is an aspect of NSZD which is not fully understood, but one that grants microorganisms the ability to alter their motion in response to a chemical concentration gradient potentially enhancing petroleum NSZD mass removal rates. This study investigates the distribution of potentially chemotactic and hydrocarbon degrading microbes (CD) across the water table of a legacy petroleum hydrocarbon site near Perth, Western Australia in areas impacted by crude oil, diesel and jet fuel. Core samples were recovered and analysed for hydrocarbon contamination using gas chromatography. Predictive metagenomic profiling was undertaken to infer functionality using a combination of 16 S rRNA sequencing and PICRUSt2 analysis. Naphthalene contamination was found to significantly increase the occurrence of potential CD microbes, including members of the Comamonadaceae and Geobacteraceae families, which may enhance NSZD. Further work to explore and define this link is important for reliable estimation of biodegradation of petroleum hydrocarbon fuels. Furthermore, the outcomes suggest that the chemotactic parameter within existing NSZD models should be reviewed to accommodate CD accumulation in areas of naphthalene contamination, thereby providing a more accurate quantification of risk from petroleum impacts in subsurface environments, and the scale of risk mitigation due to NSZD.

Towards a digital twin for characterising natural source zone depletion: A feasibility study based on the Bemidji site

Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) may be a valid long-term management option at petroleum impacted sites. However, its future long-term reliability needs to be established. NSZD includes partitioning, biotic and abiotic degradation of LNAPL components plus multiphase fluid dynamics in the subsurface. Over time, LNAPL components are depleted and those partitioning to various phases change, as do those available for biodegradation. To accommodate these processes and predict trends and NSZD over decades to centuries, for the first time, we incorporated a multi-phase multi-component multi-microbe non-isothermal approach to representatively simulate NSZD at field scale. To validate the approach we successfully mimic data from the LNAPL release at the Bemidji site. We simulate the entire depth of saturated and unsaturated zones over the 27 years of post-release measurements. The study progresses the idea of creating a generic digital twin of NSZD processes and future trends. Outcomes show the feasibility and affordability of such detailed computational approaches to improve decision-making for site management and restoration strategies. The study provided a basis to progress a computational digital twin for complex subsurface systems.

Effect of water-level fluctuation on the removal of benzene from soil by SVE

In areas that experience frequent fluctuations in groundwater levels, changes in soil void space caused by fluctuations in groundwater levels directly affect the efficacy of soil vapor extraction (SVE). Most studies of SVE technology of oil-contaminated soil remediation do not consider fluctuations in groundwater levels. Here, we constructed an experimental device to study the removal of benzene by SVE under fluctuating groundwater levels. Key parameters affecting the remediation effect of SVE, such as the extraction flow, extraction time, extraction method, initial soil moisture content and initial pollutant content, were studied to characterize their effects on the efficacy of benzene removal by SVE under stable and fluctuating groundwater levels. The removal rate of benzene by SVE was approximately 10% higher under fluctuating water levels than under stable water levels. Extraction flow can directly change the removal rate under fluctuating groundwater levels. Under fluctuating groundwater levels, the removal effect of SVE on benzene in different soils of the "stabilization zone - fluctuation zone - saturation zone" also significantly differed; specifically, the removal effect of SVE in the fluctuation zone and stabilization zone was superior to that in the saturation zone. Given the fluctuation in groundwater levels at the study site, the remediation.

Quantifying the benefits of in-time and in-place responses to remediate acute LNAPL release incidents

Acute large volume spills from storage tanks of petroleum hydrocarbons as light non aqueous phase liquids (LNAPLs) can contaminate soil and groundwater and may have the potential to pose explosive and other risks. In consideration of an acute LNAPL release scenario, we explore the value of a rapid remediation response, and the value of installing remediation infrastructure in close proximity to the spill location, in effecting greater recovery of LNAPL mass from the subsurface. For the first time, a verified three-dimensional multi-phase numerical framework and supercomputing resources was applied to explore the significance of in-time and in-place remediation actions. A sand aquifer, two release volumes and a low viscosity LNAPL were considered in key scenarios. The time of commencement of LNAPL remediation activities and the location of recovery wells were assessed requiring asymmetric computational considerations. The volume of LNAPL released considerably affected the depth of LNAPL penetration below the groundwater table, the radius of the plume over time and the recoverable LNAPL mass. The remediation efficiency was almost linearly correlated with the commencement time, but was a non-linear function of the distance of an extraction well from the spill release point. The ratio of the recovered LNAPL in a well located at the centre of the spill/release compared to a well located 5 m away was more than 3.5, for recovery starting only 7 days after the release. Early commencement of remediation with a recovery well located at the centre of the plume was estimated to recover 190 times more LNAPL mass than a one-month delayed commencement through a well 15 m away from the centre of the LNAPL plume. Optimally, nearly 40% of the initially released LNAPL could be recovered within two months of commencing LNAPL recovery actions.

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.

Three-dimensional migration and resistivity characteristics of crude oil in heterogeneous soil layers

An experimentally induced three-dimensional petroleum seepage flume was used to investigate its migration in heterogeneous soil layers and a method for monitoring resistivity was adopted, under conditions of fluctuating water levels and rainfall. The corresponding mechanisms were then analyzed based on the resistivity characteristics and combined with three-dimensional inversion images. Finally, physical and chemical property analysis was conducted to verify the results of resistivity monitoring. The results demonstrated that: (1) In the process of natural oil leakage, the variation of soil resistivity presents a concave shape in the resistivity profile. Thus, oil migration exhibited the following patterns. At first, circular migration front was dominant in a vertical direction. Subsequently, after vertical migration was impeded, lateral migration was dominant. As the crude oil gradually accumulated, the migration front broke through the limitation of lithologic interface and continued vertically. (2) By comparing the two resistivity monitoring methods, namely the Wenner and Pole-pole methods, it was demonstrated that the inversion resistivity measured by Wenner method was closer to the true resistivity, and the resistivity variations were more distinguishable. (3) The resistivity inversion profile demonstrated that the low resistivity anomaly of the crude oil leakage area was related to the low water content of the soil layer in the test area. (4) Fluctuations in water level increased the diffusion range of crude oil beyond the original pollution source area, especially horizontally. (5) Percolation of rainfall caused the water level to rise, and the crude oil was evenly distributed in the soil layers above the capillary zone. (6) Through sample analysis and verification, it was demonstrated that the resistivity method can accurately and intuitively present the characteristics of crude oil migration. These results provide theoretical support for the rapid determination of the migration range and characteristics of crude oil in heterogeneous soil layers.

Pore-by-pore modeling, analysis, and prediction of two-phase flow in mixed-wet rocks

A pore-network model is an upscaled representation of the pore space and fluid displacement, which is used to simulate two-phase flow through porous media. We use the results of pore-scale imaging experiments to calibrate and validate our simulations, and specifically to find the pore-scale distribution of wettability. We employ energy balance to estimate an average, thermodynamic, contact angle in the model, which is used as the initial estimate of contact angle. We then adjust the contact angle of each pore to match the observed fluid configurations in the experiment as a nonlinear inverse problem. The proposed algorithm is implemented on two sets of steady state micro-computed-tomography experiments for water-wet and mixed-wet Bentheimer sandstone. As a result of the optimization, the pore-by-pore error between the model and experiment is decreased to less than that observed between repeat experiments on the same rock sample. After calibration and matching, the model predictions for capillary pressure and relative permeability are in good agreement with the experiments. The proposed algorithm leads to a distribution of contact angle around the thermodynamic contact angle. We show that the contact angle is spatially correlated over around 4 pore lengths, while larger pores tend to be more oil-wet. Using randomly assigned distributions of contact angle in the model results in poor predictions of relative permeability and capillary pressure, particularly for the mixed-wet case.

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.

Impact of mineralogy and wettability on pore-scale displacement of NAPLs in heterogeneous porous media

Subsurface formations often contain multiple minerals with different wettability characteristics upon contact with nonaqueous-phase liquids (NAPLs). Constitutive relationships between microstructure heterogeneity and NAPL fate and transport in these formations are difficult to predict. Several studies have used pore-scale network models with faithful representations of rock pore space topology to predict macroscopic descriptors of two-phase flow, however wettability is usually considered as a spatially random variable. This study attempts to overcome this limitation by considering more realistic representations of rock mineralogy and wettability in these models. This is especially important for heterogeneous rocks where properties vary at the pore-scale. The work was carried out in two phases. First, pore-fluid occupancy maps during waterflooding were obtained by X-ray microtomography to elucidate the impact of pore wall mineralogy and wettability on water preferential flow paths and NAPL trapping within a heterogeneous aquifer sandstone (Arkose). Then, microtomography images of the rock were used to generate a hybrid pore network model (PNM) that incorporated both pore space topology and pore wall mineralogy. In-situ contact angles (CA) measured on the surface of different minerals were assigned to the network on a pore-by-pore basis to describe the exact wettability distribution of the rock (Pore-by-pore model). The equivalent network was used as input in a quasi-static flow model to simulate waterflooding, and the predictions of residual NAPL saturation and relative permeabilities were compared against their experimental counterparts. To examine the sensitivity of the model to the underlying fluid-solid interactions, we also used traditional methods of wettability characterization in the input data and assigned them randomly to the PNM. Wettability in this case was assessed from macroscale CA distribution of oil droplets on the surface of unpolished Arkose substrates released by spontaneous imbibition of water (Arkose model) and from pendant drop measurements on polished quartz (Quartz model). Our results revealed that the Pore-by-pore model predicted waterflooding with the highest accuracy among all three cases. The Arkose model slightly overestimated NAPL removal whereas the Quartz model failed to predict the experiments. More in-depth analysis of the Pore-by-pore and Arkose models showed that macroscopic transport quantities are less dependent to microstructure heterogeneity if minerals are distributed uniformly across the rock. The predictions herein indicate the importance of incorporating mineralogy and wettability maps to improve the prediction capabilities of PNMs especially in systems with high mineral heterogeneity, where minerals are nonuniformly distributed, or selective fluid-mineral interactions are targeted.

Numerical Modeling of Multiphase Extraction (MPE) Aiming at LNAPL Recovery in Tropical Soils

Subsurface contamination by light non-aqueous phase liquids (LNAPL) is a widespread global problem that requires appropriate techniques to remediate soil and groundwater. In this paper, the subsurface transport over multiple phases (STOMP) model was used to simulate LNAPL multiphase flow and transport during multiphase extraction (MPE) application in two Brazilian tropical soils (silty sand and oxisol) contaminated by diesel. The model was applied to a hypothetical contamination site, with the initial LNAPL thickness observed in well extraction. The first part consisted of the MPE system sensitivity analysis, varying the applied vacuum and tip tube position. The Van Genuchten α parameter and hydraulic conductivity were the properties that most affected LNAPL saturation and fluid extraction volumes. Suitable applied vacuum and tip tube position parametrization was imperative for the efficiency of LNAPL extraction. After the definition of an appropriate MPE system configuration, simulations demonstrated that the immobile LNAPL saturation affected fluid extraction and diesel oil concentrations in aqueous and gas saturation. The model applied is able to predict LNAPL contaminant behavior in porous media during MPE technique application.

Sustainable risk-based analysis towards remediation of an aquifer impacted by crude oil spills

Uncontrolled release of hydrocarbons from pipelines results in soil and groundwater contamination. However, due to the geo-environmental properties of the contaminated area, the remediation strategies might vary by light non-aqueous phase liquid (LNAPL) behaviors. In this study, a contaminated area with spilled oil from a pipeline was monitored. In the initial investigation in the 1980s, the contamination was reported in some citizen wells (CIZs), which resulted in drilling of 15 boreholes (BHs) across the Site from February 2001 to February 2007. Additionally, an area of around 0.24 km² encompassing the LNAPL was inferred. The extension of plume in 2016-7 was monitored in this study, in which new farms were impacted by the contamination. Further, a conceptual model based on the previous information and current measured data was developed to better understand the behavior of the plume. The model showed that the Site is very complex, dipping towards the south, and the groundwater contains light hydrocarbons. Pumping tests, as a part of LNAPL remedial technology, were conducted by using three pumping wells (PWs), each accompanied by a monitoring well. Accordingly, a risk-based corrective action was implemented to eliminate and control unacceptable risks in a safe and timely manner. From the remediation approach, a monitoring plan in BHs and CIZs was suggested. In the case of receptors (humans and farms), clean-up of wells, tanks, and water channels as well as replacement of contaminated soils were highly regarded. Although the recent investigation and clean up monitoring wells showed that the LNAPL was very minimal, further steps in the receptor side should be taken prior to irrigation applications.

An approach to optimize the location of LNAPL recovery wells using the concept of a LNAPL specific yield

Leakage of hydrocarbon fuel (light nonaqueous-phase liquid, LNAPL) from petroleum processing facilities and storage tanks may result in significant subsurface contamination. Remediating the contaminated areas represent considerable challenges, especially when remediation resources are limited and site data are incomplete. A reasonable management strategy under this scenario may be to identify sites where LNAPL recovery operations should be located that would provide the largest LNAPL recovery initially while minimizing the LNAPL remaining in the subsurface (entrapped and residual LNAPL), which may serve as future sources for groundwater contamination. To accomplish this objective, we use estimates of subsurface recoverable and total LNAPL specific volumes and LNAPL transmissivities to generate GIS maps that can be combined to highlight locations where to develop LNAPL recovery operations. When the approach is applied to a LNAPL-contaminated area in Iran, we were able to narrow the locations for potential LNAPL recovery operations. Specifically, we combine maps of the LNAPL specific yield, an introduced term, and the LNAPL transmissivity where the LNAPL specific yield is the ratio of the recoverable to total LNAPL specific volumes. The LNAPL specific yield is a relative measure of the amount of LNAPL that potentially can be recovered while minimizing residual LNAPL in soils. The approach can be applied to sites where the recoverable and total LNAPL specific volumes and LNAPL transmissivities can be estimated using data from boreholes in the contaminated area.

Analysis of radially convergent tracer test in a two-zone confined aquifer with vertical dispersion effect: Asymmetrical and symmetrical transports

Radially convergent tracer test (RCT) with an extraction well and a tracer injection well is commonly conducted for determining aquifer dispersivity associated with the spreading of contamination plume. A variety of analytical models for RCT have been proposed, but the effects of skin zone and vertical dispersion are rarely explored. In addition, little attention is paid to the validity of those analytical models assuming asymmetrical tracer transport (ATT) as radially symmetrical transport (RST) toward the pumping well. This study develops a new analytical model for RCT subject to the above-mentioned effects in a confined aquifer with a skin zone around the wellbore. The Laplace-domain solution of the model for a continuous or instantaneous input is developed. A finite element solution (FES) for ATT is also developed to verify the Laplace-domain solution based on RST. Results suggest ATT can be regarded as RST in predicting the breakthrough curve (BTC) at the pumping well when four quantitative conditions are met. A lumped dimensionless parameter dominates those two effects on the BTC. Both the FES and Laplace-domain solution agree with monitored concentration data from a field RCT.

Natural source zone depletion of LNAPL: A critical review supporting modelling approaches

Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) includes partitioning, transport and degradation of LNAPL components. NSZD is being considered as a site closure option during later stages of active remediation of LNAPL contaminated sites, and where LNAPL mass removal is limiting. To ensure NSZD meets compliance criteria and to design enhanced NSZD actions if required, residual risks posed by LNAPL and its long term behaviour require estimation. Prediction of long-term NSZD trends requires linking physicochemical partitioning and transport processes with bioprocesses at multiple scales within a modelling framework. Here we expand and build on the knowledge base of a recent review of NSZD, to establish the key processes and understanding required to model NSZD long term. We describe key challenges to our understanding, inclusive of the dominance of methanogenic or aerobic biodegradation processes, the potentially changeability of rates due to the weathering profile of LNAPL product types and ages, and linkages to underlying bioprocesses. We critically discuss different scales in subsurface simulation and modelling of NSZD. Focusing on processes at Darcy scale, 36 models addressing processes of importance to NSZD are investigated. We investigate the capabilities of models to accommodate more than 20 subsurface transport and transformation phenomena and present comparisons in several tables. We discuss the applicability of each group of models for specific site conditions.

Towards characterizing LNAPL remediation endpoints

Remediating sites contaminated with light non-aqueous phase liquids (LNAPLs) is a demanding and often prolonged task. It is vital to determine when it is appropriate to cease engineered remedial efforts based on the long-term effectiveness of remediation technology options. For the first time, the long term effectiveness of a range of LNAPL remediation approaches including skimming and vacuum-enhanced skimming each with and without water table drawdown was simulated through a multi-phase and multi-component approach. LNAPL components of gasoline were simulated to show how component changes affect the LNAPL's multi-phase behaviour and to inform the risk profile of the LNAPL. The four remediation approaches along with five types of soils, two states of the LNAPL specific mass and finite and infinite LNAPL plumes resulted in 80 simulation scenarios. Effective conservative mass removal endpoints for all the simulations were determined. As a key driver of risk, the persistence and mass removal of benzene was investigated across the scenarios. The time to effectively achieve a technology endpoint varied from 2 to 6 years. The recovered LNAPL in the liquid phase varied from 5% to 53% of the initial mass. The recovered LNAPL mass as extracted vapour was also quantified. Additional mass loss through induced biodegradation was not determined. Across numerous field conditions and release incidents, graphical outcomes provide conservative (i.e. more prolonged or greater mass recovery potential) LNAPL remediation endpoints for use in discussing the halting or continuance of engineered remedial efforts.

LNAPL transmissivity as a remediation metric in complex sites under water table fluctuations

Water table fluctuations affect the recoverability of light non-aqueous phase liquid (LNAPL) petroleum hydrocarbons. LNAPL transmissivity (T_n) is being applied as an improved metric for LNAPL recoverability. In this paper, the applicability of T_n as a lagging and leading metric in unconsolidated aquifers under variable water table conditions was investigated. T_n values obtained through baildown testing and recovery data-based methods (skimming) were compared in three areas of a heterogeneous gasoline contaminated site in Western Australia. High-resolution characterisation methods were applied to account for differences in the stratigraphic profile and LNAPL distribution. The results showed a range of T_n from 0 m²/day to 2.13 m²/day, exhibiting a strong spatial and temporal variability. Additionally, observations indicated that T_n reductions may be more affected by the potentiometric surface elevation (Z_aw) than by the application of mass recovery technologies. These observations reflected limitations of T_n as a lagging metric and a Remedial Endpoint. On the other hand, the consistency and accuracy of T_n as a leading metric was affected by the subsurface conditions. For instance, the area with a larger vertical LNAPL distribution and higher LNAPL saturations found T_n to be less sensitive to changes in Z_aw than the other two areas during the skimming trials. T_n values from baildown and skimming tests were generally in a close agreement (less than a factor of 2 difference), although higher discrepancies (by a factor up to 7.3) were found, probably linked to a preferential migration pathway and Z_aw. Under stable Z_aw, T_n was found to be a relatively reliable metric. However, variable water table conditions affected T_n and caution should be exercised in such scenarios. Consequently, remediation practitioners, researchers and regulators should account for the nexus between T_n, LNAPL distribution, geological setting and temporal effects for a more efficient and sustainable management of complex contaminated sites.