Refinement of the gradient method for the estimation of natural source zone depletion at petroleum contaminated sites

Intercomparison of in-well, depth profile and online soil gas measurements for estimating petroleum natural source zone depletion rates: Will in-well gas samples suffice?

Natural Source Zone Deletion (NSZD) is a viable long-term management option for sites impacted by petroleum hydrocarbon fuels. NSZD rate estimation methods for petroleum mass losses often use soil gas gradients of oxygen, carbon dioxide, methane or vapour concentrations through the vadose zone. Seeking greater efficiencies, we investigated if existing short-screened wells are reliable for representative sampling of soil gases in a vadose zone undergoing NSZD. At gasoline, diesel, aviation-gasoline, crude oil and background locations, we compared soil gas compositions from three methods: (i) multilevel samplers (VZ) on three occasions over 18 months (three per location); (ii) oxygen sensors (GP) at multiple depths and tracked over 3 years (one per location); and (iii) short-screened monitoring wells (MW) across the water table (five per location) and sampled as per VZ. GP sensors showed seasonal and other variations, yielding variable NSZD estimates. Oxygen depth profiles from GP and VZ measurements showed good correspondence despite fluctuations. Oxygen, carbon dioxide and methane concentrations from the deepest sampling ports of VZ samplers correlated strongly with concentrations from adjacent MW installations, indicating in-well sampling in the field can provide reliable data. The advantages and limitations of using MW data in accelerating NSZD field assessments are outlined.

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.

Soil gas gradient method for estimating natural source zone depletion rates of LNAPL and specific chemicals of concern

This paper presents a simplified approach for the soil gas gradient method for estimating natural source zone depletion (NSZD) rates of specific contaminants of concern (COCs) at sites contaminated by light non-aqueous phase liquids (LNAPL). Traditional approaches to quantify COC-specific NSZD rates often rely on numerical or analytical reaction-transport models that require detailed site-specific data. In contrast, the proposed method employs simple analytical solutions, making it more accessible to practitioners. Specifically, it requires only the maximum soil gas concentration, the effective diffusion coefficient, and the diffusive reaction length calculated from vertical soil gas concentration profiles. The simplified approach was validated against a reactive transport numerical model reported in the literature, showing consistent results within the same order of magnitude for BTEX NSZD rates at a gasoline spill site in South Carolina. Further validation using a larger dataset involved comparing NSZD rate estimates for benzene and total petroleum hydrocarbons (TPH) against those obtained using BioVapor, utilizing empirical soil gas data from the USEPA Petroleum Vapor Intrusion Database. Results demonstrated a strong correlation between NSZD rates and maximum soil gas concentrations, allowing the development of a rapid screening approach based only on the measured soil gas concentrations and literature values for diffusion coefficients and diffusive reaction lengths. This approach aligned well with previous modeling studies and was consistent with literature values for TPH NSZD rates. Overall, both the simplified and screening approaches offer practical, easy-to-use tools for evaluating temporal variability in natural attenuation rates, supporting baseline assessments and ongoing performance evaluations of remediation at LNAPL sites.

Sorption and attenuation of petroleum VOCs in five unsaturated soils: Microcosms and column experiments

The fate of volatile organic compounds (VOC) vapors in the unsaturated zone is the basis for evaluating the natural attenuation potential and vapor intrusion risk. Microcosm and column experiments were conducted to study the effects chemical speciation and soil types/properties on the fate of petroleum VOCs in unsaturated zone. The biodegradation and total attenuation rates of the seven VOCs obtained by microcosm experiments in black soil and yellow earth were also generally higher than those in floodplain soil, lateritic red earth, and quartz sand. The VOC vapors in floodplain soil, lateritic red earth, and quartz sand showed slow total attenuation rates (<0.3 d-1). N-pentane, methylcyclopentane, and methylcyclohexane showed lower biodegradation rates than octane and three monoaromatic hydrocarbons. Volatilization into the atmosphere and biodegradation are two important natural attenuation paths for VOCs in unsaturated soil columns. The volatilization loss fractions of different volatile hydrocarbons in all five unsaturated soils were generally in the order: n-pentane (93.5%-97.8%) > methylcyclopentane (77.2%-85.5%) > methylcyclohexane (53.5%-69.2%) > benzene (17.1%-73.3%) > toluene (0-45.7%) > octane (1.9%-34.2%) > m-xylene (0-5.7%). The fractions by volatilization into the atmosphere of all seven hydrocarbons in quartz sand, lateritic red earth, and floodplain soil were close and higher compared to the yellow earth and black soil. Overall, this study illustrated the important roles chemical speciation and soil properties in determining the vapor-phase transport and natural attenuation of VOCs in the unsaturated zone.

Impacts of LNAPL types on mechanisms and rate of natural source zone depletion

Understanding the mechanisms of natural source zone depletion (NSZD) will support an improved understanding of the long-term sustainability of NSZD as a site remedy and how NSZD rates may change over time. This is the first study that has quantified and compared the rate of three NSZD mechanisms (methanogenesis, vaporization, and aqueous biodegradation) between two chemically distinct light non-aqueous phase liquid (LNAPL) source zones (aliphatic-rich naphtha for Zone #1 vs aromatic-rich pyrolysis gasoline for Zone #2) within the same geologic and climate conditions. The rates of NSZD attributable to vaporization (400 mg C/m2/d vs. 300 mg C/m2/d) and aqueous biodegradation (92 mg C/m2/d vs. 67 mg C/m2/d) were similar for Zone #1 and #2; however, the rate of methanogenesis NSZD was 6x higher in Zone #1 (1000 mg C/m2/d vs. 170 mg C/m2/d). These results suggest that the aliphatic hydrocarbons content in an LNAPL source may be a factor in the rate of methanogenesis NSZD. For both Zone #1 and #2, total NSZD rate determined using this "three mechanism" measurement method was in reasonable agreement with two other methods used to measure total NSZD rates (CO2 Gradient Method and Dynamic Closed Chamber Method), validating the "three mechanism" method as a tool to measure the total NSZD rate at a site and to provide an improved understanding of the predominant NSZD mechanism. Overall, this study highlights the importance of LNAPL type and chemical characteristics in determining source zone natural attenuation mechanism and its total rates.

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.

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.

Incorporating the Soil Gas Gradient Method and Functional Genes to Assess the Natural Source Zone Depletion at a Petroleum-Hydrocarbon-Contaminated Site of a Purification Plant in Northwest China

An increasing number of studies have demonstrated that natural source zone depletion (NSZD) in the vadose zone accounts for the majority (90%~99%) of the natural attenuation of light non-aqueous phase liquid (LNAPL). Until now, 0.05 to 12 kg/a.m² NSZD rates at tens of petroleum LNAPL source zones have been determined in the middle or late evolution stage of LNAPL release, in which limited volatile organic compounds (VOCs) and methane (CH₄) were detected. NSZD rates are normally estimated by the gradient method, yet the associated functional microbial activity remains poorly investigated. Herein, the NSZD at an LNAPL-releasing site was studied using both soil gas gradient methods quantifying the O₂, CO₂, CH₄, and VOCs concentrations and molecular biology methods quantifying the abundance of the pmoA and mcrA genes. The results showed that the methanogenesis rates were around 4 to 40 kg/a.m². The values were greater than the rates calculated by the sum of CH₄ escaping (0.3~1.2 kg/a.m²) and O₂ consuming (3~13 kg/a.m²) or CO₂ generating rates (2~4 kg/a.m²), suggesting that the generated CH₄ was oxidized but not thoroughly to CO₂. The functional gene quantification also supported the indication of this process. Therefore, the NSZD rates at the site roughly equaled the methanogenesis rates (4~40 kg/a.m²), which were greater than most of the previously studied sites with a 90th percentile value of 4 kg/a.m². The study extended the current knowledge of the NSZD and has significant implications for LNAPL remediation management.

Natural source zone depletion (NSZD) insights from over 15 years of research and measurements: A multi-site study

Site-average Natural Source Zone Depletion (NSZD) rates measured from 40 petroleum light non-aqueous phase liquid (LNAPL) source zone sites were compiled from researchers, project reports, and scientific papers. At each site, the following data were compiled: i) general site location; ii) LNAPL fuel type; iii) measurement method, number of locations, and number of measurements per location; and iv) calculated site-average NSZD rate in liters per hectare per year (L/ha/yr) per site and the associated measurement method (i.e., Gradient Method, Carbon Traps, Dynamic Closed Chamber (DCC), or Thermal Monitoring). The resulting dataset showed site-average NSZD rates that ranged from 650 to 152,000 L/ha/yr (70 to 16,250 gallons per acre per year (gal/acre/yr)), with a median value of 9,540 L/ha/yr (1,020 gal/acre/yr). The median site-average NSZD rate by type of fuel spill did not show a statistically significant difference between fuel types. When comparing the different NSZD measurement methods applied to the same sites, the site-average NSZD rates differed by up to 4.8 times (i.e., ratio of faster rate to slower rate), with a median difference of 2.1 times. No clear bias was observed between NSZD rate measurement methods. At four sites with calculations of NSZD rates by season, NSZD rates were typically higher during summer and fall compared to winter and spring. For these sites, Q₁₀ values (a measure of the increase in NSZD rate associated with a 10 °C increase in temperature) ranged from 0.8 to 15.1, with a median of 2.2. The implications of this study suggest that increasing mean annual soil temperature at a site using engineered methods could potentially increase the biodegradation rate (e.g., an increase of 10 °C could double the NSZD rate). Finally, for five sites with site-average NSZD rates for multiple years, average NSZD rates varied by 1.1 to 4.9 times across years. Overall, the evaluation of NSZD rates measured at 40 LNAPL sites suggests that measurable NSZD occurs across a broad range of LNAPL sites. Although NSZD rates vary across sites, fuel type is not the primary factor explaining observed differences in rates.

Tracking NSZD mass removal rates over decades: Site-wide and local scale assessment of mass removal at a legacy petroleum site

Long-term estimates of natural source zone depletion (NSZD) rates for petroleum LNAPL (light non-aqueous phase liquid) sites are not available. One-off measurements are often thought valid over the lifetime of LNAPL sites. In the context of site-wide LNAPL mass estimates, we report site-specific gasoline and diesel NSZD rates spanning 21-26 years. Using depth profiles of soil gases (oxygen, carbon dioxide, methane, volatiles) above LNAPL, NSZD rates were estimated in 1994, 2006 and 2020 for diesel and 1999, 2009 and 2020 for gasoline. Each date also had soil-core mass estimates, which together with NSZD rates allow estimation of the longevity for LNAPL presence. Site-wide coring (in 1992, 2002, 2007) estimated LNAPL mass reductions of 12,000 t. For diesel NSZD, the ratio of NSZD rates for 2006 (16,000-49,000 L/ha/y) to those in 2020 (2600-14,000 L/ha/y) was ~3-6. By 2020, the 1994 diesel NSZD rates would have predicted the entire removal of measured mass (16-42 kg/m²). For gasoline, NSZD rates in 1999 were extremely high (50,000-270,000 L/ha/y) but 9-27 times lower (5800-10,000 L/ha/y) a decade later. The gasoline NSZD rates in 1999 predicted near complete mass removal in 2-12 years, but 10-11 kg/m² was measured 10 and 21 years later which is 26% of the initial mass in 1999. The outcomes substantiate the need to understand NSZD rate changes over the lifetime of LNAPL-impacted sites.

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.