Resistivity and self-potential tomography applied to groundwater remediation and contaminant plumes: Sandbox and field experiments

Geophysical methods can be used to remotely characterize contaminated sites and monitor in situ enhanced remediation processes. We have conducted one sandbox experiment and one contaminated field investigation to show the robustness of electrical resistivity tomography and self-potential (SP) tomography for these applications. In the sandbox experiment, we injected permanganate in a trichloroethylene (TCE)-contaminated environment under a constant hydraulic gradient. Inverted resistivity tomograms are able to track the evolution of the permanganate plume in agreement with visual observations made on the side of the tank. Self-potential measurements were also performed at the surface of the sandbox using non-polarizing Ag-AgCl electrodes. These data were inverted to obtain the source density distribution with and without the resistivity information. A compact horizontal dipole source located at the front of the plume was obtained from the inversion of these self-potential data. This current dipole may be related to the redox reaction occurring between TCE and permanganate and the strong concentration gradient at the front of the plume. We demonstrate that time-lapse self-potential signals can be used to track the kinetics of an advecting oxidizer plume with acceptable accuracy and, if needed, in real time, but are unable to completely resolve the shape of the plume. In the field investigation, a 3D resistivity tomography is used to characterize an organic contaminant plume (resistive domain) and an overlying zone of solid waste materials (conductive domain). After removing the influence of the streaming potential, the identified source current density had a magnitude of 0.5 A m-2. The strong source current density may be attributed to charge movement between the neighboring zones that encourage abiotic and microbially enhanced reduction and oxidation reactions. In both cases, the self-potential source current density is located in the area of strong resistivity gradient

Quantifying precision and accuracy of measurements of dissolved inorganic carbon stable isotopic composition using continuous-flow isotope-ratio mass spectrometry

RATIONALE: We describe an analytical procedure that allows sample collection and measurement of carbon isotopic composition (δ13CV-PDB value) and dissolved inorganic carbon concentration, [DIC], in aqueous samples without further manipulation post field collection. By comparing outputs from two different mass spectrometers, we quantify with the statistical rigour uncertainty associated with the estimation of an unknown measurement. This is rarely undertaken, but it is needed to understand the significance of field data and to interpret quality assurance exercises.<p></p> METHODS: Immediate acidification of field samples during collection in evacuated, pre-acidified vials removed the need for toxic chemicals to inhibit continued bacterial activity that might compromise isotopic and concentration measurements. Aqueous standards mimicked the sample matrix and avoided headspace fractionation corrections. Samples were analysed using continuous-flow isotope-ratio mass spectrometry, but for low DIC concentration the mass spectrometer response could be non-linear. This had to be corrected for.<p></p> RESULTS: Mass spectrometer non-linearity exists. Rather than estimating precision as the repeat analysis of an internal standard, we have adopted inverse linear calibrations to quantify the precision and 95% confidence intervals (CI) of the δ13CDIC values. The response for [DIC] estimation was always linear. For 0.05–0.5 mM DIC internal standards, however, changes in mass spectrometer linearity resulted in estimations of the precision in the δ13CVPDB value of an unknown ranging from ± 0.44‰ to ± 1.33‰ (mean values) and a mean 95% CI half-width of ±1.1–3.1‰.<p></p> CONCLUSIONS: Mass spectrometer non-linearity should be considered in estimating uncertainty in measurement. Similarly, statistically robust estimates of precision and accuracy should also be adopted. Such estimations do not inhibit research advances: our consideration of small-scale spatial variability at two points on a small order river system demonstrates field data ranges larger than the precision and uncertainties. However, without such statistical quantification, exercises such as inter-lab calibrations are less meaningful.<p></p&gt

Structure of the lithosphere beneath the Barotse Basin, western Zambia, from magnetotelluric data.

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Tectonics, 38(2), (2019):666-686. doi:10.1029/2018TC005246.A magnetotelluric survey in the Barotse Basin of western Zambia shows clear evidence for thinned lithosphere beneath an orogenic belt. The uppermost asthenosphere, at a depth of 60–70 km, is highly conductive, suggestive of the presence of a small amount of partial melt, despite the fact that there is no surface expression of volcanism in the region. Although the data support the presence of thicker cratonic lithosphere to the southeast of the basin, the lithospheric thickness is not well resolved and models show variations ranging from ~80 to 150 km in this region. Similarly variable is the conductivity of the mantle beneath the basin and immediately beneath the cratonic lithosphere to the southeast, although the conductivity is required to be elevated compared to normal lithospheric mantle. In a general sense, two classes of model are compatible with the magnetotelluric data: one with a moderately conductive mantle and one with more elevated conductivities. This latter class would be consistent with the impingement of a stringer of plume‐fed melt beneath the cratonic lithosphere, with the melt migrating upslope to thermally erode lithosphere beneath the orogenic belt that is overlain by the Barotse Basin. Such processes are potentially important for intraplate volcanism and also for development or propagation of rifting as lithosphere is thinned and weakened by melt. Both models show clear evidence for thinning of the lithosphere beneath the orogenic belt, consistent with elevated heat flow data in the region.Funding for MT acquisition and analysis was provided by the National Science Foundation grant EAR‐1010432 through the Continental Dynamics Program. The data used in this study are available for download at the IRIS Data Management Center through the DOI links cited in Jones et al. (2003–2008; https://doi.org/10.17611/DP/EMTF/SAMTEX) and Evans et al. (2012; https://doi.org/10.17611/DP/EMTF/PRIDE/ZAM). We would like to thank the field crew from the Geological Survey Department, Zambia, for their assistance in collecting data. Matthew Chamberlain, David Margolius, and Colin Skinner, formerly of Northeastern University, are also thanked for their field assistance. Data are available from the corresponding author pending their submission to the IRIS DMC repository at which point they will be publically available. This is Oklahoma State University, Boone Pickens School of Geology contribution number 2019‐99.2019-07-3

Field Evidence for Geophysical Detection of Subsurface Zones of Enhanced Microbial Activity

Geochemical data from closely spaced vertical intervals in a hydrocarbon-impacted aquifer were used to assess the relationship between bulk conductivity and zones of enhanced microbial activity. The bulk conductivity was measured using in situ vertical resistivity probes. Microbial activity was verified using terminal electron acceptors (nitrate, sulfate, iron, and manganese), dissolved inorganic carbon (DIC), and major ion chemistry. Peaks in bulk conductivity in the aquifer overlapped with zones where nitrates and sulfates were depleted, total petroleum hydrocarbon, iron, manganese, dissolved ions, and DIC were elevated, suggesting a link between higher electrical conductivity and zones of enhanced microbial activity stimulated by the presence of hydrocarbon. Thus the subsurface expression of microbial activity is apparently recorded in the bulk conductivity measurements. Our results argue for combining geophysics with biogeochemistry studies to delineate subsurface zones of enhanced microbial activity

Temporal Geophysical Signatures from Contaminant-Mass Remediation

We have previously documented changes in bulk electrical conductivity, self-potential (SP), and ground-penetrating-radar (GPR) reflections in a field setting caused by biogeochemical transformations of hydrocarbon-contaminated media. These transformations are associated with hydrocarbon biodegradation. The results of surface geophysical surveys acquired in 1996, 2003, and 2007 document changes in geophysical signatures associated with removing hydrocarbon mass in the contaminated zone. Initial investigations in 1996 showed that relative to background, the contaminated area was characterized by higher bulk electrical conductivity, positive SP anomaly, and attenuated GPR reflections. Repeated surveys in 2003 and 2007 over the contaminated area showed that in 2007, the bulk electrical conductivity had reverted to near-background conditions, the positive SP anomaly became more negative, and the zone of attenuated GPR reflections showed increased signal strength. Removal of hydrocarbon mass in the vadose zone over the plume by a soil vapor-extraction system installed in 2001 was primarily responsible for the changing geophysical responses. Although chemical data from groundwater showed a 3-m-thick conductive plume in 2007, the plume was not imaged by electrical resistivity. Forward modeling suggests that apparent bulk electrical conductivity of the saturated zone plume has to be three to five times higher than background values to be imaged by electrical resistivity. We suggest that removing hydrocarbon-contaminant-mass reduction by natural or engineered bioremediation can be imaged effectively by temporal geophysical surveys

Effects of Microbial Processes on Electrolytic and Interfacial Electrical Properties of Unconsolidated Sediments

The effect of microbial processes on electrical properties of unconsolidated sediments was investigated in a laboratory experiment consisting of biotic and abiotic sand columns. The biotic column (nutrient, diesel and bacteria) showed (a) temporal increase in the real, imaginary, and surface conductivity, and (b) temporal decrease in the formation factor. The abiotic columns (nutrient; and nutrient and diesel) showed no significant changes. Increase in microbial population numbers, decrease in organic carbon source, nitrate, and sulfate and increase in dissolved inorganic carbon and fluid conductivity were indicative of microbial activity in the biotic column. We also measure relative increase in the interfacial electrical properties that exceed relative increase in the electrolytic conductivity. Thus changes in the real and imaginary conductivity were induced by microbial processes. These results suggest that interpretation of geoelectrical data from near surface environments should consider effects of microbial processes

Microbial Growth and Biofilm Formation in Geologic Media Is Detected with Complex Conductivity Measurements

Complex conductivity measurements (0.1-1000 Hz) were obtained from biostimulated sand-packed columns to investigate the effect of microbial growth and biofilm formation on the electrical properties of porous media. Microbial growth was verified by direct microbial counts, pH measurements, and environmental scanning electron microscope imaging. Peaks in imaginary (interfacial) conductivity in the biostimulated columns were coincident with peaks in the microbial cell concentrations extracted from sands. However, the real conductivity component showed no discernible relationship to microbial cell concentration. We suggest that the observed dynamic changes in the imaginary conductivity (σ″) arise from the growth and attachment of microbial cells and biofilms to sand surfaces. We conclude that complex conductivity techniques, specifically imaginary conductivity measurements are a proxy indicator for microbial growth and biofilm formation in porous media. Our results have implications for microbial enhanced oil recovery, CO2 sequestration, bioremediation, and astrobiology studies

Thermal Perturbations beneath the Incipient Okavango Rift Zone, Northwest Botswana

We used aeromagnetic and gravity data to investigate the thermal structure beneath the incipient Okavango Rift Zone (ORZ) in northwestern Botswana in order to understand its role in strain localization during rift initiation. We used three-dimensional (3-D) inversion of aeromagnetic data to estimate the Curie Point Depth (CPD) and heat flow under the rift and surrounding basement. We also used two-dimensional (2-D) power-density spectrum analysis of gravity data to estimate the Moho depth. Our results reveal shallow CPD values (8-15 km) and high heat flow (60-90 mW m-2) beneath a ∼60 km wide NE-trending zone coincident with major rift-related border faults and the boundary between Proterozoic orogenic belts. This is accompanied by thin crust ( \u3c 30 km) in the northeastern and southwestern parts of the ORZ. Within the Precambrian basement areas, the CPD values are deeper (16-30 km) and the heat flow estimates are lower (30-50 mW m-2), corresponding to thicker crust (∼40-50 km). We interpret the thermal structure under the ORZ as due to upward migration of hot mantle fluids through the lithospheric column that utilized the presence of Precambrian lithospheric shear zones as conduits. These fluids weaken the crust, enhancing rift nucleation. Our interpretation is supported by 2-D forward modeling of gravity data suggesting the presence of a wedge of altered lithospheric mantle centered beneath the ORZ. If our interpretation is correct, it may result in a potential paradigm shift in which strain localization at continental rift initiation could be achieved through fluid-assisted lithospheric weakening without asthenospheric involvement

Thermal Perturbations beneath the Incipient Okavango Rift Zone, Northwest Botswana

We used aeromagnetic and gravity data to investigate the thermal structure beneath the incipient Okavango Rift Zone (ORZ) in northwestern Botswana in order to understand its role in strain localization during rift initiation. We used three-dimensional (3-D) inversion of aeromagnetic data to estimate the Curie Point Depth (CPD) and heat flow under the rift and surrounding basement. We also used two-dimensional (2-D) power-density spectrum analysis of gravity data to estimate the Moho depth. Our results reveal shallow CPD values (8-15 km) and high heat flow (60-90 mW m-2) beneath a ∼60 km wide NE-trending zone coincident with major rift-related border faults and the boundary between Proterozoic orogenic belts. This is accompanied by thin crust ( \u3c 30 km) in the northeastern and southwestern parts of the ORZ. Within the Precambrian basement areas, the CPD values are deeper (16-30 km) and the heat flow estimates are lower (30-50 mW m-2), corresponding to thicker crust (∼40-50 km). We interpret the thermal structure under the ORZ as due to upward migration of hot mantle fluids through the lithospheric column that utilized the presence of Precambrian lithospheric shear zones as conduits. These fluids weaken the crust, enhancing rift nucleation. Our interpretation is supported by 2-D forward modeling of gravity data suggesting the presence of a wedge of altered lithospheric mantle centered beneath the ORZ. If our interpretation is correct, it may result in a potential paradigm shift in which strain localization at continental rift initiation could be achieved through fluid-assisted lithospheric weakening without asthenospheric involvement

Solid State Electron Transfer via Bacterial Nanowires: Contributions Toward a Mechanistic Understanding of Geophysical Response of Biostimulated Subsurface

The degradation of organic matter by microorganisms provides a source of electrical potential or so-called 'self potential' (SP) that can be measured by using a voltmeter. During this process electrons are being produced as a waste-product and bacterial cells have to dispose of these to allow for the complete biodegradation of organic matter. Especially in anaerobic microbial communities, exo-cellular electron transfer is the most important driving force behind this process and organisms have developed different, but also similar, ways to transfer electrons to other microorganisms. Recently, it has been postulated that direct electron transfer from cell-to-cell is actually done by 'hard-wired' microorganisms. This shuttling of electrons is most likely done by certain c-type cytochromes that form the functional part of electrically conductive nanowires. In this study we investigated if nanowires can explain the geoelectrical (self potential and spectral induced polarization) signals observed at some biostimulated environments such as DOE sites. The objectives of our project are to: (1) investigate any temporal changes in the geophysical signatures (Self Potential (SP) and Induced Polarization (IP)) associated with nanowires of the bacterium Shewanella oneidensis MR-1, wild type and mtrc/omcA deletion mutant, (2) demonstrate that mutant strains of bacteria that produce nonconductive nanowires do not contribute to geoelectrical responses. We accomplished the following: (1) Provided training to students and a postdoctoral fellow that worked on the project, (2) Conducted several SP & IP measurements correlating the distribution of nanowires and SIP/SP signals in partial fulfillment of object No. 1 and 2. On the following we will report and discuss the results of our last experiment with some emphasis on the source mechanisms of both SP and IP associated with Shewanella oneidensis MR-1, wild type in sand columns