Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m,p-xylene, and o-xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o-xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m,p-xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o-xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.     

Estimation of Efficiency of Processing Soil Samples for Pesticide Residues Analysis

The efficiency of four sample processing methods was tested with eight different types of soils representing the major proportion of cultivated soils. The principle of sampling constant was applied for characterizing the efficiency of the procedures and testing the well-mixed status of the prepared soil. The test material was 14C-labeled atrazine that enabled keeping the random error of analyses ≤ about 1%. Adding water to the soil proved to be the most efficient and generally applicable procedure resulting in about 6% relative sample processing uncertainty for 20 g test portions. The expected error is inversely proportional to the mass of test portion. Smashing and manual mixing of soil resulted in about four times higher uncertainty than mixing with water. Grinding of soil is applicable for dry soils only, but the test procedure applied was not suitable for estimating a typical uncertainty of processing dry soil samples. Adding dry ice did not improve the efficiency of sample processing.     

Loss of volatile hydrocarbons from an LNAPL oil source

The light nonaqueous phase liquid (LNAPL) oil pool in an aquifer that resulted from a pipeline spill near Bemidji, Minnesota, was analyzed for volatile hydrocarbons (VHCs) to determine if the composition of the oil remains constant over time. Oil samples were obtained from wells at five locations in the oil pool in an anaerobic part of the glacial outwash aquifer. Samples covering a 21-year period were analyzed for 25 VHCs. Compared to the composition of oil from the pipeline source, VHCs identified in oil from wells sampled in 2008 were 13 to 64% depleted. The magnitude of loss for the VHCs analyzed was toluene?o-xylene, benzene, C(6) and C(10-12)n-alkanes>C(7)-C(9)n-alkanes>m-xylene, cyclohexane, and 1- and 2-methylnaphthalene>1,2,4-trimethylbenzene and ethylbenzene. Other VHCs including p-xylene, 1,3,5- and 1,2,3-trimethylbenzenes, the tetramethylbenzenes, methyl- and ethyl-cyclohexane, and naphthalene were not depleted during the time of the study. Water-oil and air-water batch equilibration simulations indicate that volatilization and biodegradation is most important for the C(6)-C(9)n-alkanes and cyclohexanes; dissolution and biodegradation is important for most of the other hydrocarbons. Depletion of the hydrocarbons in the oil pool is controlled by: the lack of oxygen and nutrients, differing rates of recharge, and the spatial distribution of oil in the aquifer. The mass loss of these VHCs in the 5 wells is between 1.6 and 7.4% in 29years or an average annual loss of 0.06-0.26%/year. The present study shows that the composition of LNAPL changes over time and that these changes are spatially variable. This highlights the importance of characterizing the temporal and spatial variabilities of the source term in solute-transport models.     

Progression of natural attenuation processes at a crude-oil spill site: I. Geochemical evolution of the plume

A 16-year study of a hydrocarbon plume shows that the extent of contaminant migration and compound-specific behavior have changed as redox reactions, most notably iron reduction, have progressed over time. Concentration changes at a small scale, determined from analysis of pore-water samples drained from aquifer cores, are compared with concentration changes at the plume scale, determined from analysis of water samples from an observation well network. The small-scale data show clearly that the hydrocarbon plume is growing slowly as sediment iron oxides are depleted. Contaminants, such as ortho-xylene that appeared not to be moving downgradient from the oil on the basis of observation well data, are migrating in thin layers as the aquifer evolves to methanogenic conditions. However, the plume-scale observation well data show that the downgradient extent of the Fe2+ and BTEX plume did not change between 1992 and 1995. Instead, depletion of the unstable Fe (III) oxides near the subsurface crude-oil source has caused the maximum dissolved iron concentration zone within the plume to spread at a rate of approximately 3 m/year. The zone of maximum concentrations of benzene, toluene, ethylbenzene and xylene (BTEX) has also spread within the anoxic plume. In monitoring the remediation of hydrocarbon-contaminated ground water by natural attenuation, subtle concentration changes in observation well data from the anoxic zone may be diagnostic of depletion of the intrinsic electron-accepting capacity of the aquifer. Recognition of these subtle patterns may allow early prediction of growth of the hydrocarbon plume.     

In situ measurements of volatile aromatic hydrocarbon biodegradation rates in groundwater

Benzene and alkylbenzene biodegradation rates and patterns were measured using an in situ microcosm in a crude-oil contaminated aquifer near Bemidji, Minnesota. Benzene-D6, toluene, ethylbenzene, o-, m- and p-xylenes and four pairs of C₃- and C₄-benzenes were added to an in situ microcosm and studied over a 3-year period. The microcosm allowed for a mass-balance approach and quantification of hydrocarbon biodegradation rates within a well-defined iron-reducing zone of the anoxic plume. Among the BTEX compounds, the apparent order of persistence is ethylbenzene > benzene > m,p-xylenes > o-xylene ≥ toluene. Threshold concentrations were observed for several compounds in the in situ microcosm, below which degradation was not observed, even after hundreds of days. In addition, long lag times were observed before the onset of degradation of benzene or ethylbenzene. The isomer-specific degradation patterns were compared to observations from a multi-year study conducted using data collected from monitoring wells along a flowpath in the contaminant plume. The data were fit with both first-order and Michaelis-Menten models. First-order kinetics provided a good fit for hydrocarbons with starting concentrations below 1 mg/L and Michaelis-Menten kinetics were a better fit when starting concentrations were above 1 mg/L, as was the case for benzene. The biodegradation rate data from this study were also compared to rates from other investigations reported in the literature.