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 X_clay ≈ 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.

Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water-bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure-generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (dv/v) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km-long telecommunication fiber-optic cable in the Sacramento Valley, CA, a so-called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily dv/v beneath each location. Near the Sacramento river, dv/v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large dv/v variations. These observations reveal short-scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed dv/v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

The detailed mechanisms of the sealing of a single fracture, from hydration to almost complete closure by increase of confining pressure, as monitored from in situ synchrotron X-ray microtomography during the flow of carbonated water, is here shown for the first time. Different mechanisms play the key role at different stages in the evolution of the fracture. Hydration mechanically weakens the surfaces of the fracture and induces a first closure due to microcracking at the asperity contacts, increasing their size and creating choke points. Increase in confining stress promptly hydraulically seals the fracture by closing the main choke point, with a relative small deformation of the sample. Finally, the more pervasive mechanical deformation observed at higher stresses almost completely seals the whole fracture. The evolution of the sample has been described and quantified using 4D image processing, focusing on the evolution of aperture and digital volume correlation. Hydraulic properties of the sample at different stages have been modeled via Stokes flow simulation, and the results compared to the morphometric analysis, finding positive correlations with the average fracture aperture variation along the flowpath in function of confining pressure. Opalinus Clay is found to be a rock markedly prone to sealing in case of flow with carbonated water; this behavior is the result of the large fraction of clays and of its microstructure, lacking both cementing phases and large stiff particles. CO2 in this sample has no evident role in the evolution of the fracture; chemically induced weathering on the surface has not been detected, in contrast with the behavior observed in samples with carbonates as cementing phase.

Surface waves are advantageous for mapping seismic structures of permafrost, in which irregular velocity gradients are common and thus the effectiveness of refraction methods are limited. Nevertheless, the complex velocity structures that are common in permafrost environments often yield unusual dispersion spectra, in which higher-order and leaky modes are dominant. Such unusual dispersion spectra were prevalent in the multichannel surface-wave data acquired fromour permafrost study site at Barrow, Alaska. Owing to the difficulties in picking and identifying dispersion curves from these dispersion spectra, conventional surface-wave inversion methods become problematic to apply. To overcome these difficulties, we adopted a full-wavefield method to invert for velocitymodels that can best fit the dispersion spectra instead of the dispersion curves. The inferred velocity models were consistent with collocated electric resistivity results and with subsequent confirmation cores, which indicated the reliability of the recovered seismic structures. The results revealed embedded low-velocity zones underlying the ice-rich permafrost at our study site - an unexpected feature considering the low ground temperatures of -10°C to -8°C. The low velocities in these zones (~70%-80% lower than the overlying ice-rich permafrost) were most likely caused by saline pore-waters that prevent the ground from freezing, and the resultant velocity structures are vivid examples of complex subsurface properties in permafrost terrain. We determined that full-wavefield inversion of surface waves, although carrying higher computational costs than conventional methods, can be an effective tool for delineating the seismic structures of permafrost.

Near-surface seismic structure, particularly the shear wave velocity (Vs), can strongly affect local site response, and should be accurately estimated for ground motion prediction during seismic hazard assessment. The Imperial Valley (California), occupying the southern end of the Salton Trough, is a seismically active basin with thick surficial lacustrine sedimentary deposits. In this study, we utilize ambient noise records and local earthquake events for high-resolution near-surface characterization and site effect estimation with an unlit fiber-optic telecommunication infrastructure (dark fiber) in Imperial Valley by using the distributed acoustic sensing (DAS) technique. We apply ambient noise interferometry to retrieve coherent surface waves from DAS records, and evaluate performances of three different surface wave methods on DAS ambient noise dispersion imaging. We develop a quality control workflow to improve the dispersion measurement of noisy portions of the DAS data set by using a data selection strategy. Using the joint inversion of both the fundamental mode and higher overtones of Rayleigh waves, a high resolution two-dimensional (2D) Vs structure down to 70 m depth is obtained. We successfully achieve an improved Vs30 (the time-averaged shear-wave velocity in the top 30 m) model with higher spatial-resolution and reliability compared to the existing community model for the area. We also explore the potential for utilizing DAS earthquake events for site amplification estimation. The preliminary results reveal a clear anti-correlation between the approximated site response and the Vs30 profile. Our results indicate the potential utility of DAS deployed on dark fiber for near-surface characterization in appropriate contexts.

Measurements of nonlinear modulus and attenuation of fractures provide the opportunity to probe their mechanical state. We have adapted a low-frequency torsional apparatus to explore the seismic signature of fractures under low normal stress, simulating low effective stress environments such as shallow or high pore pressure reservoirs. We report strain-dependent modulus and attenuation for fractured samples of Duperow dolomite (a carbon sequestration target reservoir in Montana), Blue Canyon Dome rhyolite (a geothermal analog reservoir in New Mexico), and Montello granite (a deep basement disposal analog from Wisconsin). We use a simple single effective asperity partial slip model to fit our measured stress-strain curves and solve for the friction coefficient, contact radius, and full slip condition. These observations have the potential to develop into new field techniques for measuring differences in frictional properties during reservoir engineering manipulations and estimate the stress conditions where reservoir fractures and faults begin to fully slip.

A continuous active source seismic monitoring data set was collected with crosswell geometry during CO2 injection at the Frio-II brine pilot, near Liberty, TX. Previous studies have shown that spatiotemporal changes in the P wave first arrival time reveal the movement of the injected CO2 plume in the storage zone. To further constrain the CO2 saturation, particularly at higher saturation levels, we investigate spatial-temporal changes in the seismic attenuation of the first arrivals. The attenuation changes over the injection period are estimated by the amount of the centroid frequency shift computed by local time-frequency analysis. We observe that (1) at receivers above the injection zone seismic attenuation does not change in a physical trend; (2) at receivers in the injection zone attenuation sharply increases following injection and peaks at specific points varying with distributed receivers, which is consistent with observations from time delays of first arrivals; then, (3) attenuation decreases over the injection time. The attenuation change exhibits a bell-shaped pattern during CO2 injection. Under Frio-II field reservoir conditions, White's patchy saturation model can quantitatively explain both the P wave velocity and attenuation response observed. We have combined the velocity and attenuation change data in a crossplot format that is useful for model-data comparison and determining patch size. Our analysis suggests that spatial-temporal attenuation change is not only an indicator of the movement and saturation of CO2 plumes, even at large saturations, but also can quantitatively constrain CO2 plume saturation when used jointly with seismic velocity.

Distributed acoustic sensing (DAS) is a recently developed technique that has demonstrated its utility in the oil and gas industry. Here we demonstrate the potential of DAS in teleseismic studies using the Goldstone OpticaL Fiber Seismic experiment in Goldstone, California. By analyzing teleseismic waveforms from the 10 January 2018 M7.5 Honduras earthquake recorded on ~5,000 DAS channels and the nearby broadband station GSC, we first compute receiver functions for DAS channels using the vertical-component GSC velocity as an approximation for the incident source wavelet. The Moho P-to-s conversions are clearly visible on DAS receiver functions. We then derive meter-scale arrival time measurements along the entire 20-km-long array. We are also able to measure path-averaged Rayleigh wave group velocity and local Rayleigh wave phase velocity. The latter, however, has large uncertainties. Our study suggests that DAS will likely play an important role in many fields of passive seismology in the near future.

Biosouring is commonly encountered during secondary oil recovery when seawater or another high sulfate water source is utilized for flooding; as a result, effective souring control is of great interest to the oil industry. Here we describe a laboratory study to evaluate the relative effectiveness of souring interventions through the injection of nitrate, chlorate and perchlorate, collectively (per)chlorate, solutions and whether in-situ galvanic potential measurements can be used for convenient and quantitative tracking of sulfide dynamics. Nitrate has typically been the chemical of choice for souring treatments, while the efficacy of (per)chlorate as a new candidate inhibitor has only been explored recently. (Per)chlorate is known to inhibit oil reservoir souring via mechanisms such as toxicity, bio-competitive exclusion and sulfur redox cycling. Two sets of experiments under different matrix and inoculation conditions were conducted to evaluate treatment efficiency under variable baseline physical and biogeochemical conditions. Our data demonstrated the sensitivity of the galvanic potential signals to sulfide concentrations where the sulfide-galvanic potential correlation is similar to the theoretical predictions based on the Nernst equation, demonstrating the feasibility of using galvanic potential as a quick and economical method for quantifying in situ sulfide concentrations for tracking reservoir souring processes and subsequent intervention effectiveness. Our results show that all three chemicals were effective at suppressing sulfidogenesis. A reactive transport model for perchlorate treatment was developed to simulate the reaction processes and explore the interactions between the underlying competing mechanisms of this inhibitor. A baseline simulation captured the temporal patterns of the effluent chemical species. Subsequent simulations in which individual inhibition mechanisms were systematically removed elucidated the relative role that each inhibition mechanism played at the different phases of the experiment. The simulation results complement the experimental findings. Our study supports the potential advantages of souring control with (per)chlorate treatments, and the application of galvanic signal as an economic, in-situ monitoring approach for tracking souring dynamics and treatment efficacy.