Utilizing distributed acoustic sensing and ocean bottom fiber optic cables for submarine structural characterization.

The sparsity of permanent seismic instrumentation in marine environments often limits the availability of subsea information on geohazards, including active fault systems, in both time and space. One sensing resource that provides observational access to the seafloor environment are existing networks of ocean bottom fiber optic cables; these cables, coupled to modern distributed acoustic sensing (DAS) systems, can provide dense arrays of broadband seismic observations capable of recording both seismic events and the ambient noise wavefield. Here, we report a marine DAS application which demonstrates the strength and limitation of this new technique on submarine structural characterization. Based on ambient noise DAS records on a 20 km section of a fiber optic cable offshore of Moss Landing, CA, in Monterey Bay, we extract Scholte waves from DAS ambient noise records using interferometry techniques and invert the resulting multimodal dispersion curves to recover a high resolution 2D shear-wave velocity image of the near seafloor sediments. We show for the first time that the migration of coherently scattered Scholte waves observed on DAS records can provide an approach for resolving sharp lateral contrasts in subsurface properties, particularly shallow faults and depositional features near the seafloor. Our results provide improved constraints on shallow submarine features in Monterey Bay, including fault zones and paleo-channel deposits, thus highlighting one of many possible geophysical uses of the marine cable network.

Improving Long-term Monitoring of Contaminated Groundwater at Sites where Attenuation-based Remedies are Deployed.

This study presents an effective approach to tackle the challenge of long-term monitoring of contaminated groundwater sites where remediation leaves residual contamination in the subsurface. Traditional long-term monitoring of contaminated groundwater sites focuses on measuring contaminant concentrations and is applicable to sites where contaminant mass is removed or degraded to a level below the regulatory standard. The traditional approach is less effective at sites where risk from metals or radionuclides continues to exist in the subsurface after remedial goals are achieved. We propose a long-term monitoring strategy for this type of waste site that focuses on measuring the hydrological and geochemical parameters that control attenuation or remobilization of contaminants while de-emphasizing contaminant-concentration measurements. We demonstrate how this approach would be more effective than traditional long-term monitoring, using a site in South Carolina, USA, where groundwater is contaminated by several radionuclides. A comprehensive enhanced attenuation remedy has been implemented at the site to minimize discharge of contamination to surface water. The immobilization of contaminants occurs in three locations by manipulation of hydrological and geochemical parameters, as well as by natural attenuation processes. Deployment of our proposed long-term monitoring strategy will combine subsurface and surface measurements using spectroscopic tools, geophysical tools, and sensors to monitor the parameters controlling contaminant attenuation. The advantage of this approach is that it will detect the potential for contaminant remobilization from engineered and natural attenuation zones, allowing potential adverse changes to be mitigated before contaminant attenuation is reversed.

Illuminating seafloor faults and ocean dynamics with dark fiber distributed acoustic sensing.

Distributed fiber-optic sensing technology coupled to existing subsea cables (dark fiber) allows observation of ocean and solid earth phenomena. We used an optical fiber from the cable supporting the Monterey Accelerated Research System during a 4-day maintenance period with a distributed acoustic sensing (DAS) instrument operating onshore, creating a ~10,000-component, 20-kilometer-long seismic array. Recordings of a minor earthquake wavefield identified multiple submarine fault zones. Ambient noise was dominated by shoaling ocean surface waves but also contained observations of in situ secondary microseism generation, post-low-tide bores, storm-induced sediment transport, infragravity waves, and breaking internal waves. DAS amplitudes in the microseism band tracked sea-state dynamics during a storm cycle in the northern Pacific. These observations highlight this method's potential for marine geophysics.

Distributed Acoustic Sensing Using Dark Fiber for Near-Surface Characterization and Broadband Seismic Event Detection.

We present one of the first case studies demonstrating the use of distributed acoustic sensing deployed on regional unlit fiber-optic telecommunication infrastructure (dark fiber) for broadband seismic monitoring of both near-surface soil properties and earthquake seismology. We recorded 7 months of passive seismic data on a 27 km section of dark fiber stretching from West Sacramento, CA to Woodland, CA, densely sampled at 2 m spacing. This dataset was processed to extract surface wave velocity information using ambient noise interferometry techniques; the resulting VS profiles were used to map both shallow structural profiles and groundwater depth, thus demonstrating that basin-scale variations in hydrological state could be resolved using this technique. The same array was utilized for detection of regional and teleseismic earthquakes and evaluated for long period response using records from the M8.1 Chiapas, Mexico 2017, Sep 8th event. The combination of these two sets of observations conclusively demonstrates that regionally extensive fiber-optic networks can effectively be utilized for a host of geoscience observation tasks at a combination of scale and resolution previously inaccessible.

Dynamics of geologic CO2 storage and plume motion revealed by seismic coda waves.

Quantifying the dynamics of sequestered CO2 plumes is critical for safe long-term storage, providing guidance on plume extent, and detecting stratigraphic seal failure. However, existing seismic monitoring methods based on wave reflection or transmission probe a limited rock volume and their sensitivity decreases as CO2 saturation increases, decreasing their utility in quantitative plume mass estimation. Here we show that seismic scattering coda waves, acquired during continuous borehole monitoring, are able to illuminate details of the CO2 plume during a 74-h CO2 injection experiment at the Frio-II well Dayton, TX. Our study reveals a continuous velocity reduction during the dynamic injection of CO2, a result that augments and dramatically improves upon prior analyses based on P-wave arrival times. We show that velocity reduction is nonlinearly correlated with the injected cumulative CO2 mass and attribute this correlation to the fact that coda waves repeatedly sample the heterogeneous distribution of cumulative CO2 in the reservoir zone. Lastly, because our approach does not depend on P-wave arrival times or require well-constrained wave reflections it can be used with many source-receiver geometries including those external to the reservoir, which reduces the risk introduced by in-reservoir monitoring wells. Our results provide an approach for quantitative CO2 monitoring and plume evolution that increases safety and long-term planning for CO2 injection and storage.

A new mini-triaxial cell for combined high-pressure and high-temperature in situ synchrotron X-ray microtomography experiments up to 400°C and 24†MPa.

A new experimental triaxial cell for in situ synchrotron X-ray micro-computed tomography aimed at imaging small samples of (6 mm × 19 mm) at high temperatures (up to 400°C) and pressures (up to 24 MPa confining) is presented. The system has flow-through capabilities, independent axial and radial pressure control, and has been developed and tested at the 8.3.2. beamline at the Advanced Light Source. The characteristics of this new experimental rig are described, along with the challenges, mainly concerning the combination of X-ray transparency with vessel strength at high temperature, and solutions found during the development stage. An experiment involving oil shale pyrolysis under subsurface conditions, highlighting the importance of a device able to operate in this pressure and temperature range, is also introduced. The availability of this cell enables an unprecedented range of experiments in the Earth Sciences, with a special focus on subsurface geothermal processes.

Attenuating Sulfidogenesis in a Soured Continuous Flow Column System With Perchlorate Treatment.

Hydrogen sulfide production by sulfate reducing bacteria (SRB) is the primary cause of oil reservoir souring. Amending environments with chlorate or perchlorate [collectively denoted (per)chlorate] represents an emerging technology to prevent the onset of souring. Recent studies with perchlorate reducing bacteria (PRB) monocultures demonstrated that they have the innate capability to enzymatically oxidize sulfide, thus PRB may offer an effective means of reversing souring. (Per)chlorate may be effective by (i) direct toxicity to SRB; (ii) competitive exclusion of SRB by PRB; or (iii) reversal of souring through re-oxidation of sulfide by PRB. To determine if (per)chlorate could sweeten a soured column system and assign a quantitative value to each of the mechanisms we treated columns flooded with San Francisco bay water with temporally decreasing amounts (50, 25, and 12.5 mM) of (per)chlorate. Geochemistry and the microbial community structure were monitored and a reactive transport model was developed, Results were compared to columns treated with nitrate or untreated. Souring was reversed by all treatments at 50 mM but nitrate-treated columns began to re-sour when treatment concentrations decreased (25 mM). Re-souring was only observed in (per)chlorate-treated columns when concentrations were decreased to 12.5 mM and the extent of re-souring was less than the control columns. Microbial community analyses indicated treatment-specific community shifts. Nitrate treatment resulted in a distinct community enriched in genera known to perform sulfur cycling metabolisms and genera capable of nitrate reduction. (Per)chlorate treatment enriched for (per)chlorate reducing bacteria. (Per)chlorate treatments only enriched for sulfate reducing organisms when treatment levels were decreased. A reactive transport model of perchlorate treatment was developed and a baseline case simulation demonstrated that the model provided a good fit to the effluent geochemical data. Subsequent simulations teased out the relative role that each of the three perchlorate inhibition mechanisms played during different phases of the experiment. These results indicate that perchlorate addition is an effective strategy for both souring prevention and souring reversal. It provides insight into which organisms are involved, and illuminates the interactive effects of the inhibition mechanisms, further highlighting the versatility of perchlorate as a sweetening agent.

Permafrost Degradation and Subsidence Observations during a Controlled Warming Experiment.

Global climate change has resulted in a warmer Arctic, with projections indicating accelerated modifications to permafrost in the near future. The thermal, hydrological, and mechanical physics of permafrost thaw have been hypothesized to couple in a complex fashion but data collection efforts to study these feedbacks in the field have been limited. As a result, laboratory and numerical models have largely outpaced field calibration datasets. We present the design, execution, and initial results from the first decameter-scale controlled thawing experiment, targeting coupled thermal/mechanical response, particularly the temporal sequence of surface subsidence relative to permafrost degradation at depth. The warming test was conducted in Fairbanks, AK, and utilized an array of in-ground heaters to induce thaw of a ~11 × 13 × 1.5 m soil volume over 63 days. The 4-D temperature evolution demonstrated that the depth to permafrost lowered 1 m during the experiment. The resulting thaw-induced surface deformation was ~10 cm as observed using a combination of measurement techniques. Surface deformation occurred over a smaller spatial domain than the full thawed volume, suggesting that gradients in cryotexture and ice content were significant. Our experiment provides the first large field calibration dataset for multiphysics thaw models.

The emerging role of 4D synchrotron X-ray micro-tomography for climate and fossil energy studies: five experiments showing the present capabilities at beamline 8.3.2 at the Advanced Light Source.

Continuous improvements at X-ray imaging beamlines at synchrotron light sources have made dynamic synchrotron X-ray micro-computed tomography (SXR-µCT) experiments more routinely available to users, with a rapid increase in demand given its tremendous potential in very diverse areas. In this work a survey of five different four-dimensional SXR-µCT experiments is presented, examining five different parameters linked to the evolution of the investigated system, and tackling problems in different areas in earth sciences. SXR-µCT is used to monitor the microstructural evolution of the investigated sample with the following variables: (i) high temperature, observing in situ oil shale pyrolysis; (ii) low temperature, replicating the generation of permafrost; (iii) high pressure, to study the invasion of supercritical CO2 in deep aquifers; (iv) uniaxial stress, to monitor the closure of a fracture filled with proppant, in shale; (v) reactive flow, to observe the evolution of the hydraulic properties in a porous rock subject to dissolution. For each of these examples, it is shown how dynamic SXR-µCT was able to provide new answers to questions related to climate and energy studies, highlighting the significant opportunities opened recently by the technique.

Clay, Water, and Salt: Controls on the Permeability of Fine-Grained Sedimentary Rocks.

The ability to predict the permeability of fine-grained soils, sediments, and sedimentary rocks is a fundamental challenge in the geosciences with potentially transformative implications in subsurface hydrology. In particular, fine-grained sedimentary rocks (shale, mudstone) constitute about two-thirds of the sedimentary rock mass and play important roles in three energy technologies: petroleum geology, geologic carbon sequestration, and radioactive waste management. The problem is a challenging one that requires understanding the properties of complex natural porous media on several length scales. One inherent length scale, referred to hereafter as the mesoscale, is associated with the assemblages of large grains of quartz, feldspar, and carbonates over distances of tens of micrometers. Its importance is highlighted by the existence of a threshold in the core scale mechanical properties and regional scale energy uses of shale formations at a clay content Xclay ≈ 1/3, as predicted by an ideal packing model where a fine-grained clay matrix fills the gaps between the larger grains. A second important length scale, referred to hereafter as the nanoscale, is associated with the aggregation and swelling of clay particles (in particular, smectite clay minerals) over distances of tens of nanometers. Mesoscale phenomena that influence permeability are primarily mechanical and include, for example, the ability of contacts between large grains to prevent the compaction of the clay matrix. Nanoscale phenomena that influence permeability tend to be chemomechanical in nature, because they involve strong impacts of aqueous chemistry on clay swelling. The second length scale remains much less well characterized than the first, because of the inherent challenges associated with the study of strongly coupled nanoscale phenomena. Advanced models of the nanoscale properties of fine-grained media rely predominantly on the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, a mean field theory of colloidal interactions that accurately predicts clay swelling in a narrow range of conditions (low salinity, low compaction, Na+ counterion). An important feature of clay swelling that is not predicted by these models is the coexistence, in most conditions of aqueous chemistry and dry bulk density, of two types of pores between parallel smectite particles: mesopores with a pore width of >3 nm that are controlled by long-range interactions (the osmotic swelling regime) and nanopores with a pore width <1 nm that are controlled by short-range interactions (the crystalline swelling regime). Nanogeochemical characterization and simulation techniques, including coarse-grained and all-atom molecular dynamics simulations, hold significant promise for the development of advanced constitutive relations that predict this coexistence and its dependence on aqueous chemistry.

Distributed Acoustic Sensing for Seismic Monitoring of The Near Surface: A Traffic-Noise Interferometry Case Study.

Ambient-noise-based seismic monitoring of the near surface often has limited spatiotemporal resolutions because dense seismic arrays are rarely sufficiently affordable for such applications. In recent years, however, distributed acoustic sensing (DAS) techniques have emerged to transform telecommunication fiber-optic cables into dense seismic arrays that are cost effective. With DAS enabling both high sensor counts ("large N") and long-term operations ("large T"), time-lapse imaging of shear-wave velocity (V S ) structures is now possible by combining ambient noise interferometry and multichannel analysis of surface waves (MASW). Here we report the first end-to-end study of time-lapse V S imaging that uses traffic noise continuously recorded on linear DAS arrays over a three-week period. Our results illustrate that for the top 20 meters the V S models that is well constrained by the data, we obtain time-lapse repeatability of about 2% in the model domain-a threshold that is low enough for observing subtle near-surface changes such as water content variations and permafrost alteration. This study demonstrates the efficacy of near-surface seismic monitoring using DAS-recorded ambient noise.

Alteration and Erosion of Rock Matrix Bordering a Carbonate-Rich Shale Fracture.

A novel reactive transport model has been developed to examine the processes that affect fracture evolution in a carbonate-rich shale. An in situ synchrotron X-ray microtomography experiment, flowing CO2 saturated water through a single fracture mini-core of Niobrara Shale provided the experimental observations for the development and testing of the model. The phenomena observed included the development of a porous altered layer, flow channeling, and increasingly limited calcite dissolution. The experimental observations cannot be explained by models that consider only mineral dissolution and development of an altered layer. The difference between the fracture volume change recorded by the microtomography images and what would be expected from mineral dissolution alone suggest that there is erosion of the altered layer as it develops. The numerical model includes this additional mechanism, with the erosion rate based on the thickness of the altered layer, and successfully captures the evolution of the geochemical reactions and morphology of the fracture. The findings imply that the abundance (with a threshold of approximately 35%) and reactivity of the rapidly reacting mineral control the development and erodibility of the altered layer on the fracture surfaces, and therefore fracture opening.

A 2.5D Reactive Transport Model for Fracture Alteration Simulation.

Understanding fracture alteration resulting from geochemical reactions is critical in predicting fluid migration in the subsurface and is relevant to multiple environmental challenges. Here, we present a novel 2.5D continuum reactive transport model that captures and predicts the spatial pattern of fracture aperture change and the development of an altered layer in the near-fracture region. The model considers permeability heterogeneity in the fracture plane and updates fracture apertures and flow fields based on local reactions. It tracks the reaction front of each mineral phase and calculates the thickness of the altered layer. Given this treatment, the model is able to account for the diffusion limitation on reaction rates associated with the altered layer. The model results are in good agreement with an experimental study in which a CO2-acidified brine was injected into a fracture in the Duperow Dolomite, causing dissolution of calcite and dolomite that result in the formation of a preferential flow channel and an altered layer. With an effective diffusion coefficient consistent with the experimentally observed porosity of the altered layer, the model captures the progressive decrease in the dissolution rate of the fast-reacting mineral in the altered layer.

Reactive Transport Model of Sulfur Cycling as Impacted by Perchlorate and Nitrate Treatments.

Microbial souring in oil reservoirs produces toxic, corrosive hydrogen sulfide through microbial sulfate reduction, often accompanying (sea)water flooding during secondary oil recovery. With data from column experiments as constraints, we developed the first reactive-transport model of a new candidate inhibitor, perchlorate, and compared it with the commonly used inhibitor, nitrate. Our model provided a good fit to the data, which suggest that perchlorate is more effective than nitrate on a per mole of inhibitor basis. Critically, we used our model to gain insight into the underlying competing mechanisms controlling the action of each inhibitor. This analysis suggested that competition by heterotrophic perchlorate reducers and direct inhibition by nitrite produced from heterotrophic nitrate reduction were the most important mechanisms for the perchlorate and nitrate treatments, respectively, in the modeled column experiments. This work demonstrates modeling to be a powerful tool for increasing and testing our understanding of reservoir-souring generation, prevention, and remediation processes, allowing us to incorporate insights derived from laboratory experiments into a framework that can potentially be used to assess risk and design optimal treatment schemes.

Microbial growth under supercritical CO2.

Growth of microorganisms in environments containing CO2 above its critical point is unexpected due to a combination of deleterious effects, including cytoplasmic acidification and membrane destabilization. Thus, supercritical CO2 (scCO2) is generally regarded as a sterilizing agent. We report isolation of bacteria from three sites targeted for geologic carbon dioxide sequestration (GCS) that are capable of growth in pressurized bioreactors containing scCO2. Analysis of 16S rRNA genes from scCO2 enrichment cultures revealed microbial assemblages of varied complexity, including representatives of the genus Bacillus. Propagation of enrichment cultures under scCO2 headspace led to isolation of six strains corresponding to Bacillus cereus, Bacillus subterraneus, Bacillus amyloliquefaciens, Bacillus safensis, and Bacillus megaterium. Isolates are spore-forming, facultative anaerobes and capable of germination and growth under an scCO2 headspace. In addition to these isolates, several Bacillus type strains grew under scCO2, suggesting that this may be a shared feature of spore-forming Bacillus spp. Our results provide direct evidence of microbial activity at the interface between scCO2 and an aqueous phase. Since microbial activity can influence the key mechanisms for permanent storage of sequestered CO2 (i.e., structural, residual, solubility, and mineral trapping), our work suggests that during GCS microorganisms may grow and catalyze biological reactions that influence the fate and transport of CO2 in the deep subsurface.

Isotopic insights into microbial sulfur cycling in oil reservoirs.

Microbial sulfate reduction in oil reservoirs (biosouring) is often associated with secondary oil production where seawater containing high sulfate concentrations (~28 mM) is injected into a reservoir to maintain pressure and displace oil. The sulfide generated from biosouring can cause corrosion of infrastructure, health exposure risks, and higher production costs. Isotope monitoring is a promising approach for understanding microbial sulfur cycling in reservoirs, enabling early detection of biosouring, and understanding the impact of souring. Microbial sulfate reduction is known to result in large shifts in the sulfur and oxygen isotope compositions of the residual sulfate, which can be distinguished from other processes that may be occurring in oil reservoirs, such as precipitation of sulfate and sulfide minerals. Key to the success of this method is using the appropriate isotopic fractionation factors for the conditions and processes being monitored. For a set of batch incubation experiments using a mixed microbial culture with crude oil as the electron donor, we measured a sulfur fractionation factor for sulfate reduction of -30‰. We have incorporated this result into a simplified 1D reservoir reactive transport model to highlight how isotopes can help discriminate between biotic and abiotic processes affecting sulfate and sulfide concentrations. Modeling results suggest that monitoring sulfate isotopes can provide an early indication of souring for reservoirs with reactive iron minerals that can remove the produced sulfide, especially when sulfate reduction occurs in the mixing zone between formation waters (FW) containing elevated concentrations of volatile fatty acids (VFAs) and injection water (IW) containing elevated sulfate. In addition, we examine the role of reservoir thermal, geochemical, hydrological, operational and microbiological conditions in determining microbial souring dynamics and hence the anticipated isotopic signatures.

Pore-scale controls on calcite dissolution rates from flow-through laboratory and numerical experiments.

A combination of experimental, imaging, and modeling techniques were applied to investigate the pore-scale transport and surface reaction controls on calcite dissolution under elevated pCO2 conditions. The laboratory experiment consisted of the injection of a solution at 4 bar pCO2 into a capillary tube packed with crushed calcite. A high resolution pore-scale numerical model was used to simulate the experiment based on a computational domain consisting of reactive calcite, pore space, and the capillary wall constructed from volumetric X-ray microtomography images. Simulated pore-scale effluent concentrations were higher than those measured by a factor of 1.8, with the largest component of the discrepancy related to uncertainties in the reaction rate model and its parameters. However, part of the discrepancy was apparently due to mass transport limitations to reactive surfaces, which were most pronounced near the inlet where larger diffusive boundary layers formed around grains and in slow-flowing pore spaces that exchanged mass by diffusion with fast flow paths. Although minor, the difference between pore- and continuum-scale results due to transport controls was discernible with the highly accurate methods employed and is expected to be more significant where heterogeneity is greater, as in natural subsurface materials.

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation.

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4+ production during urea hydrolysis were incorporated in the model and captured critical changes in the major metal species. The electrical phase increases were potentially due to ion exchange processes that modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to understand complex biogeochemical processes in the subsurface.

Hydrogeophysical methods for analyzing aquifer storage and recovery systems.

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site-specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ∼500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time-domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.