Three-dimensional density-dependent flow and multicomponent reactive transport modeling of chlorinated solvent oxidation by potassium permanganate

A popular method for the treatment of aquifers contaminated with chlorinated solvents is chemical oxidation based on the injection of potassium permanganate (KMnO₄). Both the high density (1025 gL^− 1) and reactivity of the treatment solution influence the fate of permanganate (MnO₄) in the subsurface and affect the degree of contaminant treatment. The MIN3P multicomponent reactive transport code was enhanced to simulate permanganate-based remediation, to evaluate the pathways of MnO₄ utilization, and to assess the role of density contrasts for the delivery of the treatment solution. The modified code (MIN3P-D) provides a direct coupling between density-dependent fluid flow, solute transport, contaminant treatment, and geochemical reactions. The model is used to simulate a field trial of TCE oxidation in a sandy aquifer that is underlain by an aquitard. Three-dimensional simulations are conducted for a coupled reactive system comprised of ten aqueous components, two mineral phases, TCE (dissolved, adsorbed, and NAPL), reactive organic matter, and including ion exchange reactions. Model parameters are constrained by literature data and a detailed data set from the field site under investigation. The general spatial and transient evolution in observed concentrations of the oxidant, dissolved TCE, and reaction products are adequately reproduced by the simulations. The model elucidates the important role of density-induced flow and transport on the distribution of the treatment solution into NAPL containing regions located at the aquifer–aquitard interface. Model results further suggest that reactions that do not directly affect the stability of MnO₄ have a negligible effect on solution density and MnO₄ delivery.     

Numerical modeling of coupled variably saturated fluid flow and reactive transport with fast and slow chemical reactions

The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.     

Inverse modeling of multicomponent reactive transport through single and dual porosity media

Compacted bentonite is foreseen as buffer material for high-level radioactive waste in deep geological repositories because it provides hydraulic isolation, chemical stability, and radionuclide sorption. A wide range of laboratory tests were performed within the framework of FEBEX (Full-scale Engineered Barrier EXperiment) project to characterize buffer properties and develop numerical models for FEBEX bentonite. Here we present inverse single and dual-continuum multicomponent reactive transport models of a long-term permeation test performed on a 2.5 cm long sample of FEBEX bentonite. Initial saline bentonite porewater was flushed with 5.5 pore volumes of fresh granitic water. Water flux and chemical composition of effluent waters were monitored during almost 4 years. The model accounts for solute advection and diffusion and geochemical reactions such as aqueous complexation, acid-base, cation exchange, protonation/deprotonation by surface complexation and dissolution/precipitation of calcite, chalcedony and gypsum. All of these processes are assumed at local equilibrium. Similar to previous studies of bentonite porewater chemistry on batch systems which attest the relevance of protonation/deprotonation on buffering pH, our results confirm that protonation/deprotonation is a key process in maintaining a stable pH under dynamic transport conditions. Breakthrough curves of reactive species are more sensitive to initial porewater concentration than to effective diffusion coefficient. Optimum estimates of initial porewater chemistry of saturated compacted FEBEX bentonite are obtained by solving the inverse problem of multicomponent reactive transport. While the single-continuum model reproduces the trends of measured data for most chemical species, it fails to match properly the long tails of most breakthrough curves. Such limitation is overcome by resorting to a dual-continuum reactive transport model.     

A field experiment on cation exchange-affected multicomponent solute transport in a sandy aquifer

A field experiment was performed in an aquifer in order to study multicomponent cation-exchange processes under natural flow conditions. The aquifer is a glacial outwash plain with sandy aquifer material having a cation-exchange capacity (CEC) of 1.0 meg/100 g. A continuous injection of groundwater spiked with sodium and potassium as chlorides was accomplished over 37 days to resemble leachate contamination from landfills. The plume was monitored by sampling in a dense spatial network (length 100 m, width 20 m) over a period of 2.5 years in order to obtain breakthrough curves and spatial contour maps of the chemical compounds. Na and especially K showed a substantial retardation caused by cation-exchange processes despite the low CEC of the aquifer material. The average velocity of K⁺ was only 10% of the velocity of chloride (0.7 m day⁻¹). The relative migration velocity of Na⁺ was not a constant in the plume, but apparently influenced by dilution. Ca²⁺ and Mg²⁺ were expelled from the cation-exchange sites of the aquifer material and subsequently transported with the same velocity as chloride. The breakthrough curves of the various compounds showed multiple peaks and low concentration zones. It was concluded by calculations with PHREEQE that changes in calcite equilibrium may occur in the lower part of the aquifer, while complexation processes seem to be of no importance. Cation exchange is then the most important process in this field experiment, and further evaluation of the data by a geochemical transport model including cation exchange is recommended.     

Analysis of a deep well recharge experiment by calibrating a reactive transport model with field data

This paper describes the modeling of the hydrogeochemical effects of deep well recharge of oxic water into an anoxic pyrite-bearing aquifer. Kinetic expressions have been used for mineral dissolution-precipitation rates and organic matter oxidation. Hydrological and chemical parameters of the model were calibrated to field measurements. The results showed that oxidation of pyrite (FeS(2)) and, to a lesser extent, organic matter dominate the changes in quality of the recharged water during its passage through the aquifer. The recharge leads to the consumption of oxygen and nitrate and the formation of sulfate and ferrihydrite. Complexation reactions, cation exchange and precipitation and dissolution of calcite, siderite and rhodochrosite were also identified through the modeling. Despite problems of non-uniqueness of the calibrated parameters, the model was used successfully to depict the geochemical processes occurring in the aquifer. Non-uniqueness can be avoided by constraining the model as much as possible to measurements and/or data from literature, although they cannot be considered always as fixed values and should be considered as stochastic variables instead.     

Multi-component reactive transport modeling of natural attenuation of an acid groundwater plume at a uranium mill tailings site

Natural attenuation of an acidic plume in the aquifer underneath a uranium mill tailings pond in Wyoming, USA was simulated using the multi-component reactive transport code PHREEQC. A one-dimensional model was constructed for the site and the model included advective-dispersive transport, aqueous speciation of 11 components, and precipitation-dissolution of six minerals. Transport simulation was performed for a reclamation scenario in which the source of acidic seepage will be terminated after 5 years and the plume will then be flushed by uncontaminated upgradient groundwater. Simulations show that successive pH buffer reactions with calcite, Al(OH)3(a), and Fe(OH)3(a) create distinct geochemical zones and most reactions occur at the boundaries of geochemical zones. The complex interplay of physical transport processes and chemical reactions produce multiple concentration waves. For SO4(2-) transport, the concentration waves are related to advection-dispersion, and gypsum precipitation and dissolution. Wave speeds from numerical simulations compare well to an analytical solution for wave propagation.     

Operator-splitting errors in coupled reactive transport codes for transient variably saturated flow and contaminant transport in layered soil profiles

One possible way of integrating subsurface flow and transport processes with (bio)geochemical reactions is to couple by means of an operator-splitting approach two completely separate codes, one for variably-saturated flow and solute transport and one for equilibrium and kinetic biogeochemical reactions. This paper evaluates the accuracy of the operator-splitting approach for multicomponent systems for typical soil environmental problems involving transient atmospheric boundary conditions (precipitation, evapotranspiration) and layered soil profiles. The recently developed HP1 code was used to solve the coupled transport and chemical equations. For steady-state flow conditions, the accuracy was found to be mainly a function of the adopted spatial discretization and to a lesser extent of the temporal discretization. For transient flow situations, the accuracy depended in a complex manner on grid discretization, time stepping and the main flow conditions (infiltration versus evaporation). Whereas a finer grid size reduced the numerical errors during steady-state flow or the main infiltration periods, the errors sometimes slightly increased (generally less than 50%) when a finer grid size was used during periods with a high evapotranspiration demand (leading to high pressure head gradients near the soil surface). This indicates that operator-splitting errors are most significant during periods with high evaporative boundary conditions. The operator-splitting errors could be decreased by constraining the time step using the performance index (the product of the grid Peclet and Courant numbers) during infiltration, or the maximum time step during evapotranspiration. Several test problems were used to provide guidance for optimal spatial and temporal discretization.     

Evaluating the sensitivity of a subsurface multicomponent reactive transport model with respect to transport and reaction parameters.

The input variables for a numerical model of reactive solute transport in groundwater include both transport parameters, such as hydraulic conductivity and infiltration, and reaction parameters that describe the important chemical and biological processes in the system. These parameters are subject to uncertainty due to measurement error and due to the spatial variability of properties in the subsurface environment. This paper compares the relative effects of uncertainty in the transport and reaction parameters on the results of a solute transport model. This question is addressed by comparing the magnitudes of the local sensitivity coefficients for transport and reaction parameters. General sensitivity equations are presented for transport parameters, reaction parameters, and the initial (background) concentrations in the problem domain. Parameter sensitivity coefficients are then calculated for an example problem in which uranium(VI) hydrolysis species are transported through a two-dimensional domain with a spatially variable pattern of surface complexation sites. In this example, the reaction model includes equilibrium speciation reactions and mass transfer-limited non-electrostatic surface complexation reactions. The set of parameters to which the model is most sensitive includes the initial concentration of one of the surface sites, the formation constant (Kf) of one of the surface complexes and the hydraulic conductivity within the reactive zone. For this example problem, the sensitivity analysis demonstrates that transport and reaction parameters are equally important in terms of how their variability affects the model results.     

Illuminating reactive microbial transport in saturated porous media: demonstration of a visualization method and conceptual transport model

We demonstrate a method to study reactive microbial transport in saturated translucent porous media using the bacteria Pseudomonas fluorescens 5RL genetically engineered to carry a plasmid with bioluminescence genes inducible by salicylate. Induced bacteria were injected into a cryolite grain filled chamber saturated with a sterile non-growth-promoting (phosphorus limited) chemical mixture containing salicylate as an aromatic hydrocarbon analogue. The amount of light produced by the bacteria serves as an estimator of the relative efficiency of aerobic biodegradation since bioluminescence is dependent on both salicylate and oxygen but only consumes oxygen. Bioluminescence was captured with a digital camera and analyzed to study the evolving spatial pattern of the bulk oxygen consuming reactions. As fluid flow transported the bacteria through the chamber, bioluminescence was observed to initially increase until an oxygen depletion zone developed behind the advective front. Bacterial transport was modeled with the advection dispersion equation and oxygen concentration was modeled assuming bacterial consumption via Monod kinetics with consideration of additional effects of rate-limited mass transfer from residual gas bubbles. Consistent with previous measurements, bioluminescence was considered proportional to oxygen consumed. Using the observed bioluminescence, model parameters were fit that were consistent with literature values and produced results in good agreement with the experimental data. These findings demonstrate potential for using this method to investigate the complex spatial and temporal dynamics of reactive microbial transport in saturated porous media.     

Process-based reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage

Reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage was used to integrate a comprehensive data set including pore water chemistry and solid phase data from several sampling events over a >3-year time period. The simulations consider the reduction of sulfate by the organic carbon-based treatment material and the removal of sulfate and iron by precipitation of reduced mineral phases including iron monosulfides and siderite. Additional parameters constraining the model include dissolved H2S, alkalinity and pH, as well as a suite of solid phase S-fractions identified by extractions. Influences of spatial heterogeneity necessitated the use of a 2-dimensional modeling approach. Simulating observed seasonal fluctuations and long-term changes in barrier reactivity required the use of temperature dependent rate coefficients and a multimodal Monod-type rate expression accounting for the variable reactivity of different organic carbon fractions. Simulated dissolved concentrations of SO4, Fe, H2S, alkalinity and pH, as well as solid phase accumulations of reduced sulfur phases generally compare well to observed trends over 23 months. Spatial variations, seasonal fluctuations and the time-dependent decline in reactivity were also captured. The modeling results generally confirm, and further strengthen, the existing conceptual model for the site. Overall sulfate reduction and S-accumulation rates are constrained with confidence within a factor of 1.5.     

Numerical simulation of reactive processes in an experiment with partially saturated bentonite

Bentonites are preferred materials for use as engineered barriers for high-level nuclear waste repositories. Simulation of geochemical processes in bentonite is therefore important for long-term safety assessment of those repositories. In this work, the porewater chemistry of a bentonite sample subject to simultaneous heating and hydration, as studied by Cuevas et al. [Cuevas, J., Villar, M., Fernández, A., Gómez, P., Martín, P., 1997. Porewaters extracted from compacted bentonite subjected to simultaneous heating and hydration. Applied Geochemistry 12, 473-481.], was assessed with a non-isothermal reactive transport model by coupling the geochemical software PHREEQC2 with the object-oriented FEM simulator GeoSys/RockFlow. Reactive transport modelling includes heat transport, two-phase flow, multicomponent transport and geochemical reactions in the liquid phase, i.e. ion exchange, mineral dissolution/precipitation and equilibrium reactions. Simulations show that the easily soluble minerals in bentonite determine the porewater chemistry. Temperature affects both two-phase flow and geochemical reactions. Porosity change due to dissolution/precipitation is low during the experiment. However, changes of the effective porosity caused by bentonite swelling can be very large. The simulated results agree well with the experimental data.     

Assessment of intrinsic bioremediation of a coal-tar-affected aquifer using Two-dimensional reactive transport and Biogeochemical mass balance approaches.

Expedited site characterization and groundwater monitoring using direct-push technology and conventional monitoring wells were conducted at a former manufactured gas plant site. Biogeochemical data and heterotrophic plate counts support the presence of microbially mediated remediation. By superimposing solutions of a two-dimensional reactive transport analytical model, first-order degradation rate coefficients ((day-1) ) of various compounds for the dissolved-phase plume were estimated (i.e., benzene [0.0084], naphthalene [0.0058], and acenaphthene [0.0011]). The total mass transformed by aerobic respiration, nitrate reduction, and sulfate reduction around the free-phase coal-tar dense-nonaqueous-phase-liquid region and in the plume was estimated to be approximately 4.5 kg/y using a biogeochemical mass-balance approach. The total mass transformed using the degradation rate coefficients was estimated to be approximately 3.6 kg/y. Results showed that a simple two-dimensional analytical model and a biochemical mass balance with geochemical data from expedited site characterization can be useful for rapid estimation of mass-transformation rates.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: dissolved-plume retardation

A natural gradient emplaced-source (ES) controlled field experiment was conducted at the Borden aquifer research site, Ontario, to study the transport of dissolved plumes emanating from residual dense nonaqueous-phase liquid (DNAPL) source zones. The specific objective of the work presented here is to determine the effects of solute and co-solute concentrations on sorption and retardation of dissolved chlorinated solvent-contaminant plumes. The ES field experiment comprised a controlled emplacement of a residual multicomponent DNAPL below the groundwater table and intensive monitoring of dissolved-phase plumes of trichloromethane (TCM), trichloroethylene (TCE), and perchloroethylene (PCE) plumes continuously generated in the aquifer down gradient from gradual source dissolution. Estimates of plume retardation (and dispersion) were obtained from 3-D numerical simulations that incorporated transient source input and flow regimes monitored during the test. PCE, the most retarded solute, surprisingly exhibited a retardation factor approximately 3 times lower than observed in a previous Borden tracer test by Mackay et al. [Water Resour. Res. 22 (1986) 2017] conducted approximately 150 m away. Also, an absence of temporal trend in PCE retardation contrasted with the previous Borden test. Supporting laboratory studies on ES site core indicated that sorption was nonlinear and competitive, i.e. reduced sorption of PCE was observed in the presence of TCE. Consideration of the effects of relatively high co-solute (TCE) concentration (competitive sorption) in addition to PCE concentration effects (nonlinear sorption) was necessary to yield laboratory-based PCE retardation estimates consistent with the field plume values. Concentration- and co-solute-based sorption and retardation analysis was also applied to the previous low-concentration pulse injection test of Mackay et al. [Water Resour. Res. 22 (1986) 2017] and was able to successfully predict the temporal field retardation trends observed in that test. While it is acknowledged that other "nonideal transport" effects may contribute, our analysis predicts differences in the PCE retardation magnitude and trend between the two experiments that are consistent with field observations based on the marked solute concentration differences that resulted from contrasting source conditions. Solute and co-solute concentration effects have heretofore received little attention, but may have wide significance in aquifers contaminated by point-source pollutants because many plumes contain mixed solutes over wide concentration ranges in strata that are likely subject to nonlinear sorption.     

Groundwater flow, multicomponent transport and biogeochemistry: development and application of a coupled process model

A research tool for modeling the reactive flow and transport of groundwater contaminants in multiple dimensions is presented. Arbitrarily complex coupled kinetic–equilibrium heterogeneous reaction networks, automatic code generation, transfer-function based solutions, parameter estimation, high-resolution methods for advection, and robust solvers for the mixed kinetic–equilibrium chemistry are some of the features of reactive flow and transport (RAFT) that make it a versatile research tool in the modeling of a wide variety of laboratory and field experiments. The treatment of reactions is quite general so that RAFT can be used to model biological, adsorption/desorption, complexation, and mineral dissolution/precipitation reactions among others. The integrated framework involving automated code generation and parameter estimation allows for the development, characterization, and evaluation of mechanistic process models. The model is described and used to solve a problem in competitive adsorption that illustrates some of these features. The model is also used to study the development of an in situ Fe(II)-zone by encouraging the growth of an iron-reducing bacterium with lactate as the electron donor. Such redox barriers are effective in sequestering groundwater contaminants such as chromate and TCE.     

Two-site kinetic modeling of bacteriophages transport through columns of saturated dune sand

Breakthrough curves, on a semi-log scale, from tests in porous media with block-input of viruses, bacteria, protozoa and colloidal particles often exhibit a typical skewness: a rather slowly rising limb and a smooth transition of a declining limb to a very long tail. One-site kinetic models fail to fit the rising and declining limbs together with the tail satisfactorily. Inclusion of an equilibrium adsorption site does not seem to improve simulation results. This was encountered in the simulation of breakthrough curves from a recent field study on the removal of bacteriophages MS2 and PRD1 by passage through dune sand. In the present study, results of laboratory experiments for the study of this issue are presented. Breakthrough curves of salt and bacteriophages MS2, PRDI, and phiX174 in 1 D column experiments have been measured. One- and two-site kinetic models have been applied to fit and predict breakthrough curves from column experiments. The two-site model fitted all breakthrough curves very satisfactorily, accounting for the skewness of the rising limb as well as for the smooth transition of the declining limb to the tail of the breakthrough curve. The one-site model does not follow the curvature of the breakthrough tail, leading to an overestimation of the inactivation rate coefficient for attached viruses. Interaction with kinetic site 1 is characterized by relatively fast attachment and slow detachment, whereas attachment to and detachment from kinetic site 2 is fast. Inactivation of viruses and interaction with kinetic site 2 provide only a minor contribution to removal. Virus removal is mainly determined by the attachment to site 1. Bacteriophage phiX174 attached more than MS2 and PRD1, which can be explained by the greater electrostatic repulsion that MS2 and PRD1 experience compared to the less negatively charged phiX174.     

Intercomparison of reactive transport models applied to UO2 oxidative dissolution and uranium migration

Oxidative dissolution of uranium dioxide (UO(2)) and the subsequent migration of uranium in a subsurface environment and an underground waste disposal have been simulated with reactive transport models. In these systems, hydrogeological and chemical processes are closely entangled and their interdependency has been analyzed in detail, notably with respect to redox reactions, kinetics of mineralogical evolution and hydrodynamic migration of species of interest. Different codes, where among CASTEM, CHEMTRAP and HYTEC, have been used as an intercomparison and verification exercise. Although the agreement between codes is satisfactory, it is shown that the discretization method of the transport equation (i.e. finite elements (FE) versus mixed-hybrid FE and finite differences) and the sequential coupling scheme may lead to systematic discrepancies.     

Colloid transport in a geochemically heterogeneous porous medium: aquifer tank experiment and modeling

To examine colloid transport in geochemically heterogeneous porous media at a scale comparable to field experiments, we monitored the migration of silica-coated zirconia colloids in a two-dimensional layered porous media containing sand coated to three different extents by ferric oxyhydroxides. Transport of the colloids was measured over 1.65 m and 95 days. Colloid transport was modeled by an advection-dispersion-deposition equation incorporating geochemical heterogeneity and colloid deposition dynamics (blocking). Geochemical heterogeneity was represented as favorable (ferric oxyhydroxide-coated) and unfavorable (uncoated sand) deposition surface areas. Blocking was modeled as random sequential adsorption (RSA). Release of deposited colloids was negligible. The time to colloid breakthrough after the onset of blocking increased with increasing ferric oxyhydroxide-coated surface area. As the ferric oxyhydroxide surface area increased, the concentration of colloids in the breakthrough decreased. Model-fits to the experimental data were made by inverse solutions to determine the fraction of surface area favorable for deposition and the deposition rate coefficients for the favorable (ferric oxyhydroxide-coated) and unfavorable sites. The favorable deposition rate coefficient was also calculated by colloid filtration theory. The model described the time to colloid breakthrough and the blocking effect reasonably well and estimated the favorable surface area fraction very well for the two layers with more than 1% ferric oxyhydroxide coating. If mica edges in the uncoated sand were considered as favorable surface area in addition to the ferric oxyhydroxide coatings, the model predicted the favorable surface area fraction accurately for the layer with less than 1% ferric oxyhydroxide coating.     

Interaction between water flow and spatial distribution of microbial growth in a two-dimensional flow field in saturated porous media

Bacterial growth and its interaction with water flow was investigated in a two-dimensional flow field in a saturated porous medium. A flow cell (56 x 44 x 1 cm) was filled with glass beads and operated under a continuous flow of a mineral medium containing nitrate as electron acceptor. A glucose solution was injected through an injection port, simulating a point source contamination. Visible light transmission was used to observe the distribution of the growing biomass and water flow during the experiment. At the end of the experiment (on day 31), porous medium samples were destructively collected and analyzed for abundance of total and active bacterial cells, bacterial cell volume and concentration of polysaccharides and proteins. Microbial growth was observed in two stripes along the length of the flow cell, starting at the glucose injection port, where highest biomass concentrations were obtained. The spatial distribution of biomass indicated that microbial activity was limited by transverse mixing between glucose and nitrate media, as only in the mixing zone between the media high biological activities were achieved. The ability of the biomass to change the flow pattern in the flow cell was observed, indicating that the biomass was locally reducing the hydraulic conductivity of the porous medium. This bioclogging effect became evident when the injection of the glucose solution was turned off and water flow still bypassed the area around the glucose injection port, preserving the flow pattern as it was during the injection of the glucose solution. As flow bypass was possible in this system, the average hydraulic properties of the flow cell were not affected by the produced biomass. Even in the vicinity of the injection port, the total volume of the bacterial cells remained below 0.01% of the pore space and was unlikely to be responsible for the bioclogging. However, the bacteria produced large amounts of extracellular polymeric substances (EPS), which likely caused the observed bioclogging effects.     

Ex situ bioremediation of a soil contaminated by mazut (heavy residual fuel oil)--a field experiment

Mazut (heavy residual fuel oil)-polluted soil was exposed to bioremediation in an ex situ field-scale (600 m(3)) study. Re-inoculation was performed periodically with biomasses of microbial consortia isolated from the mazut-contaminated soil. Biostimulation was conducted by adding nutritional elements (N, P and K). The biopile (depth 0.4m) was comprised of mechanically mixed polluted soil with softwood sawdust and crude river sand. Aeration was improved by systematic mixing. The biopile was protected from direct external influences by a polyethylene cover. Part (10 m(3)) of the material prepared for bioremediation was set aside uninoculated, and maintained as an untreated control pile (CP). Biostimulation and re-inoculation with zymogenous microorganisms increased the number of hydrocarbon degraders after 50 d by more than 20 times in the treated soil. During the 5 months, the total petroleum hydrocarbon (TPH) content of the contaminated soil was reduced to 6% of the initial value, from 5.2 to 0.3 g kg(-1) dry matter, while TPH reduced to only 90% of the initial value in the CP. After 150 d there were 96%, 97% and 83% reductions for the aliphatic, aromatic, and nitrogen-sulphur-oxygen and asphaltene fractions, respectively. The isoprenoids, pristane and phytane, were more than 55% biodegraded, which indicated that they are not suitable biomarkers for following bioremediation. According to the available data, this is the first field-scale study of the bioremediation of mazut and mazut sediment-polluted soil, and the efficiency achieved was far above that described in the literature to date for heavy fuel oil.