ALBERT

Description: 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 mid-aperture 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 beam forming 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.

Description: At the eastern end of the Azores-Gloria transform fault system to the southwest of Portugal, the plate boundary between Africa and Iberia is a region where deformation is accommodated over a wide tectonically-active area. The region has unleashed large earthquakes and tsunamis, including the Mw ~ 8.5 Great Lisbon earthquake of 1755. Although the source region of the 1755 earthquake is still disputed, most proposals include a source location in the vicinity of the Horseshoe Abyssal Plain (HAP), which is bounded by the 5000 m high Gorringe Bank (GB). In this study we characterise seismic activity in the region using data recorded by two local networks of ocean-bottom seismometers (OBS). The networks were deployed in the eastern HAP and at the GB. The dataset allowed the detection of 160 local earthquakes. These earthquakes cluster around the GB, to the SW of Cabo Sao Vicente, and in the HAP. Focal depths indicate deep-seated earthquakes, with depths increasing from 20-35 km (mean of 26.1 ± 7.2 km) at the GB to 15-45 km (mean 31.5 km ± 10.5 km) under the HAP. Seismic activity thus extends down to levels that are deeper than those mapped by active seismic profiling, with the majority of events occurring within the mantle. Thermal modelling suggests that temperatures of approximately 600 °C characterise the base of the seismogenic brittle lithosphere at ~45 km depth. The large source depth and thermal structure supports previous suggestions that catastrophic seismic rupture through the lithospheric mantle may indeed occur in the area.

Description: We present a new tomographic model of the mantle in the area of the 2010 M8.8 Maule earthquake and surrounding regions. Increased ray coverage provided by the aftershock data allows us to image the detailed subducting slab structure in the mantle, from the region of flat slab subduction north of the Maule rupture to the area of overlapping rupture between the 1960 M9.5 and the 2010 M8.8 events to the south. We have combined teleseismic primary and depth phase arrivals with available local arrivals to better constrain the teleseismic earthquake locations in the region, which we use to conduct nested regional–global tomography. The new model reveals the detailed structure of the flat slab and its transition to a more moderately dipping slab in the Maule region. South of the Maule region, a steeply dipping relic slab is imaged from ∼200 to 1000 km depth that is distinct from the moderately dipping slab above it and from the more northerly slab at similar depths. We interpret the images as revealing both horizontal and vertical tearing of the slab at ∼38°S to explain the imaged pattern of slab anomalies in the southern portion of the model. In contrast, the transition from a horizontal to moderately subducting slab in the northern portion of the model is imaged as a continuous slab bend. We speculate that the tearing was most likely facilitated by a fracture zone in the downgoing plate or alternatively by a continental scale terrane boundary in the overriding plate.

Description: The submarine Istanbul-Silivri fault segment, within 15km of Istanbul, is the only portion of the North Anatolian Fault that has not ruptured in the last 250years. We report first results of a seafloor acoustic ranging experiment to quantify current horizontal deformation along this segment and assess whether the segment is creeping aseismically or accumulating stress to be released in a future event. Ten transponders were installed to monitor length variations along 15 baselines. A joint least squares inversion for across-fault baseline changes, accounting for sound speed drift at each transponder, precludes fault displacement rates larger than a few millimeters per year during the 6month observation period. Forward modeling shows that the data better fit a locked state or a very moderate surface creep-less than 6mm/yr compared to a far-field slip rate of over 20mm/yr-suggesting that the fault segment is currently accumulating stress.

Description: Our knowledge of the absolute S wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P wave polarization (apparent P wave incidence angle) of teleseismic events to investigate the S wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent P wave incidence angle at the ocean bottom dependent on the half space S wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBS) to derive apparent S wave velocity profiles. These profiles are dependent on the S wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S wave velocity-depth models by a three step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S wave velocities. The apparent S wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 km and 0.9 km with S wave velocities between 0.7 km s−1 and 1.4 km s−1. The estimated total crustal thickness varies between 4 km and 10 km with S wave velocities between 3.5 km s−1 and 4.3 km s−1. We find a slight increase of the total crustal thickness from ∼5 km to ∼8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S wave velocities decrease from values around 5.5 km s−1 to 4.5 km s−1 towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S wave velocities for real and synthetic data to estimate S wave velocity-depth models.

Description: Seafloor seepage of hydrocarbon-bearing fluids has been identified in a number of marine forearcs. However, temporal variations in seep activity and the structural and tectonic parameters that control the seepage often remain poorly constrained. Subduction-zone earthquakes for example, are often discussed to trigger seafloor seepage but causal links that go beyond theoretical considerations have not yet been fully established. This is mainly due to the inaccessibility of offshore epicentral areas, the infrequent occurrence of large earthquakes, and challenges associated with offshore monitoring of seepage over large areas and sufficient time periods. Here, we report visual, geochemical, geophysical, and modelling results and observations from the Concepción Methane Seep Area (offshore Central Chile) located in the rupture area of the 2010 Mw. 8.8 Maule earthquake. High methane concentrations in the oceanic water column and a shallow sub-bottom depth of sulfate penetration indicate active methane seepage. The stable carbon isotope signature of the methane and hydrocarbon composition of the released gas indicate a mixture of shallow-sourced biogenic gas and a deeper sourced thermogenic component. Pristine fissures and fractures observed at the seafloor together with seismically imaged large faults in the marine forearc may represent effective pathways for methane migration. Upper-plate fault activity with hydraulic fracturing and dilation is in line with increased normal Coulomb stress during large plate-boundary earthquakes, as exemplarily modelled for the 2010 earthquake. On a global perspective our results point out the possible role of recurring large subduction-zone earthquakes in driving hydrocarbon seepage from marine forearcs over long timescales.

Description: Powerful subduction zone earthquakes rupture thousands of square kilometers along continental margins but at certain locations earthquake rupture terminates. To date detailed knowledge of the parameters that govern seismic rupture and aftershocks is still incomplete. On 16 September 2015 the Mw. 8.3 Illapel earthquake ruptured a 200 km long stretch of the Central Chilean subduction zone, triggering a tsunami and causing significant damage. Here we analyze the temporal and spatial pattern of the co-seismic rupture and aftershocks in relation to the tectonic setting in the earthquake area. Aftershocks cluster around the area of maximum coseismic slip, in particular in lateral and downdip direction. During the first 24 hours after the mainshock, aftershocks migrated in both lateral directions with velocities of approximately 2.5 and 5 km/h. At the southern rupture boundary aftershocks cluster around individual subducted seamounts that are related to the downthrusting Juan Fernández Ridge. In the northern part of the rupture area aftershocks separate into an upper cluster (above 25 km depth) and a lower cluster (below 35 km depth). This dual seismic-aseismic transition in downdip direction is also observed in the interseismic period suggesting that it may represent a persistent feature for the Central Chilean subduction zone.

Description: The postseismic deformation captured with continuous Global Positioning System (cGPS) monitoring following many recent mega-thrust events has been shown to be a signal composed of two dominant processes: afterslip on the plate interface and viscoelastic relaxation of the continental and oceanic mantles in response to the coseismic stress perturbation. Following the south-central Chile 2010 Maule Mw 8.8 earthquake, the time series from the regional cGPS network show a distinct curvature in the pathway of the horizontal motion that is not easily fit by a stationary decaying pattern of afterslip in combination with viscoelastic relaxation. Here we show that with realistic assumptions about the long-term decay of the afterslip signal, the postseismic signal can be decomposed into three first-order contributing processes: plate interface re-locking, plate interface afterslip, and mantle viscoelastic relaxation. From our analyses we conclude that the plate interface recovers its interseismic locking state rapidly (model space ranges between an instant recovery and a period of 1 year); a finding that supports laboratory experimental evidence as well as some recent studies of aftershocks and postseismic surface deformation. Furthermore, re-locking is the main cause of the curvature in the cGPS signal, and this study presents a plausible range of geodetic re-locking rates following a megathrust earthquake.

Description: An important tool for understanding deformation occurring within a subduction zone is the measurement of seismic anisotropy through observations of shear wave splitting (SWS). In Sumatra, two temporary seismic networks were deployed between December 2007 and February 2009, covering the fore arc between the fore-arc islands to the back arc. We use SKS and local SWS measurements to determine the type, amount, and location of anisotropy. Local SWS measurements from the fore-arc islands exhibit trench-parallel fast directions which can be attributed to shape preferred orientation of cracks/fractures in the overriding sediments. In the Sumatran Fault region, the predominant fast direction is fault/trench parallel, while in the back-arc region it is trench perpendicular. The trench-perpendicular measurements exhibit a positive correlation between delay time and raypath length in the mantle wedge, while the fault-parallel measurements are similar to the fault-parallel fast directions observed for two crustal events at the Sumatran Fault. This suggests that there are two layers of anisotropy: one due to entrained flow within the mantle wedge and a second layer within the overriding crust due to the shear strain caused by the Sumatran Fault. SKS splitting results show a NNW-SSE fast direction with delay times of 0.8–3.0 s. The fast directions are approximately parallel to the absolute plate motion of the subducting Indo-Australian Plate. The small delay times exhibited by the local SWS (0.05–0.45 s), in combination with the large SKS delay times, suggest that the anisotropy generating the teleseismic SWS is dominated by entrained flow in the asthenosphere below the slab.

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.