3D architecture and complex behavior along the simple central San Andreas fault

The central San Andreas Fault (CSAF) exhibits a simple linear large-scale fault geometry, yet seismic and aseismic deformation features vary in a complex way along the fault. Here we investigate fault zone behaviors using geodetic observation, seismicity and microearthquake focal mechanisms. We employ an improved focal-mechanism characterization method using relative earthquake radiation patterns on 75,164 Ml ≥ 1 earthquakes along a 2-km-wide, 190-km-long segment of the CSAF, from 1984 to 2015. The data reveal the 3D fine-scale structure and interseismic kinematics of the CSAF. Our findings indicate that the first-order spatial variations in interseismic fault creep rate, creep direction, and the fault zone stress field can be explained by a simple fault coupling model. The inferred 3D mechanical properties of a mechanically weak and poorly coupled fault zone provide a unified understanding of the complex fine-scale kinematics, indicating distributed slip deficits facilitating small-to-moderate earthquakes, localized stress heterogeneities, and complex multi-scale ruptures along the fault. Through this detailed mapping, we aim to relate the fine-scale fault architecture to potential future faulting behavior along the CSAF.

De-risking fault leakage risk and containment integrity for subsurface storage applications

The subsurface is pivotal in the energy transition, for the sequestration of CO2 and energy storage. It is crucial to understand to what extent geological faults may form leakage pathways that threaten the containment integrity of these projects. Fault flow behavior has been studied in the context of hydrocarbon development, supported by observations from wells drilled through faults, but such observations are rare in geoenergy projects. Focusing on mechanical behavior as early indicator of potential leakage risks, a probabilistic Coulomb Failure Stress workflow is developed and demonstrated using data from the Decatur CO2 sequestration project to rank faults based on their containment risk. The analysis emphasizes the importance of fault throw relative to reservoir thickness and pore pressure change in assessing reactivation risks. Integrating this mechanical assessment with geological and dynamic fault analyses contributes to derisking fault containment for geoenergy applications, providing valuable insights for the successful development of subsurface storage projects.

Hydro-mechanical coupling characteristics and weakening mechanisms of filling joint resulting from water injection

The injection of fluids into fault gouges in rock formations disturbs the in situ stress conditions, leading to fault slip and increasing the risk of inducing earthquakes. The weakening effect and the permeation of the injected fluid can be influenced significantly by the presence of fault gouges. To investigate this issue, the hydro-mechanical characteristics of fault gouges were evaluated using physical tests to study the combined effects of coupling injecting water and shear deformation. We propose a new experimental apparatus that allows us to measure the spatial distribution of the thickness of a gouge layer sample under combined conditions of shearing and water injection, using 3D scanning technology to evaluate the primary flow path. The test results showed that injecting water had a significant effect in reducing the maximum shear strength, but the degree to which the strength was affected varied with the gouge fill material. The effect of shear deformation is that it will increase the inhomogeneity of the thickness distribution and therefore the distribution of injected water along the fault is not uniformly radial. Additionally, the properties of gouge fill material have an important influence on flow characteristics during fault slipping.

Effects of breccia and water contents on the mechanical properties of fault-core-zone materials

Determining the mechanical properties of fault-core-zone materials is challenging because of the low strength of such materials, which affects field sampling, specimen preparation, and laboratory testing. We overcame this problem by preparing and testing mechanical properties of 132 artificial fault-core-zone specimens consisting of mixtures of breccia, sand, clay, and water. The unconfined compressive strength (UCS), elastic modulus (E), and penetration resistance value (PRV) of these fault-core-zone materials were measured, and the effects of breccia content and water content on mechanical properties were assessed. Results show that UCS is inversely proportional to breccia content and water content, and that E is inversely proportional to water content. Furthermore, the inverse relationship of UCS with water content varies with breccia content. UCS is proportional to both PRV and E, and the relationship for each varies with breccia content. High coefficients of determination (R² = 0.62–0.88) between the parameters suggest that breccia content, water content, and PRV are potentially useful parameters for estimating the mechanical properties of fault core zones.

Signature of transition to supershear rupture speed in the coseismic off-fault damage zone

Most earthquake ruptures propagate at speeds below the shear wave velocity within the crust, but in some rare cases, ruptures reach supershear speeds. The physics underlying the transition of natural subshear earthquakes to supershear ones is currently not fully understood. Most observational studies of supershear earthquakes have focused on determining which fault segments sustain fully grown supershear ruptures. Experimentally cross-validated numerical models have identified some of the key ingredients required to trigger a transition to supershear speed. However, the conditions for such a transition in nature are still unclear, including the precise location of this transition. In this work, we provide theoretical and numerical insights to identify the precise location of such a transition in nature. We use fracture mechanics arguments with multiple numerical models to identify the signature of supershear transition in coseismic off-fault damage. We then cross-validate this signature with high-resolution observations of fault zone width and early aftershock distributions. We confirm that the location of the transition from subshear to supershear speed is characterized by a decrease in the width of the coseismic off-fault damage zone. We thus help refine the precise location of such a transition for natural supershear earthquakes.

Weakening Mechanisms in a Basalt-Hosted Subduction Megathrust Fault Segment, Southern Alaska

Basaltic and gabbroic rocks that define the seafloor have been suggested to act as sources of rheological heterogeneity during subduction, with the capacity to enhance or dampen seismicity. Despite this, relatively little is known from the rock record regarding the progression and conditions of mafic oceanic crust deformation during subduction, particularly in the shallow megathrust region of the seismogenic zone. We describe subduction-related deformation structures and characterize deformation conditions from an exhumed, basalt-hosted megathrust in the Chugach accretionary complex of south-central Alaska. Rocks in the Chugach preserve a record of seafloor mineralogical changes from pre-subduction, hydrothermal circulation that produced sheet silicates with a lower frictional strength than intact basalt. Pre-subduction alteration also served to introduce hydrous phases that can expel water during deformation and raise the pore fluid pressure. Once strain localized within basalts onto a megathrust fault plane at lithostatic pore fluid pressures, the basalt weakened further through a combination of cataclasis, dilatational shear fracturing, and slip on chlorite-rich shear bands. This process occurred in a narrower fault zone, and at higher maximum differential stress and greater pore fluid pressure fluctuations than recorded in some sediment-hosted megathrusts at similar pressure and temperature conditions. Our data indicate that when the lower plate contains basalt bathymetric features, basalt dismembers during subduction into a chlorite-rich fault gouge that surrounds lenses or slices of intact, less-altered basalt.

Seismic Slip-Pulse Experiments Simulate Induced Earthquake Rupture in the Groningen Gas Field

Rock materials show dramatic dynamic weakening in large-displacement (m), high-velocity (∼1 m/s) friction experiments, providing a mechanism for the generation of large, natural earthquakes. However, whether such weakening occurs during induced M3-4 earthquakes (dm displacements) is unknown. We performed rotary-shear experiments on simulated fault gouges prepared from the source-, reservoir- and caprock formations present in the seismogenic Groningen gas field (Netherlands). Water-saturated gouges were subjected to a slip pulse reaching a peak circumferential velocity of 1.2-1.7 m/s and total displacements of 13-20 cm, at 2.5-20 MPa normal stress. The results show 22%-81% dynamic weakening within 5-12 cm of slip, depending on normal stress and gouge composition. At 20 MPa normal stress, dynamic weakening from peak friction coefficients of 0.4-0.9 to 0.19-0.27 was observed, probably through thermal pressurization. We infer that similar effects play a key role during induced seismic slip on faults in the Groningen and other reservoir systems.

A unified first-order hyperbolic model for nonlinear dynamic rupture processes in diffuse fracture zones

Earthquake fault zones are more complex, both geometrically and rheologically, than an idealized infinitely thin plane embedded in linear elastic material. To incorporate nonlinear material behaviour, natural complexities and multi-physics coupling within and outside of fault zones, here we present a first-order hyperbolic and thermodynamically compatible mathematical model for a continuum in a gravitational field which provides a unified description of nonlinear elasto-plasticity, material damage and of viscous Newtonian flows with phase transition between solid and liquid phases. The fault geometry and secondary cracks are described via a scalar function ξ ∈ [0, 1] that indicates the local level of material damage. The model also permits the representation of arbitrarily complex geometries via a diffuse interface approach based on the solid volume fraction function α ∈ [0, 1]. Neither of the two scalar fields ξ and α needs to be mesh-aligned, allowing thus faults and cracks with complex topology and the use of adaptive Cartesian meshes (AMR). The model shares common features with phase-field approaches, but substantially extends them. We show a wide range of numerical applications that are relevant for dynamic earthquake rupture in fault zones, including the co-seismic generation of secondary off-fault shear cracks, tensile rock fracture in the Brazilian disc test, as well as a natural convection problem in molten rock-like material. This article is part of the theme issue 'Fracture dynamics of solid materials: from particles to the globe'.

Nucleation of Stick-Slip Instability Within a Large-Scale Experimental Fault: Effects of Stress Heterogeneities Due to Loading and Gouge Layer Compaction

Geodetic observations and large-scale laboratory experiments show that seismic instability is preceded by slow slip within a finite nucleation zone. In laboratory experiments rupture nucleation is studied mostly using bare (rock) interfaces, whereas upper crustal faults are typically filled with gouge. To investigate effects of gouge on rupture nucleation, we performed a biaxial shearing experiment on a 350 mm long saw-cut fault filled with gypsum gouge, at room temperature and a minimum horizontal stress σ 2 = 0.3-5 MPa. The gouge layer was sandwiched between polymethylmethacrylate (PMMA) plates For reference also a fault without gouge was deformed. Strain gauges and Digital Image Correlation were used to monitor the deformation field along the fault zone margins. Stick-slip behavior occurred on both the gouge-filled fault and the PMMA fault. Nucleation of instability on the PMMA fault persistently occurred from one location 2/3 to 3/4 along the fault adjacent to a slow slip zone at the fault end, but nucleation on the gouge-filled fault was more variable, nucleating at the ends and/or at approximately 2/3 along the fault, with precursory slip occurring over a large fraction of the fault. Nucleation correlated to regions of high average fault stress ratio τ/σ n , which was more variable for the gouge-filled fault due to small length scale variations in normal stress caused by heterogeneous gouge compaction. Rupture velocities and slip rates were lower for the gouge-filled fault than for the bare PMMA fault. Stick-slip persisted when σ 2 was lowered and the nucleation zone length increased, expanding from the center to the sample ends before transitioning into instability.

Effect of Mineral Processes and Deformation on the Petrophysical Properties of Soft Rocks during Active Faulting

We have studied damage zones of two active faults, Baza and Padul faults in Guadix-Baza and Granada basins, respectively, in South Spain. Mineral and microstructural characterization by X-ray diffraction and field emission electron microscopy studies have been combined with structural fieldwork and in situ measurements of rock properties (permeability and Young’s modulus) to find out the relation between deformation behavior, mineral processes, and changes in the soft rock and sediment properties produced by fluid flow during seismic cycles. Our results show that microsealing produced by precipitation of dolomite and aragonite along fractures in the damage zone of Baza Fault reduces the permeability and increases the Young’s modulus. In addition, deformation bands formed in sediments richer in detrital silicates involved cataclasis as deformation mechanism, which hamper permeability of the sediments. In the Granada Basin, the calcarenitic rocks rich in calcite and clays in the damage zone of faults associated to the Padul Fault are characterized by the presence of stylolites without any carbonate cement. On the other hand, marly lithofacies affected by faults are characterized by the presence of disaggregation bands that involve cracking and granular flow, as well as clay smear. The presence of stylolites and deformation bands in these rocks reduces permeability.

Transformation Mechanism of a Fault and Its Associated Microstructures in Low-Porosity Rocks: A Case Study of the Tanan Depression in the Hailar-Tamtsag Basin

Faults have complex internal structures, which can be divided into the fault core and the damaged zone. During the fault formation process, the damaged zone will develop a large number of fractures cross cutting the adjacent host rocks, where the size and density of fractures would decrease as a function of distance from the fault core. Statistics show that compared with the host rocks, the rock’s secondary porosity may significantly get improved until it is 5%–10% higher than the host rock’s, and the permeability could be increased from zero to six orders of magnitude higher where fractures are developed. Based on the distribution of the fault core and damaged zone, we established an idealized geologic model to analyze the influence of the fault associated microstructures on the fluid flow efficiency. The results demonstrated that the reservoir property will be effectively improved when the fractures are developed to a certain magnitude, which will provide an advantageous conduit for fluid flow. The physical properties of wells are significantly different between the reservoirs transformed by the fractures and the reservoirs that have not been affected. The reservoir unit near faults has been modified by the associated fractures improving the reservoir quality. In addition, the portions of the reservoir farther from the fault core are less influenced and retain their initial poor characteristics. In order to evaluate the enhanced reservoir properties caused by faults and associated fractures, we performed statistical analysis of valid porosity ratio of the Nantun reservoirs as a function of the distance from 49 wells to faults in the research area. The results of this study demonstrated that connected fractures enhance the properties of the reservoirs and there is a distinct range that separates oil producing wells and water producing ones. This phenomenon verifies that faults are important and must be considered carefully during the exploration and production for hydrocarbon to provide higher quality reservoirs.

The Nanoscale Basis of CO2 Trapping for Geologic Storage

Carbon capture and storage (CCS) is likely to be a critical technology to achieve large reductions in global carbon emissions over the next century. Research on the subsurface storage of CO2 is aimed at reducing uncertainties in the efficacy of CO2 storage in sedimentary rock formations. Three key parameters that have a nanoscale basis and that contribute uncertainty to predictions of CO2 trapping are the vertical permeability kv of seals, the residual CO2 saturation Sg,r in reservoir rocks, and the reactive surface area ar of silicate minerals. This review summarizes recent progress and identifies outstanding research needs in these areas. Available data suggest that the permeability of shale and mudstone seals is heavily dependent on clay fraction and can be extremely low even in the presence of fractures. Investigations of residual CO2 trapping indicate that CO2-induced alteration in the wettability of mineral surfaces may significantly influence Sg,r. Ultimately, the rate and extent of CO2 conversion to mineral phases are uncertain due to a poor understanding of the kinetics of slow reactions between minerals and fluids. Rapidly improving characterization techniques using X-rays and neutrons, and computing capability for simulating chemical interactions, provide promise for important advances.

Fault seal calibration: a brief review

Calibration is a necessary step in the workflow for prediction of fault seal because there is no direct way to detect the hydraulic behaviour of a fault at the scale of a hydrocarbon trap. Over the last 20 years two general approaches have been developed: