Modelling of iron cycling and its impact on the electron balance at a petroleum hydrocarbon contaminated site in Hnevice, Czech Republic

Enhancing in situ biodegradation in groundwater using pump and treat remediation: a proof of concept and modelling analysis of controlling variables

A remediation approach which uses pump and treatment (PAT) to enhance the biodegradation of organic contaminants by increasing dispersive mixing between plumes and groundwater was evaluated for a phenol-contaminated aquifer, using a reactive transport model which simulates kinetic reactions between an electron donor (ED) in the plume and electron acceptor (EA) in the groundwater. The influence of system design and operation was examined in six modelling scenarios. Injection or extraction of groundwater increases biodegradation above no action and the location, pumping rate, and distance between well(s) are important variables which influence biodegradation. An increase in pumping rate, distance of the wells from the plume centreline, and changing the flow direction increase dispersive mixing between the plume and groundwater. This increases plume spreading and the plume fringe interface, providing a greater flux of dissolved EAs for biodegradation. In general, injection of groundwater containing natural EAs enhances biodegradation more than extraction. The enhancement of biodegradation is sensitive to the relative fluxes of ED and EA, as controlled by the arrangement of the wells. In the best performing scenario, biodegradation was enhanced by 128%, compared with no action.

Effects of freeze-thaw cycles on methanogenic hydrocarbon degradation: Experiment and modeling

Cold regions are warming much faster than the global average, resulting in more frequent and intense freeze-thaw cycles (FTCs) in soils. In hydrocarbon-contaminated soils, FTCs modify the biogeochemical and physical processes controlling petroleum hydrocarbon (PHC) biodegradation and the associated generation of methane (CH₄) and carbon dioxide (CO₂). Thus, understanding the effects of FTCs on the biodegradation of PHCs is critical for environmental risk assessment and the design of remediation strategies for contaminated soils in cold regions. In this study, we developed a diffusion-reaction model that accounts for the effects of FTCs on toluene biodegradation, including methanogenic biodegradation. The model is verified against data generated in a 215 day-long batch experiment with soil collected from a PHC contaminated site in Ontario, Canada. The fully saturated soil incubations with six different treatments were exposed to successive 4-week FTCs, with temperatures oscillating between -10 °C and +15 °C, under anoxic conditions to stimulate methanogenic biodegradation. We measured the headspace concentrations and ¹³C isotope compositions of CH₄ and CO₂ and analyzed the porewater for pH, acetate, dissolved organic and inorganic carbon, and toluene. The numerical model represents solute diffusion, volatilization, sorption, as well as a reaction network of 13 biogeochemical processes. The model successfully simulates the soil porewater and headspace concentration time series data by representing the temperature dependencies of microbial reaction and gas diffusion rates during FTCs. According to the model results, the observed increases in the headspace concentrations of CH₄ and CO₂ by 87% and 136%, respectively, following toluene addition are explained by toluene fermentation and subsequent methanogenesis reactions. The experiment and the numerical simulation show that methanogenic degradation is the primary toluene attenuation mechanism under the electron acceptor-limited conditions experienced by the soil samples, representing 74% of the attenuation, with sorption contributing to 11%, and evaporation contributing to 15%. Also, the model-predicted contribution of acetate-based methanogenesis to total produced CH₄ agrees with that derived from the ¹³C isotope data. The freezing-induced soil matrix organic carbon release is considered as an important process causing DOC increase following each freezing period according to the calculations of carbon balance and SUVA index. The simulation results of a no FTC scenario indicate that, in the absence of FTCs, CO₂ and CH₄ generation would decrease by 29% and 26%, respectively, and that toluene would be biodegraded 23% faster than in the FTC scenario. Because our modeling approach represents the dominant processes controlling PHC biodegradation and the associated CH₄ and CO₂ fluxes, it can be used to analyze the sensitivity of these processes to FTC frequency and duration driven by temperature fluctuations.

Fermentation, methanotrophy and methanogenesis influence sedimentary Fe and As dynamics in As-affected aquifers in Vietnam

High arsenic (As) concentrations in groundwater are a worldwide problem threatening the health of millions of people. Microbial processes are central in the (trans)formation of the As-bearing ferric and ferrous minerals, and thus regulate dissolved As levels in many aquifers. Mineralogy, microbiology and dissolved As levels can vary sharply within aquifers, making high-resolution measurements particularly valuable in understanding the linkages between them. We conducted a high spatial resolution geomicrobiological study in combination with analysis of sediment chemistry and mineralogy in an alluvial aquifer system affected by geogenic As in the Red River delta in Vietnam. Microbial community analysis revealed a dominance of fermenters, methanogens and methanotrophs whereas sediment mineralogy along a 46 m deep core showed a diversity of Fe minerals including poorly crystalline Fe (II/III) and Fe(III) (oxyhydr)oxides such as goethite, hematite, and magnetite, but also the presence of Fe(II)-bearing carbonates and sulfides which likely formed as a result of microbially driven organic carbon (OC) degradation. A potential important role of methane (CH₄) as electron donor for reductive Fe mineral (trans)formation was supported by the high abundance of Candidatus Methanoperedens, a known Fe(III)-reducing methanotroph. Overall, these results imply that OC turnover including fermentation, methanogenesis and CH₄ oxidation are important mechanisms leading to Fe mineral (trans)formation, dissolution and precipitation, and thus indirectly affecting As mobility by changing the Fe-mineral inventory.

Spatial and temporal evolution of groundwater arsenic contamination in the Red River delta, Vietnam: Interplay of mobilisation and retardation processes

Geogenic arsenic (As) contamination of groundwater poses a major threat to global health, particularly in Asia. To mitigate this exposure, groundwater is increasingly extracted from low-As Pleistocene aquifers. This, however, disturbs groundwater flow and potentially draws high-As groundwater into low-As aquifers. Here we report a detailed characterisation of the Van Phuc aquifer in the Red River Delta region, Vietnam, where high-As groundwater from a Holocene aquifer is being drawn into a low-As Pleistocene aquifer. This study includes data from eight years (2010-2017) of groundwater observations to develop an understanding of the spatial and temporal evolution of the redox status and groundwater hydrochemistry. Arsenic concentrations were highly variable (0.5-510 μg/L) over spatial scales of <200 m. Five hydro(geo)chemical zones (indicated as A to E) were identified in the aquifer, each associated with specific As mobilisation and retardation processes. At the riverbank (zone A), As is mobilised from freshly deposited sediments where Fe(III)-reducing conditions occur. Arsenic is then transported across the Holocene aquifer (zone B), where the vertical intrusion of evaporative water, likely enriched in dissolved organic matter, promotes methanogenic conditions and further release of As (zone C). In the redox transition zone at the boundary of the two aquifers (zone D), groundwater arsenic concentrations decrease by sorption and incorporations onto Fe(II) carbonates and Fe(II)/Fe(III) (oxyhydr)oxides under reducing conditions. The sorption/incorporation of As onto Fe(III) minerals at the redox transition and in the Mn(IV)-reducing Pleistocene aquifer (zone E) has consistently kept As concentrations below 10 μg/L for the studied period of 2010-2017, and the location of the redox transition zone does not appear to have propagated significantly. Yet, the largest temporal hydrochemical changes were found in the Pleistocene aquifer caused by groundwater advection from the Holocene aquifer. This is critical and calls for detailed investigations.

Redox buffering and de-coupling of arsenic and iron in reducing aquifers across the Red River Delta, Vietnam, and conceptual model of de-coupling processes

Analysis of over 500 groundwater samples from throughout the Red River Delta indicates de-coupling of dissolved arsenic (As) and dissolved iron (Fe). Sorting of all data along the redox potentials suggests re-adsorption of As released initially from Mn(IV)-oxyhydroxides and later from Fe(III)-oxyhydroxides on remaining ferric phases at moderate redox levels. A gradually decreasing specific surface area available for re-adsorption of As probably plays a role as a consequence of limited reactivity of more crystalline phases such as goethite and hematite. At low redox levels, concentrations of Fe and phosphate decrease, but As concentrations keep increasing and most As is present as As(III) with limited adsorption affinity. Based on the results of speciation modeling, the water is supersaturated with respect to siderite and vivianite. A general conceptual model of As and Fe behavior is presented, suggesting that coupled behavior is possible in two geochemical "windows", i.e., 1: between saturation of remaining adsorption sites and the onset of siderite and vivianite precipitation, and 2: after the beginning of secondary sulfide phases precipitation and during methanogenesis. The de-coupling of As from Fe is common and has been observed at many sites around the world where As is released as a consequence of redox processes, e.g., in Bangladesh, West Bengal and Assam in India, the Mekong Delta in Cambodia and Vietnam, and Taiwan. The presented general conceptual model of de-coupling processes can be applied to the interpretation of As and Fe data, and, thus, it can help in the preparation of a site conceptual model which is a necessary prerequisite for reactive transport modeling.

Anaerobic biodegradation of dissolved ethanol in a pilot-scale sand aquifer: Variability in plume (redox) biogeochemistry

The use of ethanol in alternative fuels has led to contamination of groundwater with high concentrations of this easily biodegradable organic compound. Previous laboratory and field studies have shown vigorous biodegradation of ethanol plumes, with prevalence of reducing conditions and methanogenesis. The objective of this study was to further our understanding of the dynamic biogeochemistry processes, especially dissolved gas production, that may occur in developing and aging plume cores at sites with ethanol or other organic contamination of groundwater. The experiment performed involved highly-detailed spatial and temporal monitoring of ethanol biodegradation in a 2-dimensional (175cm high×525cm long) sand aquifer tank for 330days, with a vertical shift in plume position and increased nutrient inputs occurring at ~Day 100. Rapid onset of fermentation, denitrification, sulphate-reduction and iron(III)-reduction occurred following dissolved ethanol addition, with the eventual widespread development of methanogenesis. The detailed observations also demonstrate a redox zonation that supports the plume fringe concept, secondary reactions resulting from a changing/moving plume, and time lags for the various biodegradation processes. Additional highlights include: i) the highest dissolved H₂ concentrations yet reported for groundwater, possibly linked to vigorous fermentation in the absence of common terminal electron-acceptors (i.e., dissolved oxygen, nitrate, and sulphate, and iron(III)-minerals) and methanogenesis; ii) evidence of phosphorus nutrient limitation, which stalled ethanol biodegradation and perhaps delayed the onset of methanogenesis; and iii) the occurrence of dissimilatory nitrate reduction to ammonium, which has not been reported for ethanol biodegradation to date.

Elucidating the fate of a mixed toluene, DHM, methanol, and i-propanol plume during in situ bioremediation

Organic pollutants such as solvents or petroleum products are widespread contaminants in soil and groundwater systems. In-situ bioremediation is a commonly used remediation technology to clean up the subsurface to eliminate the risks of toxic substances to reach potential receptors in surface waters or drinking water wells. This study discusses the development of a subsurface model to analyse the performance of an actively operating field-scale enhanced bioremediation scheme. The study site was affected by a mixed toluene, dihydromyrcenol (DHM), methanol, and i-propanol plume. A high-resolution, time-series of data was used to constrain the model development and calibration. The analysis shows that the observed failure of the treatment system is linked to an inefficient oxygen injection pattern. Moreover, the model simulations also suggest that additional contaminant spillages have occurred in 2012. Those additional spillages and their associated additional oxygen demand resulted in a significant increase in contaminant fluxes that remained untreated. The study emphasises the important role that reactive transport modelling can play in data analyses and for enhancing remediation efficiency.

Microbial iron-redox cycling in subsurface environments

In addition to its central role in mediating electron-transfer reactions within all living cells, iron undergoes extracellular redox transformations linked to microbial energy generation through utilization of Fe(II) as a source of chemical energy or Fe(III) as an electron acceptor for anaerobic respiration. These processes permit microbial populations and communities to engage in cyclic coupled iron oxidation and reduction within redox transition zones in subsurface environments. In the present paper, I review and synthesize a few case studies of iron-redox cycling in subsurface environments, highlighting key biochemical aspects of the extracellular iron-redox metabolisms involved. Of specific interest are the coupling of iron oxidation and reduction in field and experimental systems that model redox gradients and fluctuations in the subsurface, and novel pathways and organisms involved in the redox cycling of insoluble iron-bearing minerals. These findings set the stage for rapid expansion in our knowledge of the range of extracellular electron-transfer mechanisms utilized by subsurface micro-organisms. The observation that closely coupled oxidation and reduction of iron can take place under conditions common to the subsurface motivates this expansion in pursuit of molecular tools for studying iron-redox cycling communities in situ.

Biogeochemical interactions among the arsenic, iron, humic substances, and microbes in mud volcanoes in southern Taiwan

Fluid and mud samples collected from Hsiaokunshui (HKS), Wushanting (WST), Yenshuikeng (YSK), Kunshuiping (KSP), Liyushan (LYS), and Sinyangnyuhu (SYNH) mud volcanoes of southwestern Taiwan were characterized for major ions, humic substances (HS) and trace elements concentrations. The relationship between the release of arsenic (As) and activities of sulfate-reducing bacteria has been assessed to understand relevant geochemical processes in the mud volcanoes. Arsenic (0.02-0.06 mg/L) and humic substances (4.13 × 10(-4) to 1.64 × 10(-3) mM) in the fluids of mud volcanoes showed a positive correlation (r = 0.99, p < 0.05) except in Liyushan mud volcano. Arsenic and iron in mud sediments formed two separate groups i) high As, but low Fe in HKS, WST, and SYNH; and ii) low As, but high Fe in the YSK, KSP, and LYS mud volcanoes. The Eh(S.H.E.) values of the mud volcano liquids were characterized by mild to strongly reducing conditions. The HKS, SYNH, and WST mud volcanoes (near the Chishan Fault) belongs to strong reducing environment (-33 to -116 mV), whereas the LYS, YSK, and KSP mud volcanoes located near the coastal plain are under mild reducing environment (-11 to 172 mV). At low Eh values mud volcanoes, saturation index (SI) values of poorly crystalline phases such as amorphous ferric hydroxide indicate understaturation, whereas saturation is reached in relatively high Eh(S.H.E.) values mud volcanoes. Arsenic contents in sediments are low, presumably due to its release to fluids (As/Fe ratio in YSK, KSP, and LYS sediment: 4.86 × 10(-4)-6.20 × 10(-4)). At low Eh(S.H.E.) values (mild to strong reducing environment), arsenic may co-precipitate with sulfides as a consequence of sulfate reduction (As/Fe ratios in WST, HKS, and SYNH sediments: 0.42-0.69).

Groundwater chemistry and arsenic mobilization in the Holocene flood plains in south-central Bangladesh

A comparative study of arsenic enrichment in the Bengal Delta (BD) was carried out in three alluvial aquifers in south-central Bangladesh. Investigated sites included Sonargaon in Narayanganj, Chandina in Comilla and Sirajdikhan in Munshiganj districts. At all sites samples from different depths were collected, and water chemistry and redox status vs. depth trends were determined. The concentrations of DOC and HCO(3)(-) were highest at Sirajdikhan site, while at the Sonargaon and Chandina sites the concentrations were lower. On the contrary, the NH(4)(+) concentration was high at the Chandina site as compared to the other sites. There was a good match between dissolved As and Fe at the Sirajdikhan and Sonargaon sites, but not at the Chandina site. The dissolved aqueous concentration of Mn was low at the Chandina site, which suggested that the Mn(IV) redox buffering step was missing. Speciation modeling indicated a possibility of siderite precipitation at all sites, but precipitation of rhodochrosite only at the Sonargaon and Sirajdikhan sites. At the Sirajdikhan site, the log P(CO2) values were very high (-1.37), which revealed the production of CO(2) in redox processes. Principal component analysis (PCA) indicated an impact of sea water and redox status of different samples. These results suggest that the dissolved As is de-coupled from dissolved Mn because when released, As is re-adsorbed onto the Fe(III) minerals in solid phase, as well as from dissolved Fe when precipitation of Fe(II) minerals controls the aqueous concentrations of Fe. In addition, several other concurrent redox processes may exert kinetic constraints depending on refractory characteristics of Fe(III) minerals.

Multicomponent simulations of contrasting redox environments at an LNAPL site

A two-dimensional multicomponent reactive transport modeling approach was used to simulate contaminant transport and the evolution of redox processes at a large-scale kerosene-contaminated site near Berlin, Germany. In contrast to previous site-scale modeling studies that focused either on one or two contaminants or on steady-state redox conditions, multiple contaminants and electron acceptors, including mineral phase Iron (III), were considered with an evolving redox zonation. Inhibition terms were used to switch between the different electron acceptor processes in the reaction scheme. The transient evolution of redox zones and contaminant plumes was simulated for two separate transects of the site, which have different geology and ground water recharge distributions and where quite different downstream contaminant and terminal electron-accepting process (TEAP) distributions are observed. The same reaction system, calibrated to measured concentrations along one of the transects, was used in both cases, achieving a reasonable match with observed concentrations. The differences between the two transects could thus to some extent be attributed to the different hydrological and hydrogeological conditions, in particular ground water recharge distributions. Long-term simulations showed that the distribution of TEAPs evolves as Fe(III) becomes depleted, with conditions becoming increasingly methanogenic, leading to changes in contaminant plume lengths. The models were applied to assess the potential effects of planned changes in land use at the site that may affect the ground water recharge distribution. The simulated redox zonation responded strongly to changes in recharge, which in turn led to changes in the contaminant plume lengths.