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

Over a period of several decades multiple leaks of large volumes from storage facilities located near Hnevice (Czech Republic) have caused the underlying Quaternary aquifer to be severely contaminated with nonaqueous phase liquid (NAPL) petroleum hydrocarbons. Beginning in the late 1980's the NAPL plume started to shrink as a consequence of NAPL dissolution exceeding replenishment and due to active remediation. The subsurface was classified geochemically into four different zones, (i) a contaminant-free zone never occupied by NAPL or dissolved contaminants, (ii) a re-oxidation zone formerly occupied by NAPL, (iii) a zone currently occupied by NAPL, and (iv) a lower fringe zone between the overlying NAPL and the deeper underlying contaminant-free zone. The study investigated the spatial and temporal variability of the redox zonation at the Hnevice site and quantified the influence of iron-cycling on the overall electron balance. As a first step inverse geochemical modelling was carried out to identify possible reaction models and mass transfer processes. In a subsequent step, two-dimensional (forward) multi-component reactive transport modelling was performed to evaluate and quantify the major processes that control the geochemical evolution at the site. The study explains the observed enrichment of the lower fringe zone with ferrihydrite as a result of the re-oxidation of ferrous iron. It suggests that once the NAPL zone started to shrink the dissolution of previously formed siderite and FeS by oxygen and nitrate consumed a significant part of the oxidation capacity for a considerable time period and therefore limited the penetration of electron acceptors into the NAPL contaminated zone.