Open Access

There is a long-standing debate over how many and what types of plumes underlie the East African Rift and whether they do or do not drive its extension and consequent magmatism and seismicity. Here we present a new tomographic study of relative teleseismic S and SKS residuals that expands the resolution from previous regional studies below the northern East African Rift to image structure from the surface to the base of the transition zone. The images reveal two low-velocity clusters, below Afar and west of the Main Ethiopian Rift, that extend throughout the upper mantle and comprise several smaller-scale (about 100 km diameter), low-velocity features. These structures support those of our recent P tomographic study below the region. The relative magnitude of S to P residuals is around 3.5, which is consistent with a predominantly thermal nature of the anomalies. The S and P velocity anomalies in the low-velocity clusters can be explained by similar excess temperatures in the range of 100–200°C, consistent with temperatures inferred from other seismic, geochemical, and petrological studies. Somewhat stronger VS anomalies below Afar than west of the Main Ethiopian Rift may include an expression of volatiles and/or melt in this region. These results, together with a comparison with previous larger-scale tomographic models, indicate that these structures are likely small-scale upwellings with mild excess temperatures, rising from a regional thermal boundary layer at the base of the upper mantle.

During the first two days of August 2016 a seismic crisis occurred on Brava, Cape Verde, which – according to observations based on a local seismic network – was characterized by more than 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, 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 migration of 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.

We have analysed the microseismic activity within the Rwenzori Mountains area in the western branch of the East African Rift. Seismogram recordings from a temporary array of up to 27 stations reveal approximately 800 events per month with local magnitudes ranging from –0.5 to 5.1. The earthquake distribution is highly heterogeneous. The majority of located events lie within faults zones to the east and west of the Rwenzoris with the highest seismic activity observed in the northeastern area, where the mountains are in contact with the rift shoulders. The hypocentral depth distribution exhibits a pronounced peak of seismic energy release at 15 km depth. The maximum extent of seismicity ranges from 20 to 32 km and correlates well with Moho depths that were derived from teleseismic receiver functions. We observe two general features: (i) beneath the rift shoulders, seismicity extends from the surface down to ca. 30 km depth; (ii) beneath the rift valley, seismicity is confined to depths greater than 10 km. From the observations there is no indication for a crustal root beneath the Rwenzori Mountains. The magnitude frequency distribution reveals a b-value of 1.1, which is consistent with the hypothesis that part of the seismicity is caused by magmatic processes within the crust. Fault plane solutions of 304 events were derived from P-polarities and SV/P amplitude ratios. More than 70 % of the source mechanisms exhibit pure or predominantly normal faulting. T-axis trends are highly uniform and oriented WNW–ESE, which is perpendicular to the rift axis and in good agreement with kinematic rift models. At the northernmost part of the region we observe a rotation of the T-axis trends to NEN–SWS, which may be indicative of a local perturbation of the regional stress field.

Seismic arrays provide tools for the localization of events without clear phases or events outside the network, where the station coverage prohibits classical localization techniques. Beam forming allows the determination of the direction (back azimuth) and horizontal (apparent) velocity of an incoming wavefront. Here we combine multiple arrays to retrieve event epicentres from the area of intersecting beams without the need to specify a velocity model. The analysis is performed in the time domain, which allows selecting a relatively narrow time window around the phase of interest while preserving frequency bandwidth. This technique is applied to earthquakes and hybrid events in the region of Fogo and Brava, two islands of the southern chain of the Cape Verde archipelago. The results show that the earthquakes mainly originate near Brava, whereas the hybrid events are located on Fogo. By multiple-event beam stacking we are able to further constrain the epicentral locations of the hybrid events in the northwestern part of the collapse scar of Fogo. In previous studies, these events were attributed to shallow hydrothermal processes. However, we obtain relatively high apparent velocities at the arrays, pointing to either deeper sources or complex ray paths. For a better understanding of possible errors of the multi-array analysis, we also compare slowness values obtained from the array analysis with those derived from earthquake locations from classical (local network) localizations. In general, the results agree well. Nevertheless, some systematic deviations of the array-derived back-azimuth and slowness values occur that can be quantified for certain event locations.

In this study, we determine spectral characteristics and amplitude decays of wind turbine induced seismic signals in the far field of a wind farm (WF) close to Uettingen, Germany. Average power spectral densities (PSDs) are calculated from 10 min time segments extracted from (up to) 6 months of continuous recordings at 19 seismic stations, positioned along an 8 km profile starting from the WF. We identify seven distinct PSD peaks in the frequency range between 1 and 8 Hz that can be observed to at least 4 km distance; lower-frequency peaks are detectable up to the end of the profile. At distances between 300 m and 4 km the PSD amplitude decay can be described by a power law with exponent b. The measured b values exhibit a linear frequency dependence and range from b=0.39 at 1.14 Hz to b=3.93 at 7.6 Hz. In a second step, the seismic radiation and amplitude decays are modeled using an analytical approach that approximates the surface wave field. Since we observe temporally varying phase differences between seismograms recorded directly at the base of the individual wind turbines (WTs), source signal phase information is included in the modeling approach. We show that phase differences between source signals have significant effects on the seismic radiation pattern and amplitude decays. Therefore, we develop a phase shift elimination method to handle the challenge of choosing representative source characteristics as an input for the modeling. To optimize the fitting of modeled and observed amplitude decay curves, we perform a grid search to constrain the two model parameters, i.e., the seismic shear wave velocity and quality factor. The comparison of modeled and observed amplitude decays for the seven prominent frequencies shows very good agreement and allows the constraint of shear velocities and quality factors for a two-layer model of the subsurface. The approach is generalized to predict amplitude decays and radiation patterns for WFs of arbitrary geometry.