Spectral Induced Polarization Signatures of Hydroxide Adsorption and

Mar 15, 2012 - Idaho National Laboratory, Idaho Falls, Idaho 83415, United States. §. Pacific Northwest National Laboratory, Richland, Washington 993...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/est

Spectral Induced Polarization Signatures of Hydroxide Adsorption and Mineral Precipitation in Porous Media Chi Zhang,*,† Lee Slater,† George Redden,‡ Yoshiko Fujita,‡ Timothy Johnson,§ and Don Fox‡ †

Rutgers University-Newark, Newark, New Jersey 07102, United States Idaho National Laboratory, Idaho Falls, Idaho 83415, United States § Pacific Northwest National Laboratory, Richland, Washington 99352, United States ‡

S Supporting Information *

ABSTRACT: The spectral induced polarization (SIP) technique is a promising approach for delineating subsurface physical and chemical property changes in a minimally invasive manner. To facilitate the understanding of position and chemical properties of reaction fronts that involve mineral precipitation in porous media, we investigated spatiotemporal variations in complex conductivity during evolution of urea hydrolysis and calcite precipitation reaction fronts within a silica gel column. The real and imaginary parts of complex conductivity were shown to be sensitive to changes in both solution chemistry and calcium carbonate precipitation. Distinct changes in imaginary conductivity coincided with increased hydroxide ion concentration during urea hydrolysis. In a separate experiment focused on the effect of hydroxide concentration on interfacial polarization of silica gel and well-sorted sand, we found a significant dependence of the polarization response on pH changes of the solution. We propose a conceptual model describing hydroxide ion adsorption behavior in silica gel and its control on interfacial polarizability. Our results demonstrate the utility of SIP for noninvasive monitoring of reaction fronts, and indicate its potential for quantifying geochemical processes that control the polarization responses of porous media at larger spatial scales in the natural environment.



INTRODUCTION Proposed approaches for in situ remediation of groundwater contamination include induced mineral precipitation schemes where contaminants are immobilized via coprecipitation or retarded by flowpath diversion. Successful application of such remediation strategies requires understanding of evolution of the relevant reaction fronts, and of changes in media properties due to precipitation. Given the transience of precipitation reaction fronts, the complexity of the reaction networks, and the challenge of accurately capturing potentially widely dispersed changes using traditional point sampling, novel monitoring approaches that sense both media property changes and temporal reaction fronts are needed to support both laboratory and field investigations. © 2012 American Chemical Society

Spectral induced polarization (SIP) is a promising geophysical technique for delineating variations in the properties of solid−fluid and fluid−fluid interfaces in porous media in a minimally invasively manner.1−4 Recently SIP has been applied for spatiotemporal monitoring of biogeochemical processes in both laboratory and field settings.2,3,5−7 Applications include monitoring of metal and nutrient cycling by sulfate-reducing microbes,2,6,8 iron corrosion,9,10 and calcite precipitation.11−13 A recent review4 describes how SIP has the potential to provide information on microbial growth, biofilm formation, and Received: Revised: Accepted: Published: 4357

December March 13, March 15, March 15,

9, 2011 2012 2012 2012

dx.doi.org/10.1021/es204404e | Environ. Sci. Technol. 2012, 46, 4357−4364

Environmental Science & Technology

Article

interface.21 It is usually assumed22,23 that the electrolytic and surface conduction paths add in parallel such that

biomineralization in porous media. The strong polarization of the metal/electrolyte interface results in large SIP signatures when metallic mineral phases (e.g., FeS) are present.2,3,6 However, the electrical properties of nonmetallic minerals and associated reaction fronts have been less studied. In one recent experiment,11 significant SIP signals were recorded during calcite precipitation on glass beads induced by mixing of CaCl2 and Na2CO3 solutions. These previous studies suggest that SIP may be a valuable technique for noninvasive monitoring of chemical reaction fronts associated with in situ mineral precipitation. We describe an experiment to evaluate the ability of SIP to monitor the spatiotemporal evolution of processes associated with calcite precipitation driven by urea hydrolysis within a silica gel column, catalyzed by immobilized extracellular urease enzyme. This model system is an abiotic analog to an approach that has been proposed for remediation of contaminants that can be coprecipitated in calcite, where ureolysis is mediated by indigenous subsurface microorganisms.14−16 The two primary reactions capable of changing interfacial electrical properties in this system are

σ′ = σele + σ′surf

and σ″ = σ″surf

urease

(1)

and (2)

To aid interpretation of the urease column experimental results, we conducted an additional study focused on the reactivity of hydroxide with silica gel, and specifically to examine whether (1) hydroxide adsorption affects SIP signals, and/or (2) changes in SIP signals can be used to discern between varying reaction fronts at solid−solution interfaces. Our study could improve interpretation of geophysical data collected for monitoring of in situ mineral precipitation. Spectral Induced Polarization. SIP measures the impedance magnitude and phase shift (φ) of a received sinusoidal voltage across a sample relative to the current waveform over a range of frequencies (typically 8.2 in the urease and posturease zones.

The column was operated in an upflow configuration. It was first flushed with 50 mM NaCl background solution (all influent solutions mentioned hence contain 50 mM NaCl) for two pore volumes (PV; 1 PV = 916 mL). This was followed by a bromide tracer test (22.5 mM LiBr; 20 mL/min). After removal of the bromide tracer, a 10 mM urea solution was introduced into the column at 10 mL/min for 544 min (Stage 1). Next a solution containing urea (10 mM), CaCl2 (10 mM), and SrCl2 (0.1 mM) was injected at 10 mL/min for 4400 min (Stage 2). The pH indicator phenol red (0.005 mM) was added to all influent solutions to facilitate visual observation of pH changes using an automated image capture system. The dye 4359

dx.doi.org/10.1021/es204404e | Environ. Sci. Technol. 2012, 46, 4357−4364

Environmental Science & Technology

Article

dynamic signal analyzer (National Instruments 4461) was employed to record SIP signatures over a frequency range from 0.1 Hz to 1 kHz at 100 logarithmic intervals on a daily basis. The measurement error is less than 1% for |σ| and ≈0.1 mrad for φ below 100 Hz. The current was injected through the coiled Ag-AgCl electrodes at both ends of the column and the resulting potential differences were recorded between Ag-AgCl electrode pairs spanning both sand zones (E1−E3, E3−E5, E13−E15, and E15−E17) and the silica gel zone (E6−E9 and E9−E12). The magnitude and φ were measured relative to a known high-precision resistor on channel 1. Spectral Induced Polarization Data Modeling. The SIP data sets from Experiment 2 were modeled with eq 8 using a two-step procedure. First, a Bayesian model using a Markovchain Monte Carlo (MCMC) based method was used to obtain the medians of unknown parameters by starting from an arbitrary set of initial values. This step provides marginal probability distributions but not explicitly optimal solutions.39 The medians of the parameter estimates were then used as the initial values for an iterative least-squares deterministic method with Marquardt regularization.40 However, there were uncertainties and nonuniqueness in the parameter estimation. For example, the correlation matrix reveals that m is highly correlated with τ, the average correlation coefficient between these two parameters being 0.67 for measurements in silica gel. This problem is frequently experienced, especially for small measurement signals, and is not easily overcome.

Figure 2. Experiment 1. Mineralized Ca and Sr in the column from destructive sampling of the solid media at the end of the experiment.

decreased later at the end of Stage 1 due to possible downstream migration of enzyme, which would have lowered the activity within the original urease zones. Urea arrival in Zone 5 during Stage 1 was evident in σ′ data as a consequence of the production of ionic species at the beginning of Stage 1. Afterward, the σ′ in Zone 5 increased significantly and reached a plateau at ≈0.52 S/m toward the end of Stage 1. With the onset of Stage 2, σf initially increased then dropped thereafter. The increase in σ′ at the beginning of Stage 2 was associated with the arrival of Ca2+ and Sr2+. Although changes in pore space and σf in the urease zones during Stage 2 were very small, apparent formation factor calculated from the ratio of σf to σ′ increased after the precipitation, consistent with a decrease in σ*surf as confirmed by σ″ measurements. This was presumably due to lower surface charge and lower surface area of CaCO 3 precipitation relative to silica gel.41 In posturease zones (using Zone 7 as representative), the σf changes were consistent with those observed in urease zones during Stage 1 except that the σf values were higher (≈ 0.64 S/ m) due to arrival of ions generated from upstream ureolysis. The correlation between σf and σ′ is strong (R2 = 0.96), and the σ′ first increased when Stage 1 was initiated and remained relatively stable initially, although it started gradually increasing during the latter part of Stage 1. This may have been due to migration/spreading of the urease zone as a consequence of incomplete enzyme immobilization. With the onset of Stage 2, the σ′ increase was much steeper before decreasing after 860 min. Qualitatively, the trends in σf, σ′, and pH were consistent with the migration of ionic species (background electrolyte and reaction products) coupled with ion depletion by CaCO3 precipitation. The pH and Imaginary Conductivity (σ″). The σ″ data compared with pH changes in Zones 1, 5, and 7 are plotted in Figure 3b, and the data in all zones are shown in Figure SI3 (the φ data are also plotted in the same manner in Figure SI4). The pH in Zone 1 was steady around 7.5, and σ″ was ≈0.00072 S/m and remained steady during the urea injection. With Stage 2, the pH dropped in the preurease zones due to pH changes in the influent, and the σ″ increased slightly in Zone 1.



RESULTS AND DISCUSSION Experiment 1: Urease Mediated Mineral Precipitation. Scanning electron microscopy (SEM) images of solids excavated from another column with similar experimental conditions after the conclusion of the experiment (Figure SI2) showed that crystals (CaCO3 composition inferred by energy dispersive spectroscopy) were produced at grain−grain contacts with a length scale of 100−200 μm. The ICP-MS analyses showed that mineralized Ca2+ reached 0.15 mmol/g in the urease zones, and that precipitation continued downstream, with decreased mineralized Ca2+ concentrations at the end of experiment (Figure 2). Fluid Conductivity (σf) and Real Conductivity (σ′). The hydrolysis of urea, with accompanying changes in solution chemistry and the subsequent formation of carbonate precipitates, was reflected in the SIP measurements. The σ′ and the σf data at different locations (Zones 1, 5, and 7) along the column are shown in Figure 3a, and the data for all zones (Zones 1−8) are plotted in Figure SI3. Note that for the individual zones, the measured pH and σf values at the two adjacent sampling ports were averaged to assist comparison with σ*(ω). In preurease zones (Zone 1 as representative), σf was steady at 0.531 S/m during Stage 1 and σ′ also remained steady (≈ 0.33 S/m). During Stage 2, the σf increased to 0.7 S/m due to the introduction of Ca2+ and Sr2+ into the column, and remained elevated through the experiment. Similarly, the σ′ also increased during Stage 2 to ≈0.42 S/m. The changes in σ′ and σf in the preurease zone were directly proportional throughout the experiment, consistent with a dominant control of electrolytic conduction (σele) on σ′, and it resulted in a formation factor of 1.64. In the urease zones ureolysis produces NH4+, HCO3−, and OH− (eq 1), which will increase both σf and σ′. In Zone 5 (as representative of the urease zone), σf increased to 0.62 S/m but 4360

dx.doi.org/10.1021/es204404e | Environ. Sci. Technol. 2012, 46, 4357−4364

Environmental Science & Technology

Article

Figure 3. Experiment 1. (a) Real (σ′) and fluid conductivities (σf) and (b) imaginary conductivities (σ″) and pH changes in the preurease (Zone 1 as representative), urease (Zone 5 as representative), and posturease zones (Zone 7 as representative) during ureolytically driven mineral precipitation experiment. The urea injection stage occurred from 0 to 544 min (Stage 1), followed by the coinjection of urea, Ca2+, and Sr2+ (Stage 2) began until the end of the experiment, and the two stages are separated by the vertical line. The σ*surf data are from 100 Hz.

Mineral precipitation can result in changes in surface area, formation porosity, and pore throat diameter, with significant impact on SIP signals.4 The electrical signals from precipitation of nonmetallic minerals such as calcite are small and associated with changes in pore volume/pore tortuosity and/or surface area/surface roughness.4 However, both SEM image and tracer test suggest porosity change during Stage 2 was insignificant (≈1%) as the precipitation of CaCO3 was localized at the grain−grain contacting points, indicating the contribution of pore space to decreased σ*surf is negligible. Instead, the CaCO3 precipitation has lower surface area relative to silica gel, and presumably reduced the total surface charge by forming a less negatively charged surface on the silica gel in the pH range 6− 10.41 In addition, calcium adsorption on silica gel is strong and would also lower the total surface charge. Our observations can be compared to those of Wu et al.,11 who conducted a column experiment with 3-mm glass beads and observed increased polarization responses during calcite precipitation. They suggested that the increase of σ″ during a discrete precipitation phase was caused by an increase in total surface area, due to the accumulation of small, well distributed, discrete calcite particles (average size 0.98). The pH increases due to ureolysis were observed first in Zone 3, then 4 and 5, and the σ″ also increased in Zones 3 to 5 in that order. The polarization magnitude in urease zones was much higher than that in the pre- and posturease zones. The σ″ in Zone 5 increased to a peak of 0.003 S/m close to the end of Stage 1. Early in Stage 2, the pH decreased due to proton production from calcium carbonate precipitation. The σ″ dropped in all urease zones during Stage 2, coincident with the carbonate mineral precipitation. The variations of phase data (Figure SI4) were similar to σ″ data, and also associated with pH changes in urease zones. In posturease zones, the pH was stable (≈ 7.5) during Stage 1, but increased during Stage 2. The pH in Zone 7 peaked at 950 min and then decreased by the end of the experiment. This can be explained by the migration of OH− formed during Stage 1 that continued after the start of stage 2. During Stage 1, the σ″ was much smaller (≈ 0.0004 S/m) in the posturease zones than in the urease zones. It increased slightly in Zone 7 near the end of Stage 1. After Stage 2 started, the σ″ in Zones 7 first increased, but then fell again, although it was still higher than the initial Stage 1 value at the end of the experiment. The σ″ changes coincided with changes in σ′ in these three posturease zones. Polarization Mechanisms. In this highly porous and high surface area medium, adsorption of ions on the surface of silica gel is significant. Ion or charge accumulation at the surface would presumably contribute to an increase in surface polarization, and therefore the σ″. The OH− adsorption, for example, would increase the negative surface charge, which would then attract counterions (i.e., H+, Na+, NH4+, Ca2+, Sr2+) into the Stern layer. During Stage 1, σ*surf has risen as indicated by increases in σ″. This could be attributed to an accumulation of surface charge resulting from adsorption of OH− or other ions. Similarly, both σ′ and σ″ gradually reduced without significant changes in σf in the downstream direction during Stage 2, suggesting that σ*surf decreased as CaCO3 precipitation occurred and pH decreased. These observations indicate that σ*surf (not σele) primarily controls the SIP responses to CaCO3 precipitation. 4361

dx.doi.org/10.1021/es204404e | Environ. Sci. Technol. 2012, 46, 4357−4364

Environmental Science & Technology

Article

increased, but then decreased on Days 7 and 8. There were no consistent variations in the characteristic critical frequency. The spectral phase (φ) changes and pH dependence in both silica gel and sand are shown in Figure SI5. Figure 5b shows the σ″ changes in silica gel and sand over time at the characteristic critical frequency (0.9 Hz), normalized chargeability mn (mn = m × σ0) in silica gel over time and compared with influent pH. The σ″ was stable at 0.0009 S/m for the first 2 days, increasing by about 10% when the influent pH increased from 7 to 8. When the influent pH changed to 10, the σ″ increased and reached a plateau by Day 4, while the pH of the effluent stabilized around 6.8 throughout the rest of the experiment. The behavior of σ″ suggests that OH− adsorbed onto silica gel after Day 4 did not contribute further to the polarization. The trend for σ″ (0.9 Hz) was very similar to the mn changes over time in silica gel, but the value of mn was 4−5 fold higher than σ″, since the σ″ here represents a single frequency whereas mn is the integrated measurement of charge storage over the investigated frequency range. The increased adsorption of OH− on the silica gel presumably caused the increase of mn, which was apparent from Day 3 when the pH of the influent increased to 10. The mn decreased on Day 6 and stayed relatively stable thereafter. The σ″ responses in the sand were much smaller than in the silica gel, and differences between the σ″ values over time were insignificant. We attribute the much weaker signals observed in the sand to the much smaller total surface area to pore volume ratio compared to silica gel. With similar particle size (300−212 μm) but no internal porous structure, the sands have surface area pzc) at salinity between 10−1 and 10−3 M. Skold et al.29 also showed increased σ″ as pH increased above the pzc in sand saturated by 0.01 M NaCl. Cole−Cole Parameters m and τ. Figure SI6a shows the change of the mn and τ over time in Experiment 2. The mn was well resolved from the model in both sand and silica gel. The mn for sand was significantly smaller than that for silica gel, and it tracked the change in σ″ as expected (both being measures of the polarization magnitude). The τ in the silica gel had a pattern similar to mn, slowly increasing to a maximum on Day 5, decreasing thereafter. It is reasonable to interpret changes in τ to changes in βs (eq 9). However, this would imply an increase in βs after Day 5, which is unreasonable given continued accumulation of charge at the surface. Furthermore, the σ″ stayed stable after Day 5, consistent with a maximum surface charge density on the surface sites. Therefore, the significance of changes in τ is uncertain, and possibly complicated by the nonuniqueness of this parameter. In the following section, we propose a conceptual model for OH− adsorption and polarization behavior in silica gel that could explain this plateau. Model for Hydroxide Adsorption on Silica Gel. Figure SI7 is a conceptual representation for the formation of surface charge on silica gel resulting from reactions with OH−. The silica gel matrix consists of silicon atoms joined by oxygen atoms in siloxane bonds. On contact with water, hydration of the silica surface initially produces multiple types of surface hydroxyl groups, and the surface of SiO2 has a net negative change above pH 2.5−3. Addition of OH− drives the deprotonation reaction forward and increases the concentration of “fixed” negative surface (often described as OH− adsorption). Formation of these fixed charge results in the formation of the EDL on the surface of primary backbones of silica gel, which can be modeled as a consisting of a Stern layer of mobile counterions that move tangentially along the surface and an outer diffuse layer. The arrangement of surface charges in the EDL with varying mobilities contributes to the SIP effect in the presence of an external electric field. We neglect the specific adsorption of Cl− anions at the surface since that usually occurs at pH < 343 where the net charge is close to zero, but where positively charged (protonated) sites can be present. The OH− ions are absorbed not only on the surface of the primary backbone, but possibly also in the numerous inner micropore structures. However, due to the small size of these pore structures and subsequent condensation reactions in silica gel, transport of Na+ (a counterion of OH−) may be restricted (Na+ hydrated radius 0.79 × 10−9 m) and will have less impact on the SIP signatures. We assume that the accumulation of OH− happened at a later stage (after Day 4) in our study, and this explains the plateau in σ″ since the excess negative charges are being stored in the internal structure. Because the chargeability mn relates to the polarizable interfacial charge storage, the behavior of mn (Figure 5b) is consistent with our conceptual model that the accumulation of OH− at the later stage did not contribute to the polarizability. Similarly to the observations in Experiment 1, we attribute the polarization responses observed in Experiment 2 to electrochemical polarization occurring at the EDL. The EDL polarization involves charge movement in both the fixed and diffused layers, and pH has an important influence on the polarization. The size of the mineral particle is fixed, so the

relaxation time is determined primarily by charge density and charge mobility in the EDL.17,35 We assert that the observed dependence of polarization on pH may reflect the relative contributions of charge density to the net surface conductivity response of the sample. Our conceptual model proposes that the increase of OH−onto the primary silica gel surface will increase accumulation of Na+, increasing the σ*surf. In reality, various surface conduction and competing physicochemical mechanisms likely combine to yield the net surface conductivity of the sample. To test the relative importance of the proposed OH− adsorption behavior on the silica gel as illustrated in the conceptual model, further studies should analyze the inner structure of the material and characterize ion adsorption behavior. We have shown that the SIP technique is sensitive enough to capture surface chemistry and associated interfacial properties of a solid−fluid surface. Our SIP measurements in the column experiments demonstrate the sensitivity of electrical geophysical methods to subtle and transient changes in porous media due to transport of ions and/or reactions such as precipitation and adsorption. The results suggest that pH controls the polarization response in silica gel. Although silica gel is a synthetic material, and the pH dependence of SIP is likely weaker for natural porous media, we believe our findings are still relevant to the assessment of methods for sensing and characterizing physical and chemical/biogeochemical events linked with subsurface remediation at larger spatial scales of natural environment. SIP can play an important role in improving the understanding of complex processes in porous media in both laboratory and field settings.



ASSOCIATED CONTENT

S Supporting Information *

Figures SI1−SI7. This information is available free of charge via the Internet at http://pubs.acs.org



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.

■ ■

ACKNOWLEDGMENTS This work was funded by the Department of Energy under contract DE-AC07-05ID14517. REFERENCES

(1) Atekwana, E. A.; Sauck, W. A.; Werkema, D. D. Investigations of geoelectrical signatures at a hydrocarbon contaminated site. J. Appl. Geophys. 2000, 44 (2−3), 167−180. (2) Personna, Y. R.; Ntarlagiannis, D.; Slater, L.; Yee, N.; O’Brien, M.; Hubbard, S. Spectral induced polarization and electrodic potential monitoring of microbially mediated iron sulfide transformations. J. Geophys. Res. 2008, 113 (G2), G02020. (3) William, K. H.; Kemna, A.; Wilkins, M. J.; Druhan, J.; Amtzen, E.; N’Guessan, A. L.; Hubbard, S. S.; Banfield, J. F. Geophysical monitoring of coupled microbial and geochemical processes during stimulated subsurface bioremediation. Environ. Sci. Technol. 2009, 43 (17), 6717−6723. (4) Atekwana, E. A.; Slater, L. D. Biogeophysics: A new frontier in Earth science research. Rev. Geophys. 2009, 47 (4), RG4004. (5) Williams, K. H.; Ntarlagiannis, D.; Slater, L. D.; Dohnalkova, A.; Hubbard, S. S.; Banfield, J. F. Geophysical Imaging of Stimulated

4363

dx.doi.org/10.1021/es204404e | Environ. Sci. Technol. 2012, 46, 4357−4364

Environmental Science & Technology

Article

Microbial Biomineralization. Environ. Sci. Technol. 2005, 39 (19), 7592−7600. (6) Ntarlagiannis, D.; Williams, K.; Slater, L.; Hubbard, S. Lowfrequency electrical response to microbial induced sulfide precipitation. J. Geophys. Res. 2005, 110 (G2), 12. (7) Davis, C. A.; Atekwana, E.; Atekwana, E.; Slater, L. D.; Rossbach, S.; Mormile, M. R. Microbial growth and biofilm formation in geologic media is detected with complex conductivity measurements. Geophys. Res. Lett. 2006, 33 (18), L18403. (8) Slater, L.; Ntarlagiannis, D.; Personna, Y. R.; Hubbard, S. Porescale spectral induced polarization signatures associated with FeS biomineral transformations. Geophys. Res. Lett. 2007, 34 (21), L21404. (9) Wu, Y.; D Slater, L.; Korte, N. Low Frequency Electrical Properties of Corroded Iron Barrier Cores. Environ. Sci. Technol. 2006, 40 (7), 2254−2261. (10) Wu, Y.; Slater, L. D.; Korte, N. Effect of Precipitation on Low Frequency Electrical Properties of Zerovalent Iron Columns. Environ. Sci. Technol. 2005, 39 (23), 9197−9204. (11) Wu, Y.; Hubbard, S.; Williams, K. H.; Ajo-Franklin, J. On the complex conductivity signatures of calcite precipitation. J. Geophys. Res. 2010, 115, G00G04. (12) Wu, Y.; Ajo-Franklin, J.; Spycher, N.; Hubbard, S.; Zhang, G.; Williams, K.; Taylor, J.; Fujita, Y.; Smith, R. Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation. Geochem. Trans. 2011, 12 (1), 7. (13) Wu, Y.; Versteeg, R.; Slater, L.; LaBrecque, D. Calcite precipitation dominates the electrical signatures of zero valent iron columns under simulated field conditions. J. Contam. Hydrol. 2009, 106 (3−4), 131−143. (14) Warren, L. A.; Maurice, P. A.; Parmar, N.; Ferris, F. G. Microbially Mediated Calcium Carbonate Precipitation: Implications for Interpreting Calcite Precipitation and for Solid-Phase Capture of Inorganic Contaminants. Geomicrobiol. J. 2001, 18 (1), 93−115. (15) Fujita, Y.; Ferris, F. G.; Lawson, R. D.; Colwell, F. S.; Smith, R. W. Calcium Carbonate Precipitation by Ureolytic Subsurface Bacteria. Geomicrobiol. J. 2000, 17 (4), 305−318. (16) Fujita, Y.; Taylor, J. L.; Gresham, T. L. T.; Delwiche, M. E.; Colwell, F. S.; McLing, T. L.; Petzke, L. M.; Smith, R. W. Stimulation Of Microbial Urea Hydrolysis In Groundwater To Enhance Calcite Precipitation. Environ. Sci. Technol. 2008, 42 (8), 3025−3032. (17) Lesmes, D. P.; Morgan, F. D. Dielectric spectroscopy of sedimentary rocks. J. Geophys. Res. 2001, 106 (B7), 13329−13346. (18) Slater, L.; Lesmes, D. P. Electrical-hydraulic relationships observed for unconsolidated sediments. Water Resour. Res. 2002, 38 (10), 1213. (19) Slater, L.; Glaser, D. Controls on induced polarization in sandy unconsolidated sediments and application to aquifer characterization. Geophysics 2003, 68 (5), 1547−1558. (20) Börner, F. D.; Schön, J. H. A relation between the quadrature component of electrical conductivity and the specific surface area of sedimentary rocks. Log Anal. 1991, 32, 612−613. (21) Lesmes, D. P.; Frye, K. M. Influence of pore fluid chemistry on the complex conductivity and induced polarization responses of Berea sandstone. J. Geophys. Res. 2001, 106 (B3), 4079−4090. (22) Waxman, M. H.; Smits, L. J. Electrical conductivities in oilbearing shaly sands. Trans. Am. Inst. Miner., Metall. Pet. Eng. 1968, 243, 107−122. (23) Vinegar, H. J.; Waxman, M. H. Induced polarization of shaly sands. Geophysics 1984, 49 (8), 1267−1287. (24) Archie, G. E. The electrical resistivity log as an aid in determining some reservoir characteristics. Trans. Am. Inst. Miner., Metall. Pet. Eng. 1942, 146, 54−62. (25) Scott, J. B. T.; Barker, R. D. Determining pore-throat size in Permo-Triassic sandstones from low-frequency electrical spectroscopy. Geophys. Res. Lett. 2003, 30 (9), 1450. (26) Binley, A.; Slater, L. D.; Fukes, M.; Cassiani, G. Relationship between spectral induced polarization and hydraulic properties of saturated and unsaturated sandstone. Water Resour. Res. 2005, 41 (12), W12417.

(27) Scott, J. B. T.; Barker, R. D. Characterization of sandstone by electrical spectroscopy for stratigraphical and hydrogeological investigations. Q. J. Eng. Geol. Hydrogeol. 2005, 38 (2), 143−154. (28) Vaudelet, P.; Revil, A.; Schmutz, M.; Franceschi, M.; Bégassat, P. Induced polarization signatures of cations exhibiting differential sorption behaviors in saturated sands. Water Resour. Res. 2011, 47 (2), W02526. (29) Skold, M.; Revil, A.; Vaudelet, P. The pH dependence of spectral induced polarization of silica sands: Experiment and modeling. Geophys. Res. Lett. 2011, 38 (12), L12304. (30) Cole, K. S.; Cole, R. H. Dispersion and Absorption in Dielectrics I. Alternating Current Characteristics. J. Chem. Phys. 1941, 9 (4), 341− 351. (31) Titov, K.; Komarov, V.; Tarasov, V.; Levitski, A. Theoretical and experimental study of time domain-induced polarization in watersaturated sands. J. Appl. Geophys. 2002, 50 (4), 417−433. (32) Pelton, W. H.; Ward, S. H.; Hallof, P. G.; Sill, W. R.; Nelson, P. H. Mineral discrimination and removal of inductive coupling with multifrequency IP. Geophysics 1978, 43 (3), 588−609. (33) Kruschwitz, S.; Binley, A.; Lesmes, D.; Elshenawy, A. Textural controls on low-frequency electrical spectra of porous media. Geophysics 2010, 75 (4), WA113−WA123. (34) Weller, A.; Breede, K.; Slater, L.; Nordsiek, S. Effect of changing water salinity on complex conductivity spectra of sandstones. Geophysics 2011, 76 (5), F315−F327. (35) Schwarz, G. A theory of the low-frequency dielectric dispersion of colloidal particles in electrolyte solution1,2. J. Phys. Chem. 1962, 66 (12), 2636−2642. (36) Benini, S.; Ciurli, S.; Rypniewski, W. R.; Wilson, K. S.; Mangani, S. Crystallization and preliminary high-resolution X-ray diffraction analysis of native and β-mercaptoethanol-inhibited urease from Bacillus pasteurii. Acta Crystallogr., D 1998, 54 (3), 409−412. (37) Wu, Y.; Slater, L.; Versteeg, R.; LaBrecque, D. A comparison of the low frequency electrical signatures of iron oxide versus calcite precipitation in granular zero valent iron columns. J. Contam. Hydrol. 2008, 95 (3−4), 154−167. (38) Versteeg, R. J.; Richardson, A. N.; Rowe, T. Web-Accessible Scientific Workflow System for Performance Monitoring. Environ. Sci. Technol. 2006, 40 (8), 2692−2698. (39) Chen, J.; Kemna, A.; Hubbard, S. S. A comparison between Gauss-Newton and Markov-chain Monte Carlo-based methods for inverting spectral induced-polarization data for Cole-Cole parameters. Geophysics 2009, 74 (4), Y7−Y7. (40) Kemna, A. Tomographic Inversion of Complex ResistivityTheory and Application; Der Andere Verlag: Osnabrück, Germany, 2000. (41) Eriksson, R.; Merta, J.; Rosenholm, J. B. The calcite/water interface: I. Surface charge in indifferent electrolyte media and the influence of low-molecular-weight polyelectrolyte. J. Colloid Interface Sci. 2007, 313 (1), 184−193. (42) Bolt, G. H. Determination of the Charge Density of Silica Sols. J. Phys. Chem. 1957, 61 (9), 1166−1169. (43) Leroy, P.; Revil, A.; Kemna, A.; Cosenza, P.; Ghorbani, A. Complex conductivity of water-saturated packs of glass beads. J. Colloid Interface Sci. 2008, 321 (1), 103−117. (44) Revil, A. Spectral induced polarization of shaly sands: Influence of the electrical double layer. Water Resour. Res. 2012, 48 (2), W02517.

4364

dx.doi.org/10.1021/es204404e | Environ. Sci. Technol. 2012, 46, 4357−4364