Search

Imaging the structure of major fault zones is essential for our understanding of crustal deformations and their implications on seismic hazards. Investigating such complex regions presents several issues, including the variation of seismic velocity due to the diversity of geological units and the cumulative damage caused by earthquakes. Conventional migration techniques are in general strongly sensitive to the available velocity model. Here we apply a passive matrix imaging approach which is robust to the mismatch between this model and the real seismic velocity distribution. This method relies on the cross‐correlation of ambient noise recorded by a geophone array. The resulting set of impulse responses form a reflection matrix that contains all the information about the subsurface. In particular, the reflected body waves can be leveraged to: (a) determine the transmission matrix between the Earth's surface and any point in the subsurface; (b) build a confocal image of the subsurface reflectivity with a transverse resolution only limited by diffraction. As a study case, we consider seismic noise (0.1–0.5 Hz) recorded by the Dense Array for Northern Anatolia that consists of 73 stations deployed for 18 months in the region of the 1999 Izmit earthquake. Passive matrix imaging reveals the scattering structure of the crust and upper mantle around the North Anatolian Fault zone over a depth range of 60 km. The results show that most of the scattering is associated with the Northern branch that passes throughout the crust and penetrates into the upper mantle. Plain Language Summary: Investigating the structure of major fault zones is important to understand the deformations of the Earth's crust and their potential impact on future earthquakes. However, the large variations of seismic velocity between different geological units drastically hampers the ability to image those areas. The North Anatolian Fault is no exception since it splits into two branches separating three major geological blocks in the region of the 1999 Izmit earthquake. To image the deep structure of this fault, seismic noise is here exploited to retrieve information on the body waves reflected by the underground heterogeneities. A three‐dimensional map of the subsurface reflectivity is then built by applying methods originally developed in optical microscopy and ultrasound for deep imaging inside complex media. The inner structure of the Earth around the North Anatolian Fault is revealed. In particular, the Northern branch is shown to exhibit a strong damage pattern and a deep penetration inside the upper mantle. Key Points: A matrix approach of passive seismology provides 3D imaging of the North Anatolian Fault, Turkey, over a depth range of 60 kmTransmission matrix between the Earth surface and the underground is extracted by means of spatio‐temporal correlations of seismic noiseTime reversal of this transmission matrix enables a sharp compensation of the seismic velocity variations in the fault area [ABSTRACT FROM AUTHOR]

• A matrix approach allows resolving small-scale heterogeneities such as cracks and fractures with an optimal resolution. • Scatterers and reflectors are imaged at the SGB site on the San Jacinto Fault Zone, California. • Clear differences are found among different sections of the fault zone in the study area. Fault zones are associated with multi-scale heterogeneities of rock properties. Large scale variations may be imaged with conventional seismic reflection methods that detect offsets in geological units, and tomographic techniques that provide average seismic velocities in resolved volumes. However, characterizing elementary localized inhomogeneities of fault zones, such as cracks and fractures, constitutes a challenge for conventional techniques. Resolving these small-scale heterogeneities can provide detailed information for structural and mechanical models of fault zones. Recently, the reflection matrix approach utilizing body wave reflections in ambient noise cross-correlations was extended with the introduction of aberration corrections to handle the actual lateral velocity variations in the fault zone (Touma et al., 2021). Here this method is applied further to analyze the distribution of scatterers in the first few kilometers of the crust in the San Jacinto Fault Zone at the Sage Brush Flat (SGB) site, southeast of Anza, California. The matrix approach allows us to image not only specular reflectors but also to resolve the presence, location and reflectivity of scatterers for seismic waves starting with a simple homogeneous background velocity model of the medium. The derived three-dimensional image of the fault zone resolves lateral variations of scattering properties in the region within and around the surface fault traces, as well as differences between the Northwest (NW) and the Southeast (SE) parts of the study area. A localized intense damage zone at depth is observed in the SE section, suggesting that a geometrical complexity of the fault zone at depth induces ongoing generation of rock damage. [ABSTRACT FROM AUTHOR]