Dissolution of multi-component LNAPL gasolines: the effects of weathering and composition

Natural Source Zone Depletion of crude oil, gasoline, diesel and aviation gasoline petroleum in the same soil/aquifer system - An intensive intercomparison field study and simulation

Petroleum biodegrades and naturally depletes. Natural Source Zone Depletion (NSZD) quantifies this at petroleum affected sites in support of management decisions for cessation of active remediation efforts. Whilst a range of NSZD estimates and methods are available, side by side comparison of NSZD rates across petroleum types in the same soil/groundwater system are lacking, especially linked to the weathering status of petroleum. At a former refinery site near Perth Western Australia, locations contaminated by crude oil, gasoline, diesel and aviation gasoline, have been intensively instrumented to enable (i) measurement of vadose zone major gas (O2, CO2, CH4), volatile organic compounds (VOCs) and temperature depth profiles, (ii) online near-continuous temperature, water level, O2 and VOC concentrations, (iii) depth profiles of groundwater parameters, and (iv) in-well gases, temperatures, water and LNAPL thickness. These measurements were compared to a background location with no history of contamination. Multiple coring events were also undertaken to determine LNAPL mass and its vertical distribution at each location. Additionally, LiCor and Eflux was conducted to measure CO2 fluxes at ground surface. NSZD rates were estimated from the measurement methods across the four petroleum types and the background site. Despite NSZD estimates that differed across some methods at sites (for example at the gasoline and diesel sites Eflux/LiCor estimates were consistently lower by a factor of 3-4 than those obtained using oxygen/temperature depth profile data) the minimum-maximum range of mean NSZD rates showed a distinct decreasing order across the fuel types: highest rates being aviation gasoline (69,000-91,000 L/ha/y), then gasoline, diesel and crude oil (2700-6200 L/ha/y). Reasons for differences are explored. Analysis of LNAPL in cores and from wells, historical data comparisons and simulations over 50 years, show that composition, age and weathering of the releases are critical to current and long-term NSZD mass losses and rate estimates.

An experimental multi-method approach to better characterize the LNAPL fate in soil under fluctuating groundwater levels

Light-Non-Aqueous phase liquids (LNAPLs) are important soil contamination sources, and groundwater fluctuations may significantly affect their migration and release. However, the risk assessment remains complex due to the continuous three-phase fluid redistribution caused by water table level variations. Hence, monitoring methods must be improved to integrate better the LNAPL multi-compound and multi-phase aspects tied to the groundwater level dynamics. For this purpose, a lysimetric contaminated soil column (2 m3) combining in-situ monitoring (electrical permittivity, soil moisture, temperature, pH, Eh), direct water and gas sampling and analyses (GC/MS-TQD, μGC) in monitoring well, gas collection chambers, and suction probes) were developed. This experiment assesses in an integrated way how controlled rainfalls and water table fluctuation patterns may affect LNAPL vertical soil saturation distribution and release. Coupling these methods permitted the investigation of the effects of rainwater infiltration and water table level fluctuation on contaminated soil oxygen turnover, LNAPL contaminants' soil distribution and remobilization towards the dissolved and the gaseous phase, and the estimate of the LNAPL source attenuation rate. Hence, 7.5% of the contamination was remobilized towards the dissolved and gaseous phase after 120 days. During the experiment, groundwater level variations were responsible for the free LNAPL soil spreading and trapping, modifying dissolved LNAPL concentrations. Nevertheless, part of the dissolved contamination was rapidly biodegraded, leaving only the most bio-resistant components in water. This result highlights the importance of developing new experimental devices designed to assess the effect of climate-related parameters on LNAPL fate at contaminated sites.

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.

Emulsion-based recovery of a multicomponent petroleum hydrocarbon NAPL using nonionic surfactant formulations

Surfactants can aid subsurface remediation through three primary mechanisms - solubilization, mobilization and/or emulsification. Among these mechanisms, emulsification in porous media is generally not well studied or well understood; particularly in the context of treating sources containing multicomponent NAPL. The objective of this research was to elucidate the processes responsible for recovery of a multicomponent hydrocarbon NAPL when surfactant solutions are introduced within a porous medium to promote the formation of kinetically-stable oil-in-water emulsions. Emulsifier formulations considered here were selected to offer similar performance characteristics while relying on different families of non-ionic surfactants - nonylphenol ethoxylates or alcohol ethoxylates - for emulsification. The families of surfactants have particular environment relevance, as alcohol ethoxylates are often used where replacement of nonylphenol content is necessary. Results from batch and column studies suggest performance of the two formulations was similar. With both, a synergistic combination of emulsification and mobilization led to recovery of a synthetic gasoline NAPL. The relative contribution of solubilization to the recovery was found to be minor. Moreover, the physical processes associated with emulsification and mobilization acted to limit the amount of preferential recovery (or fractionation) of the multicomponent NAPL.

A critical review of the influence of groundwater level fluctuations and temperature on LNAPL contaminations in the context of climate change

The intergovernmental panel on climate change (IPCC) predicts significant changes in precipitation patterns, an increase in temperature, and groundwater level variations by 2100. These changes are expected to alter light non-aqueous phase liquid (LNAPL) impacts since groundwater level fluctuations and temperature are known to influence both the mobility and release of LNAPL compounds to air and groundwater. Knowledge of these potential effects is currently dispersed in the literature, hindering a clear vision of the processes at play. This review aims to synthesize and discuss the possible effects of the increase in temperature and groundwater level fluctuations on the behavior of LNAPL and its components in a climate change context. In summary, a higher amplitude of groundwater table variations and higher temperatures will probably increase biodegradation processes, the LNAPL mobility, and spreading across the smear zone, favoring the release of LNAPL compounds to the atmosphere and groundwater but decreasing the LNAPL mass and its longevity. Outcomes will, nevertheless, vary greatly across arid, cold, or humid coastal environments, where different effects of climate change are expected. The effects of the climate change factors linked to soil heterogeneities, local conditions, and weathering processes will govern LNAPL behavior and need to be further clarified.

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.

Unravelling Microbial Communities Associated with Different Light Non-Aqueous Phase Liquid Types Undergoing Natural Source Zone Depletion Processes at a Legacy Petroleum Site

Petroleum contaminants are exposed to weathering when released into environment, resulting in the alteration of their chemical composition. Here, we investigated microbial communities through the soil profile at an industrial site, which was exposed to various petroleum products for over 50 years. The petroleum is present as light non-aqueous phase liquid (LNAPL) and is undergoing natural source zone depletion (NSZD). Microbial community composition was compared to the contaminant type, concentration, and its depth of obtained soil cores. A large population of Archaea, particularly Methanomicrobia and Methanobacteria and indication of complex syntrophic relationships of methanogens, methanotrophs and bacteria were found in the contaminated cores. Different families were enriched across the LNAPL types. Results indicate methanogenic or anoxic conditions in the deeper and highly contaminated sections of the soil cores investigated. The contaminant was highly weathered, likely resulting in the formation of recalcitrant polar compounds. This research provides insight into the microorganisms fundamentally associated with LNAPL, throughout a soil depth profile above and below the water table, undergoing NSZD processes at a legacy petroleum site. It advances the potential for integration of microbial community effects on bioremediation and in response to physicochemical partitioning of LNAPL components from different petroleum types.

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.

Magnetic resonance imaging of enhanced mobility of light non aqueous phase liquid (LNAPL) during drying of water wet porous media

Visualization of NAPLs in multiphase systems in porous media is important for determining contaminant transport in the environment. In this study, magnetic resonance imaging (MRI) was used to confirm the recent observations of mobilisation of a light non aqueous phase liquid (LNAPL) trapped in wet sand under natural drying conditions of the wet porous medium. Visualization of LNAPL (motor oil) and water mobility during the drying of wet glass beads (0.5 mm) in a cylindrical glass column (15 mm ID, 45 mm long) was obtained using spin echo-based NMR microimaging performed at 500 MHz, corresponding to a field of ca. 11.75 T. Sagittal and axial images of LNAPL and water in the porous medium were obtained at a spatial resolution of 59 μm/pixel at different time intervals. A rise of 15-20 mm was observed in the presence of evaporation of water as compared to a 2-3 mm rise in the absence of evaporation in a time span of about 1400 min. The spatio-temporal MRI scans of the water and LNAPL in the glass column reveals that LNAPL rise occurs when the water evaporation front reaches the LNAPL layer. This implied that the enhanced LNAPL rise was strongly linked to the process of water evaporation. A linear correlation of the MRI signal intensities of LNAPL and water with reference to different saturation levels of LNAPL and water in the porous media was obtained. This calibration information was used to quantify the saturation levels of the LNAPL and water during the drying process. These findings show the application of non-invasive techniques such as MRI in quantifying and understanding the mechanism of fate and transport of LNAPLs in porous media, towards effective environmental quality assessment.

Intra-NAPL diffusion and dissolution of a MGP NAPL exposed to persulfate in a flow-through system

Intra-NAPL diffusion is a critical process that can influence NAPL/water mass transfer. A series of physical model experiments was performed to investigate the role of intra-NAPL diffusion on the transient dissolution of a complex multicomponent NAPL subjected to persulfate treatment. To support these observations, a diffusion-based model was developed and calibrated using the experimental data. The experimental results indicated that while persulfate was able to completely degrade dissolved phase components, mass loss after ∼410 pore volumes of persulfate flushing was less than the no-treatment system. Intra-NAPL diffusion limitations were not observed in the physical model experiments. A comparison of experimental and simulated results indicated that processes related to persulfate/NAPL interactions restricted mass transfer, and yielded multicomponent mass transfer rate coefficients that were ∼30% of those estimated from an equivalent water-flushing experiment. Simulation results showed that a combination of NAPL composition and geometry, and interphase mass transfer rate can yield intra-NAPL diffusion limitations. Remedial technologies that rely on the aggressive flushing of reagents into NAPL zones may give rise to intra-NAPL diffusion limitations, which will directly affect treatment efficiency.

Weathering of Oil in a Surficial Aquifer

The composition of crude oil in a surficial aquifer was determined in two locations at the Bemidji, MN, spill site. The abundances of 71 individual hydrocarbons varied within 16 locations sampled. Little depletion of these hydrocarbons (relative to the pipeline oil) occurred in the first 10 years after the spill, whereas losses of 25% to 85% of the total measured hydrocarbons occurred after 30 years. The C_6-30 n-alkanes, toluene, and o-xylene were the most depleted hydrocarbons. Some hydrocarbons, such as the n-C_10-24 cyclohexanes, tri- and tetra- methylbenzenes, acyclic isoprenoids, and naphthalenes were the least depleted. Benzene was detected at every sampling location 30 years after the spill. Degradation of the oil led to increases in the percent organic carbon and in the δ ¹³ C of the oil. Another method of determining hydrocarbon loss was by normalizing the total measured hydrocarbon concentrations to that of the most conservative analytes. This method indicated that the total measured hydrocarbons were depleted by 47% to 77% and loss of the oil mass over 30 years was 18% to 31%. Differences in hydrocarbon depletion were related to the depth of the oil in the aquifer, local topography, amount of recharge reaching the oil, availability of electron acceptors, and the presence of less permeable soils above the oil. The results from this study indicate that once crude oil has been in the subsurface for a number of years there is no longer a "starting oil concentration" that can be used to understand processes that affect its fate and the transport of hydrocarbons in groundwater.

The role of intra-NAPL diffusion on mass transfer from MGP residuals

An experimental and computational study was performed to investigate the role of multi-component intra-NAPL diffusion on NAPL-water mass transfer. Molecular weight and the NAPL component concentrations were determined to be the most important parameters affecting intra-NAPL diffusion coefficients. Four NAPLs with different viscosities but the same quantified mass were simulated. For a spherical NAPL body, a combination of NAPL properties and interphase mass transfer rate can result in internal diffusion limitations. When the main intra-NAPL diffusion coefficients are in the range of self-diffusion coefficients (10^-5 to 10^-6 cm²/s), dissolution is not limited by internal diffusion except for high mass transfer rate coefficients (>180 cm/day). For a complex and relatively high viscous NAPL (>50 g/(cm s)), smaller intra-NAPL diffusion coefficients (<10^-8) are expected and even low mass transfer rate coefficients (~6 cm/day) can result in diffusion-limited dissolution.

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

Remediation of subsurface systems, including groundwater, soil and soil gas, contaminated with light non-aqueous phase liquids (LNAPLs) is challenging. Field-scale pilot trials of multi-phase remediation were undertaken at a site to determine the effectiveness of recovery options. Sequential LNAPL skimming and vacuum-enhanced skimming, with and without water table drawdown were trialled over 78days; in total extracting over 5m³ of LNAPL. For the first time, a multi-component simulation framework (including the multi-phase multi-component code TMVOC-MP and processing codes) was developed and applied to simulate the broad range of multi-phase remediation and recovery methods used in the field trials. This framework was validated against the sequential pilot trials by comparing predicted and measured LNAPL mass removal rates and compositional changes. The framework was tested on both a Cray supercomputer and a cluster. Simulations mimicked trends in LNAPL recovery rates (from 0.14 to 3mL/s) across all remediation techniques each operating over periods of 4-14days over the 78day trial. The code also approximated order of magnitude compositional changes of hazardous chemical concentrations in extracted gas during vacuum-enhanced recovery. The verified framework enables longer term prediction of the effectiveness of remediation approaches allowing better determination of remediation endpoints and long-term risks.

Enhanced mobility of non aqueous phase liquid (NAPL) during drying of wet sand

Enhanced upward mobility of a non aqueous phase liquid (NAPL) present in wet sand during natural drying, and in the absence of any external pressure gradients, is reported for the first time. This mobility was significantly higher than that expected from capillary rise. Experiments were performed in a glass column with a small layer of NAPL-saturated sand trapped between two layers of water-saturated sand. Drying of the wet sand was induced by flow of air across the top surface of the wet sand. The upward movement of the NAPL, in the direction of water transport, commenced when the drying effect reached the location of the NAPL and continued as long as there was significant water evaporation in the vicinity of NAPL, indicating a clear correlation between the NAPL rise and water evaporation. The magnitude and the rate of NAPL rise was measured at different water evaporation rates, different initial locations of the NAPL, different grain size of the sand and the type of NAPL (on the basis of different NAPL-glass contact angle, viscosity and density). A positive correlation was observed between average rate of NAPL rise and the water evaporation while a negative correlation was obtained between the average NAPL rise rate and the NAPL properties of contact angle, viscosity and density. There was no significant correlation of average NAPL rise rate with variation of sand grain size between 0.1 to 0.5mm. Based on these observations and on previous studies reported in the literature, two possible mechanisms are hypothesized -a) the effect of the spreading coefficient resulting in the wetting of NAPL on the water films created and b) a moving water film due to evaporation that "drags" the NAPL upwards. The NAPL rise reported in this paper has implications in fate and transport of chemicals in NAPL contaminated porous media such as soils and exposed dredged sediment material, which are subjected to varying water saturation levels due to drying and rewetting.

In situ stabilization of NAPL contaminant source-zones as a remediation technique to reduce mass discharge and flux to groundwater

Widely used flushing and in-situ destruction based remediation techniques (i.e. pump-and treat, enhanced-solubilization, and chemical oxidation/reduction) for sites contaminated by nonaqueous phase liquid (NAPL) contaminant sources have been shown to be ineffective at complete mass removal and reducing aqueous-phase contaminant of concern (COC) concentrations to levels suitable for site closure. A remediation method was developed to reduce the aqueous solubility and mass-flux of COCs within NAPL through the in-situ creation of a NAPL mixture source-zone. In contrast to remediation techniques that rely on the rapid removal of contaminant mass, this technique relies on the stabilization of difficult-to-access NAPL sources to reduce COC mass flux to groundwater. A specific amount (volume) of relatively insoluble n-hexadecane (HEXDEC) or vegetable oil (VO) was injected into a trichloroethene (TCE) contaminant source-zone through a bench-scale flow cell port (i.e. well) to form a NAPL mixture of targeted mole fraction (TCE:HEXDEC or TCE:VO). NAPL-aqueous phase batch tests were conducted prior to the flow-cell experiments to evaluate the effects of various NAPL mixture ratios on equilibrium aqueous-phase concentrations of TCE to design optimal NAPL (HEXDEC or VO) injection volumes for the flow-cell experiments. The NAPL-stabilization flow-cell experiments initiated and sustained significant reductions in COC concentration and mass flux due to a combination of both reduced relative permeability (increased NAPL-saturation) and via modification of NAPL composition (decreased TCE mole fraction). Variations in remediation performance (i.e. impacts on TCE concentration and mass flux reduction) between the different HEXDEC injection volumes were relatively minor, and therefore inconsistent with Raoult's Law predictions. This phenomenon likely resulted from non-uniform mixing of the injected HEXDEC with TCE in the source-zone. VO injection caused TCE concentrations and mass-flux to decrease more rapidly than with HEXDEC injections. This phenomenon occurred because the injected VO was observed to mix more uniformly with TCE in the source-zone due to a lower mobilization potential. The relative lower density differences (buoyancy effects) between VO and the flushing solution (water) was the primary factor contributing to the lower mobilization potential for VO. Overall, this study indicated that the delivery of HEXDEC or VO into the toxic TCE source-zone was effective in significantly reducing contaminant aqueous-phase concentration and mass-flux. However, the effectiveness of this in-situ NAPL stabilization technique depends on source delivery, uniform mixing of amendment, and that the amendment remains immobilized within and around the NAPL contaminant source.

Chemical structure influence on NAPL mixture nonideality evolution, rate-limited dissolution, and contaminant mass flux

The influence of chemical structure on NAPL mixture nonideality evolution, rate-limited dissolution, and contaminant mass flux was examined. The variability of measured and UNIFAC modeled NAPL activity coefficients as a function of mole fraction was compared for two NAPL mixtures containing structurally-different contaminants of concern including toluene (TOL) or trichloroethene (TCE) within a hexadecane (HEXDEC) matrix. The results showed that dissolution from the NAPL mixtures transitioned from ideality for mole fractions >0.05 to nonideality as mole fractions decreased. In particular, the TCE generally exhibited more ideal dissolution behavior except at lower mole fractions, and may indicate greater structural/polarity similarity between the two compounds. Raoult's Law and UNIFAC generally under-predicted the batch experiment results for TOL:HEXDEC mixtures especially for mole fractions ≤0.05. The dissolution rate coefficients were similar for both TOL and TCE over all mole fractions tested. Mass flux reduction (MFR) analysis showed that more efficient removal behavior occurred for TOL and TCE with larger mole fractions compared to the lower initial mole fraction mixtures (i.e. <0.2). However, compared to TOL, TCE generally exhibited more efficient removal behavior over all mole fractions tested and may have been the result of structural and molecular property differences between the compounds. Activity coefficient variability as a function of mole fraction was quantified through regression analysis and incorporated into dissolution modeling analyses for the dynamic flushing experiments. TOL elution concentrations were modeled (predicted) reasonable well using ideal and equilibrium assumptions, but the TCE elution concentrations could not be predicted using the ideal model. Rather, the dissolution modeling demonstrated that TCE elution was better described by the nonideal model whereby NAPL-phase activity coefficient varied as a function of COC mole fraction. For dynamic column flushing experiments, dissolution rate kinetics can vary significantly with changes in NAPL volume and surface area. However, under conditions whereby NAPL volume and area are not significantly altered during dissolution, mixture nonideality effects may have a greater relative control on dissolution (elution) and MFR behavior compared to kinetic rate limitations.

Field data and numerical simulation of btex concentration trends under water table fluctuations: Example of a jet fuel-contaminated site in Brazil

Mass transfer of light non-aqueous phase liquids (LNAPLs) trapped in porous media is a complex phenomenon. Water table fluctuations have been identified as responsible for generating significant variations in the concentration of dissolved hydrocarbons. Based on field evidence, this work presents a conceptual model and a numerical solution for mass transfer from entrapped LNAPL to groundwater controlled by both LNAPL saturation and seasonal water table fluctuations within the LNAPL smear zone. The numerical approach is capable of reproducing aqueous BTEX concentration trends under three different scenarios - water table fluctuating within smear zone, above the smear zone and partially within smear zone, resulting in in-phase, out-of-phase and alternating in-phase and out-of-phase BTEX concentration trend with respect to water table oscillation, respectively. The results demonstrate the model's applicability under observed field conditions and its ability to predict source zone depletion.

Evaluating the reliability of equilibrium dissolution assumption from residual gasoline in contact with water saturated sands

Understanding dissolution dynamics of hazardous compounds from complex gasoline mixtures is a key to long-term predictions of groundwater risks. The aim of this study was to investigate if the local equilibrium assumption for BTEX and TMBs (trimethylbenzenes) dissolution was valid under variable saturation in two dimensional flow conditions and evaluate the impact of local heterogeneities when equilibrium is verified at the scale of investigation. An initial residual gasoline saturation was established over the upper two-thirds of a water saturated sand pack. A constant horizontal pore velocity was maintained and water samples were recovered across 38 sampling ports over 141days. Inside the residual NAPL zone, BTEX and TMBs dissolution curves were in agreement with the TMVOC model based on the local equilibrium assumption. Results compared to previous numerical studies suggest the presence of small scale dissolution fingering created perpendicular to the horizontal dissolution front, mainly triggered by heterogeneities in the medium structure and the local NAPL residual saturation. In the transition zone, TMVOC was able to represent a range of behaviours exhibited by the data, confirming equilibrium or near-equilibrium dissolution at the scale of investigation. The model locally showed discrepancies with the most soluble compounds, i.e. benzene and toluene, due to local heterogeneities exhibiting that at lower scale flow bypassing and channelling may have occurred. In these conditions mass transfer rates were still high enough to fall under the equilibrium assumption in TMVOC at the scale of investigation. Comparisons with other models involving upscaled mass transfer rates demonstrated that such approximations with TMVOC could lead to overestimate BTEX dissolution rates and underestimate the total remediation time.

Effect of compositional heterogeneity on dissolution of non-ideal LNAPL mixtures

The extent of dissolution of petroleum hydrocarbon fuels into groundwater depends greatly on fuel composition. Petroleum fuels can consist of thousands of compounds creating different interactions within the non-aqueous phase liquid (NAPL), thereby affecting the relative dissolution of the components and hence a groundwater plume's composition over long periods. Laboratory experiments were conducted to study the variability in the effective solubilities and activity coefficients for common constituents of gasoline fuels (benzene, toluene, p-xylene and 1,2,4-trimethylbenzene) (BTX) in matrices with an extreme range of molar volumes and chemical affinities. Four synthetic mixtures were investigated comprising BTX with the bulk of the NAPL mixtures made up of either, ethylbenzene (an aromatic like BTX with similar molar volume); 1,3,5-trimethylbenzene (an aromatic with a greater molar volume); n-hexane (an aliphatic with a low molar volume); and n-decane (an aliphatic with a high molar volume). Equilibrium solubility values for the constituents were under-predicted by Raoult's law by up to 30% (higher experimental concentrations) for the mixture with n-hexane as a filler and over-predicted by up to 12% (lower experimental concentrations) for the aromatic mixtures with ethylbenzene and 1,3,5-trimethylbenzene as fillers. Application of PP-LFER (poly-parameter linear free energy relationship) model for non-ideal mixtures also resulted in poor correlation between experimentally measured and predicted concentrations, indicating that differences in chemical affinities can be the major cause of deviation from ideal behavior. Synthetic mixtures were compared with the dissolution behavior of fresh and naturally weathered unleaded gasoline. The presence of lighter aliphatic components in the gasoline had a profound effect on estimating effective solubility due to chemical affinity differences (estimated at 0.0055 per percentage increase in the molar proportion of aliphatic) as well as reduced molar volumes (estimated at -0.0091 in the activity coefficient per unit increase in molar volume, mL/mol). Previously measured changes in activity coefficients due to natural weathering of 0.25 compares well to 0.27 calculated here based on changes in the chemical affinity and molar volumes. The study suggests that the initial estimation of the composition of a fuel is crucial in evaluating dissolution processes due to ideal and non-ideal dissolution, and in predicting long term dissolution trends and the longevity of NAPL petroleum plume risks.

Scenario-based modelling of mass transfer mechanisms at a petroleum contaminated field site-numerical implications

Knowledge about distribution of dissolved plumes and their influencing factors is essential for risk assessment and remediation of light non-aqueous phase liquid contamination in groundwater. Present study deals with the applicability of numerical model for simulating various hydro-geological scenarios considering non-uniform source distribution at a petroleum contaminated site in Chennai, India. The complexity associated with the hydrogeology of the site has limited scope for on-site quantification of petroleum pipeline spillage. The change in fuel composition under mass-transfer limited conditions was predicted by simultaneously comparing deviations in aqueous concentrations and activity coefficients (between Raoult's law and analytical approaches). The effects of source migration and weathering on the dissolution of major soluble fractions of petroleum fuel were also studied in relation to the apparent change in their activity coefficients and molar fractions. The model results were compared with field observations and found that field conditions were favourable for biodegradation, especially for the aromatic fraction (benzene and toluene (nearly 95% removal), polycyclic aromatic hydrocarbons (up to 65% removal) and xylene (nearly 45% removal). The results help to differentiate the effect of compositional non-ideality from rate-limited dissolution towards tailing of less soluble compounds (alkanes and trimethylbenzene). Although the effect of non-ideality decreased with distance from the source, the assumption of spatially varying residual saturation could effectively illustrate post-spill scenario by estimating the consequent decrease in mass transfer rate.