Biological remediation of explosives and related nitroaromatic compounds

Nitroaromatics form an important group of recalcitrant xenobiotics. Only few aromatic compounds, bearing one nitro group as a substituent of the aromatic ring, are produced as secondary metabolites by microorganisms. The majority of nitroaromatic compounds in the biosphere are industrial chemicals such as explosives, dyes, polyurethane foams, herbicides, insecticides and solvents. These compounds are generally recalcitrant to biological treatment and remain in the biosphere, where they constitute a source of pollution due to both toxic and mutagenic effects on humans, fish, algae and microorganisms. However, relatively few microorganisms have been described as being able to use nitroaromatic compounds as nitrogen and/or carbon and energy source. The best-known nitroaromatic compound is the explosive TNT (2,4,6-trinitrotoluene). This article reviews the bioremediation strategies for TNT-contaminated soil and water. It comes to the following conclusion: The optimal remediation strategy for nitroaromatic compounds depends on many site-specific factors. Composting and the use of reactor systems lend themselves to treating soils contaminated with high levels of explosives (e.g. at former ammunition production facilities, where areas with a high contamination level are common). Compared to composting systems, bioreactors have the major advantage of a short treatment time, but the disadvantage of being more labour intensive and more expensive. Studies indicate that biological treatment systems, which are based on the activity of the fungus Phanerochaete chrysosporium or on Pseudomonas sp. ST53, might be used as effective methods for the remediation of highly contaminated soil and water. Phytoremediation, although not widely used now, has the potential to become an important strategy for the remediation of soil and water contaminated with explosives. It is best suited where contaminant levels are low (e.g. at military sites where pollution is rather diffuse) and where larger contaminated surfaces or volumes have to be treated. In addition, phytoremediation can be used as a polishing method after other remediation treatments, such as composting or bioslurry, have taken place. This in-situ treatment method has the advantage of lower treatment costs, but has the disadvantage of a considerable longer treatment time. In order to improve the cost-efficiency, phytoremediation of nitroaromatics (and other organic xenobiotics) could be combined with bio-energy production. This requires, however, detailed knowledge on the fate of the contaminants in the plants as well as the development of efficient treatment methods for the contaminated biomass that minimise the spreading of the contaminants into the environment during post harvest treatment.     

The effect of particle size reduction by grinding on subsampling variance for explosives residues in soil

Efforts to characterize the surface soil contamination on military training ranges have been compromised by the inability to obtain representative subsamples of soils submitted to analytical laboratories for determination of explosives residues. Two factors affecting subsampling error for explosives residues were examined using soils collected from hand grenade and anti-tank ranges. These factors were increased subsample size and particle size reduction prior to subsampling of soils. Increasing the subsample size from 2 to 50 g did not reduce the soil subsampling error because of the extreme heterogeneous distribution of the solid contaminants. Alternatively, particle size reduction by machine grinding on a ring mill reduced subsampling error to less than 10% relative standard deviation for replicate analyses using 10-g subsamples.     

Comparison of biotic and abiotic treatment approaches for co-mingled perchlorate, nitrate, and nitramine explosives in groundwater

Biological and abiotic approaches for treating co-mingled perchlorate, nitrate, and nitramine explosives in groundwater were compared in microcosm and column studies. In microcosms, microscale zero-valent iron (mZVI), nanoscale zero-valent iron (nZVI), and nickel catalyzed the reduction of RDX and HMX from initial concentrations of 9 and 1 mg/L, respectively, to below detection (0.02 mg/L), within 2 h. The mZVI and nZVI also degraded nitrate (3 mg/L) to below 0.4 mg/L, but none of the metal catalysts were observed to appreciably reduce perchlorate ( approximately 5 mg/L) in microcosms. Perchlorate losses were observed after approximately 2 months in columns of aquifer solids treated with mZVI, but this decline appears to be the result of biodegradation rather than abiotic reduction. An emulsified vegetable oil substrate was observed to effectively promote the biological reduction of nitrate, RDX and perchlorate in microcosms, and all four target contaminants in the flow-through columns. Nitrate and perchlorate were biodegraded most rapidly, followed by RDX and then HMX, although the rates of biological reduction for the nitramine explosives were appreciably slower than observed for mZVI or nickel. A model was developed to compare contaminant degradation mechanisms and rates between the biotic and abiotic treatments.     

Comparing solid phase extraction and direct injection for the analysis of ultra-trace levels of relevant explosives in lake water and tributaries using liquid chromatography-electrospray tandem mass spectrometry

Off-line solid phase extraction and direct injection analysis were evaluated for the determination of traces of explosives such as TNT and its mono and diamino metabolites, HMX, RDX, nitroglycerin and PETN in lake water and tributaries applying liquid chromatography-electrospray tandem mass spectrometry. Improved chromatographic separation was achieved on a phenyl based stationary phase with baseline resolution of the mono- and diamino metabolites of TNT. Identification and quantification of the target compounds was performed by multiple reaction monitoring applying electrospray ionization in either the positive mode for the diaminometabolites of TNT or the negative mode for all other compounds. An extensive method validation was performed and limits of quantification were obtained for the explosives in preconcentrated lake water samples from 0.03 to 1 ng l(-1) and 0.1 to 5 ng l(-1) in river water. Direct injection analysis revealed comparable results to preconcentrated water samples for the most persistent explosives. Analysis of lake water samples collected at different depths showed the presence of HMX, RDX and PETN at concentrations from 0.1 to 0.4 ng l(-1). The analysis of main tributaries revealed concentrations from 0.1 to 0.9 ng l(-1) of the same compounds. They seem to be responsible for the contamination of the explosives in the lakes.     

Determination of mono-, di- and trinitronaphthalenes in soil samples contaminated by explosives

We investigated the extraction of eight mono-, di- and trinitronaphthalenes from soil samples by ultrasonic, respectively, soxhlet extraction. Sonication was the more convenient procedure with a usually higher repeatability than the soxhlet method. All extracts were analysed by high-performance liquid chromatography with diode array detection. These methods were used to determine the concentration of nitronaphthalenes in five soil samples polluted by explosives. The most abundant contaminants were 1,5- and 1,8-dinitronaphthalene as well as 1,3,5-, 1,3,8- and 1,4,5-trinitronaphthalene. In all samples the sum of trinitronaphthalenes exceeded the sum of mono- and dinitronaphthalenes. It is recommendable to include these eight compounds in the analytical exploration of sites contaminated by explosives.     

Anaerobic biotransformation of explosives in aquifer slurries amended with ethanol and propylene glycol

Hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX), and 2,4,6-trinitrotoluene (TNT) are explosives that are frequently found as environmental contaminants on military installations. Hydrogen has been shown to support the anaerobic transformation of these explosives. We investigated ethanol and propylene glycol as electron donors for providing syntrophically produced H2 for stimulating the anaerobic biodegradation of explosives in contaminated soil. The study was conducted using anoxic microcosms constructed with slurries of the contaminated soil and groundwater. The addition of 5mM ethanol and propylene glycol enhanced the biodegradation of RDX and HMX relative to the control bottles. Ethanol was depleted within about 20 days, resulting in the transient formation of hydrogen, acetate, and methane. The hydrogen headspace concentration increased from 8 ppm to 1838 ppm before decreasing to background concentrations. Propylene glycol was completely degraded after 15 days, forming hydrogen, propionate, and acetate as end-products. The hydrogen headspace concentrations increased from 56 ppm to 628 ppm before decreasing to background concentrations. No methane formation was observed during the incubation period of 48 days. Our findings indicate the addition of ethanol and propylene to the aquifer slurries increased the hydrogen concentrations and enhanced the biotransformation of RDX and HMX in the explosive-contaminated soil.     

A Probability Model for Evaluating Building Contamination from an Environmental Event

ABSTRACT

Asbestos dust and bioaerosol sampling data from suspected contaminated zones in buildings allowed development of an environmental data evaluation protocol based on the differences in frequency of detection of a target contaminant between zones of comparison. Under the assumption that the two test zones of comparison are similar, application of population proportion probability calculates the significance of observed differences in contaminant levels. This was used to determine whether levels of asbestos dust contamination detected after a fire were likely the result of smoke-borne contamination, or were caused by pre-existing/background conditions.

Bioaerosol sampling from several sites was also used to develop the population proportion probability protocol. In this case, significant differences in indoor air contamination relative to the ambient conditions were identified that were consistent with the visual observations of contamination. Implicit in this type of probability analysis is a definition of “contamination” based on significant differences in contaminant levels relative to a control zone. Detection of a suspect contaminant can be assessed as to possible sources(s) as well as the contribution made by pre-existing (i.e., background) conditions, provided the test and control zones are subjected to the same sampling and analytical methods.     

A probability model for evaluating building contamination from an environmental event

Asbestos dust and bioaerosol sampling data from suspected contaminated zones in buildings allowed development of an environmental data evaluation protocol based on the differences in frequency of detection of a target contaminant between zones of comparison. Under the assumption that the two test zones of comparison are similar, application of population proportion probability calculates the significance of observed differences in contaminant levels. This was used to determine whether levels of asbestos dust contamination detected after a fire were likely the result of smoke-borne contamination, or were caused by pre-existing/background conditions. Bioaerosol sampling from several sites was also used to develop the population proportion probability protocol. In this case, significant differences in indoor air contamination relative to the ambient conditions were identified that were consistent with the visual observations of contamination. Implicit in this type of probability analysis is a definition of "contamination" based on significant differences in contaminant levels relative to a control zone. Detection of a suspect contaminant can be assessed as to possible sources(s) as well as the contribution made by pre-existing (i.e., background) conditions, provided the test and control zones are subjected to the same sampling and analytical methods.     

Redistribution of contaminants by a fluctuating water table in a micro-porous, double-porosity aquifer: field observations and model simulations

Large seasonal fluctuations of the water table are characteristic of aquifers with a low specific yield, including those fractured, double-porosity aquifers that have significant matrix porosity containing virtually immobile porewater, such as the Chalk of northern Europe. Where these aquifers are contaminated, a strong relationship between water table elevation and contaminant concentration in groundwater is commonly observed, of significance to the assessment, monitoring, and remediation of contaminated groundwater. To examine the processes governing contaminant redistribution by a fluctuating water table within the 'seasonally unsaturated zone', or SUZ, profiles of porewater solute concentrations have been established at a contaminated site in southern England. These profiles document the contaminant distribution in porewater of the Chalk matrix over the SUZ at a greater level of detail than recorded previously. A novel double-porosity solute transport code has been developed to simulate the evolution of the SUZ matrix porewater contaminant profiles, given a fluctuating water table, when the groundwater is initially contaminated and the SUZ is initially free of contamination. The model is simply characterised by: the matrix-fracture porosity ratio, the matrix block geometry, and a characteristic diffusion time. De-saturation and re-saturation of fractures is handled by a new approximation method. Contaminant accumulates in the upper levels of the SUZ, where it is less accessible to mobile groundwater, and acts as a persistent secondary source of contamination once the original source of contamination has been removed or has become depleted. The 'SUZ process' first attenuates the progress of contaminants in groundwater, and subsequently controls the slow release of contamination back to the mobile groundwater, thus prolonging the duration of groundwater contamination by many years. The SUZ process should operate in any fractured, micro-porous lithology e.g. fractured clays and mudstones, making this approach widely applicable.     

Studies on plant-mediated fate of the explosives RDX and HMX

The fate of the explosives RDX and HMX on exposure to plants was investigated in 'natural' aquatic systems of Myriophyllum aquaticum for 16 days, and in axenic hairy root cultures of Catharanthus roseus for > or = 9 weeks. Exposure levels were: HMX, 5 mg/l; and RDX, approximately 8 mg/l. Exposure outcomes observed include: HMX, no transformation by aquatic plants, and minimal biological activity by axenic roots; and RDX, removal by both plant systems. In the case of RDX exposure to axenic roots, since 14C-RDX was included, removal was confirmed by the accumulation of 14C-label in the biomass. The intracellular 14C-label in these RDX studies was detected in two forms: intact RDX and bound unknown(s).     

UV-VIS spectroscopy of 2,4,6-trinitrotoluene-hydroxide reaction

Contamination of groundwater, surface water and soil by explosives has occurred at military sites throughout the world as a result of manufacture of explosive compounds, assembly of munitions, and deployment of explosives containing devices. Due to the adverse effects of explosives on humans and other natural receptors, a low cost means of decontaminating these areas of contamination is needed. Base-induced transformation of explosives has shown promise as a rapid, low cost, and minimally resource-intensive technology for detoxifying explosives in soil and water. In order to understand the reaction mechanism, a reaction mixture of 2:1:1 (water:2,4,6-trinitrotoluene (TNT):1 N KOH) was analyzed by UV/VIS spectrometry from 190 to 1,100 nm. Time course measurements were conducted at 25, 20, 15, and 12 degrees C. A factor analysis program was used to analyze the spectral data. Principal component analysis indicated that six principal components explained the spectra to within experimental error, with four factors explaining the majority of the variance. Test spectral vectors for four components were developed, including TNT, two intermediates, and the final product, and were tested against the abstract vectors. Two possible reaction mechanisms were suggested and tested to explain the spectral data.     

Assessing performance and closure for soil vapor extraction: integrating vapor discharge and impact to groundwater quality

Soil vapor extraction (SVE) is typically effective for removal of volatile contaminants from higher-permeability portions of the vadose zone. However, contamination in lower-permeability zones can persist due to mass transfer processes that limit the removal effectiveness. After SVE has been operated for a period of time and the remaining contamination is primarily located in lower-permeability zones, the remedy performance needs to be evaluated to determine whether the SVE system should be optimized, terminated, or transitioned to another technology to replace or augment SVE. Numerical modeling of vapor-phase contaminant transport was used to investigate the correlation between measured vapor-phase mass discharge, MF(r), from a persistent, vadose-zone contaminant source and the resulting groundwater contaminant concentrations. This relationship was shown to be linear, and was used to directly assess SVE remediation progress over time and to determine the level of remediation in the vadose zone necessary to protect groundwater. Although site properties and source characteristics must be specified to establish a unique relation between MF(r) and the groundwater contaminant concentration, this correlation provides insight into SVE performance and support for decisions to optimize or terminate the SVE operation or to transition to another type of treatment.     

Analysis of groundwater contamination using concentration–time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume

When only few monitoring wells are available to assess the extent and level of groundwater contamination, inversion of concentration breakthrough curves acquired during an integral pumping test can be used as an alternative quantification method. The idea is to use concentration–time series recorded during integral pumping tests through an inversion technique to estimate contaminant mass fluxes crossing a control plane. In this paper, we examine how a longitudinal concentration gradient along a contaminant plume length scale affects the estimated inversed-concentration distribution and its associated mass flux. The analytically inversed-concentration distribution at the imaginary control plane (ICP) is compared to a numerically generated concentration distribution, treating the latter one as a “real contaminant plume” characterized by the presence of a longitudinal concentration gradient. It is found that the analytically inversed-concentration can lead to overestimation or underestimation of concentration distribution values depending on the transport time period and dispersivity values. At lower dispersivity values, with shorter transport time periods, the analytically inversed-concentration distribution overestimates the “real” concentration distribution.A better fit of the estimated concentration distribution to the “real” one is observed when the transport time period increases, i.e. when the advective front has already crossed the ICP. However, for higher dispersivity values, underestimation of the real concentration distribution is observed. Deviation of the inversed-concentration distribution from the “real” one is assessed for a site-specific concentration gradient term. A concentration gradient adjusted contaminant mass flux is thus formulated to evaluate groundwater contamination levels at a given time period through an ICP. This concentration gradient ratio can indicate whether the ICP is well positioned to evaluate accurately contaminant mass fluxes which are representative of groundwater contamination levels.     

The mathematical modelling of total nitrogen contamination transport via groundwater from a sugar factory in Eskisehir, Turkey

In this study, mathematical modelling of the total nitrogen contamination transport in a porous medium was evaluated in order to determine the potential groundwater pollution caused by a sugar factory in the Eskisehir region of Turkey. Analytical solutions of mathematical modelling were performed to show graphically the distributions of contaminant concentrations. Multiflow computer programming was used to determine the distribution of contaminant concentrations with respect to time and distance. The distribution distance of contaminant concentrations was determined at any time interval. From this study, the potential pollution of groundwater can be effectively estimated. Prediction of pollution by means of the model will help to form the future predictions of water resource management.     

Analysis of groundwater contamination using concentration-time series recorded during an integral pumping test: bias introduced by strong concentration gradients within the plume

When only few monitoring wells are available to assess the extent and level of groundwater contamination, inversion of concentration breakthrough curves acquired during an integral pumping test can be used as an alternative quantification method. The idea is to use concentration-time series recorded during integral pumping tests through an inversion technique to estimate contaminant mass fluxes crossing a control plane. In this paper, we examine how a longitudinal concentration gradient along a contaminant plume length scale affects the estimated inversed-concentration distribution and its associated mass flux. The analytically inversed-concentration distribution at the imaginary control plane (ICP) is compared to a numerically generated concentration distribution, treating the latter one as a "real contaminant plume" characterized by the presence of a longitudinal concentration gradient. It is found that the analytically inversed-concentration can lead to overestimation or underestimation of concentration distribution values depending on the transport time period and dispersivity values. At lower dispersivity values, with shorter transport time periods, the analytically inversed-concentration distribution overestimates the "real" concentration distribution. A better fit of the estimated concentration distribution to the "real" one is observed when the transport time period increases, i.e. when the advective front has already crossed the ICP. However, for higher dispersivity values, underestimation of the real concentration distribution is observed. Deviation of the inversed-concentration distribution from the "real" one is assessed for a site-specific concentration gradient term. A concentration gradient adjusted contaminant mass flux is thus formulated to evaluate groundwater contamination levels at a given time period through an ICP. This concentration gradient ratio can indicate whether the ICP is well positioned to evaluate accurately contaminant mass fluxes which are representative of groundwater contamination levels.     

Transport of reactive colloids and contaminants in groundwater: effect of nonlinear kinetic interactions

Transport of reactive colloids in groundwater may enhance the transport of contaminants in groundwater. Often, the interpretation of results of transport experiments is not a simple task as both reactions of colloids with the solid matrix and reactions of contaminants with the solid matrix and mobile and immobile colloids may be time dependent and nonlinear. Further colloid transport properties may differ from solute transport properties. In this paper, a one-dimensional model for coupled and contaminant in a porous medium (COLTRAP) is presented together with simulation results. Calculated breakthrough curves (BTC's) during contamination and decontamination show systematically the effect of nonlinear and kinetic interactions on contaminant transport in the presence of reactive colloids, and the effect of colloid transport properties that differ from solute transport properties. It is shown that in case of linear kinetic reactions, the rate of exchange of mobile and immobile colloids have a large impact on the shape of BTC's even if the solid matrix is saturated with respect to colloids. BTC's during the contamination and decontamination phase have identical shapes in this case. Moreover, the slow reactions of contaminants and colloids may lead to unretarded breakthrough of contaminants. Independent of reaction rates, nonlinear reactions lead to BTC's that are steeper during contamination than in the linear case. A characteristic aspect of nonlinear sorption is that shapes of BTC's differ during the contamination and decontamination phase. It has been observed that shapes of some of the simulated adsorption and desorption curves are similar as shapes found in experiments reported in literature. This stresses the importance of incorporating both kinetics and nonlinearity in models for coupled colloid and contaminant transport and the capability of COLTRAP to interpret experimental results. Finally, to figure out whether nonlinear processes play a role, it is very important to consider both contamination and decontamination in transport experiments.     

Optimization-Based Atmospheric Plume Source Identification

This work applies optimization and an Eulerian inversion approach presented by Bagtzoglou and Baun in 2005 in order to reconstruct contaminant plume time histories and to identify the likely source of atmospheric contamination using data from a real test site for the first time. Present-day distribution of an atmospheric contaminant plume as well as data points reflecting the plume history allow the reconstruction and provide the plume velocity, distribution, and probable source. The method was tested to a hypothetical case and with data from the Forest Atmosphere Transfer and Storage (FACTS) experiment in the Duke experimental forest site. In the scenarios presented herein, as well as in numerous cases tested for verification purposes, the model conserved mass, successfully located the peak of the plume, and managed to capture the motion of the plume well but underestimated the contaminant peak.     

One-electron standard reduction potentials of nitroaromatic and cyclic nitramine explosives

Extensive studies have been conducted in the past decades to predict the environmental abiotic and biotic redox fate of nitroaromatic and nitramine explosives. However, surprisingly little information is available on one-electron standard reduction potentials (E^o(R-NO₂/R-NO₂⁻)). The E^o(R-NO₂/R-NO₂⁻) is an essential thermodynamic parameter for predicting the rate and extent of reductive transformation for energetic residues. In this study, experimental (linear free energy relationships) and theoretical (ab initio calculation) approaches were employed to determine E^o(R-NO₂/R-NO₂⁻) for nitroaromatic, (caged) cyclic nitramine, and nitroimino explosives that are found in military installations or are emerging contaminants. The results indicate a close agreement between experimental and theoretical E^o(R-NO₂/R-NO₂⁻) and suggest a key trend: E^o(R-NO₂/R-NO₂⁻) value decreases from di- and tri-nitroaromatic (e.g., 2,4-dinitroanisole) to nitramine (e.g., RDX) to nitroimino compound (e.g., nitroguanidine). The observed trend in E^o(R-NO₂/R-NO₂⁻) agrees with reported rate trends for reductive degradation, suggesting a thermodynamic control on the reduction rate under anoxic/suboxic conditions.     

Trends in spatial and temporal levels of persistent organic pollutants in Lake Erie sediments

A Lake Erie sediment survey was conducted in 1997 to characterize spatial trends in contamination, and for comparison with a 1971 survey to assess any changes in environmental quality since the advent of measures to reduce contaminant sources. Contaminant data for some contaminant classes in 1971 was based on analysis of frozen archived samples, thereby allowing a direct comparison between surficial sediment contamination in 1971 and 1997 based on modern analytical methods. Lake-wide contaminant concentrations of polychlorinated biphenyls (PCBs) and organochlorine compounds including hexachlorobenzene and the DDT compounds in 1997 were substantially lower, compared to levels in samples collected in 1971. Lake-wide average sediment PCB levels were found to have decreased roughly 70% from 136 ng/g in 1971 to 43 ng/g in 1997. Similarly, reductions in other compound classes ranged from 40% to 80%. In 1997, some individual contaminants classes including PCBs exhibited a spatial trend toward increasing sediment contamination from the eastern basin to the western basin, and from the north-central basin to the south-central basin. Levels of organic contaminants in sediments in some areas of Lake Erie still exceeded the strictest Canadian Federal and Ontario Provincial guidelines. However, exceedances of guidelines describing contaminated environments in 1997 were predominantly restricted to the western basin and near-shore sites in the southern part of the central basin.