Search

We deployed a dense geodetic and seismological network in the Atacama seismic gap in Chile. We derive a microseismicity catalog of >30,000 events, time series from 70 GNSS stations, and utilize a transdimensional Bayesian inversion to estimate interplate locking. We identify two highly locked regions of different sizes whose geometries appear to control seismicity patterns. Interface seismicity concentrates beneath the coastline, just downdip of the highest locking. A region with lower locking (27.5°S–27.7°S) coincides with higher seismicity levels, a high number of repeating earthquakes and events extending toward the trench. This area is situated where the Copiapó Ridge is subducted and has shown previous indications of both seismic and aseismic slip, including an earthquake sequence in 2020. While these findings suggest that the structure of the downgoing oceanic plate prescribes patterns of interplate locking and seismicity, we note that the Taltal Ridge further north lacks a similar signature. Plain Language Summary: Deformation along plate boundaries can occur seismically (i.e. through earthquakes) as well as aseismically (i.e. slipping slowly), and it is important to understand where each of these modes is dominant. Along the Chilean subduction contact, North‐Central Chile is the only place where aseismic deformation episodes have been observed so far. In order to study these processes in detail, we deployed and operated dense geodetic and seismological networks in this region. Analyzing the data collected by these networks, we find notable relationships between seismic and aseismic processes. Thousands of small earthquakes are found at the boundaries of locked regions, whereas no small earthquakes are found at their interior. Thus, implying such regions are mechanically coupled, that is, currently accumulating elastic deformation energy that will 1 day be released during a large earthquake. Along the North‐Central Chilean plate boundary, there is one region (around 27.5°S) that shows many signs of aseismic deformation. It is located where a chain of seamounts is being subducted, which is likely responsible for the different behavior of this segment. Key Points: Microseismicity catalog and map of interplate locking derived for the Atacama 1922 seismic gap in North‐Central ChileSeismicity in vicinity of plate interface coincides with downdip edge of high couplingSeismo‐geodetic signals due to the subduction of the Copiapó ridge are prominent but negligible for the subducting Taltal Ridge [ABSTRACT FROM AUTHOR]

At the Blanco transform fault system (BTFS) off Oregon, 138 local earthquakes and 84 double‐couple focal mechanisms from ocean‐bottom‐seismometer recordings jointly discussed with bathymetric features reveal a highly segmented transform system without any prominent fracture zone traces longer than 100 km. In the west, seismicity is focused at deep troughs (i.e., the West and East Blanco, and Surveyor Depressions). In the east, the BTFS lacks a characteristic transform valley and instead developed the Blanco Ridge, which is the most seismically active feature, showing strike‐slip and dip‐slip faulting. Sandwiched between the two main segments of the BTFS is the Cascadia Depression, representing a short intra‐transform spreading segment. Seismic slip vectors reveal that stresses at the eastern BTFS are roughly in line with plate motion. In contrast, stresses to the west are clockwise skewed, indicating ongoing reorganization of the OTF system. As we observed no prominent fracture zones at the BTFS, plate tectonic reconstructions suggest that the BTFS developed from non‐transform offsets rather than pre‐existing transform faults during a series of ridge propagation events. Our observations suggest that the BTFS can be divided into two oceanic transform systems. The eastern BTFS is suggested to be a mature transform plate boundary since ∼0.6 Ma. In contrast, the western BTFS is an immature transform system, which is still evolving to accommodate far‐field stress change. The BTFS acts as a natural laboratory to yield processes governing the development of oceanic transform faults. Plain Language Summary: The Blanco transform fault system (BTFS) northwest off the coast of Oregon is seismically very active. We used 1 year of ocean bottom seismometer data collected between September 2012 and October 2013 to locate 138 local earthquakes. The events align perfectly with the morphologic features of the BTFS, dividing the BTFS into five transform segments and two short intra‐transform spreading centers. Furthermore, we observe different seismotectonic behaviors of the western and eastern BTFS based on the along‐strike variation in morphology, magnetization, focal depth distribution, and strain partitioning. Although many segmented oceanic transform systems were formed from a single transform fault in response to rotations in plate motion, the BTFS turns out to be originated from non‐transform offsets between ridge segments, as we observed no prominent fracture zone traces neither in morphology nor gravity field data. A clockwise shift in the Juan de Fuca/Pacific pole of rotation at ∼5 Ma followed by a series of ridge propagation events initiated the formation of the BTFS, integrated each segment of the BTFS by shortening the ridge segments in between. Our observations suggest that the Blanco Ridge and the Gorda transform segment in the eastern BTFS were formed at ∼1.6 and 0.6 Ma, respectively, and ever since, the eastern BTFS became a mature transform boundary. In contrast, seismic slip vectors comparing to plate motion directions reveal that stresses in the western BTFS are systematically skewed, suggesting the immature transform plate boundary is still adjusting to the new stress regime. Key Points: Local seismicity of the Blanco transform fault system (BTFS) reveals along‐strike variations dominated by strike‐slip and oblique dip‐slipThe BTFS developed from non‐transform offsets rather than discrete transform faults in response to plate rotation and ridge propagationThe BTFS consists of a mature plate boundary in the east and an immature system in the west, separated by a central spreading center [ABSTRACT FROM AUTHOR]

Seismic event detection and phase picking are the base of many seismological workflows. In recent years, several publications demonstrated that deep learning approaches significantly outperform classical approaches, achieving human‐like performance under certain circumstances. However, as studies differ in the datasets and evaluation tasks, it is unclear how the different approaches compare to each other. Furthermore, there are no systematic studies about model performance in cross‐domain scenarios, that is, when applied to data with different characteristics. Here, we address these questions by conducting a large‐scale benchmark. We compare six previously published deep learning models on eight data sets covering local to teleseismic distances and on three tasks: event detection, phase identification and onset time picking. Furthermore, we compare the results to a classical Baer‐Kradolfer picker. Overall, we observe the best performance for EQTransformer, GPD and PhaseNet, with a small advantage for EQTransformer on teleseismic data. Furthermore, we conduct a cross‐domain study, analyzing model performance on data sets they were not trained on. We show that trained models can be transferred between regions with only mild performance degradation, but models trained on regional data do not transfer well to teleseismic data. As deep learning for detection and picking is a rapidly evolving field, we ensured extensibility of our benchmark by building our code on standardized frameworks and making it openly accessible. This allows model developers to easily evaluate new models or performance on new data sets. Furthermore, we make all trained models available through the SeisBench framework, giving end‐users an easy way to apply these models. Plain Language Summary: The first step in many seismological workflows is identifying if a signal contains an earthquake, and at which time which type of seismic wave arrived. These steps are known as event detection, phase identification and phase picking. In recent years, machine learning methods, in particular deep learning methods have been developed, showing promising performance on these tasks. However, so far these models have not been compared systematically in a quantitative way. Here we evaluate the performance of six deep learning models on eight datasets. Additionally, we compare them to a traditional picking algorithm not using machine learning. From our results we identify that the models EQTransformer, GPD and PhaseNet perform best. As in many use cases no picker trained on the target region will be available, we further evaluated how well models are transferable across regions. We identified that transfer across regions works well as long as the distance ranges stay similar. To foster application of the results, we make all our trained models available through the SeisBench framework. Key Points: We conducted a large scale benchmark of machine learning pickers using six models and eight datasetsBest overall performance is observed for EQTransformer, GPD and PhaseNet, with advantages for EQTransformer on teleseismic distancesModels transfer well between different regions with similar distances, but not between regional and teleseismic distances [ABSTRACT FROM AUTHOR]

The aftershock distribution of the 2014 Mw 8.1 Iquique earthquake offshore northern Chile, identified from a long‐term deployment of ocean bottom seismometers installed eight months after the mainshock, in conjunction with seismic reflection imaging, provides insights into the processes regulating the updip limit of coseismic rupture propagation. Aftershocks updip of the mainshock hypocenter frequently occur in the upper plate and are associated with normal faults identified from seismic reflection data. We propose that aftershock seismicity near the plate boundary documents subduction erosion that removes mass from the base of the wedge and results in normal faulting in the upper plate. The combination of very little or no sediment accretion and subduction erosion over millions of years has resulted in a very weak and aseismic frontal wedge. Our observations thus link the shallow subduction zone seismicity to subduction erosion processes that control the evolution of the overriding plate. Plain Language Summary: To better understand the controls on shallow seismicity and subduction erosion following large subduction earthquakes, we use marine recordings of the Mw 8.1 2014 Iquique earthquake aftershocks and long‐offset multi‐channel seismic data. By comparing the aftershock locations and seismic imaging, we observe that most aftershocks occurred in the upper continental plate and abruptly stopped in the frontal forearc. The amplitude characteristics of upper‐crust reflections indicate a fractured and fluid‐filled outer forearc, which is associated with the absence of aftershocks. Large‐scale faulting, as evidenced by disrupted reflections in the seismic image, can be correlated to upper plate seismicity. We propose that the aftershocks updip of the main earthquake area reflect active subduction erosion processes. Key Points: We investigate structure and seismicity at the updip end of the 2014 Iquique earthquake rupture using amphibious seismic dataSeismicity updip of the 2014 Iquique earthquake occurs over a broad range likely interpreted to be related to the basal erosion processesCoseismic stress changes and aftershocks activate extensional faulting of the upper plate and subduction erosion [ABSTRACT FROM AUTHOR]

What controls subduction megathrust seismogenesis downdip of the mantle wedge corner (MWC)? We propose that, in the region of the 2010 Mw = 8.8 Maule, Chile, earthquake, serpentine minerals derived from the base of the hydrated mantle wedge exert a dominant control. Based on modeling, we predict that the megathrust fault zone near the MWC contains abundant lizardite/chrysotile‐rich serpentinite that transforms to antigorite‐rich serpentinite at greater depths. From the MWC at 32–40 km depth to at least 55 km, the predominantly velocity‐strengthening megathrust accommodated dynamic propagation of the 2010 rupture but with small slip and negative stress drop. The downdip distribution of interplate aftershocks exhibits a gap around the MWC that can be explained by the velocity‐strengthening behavior of lizardite/chrysotile. Interspersed velocity‐weakening and dynamic weakening antigorite‐rich patches farther downdip may be responsible for increased abundance of aftershocks and possibly for some of the high‐frequency energy radiation during the 2010 rupture. Plain Language Summary: A subduction megathrust rupture may extend from the trench to beneath populated coastal area. To understand what controls the megathrust seismogenic behavior of its deeper part, we conduct a case study of the 2010 magnitude 8.8 Maule, Chile, earthquake and its aftershocks. With numerical thermal modeling, we find that the deep seisomogenic behavior may be controlled by serpentine materials derived from the base of the overlying hydrated mantle wedge. The deep part of the megathrust ruptured "passively" against increasing resistance (negative stress drop). Lower‐temperature serpentines lizardite and chrysotile that are known to facilitate aseismic slip may be responsible for an observed gap in aftershock distribution in the downdip direction. Higher‐temperature serpentine antigorite that is known to facilitate seismic slip may be responsible for increased aftershocks farther downdip and possibly radiation of high‐frequency energy during the 2010 Maule earthquake. Key Points: Seismogenic behavior of subduction megathrust may be affected by materials derived from the base of the serpentinized forearc mantle wedgeThermopetrologic modeling suggests abundance of lizardite and chrysotile around the mantle wedge corner but antigorite farther downdipSerpentine frictional behavior may explain distributions of stress drop and seismic energy radiation in 2010 rupture and aftershocks [ABSTRACT FROM AUTHOR]

We examine the intra‐arc crustal seismicity of the Andean Southern Volcanic Zone. Our aim is to resolve interseismic deformation in an active magmatic arc dominated by both margin‐parallel (Liquiñe‐Ofqui fault system, LOFS) and Andean transverse faults. Crustal seismicity provides information about the schizosphere tectonic state, delineating the geometry and kinematics of high strain domains driven by oblique‐subduction. Here, we present local seismicity based on 16‐month data collected from 34 seismometers monitoring a ~200‐km‐long section of the Southern Volcanic Zone, including the Lonquimay and Villarrica volcanoes. We located 356 crustal events with magnitudes between Mw 0.6 and Mw 3.6. Local seismicity occurs at depths down to 40 km in the forearc and consistently shallower than 12 km beneath the volcanic chain, suggesting a convex shape of the crustal seismogenic layer bottom. Focal mechanisms indicate strike‐slip faulting consistent with ENE‐WSW shortening in line with the long‐term deformation history revealed by structural geology studies. However, we find regional to local‐scale variations in the shortening axes orientation as revealed by the nature and spatial distribution of microseismicity, within three distinctive latitudinal domains. In the northernmost domain, seismicity is consistent with splay faulting at the northern termination of the LOFS; in the central domain, seismicity distributes along ENE‐ and WNW‐striking discrete faults, spatially associated with, hitherto seismic Andean transverse faults. The southernmost domain, in turn, is characterized by activity focused along a N15°E striking master branch of the LOFS. These observations indicate a complex strain compartmentalization pattern within the intra‐arc crust, where variable strike‐slip faulting dominates over dip‐slip movements. Plain Language Summary: In active volcanic chains, there is a strong interplay between deformation and volcanism. In this research, we take the "pulse" of active tectonics in a volcanic arc setting by measuring natural seismicity. We installed a network of 34 seismometers over 200 km along the volcanic arc in Southcentral Chile. Our results indicate active faulting in coherence with long‐lived faults and the overall regional‐scale stress regime. In fact, crustal regions in which faults are more active are spatially associated with inherited faults and with higher temperature gradient domains, as inferred from volcanic and geothermal activity. The maximum depth of seismicity is shallow (<12 km) in the central part of the volcanic arc, whereas is much deeper (~40 km) toward the forearc and back‐arc regions. This observation suggests that elevated geothermal gradients lead to a thinner brittle portion of crust, which is then mechanically easier to (re‐)break up. Our results can be used for understanding the way by which faulting interacts with crustal fluids, which in turn may help improving geothermal exploration strategies and seismic hazard assessment. Furthermore, seismicity reveals detailed information on how ongoing deformation accommodates within complex fault systems characterized by different orientations and kinematics. Key Points: Crustal seismicity in the SVZ of the Andes reveals the nature of active upper crustal faulting consistent with long‐term arc tectonicsStrain partitioning is compartmentalized into arc‐parallel (LOFS) and Andean transverse faults (ATF)Crustal seismicity occurs at all depths down to 40 km depth in the forearc, but shallower than 12 km along the volcanic chain [ABSTRACT FROM AUTHOR]

Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (∼1 year), resulting in a small number of usable records of teleseismic earthquakes. Here we use OBSs deployed as midaperture array in the deep ocean (4.5-5.5 km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beamforming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8 km) of oceanic crust. Observations at single stations with thin sediments (300-400 m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of ∼70-80 km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at ∼410 km and ∼660 km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam-formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events. [ABSTRACT FROM AUTHOR]