Search

The contribution of mantle upwellings of varying spatial extent to Cenozoic magmatism across Africa is debated because geochemical and seismological tools used to interrogate them are primarily sensitive to either composition or temperature. Thermochemical conditions control the depth at which mantle materials undergo phase changes, which cause seismic discontinuities. Mapping seismic discontinuities across the mantle transition zone (MTZ) and below provides insight into the variable thermochemical nature of upwellings. We present observations of seismic discontinuities beneath Africa obtained from a compilation of P-to-s receiver functions (using Pds, PPds, and PKPds phases), recorded at seismograph networks across Africa between 1990-2019. We exploit a recent high-resolution African continental P-wavespeed model to migrate our receiver functions to depth in a common conversion point stack. Cenozoic magmatism along the East African Rift is largely underlain by a thin MTZ implying a contribution to rift magmatism from sources at or below MTZ depths. The Ethiopian rift is underlain by a depressed d410 and uplifted d660 indicating a moderate positive thermal anomaly at MTZ depths (~100-150K). The southern East African Rift displays a greater d410 depression and a regional d660 depression, suggesting a stronger thermochemical anomaly at MTZ depths. Here, seismic conversions at ~1025km depth are collocated with slow wavespeeds within the African Superplume, corroborating evidence for a compositional anomaly. We suggest that the contribution of a purely thermal plume directly below Ethiopia augments conditions for mantle melting and rifting. Distinct upwellings may also affect the MTZ below Cenozoic magmatism in Cameroon and Madagascar.

©2020. American Geophysical Union. All Rights Reserved. Geochemical observations of ocean island and mid-ocean ridge basalts suggest that abundances of heat-producing elements (HPEs: U, Th, and K) vary within the mantle. Combined with bulk silicate Earth models and constraints on the Earth's heat budget, these observations suggest the presence of a more enriched (potentially deep and undepleted) reservoir in the mantle. Such a reservoir may be related to seismically observed deep mantle structures known as large low shear velocity provinces (LLSVPs). LLSVPs might represent thermochemical piles of an intrinsically denser composition, and many studies have shown such piles to remain stable over hundreds of Myr or longer. However, few studies have examined if thermochemical piles can remain stable if they are enriched in HPEs, a necessary condition for them to constitute an enriched HPE reservoir. We conduct a suite of mantle convection simulations to examine the effect of HPE enrichment up to 25× the ambient mantle on pile stability. Model results are evaluated against present-day pile morphology and tested for resulting seismic signatures using self-consistent potential pile compositions. We find that stable piles can form from an initial basal layer of dense material even if the layer is enriched in HPEs, depending on the density of the layer and degree of HPE enrichment, with denser basal layers requiring increased HPE enrichment to form pile-like morphology instead of a stable layer. Thermochemical piles or LLSVPs may therefore constitute an enriched reservoir in the deep mantle.

Iceland sits astride a mid-ocean ridge underlain by a {mantle} hotspot. The interplay of these two geological processes has the potential to generate a complex and laterally variable crustal structure. The thickness of the Icelandic crust is a long running and controversial debate, with estimates ranging from a "thin'' 20 km crust to a "thick'' 40 km crust. We present new images of the first order seismic discontinuity structure of the Icelandic crust based on a joint inversion of receiver function and ambient noise derived surface wave dispersion data. Inversion results are validated through comparison to receiver functions multi-phase common conversion point stacks across the densely instrumented Northern Volcanic Zone. We find a multi-layered crustal structure consisting of a 6-10 km deep upper crust underlain by either one or two discontinuities. The shallower discontinuity is found at depths of ~20 km throughout Iceland. The deeper discontinuity is only present in some regions, defining the base of a lens-like lower layer with maximum depths of 44 km above the center of the mantle plume. Either of these two discontinuities could be interpreted as the seismic Moho, providing an explanation why previous estimates of crustal thickness have diverged. Such structure may form via underplating of a pre-existing oceanic crust as has been hypothesized in other ocean island plume settings. However we demonstrate with a simple petrological model that variability in seismic discontinuity structure can also be understood as a consequence of compositional variation in melts generated with distance from the plume center.

Constraints on chemical heterogeneities in the upper mantle may be derived from studying the seismically observable impedance contrasts that they produce. Away from subduction zones, several causal mechanisms are possible to explain the intermittently observed X-discontinuity (X) at 230--350\,km depth: the coesite-stishovite phase transition, the enstatite to clinoenstatite phase transition and/or carbonated silicate melting, all requiring a local enrichment of basalt. Africa hosts a broad range of terranes, from Precambrian cores to Cenozoic hotspots with or without lowermost mantle origins. With the absence of subduction below the margins of the African plate for $>$0.5\,Ga, Africa presents an ideal study locale to explore the origins of the X. Traditional receiver function (RF) approaches used to map seismic discontinuities, like common conversion-point stacking, ignore slowness information crucial for discriminating converted upper mantle phases from surface multiples. By manually assessing depth and slowness stacks for 1$^\circ$ radius overlapping bins, normalized vote mapping of RF stacks is used to robustly assess the spatial distribution of converted upper mantle phases. The X is mapped beneath Africa at 233--340\,km depth, revealing patches of heterogeneity proximal to mantle upwellings in Afar, Canaries, Cape Verde, East Africa, Hoggar, and Réunion with further observations beneath Cameroon, Madagascar, and Morocco. There is a lack of an X beneath southern Africa, and strikingly, the magmatic eastern rift branch of the southern East African Rift. With no relationships existing between depth and amplitudes of observed X and estimated mantle temperatures, multiple causal mechanisms are required across a range of continental geodynamic settings.

Africa’s Cenozoic tectonism is often attributed to mantle plumes, particularly below East Africa, but their morphology, number, location, and impact on the African lithosphere are debated. The broad slow wavespeed African Superplume, ubiquitous in large-scale tomographic models, originates below South Africa, reaching the surface somewhere below East Africa. However, whether the diverse East African mantle geochemistry is best reconciled with one heterogeneous upwelling, or current tomographic models lack the resolution to image multiple distinct plumes, remains enigmatic. S-wavespeed tomographic images of Africa are legion, but higher-frequency P-wavespeed whole-mantle models possessing complementary diagnostic capabilities are comparatively lacking. This hinders attempts to disentangle the effects of Cenozoic hotspot tectonism and Pan African (and older) tectonic events on the East African lithosphere. Here we develop a continental-scale P-wave tomographic model capable of resolving structure from upper-to-lower mantle depths using a recently-developed technique to extract absolute arrival-times from noisy, temporary African seismograph deployments. Shallow-mantle wavespeeds are δVP ≈–4% below Ethiopia, but less anomalous (δVP ≥–2%) below other volcanic provinces. The heterogeneous African Superplume reaches the upper mantle below the Kenyan plateau. Below Ethiopia/Afar we image a second sub-vertical slow wavespeed anomaly rooted near the core-mantle boundary outside the African LLVP, meaning multiple disparately sourced whole-mantle plumes may influence East African magmatism. In contrast to other African cratons, wavespeeds below Tanzania are only fast to 90–135km depth. When interpreted alongside Lower Eocene on-craton kimberlites, our results support pervasive metasomatic lithospheric modification caused by subduction during the Neoproterozoic Pan-African orogeny.