CO₂ injection into porous sandstone reservoirs offers a promising pathway to curb anthropogenic carbon emissions, but poses risks of leakage and induced seismicity from stress perturbations and fault reactivation without meticulous monitoring. Here, we present a time-lapse monitoring approach based on laboratory measurements of ultrasonic Vp, Vs and corresponding peak amplitudes in critically stressed, partially saturated North Sea sandstones (porosity 9–23%). Our experiments show that Vp and Vs exhibit higher sensitivity (4–15%) to stress changes compared to fluid saturation changes (0.8–1%), whereas amplitudes are more responsive (30–500%) to saturation, showing staggered change when brine is displaced by CO₂. Under pure stress perturbation, amplitude variations are smaller (10–50%). During elastic deformation, the Vp/Vs ratio decreases while the ratio of their corresponding amplitudes increases, underscoring the need for both P- and S-wave measurements. Velocity and amplitude changes are more pronounced in high-porosity rocks. In a critically stressed state (beyond yield/before failure), the rise in pore fluid density from CO₂ injection boosts shear wave amplitudes, offsetting attenuation from inelastic deformation. Knowing the pre-injection stress state enables these velocity and amplitude trends to serve as robust indicators of reservoir conditions during and after CO₂ injection. This cost-effective approach can be adapted to reservoir-scale monitoring and extends beyond CCS, supporting enhanced detection of stress and fluid-induced changes in subsurface formations.

The global expansion of subsurface CO₂ and hydrogen storage, alongside geothermal energy development, offers promising pathways for gigaton-scale CO₂ abatement. However, fluid injections and associated thermal effects can significantly alter reservoir stress states, risking fault reactivation and compromising caprock integrity. Direct stress measurements in the subsurface remain technically challenging, particularly beyond the near-wellbore zone. This study investigates how stress-induced changes in ultrasonic P- and S-wave velocities and amplitudes can serve as early indicators of irreversible rock deformation. Using triaxial cyclic and failure experiments on core samples from offshore Netherlands (depths: 3.1–4.2 km; porosity: 8–23 %), we demonstrate that wave velocities and amplitudes increase with axial loading in the elastic regime but decline progressively following crack initiation—well before mechanical failure. This trend reversal provides a reliable sonic precursor to failure. We propose a field-applicable traffic-light monitoring framework using sonic parameters to infer stress changes during injection operations. The observed inverse relationships between porosity and both mechanical strength and sonic velocity, along with the porosity-dependent velocity enhancement under confinement, present a novel opportunity to develop constitutive geomechanical models directly from reservoir sonic logs. This work advances non-invasive stress monitoring approaches and provides engineering geologists with robust tools to improve safety and predictability in subsurface energy storage projects. Moreover, such techniques can also be translated to integrity monitoring for underground mines and engineered structures.

Forecasting the occurrence of natural hazards, such as earthquakes or landslides, remain very challenging. These hazards are often caused by stress changes in the subsurface, therefore detecting and monitoring these changes can help the prediction and mitigation. Active ultrasonic transmission experiments were performed on Red Pfaelzer sandstones to investigate the monitoring and forecasting potential of these measurements. The sandstone samples were loaded until failure at different initial confining stress conditions. The forecasting potential to failure of different analysis methods, such as coda wave interferometry or wave attenuation, is investigated and compared. Our results show we can detect the forecast the upcoming failure of the samples from 40 to 70% of its failure point. Small differences between each analysis method are visible, but the trend of the signal is leading and therefore a robust prediction of failure can be made by combining analysis methods. In this paper, we propose a traffic light forecasting system using the precursory signals from ultrasonic monitoring. This system is applicable for monitoring failure at various depths and or stress conditions, for a better prediction of small stress-induced changes in the subsurface and thus mitigation of failure (natural hazards) in the subsurface.

A central goal of laboratory seismology is to infer large-scale seismic processes from small-scale experiments, with acoustic emissions (AE) being a common observable. These signals, indicative of microfracturing, slip localization, and damage evolution, are often paralleled with earthquakes to understand seismic behaviors. This study challenges traditional perspectives by applying Coulomb rate-and-state seismicity theory, originally developed for earthquake clustering, to AE experiments. This theory maps stressing history to seismicity rates using rate-and-state friction, however, its validity under controlled experimental conditions remains an open question. We conducted four experiments on a sawcut sample of red felser sandstone, representing a fault under variable stress conditions. Adjustments in loading rates and initial conditions revealed that, while a single free parameter A—related to the direct effect—should suffice, a rescaling of the model by 1.5 to 2.2 was necessary for fitting the data. Differences in values across experiments appeared mostly non-systematic, and partial data usage did not yield consistently systematic parameter migrations. These findings suggest that fault microstructure may complexly alter parameter values during loading beyond what is accounted for in the Coulomb rate-and-state theory. Nonetheless, with the introduction of the scaling parameter, the Coulomb rate-and-state theory effectively captures the fundamental aspects of AE responses to complex controlled loading histories.

Induced earthquakes are still highly unpredictable, and often caused by variations in pore fluid pressure. Monitoring and understanding the mechanisms of fluid-induced fault slip is essential for seismic risk mitigation and seismicity forecasting. Fluid-induced slip experiments were performed on critically stressed faulted sandstone samples, and the evolution of the actively sent ultrasonic waves throughout the experiment was measured. Two different fault types were used: smooth saw-cut fault samples at a 35° angle, and a rough fault created by in situ faulting of the samples. Variations in the seismic slip velocity and friction along the fault plane were identified by the coda of the ultrasonic waves. Additionally, ultrasonic amplitudes show precursory signals to laboratory fault reactivation. Our results show that small and local variations in stress before fault failure can be inferred using coda wave interferometry for time-lapse monitoring, as coda waves are more sensitive to small perturbations in a medium than direct waves. Hence, these signals can be used as precursors to laboratory fault slip and to give insight into reactivation mechanisms. Our results show that time-lapse monitoring of coda waves can be used to monitor local stress changes associated with fault reactivation in this laboratory setting of fluid-induced fault reactivation. This is a critical first step toward a method for continuous monitoring of natural fault zones, contributing to seismic risk mitigation of induced and natural earthquakes.

As part of the Synergetic Utilisation of CO (Formula presented.) storage Coupled with geothermal EnErgy Deployment project, investigating CO (Formula presented.) reinjection with different seismic methods, both passive and active seismic surveys have been conducted at the geothermal power plant at Hellisheiði, Iceland. During the 2021 survey, two geophone lines recorded noise for a week. We process the passive-source data with seismic interferometry to image the subsurface structure around the CarbFix2 reinjection reservoir. To improve image quality, we perform an illumination analysis to select only noise panels dominated by body-wave energy. The results show that most noise panels are dominated by air-wave energy arriving from the direction of the power plant. We use panels with a near-vertical incidence to create a zero-offset image and a larger selection of body-wave-dominated panels to create virtual common-shot gathers. We process the gathers with a simple reflection seismology processing workflow to obtain stacked images. The zero-offset images show a relatively lower signal-to-noise ratio and only horizontal reflectors. The stacked images show slightly dipping reflectors and possibly lateral amplitude variations around the expected injection region. This could indicate a region of interest for future research into the reinjection reservoir.

Seismic interferometry (SI) retrieves new seismic responses, for example reflections, between either receivers or sources. When SI is applied to a reflection survey with active sources and receivers at the surface, non-physical (ghost) reflections are retrieved as well. Ghost reflections, retrieved from the correlation of two primary reflections or multiples from two different depth levels, are only sensitive to the properties in the layer that cause them to appear in the result of SI, such as velocity, density and thickness. We aim to use these ghost reflections for monitoring subsurface changes, to address challenges associated with detecting and isolating changes within the target layer in monitoring. We focus on the feasibility of monitoring pore-pressure changes in the Groningen gas field in the Netherlands using ghost reflections. To achieve this, we utilize numerical modelling to simulate scalar reflection data, deploying sources and receivers at the surface. To build up subsurface models for monitoring purposes, we perform an ultrasonic transmission laboratory experiment to measure S-wave velocities at different pore pressures. Applying SI by autocorrelation to the modelled data sets, we retrieve zero-offset ghost reflections. Using a correlation operator, we determine time differences between a baseline survey and monitoring surveys. To enhance the ability to detect small changes, we propose subsampling the ghost reflections before the correlation operator and using only virtual sources with a complete illumination of receivers. We demonstrate that the retrieved time differences between the ghost reflections exhibit variations corresponding to velocity changes inside the reservoir. This highlights the potential of ghost reflections as valuable indicators for monitoring even small changes. We also investigate the effect of the sources and receivers’ geometry and spacing and the number of virtual sources and receivers in retrieving ghost reflections with high interpretability resolution.

Many rocks contain planar heterogeneities, in the form of open fractures, veins and/or stylolites, but scarce data exist on how strength and fracture pattern formation is affected by the presence of a singular planar heterogeneity in an otherwise uniform matrix. The mechanics of stylolite-bearing and/or fractured limestone is of interest to several engineering applications, from quarries to subsurface gas or geothermal reservoirs. We have performed Brazilian Disc tests on pre-fractured Indiana limestone samples and Treuchtlinger Marmor discs which contain cohesive stylolites, investigating Brazilian test Strength and the resulting fracture pattern. All experiments were filmed, and where possible analyzed with particle image velocimetry. When viewed in 2D, the planar discontinuity was set at different rotation angles compared to the principal loading direction, where perpendicular to the loading direction is defined as 0⁰. The results show that all samples are weaker than their intact counterparts. For the pre-fractured Indiana limestone, there is 10–75% angle-dependent weakening. However, in the samples with a stylolite, strength is weakened by 35–75%, independent of direction. Several new cracks appeared when fracturing a stylolite-sample, where the orientation is heavily influenced by the stylolite orientation. The fracture pattern and associated stress drops are more complex for high angles. In these samples always more than one fracture formed, whereas in pre-fractured samples usually only one new fracture formed. This suggests a potential for more permeability increase when hydrofracturing a stylolite-rich interval. Comparison with Finite Element Models indicates that this difference in fracture pattern is caused by the strength contrast between the anastomosing stylolite zone and the matrix material, leading to stress concentrations effects. This causes (micro-) fracture nucleation to occur locally, promotes fracture coalescence and fracture growth at lower overall sample-load conditions compared to intact samples.

Geothermal energy production often involves use of corrosion inhibitors. We performed rock mechanical experiments (room temperature; confining pressure of 10/20/30 MPa) on typical reservoir rocks (Bentheim sandstone and Treuchtlinger limestone) in contact with two different inhibitor solutions or with demineralized water. The sandstone experiments show no discernible difference in rock strength between inhibitors or water, attributed to low quartz reactivity. The limestone experiments show a significant difference in rock strength (and Mohr–Coulomb envelope), dependent on inhibitor type, attributed to high carbonate reactivity. This implies that, depending on the reactivity of the rocks and local stress conditions, inhibitor leakage may lead to unpredicted reservoir failure.