Fault Observatory Research Articles

Noise attenuation in distributed acoustic sensing data using a guided unsupervised deep learning network

Distributed acoustic sensing (DAS) is a promising technology introducing a new paradigm in the acquisition of high-resolution seismic data. However, DAS data often show weak signals compared with the background noise, especially in challenging installation environments. In this study, we develop a new approach to denoise DAS data that leverages an unsupervised deep learning (DL) model, eliminating the need for labeled training data. The input DAS data undergo band-pass filtering to eliminate high-frequency content. Subsequently, a continuous wavelet transform (CWT) is performed, and the finest scale is used to guide the DL model in reconstructing the DAS signal. First, we extract 2D patches from the band-pass filtered data and the CWT scale of the data. Then, these patches are converted using an unrolling mechanism into 1D vectors to form the input of the DL model. A self-attention layer is included in each layer to extract the spatial relation between the band-pass filtered data and the CWT scale. Through an iterative process, the DL model tunes its parameters to suppress DAS noise, with the band-pass filtered data serving as the target for the network. The denoising performance of our framework is validated using field examples from the San Andreas Fault Observatory at Depth and Frontier Observatory for Research in Geothermal Energy data sets, where the data are recorded by a fiber-optic cable. Comparative analyses against three benchmark methods reveal the robust denoising performance of our framework.

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Abstract. Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr−1. The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

A study of fluid overpressure microstructures from the creeping segment of the San Andreas fault

Evidence of episodic fluid overpressure events noted in samples from the San Andreas Fault Observatory at Depth (SAFOD) have remained largely uncorrelated in terms of their collective significance for seismic history of the fault zone. The compositional and microstructural correlations sought in this study could shed light on questions about potential for major seismic events in the creeping segment of the SAF in central California. We used quantitative energy dispersive spectroscopy (EDS), Cathodoluminescence (CL) and Scanning Electron Microscope (SEM) imaging, and electron backscatter diffraction (EBSD) analysis to acquire geochemical and microstructural data from a suite of twenty SAFOD core samples including the damage zone and the active core of the fault. The results indicate intermittent coseismic fluid overpressure events that overprint the background aseismic creep across the fault. Analysis of trace elements and deformation in the coseismic calcite vein generations and their associated hydrothermal mineral phases indicate progressive uplift and exhumation followed by an asymmetric incursion of meteoric water into the damage zone. The same analysis suggests that the actively creeping intervals act as permeability barriers. Our results are in overall agreement with recent studies of the SAF in central California that indicate large seismic events have occurred intermittent with aseismic creep in recent geological time or suggest future potential for such events.

A novel infrastructure for the continuous monitoring of soil CO2 emissions: a case study at the alto Tiberina near fault observatory in Italy

Static and dynamic stress, along with earthquakes, can trigger the emission and migration of crustal fluids, as frequently observed on the surface and within the upper crust of tectonically active areas such as the northern Apennines of Italy. To investigate the origin of these fluids and their interconnection with the seismogenic process, we complemented The Alto Tiberina Near Fault Observatory (TABOO-NFO), a multidisciplinary monitoring infrastructure composed of a dense array of seismic, geodetic, strain, and radon sensors, with a proper geochemical network grounded on four soil CO2 flux monitoring stations and weather sensors, placed near the main vents of the superficial manifestations. The TABOO-NFO is a state-of-the-art monitoring infrastructure, which allows for studying various geophysical parameters connected to the deformation processes active along a crustal fault system dominated by the Alto Tiberina fault (ATF), which is a 60 km long normal fault dipping at a low angle (<15°–20°). The region is favourable for conducting geochemical studies, as it is characterised by the presence of over-pressurised fluids trapped at certain depths and superficial manifestations associated with the emission of large quantities of fluids. After describing the theoretical framework and the technological aspects based on which we developed the geochemical monitoring network, we described the data recorded in the first months. Over the studied period, the results showed that soil CO2 flux was primarily influenced by environmental parameters, and that the selected sites received a regular supply of deep-origin CO2.

Enhancing earthquake detection from distributed acoustic sensing data by coherency measure and moving-rank-reduction filtering

Distributed acoustic sensing (DAS) enables the recording of earthquake signals at an unprecedented low cost with high durability using fiber optic cables. It is often compromised by relatively lower data quality in single-channel measurements compared with traditional seismic receivers but compensated by high-density recordings from hundreds of channels per earthquake event. The multichannel nature of DAS data sets facilitates the applications of well-developed coherency-based denoising methods arising from reflection seismology for detecting more earthquakes. We first introduce a coherency measure for detecting earthquake signals from DAS data sets. Then, we apply a moving-rank-reduction (MRR) filter to enhance the DAS data quality so as to improve the earthquake detectability. The MRR filter is tailored from a rank-reduction filter that is widely used in processing multichannel reflection seismic data. We find that a simple band-pass or median filter is incapable of revealing weak signals generated from small-magnitude or far-away earthquake events, whereas the MRR filter significantly improves the signal-to-noise ratio that enables the detection of those weak signals. We apply the MRR method and the coherency measure on the San Andreas Fault Observatory at Depth DAS data sets for denoising and earthquake detection. As a result, our framework detects all 31 cataloged events, outperforming the previous detection of 25 events using the same data set.

Comparing and integrating artificial intelligence and similarity search detection techniques: application to seismic sequences in Southern Italy

SUMMARYUnderstanding mechanical processes occurring on faults requires detailed information on the microseismicity that can be enhanced today by advanced techniques for earthquake detection. This problem is challenging when the seismicity rate is low and most of the earthquakes occur at depth. In this study, we compare three detection techniques, the autocorrelation FAST, the machine learning EQTransformer, and the template matching EQCorrScan, to assess their ability to improve catalogues associated with seismic sequences in the normal fault system of Southern Apennines (Italy) using data from the Irpinia Near Fault Observatory (INFO). We found that the integration of the machine learning and template matching detectors, the former providing templates for the cross-correlation, largely outperforms techniques based on autocorrelation and machine learning alone, featuring an enrichment of the automatic and manual catalogues of factors 21 and 7, respectively. Since output catalogues can be polluted by many false positives, we applied refined event selection based on the cumulative distribution of their similarity level. We can thus clean up the detection lists and analyse final subsets dominated by real events. The magnitude of completeness decreases by more than one unit compared to the reference value for the network. We report b-values associated with sequences smaller than the average, likely corresponding to larger differential stresses than for the background seismicity of the area. For all the analysed sequences, we found that main events are anticipated by foreshocks, indicating a possible preparation process for main shocks at subkilometric scales.

The Near Fault Observatory community in Europe: a new resource for faulting and hazard studies

The Near Fault Observatories (NFOs) community is one of the European Plate Observing System (EPOS, http://www.epos-eu.org) Thematic Communities, today consisting of six research infrastructures that operate in regions characterised by high seismic hazard originating from different tectonic regimes. Earthquakes respond to complex natural systems whose mechanical properties evolve over time. Thus, in order to understand the multi-scale, physical/chemical processes responsible for the faulting that earthquakes occur on, it is required to consider phenomena that intersect different research fields, i.e., to put in place multidisciplinary monitoring. Hence, NFOs are grounded on modern and multidisciplinary infrastructures, collecting near fault high resolution raw data that allows generation of innovative scientific products. The NFOs usually complement regional backbone networks with a higher density distribution of seismic, geodetic, geochemical and other geophysical sensors, at surface and sometimes below grade. These dense and modern networks of multi-parametric sensors are sited at and around active faults, where moderate to large earthquakes have occurred in the past and are expected in the future. They continuously monitor the underlying Earth instability processes over a broad time interval. Data collected at each NFO results in an exceptionally high degree of knowledge of the geometry and parameters characterizing the local geological faults and their deformation pattern. The novel data produced by the NFO community is aggregated in EPOS and is made available to a diverse set of stakeholders through the NFO Federated Specific Data Gateway (FRIDGE). In the broader domain of the Solid Earth sciences, NFOs meet the growing expectations of the learning and communication sectors by hosting a large variety of scientific information about earthquakes as a natural phenomenon and a societal issue. It represents the EPOS concept and objective of aggregating and harmonising the European research infrastructures capabilities to facilitate broader scientific opportunity. The NFOs are at the cutting edge of network monitoring. They conduct multidisciplinary experiments for testing multi-sensor stations, as well as realise robust and ultra-low latency, transmission systems that can routinely accommodate temporary monitoring densification. The effort to continuously upgrade the technological efficiency of monitoring systems positions the NFO at the centre of marketing opportunities for the European enterprises devoted to new sensor technology. The NFOs constitute ideal test beds for generating expertise on data integration, creating tools for the next generation of multidisciplinary research, routine data analysis and data visualization. In particular focus is often on near-real time tools and triggering alarms at different levels are tested and implemented, strengthening the cooperation with the Agencies for risk management. NFOs have developed innovative operational actions such as the Testing Centre for Earthquake Early Warning and Source Characterisation (CREW) and detailed fast ground shaking and damage characterization. Complementing the recent growth of modern laboratory and computational models, the NFOs can provide interdisciplinary observations of comparable high resolution to describe the behaviour of fault slip over a vast range of spatial and temporal scales and aiding to provide more accurate earthquake hazard characterizations.

History of earthquakes along the creeping section of the San Andreas fault, California, USA

Abstract Creeping faults are difficult to assess for seismic hazard because they may participate in rupture even though they likely cannot nucleate large earthquakes. The creeping central section of the San Andreas fault in California (USA) has not participated in a historical large earthquake; however, earthquake ruptures nucleating in the locked northern and southern sections may propagate through the creeping section. We used biomarker thermal maturity and K/Ar dating on samples from the San Andreas Fault Observatory at Depth to look for evidence of earthquakes. Biomarkers show evidence of many earthquakes with displacements >1.5 m in and near a 3.5-m-wide patch of the fault. We show that K/Ar ages decrease with thermal maturity, and partial resetting occurs during coseismic heating. Therefore, measured ages provide a maximum constraint on earthquake age, and the youngest earthquakes here are younger than 3 Ma. Our results demonstrate that creeping faults may host large earthquakes over longer time scales.

Spatial-temporal characterization of the San Andreas Fault by fault-zone trapped waves at seismic experiment site, Parkfield, California

In this article, we review our previous research for spatial and temporal characterizations of the San Andreas Fault (SAF) at Parkfield, using the fault-zone trapped wave (FZTW) since the middle 1980s. Parkfield, California has been taken as a scientific seismic experimental site in the USA since the 1970s, and the SAF is the target fault to investigate earthquake physics and forecasting. More than ten types of field experiments (including seismic, geophysical, geochemical, geodetic and so on) have been carried out at this experimental site since then. In the fall of 2003, a pair of scientific wells were drilled at the San Andreas Fault Observatory at Depth (SAFOD) site; the main-hole (MH) passed a ~200-m-wide low-velocity zone (LVZ) with highly fractured rocks of the SAF at a depth of ~3.2 km below the wellhead on the ground level (Hickman et al., 2005; Zoback, 2007; Lockner et al., 2011). Borehole seismographs were installed in the SAFOD MH in 2004, which were located within the LVZ of the fault at ~3-km depth to probe the internal structure and physical properties of the SAF. On September 28 2004, a M6 earthquake occurred ~15 km southeast of the town of Parkfield. The data recorded in the field experiments before and after the 2004 M6 earthquake provided a unique opportunity to monitor the co-mainshock damage and post-seismic heal of the SAF associated with this strong earthquake. This retrospective review of the results from a sequence of our previous experiments at the Parkfield SAF, California, will be valuable for other researchers who are carrying out seismic experiments at the active faults to develop the community seismic wave velocity models, the fault models and the earthquake forecasting models in global seismogenic regions.

Geospatial Management and Analysis of Microstructural Data from San Andreas Fault Observatory at Depth (SAFOD) Core Samples

Core samples obtained from scientific drilling could provide large volumes of direct microstructural and compositional data, but generating results via the traditional treatment of such data is often time-consuming and inefficient. Unifying microstructural data within a spatially referenced Geographic Information System (GIS) environment provides an opportunity to readily locate, visualize, correlate, and apply remote sensing techniques to the data. Using 26 core billet samples from the San Andreas Fault Observatory at Depth (SAFOD), this study developed GIS-based procedures for: 1. Spatially referenced visualization and storage of various microstructural data from core billets; 2. 3D modeling of billets and thin section positions within each billet, which serve as a digital record after irreversible fragmentation of the physical billets; and 3. Vector feature creation and unsupervised classification of a multi-generation calcite vein network from cathodluminescence (CL) imagery. Building on existing work which is predominantly limited to the 2D space of single thin sections, our results indicate that a GIS can facilitate spatial treatment of data even at centimeter to nanometer scales, but also revealed challenges involving intensive 3D representations and complex matrix transformations required to create geographically translated forms of the within-billet coordinate systems, which are suggested for consideration in future studies.

Passive Seismic Imaging of Near Vertical Structures Around the SAFOD Site, California, Jointly Using Scattered P and SH Waves

AbstractTo better illuminate structural discontinuities around the San Andreas Fault Observatory at Depth (SAFOD) site, following the previous study of Zhang et al. (2009), we have extended the Generalized Radon Transform (GRT) method to jointly use scattered P and SH waves from abundant local earthquakes recorded by a local dense seismic network. A hybrid imaging condition is applied to extract consistent structures from scattered P and SH images. In this way, more robust structure information can be determined from separate images with potential artifacts. For the transverse component waveforms, coherent S‐S scattered waves after the direct SH waves can be clearly identified, and they are less interfered with by other scattered and converted waves compared to scattered P‐P waves on the vertical component. Similar to Zhang, Wang, et al. (2009), near‐vertical reflectors are imaged on both sides of the San Andreas Fault (SAF) with scattered SH waves, which is generally consistent with the imaging results using scattered P‐P waves from local earthquakes and a wide‐angle active seismic reflection profile. Compared to the P‐P scattering imaging result, the imaging result using scattered SH waves has higher resolution due to shorter S wavelength and cleaner S‐S scattered waveforms, and the SAF, as well as other reflectors around it, is better resolved. Although the resolution of the joint image obtained by combining separate imaging results from different scattered waves may be degraded, it is able to more robustly characterize strong structure discontinuities using passive seismic sources.

Core surprise: What's inside a plate boundary?

Despite the fact that 90% of global seismicity occurs at plate boundary faults, our understanding of their internal structure is lacking. It’s not easy to see inside a plate boundary fault – typically composed of a high-strain fault core surrounded by a fractured damage zone – and when we can, it often requires expensive drilling projects that yield limited information on the internal structure of the whole fault. Understanding the internal structure of large faults is crucial, because their chemical and mechanical properties control how and where earthquakes rupture, nucleate and propagate. This in turn limits the size of the earthquake or the amount of radiated seismic energy, and consequently the severity of surface damage. The 1999 magnitude 7.7 earthquake along the Chelungpu plate boundary fault, for example – the second deadliest earthquake in Taiwan’s recorded history – saw significant variations in slip and ground motion at different locations along the fault which resulted in large local variations in casualties and damage. Subsequent field investigations related these variations to changes in the fault’s structure (i.e., clay width, geometry), which in turn controlled how the fault moved.

Velocity‐Based Earthquake Detection Using Downhole Distributed Acoustic Sensing—Examples from the San Andreas Fault Observatory at Depth

AbstractConventional seismographic networks sparsely sample the wavefields excited by earthquakes. Thus, standard event detection is conducted by analyzing separate stations and merging their results. Emerging distributed acoustic sensing recording technologies allow for unbiased spatial sampling of the wavefield and, as a result, array‐based processing of the recorded signals. Using a cemented fiber in the San Andreas Fault Observatory at Depth main hole, 800 virtual receivers are sampled at a 1 m interval from the surface to 800 m depth. Recorded earthquakes are approximated as plane waves reaching the bottom of the array first. Following this assumption, the relative travel times of the recorded event depend on the local velocity at the array location and the angle of incidence at which the planar wavefront reaches it. Given the seismic velocity, a newly proposed detection algorithm amounts to a single‐parameter scan of the incidence angle and measurement of data coherency along the different possible travel‐time curves. Using the entire recording array, a much higher effective signal‐to‐noise ratio can be obtained when compared to individual channel processing. About 20 days of recorded seismic activity from the San Andreas Fault is analyzed. Using a downhole single array, the majority of cataloged events in the area are detected. In addition, a previously unknown event is unveiled. We estimate its magnitude at roughly −0.5.

Detection of weak seismic sequences based on arrival time coherence and empiric network detectability: an application at a near fault observatory

SUMMARYMicroseismic monitoring is a primary tool for understanding and tracking the progress of mechanical processes occurring in active rock fracture systems. In geothermal or hydrocarbon fields or along seismogenic fault systems, the detection and location of microseismicity facilitates resolution of the fracture system geometry and the investigation of the interaction between fluids and rocks, in response of stress field perturbations. Seismic monitoring aims to detect locate and characterize seismic sources. The detection of weak signals is often achieved at the cost of increasing the number of false detections, related to transient signals generated by a range of noise sources, or related to instrumental problems, ambient conditions or human activity that often affect seismic records. A variety of fast and automated methods has been recently proposed to detect and locate microseismicity based on the coherent detection of signal anomalies, such as increase in amplitude or coherent polarization, at dense seismic networks. While these methods proved to be very powerful to detect weak events and to reduce the magnitude of completeness, a major problem remains to discriminate among weak seismic signals produced by microseismicity and false detections. In this work, the microseimic data recorded along the Irpinia fault zone (Southern Apennines, Italy) are analysed to detect weak, natural earthquakes using one of such automated, migration-based, method. We propose a new method for the automatic discrimination of real vs false detections, which is based on empirical data and information about the detectability (i.e. detection capability) of the seismic network. Our approach allows obtaining high performances in detecting earthquakes without requiring a visual inspection of the seismic signals and minimizing analyst intervention. The proposed methodology is automated, self-updating and can be tuned at different success rates.

Seismic Velocity Estimation Using Passive Downhole Distributed Acoustic Sensing Records: Examples From the San Andreas Fault Observatory at Depth

AbstractStructural imaging and event location require an accurate estimation of the seismic velocity. However, active seismic surveys used to estimate it are expensive and time‐consuming. During the last decade, fiber‐optic‐based distributed acoustic sensing has emerged as a reliable, enduring, and high‐resolution seismic sensing technology. We show how downhole distributed acoustic sensing passive records from the San Andreas Fault Observatory at Depth can be used for seismic velocity estimation. Using data recorded from earthquakes propagating near‐vertically, we compute seismic velocities using first‐break picking as well as slant stack decomposition. This methodology allows for the estimation of both P and S wave velocity models. We also use records of the ambient seismic field for interferometry and P wave velocity model extraction. Results are compared to a regional model obtained from surface seismic as well as a conventional downhole geophone survey. We find that using recorded earthquakes, we obtain the highest P wave model resolution. In addition, it is the only method that allows for S wave velocity estimation. Computed P and S models unravel three distinct areas at the depth range of 50‐750 m, which were not present in the regional model. In addition, we find high VP/VS values near the surface and a possible VP/VS anomaly about 500 m deep. We confirm its existence by observing a strong S‐P mode conversion at that depth.

Serpentinite‐Rich Gouge in a Creeping Segment of the Bartlett Springs Fault, Northern California: Comparison With SAFOD and Implications for Seismic Hazard

AbstractAn exposure of a creeping segment of the Bartlett Springs Fault (BSF), part of the San Andreas Fault system in northern California, is a ~1.5‐m‐wide zone of serpentinite‐bearing fault gouge cutting through Late Pleistocene fluvial deposits. The fault gouge consists of porphyroclasts of antigorite serpentinite, talc, chlorite, and tremolite‐actinolite, along with some Franciscan metamorphic rocks, in a matrix of the same materials. The Mg‐mineral assemblage is stable at temperatures above 250–300 °C. The BSF gouge is interpreted to have been tectonically incorporated into the fault from depths near the base of the seismogenic zone and to have risen buoyantly to the surface where it is now undergoing right‐lateral displacement. The ultramafic‐rich composition, frictional properties, and inferred mode of emplacement of the BSF serpentinitic gouge correspond to those of the creeping traces of the San Andreas Fault identified in the SAFOD (San Andreas Fault Observatory at Depth) drill hole. This suggests a common origin for creep at both locations. A tectonic model for the source of the ultramafic‐rich materials in the BSF is proposed that potentially could explain the distribution of creep throughout the northernmost San Andreas Fault system.

Continuous Measurement of Stress‐Induced Travel‐Time Variations at SAFOD

Author(s): Yang, C; Niu, F; Daley, TM; Taira, T | Abstract: In situ stress measurement at seismogenic depth is critically important for deciphering fault zone processes. In this study, we conducted a second active-source crosswell field experiment at the Parkfield San Andreas Fault Observatory at Depth (SAFOD) drill site to investigate the detectability of stress-induced seismic velocity changes at the top part of the seismogenic zone. We employed the same configuration of our previous experiments, which deployed a piezoelectric source and a three-component (3C) accelerometer at 1 km deep inside the pilot and main holes, respectively. We also added a hydrophone, which is attached to the source, to monitor the repeatability of the source waveforms. Over a 40-day recording period, we confirmed an ∼0:04% travel-time variation in S wave and coda that roughly follows the fluctuation of barometric pressure. We attributed this correlation to stress sensitivity of seismic velocity and the stress sensitivity is estimated to be 2:0 × 10 −7 Pa −1 , which is approximately two orders of magnitude higher than those measured in laboratory with dry rock samples, but is consistent with our previous results. Our results confirm the hypothesis that substantial cracks and/or pore spaces exist at seismogenic depths and thus may be used to monitor the subsurface stress field with active-source crosswell seismic.

The NextM~6 Event in Parkfield Implied by a Physical Model Linking Interseismic, Coseismic, and Postseismic Phase

AbstractDetermination of the constitutive frictional parameters is crucial for describing the dynamic behaviors of natural faults. In this study, we explore the along‐strike variation in fault friction along the San Andreas Fault in Parkfield, where characteristicM6 earthquakes occur nearly every 25 years. Using displacement data measured over decades encompassing the interseismic, coseismic, and postseismic phase of the 2004M6 event, which involve inherent fault rheological responses to tectonic loading, we construct a physical model based on the rate‐dependent friction law to describe the fault slip and stress evolution. Our modeled friction rate parameter is consistent with available laboratory measurements using fault cores collected from the San Andreas Fault Observatory at Depth. Based on the physical model linking interseismic, coseismic, and postseismic phase, we predict that the 25‐year slip deficit following the 2004M6 earthquake will predominantly accumulate in the area to the south of the 2004 event hypocenter, with an energy equivalent to aM6.0 event. Together with the slip deficit accumulated over the 25 years prior to the 2004 event that was unreleased by this event, the total slip deficit amounts to a moment magnitude ofM6.35 ± 0.11, suggesting the upper limit of the seismic release level until 2029. The proposed model (1) involves only two free quantities (the friction rate parameter and velocity coefficient used to determine the effective normal stress); (2) is well constrained by available interseismic, coseismic, and postseismic displacement data; and (3) highlights along‐strike variations in the fault friction.

A study of secondary pyrite deformation and calcite veins in SAFOD damage zone with implications for aseismic creep deformation mechanism at depths >3 km

Previous studies of the San Andreas Fault damage-zone samples from the San Andreas Fault Observatory at Depth (SAFOD) have identified a variety of tectonic microstructures including pressure solution cleavage, calcite-sealed fractures vein fabric, and pyrite and anhydrite hydrothermal fracture sealing. Understanding the deformation provenance of the damage zone rocks and operative deformation mechanism(s) based on preserved microstructures provide insight into overall deformation behavior of the entire seismogenic zone in the creeping section of this transform fault. We analysed the deformation of hydrothermal secondary pyrite in connection with network of calcite veins in a sample of foliated ultracataclasites bordering the actively creeping Southwestern Deforming Zone (SDZ), using SEM, EBSD and CL microscopy. The results show that calcite veins associated with the pressure solution cleavage are crosscut by the secondary pyrite deformed under a range of P-T conditions. Relatively undeformed secondary pyrite is found sealing implosion microbreccia. Our review of previously available data indicates that the damage zone rocks may represent a collage of structural and compositional domains from both locked and creeping sections of the SAF. This interpretation together with results of this study suggest that weak-clay frictional deformation mechanism(s) is likely to be the predominant aseismic creep mechanism at depths below the SAFOD.

Grain size‐dependent strength of phyllosilicate‐rich gouges in the shallow crust: Insights from the SAFOD site

AbstractThe San Andreas Fault Observatory at Depth (SAFOD) drilling project directly sampled a transitional (between creeping and locked) segment of the San Andreas Fault at 2.7 km depth. At the site, changes in strain rate occur between periods of coseismic slip (>10−7 s−1) and interseismic creep (10−10 s−1) over decadal scales (~30 years). Microstructural observations of core retrieved from the SAFOD site show throughgoing fractures and gouge‐rich cores within the fractures, evidence of predominantly brittle deformation mechanisms. Within the gouge‐rich cores, strong phases show evidence of deformation by pressure solution once the grain size is reduced to a critical effective grain size. Models of pressure solution‐accommodated creep for quartz‐phyllosilicate mixtures indicate that viscous weakening of quartz occurs during the interseismic period once a critical effective grain size of 1 μm is achieved, consistent with microstructural observations. This causes pronounced weakening, as the strength of the mixture is then controlled by the frictional properties of the phyllosilicate phases. These results have pronounced implications for the internal deformation of fault zones in the shallow crust, where at low strain rates, deformation is accommodated by both viscous and brittle deformation mechanisms. As strain rates increase, the critical effective grain size for weakening decreases, localizing deformation into the finest‐grained gouges until deformation can no longer be accommodated by viscous processes and purely brittle failure occurs.