Rift-induced disruption of cratonic keels drives kimberlite volcanism.

Kimberlites are volatile-rich, occasionally diamond-bearing magmas that have erupted explosively at Earth's surface in the geologic past1-3. These enigmatic magmas, originating from depths exceeding 150 km in Earth's mantle1, occur in stable cratons and in pulses broadly synchronous with supercontinent cyclicity4. Whether their mobilization is driven by mantle plumes5 or by mechanical weakening of cratonic lithosphere4,6 remains unclear. Here we show that most kimberlites spanning the past billion years erupted about 30 million years (Myr) after continental breakup, suggesting an association with rifting processes. Our dynamical and analytical models show that physically steep lithosphere-asthenosphere boundaries (LABs) formed during rifting generate convective instabilities in the asthenosphere that slowly migrate many hundreds to thousands of kilometres inboard of rift zones. These instabilities endure many tens of millions of years after continental breakup and destabilize the basal tens of kilometres of the cratonic lithosphere, or keel. Displaced keel is replaced by a hot, upwelling mixture of asthenosphere and recycled volatile-rich keel in the return flow, causing decompressional partial melting. Our calculations show that this process can generate small-volume, low-degree, volatile-rich melts, closely matching the characteristics expected of kimberlites1-3. Together, these results provide a quantitative and mechanistic link between kimberlite episodicity and supercontinent cycles through progressive disruption of cratonic keels.

Spatial forecasting of seismicity provided from Earth observation by space satellite technology.

Understanding the controls on the distribution and magnitude of earthquakes is required for effective earthquake forecasting. We present a study that demonstrates that the distribution and size of earthquakes in Italy correlates with the steady state rate at which the Earth's crust moves. We use a new high-resolution horizontal strain rate (S) field determined from a very dense velocity field derived from the combination of Global Navigation Satellite System (GNSS) and satellite radar interferometry from two decades of observations. Through a statistical approach we study the correlation between the S and the magnitude of M ≥ 2.5 earthquakes that occurred in the same period of satellite observations. We found that the probability of earthquakes occurring is linked to S by a linear correlation, and more specifically the probability that a strong seismic event occurs doubles with the doubling of S. It also means that lower horizontal strain rate zone can have as large earthquakes as high horizontal strain rate zones, just with a reduced probability. The work demonstrates an independent and quantitative tool to spatially forecast seismicity.

3D Extension at Plate Boundaries Accommodated by Interacting Fault Systems.

Complex patterns of normal faults with multiple orientations and/or highly curved shapes have been traditionally explained by successive tectonic phases of 2-dimensional deformation. Alternatively, multiple fault sets have been proposed to develop simultaneously and in orthorhombic symmetry during a single phase of 3-dimensional deformation. We use analogue models of normal faults to demonstrate that, without the influence of pre-existing structures, 3D extension is preferentially accommodated by the alternate, rather than simultaneous, development of faults with different trends. By means of stress-driven interactions, 3D deformation can be partitioned into coupled systems of normal faults, which display geometries commonly observed in tectonic settings affected by interacting plate boundaries. Under radial extension, deformation is accommodated by major curvilinear grabens coupled with minor perpendicular faults, resulting in the triple junctions of grabens observed in Afar. On the other hand, the alternate development of perpendicular faults accommodates synchronous bi-directional and mutually perpendicular extension, giving the same fault pattern observed in the Barents Sea rift-shear margin.

Aborted propagation of the Ethiopian rift caused by linkage with the Kenyan rift.

Continental rift systems form by propagation of isolated rift segments that interact, and eventually evolve into continuous zones of deformation. This process impacts many aspects of rifting including rift morphology at breakup, and eventual ocean-ridge segmentation. Yet, rift segment growth and interaction remain enigmatic. Here we present geological data from the poorly documented Ririba rift (South Ethiopia) that reveals how two major sectors of the East African rift, the Kenyan and Ethiopian rifts, interact. We show that the Ririba rift formed from the southward propagation of the Ethiopian rift during the Pliocene but this propagation was short-lived and aborted close to the Pliocene-Pleistocene boundary. Seismicity data support the abandonment of laterally offset, overlapping tips of the Ethiopian and Kenyan rifts. Integration with new numerical models indicates that rift abandonment resulted from progressive focusing of the tectonic and magmatic activity into an oblique, throughgoing rift zone of near pure extension directly connecting the rift sectors.

Imaging Lithospheric Discontinuities Beneath the Northern East African Rift Using S-to-P Receiver Functions.

Imaging the lithosphere is key to understand mechanisms of extension as rifting progresses. Continental rifting results in a combination of mechanical stretching and thinning of the lithosphere, decompression upwelling, heating, sometimes partial melting of the asthenosphere, and potentially partial melting of the mantle lithosphere. The northern East African Rift system is an ideal locale to study these processes as it exposes the transition from tectonically active continental rifting to incipient seafloor spreading. Here we use S-to-P receiver functions to image the lithospheric structure beneath the northernmost East African Rift system where it forms a triple junction between the Main Ethiopian rift, the Red Sea rift, and the Gulf of Aden rift. We image the Moho at 31 ± 6 km beneath the Ethiopian plateau. The crust is 28 ± 3 km thick beneath the Main Ethiopian rift and thins to 23 ± 2 km in northern Afar. We identify a negative phase, a velocity decrease with depth, at 67 ± 3 km depth beneath the Ethiopian plateau, likely associated with the lithosphere-asthenosphere boundary (LAB), and a lack of a LAB phase beneath the rift. Using observations and waveform modeling, we show that the LAB phase beneath the plateau is likely defined by a small amount of partial melt. The lack of a LAB phase beneath the rift suggests melt percolation through the base of the lithosphere beneath the northernmost East African Rift system.

The initiation of segmented buoyancy-driven melting during continental breakup.

Melting of the mantle during continental breakup leads to magmatic intrusion and volcanism, yet our understanding of the location and dominant mechanisms of melt generation in rifting environments is impeded by a paucity of direct observations of mantle melting. It is unclear when during the rifting process the segmented nature of magma supply typical of seafloor spreading initiates. Here, we use Rayleigh-wave tomography to construct a high-resolution absolute three-dimensional shear-wave velocity model of the upper 250 km beneath the Afar triple junction, imaging the mantle response during progressive continental breakup. Our model suggests melt production is highest and melting depths deepest early during continental breakup. Elevated melt production during continental rifting is likely due to localized thinning and melt focusing when the rift is narrow. In addition, we interpret segmented zones of melt supply beneath the rift, suggesting that buoyancy-driven active upwelling of the mantle initiates early during continental rifting.

First recorded eruption of Nabro volcano, Eritrea, 2011.

We present a synthesis of diverse observations of the first recorded eruption of Nabro volcano, Eritrea, which began on 12 June 2011. While no monitoring of the volcano was in effect at the time, it has been possible to reconstruct the nature and evolution of the eruption through analysis of regional seismological and infrasound data and satellite remote sensing data, supplemented by petrological analysis of erupted products and brief field surveys. The event is notable for the comparative rarity of recorded historical eruptions in the region and of caldera systems in general, for the prodigious quantity of SO2 emitted into the atmosphere and the significant human impacts that ensued notwithstanding the low population density of the Afar region. It is also relevant in understanding the broader magmatic and tectonic significance of the volcanic massif of which Nabro forms a part and which strikes obliquely to the principal rifting directions in the Red Sea and northern Afar. The whole-rock compositions of the erupted lavas and tephra range from trachybasaltic to trachybasaltic andesite, and crystal-hosted melt inclusions contain up to 3,000 ppm of sulphur by weight. The eruption was preceded by significant seismicity, detected by regional networks of sensors and accompanied by sustained tremor. Substantial infrasound was recorded at distances of hundreds to thousands of kilometres from the vent, beginning at the onset of the eruption and continuing for weeks. Analysis of ground deformation suggests the eruption was fed by a shallow, NW-SE-trending dike, which is consistent with field and satellite observations of vent distributions. Despite lack of prior planning and preparedness for volcanic events in the country, rapid coordination of the emergency response mitigated the human costs of the eruption.

Mapping the evolving strain field during continental breakup from crustal anisotropy in the Afar Depression.

Rifting of the continents leading to plate rupture occurs by a combination of mechanical deformation and magma intrusion, yet the spatial and temporal scales over which these alternate mechanisms localize extensional strain remain controversial. Here we quantify anisotropy of the upper crust across the volcanically active Afar Triple Junction using shear-wave splitting from local earthquakes to evaluate the distribution and orientation of strain in a region of continental breakup. The pattern of S-wave splitting in Afar is best explained by anisotropy from deformation-related structures, with the dramatic change in splitting parameters into the rift axis from the increased density of dyke-induced faulting combined with a contribution from oriented melt pockets near volcanic centres. The lack of rift-perpendicular anisotropy in the lithosphere, and corroborating geoscientific evidence of extension dominated by dyking, provide strong evidence that magma intrusion achieves the majority of plate opening in this zone of incipient plate rupture.

Magma-maintained rift segmentation at continental rupture in the 2005 Afar dyking episode.

Seafloor spreading centres show a regular along-axis segmentation thought to be produced by a segmented magma supply in the passively upwelling mantle. On the other hand, continental rifts are segmented by large offset normal faults, and many lack magmatism. It is unclear how, when and where the ubiquitous segmented melt zones are emplaced during the continental rupture process. Between 14 September and 4 October 2005, 163 earthquakes (magnitudes greater than 3.9) and a volcanic eruption occurred within the approximately 60-km-long Dabbahu magmatic segment of the Afar rift, a nascent seafloor spreading centre in stretched continental lithosphere. Here we present a three-dimensional deformation field for the Dabbahu rifting episode derived from satellite radar data, which shows that the entire segment ruptured, making it the largest to have occurred on land in the era of satellite geodesy. Simple elastic modelling shows that the magmatic segment opened by up to 8 m, yet seismic rupture can account for only 8 per cent of the observed deformation. Magma was injected along a dyke between depths of 2 and 9 km, corresponding to a total intrusion volume of approximately 2.5 km3. Much of the magma appears to have originated from shallow chambers beneath Dabbahu and Gabho volcanoes at the northern end of the segment, where an explosive fissural eruption occurred on 26 September 2005. Although comparable in magnitude to the ten year (1975-84) Krafla events in Iceland, seismic data suggest that most of the Dabbahu dyke intrusion occurred in less than a week. Thus, magma intrusion via dyking, rather than segmented normal faulting, maintains and probably initiated the along-axis segmentation along this sector of the Nubia-Arabia plate boundary.