Measurement of Natural Losses of LNAPL Using CO₂ Traps

Evaluating Natural Source Zone Depletion and Enhanced Source Zone Depletion in laboratory columns via soil redox continuous sensing and microbiome characterization

To optimally employ Natural Source Zone Depletion (NSZD) and Enhanced Source Zone Depletion (ESZD) at sites impacted by light non-aqueous phase liquids (LNAPL), monitoring strategies are required. Emerging use of subsurface oxidation-reduction potential (ORP) sensors shows promise for tracking redox evolution, which reflects ongoing biogeochemical processes. However, further understanding of how soil redox dynamics relate to subsurface microbial activity and LNAPL degradation pathways is needed. In this work, soil ORP sensors and DNA and RNA sequencing-based microbiome analysis were combined to elucidate NSZD and ESZD (biostimulation via periodic sulfate addition and biosparging) processes in columns containing LNAPL-impacted soils from a former petroleum refinery. Results show expected relationships between continuous soil redox and active microbial communities. Continuous data revealed spatial and temporal detail that informed interpretation of the hydrocarbon biodegradation data. Redox increases were transient for sulfate addition, and sequencing revealed how hydrocarbon concentration and composition impacted microbiome structure and naphthalene degradation. Periodic biosparging did not result in fully aerobic conditions suggesting observed biodegradation improvements could be explained by alternative anaerobic metabolisms (e.g., iron reduction due to air oxidizing reduced iron). Collectively, data suggest combining continuous redox sensing with microbiome analysis provides insights beyond those possible with either monitoring tool alone.

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.

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.

Natural Source Zone Depletion (NSZD) Quantification Techniques: Innovations and Future Directions

Natural source zone depletion (NSZD) is an emerging technique for sustainable and cost-effective bioremediation of light non-aqueous phase liquid (LNAPL) in oil spill sites. Depending on regulatory objectives, NSZD has the potential to be used as either the primary or sole LNAPL management technique. To achieve this goal, NSZD rate (i.e., rate of bulk LNAPL mass depletion) should be quantified accurately and precisely. NSZD has certain characteristic features that have been used as surrogates to quantify the NSZD rates. This review highlights the most recent trends in technology development for NSZD data collection and rate estimation, with a focus on the operational and technical advantages and limitations of the associated techniques. So far, four principal techniques are developed, including concentration gradient (CG), dynamic closed chamber (DCC), CO2 trap and thermal monitoring. Discussions revolving around two techniques, “CO2 trap” and “thermal monitoring”, are expanded due to the particular attention to them in the current industry. The gaps of knowledge relevant to the NSZD monitoring techniques are identified and the issues which merit further research are outlined. It is hoped that this review can provide researchers and practitioners with sufficient information to opt the best practice for the research and application of NSZD for the management of LNAPL impacted sites.

Advanced methods for RNA recovery from petroleum impacted soils

Microbially-mediated hydrocarbon degradation is well documented. However, how these microbial processes occur in complex subsurface petroleum impacted systems remains unclear, and this knowledge is needed to guide technologies to enhance microbial degradation effectively. Analysis of RNA derived from soils impacted by petroleum liquids would allow for analysis of active microbial communities, and a deeper understanding of the dynamic biochemistry occurring during site remediation. However, RNA analysis in soils impacted with petroleum liquids is challenging due to: (A) RNA being inherently unstable, and (B) petroleum impacted soils containing problematic levels of polymerase chain reaction (PCR) inhibitors that must be removed to yield high-purity RNA for downstream analysis. A previously published soil wash pretreatment step and a commercially available DNA extraction kit protocol were combined and modified to be able to purify RNA from soils containing petroleum liquids.•

A key modification involved reformulation of the pretreatment solution via replacing water as the diluent with a commercially-available RNA preservation solution.

•

Methods were developed and demonstrated using cryogenically preserved soils from three former petroleum refineries. Results showed the new soil washing approach had no adverse effects on RNA recovery but did improve RNA quality, by PCR inhibitor removal, which in turn allows for characterization of active microbial communities present in petroleum impacted soils.

•

In summary, our method for extracting RNA from petroleum-impacted soils provides a promising new tool for resolving metabolic processes at sites as they progress toward restoration via natural and/or engineered remediation.

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

Rates of natural source zone depletion (NSZD) are increasingly being used to aid remedial decision making and light non-aqueous phase liquid (LNAPL) longevity estimates at petroleum release sites. Current NSZD estimate methods, based on analyses of carbon dioxide (CO₂) and oxygen (O₂) soil-gas concentration gradients ("gradient method") assume linear concentration profiles with depth. This assumption can underestimate the concentration gradients especially above LNAPL sources that are typically characterized by curvilinear or semi-curvilinear O₂ and CO₂ concentration profiles. In this work, we proposed a new method that relies on calculating the O₂ and CO₂ concentration gradient using a first-order reaction model. The method requires an estimate of the diffusive reaction length that can be easily derived from soil-gas concentration data. A simple step-by-step guide for applying the new method is provided. Nomographs were also developed to facilitate method application. Application of the nomographs using field data from published literature showed that NSZD rates could be underestimated by nearly an order of magnitude if reactivity in the vadose zone is not accounted for. The new method helps refine NSZD rates estimation and improve risk-based decision making at certain petroleum contaminated sites.

Quantifying natural source zone depletion at petroleum hydrocarbon contaminated sites: A comparison of 14C methods

Surficial CO₂ efflux surveys have been used to delineate hydrocarbon source zones in contaminated aquifers and provide estimates of hydrocarbon biodegradation rates. This approach requires distinguishing between CO₂ derived from petroleum degradation and CO₂ produced from natural soil respiration. To this end, radiocarbon has been used to differentiate between ¹⁴C-depleted CO₂ from hydrocarbon degradation and ¹⁴C-enriched CO₂ from natural soil respiration to effectively quantify the contribution of each source to total CO₂ efflux, and by deduction natural source zone depletion (NSZD) rates. In this study, a systematic method comparison has been conducted to evaluate available approaches for collecting CO₂ gas samples for radiocarbon analysis used to correct total CO₂ efflux measurements for quantifying natural source zone depletion rates. Gas samples for radiocarbon analysis were sampled from (i) dynamic closed chambers (located at ground surface), (ii) static chambers (also at ground surface), (iii) shallow soil gas probes (0.3 m bgs), and (iv) soil gas monitoring wells (~0.6 m below ground surface) during a CO₂ efflux survey conducted at the site of a historical pipeline rupture near Bemidji, MN. The mean fraction of radiocarbon (F¹⁴C) obtained from samples overlying the source zone were (i) 0.93 ± 0.01, (ii) 0.73 ± 0.03, (iii) 0.71 ± 0.04, and (iv) 0.41 ± 0.06, for the four methods respectively. These F¹⁴C values were used to apportion total CO₂ efflux measurements into contributions of contaminant-derived CO₂ efflux and natural soil respiration to evaluate natural source zone depletion processes. Results suggest that the method of radiocarbon sampling has a significant effect on the calculated fraction of the CO₂ efflux originating from contaminant-related soil respiration, with contributions ranging between 27% and 59% of total soil respiration. Results indicate that radiocarbon sampled from static chambers and shallow soil gas probes methods offer the best compromise between CO₂ sample yield and sample representativeness, providing the most reliable estimates of CO₂ effluxes originating from contaminant degradation. However, the results also show that for this study, all methods agree within a factor of <2.3 regarding the inferred NSZD rates.

Real-time soil and groundwater monitoring via spatial and temporal resolution of biogeochemical potentials

Oxidation-reduction potentials (ORP) govern the transformation of organic compounds in water and soils. Standard methods for measurements of ORPs in subsurface setting are deeply flawed due to heterogeneous samples from wells, failure to capture weakly poised redox couples, and biases with ex-situ measurements. In this study, we developed a real-time in-situ ORP sensor system that continuously measures biogeochemical electrical potentials using vertically distributed point sensing electrodes in direct contact with the soil. Three hundred thousand data points, providing a full range of aqueous ORP values (+ 600 to - 600 mV vs. Ag/AgCl) were collected over 513 days to spatially and temporally resolve subsurface biogeochemical processes at a former petroleum refinery. Water quality and microbial community data support the validity of the ORP data. In locations impacted by petroleum light non-aqueous phase liquids (LNAPLs), barometric pumping and ebullition events drive near-daily cycles of ORP changes in the vadose zone of 400 mV. When only dissolved phase hydrocarbons are present, near-daily redox cycles are absent and values for ORP indicate methanogenic conditions immediately about the water table. When hydrocarbons are not present, redox conditions are more oxidizing by + 400 to + 700 mV. The embedded electrodes revealed variations in hydrocarbon biodegradation in time and space that cannot be resolved by collection and analysis of conventional samples of groundwater and soil gas.

Direct aerobic NSZD of a basalt vadose zone LNAPL source in Hawaii

In recent years, a number of methods have been used to measure the biodegradation of petroleum light non-aqueous phase liquids (LNAPL) at petroleum release sites, a process known as natural source zone depletion (NSZD). Most commonly, NSZD rates have been measured at sites with unconsolidated geology and relatively shallow groundwater (<50 ft. bgs, <15 m bgs). For this study, we have used two methods (1. carbon dioxide flux measured using carbon traps and 2. heat flux based on subsurface temperature gradients) to measure NSZD rates at a petroleum release site in Hawaii with basalt geology and deep groundwater (>300 ft. bgs, >100 m bgs). Both methods documented the occurrence of NSZD at the facility and the two methods yield estimates of the NSZD rate that agreed within a factor of 2 (4600 to 7400 gal/yr; 17,000 to 28,000 L/yr for the flux method and 8600 to 13,000 gal/yr; 33,000 to 49,000 L/yr for the temperature method). Soil gas samples collected directly above the water table and at shallower depths within the vadose zone indicated aerobic conditions throughout the vadose zone (oxygen >13%) and no detectable methane. These results indicate that NSZD occurs at this site through the direct aerobic biodegradation of LNAPL rather than the two-step process of anaerobic methanogenesis followed by methane oxidation at a shallow depth interval documented at other sites.

Thermal estimation of natural source zone depletion rates without background correction

Real-time monitoring of subsurface temperature profiles is a promising approach to resolving natural source zone depletion (NSZD) rates for shallow petroleum liquids. Herein, a new "single stick" computational method for transforming temperature data into NSZD rates is advanced. The method is predicated on subsurface temperatures being a function of surface heating and cooling, and the heat associated with NSZD. Given subsurface temperature at two points, a system of two-equation two-unknown is used to resolve NSZD rates. Mathematical formulations and computational algorithms are validated through computational tests showing near perfect agreement between prescribed and predicted NSZD heating, and observed and predicted subsurface temperatures. The method is applied to temperature data from five field sites. Results include lower NSZD rates in areas where petroleum liquids are absent by a factor of 0.5-7.5 as compared to background correction methods. While the single stick method provides the lowest rate in unimpacted areas, it also provides reasonable rates as compared to two of three background correction methods. In addition, the single stick method yields a coefficient of variation equal to 6% in a triplicate analysis and reasonable estimates of NSZD in LNAPL-impacted areas without background corrected data.

A Screening Model to Predict Entrapped LNAPL Depletion

Accidental leakage of hydrocarbons is a common subsurface contamination scenario. Once released, the hydrocarbons migrate until they reach the vicinity of the uppermost portion of the saturated zone, where it accumulates. Whenever the amplitude of the water table fluctuation is high, the light non-aqueous phase liquid (LNAPL) may be completely entrapped in the saturated zone. The entrapped LNAPL, comprised of multicomponent products (e.g., gasoline, jet fuel, diesel), is responsible for the release of benzene, toluene, ethylbenzene, and xylenes (BTEX) into the water, thus generating the dissolved phase plumes of these compounds. In order to estimate the time required for source-zone depletion, we developed an algorithm that calculates the mass loss of BTEX compounds in LNAPL over time. The simulations performed with our algorithm provided results akin to those observed in the field and demonstrated that the depletion rate will be more pronounced in regions with high LNAPL saturation. Further, the LNAPL depletion rate is mostly controlled by flow rate and is less sensible to the biodegradation rate in the aqueous phase.

Using dynamic flux chambers to estimate the natural attenuation rates in the subsurface at petroleum contaminated sites

In this work, we introduce a screening method for the evaluation of the natural attenuation rates in the subsurface at sites contaminated by petroleum hydrocarbons. The method is based on the combination of the data obtained from standard source characterization with dynamic flux chambers measurements. The natural attenuation rates are calculated as difference between the flux of contaminants estimated with a non-reactive diffusive model starting from the concentrations of the contaminants detected in the source (soil and/or groundwater) and the effective emission rate of the contaminants measured using dynamic flux chambers installed at ground level. The reliability of this approach was tested in a contaminated site characterized by the presence of BTEX in soil and groundwater. Namely, the BTEX emission rates from the subsurface were measured in 4 seasonal campaigns using dynamic flux chambers installed in 14 sampling points. The comparison of measured fluxes with those predicted using a non-reactive diffusive model, starting from the source concentrations, showed that, in line with other recent studies, the modelling approach can overestimate the expected outdoor concentration of petroleum hydrocarbons even up to 4 orders of magnitude. On the other hand, by coupling the measured data with the fluxes estimated with the diffusive non-reactive model, it was possible to perform a mass balance to evaluate the natural attenuation loss rates of petroleum hydrocarbons during the migration from the source to ground level. Based on this comparison, the estimated BTEX loss rates in the test site were up to almost 0.5kg/year/m². These rates are in line with the values reported in the recent literature for natural source zone depletion. In short, the method presented in this work can represent an easy-to-use and cost-effective option that can provide a further line of evidence of natural attenuation rates expected at contaminated sites.

Understanding complex LNAPL sites: Illustrated handbook of LNAPL transport and fate in the subsurface

The goal of the paper is to highlight the management of the complexities and risks for light non-aqueous phase liquid (LNAPL) sites, and how the Illustrated Handbook of LNAPL Transport and Fate in the Subsurface (CL:AIRE, London. ISBN 978-1-905046-24-9. http://www.claire.co.uk/LNAPL; "LNAPL illustrated handbook") is useful guidance and a tool for professionals to understand these complexities and risks. The LNAPL illustrated handbook provides a clear and concise best-practice guidance document, which is a valuable decision support tool for use in discussions and negotiations regarding LNAPL impacted sites with respect to the risks of LNAPL. The LNAPL illustrated handbook is a user-friendly overview of the nature of LNAPL contamination in various geological settings including unconsolidated, consolidated, and fractured rock environments to best understand its fate and behavior leading to the appropriate management and/or remedial approach of the two major risks associated with a LNAPL source. As a source term, LNAPL has chemicals that form dissolved- and vapor-phase plumes, which are referred to as composition-based risks; and being a liquid there is the risk that the source may expand impacting a greater volume of the aquifer, which are referred to as saturation-based risks. There have been significant developments in recent years on the understanding of the complex behavior of LNAPL and associated groundwater and vapor plumes; however, the state of practice has often lagged these improvements in knowledge. The LNAPL illustrated handbook aids the site investigator, site owners, and regulators to understand these risks, and understand how these risks behave through better conceptual understanding of LNAPL transport and fate in the subsurface.

Measuring Carbon-based Contaminant Mineralization Using Combined CO2 Flux and Radiocarbon Analyses

A method is described which uses the absence of radiocarbon in industrial chemicals and fuels made from petroleum feedstocks which frequently contaminate the environment. This radiocarbon signal - or rather the absence of signal - is evenly distributed throughout a contaminant source pool (unlike an added tracer) and is not impacted by biological, chemical or physical processes (e.g., the ¹⁴C radioactive decay rate is immutable). If the fossil-derived contaminant is fully degraded to CO2, a harmless end-product, that CO2 will contain no radiocarbon. CO2 derived from natural organic matter (NOM) degradation will reflect the NOM radiocarbon content (usually <30,000 years old). Given a known radiocarbon content for NOM (a site background), a two end-member mixing model can be used to determine the CO2 derived from a fossil source in a given soil gas or groundwater sample. Coupling the percent CO2 derived from the contaminant with the CO2 respiration rate provides an estimate for the total amount of contaminant degraded per unit time. Finally, determining a zone of influence (ZOI) representing the volume from which site CO2 is collected allows determining the contaminant degradation per unit time and volume. Along with estimates for total contaminant mass, this can ultimately be used to calculate time-to-remediate or otherwise used by site managers for decision-making.

Comparison of bacterial and archaeal communities in depth-resolved zones in an LNAPL body

Advances in our understanding of the microbial ecology at sites impacted by light non-aqueous phase liquids (LNAPLs) are needed to drive development of optimized bioremediation technologies, support longevity models, and develop culture-independent molecular tools. In this study, depth-resolved characterization of geochemical parameters and microbial communities was conducted for a shallow hydrocarbon-impacted aquifer. Four distinct zones were identified based on microbial community structure and geochemical data: (i) an aerobic, low-contaminant mass zone at the top of the vadose zone; (ii) a moderate to high-contaminant mass, low-oxygen to anaerobic transition zone in the middle of the vadose zone; (iii) an anaerobic, high-contaminant mass zone spanning the bottom of the vadose zone and saturated zone; and (iv) an anaerobic, low-contaminant mass zone below the LNAPL body. Evidence suggested that hydrocarbon degradation is mediated by syntrophic fermenters and methanogens in zone III. Upward flux of methane likely contributes to promoting anaerobic conditions in zone II by limiting downward flux of oxygen as methane and oxygen fronts converge at the top of this zone. Observed sulfate gradients and microbial communities suggested that sulfate reduction and methanogenesis both contribute to hydrocarbon degradation in zone IV. Pyrosequencing revealed that Syntrophus- and Methanosaeta-related species dominate bacterial and archaeal communities, respectively, in the LNAPL body below the water table. Observed phylotypes were linked with in situ anaerobic hydrocarbon degradation in LNAPL-impacted soils.