Relating subsurface temperature changes to microbial activity at a crude oil-contaminated site

Experimental research on the transport-transformation of organic contaminants under the influence of multi-field coupling at a site scale

Organic-contaminated shallow aquifers have become a global concern of groundwater contamination, yet little is known about the coupled effects of hydrodynamic-thermal-chemical-microbial (HTCM) multi-field on organic contaminant transport and transformation over a short time in aquifers. Therefore, this study proposed a quick and efficient field experimental method for the transport-transformation of contaminants under multi-field coupling to explore the relationship between organic contaminants (total petroleum hydrocarbon (TPH), polycyclic aromatic hydrocarbons (PAHs), benzene-toluene-ethylbenzene-xylene (BTEX) and phthalates acid esters (PAEs)) and multi-field factors. The results showed that hydrodynamics (affecting pH, p < 0.001) and temperature (affecting dissolved oxygen, pH and HCO3-, p < 0.05) mainly affected the organic contaminants indirectly by influencing the hydrochemistry to regulate redox conditions in the aquifer. The main degradation reactions of the petroleum hydrocarbons (TPH, PAHs and BTEX) and PAEs in the aquifer were sulfate reduction and nitrate reduction, respectively. Furthermore, the organic contamination was directly influenced by microbial communities, whose spatial patterns were shaped by the combined effects of the spatial pattern of hydrochemistry (induced by the organic contamination pressure) and other multi-field factors. Overall, our findings imply that the spatiotemporal patterns of organic contaminants are synergistically regulated by HTCM, with distinct mechanisms for petroleum hydrocarbons and PAEs.

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.

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.

Author response to technical commentary: Direct aerobic NSZD of a basalt vadose zone LNAPL source in Hawaii, McHugh et al., JCH #235 (2020)

Our Paper uses two independent methods (carbon dioxide flux and heat flux) to measure the rates of natural source zone depletion (NSZD) at a petroleum release site in Hawaii. The two methods yielded estimates of the NSZD rate that agreed within a factor of 2. We disagree with the technical commentors (Beckett et al., 2022). Specifically, the available data indicate that the observed NSZD is not occurring through a two-stage process of methane generation at the water table followed by methane oxidation in the vadose zone; rather direct aerobic oxidation of LNAPL in the vadose zone is the simplest and most likely explanation for the observed heat generation. In addition, the agreement between the two independent NSZD rate measurement methods and the temporal consistency in the measured heat flux across two field events provides confidence that the NSZD rates presented in the paper are not greatly over or under-estimated.

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.

Influence of Temperature on Denitrification and Microbial Community Structure and Diversity: A Laboratory Study on Nitrate Removal from Groundwater

Temperature is an extremely important environmental condition in the application of microbial denitrification for nitrate removal from groundwater. Understanding the nitrate removal efficiency of groundwater and the diversity, composition, and structure of microbial communities under different temperature conditions is of great significance for effective mitigation of groundwater nitrate pollution. This study investigated the effects of temperature on denitrification at 15 °C, 25 °C, 40 °C, and 45 °C. Moreover, the characteristics of microbial community structure and diversity were analyzed by combining high-throughput sequencing and polymerase chain reaction methods in order to fully clarify the denitrification efficiency under different temperature conditions. According to laboratory batch experiments and the findings of previous research, glucose was set as the carbon source and changes in “three nitrogen” indicators of the four temperature systems were mainly tested to clarify the effectiveness of nitrate removal. The maximum removal rates of nitrate nitrogen at 15 °C, 25 °C, 40 °C, and 45 °C were 44.05%, 87.03%, 99.26%, and 92.79%, respectively. Therefore, the most efficient nitrate removal can be achieved at 40℃. The Chao abundance indexes in the denitrification systems at 15 °C, 25 °C, 40 °C, and 45 °C were 1873, 352, 466, and 640, respectively. Therefore, the highest species richness was observed at 15 °C, but there were only a few dominant bacteria species. The composition of the bacterial community and the most dominant phylum varied at different temperatures. Among them, Gammaproteobacteria in Proteobacteria phylum plays an important role in the degradation of nitrate nitrogen. The relative abundance of Gammaproteobacteria at 15 °C, 25 °C, 40 °C, and 45 °C were 25.32%, 66.56%, 72.83%, and 3.47%. Tolumonas belongs to Gammaproteobacteria. The relative abundance of Tolumonas at 15 °C, 25 °C, 40 °C, and 45 °C were 9.41%, 65.47%, 62.49%, and 0.03%, respectively. The results of this study show that different temperature conditions affect the diversity, composition, and structure of the microbial community, thereby affecting the efficiency of denitrification for nitrate removal from groundwater.

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.

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.

Methane in groundwater from a leaking gas well, Piceance Basin, Colorado, USA

Site-specific and regional analysis of time-series hydrologic and geochemical data collected from 15 monitoring wells in the Piceance Basin indicated that a leaking gas well contaminated shallow groundwater with thermogenic methane. The gas well was drilled in 1956 and plugged and abandoned in 1990. Chemical and isotopic data showed the thermogenic methane was not from mixing of gas-rich formation water with shallow groundwater or natural migration of a free-gas phase. Water-level and methane-isotopic data, and video logs from a deep monitoring well, indicated that a shale confining layer ~125m below the zone of contamination was an effective barrier to upward migration of water and gas. The gas well, located 27m from the contaminated monitoring well, had ~1000m of uncemented annular space behind production casing that was the likely pathway through which deep gas migrated into the shallow aquifer. Measurements of soil gas near the gas well showed no evidence of methane emissions from the soil to the atmosphere even though methane concentrations in shallow groundwater (16 to 20mg/L) were above air-saturation levels. Methane degassing from the water table was likely oxidized in the relatively thick unsaturated zone (~18m), thus rendering the leak undetectable at land surface. Drilling and plugging records for oil and gas wells in Colorado and proxies for depth to groundwater indicated thousands of oil and gas wells were drilled and plugged in the same timeframe as the implicated gas well, and the majority of those wells were in areas with relatively large depths to groundwater. This study represents one of the few detailed subsurface investigations of methane leakage from a plugged and abandoned gas well. As such, it could provide a useful template for prioritizing and assessing potentially leaking wells, particularly in cases where the leakage does not manifest itself at land surface.

Relative contributions of microbial and infrastructure heat at a crude oil-contaminated site

Biodegradation of contaminants can increase the temperature in the subsurface due to heat generated from exothermic reactions, making temperature observations a potentially low-cost approach for determining microbial activity. For this technique to gain more widespread acceptance, it is necessary to better understand all the factors affecting the measured temperatures. Biodegradation has been occurring at a crude oil-contaminated site near Bemidji, Minnesota for 39 years, creating a quasi-steady-state plume of contaminants and degradation products. A model of subsurface heat generation and transport helps elucidate the contribution of microbial and infrastructure heating to observed temperature increases at this site. We created a steady-state, two-dimensional, heat transport model using previous-published parameter values for physical, chemical and biodegradation properties. Simulated temperature distributions closely match the observed average annual temperatures measured in the contaminated area at the site within less than 0.2 °C in the unsaturated zone and 0.4 °C in the saturated zone. The model results confirm that the observed subsurface heat from microbial activity is due primarily to methane oxidation in the unsaturated zone resulting in a 3.6 °C increase in average annual temperature. Another important source of subsurface heat is from the active, crude-oil pipelines crossing the site. The pipelines impact temperatures for a distance of 200 m and contribute half the heat. Model results show that not accounting for the heat from the pipelines leads to overestimating the degradation rates by a factor of 1.7, demonstrating the importance of identifying and quantifying all heat sources. The model results also highlighted a zone where previously unknown microbial activity is occurring at the site.

Long-term ERT monitoring of biogeochemical changes of an aged hydrocarbon contamination

Adequate management of contaminated sites requires information with improved spatio-temporal resolution, in particular to assess bio-geochemical processes, such as the transformation and degradation of contaminants, precipitation of minerals or changes in groundwater geochemistry occurring during and after remediation procedures. Electrical Resistivity Tomography (ERT), a geophysical method sensitive to pore-fluid and pore-geometry properties, permits to gain quasi-continuous information about subsurface properties in real-time and has been consequently widely used for the characterization of hydrocarbon-impacted sediments. However, its application for the long-term monitoring of processes accompanying natural or engineered bioremediation is still difficult due to the poor understanding of the role that biogeochemical processes play in the electrical signatures. For in-situ studies, the task is further complicated by the variable signal-to-noise ratio and the variations of environmental parameters leading to resolution changes in the electrical images. In this work, we present ERT imaging results for data collected over a period of two years on a site affected by a diesel fuel contamination and undergoing bioremediation. We report low electrical resistivity anomalies in areas associated to the highest contaminant concentrations likely due transformations of the contaminant due to microbial activity and accompanying release of metabolic products. We also report large seasonal variations of the bulk electrical resistivity in the contaminated areas in correlation with temperature and groundwater level fluctuations. However, the amplitude of bulk electrical resistivity variations largely exceeds the amplitude expected given existing petrophysical models. Our results suggest that the variations in electrical properties are mainly controlled by microbial activity which in turn depends on soil temperature and hydrogeological conditions. Therefore, ERT can be suggested as a promising tool to track microbial activity during bioremediation even though further research is still needed to completely understand the bio-geochemical processes involved and their impact on electrical signatures.