ALBERT

Description: SUMMARY We present a technique to derive robust estimates for the crustal thickness and elastic properties, including anisotropy, from shear wave splitting of converted phases in receiver functions. We combine stacking procedures with a correction scheme for the splitting effect of the crustal converted Ps-phase and its first reverberation, the PpPs-phase, where we also allow for a predefined dipping Moho. The incorporation of two phases stabilizes the analysis procedure and allows to simultaneously solve for the crustal thickness, the ratio of average P- to S-wave velocities, the percentage of anisotropy and the fast-axis direction. The stacking is based on arrival times and polarizations computed using a ray-based algorithm. Synthetic tests show the robustness of the technique and its applicability to tectonic settings where dip of the Moho is significant. These tests also demonstrate that the effects of a dipping layer boundary may overprint a possible anisotropic signature. To constrain the uncertainty of our results we perform statistical tests based on a bootstrapping approach. We distinguish between different model classes by comparing the coherency of the stacked amplitudes after moveout correction. We apply the new technique to real-data examples from different tectonic regimes and show that coherency of the stacked receiver functions can be improved, when anisotropy and a dipping Moho are included in the analysis. The examples underline the advantages of statistical analyses when dealing with stacking procedures and potentially ambiguous solutions.

Description: We have constructed a 3-D shear wave velocity (Vs) model for the crust and uppermost mantle beneath the Middle East using Rayleigh wave records obtained from ambient-noise cross-correlations and regional earthquakes. We combined one decade of data collected from 852 permanent and temporary broad-band stations in the region to calculate group-velocity dispersion curves. A compilation of >54 000 ray paths provides reliable group-velocity measurements for periods between 2 and 150 s. Path-averaged group velocities calculated at different periods were inverted for 2-D group-velocity maps. To overcome the problem of heterogeneous ray coverage, we used an adaptive grid parametrization for the group-velocity tomographic inversion. We then sample the period-dependent group-velocity field at each cell of a predefined grid to generate 1-D group-velocity dispersion curves, which are subsequently inverted for 1-D Vs models beneath each cell and combined to approximate the 3-D Vs structure of the area. The Vs model shows low velocities at shallow depths (5–10 km) beneath the Mesopotamian foredeep, South Caspian Basin, eastern Mediterranean and the Black Sea, in coincidence with deep sedimentary basins. Shallow high-velocity anomalies are observed in regions such as the Arabian Shield, Anatolian Plateau and Central Iran, which are dominated by widespread magmatic exposures. In the 10–20 km depth range, we find evidence for a band of high velocities (>4.0 km s–1) along the southern Red Sea and Arabian Shield, indicating the presence of upper mantle rocks. Our 3-D velocity model exhibits high velocities in the depth range of 30–50 km beneath western Arabia, eastern Mediterranean, Central Iranian Block, South Caspian Basin and the Black Sea, possibly indicating a relatively thin crust. In contrast, the Zagros mountain range, the Sanandaj-Sirjan metamorphic zone in western central Iran, the easternmost Anatolian plateau and Lesser Caucasus are characterized by low velocities at these depths. Some of these anomalies may be related to thick crustal roots that support the high topography of these regions. In the upper mantle depth range, high-velocity anomalies are obtained beneath the Arabian Platform, southern Zagros, Persian Gulf and the eastern Mediterranean, in contrast to low velocities beneath the Red Sea, Arabian Shield, Afar depression, eastern Turkey and Lut Block in eastern Iran. Our Vs model may be used as a new reference crustal model for the Middle East in a broad range of future studies.

Description: Rodrigues Ridge connects the Réunion hotspot track with the Central Indian Ridge (CIR) and has been suggested to represent the surface expression of a sub-lithospheric flow channel. From global earthquake catalogues, the seismicity in the region has been associated mainly with events related to the fracture zones at the CIR. However, some segments of the CIR appear void of seismic events. Here, we report on the seismicity recorded at a temporary array of 10 seismic stations operating on Rodrigues Island from September 2014 to June 2016. The array analysis was performed in the time domain by time shifting and stacking the complete waveforms. Event distances were estimated based on a 1-D velocity model and the travel time differences between S and P wave arrivals. We detected and located 63 new events that were not reported by the global networks. Most of the events (51) are located off the CIR and can be classified as intraplate earthquakes. Local magnitudes varied between 1.6 and 3.7. Four seismic clusters were observed that occurred to the west of the spreading segment of the CIR. The Rodrigues Ridge appeared to be aseismic during the period of operation. The lack of seismic activity along both Rodrigues Ridge and the sections of the CIR to the east of Rodrigues may be explained by partially molten upper-mantle material, possibly in relation to the proposed material flow between the Réunion plume and the CIR.

Description: Constraining the architecture of complex 3D volcanic plumbing systems within active rifts, and their impact on rift processes, is critical for examining the interplay between faulting, magmatism and magmatic fluids in developing rift segments. The Natron basin of the East African Rift System provides an ideal location to study these processes, owing to its recent magmatic-tectonic activity and ongoing active carbonatite volcanism at Oldoinyo Lengai. Here, we report seismicity and fault plane solutions from a 10 month-long temporary seismic network spanning Oldoinyo Lengai, Naibor Soito volcanic field and Gelai volcano. We locate 6,827 earthquakes with ML −0.85 to 3.6, which are related to previous and ongoing magmatic and volcanic activity in the region, as well as regional tectonic extension. We observe seismicity down to ∼17 km depth north and south of Oldoinyo Lengai and shallow seismicity (3–10 km) beneath Gelai, including two swarms. The deepest seismicity (∼down to 20 km) occurs above a previously imaged magma body below Naibor Soito. These seismicity patterns reveal a detailed image of a complex volcanic plumbing system, supporting potential lateral and vertical connections between shallow- and deep-seated magmas, where fluid and melt transport to the surface is facilitated by intrusion of dikes and sills. Focal mechanisms vary spatially. T-axis trends reveal dominantly WNW-ESE extension near Gelai, while strike-slip mechanisms and a radial trend in P-axes are observed in the vicinity of Oldoinyo Lengai. These data support local variations in the state of stress, resulting from a combination of volcanic edifice loading and magma-driven stress changes imposed on a regional extensional stress field. Our results indicate that the southern Natron basin is a segmented rift system, in which fluids preferentially percolate vertically and laterally in a region where strain transfers from a border fault to a developing magmatic rift segment.

Description: Previous investigation of seismic anisotropy indicates the presence of a simple mantle flow regime beneath the Turkish-Anatolian Plateau and Arabian Plate. Numerical modeling suggests that this simple flow is a component of a large-scale global mantle flow associated with the African superplume, which plays a key role in the geodynamic framework of the Arabia-Eurasia continental collision zone. However, the extent and impact of the flow pattern farther east beneath the Iranian Plateau and Zagros remains unclear. While the relatively smoothly varying lithospheric thickness beneath the Anatolian Plateau and Arabian Plate allows progress of the simple mantle flow, the variable lithospheric thickness across the Iranian Plateau is expected to impose additional boundary conditions on the mantle flow field. In this study, for the first time, we use an unprecedented data set of seismic waveforms from a network of 245 seismic stations to examine the mantle flow pattern and lithospheric deformation over the entire region of the Iranian Plateau and Zagros by investigation of seismic anisotropy. We also examine the correlation between the pattern of seismic anisotropy, plate motion using GPS velocities and surface strain fields. Our study reveals a complex pattern of seismic anisotropy that implies a similarly complex mantle flow field. The pattern of seismic anisotropy suggests that the regional simple mantle flow beneath the Arabian Platform and eastern Turkey deflects as a circular flow around the thick Zagros lithosphere. This circular flow merges into a toroidal component beneath the NW Zagros that is likely an indicator of a lateral discontinuity in the lithosphere. Our examination also suggests that the main lithospheric deformation in the Zagros occurs as an axial shortening across the belt, whereas in the eastern Alborz and Kopeh-Dagh a belt-parallel horizontal lithospheric deformation plays a major role.

Description: During the first two days of August 2016 a seismic crisis occurred on Brava, Cabo Verde, which – according to observations based on a local seismic network – was characterized by more than a thousand volcano-seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic-array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes in the region of Fogo and Brava; 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a significant variation in seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.

Description: The Alpine orogeny is characterized by tectonic sequences of subduction and collision accompanied by break-off events and possibly preceded by a flip of subduction polarity. The tectonic evolution of the transition to the Eastern Alps has thus been under debate. The dense SWATH-D seismic network as a complementary experiment to the AlpArray seismic network provides unprecedented lateral resolution to address this ongoing discussion. We analyze the shear-wave splitting of this data set including stations of the AlpArray backbone in the region to obtain new insights into the deformation at depth from seismic anisotropy. Previous studies indicate two-layer anisotropy in the Eastern Alps. This is supported by the azimuthal pattern of the measured fast axis direction across all analyzed stations. However, the temporary character of the deployment requires a joint analysis of multiple stations to increase the number of events adding complementary information of the anisotropic properties of the mantle. We, therefore, perform a cluster analysis based on a correlation of energy tensors between all stations. The energy tensors are assembled from the remaining transverse energy after the trial correction of the splitting effect from two consecutive anisotropic layers. This leads to two main groups of different two-layer properties, separated approximately at 13°E. We identify a layer with a constant fast axis direction (measured clockwise with respect to north) of about 60° over the whole area, with a possible dip from west to east. The lower layer in the west shows N–S fast direction and the upper layer in the east shows a fast axis of about 115°. We propose two likely scenarios, both accompanied by a slab break-off in the eastern part. The continuous layer can either be interpreted as frozen-in anisotropy with a lithospheric origin or as an asthenospheric flow evading the retreat of the European slab that would precede the break-off event. In both scenarios, the upper layer in the east is a result of a flow through the gap formed in the slab break-off. The N–S direction can be interpreted as an asthenospheric flow driven by the retreating European slab but might also result from a deep-reaching fault-related anisotropy.

Description: SUMMARY To constrain seismic anisotropy under and around the Alps in Europe, we study SKS shear wave splitting from the region densely covered by the AlpArray seismic network. We apply a technique based on measuring the splitting intensity, constraining well both the fast orientation and the splitting delay. Four years of teleseismic earthquake data were processed, from 723 temporary and permanent broad-band stations of the AlpArray deployment including ocean-bottom seismometers, providing a spatial coverage that is unprecedented. The technique is applied automatically (without human intervention), and it thus provides a reproducible image of anisotropic structure in and around the Alpine region. As in earlier studies, we observe a coherent rotation of fast axes in the western part of the Alpine chain, and a region of homogeneous fast orientation in the Central Alps. The spatial variation of splitting delay times is particularly interesting though. On one hand, there is a clear positive correlation with Alpine topography, suggesting that part of the seismic anisotropy (deformation) is caused by the Alpine orogeny. On the other hand, anisotropic strength around the mountain chain shows a distinct contrast between the Western and Eastern Alps. This difference is best explained by the more active mantle flow around the Western Alps. The new observational constraints, especially the splitting delay, provide new information on Alpine geodynamics.