Impact of transverse and longitudinal dispersion on first-order degradation rate constant estimation

A two-dimensional analytical model is employed for estimating the first-order degradation rate constant of hydrophobic organic compounds (HOCs) in contaminated groundwater under steady-state conditions. The model may utilize all aqueous concentration data collected downgradient of a source area, but does not require that any data be collected along the plume centerline. Using a least squares fit of the model to aqueous concentrations measured in monitoring wells, degradation rate constants were estimated at a former manufactured gas plant (FMGP) site in the Midwest U.S. The estimated degradation rate constants are 0.0014, 0.0034, 0.0031, 0.0019, and 0.0053 day(-1) for acenaphthene, naphthalene, benzene, ethylbenzene, and toluene, respectively. These estimated rate constants were as low as one-half those estimated with the one-dimensional (centerline) approach of Buscheck and Alcantar [Buscheck, T.E., Alcantar, C.M., 1995. Regression techniques and analytical solutions to demonstrate intrinsic bioremediation. In: Hinchee, R.E., Wilson, J.T., Downey, D.C. (Eds.), Intrinsic Bioremediation, Battelle Press, Columbus, OH, pp. 109-116] which does not account for transverse dispersivity. Varying the transverse and longitudinal dispersivity values over one order of magnitude for toluene data obtained from the FMGP site resulted in nearly a threefold variation in the estimated degradation rate constant-highlighting the importance of reliable estimates of the dispersion coefficients for obtaining reasonable estimates of the degradation rate constants. These results have significant implications for decision making and site management where overestimation of a degradation rate may result in remediation times and bioconversion factors that exceed expectations. For a complex source area or non-steady-state plume, a superposition of analytical models that incorporate longitudinal and transverse dispersion and time may be used at sites where the centerline method would not be applicable.     

Flux-based assessment at a manufacturing site contaminated with trichloroethylene

Groundwater and contaminant fluxes were measured, using the passive flux meter (PFM) technique, in wells along a longitudinal transect passing approximately through the centerline of a trichloroethylene (TCE) plume at a former manufacturing plant located in the Midwestern US. Two distinct zones of hydraulic conductivity were identified from the measured groundwater fluxes; a 6-m-thick upper zone ( approximately 7 m to 13 m below the ground surface or bgs) with a geometric mean Darcy flux (q(0)) of 2 cm/day, and a lower zone ( approximately 13 m to 16.5m bgs) with a q(0) approximately 15 cm/day; this important hydrogeologic feature significantly impacts any remediation technology used at the site. The flux-averaged TCE concentrations estimated from the PFM results compared well with existing groundwater monitoring data. It was estimated that at least 800 kg of TCE was present in the source zone. The TCE mass discharge across the source control plane (85 m x 38 m) was used to estimate the "source strength" ( approximately 365 g/day), while mass discharges across multiple down-gradient control planes were used to estimate the plume-averaged, TCE degradation rate constant (0.52 year(-1)). This is close to the rate estimated using the conventional centerline approach (0.78 year(-1)). The mass discharge approach provides a more robust and representative estimate than the centerline approach since the latter uses only data from wells along the plume centerline while the former uses all wells in the plume.     

Estimation of Methyl tert -Butyl Ether Plume Length Using the Domenico Analytical Model

Methyl tert -butyl ether (MTBE) plume is controlled by many factors, primarily by groundwater flow velocity, dispersion, natural attenuation. This study employed an analytical model introduced by Domemico (1987, J. Hydrol 91 , 49-58.) to describe the MTBE concentration distribution horizontal pattern and estimated the MTBE plume length. The model was applied to 90 leaking underground storage tank cases in Los Angeles, CA, U.S.A. The analytical model was calibrated with field data for each ease using a Microsoft Excel spreadsheet program. Methyl tert -butyl ether concentrations in one source monitoring well and one to two downgradient centerline monitoring wells were used for each case study. When the centerline well is not available, the closest off-centerline wells were projected to the centerline using an ellipse trigonometry method. The model parameter values for longitudinal dispersivity, groundwater velocity, and degradation rate constant were calibrated using the field data and then used to estimate the maximum distance between source well and the plume edge. This study demonstrates that the Domenico model can be applied to MTBE plume investigation when adequate field data are available. The correlation coefficients calculated based on the results of the 90 case studies indicate that MTBE plume length has a poor correlation with MTBE concentration at the source well, and a moderate negative correlation with the degradation rate constant ( m 0.65) and u / v ratio (0.64). Furthermore, MTBE plume length has a poor correlation with the longitudinal dispersivity ( m 0.4), hydraulic gradient ( m 0.1), and groundwater velocity (0.17).     

Assessing the natural attenuation of organic contaminants in aquifers using plume-scale electron and carbon balances: model development with analysis of uncertainty and parameter sensitivity

A quantitative methodology is described for the field-scale performance assessment of natural attenuation using plume-scale electron and carbon balances. This provides a practical framework for the calculation of global mass balances for contaminant plumes, using mass inputs from the plume source, background groundwater and plume residuals in a simplified box model. Biodegradation processes and reactions included in the analysis are identified from electron acceptors, electron donors and degradation products present in these inputs. Parameter values used in the model are obtained from data acquired during typical site investigation and groundwater monitoring studies for natural attenuation schemes. The approach is evaluated for a UK Permo-Triassic Sandstone aquifer contaminated with a plume of phenolic compounds. Uncertainty in the model predictions and sensitivity to parameter values was assessed by probabilistic modelling using Monte Carlo methods. Sensitivity analyses were compared for different input parameter probability distributions and a base case using fixed parameter values, using an identical conceptual model and data set. Results show that consumption of oxidants by biodegradation is approximately balanced by the production of CH4 and total dissolved inorganic carbon (TDIC) which is conserved in the plume. Under this condition, either the plume electron or carbon balance can be used to determine contaminant mass loss, which is equivalent to only 4% of the estimated source term. This corresponds to a first order, plume-averaged, half-life of > 800 years. The electron balance is particularly sensitive to uncertainty in the source term and dispersive inputs. Reliable historical information on contaminant spillages and detailed site investigation are necessary to accurately characterise the source term. The dispersive influx is sensitive to variability in the plume mixing zone width. Consumption of aqueous oxidants greatly exceeds that of mineral oxidants in the plume, but electron acceptor supply is insufficient to meet the electron donor demand and the plume will grow. The aquifer potential for degradation of these contaminants is limited by high contaminant concentrations and the supply of bioavailable electron acceptors. Natural attenuation will increase only after increased transport and dilution.     

Rebound of a coal tar creosote plume following partial source zone treatment with permanganate

The long-term management of dissolved plumes originating from a coal tar creosote source is a technical challenge. For some sites stabilization of the source may be the best practical solution to decrease the contaminant mass loading to the plume and associated off-site migration. At the bench-scale, the deposition of manganese oxides, a permanganate reaction byproduct, has been shown to cause pore plugging and the formation of a manganese oxide layer adjacent to the non-aqueous phase liquid creosote which reduces post-treatment mass transfer and hence mass loading from the source. The objective of this study was to investigate the potential of partial permanganate treatment to reduce the ability of a coal tar creosote source zone to generate a multi-component plume at the pilot-scale over both the short-term (weeks to months) and the long-term (years) at a site where there is >10 years of comprehensive synoptic plume baseline data available. A series of preliminary bench-scale experiments were conducted to support this pilot-scale investigation. The results from the bench-scale experiments indicated that if sufficient mass removal of the reactive compounds is achieved then the effective solubility, aqueous concentration and rate of mass removal of the more abundant non-reactive coal tar creosote compounds such as biphenyl and dibenzofuran can be increased. Manganese oxide formation and deposition caused an order-of-magnitude decrease in hydraulic conductivity. Approximately 125 kg of permanganate were delivered into the pilot-scale source zone over 35 days, and based on mass balance estimates <10% of the initial reactive coal tar creosote mass in the source zone was oxidized. Mass discharge estimated at a down-gradient fence line indicated >35% reduction for all monitored compounds except for biphenyl, dibenzofuran and fluoranthene 150 days after treatment, which is consistent with the bench-scale experimental results. Pre- and post-treatment soil core data indicated a highly variable and random spatial distribution of mass within the source zone and provided no insight into the mass removed of any of the monitored species. The down-gradient plume was monitored approximately 1, 2 and 4 years following treatment. The data collected at 1 and 2 years post-treatment showed a decrease in mass discharge (10 to 60%) and/or total plume mass (0 to 55%); however, by 4 years post-treatment there was a rebound in both mass discharge and total plume mass for all monitored compounds to pre-treatment values or higher. The variability of the data collected was too large to resolve subtle changes in plume morphology, particularly near the source zone, that would provide insight into the impact of the formation and deposition of manganese oxides that occurred during treatment on mass transfer and/or flow by-passing. Overall, the results from this pilot-scale investigation indicate that there was a significant but short-term (months) reduction of mass emanating from the source zone as a result of permanganate treatment but there was no long-term (years) impact on the ability of this coal tar creosote source zone to generate a multi-component plume.     

Uncertainty assessment of contaminant plume length estimates in heterogeneous aquifers

The Virtual Aquifer approach is used in this study to assess the uncertainty involved in the estimation of contaminant plume lengths in heterogeneous aquifers. Contaminant plumes in heterogeneous two-dimensional conductivity fields and subject to first order and Michaelis-Menten (MM) degradation kinetics are investigated by the center line method. First order degradation rates and plume lengths are estimated from point information obtained along the plume center line. Results from a Monte-Carlo investigation show that the estimated rate constant is highly uncertain and biased towards overly high values. Uncertainty and bias amplify with increasing heterogeneity up to maximum values of one order of magnitude. Calculated plume lengths reflect this uncertainty and bias. On average, plume lengths are estimated to about 50% of the true plume length. When plumes subject to MM degradation kinetics are investigated by using a first order rate law, an additional error is introduced and uncertainty as well as bias increase, causing plume length estimates to be less than 40% of the true length. For plumes with MM degradation kinetics, therefore, a regression approach is used which allows the determination of the MM parameters from center line data. Rate parameters are overestimated by a factor of two on average, while plume length estimates are about 80% of the true length. Plume lengths calculated using the MM parameters are thus closer to the correct length, as compared to the first order approximation. This approach is therefore recommended if field data collected along the center line of a plume give evidence of MM kinetics.     

Natural attenuation of trichloroethene and its degradation products at a lake-shore site

Subsurface contamination by trichloroethene (TCE) was detected at a Michigan National Priorities List (NPL) site in 1982. The TCE plume resulted from the disposal of spent solvent and other chemicals at an industrial facility located in the eastern shore of Lake Michigan. TCE degradation products of three dichloroethene (DCE) isomers, vinyl chloride (VC) and ethene were present. The plume was depleted of oxygen and methanogenic at certain depths. Transects of the plume were sampled by slotted auger borings the year after the TCE plume was first discovered. Water samples were also taken from lake sediments to a depth of 12 m about 100 m offshore. Later samples were taken along the shoreline of the lake with a hand-driven probe. Later in 1998 water was taken from sediments about 3-m from the shoreline. The average concentration of each chemical and net apparent base coefficient between appropriate pairs of transects between the lower site and lakeshore were calculated. Loss rates were then calculated from an analytical solution of the two-dimensional advective-dispersive-reactive transport equation. Net apparent rate coefficients and a set of coupled reaction rate equations were used to extract the apparent loss coefficients. This study showed the field evidence for natural attenuation of TCE.     

Field demonstration and evaluation of the Passive Flux Meter on a CAH groundwater plume

This study comprises the first application of the Passive Flux Meter (PFM) for the measurement of chlorinated aliphatic hydrocarbon (CAH) mass fluxes and Darcy water fluxes in groundwater at a European field site. The PFM was originally developed and applied to measurements near source zones. The focus of the PFM is extended from near source to plume zones. For this purpose, 48 PFMs of 1.4 m length were constructed and installed in eight different monitoring wells in the source and plume zone of a CAH-contaminated field site located in France. The PFMs were retrieved, sampled, and analyzed after 3 to 11 weeks of exposure time, depending on the expected contaminant flux. PFM evaluation criteria include analytical, technical, and practical aspects as well as conditions and applicability. PFM flux data were compared with so-called traditional soil and groundwater concentration data obtained using active sampling methods. The PFMs deliver reasonable results for source as well as plume zones. The limiting factor in the PFM applicability is the exposure time together with the groundwater flux. Measured groundwater velocities at the field site range from 2 to 41 cm/day. Measured contaminant flux data raise up to 13 g/m²/day for perchloroethylene in the plume zone. Calculated PFM flux averaged concentration data and traditional concentration data were of similar magnitude for most wells. However, both datasets need to be compared with reservation because of the different sampling nature and time. Two important issues are the PFM tracer loss during installation/extraction and the deviation of the groundwater flow field when passing the monitoring well and PFM. The demonstration of the PFM at a CAH-contaminated field site in Europe confirmed the efficiency of the flux measurement technique for source as well as plume zones. The PFM can be applied without concerns in monitoring wells with European standards. The acquired flux data are of great value for the purpose of site characterization and mass discharge modeling, and can be used in combination with traditional soil and groundwater sampling methods.     

Natural attenuation of chlorinated ethene compounds: model development and field-scale application at the Dover site

A multi-dimensional and multi-species reactive transport model was developed to aid in the analysis of natural attenuation design at chlorinated solvent sites. The model can simulate several simultaneously occurring attenuation processes including aerobic and anaerobic biological degradation processes. The developed model was applied to analyze field-scale transport and biodegradation processes occurring at the Area-6 site in Dover Air Force Base, Delaware. The model was calibrated to field data collected at this site. The calibrated model reproduced the general groundwater flow patterns, and also, it successfully recreated the observed distribution of tetrachloroethene (PCE), trichloroethene (TCE), dichloroethylene (DCE), vinyl chloride (VC) and chloride plumes. Field-scale decay rates of these contaminant plumes were also estimated. The decay rates are within the range of values that were previously estimated based on lab-scale microcosm and field-scale transect analyses. Model simulation results indicated that the anaerobic degradation rate of TCE, source loading rate, and groundwater transport rate are the important model parameters. Sensitivity analysis of the model indicated that the shape and extent of the predicted TCE plume is most sensitive to transmissivity values. The total mass of the predicted TCE plume is most sensitive to TCE anaerobic degradation rates. The numerical model developed in this study is a useful engineering tool for integrating field-scale natural attenuation data within a rational modeling framework. The model results can be used for quantifying the relative importance of various simultaneously occurring natural attenuation processes.     

Cost-effective sampling network design for contaminant plume monitoring under general hydrogeological conditions

A new simulation-optimization methodology is developed for cost-effective sampling network design associated with long-term monitoring of large-scale contaminant plumes. The new methodology is similar in concept to the one presented by Reed et al. (Reed, P.M., Minsker, B.S., Valocchi, A.J., 2000a. Cost-effective long-term groundwater monitoring design using a genetic algorithm and global mass interpolation. Water Resour. Res. 36 (12), 3731-3741) in that an optimization model based on a genetic algorithm is coupled with a flow and transport simulator and a global mass estimator to search for optimal sampling strategies. However, this study introduces the first and second moments of a three-dimensional contaminant plume as new constraints in the optimization formulation, and demonstrates the proposed methodology through a real-world application. The new moment constraints significantly increase the accuracy of the plume interpolated from the sampled data relative to the plume simulated by the transport model. The plume interpolation approaches employed in this study are ordinary kriging (OK) and inverse distance weighting (IDW). The proposed methodology is applied to the monitoring of plume evolution during a pump-and-treat operation at a large field site. It is shown that potential cost savings up to 65.6% may be achieved without any significant loss of accuracy in mass and moment estimations. The IDW-based interpolation method is computationally more efficient than the OK-based method and results in more potential cost savings. However, the OK-based method leads to more accurate mass and moment estimations. A comparison of the sampling designs obtained with and without the moment constraints points to their importance in ensuring a robust long-term monitoring design that is both cost-effective and accurate in mass and moment estimations. Additional analysis demonstrates the sensitivity of the optimal sampling design to the various coefficients included in the objective function of the optimization model.     

Groundwater-surface water interaction and its role on TCE groundwater plume attenuation

A field investigation of a TCE plume in a surficial sand aquifer shows that groundwater-surface water interactions strongly influence apparent plume attenuation. At the site, a former industrial facility in Connecticut, depth-discrete monitoring along three cross-sections (transects) perpendicular to groundwater flow shows a persistent VOC plume extending 700 m from the DNAPL source zone to a mid-size river. Maximum TCE concentrations along a transect 280 m from the source were in the 1000s of microg/L with minimal degradation products. Beyond this, the land surface drops abruptly to a lower terrace where a shallow pond and small streams occur. Two transects along the lower terrace, one midway between the facility and river just downgradient of the pond and one along the edge of the river, give the appearance that the plume has strongly attenuated. At the river, maximum TCE concentrations in the 10s of microg/L and similar levels of its degradation product cis-DCE show direct plume discharge from groundwater to the river is negligible. Although degradation plays a role in the strong plume attenuation, the major attenuation factor is partial groundwater plume discharge to surface water (i.e. the pond and small streams), where some mass loss occurs via water-air exchange. Groundwater and stream mass discharge estimates show that more than half of the plume mass discharge crossing the first transect, before surface water interactions occur, reaches the river directly via streamflow, although river concentrations were below detection due to dilution. This study shows that groundwater and surface water concentration measurements together provide greater confidence in identifying and quantifying natural attenuation processes at this site, rather than groundwater measurements alone.     

Assessing impacts of partial mass depletion in DNAPL source zones: II. Coupling source strength functions to plume evolution

Analytical solutions, describing the time-dependent DNAPL source-zone mass and contaminant discharge rate, derived previously in Part I [Falta, R.W., Rao, P.S., Basu, N., this issue. Assessing the impacts of partial mass depletion in DNAPL source zones: I. Analytical modeling of source strength functions and plume response. J. Contam. Hydrol.] are used as a flux-boundary condition in a semi-analytical contaminant transport model. These analytical solutions assume a power relationship between the flow-averaged source concentration, and the source DNAPL mass; the empirical exponent (gamma) is a function of the flow field heterogeneity, DNAPL architecture, and the correlation between them. The DNAPL source strength terms can account for partial source remediation, either at time zero, or at some later time after the DNAPL release. The transport model considers advection, retardation, three-dimensional dispersion, and sequential first-order decay/production of several species. A separate solution is used to compute the time-dependent mass of each contaminant in the plume. A series of examples using different values of gamma shows how the benefits of partial DNAPL source remediation can vary with site conditions. In general, when gamma>1, relatively large short-term reductions in the plume concentrations and mass occur, but the source longevity is not strongly affected. Conversely, when gamma<1, the short-term reductions in the plume concentrations and mass are smaller, but the source longevity can be greatly reduced. In either case, the source remediation effort is much more effective if it is undertaken at an early time, before much contaminant mass has entered the plume. If the remediation effort is significantly delayed, the leading parts of the plume are not affected by the source remediation, and additional control or remediation of the plume itself is required.     

Anaerobic degradation of naphthalene in a fluvial aquifer: a radiotracer study

A radiotracer study was conducted in a creosote-contaminated aquifer beneath the Fraser River, British Columbia Canada to investigate the in situ degradation of naphthalene. The groundwater is anaerobic, with abundant methane, ferrous iron and carbon dioxide. This study followed earlier work at the site where the contaminant distribution could only be explained by invoking a mass loss through degradation, even though extensive field and laboratory microcosm studies closer to the source zone onshore could not confirm degradation. Accordingly, 14C-naphthalene was injected into the aquifer offshore, further from the source zone where modeling suggested degradation was occurring. During the 230-day monitoring period, 14CO2 was detected, confirming the degradation of the radio-labeled naphthalene tracer. A zero-order degradation rate of naphthalene of 5 microg/L-day was estimated based on the decrease in 14C-naphthalene concentration with time. While the degradation pathway could not be determined from the radiotracer study alone, the geochemistry of the site suggests that either iron reduction or methanogenesis is the terminal electron accepting processes responsible for naphthalene oxidation.     

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.     

Processes controlling the distribution and natural attenuation of dissolved phenolic compounds in a deep sandstone aquifer

Processes controlling the distribution and natural attenuation (NA) of phenol, cresols and xylenols released from a former coal-tar distillation plant in a deep Triassic sandstone aquifer are evaluated from vertical profiles along the plume centerline at 130 and 350 m from the site. Up to four groups of contaminants (phenols, mineral acids, NaOH, NaCl) form discrete and overlapping plumes in the aquifer. Their distribution reflects changing source history with releases of contaminants from different locations. Organic contaminant distribution in the aquifer is determined more by site source history than degradation. Contaminant degradation at total organic carbon (TOC) concentrations up to 6500 mg l(-1) (7500 mg l(-1) total phenolics) is occurring by aerobic respiration NO3-reduction, Mn(IV)-/Fe(III)-reduction, SO4-reduction, methanogenesis and fermentation, with the accumulation of inorganic carbon, organic metabolites (4-hydroxybenzaldehyde, 4-hydroxybenzoic acid), acetate, Mn(II), Fe(II), S(-II), CH4 and H2 in the plume. Aerobic and NO3-reducing processes are restricted to a 2-m-thick plume fringe but Mn(IV)-/Fe(II)-reduction, SO4-reduction, methanogenesis and fermentation occur concomitantly in the plume. Dissolved H2 concentrations in the plume vary from 0.7 to 110 nM and acetate concentrations reach 200 mg l(-1). The occurrence of a mixed redox system and concomitant terminal electron accepting processes (TEAPs) could be explained with a partial equilibrium model based on the potential in situ free energy (deltaGr) yield for oxidation of H2 by specific TEAPs. Respiratory processes rather than fermentation are rate limiting in determining the distribution of H2 and TEAPs and H2 dynamics in this system. Most (min. 90%) contaminant degradation has occurred by aerobic and NO3-reducing processes at the plume fringe. This potential is determined by the supply of aqueous O2 and NO3 from uncontaminated groundwater, as controlled by transverse mixing, which is limited in this aquifer by low dispersion. Consumption to date of mineral oxides and SO4 is, respectively, <0.15% and 0.4% of the available aquifer capacity, and degradation using these oxidants is <10%. Fermentation is a significant process in contaminant turnover, accounting for 21% of degradation products present in the plume, and indicating that microbial respiration rates are slow in comparison with fermentation. Under present conditions, the potential for degradation in the plume is very low due to inhibitory effects of the contaminant matrix. Degradation products correspond to <22% mass loss over the life of the plume, providing a first-order plume scale half-life >140 years. The phenolic compounds are biodegradable under the range of redox conditions in the aquifer and the aquifer is not oxidant limited, but the plume is likely to be long-lived and to expand. Degradation is likely to increase only after contaminant concentrations are reduced and aqueous oxidant inputs are increased by dispersion of the plume. The results imply that transport processes may exert a greater control on the natural attenuation of this plume than aquifer oxidant availability.     

Comparing estimates of fugitive landfill methane emissions using inverse plume modeling obtained with Surface Emission Monitoring (SEM), Drone Emission Monitoring (DEM), and Downwind Plume Emission Monitoring (DWPEM)

ABSTRACT

As part of the global effort to quantify and manage anthropogenic greenhouse gas emissions, there is considerable interest in quantifying methane emissions in municipal solid waste landfills. A variety of analytical and experimental methods are currently in use for this task. In this paper, an optimization-based estimation method is employed to assess fugitive landfill methane emissions. The method combines inverse plume modeling with ambient air methane concentration measurements. Three different measurement approaches are tested and compared. The method is combined with surface emission monitoring (SEM), above ground drone emission monitoring (DEM), and downwind plume emission monitoring (DWPEM). The methodology is first trialed and validated using synthetic datasets in a hand-generated case study. A field study is also presented where SEM, DEM and DWPEM are tested and compared. Methane flux during two-days measurement campaign was estimated to be between 228 and 350 g/s depending on the type of measurements used. Compared to SEM, using unmanned aerial systems (UAS) allows for a rapid and comprehensive coverage of the site. However, as showed through this work, advancement of DEM-based methane sampling is governed by the advances that could be made in UAS-compatible measurement instrumentations. Downwind plume emission monitoring led to a smaller estimated flux compared with SEM and DEM without information about positions of major leak points in the landfill. Even though, the method is simple and rapid for landfill methane screening. Finally, the optimization-based methodology originally developed for SEM, shows promising results when it is combined with the drone-based collected data and downwind concentration measurements. The studied cases also discovered the limitations of the studied sampling strategies which is exploited to identify improvement strategies and recommendations for a more efficient assessment of fugitive landfill methane emissions.

Implications: Fugitive landfill methane emission estimation is tackled in the present study. An optimization-based method combined with inverse plume modeling is employed to treat data from surface emission monitoring, drone-based emission monitoring and downwind plume emission monitoring. The study helped revealing the advantages and the limitations of the studied sampling strategies. Recommendations for an efficient assessment of landfill methane emissions are formulated. The method trialed in this study for fugitive landfill methane emission could also be appropriate for rapid screening of analogous greenhouse gas emission hotspots.     

Simultaneous optimization of dense non-aqueous phase liquid (DNAPL) source and contaminant plume remediation

A framework is developed for simultaneous, optimal design of groundwater contaminant source removal and plume remediation strategies. The framework allows for varying degrees of effort and cost to be dedicated to source removal versus plume remediation. We have accounted for the presence of physical heterogeneity in the DNAPL source, since source heterogeneity controls mass release into the plume and the efficiency of source removal efforts. We considered high and low estimates of capital and operating costs for chemical flushing removal of the source, since these are expected to vary form site to site. Using the lower chemical flushing cost estimates, it is found that the optimal allocation of funds to source removal or plume remediation is sensitive to the degree of heterogeneity in the source. When the time elapsed between the source release and the implementation of remediation was varied, it was found that, except for the longest elapsed time (50,000 days), a combination of partial source removal and plume remediation was most efficient. When first-order, dissolved contaminant degradation was allowed, source removal was found to be unnecessary for the cases where the degradation rate exceeded intermediate values of the first-order rate constant. Finally, it was found that source removal became more necessary as the degree of aquifer heterogeneity increased.     

Microbiological analysis of multi-level borehole samples from a contaminated groundwater system

A range of bacteriological, geochemical process-related and molecular techniques have been used to assess the microbial biodegradative potential in groundwater contaminated with phenol and other tar acids. The contaminant plume has travelled 500 m from the pollutant source over several decades. Samples were obtained from the plume using a multi-level sampler (MLS) positioned in two boreholes (boreholes 59 and 60) which vertically transected two areas of the plume. Activity of the microbial community, as represented by phenol degradation potential and ability to utilise a range of substrates, was found to be influenced by the plume. Phenol degradation potential appeared to be influenced more by the concentration of the contaminants than the total bacterial cell numbers. However, in the areas of highest phenol concentration, the depression of cell numbers clearly had an effect. The types of bacteria present were assessed by culture and DNA amplification by polymerase chain reaction (PCR). Bacterial groups or processes associated with major geochemical processes, such as methanogenesis, sulphate reduction and denitrification, that have the potential to drive contaminant degradation, were detected at various borehole levels. A comparative molecular analysis of the microbial community between samples obtained from the MLS revealed the microbial community was diverse. The examination of microbial activity complemented those results obtained through chemical analysis, and when combined with hydrological data, showed that MLS samples provided a realistic profile of plume effects and could be related to the potential for natural attenuation of the site.     

Plume persistence caused by back diffusion from thin clay layers in a sand aquifer following TCE source-zone hydraulic isolation

This paper concludes that back diffusion from one or a few thin clayey beds in a sand aquifer can cause contaminant persistence above MCLs in a sand aquifer long after the source zone initially causing the plume is isolated or removed. This conclusion is based on an intensive case study of a TCE contaminated site in Florida, with the processes evaluated using numerical modeling. At this site, the TCE DNAPL zone formed decades ago, and was hydraulically isolated by means of an innovative system performing groundwater extraction, treatment and re-injection. Treated water is re-injected in a row of injection wells situated a short distance downgradient of the extraction wells, creating a clean-water displacement front to efficiently flush the downgradient plume. This scheme avoids the creation of stagnation zones typical of most groundwater pump-and-treat systems, thereby minimizing the time for aquifer flushing and therefore downgradient cleanup. The system began operation in August 2002 and although the performance monitoring shows substantial declines in concentrations, detectable levels of TCE and degradation products persist downgradient of the re-injection wells, long after the TCE should have disappeared based on calculations assuming a nearly homogenous sand aquifer. Three hypotheses were assessed for this plume persistence: 1) incomplete source-zone capture, 2) DNAPL occurrence downgradient of the re-injection wells, and 3) back diffusion from one or more thin clay beds in the aquifer. After careful consideration, the first two hypotheses were eliminated, leaving back diffusion as the only plausible hypothesis, supported by detailed measurements of VOC concentrations within and near the clay beds and also by numerical model simulations that closely represent the field site hydrogeologic conditions. The model was also used to simulate a more generalized, hypothetical situation where more thin clayey beds occur in a sand aquifer with an underlying aquitard. While there is no doubt that DNAPL source mass reduction can eventually improve downgradient groundwater quality, the magnitude and time scale over which the improvement occurs is the major uncertainty given current characterization approaches. This study shows that even one thin clay bed, less than 0.2 m thick, can cause plume persistence due to back diffusion for several years or even decades after the flux from the source is completely isolated. Thin clay beds, which have a large storage capacity for dissolved and sorbed contaminant mass, are common in many types of sandy aquifers. However, without careful inspection of continuous cores and sampling, such thin clay beds, and their potential for causing long-term back-diffusion effects, can easily go unnoticed during site characterization.     

An analytical quantification of mass fluxes and natural attenuation rate constants at a former gasworks site

A new integral groundwater investigation approach was used for the first time to quantify natural attenuation rates at field scale. In this approach, pumping wells positioned along two control planes were operated at distances of 140 and 280 m downstream of a contaminant source zone at a former gasworks site polluted with BTEX- (benzene, toluene, ethyl-benzene, o-, p-xylene) and PAH- (polycyclic aromatic hydrocarbons) compounds. Based on the quantified changes in total contaminant mass fluxes between the control planes, first-order natural attenuation rate constants could be estimated. For BTEX-compounds, these ranged from 1.4e-02 to 1.3e-01 day(-1) whereas for PAH-compounds natural attenuation rate constants of 3.7e-04 to 3.1e-02 day(-1) were observed. Microbial degradation activity at the site was indicated by an increase in dissolved iron mass flux and a reduction in sulphate mass flux between the two investigated control planes. In addition to information about total contaminant mass fluxes and average concentrations, an analysis of the concentration-time series measured at the control planes also allowed to semi-quantitatively delineate the aquifer regions most likely contaminated by the BTEX- and PAH-compounds.