Transverse Hydrodynamic Dispersion Effects on ... - ACS Publications

8 Jun 2012 - Department of Civil and Environmental Engineering, Stanford ..... In summary, ongoing degradation processes at plume fringes may be stron...
5 downloads 0 Views 2MB Size
Article pubs.acs.org/est

Transverse Hydrodynamic Dispersion Effects on Isotope Signals in Groundwater Chlorinated Solvents’ Plumes Boris M. Van Breukelen*,† and Massimo Rolle‡,# †

Faculty of Earth and Life Sciences, VU University Amsterdam, De Boelelaan 1085, NL-1081 HV Amsterdam, The Netherlands Center for Applied Geosciences, University of Tübingen, Sigwartstrasse 10, D-72076 Tübingen, Germany # Department of Civil and Environmental Engineering, Stanford University, 473 Via Ortega, 94305 Stanford, California, United States ‡

S Supporting Information *

ABSTRACT: The effects of transverse hydrodynamic dispersion on altering transformation-induced compound-specific isotope analysis (CSIA) signals within groundwater pollution plumes have been assessed with reactive transport modeling accommodating diffusion-induced isotope fractionation (DIF) and implementing different parameterizations of local transverse dispersion. The model reproduced previously published field data showing a negative carbon isotope pattern (−2 ‰) at the fringes of a nondegrading PCE plume. We extended the study to reactive transport scenarios considering vinyl chloride as a model compound and assessing, through a detailed sensitivity analysis, the coupled effects of transverse hydrodynamic dispersion (with and without DIF) and aerobic fringe degradation on the evolution of carbon and chloride isotope ratios. Transformation-induced positive isotope signals were increasingly attenuated with distance from the source and higher degradation rate. The effect of DIF on the overall isotope signal attenuation was greatest near the source and for low values of groundwater flow velocity, transverse dispersion coefficient, molecular weight, rate constant, and isotope fractionation factor, α, of the degradation reaction. Models disregarding DIF underestimate the actual α. The approximately twice larger DIF effect for chlorine than for carbon together with the low α for oxidation resulted in strong chlorine CSIA depletions for VC at the plume fringe.



INTRODUCTION Fringe-controlled biodegradation plays a pivotal role in determining the fate of groundwater pollution plumes. In particular, for plumes approaching steady-state conditions, transverse hydrodynamic dispersion determines the extent of mixing of the organic pollutants with dissolved oxidants (oxygen, nitrate, and sulfate) in the ambient groundwater and, thereby, is a key factor in plume attenuation, as shown by a number of model-based and field studies (e.g., refs 1−7). Therefore, a correct quantification of the interplay between physical transport processes and (bio)degradation reactions occurring at the plume fringe is of utmost importance for understanding and assessing natural attenuation of pollutants in aquifer systems. The analysis of the spatial and temporal evolution of isotope ratios during contaminant transport has the potential of providing additional and often decisive information on contaminant degradation. Compound-specific stable isotope analysis (CSIA) has been established as a monitoring tool in the past decade and provides a more direct approach than concentration-based monitoring to identify the occurrence of degradation of organic pollutants.8−11 Isotope ratios of widespread organic groundwater contaminants like chlorinated ethenes increase considerably over the course of degradation for most transformation pathways. In particular, for © 2012 American Chemical Society

oxidative processes of chlorinated ethenes at the plume fringe, CSIA can be applied to obtain direct information on biodegradation, and, under ideal circumstances, isotope enrichments can be linked directly to mass destruction as is typically required for monitored natural attenuation to be accepted as remediation technology. In many practical applications at contaminated field sites, the interpretation of isotopic data needs to take into account the effects of depth-averaged sampling from long well screens which attenuate degradation-induced isotope signals.12−14 Moreover, recent studies, showed that CSIA data are more affected than previously thought by physical processes like mass transfer from nonaqueous phase liquids to the aqueous phase,15 sorption,12,14,16,17 volatilization,18 dispersion,14,19,20 and diffusion,20,21 thus calling for comprehensive interpretation approaches taking into account the interplay between physical and reactive processes. In this contribution, we make a detailed, model-assisted analysis of the effects of transverse diffusive/ dispersive processes on carbon and chlorine isotope ratios at Received: Revised: Accepted: Published: 7700

March 17, 2012 June 6, 2012 June 8, 2012 June 8, 2012 dx.doi.org/10.1021/es301058z | Environ. Sci. Technol. 2012, 46, 7700−7708

Environmental Science & Technology

Article

homogeneous porous medium.36 As shown in the Supporting Information, SI, a two-dimensional steady state transport problem can be reduced to a 1-D transient one by a simple transformation of coordinates. The comparison of both numerical codes with the analytical solution (eq S4 of the SI) showed an excellent agreement over a wide range (i.e., 6 orders of magnitudes) of contaminant concentrations. The numerical simulations performed include: (i) conservative transport to interpret the high-resolution concentrations and carbon isotope vertical profiles of a PCE plume observed in a shallow sandy aquifer, (ii) a reactive transport modeling scenario investigating the effect of transverse hydrodynamic dispersion and fringe biodegradation on C and Cl isotopic signatures, and (iii) a detailed sensitivity analysis to systematically assess the influence of changes in transport and reaction parameters on isotope data and their interpretation to quantify biodegradation. Diffusion and Dispersion Coefficients. The modeling tools used in this study are able to accommodate isotopologuespecific diffusion and dispersion coefficients. Few literature data are available on the effect of different isotopic masses on the diffusion of neutral species in water. To the best of our knowledge, these properties should still be well documented for organic contaminants such as chlorinated ethenes. As a working hypothesis for the present study, we calculated the aqueous diffusion coefficients, Daq, of the different contaminants’ isotopologues according to the empirical correlation proposed by Worch 37 for organic compounds:

the fringes of chlorinated solvents’ plumes and their implications for CSIA data evaluation. Isotopic fractionation by diffusion has been studied mainly in diffusion-dominated aqueous environments for cases such as methane oxidation in marine sediments,22 denitrification in the hypolimnion of a lake,23 and chloride diffusion through clay layers for paleo groundwater studies,24,25 and in unsaturated soils for cases including gaseous CO2 transport,26 and diffusion of volatile petroleum hydrocarbons27,28 and chlorinated ethenes.29 Diffusion-induced isotope fractionation (DIF) has not yet been quantitatively investigated in advection-dominated saturated environments or more specifically groundwater pollution plume studies with the exception of LaBolle et al.21 A few recent contributions inspired us for the present work. LaBolle et al. 21 simulated significant carbon isotope enrichments of organic pollutants in thin aquifers ( 1 implying an overestimation of degradation applying the Rayleigh equation. Since we focused especially on the behavior at the plume fringe we excluded these points from the representation of Figure 3. Figure 3 shows that diffusion-induced isotope fractionation further attenuates the biodegradation-induced isotope signal (lower θ for HDIF versus HD), whereas the new parameterization for hydrodynamic dispersion, NP, gives typically lower θ values than the classical one, CP. Distance from the source and kox largely control θ for HD simulations, whereas θ is nearly insensitive to variations of the other parameters for HD simulations (Figure 3). The two aforementioned parameters show a nearly identical effect to θ suggesting a similar mechanistic explanation (Figure 3A,E). Higher degradation rates decrease θ as discussed before (Figure 3E; SI;14,19). Similarly, at larger distance from the source, mixing becomes progressively more limiting 51 and its effect on the overall transformation is well captured by the computed θ, which decreases indicating more substantial deviation from the 7705

dx.doi.org/10.1021/es301058z | Environ. Sci. Technol. 2012, 46, 7700−7708

Environmental Science & Technology

Article

Chlorine isotope fractionation factors for cDCE and VC oxidation are particularly low,45 whereas the DIF effect is approximately twice as large for chlorine than for carbon and may consequently result in negative isotope patterns in plume fringe areas. However, as the carbon isotope fractionation factors for oxidation of these lower chlorinated ethenes are relatively high, the dual carbon−chlorine isotope approach to detect possible cDCE and VC oxidation33,45 should still be useful. The chlorine versus carbon dual isotope slope for fringedegradation may, however, not show a gentle incline,33 but rather moves into the fourth quadrant on a dual isotope plot (Figure 2F). Underestimation or nondetecting biodegradation at pollution fringes based on CSIA data may occur for pollutants which have low fractionation factors for oxidation depending on the microorganisms involved (e.g., MTBE: εC = −0.2 to −2.5 ‰, ε H ≈ 0 to −100 ‰; 52−54 1,2,4trichlorobenzene: εC ≤ −0.4 ‰;55 toluene: εC = −0.4 to −6.2 ‰, but εH ≫ −2 ‰),56 especially under conditions favoring the HDIF effect. This study also pointed out the key role of diffusion on contaminant transport and degradation. Diffusive processes take place at the small pore scale but they are relevant and manifest themselves at larger macroscopic scales.57 Our results shows that a reactive transport model approach which considers DIF was required to reproduce the negative isotope patterns at the investigated field site. Moreover, a reactive transport model accounting for DIF also allows obtaining better estimates on actual degradation rates occurring within groundwater pollution plumes. Finally, this study calls for further experimental research on diffusive/ dispersive processes and their importance for an improved interpretation of isotopic signatures at contaminated field sites.

Rayleigh equation which assumes a well-mixed system. The additional effect of HDIF becomes dominant when degradation is sluggish or has progressed to limited extent when the point of observation is near the source area. For low extents of degradation, negative isotope patterns produced by HDIF are only marginally turned toward positive isotope values indicative of biodegradation. At high rates of degradation and with large travel distance all methods converge to low θ (Figure 3A,E). Variations in the degree of heterogeneity of the aquifer (σ2lnK) and v also show comparable control on θ which follows from the influence of these parameters on transverse hydrodynamic dispersion. The importance of pore diffusion over mechanical dispersion is higher at low flow velocities and in mildly heterogeneous media causing strong reductions in θ for HDIF simulations. For higher values of σ2lnK and v, CP gradually increases to the invariant θ of HD, whereas θ remains low for NP. Low ε values produce a positive but small isotope signal which easily becomes reduced due to HDIF (Figure 3F). For the base case settings, CP and NP inhibit the development of positive isotope signals for εC below roughly −2 and −4 ‰, respectively (Figure 3F; SI). Figure 3D shows that the molecular weight has no effect on HD, whereas θ becomes lower with lower MW for CP and NP due to enhanced HDIF. In summary, ongoing degradation processes at plume fringes may be strongly underestimated by the Rayleigh equation if the distance from the source zone to the monitoring well is long, degradation rates are fast, and the groundwater flow velocity, heterogeneity of the hydraulic conductivity field, and ε are low. Diffusion-induced isotope fractionation effects have the most impact at low values of the following factors: distance from the source, groundwater flow velocity, heterogeneity of the hydraulic conductivity field, degradation rate constant, kinetic isotope enrichment factor, and molecular weight. Implications for CSIA-Based Assessment and Modeling of Biodegradation. The field site model demonstration showed that diffusion-induced isotope fractionation (DIF) is a relevant process at plume fringes and could be satisfactorily simulated with the developed model, which could capture the diffusive/dispersive effects on the measured δ13C signals. Neither of the two parameterizations of hydrodynamic dispersion, CP and NP, always gave a good fit for the three investigated field locations but they allowed an explanation of the observed negative isotope patterns at the PCE plume fringe and their combination indicated the potential range of effects. The sensitivity analysis indicates that the DIF effect needs to be considered for CSIA data interpretation under a wide range of conditions. CSIA as a tool to assess degradation is limited to low molecular weight organic pollutants, which are inherently sensitive to DIF. Beside chlorine CSIA, oxygen (δ18O−NO3, δ18O −SO4), sulfur (δ34S−SO4), and bromine CSIA are more sensitive to DIF because of the isotope mass difference of two compared to one for carbon, hydrogen, and nitrogen. The DIF effect may have gone unnoticed at contaminated sites because of the lack of high-resolution monitoring devices necessary to capture the steep vertical diffusive/dispersive gradients and, where degradation does occur, as biodegradation-induced positive isotope patterns have only become reduced to less positive values. Field calibration of models which disregard DIF results in underestimation of actual isotope enrichment factors operative at field sites. This tendency becomes more important for low actual ε values.



ASSOCIATED CONTENT

S Supporting Information *

Model validation to analytical solution; outline of the fluxrelated moment analysis; parameters and results of all simulations performed for the sensitivity analysis. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*Phone: +31-20-5987393; fax: +31-20-5989940; e-mail: b.m. [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Both authors, B.M.V.B. and M.R., contributed equally to this work. B.M.V.B. was financially supported by the European Union under the 7th Framework Programme (project acronym CSI:ENVIRONMENT, contract number PITN-GA-2010264329), ESTCP project ER-201029, and Direct Funding. M.R. acknowledges the support of the Marie Curie International Outgoing Fellowship (DILREACT project) within the 7th European Community Framework Programme. We thank three anonymous reviewers for their helpful and constructive comments.



REFERENCES

(1) Chu, M.; Kitanidis, P. K.; McCarty, P. L. Modeling microbial reactions at the plume fringe subject to transverse mixing in porous media: When can the rates of microbial reaction be assumed to be

7706

dx.doi.org/10.1021/es301058z | Environ. Sci. Technol. 2012, 46, 7700−7708

Environmental Science & Technology

Article

instantaneous? Water Resour. Res. 2005, 41, W06002 DOI: 10.1029/ 2004wr003495. (2) Cirpka, O. A.; Frind, E. O.; Helmig, R. Numerical simulation of biodegradation controlled by transverse mixing. J. Contam. Hydrol. 1999, 40, 159−182. (3) Davis, G. B.; Barber, C.; Power, T. R.; Thierrin, J.; Patterson, B. M.; Rayner, J. L.; Wu, Q. L. The variability and intrinsic remediation of a BTEX plume in anaerobic sulphate-rich groundwater. J. Contam. Hydrol. 1999, 36, 265−290. (4) Prommer, H.; Anneser, B.; Rolle, M.; Einsiedl, F.; Griebler, C. Biogeochemical and isotopic gradients in a BTEX/PAH contaminant plume: Model-based interpretation of a high-resolution field data set. Environ. Sci. Technol. 2009, 43, 8206−8212. (5) Rolle, M.; Clement, T. P.; Sethi, R.; Di Molfetta, A. A kinetic approach for simulating redox-controlled fringe and core biodegradation processes in groundwater: Model development and application to a landfill site in Piedmont, Italy. Hydrol. Process. 2008, 22, 4905−4921. (6) Tuxen, N.; Albrechtsen, H. J.; Bjerg, P. L. Identification of a reactive degradation zone at a landfill leachate plume fringe using high resolution sampling and incubation techniques. J. Contam. Hydrol. 2006, 85, 179−194. (7) van Breukelen, B. M.; Griffioen, J. Biogeochemical processes at the fringe of a landfill leachate pollution plume: Potential for dissolved organic carbon, Fe(II), Mn(II), NH4, and CH4 oxidation. J. Contam. Hydrol. 2004, 73, 181−205. (8) Elsner, M. Stable isotope fractionation to investigate natural transformation mechanisms of organic contaminants: Principles, prospects and limitations. J. Environ. Monit. 2010, 12, 2005−2031. (9) Hofstetter, T. B.; Schwarzenbach, R. P.; Bernasconi, S. M. Assessing transformation processes of organic compounds using stable isotope fractionation. Environ. Sci. Technol. 2008, 42, 7737−7743. (10) Hunkeler, D.; Meckenstock, R. U.; Lollar, B. S.; Schmidt, T. C.; Wilson, J. T. a Guide for Assessing Biodegradation and Source Identification of Organic Ground Water Contaminants Using Compound Specific Isotope Analysis (CSIA); US EPA: 2008. (11) Meckenstock, R. U.; Morasch, B.; Griebler, C.; Richnow, H. H. Stable isotope fractionation analysis as a tool to monitor biodegradation in contaminated acquifers. J. Contam. Hydrol. 2004, 75, 215−255. (12) Kopinke, F. D.; Georgi, A.; Voskamp, M.; Richnow, H. H. Carbon isotope fractionation of organic contaminants due to retardation on humic substances: Implications for natural attenuation studies in aquifers. Environ. Sci. Technol. 2005, 39, 6052−6062. (13) Thullner, M.; Centler, F.; Richnow, H. H.; Fischer, A. Quantification of organic pollutant degradation in contaminated aquifers using compound specific stable isotope analysisReview of recent developments. Org. Geochem. 2012, 42, 1440−1460. (14) Van Breukelen, B. M.; Prommer, H. Beyond the Rayleigh equation: Reactive transport modeling of isotope fractionation effects to improve quantification of biodegradation. Environ. Sci. Technol. 2008, 42, 2457−2463. (15) Aeppli, C.; Berg, M.; Cirpka, O. A.; Holliger, C.; Schwarzenbach, R. P.; Hofstetter, T. B. Influence of mass-transfer limitations on carbon isotope fractionation during microbial dechlorination of trichloroethene. Environ. Sci. Technol. 2009, 43, 8813−8820. (16) Hohener, P.; Atteia, O. Multidimensional analytical models for isotope ratios in groundwater pollutant plumes of organic contaminants undergoing different biodegradation kinetics. Adv. Water Resour 2010, 33, 740−751. (17) Höhener, P.; Yu, X. Stable carbon and hydrogen isotope fractionation of dissolved organic groundwater pollutants by equilibrium sorption. J. Contam. Hydrol 2012, 129−130, 54−61. (18) Kuder, T.; Philp, P.; Allen, J. Effects of volatilization on carbon and hydrogen isotope ratios of MTBE. Environ. Sci. Technol. 2009, 43, 1763−1768. (19) Abe, Y.; Hunkeler, D. Does the Rayleigh equation apply to evaluate field isotope data in contaminant hydrogeology? Environ. Sci. Technol. 2006, 40, 1588−1596.

(20) Rolle, M.; Chiogna, G.; Bauer, R.; Griebler, C.; Grathwohl, P. Isotopic fractionation by transverse dispersion: Flow-through microcosms and reactive transport modeling study. Environ. Sci. Technol. 2010, 44, 6167−6173. (21) LaBolle, E. M.; Fogg, G. E.; Eweis, J. B.; Gravner, J.; Leaist, D. G. Isotopic fractionation by diffusion in groundwater. Water Resour. Res. 2008, 44, W07405 DOI: 10.1029/2006wr005264. (22) Alperin, M. J.; Reeburgh, W. S.; Whiticar, M. J. Carbon and hydrogen isotope fractionation resulting from anaerobic methane oxidation. Global Biogeochem. Cycles 1988, 2, 279−288. (23) Lehmann, M. F.; Reichert, P.; Bernasconi, S. M.; Barbieri, A.; McKenzie, J. A. Modelling nitrogen and oxygen isotope fractionation during denitrification in a lacustrine redox-transition zone. Geochim. Cosmochim. Acta 2003, 67, 2529−2542. (24) Beekman, H. E.; Eggenkamp, H. G. M.; Appelo, C. A. J. An integrated modelling approach to reconstruct complex solute transport mechanismsCl and delta(37)Cl in pore water of sediments from a former brackish lagoon in The Netherlands. Appl. Geochem. 2011, 26, 257−268. (25) Groen, J.; Velstra, J.; Meesters, A. Salinization processes in paleowaters in coastal sediments of Suriname: evidence from delta Cl37 analysis and diffusion modelling. J. Hydrol. 2000, 234, 1−20. (26) Risk, D.; Kellman, L. Isotopic fractionation in non-equilibrium diffusive environments. Geophys. Res. Lett. 2008, 35, L02403 DOI: 10.1029/2007GL032374. (27) Bouchard, D.; Hohener, P.; Hunkeler, D. Carbon isotope fractionation during volatilization of petroleum hydrocarbons and diffusion across a porous medium: A column experiment. Environ. Sci. Technol. 2008, 42, 7801−7806. (28) Bouchard, D.; Hunkeler, D.; Gaganis, P.; Aravena, R.; Hohener, P.; Broholm, M. M.; Kjeldsen, P. Carbon isotope fractionation during diffusion and biodegradation of petroleum hydrocarbons in the unsaturated zone: Field experiment at Vaerlose airbase, Denmark, and modeling. Environ. Sci. Technol. 2008, 42, 596−601. (29) Hunkeler, D.; Aravena, R.; Shouakar-Stash, O.; Weisbrod, N.; Nasser, A.; Netzer, L.; Ronen, D. Carbon and chlorine isotope ratios of chlorinated ethenes migrating through a thick unsaturated zone of a sandy aquifer. Environ. Sci. Technol. 2011, 45, 8247−8253. (30) Chiogna, G.; Eberhardt, C.; Grathwohl, P.; Cirpka, O. A.; Rolle, M. Evidence of compound-dependent hydrodynamic and mechanical transverse dispersion by multitracer laboratory experiments. Environ. Sci. Technol. 2010, 44, 688−693. (31) Hunkeler, D.; Chollet, N.; Pittet, X.; Aravena, R.; Cherry, J. A.; Parker, B. L. Effect of source variability and transport processes on carbon isotope ratios of TCE and PCE in two sandy aquifers. J. Contam. Hydrol. 2004, 74, 265−282. (32) Cirpka, O. A.; Frind, E. O.; Helmig, R. Streamline-oriented grid generation for transport modelling in two-dimensional domains including wells. Adv. Water Resour. 1999, 22, 697−710. (33) Hunkeler, D.; Van Breukelen, B. M.; Elsner, M. Modeling chlorine isotope trends during sequential transformation of chlorinated ethenes. Environ. Sci. Technol. 2009, 43, 6750−6756. (34) van Breukelen, B. M.; Griffioen, J.; Roling, W. F. M.; van Verseveld, H. W. Reactive transport modelling of biogeochemical processes and carbon isotope geochemistry inside a landfill leachate plume. J. Contam. Hydrol. 2004, 70, 249−269. (35) van Breukelen, B. M.; Hunkeler, D.; Volkering, F. Quantification of sequential chlorinated ethene degradation by use of a reactive transport model incorporating isotope fractionation. Environ. Sci. Technol. 2005, 39, 4189−4197. (36) Domenico, P. A.; Palciauskas, V. V. Alternative boundaries in solid-waste management. Ground Water 1982, 20, 303−311. (37) Worch, E. A. A new equation for the calculation of diffusion coefficients for dissolved substances. Vom Wasser 1993, 81, 289−297. (38) Mills, R. Diffusion relationships in binary-system benzeneperdeuteriobenzene at 25°C. J. Phys. Chem. 1976, 80, 888−890. (39) Scheidegger, A. General theory of dispersion in porous media. J. Geophys. Res. 1961, 66, 3273−3278. 7707

dx.doi.org/10.1021/es301058z | Environ. Sci. Technol. 2012, 46, 7700−7708

Environmental Science & Technology

Article

(40) Chiogna, G.; Cirpka, O. A.; Grathwohl, P.; Rolle, M. Transverse mixing of conservative and reactive tracers in porous media: Quantification through the concepts of flux-related and critical dilution indices. Water Resour. Res. 2011, 47, W02505 DOI: 10.1029/ 2010wr009608. (41) Rolle, M.; Eberhardt, C.; Chiogna, G.; Cirpka, O. A.; Grathwohl, P. Enhancement of dilution and transverse reactive mixing in porous media: Experiments and model-based interpretation. J. Contam. Hydrol. 2009, 110, 130−142. (42) Chiogna, G.; Cirpka, O. A.; Grathwohl, P.; Rolle, M. Relevance of local compound-specific transverse dispersion for conservative and reactive mixing in heterogeneous porous media. Water Resour. Res. 2011, 47, W07540 DOI: 10.1029/2010wr010270. (43) Cirpka, O. A.; de Barros, F. P. J.; Chiogna, G.; Rolle, M.; Nowak, W. Stochastic flux-related analysis of transverse mixing in twodimensional heterogeneous porous media. Water Resour. Res. 2011, 47, W06515 DOI: 10.1029/2010wr010279. (44) Bauer, R. D.; Rolle, M.; Bauer, S.; Eberhardt, C.; Grathwohl, P.; Kolditz, O.; Meckenstock, R. U.; Griebler, C. Enhanced biodegradation by hydraulic heterogeneities in petroleum hydrocarbon plumes. J. Contam. Hydrol. 2009, 105, 56−68. (45) Abe, Y.; Aravena, R.; Zopfi, J.; Shouakar-Stash, O.; Cox, E.; Roberts, J. D.; Hunkeler, D. Carbon and chlorine isotope fractionation during aerobic oxidation and reductive dechlorination of vinyl chloride and cis-1,2-dichloroethene. Environ. Sci. Technol. 2009, 43, 101−107. (46) Hofstetter, T. B.; Reddy, C. M.; Heraty, L. J.; Berg, M.; Sturchio, N. C. Carbon and chlorine isotope effects during abiotic reductive dechlorination of polychlorinated ethanes. Environ. Sci. Technol. 2007, 41, 4662−4668. (47) Guilbeault, M. A.; Parker, B. L.; Cherry, J. A. Mass and flux distributions from DNAPL zones in sandy aquifers. Ground Water 2005, 43, 70−86. (48) Sudicky, E. A. A natural gradient experiment on solute transport in a sand aquiferSpatial variability of hydraulic conductivity and its role in the dispersion process. Water Resour. Res 1986, 22, 2069−2082. (49) Kampara, M.; Thullner, M.; Richnow, H. H.; Harms, H.; Wick, L. Y. Impact of bioavailability restrictions on microbially induced stable isotope fractionation. 2. Experimental evidence. Environ. Sci. Technol. 2008, 42, 6552−6558. (50) Thullner, M.; Kampara, M.; Richnow, H. H.; Harms, H.; Wick, L. Y. Impact of bioavailability restrictions on microbially induced stable isotope fractionation. 1. Theoretical calculation. Environ. Sci. Technol. 2008, 42, 6544−6551. (51) Knutson, C.; Valocchi, A.; Werth, C. Comparison of continuum and pore-scale models of nutrient biodegradation under transverse mixing conditions. Adv. Water Resour. 2007, 30, 1421−1431. (52) Thornton, S. F.; Bottrell, S. H.; Spence, K. H.; Pickup, R.; Spence, M. J.; Shah, N.; Mallinson, H. E. H.; Richnow, H. H. Assessment of MTBE biodegradation in contaminated groundwater using (13)C and (14)C analysis: Field and laboratory microcosm studies. Appl. Geochem. 2011, 26, 828−837. (53) Youngster, L. K. G.; Rosell, M.; Richnow, H. H.; Haggblom, M. M. Assessment of MTBE biodegradation pathways by two-dimensional isotope analysis in mixed bacterial consortia under different redox conditions. Appl. Microbiol. Biotechnol. 2010, 88, 309−317. (54) Rosell, M.; Barcelo, D.; Rohwerder, T.; Breuer, U.; Gehre, M.; Richnow, H. H. Variations in C-13/C-12 and D/H enrichment factors of aerobic bacterial fuel oxygenate degradation. Environ. Sci. Technol. 2007, 41, 2036−2043. (55) Liang, X.; Howlett, M. R.; Nelson, J. L.; Grant, G.; Dworatzek, S.; Lacrampe-Couloume, G.; Zinder, S. H.; Edwards, E. A.; Lollar, B. S. Pathway-dependent isotope fractionation during aerobic and anaerobic degradation of monochlorobenzene and 1,2,4-trichlorobenzene. Environ. Sci. Technol. 2011, 45, 8321−8327. (56) Vogt, C.; Cyrus, E.; Herklotz, I.; Schlosser, D.; Bahr, A.; Herrmann, S.; Richnow, H. H.; Fischer, A. Evaluation of toluene degradation pathways by two-dimensional stable isotope fractionation. Environ. Sci. Technol. 2008, 42, 7793−7800.

(57) Rolle, M.; Hochstetler, D.; Chiogna, G.; Kitanidis, P.; Grathwohl, P. Experimental investigation and pore-scale modeling interpretation of compound-specific transverse dispersion in porous media. Trans. Por. Med. 2012, DOI: 10.1007/s11242-012-9953-8.

7708

dx.doi.org/10.1021/es301058z | Environ. Sci. Technol. 2012, 46, 7700−7708