ALBERT

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel-time data of local earthquakes for vp and vp∕vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean-bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). The 2-D and 3-D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3-D velocity model shows prevailing trench-parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at  ∼ 180 km distance to the trench. vp velocities beneath the magmatic arc and the Sumatran fault zone (SFZ) are around 5 km s−1 at 10 km depth and the vp∕vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp∕vs values of 1.85 at  ∼ 150 km distance to the trench in the region of the Mentawai Fault. vp∕vs ratios suggest an absence of large-scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between the Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: We present results of a seismic refraction experiment which determines the crustal and upper-mantle structure of an oceanic core complex (OCC) and its conjugate side located south of the 5°S ridge–transform intersection at the Mid-Atlantic Ridge. The core complex with a corrugated surface has been split by a change in location of active seafloor spreading, resulting in two massifs on either side of the current spreading axis. We applied a joint tomographic inversion of wide-angle reflected and refracted phases for five intersecting seismic profiles. The obtained velocity models are used to constrain the magmatic evolution of the core complex from the analysis of seismic layer 3 and crustal thickness. An abrupt increase of crustal velocities at shallow depth coincides with the onset of the seafloor corrugations at the exposed footwall. The observed velocity structure is consistent with the presence of gabbros directly beneath the corrugated fault surface. The thickness of the high-velocity body is constrained by PmP reflections to vary along and across axis between 〈3 and 5 km. The thickest crust is associated with the central phase of detachment faulting at the higher-elevated northern portion of the massif. Beneath the breakaway of the OCC the crust is 2.5 km thick and reveals significantly lower velocities. This implies that the fault initially exhumed low-velocity material overlying the gabbro plutons. In contrast, crust formed at the conjugate side during OCC formation is characterized by an up to 2-km-thick seismic layer 2 overlying a 1.7-km-thick seismic layer 3. Obtained upper-mantle velocities range from 7.3 to 7.9 km s−1 and seem to increase with distance from the median valley. However, velocities of 7.3–7.5 km s−1 beneath the older portions of the OCC may derive from deep fluid circulation and related hydrothermal alteration, which may likely be facilitated by the subsequent rifting. Our velocity models reveal a strongly asymmetric velocity structure across the ridge axis, associated with the accretion of gabbros into the footwall of the detachment fault and upper-crustal portions concentrated at the conjugate side. Our results do not support a substantial increase in the axial ridge's melt supply related to the final phase of detachment faulting. Hence, the footwall rifting at 5°S may be a generic mechanism of detachment termination under very low melt conditions, as predicted by recent numerical models of Tucholke et al.

Description: We present a detailed 3-D P-wave velocity model obtained by first-arrival travel-time tomography with seismic refraction data in the segment boundary of the Sumatra subduction zone across Simeulue Island, and an image of the top of the subducted oceanic crust extracted from depth-migrated multi-channel seismic reflection profiles. We have picked P-wave first arrivals of the air-gun source seismic data recorded by local networks of ocean-bottom seismometers, and inverted the travel-times for a 3-D velocity model of the subduction zone. This velocity model shows an anomalous zone of intermediate velocities between those of oceanic crust and mantle that is associated with raised topography on the top of the oceanic crust. We interpret this feature as a thickened crustal zone in the subducting plate with compositional and topographic variations, providing a primary control on the upper plate structure and on the segmentation of the 2004 and 2005 earthquake ruptures.

Description: We present a new 3-D shear-velocity model for the top 30 km of the crust in the wider Vienna Basin region based on surface waves extracted from ambient-noise cross-correlations. We use continuous seismic records of 63 broad-band stations of the AlpArray project to retrieve interstation Green’s functions from ambient-noise cross-correlations in the period range from 5 to 25 s. From these Green’s functions, we measure Rayleigh group traveltimes, utilizing all four components of the cross-correlation tensor, which are associated with Rayleigh waves (ZZ, RR, RZ and ZR), to exploit multiple measurements per station pair. A set of selection criteria is applied to ensure that we use high-quality recordings of fundamental Rayleigh modes. We regionalize the interstation group velocities in a 5 km × 5 km grid with an average path density of ∼20 paths per cell. From the resulting group-velocity maps, we extract local 1-D dispersion curves for each cell and invert all cells independently to retrieve the crustal shear-velocity structure of the study area. The resulting model provides a previously unachieved lateral resolution of seismic velocities in the region of ∼15 km. As major features, we image the Vienna Basin and Little Hungarian Plain as low-velocity anomalies, and the Bohemian Massif with high velocities. The edges of these features are marked with prominent velocity contrasts correlated with faults, such as the Alpine Front and Vienna Basin transfer fault system. The observed structures correlate well with surface geology, gravitational anomalies and the few known crystalline basement depths from boreholes. For depths larger than those reached by boreholes, the new model allows new insight into the complex structure of the Vienna Basin and surrounding areas, including deep low-velocity zones, which we image with previously unachieved detail. This model may be used in the future to interpret the deeper structures and tectonic evolution of the wider Vienna Basin region, evaluate natural resources, model wave propagation and improve earthquake locations, among others.

Description: The dense AlpArray network allows studying seismic wave propagation with high spatial resolution. Here we introduce an array approach to measure arrival angles of teleseismic Rayleigh waves. The approach combines the advantages of phase correlation as in the two-station method with array beamforming to obtain the phase-velocity vector. 20 earthquakes from the first two years of the AlpArray project are selected, and spatial patterns of arrival-angle deviations across the AlpArray are shown in maps, depending on period and earthquake location. The cause of these intriguing spatial patterns is discussed. A simple wave-propagation modelling example using an isolated anomaly and a Gaussian beam solution suggests that much of the complexity can be explained as a result of wave interference after passing a structural anomaly along the wave paths. This indicates that arrival-angle information constitutes useful additional information on the Earth structure, beyond what is currently used in inversions.

Description: 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.

Description: 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 and even achieve human-like performance under certain circumstances. However, as most studies differ in the datasets and exact evaluation tasks studied, it is yet unclear how the different approaches compare to each other. Furthermore, there are no systematic studies how the models perform in a cross-domain scenario, i.e., when applied to data with different characteristics. Here, we address these questions by conducting a large-scale benchmark study. We compare six previously published deep learning models on eight datasets 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 EQTransformer having a small advantage for teleseismic data. Furthermore, we conduct a cross-domain study, in which we analyze model performance on datasets they were not trained on. We show that trained models can be transferred between regions with only mild performance degradation, but not from regional to teleseismic data or vice versa. 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 compare new models or evaluate performance on new datasets, beyond those presented here. Furthermore, we make all trained models available through the SeisBench framework, giving end-users an easy way to apply these models in seismological analysis.

Description: 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 apply a transdimensional Bayesian inversion to estimate interplate locking degree. 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 of lower interplate locking around 27.5ºS coincides with higher seismicity levels, a high number of repeating earthquakes and events extending further towards the trench. Having shown numerous signs of aseismic deformation (slow-slip events and earthquake swarms), this area is situated where the Copiapó Ridge is subducted. 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.