Body-wave tomographic imaging of the Turkana Depression: Implications for rift development and plume-lithosphere interactions

The Turkana Depression, a topographically-subdued, broadly-rifted zone between the elevated East African and Ethiopian plateaus, disrupts the N–S, fault-bounded rift basin morphology that characterizes most of the East African Rift. The unusual breadth of the Turkana Depression leaves unanswered questions about the initiation and evolution of rifting between the Main Ethiopian and Eastern rifts. Hypotheses explaining the unusually broad, low-lying area include superposed Mesozoic and Cenozoic rifting and a lack of mantle lithospheric thinning and dynamic support. To address these issues, we have carried out the first body-wave tomographic study of the Depression’s upper mantle. Seismically-derived temperatures at 100 km depth exceed petrological estimates, suggesting the presence of mantle melt, although not as voluminous as the Main Ethiopian Rift, contributes to velocity anomalies. A NW–SE-trending high wavespeed band in southern Ethiopia at urn:x-wiley:15252027:media:ggge22580:ggge22580-math-0001200 km depth is interpreted as refractory Proterozoic lithosphere which has likely influenced the localization of both Mesozoic and Cenozoic rifting. At urn:x-wiley:15252027:media:ggge22580:ggge22580-math-0002100 km depth below the central Depression, a single localized low wavespeed zone is lacking. Only in the northernmost Eastern Rift and southern Lake Turkana is there evidence for focused low wavespeeds resembling the Main Ethiopian Rift, that bifurcate below the Depression and broaden approaching southern Ethiopia further north. These low wavespeeds may be attributed to melt-intruded mantle lithosphere or ponded asthenospheric material below lithospheric thin-spots induced by the region's multiple rifting phases. Low wavespeeds persist to the mantle transition zone suggesting the Depression may not lack mantle dynamic support in comparison to the two plateaus

The development of multiple phases of superposed rifting in the Turkana Depression, East Africa: evidence from receiver functions

The Turkana Depression in Eastern Africa separates the elevated plateaus of East Africa to the south and Ethiopia-Yemen to the north. It remains unclear whether the Depression lacks dynamic mantle support, or if the entire East Africa region is dynamically supported and the Depression compensated isostatically by thinned crust. Also poorly understood is how Miocene-Recent extension has developed across the Depression, connecting spatially separated magmatic rift zones in Ethiopia and Kenya. Receiver function analysis is used to constrain Moho depth and bulk-crustal V P /V S ratio below new seismograph networks in the Depression, and on the northern Tanzania craton. Crustal thickness is ∼40 km below northern Uganda and 30–35 km below southern Ethiopia, but 20–30 km below most of the Depression, where mass-balance calculations reveal low elevations can be explained adequately by crustal thinning alone. Despite the fact that magmatism has occurred for 45 Ma across the Depression, more than 15 Ma before East African Rift (EAR) extension initiated, bulk crustal V P /V S across southern Ethiopia and the Turkana Depression (∼1.74) is similar to that observed in areas unaffected by Cenozoic rifting and magmatism. Evidence for voluminous lower crustal intrusions and/or melt, widespread below the Ethiopian rift and Ethiopian plateau to the north, is therefore lacking. These observations, when reviewed in light of high stretching factors (β ≤ 2.11), suggest Cenozoic extension has been dominated until recently by faulting and plate stretching, rather than magma intrusion, which is likely an incipient process, operating directly below seismically-active Lake Turkana. Early-stage EAR basins to the west of Lake Turkana, with associated stretching factors of β ≈ 2, formed in crust only moderately thinned during earlier rifting episodes. Conversely, ∼23 km-thick crust beneath the Kino Sogo Fault Belt (KSFB) has small offset faults and thin sedimentary strata, suggesting almost all of the observed stretching occurred in Mesozoic times. Despite the KSFB marking the shortest path between focused extensional zones to the north and south, seismicity and GPS data show that modern extension is localized below Lake Turkana to the west. Failed Mesozoic rift zones, now characterized by thinned crust and relatively refractory mantle lithosphere, are being circumnavigated, not exploited by EAR rifting

Subsurface images of the Eastern Rift, Africa, from the joint inversion of body waves, surface waves and gravity: investigating the role of fluids in early-stage continental rifting

International audienceThe Eastern Rift System (ERS) of northern Tanzania and southern Kenya, where a cratonic lithosphere is in the early stages of rifting, offers an ideal venue for investigating the roles of magma and other fluids in such an environment. To illuminate these roles, we jointly invert arrival times of locally recorded P and S body waves, phase delays of ambient noise generated Rayleigh waves and Bouguer anomalies from gravity observations to generate a 3-D image of P and S wave speeds in the upper 25 km of the crust. While joint inversion of gravity and arrival times requires a relationship between density and wave speeds, the improvement in resolution obtained by the combination of these disparate data sets serves to further constrain models, and reduce uncertainties. The most significant features in the 3-D model are (1) P and S wave speeds that are 10–15 per cent lower beneath the rift zone than in the surrounding regions, (2) a relatively high wave speed tabular feature located along the western edge of the Natron and Manyara rifts, and (3) low (∼1.71) values of Vp/Vs throughout the upper crust, with the lowest ratios along the boundaries of the rift zones. The low P and S wave speeds at mid-crustal levels beneath the rift valley are an expected consequence of active volcanism, and the tabular, high-wave speed feature is interpreted to be an uplifted footwall at the western edge of the rift. Given the high levels of CO2 outgassing observed at the surface along border fault zones, and the sensitivity of Vp/Vs to pore-fluid compressibility, we infer that the low Vp/Vs values in and around the rift zone are caused by the volcanic plumbing in the upper crust being suffused by a gaseous CO2 froth on top of a deeper, crystalline mush. The repository for molten rock is likely located in the lower crust and upper mantle, where the Vp/Vs ratios are significantly higher

Lithospheric modification by extension and magmatism at the craton-orogenic boundary: North Tanzania Divergence, East Africa

We present a joint analysis of newly acquired gravity and teleseismic data in the North Tanzanian Divergence, where the lithospheric break-up is at its earliest stage. The impact of a mantle upwelling in more mature branches of the East African Rift has been extensively studied at a lithospheric scale. However, few studies have been completed that relate the deep-seated mantle anomaly detected in broad regional seismic tomography with the surface deformation observed in the thick Archaean Pan-African suture zone located in North Tanzania. Our joint inversion closes the gap between local and regional geophysical studies, providing velocity and density structures from the surface down to ca. 250 km depth with new details. Our results support the idea of a broad mantle upwelling rising up to the lithosphere and creating a thermal modification along its path. However, our study clearly presents an increasing amplitude of the associated anomaly both in velocity and density above 200 km depth, which cannot be solely explained by a temperature rise. We infer from our images the combined impact of melt (2-3 per cent), composition and hydration that accompany the modification of a thick heterogenous cratonic lithosphere are a response to the hot mantle rising. The detailed images we obtained in density and velocity assert that Archaean and Proterozoic units interact with the mantle upwelling to restrict the lithosphere modifications within the Magadi-Natron-Manyara rift arm. The composition and hydration variations associated with those units equilibrate the thermal erosion of the craton root and allow for its stability between 100 and 200 km depth. Above 80 km depth, the crustal part is strongly affected by intruding bodies (melt and gas) which produces large negative anomalies in both velocity and density beneath the main magmatic centres. In addition to the global impact of a superplume, the velocity and density anomaly pattern suggests a 3-D distribution of the crust and mantle lithospheric stretching, which is likely to be controlled by inherited fabrics and enhanced by lateral compositional and hydration variations at the Tanzanian craton-orogenic belt boundary

Active deformation constraints on the Nubia-Somalia plate boundary through heterogenous lithosphere of the Turkana depression

Abstract The role of lithospheric heterogeneities, presence or absence of melt, local and regional stresses, and gravitational potential energy in strain localization in continental rifts remains debated. We use new seismic and geodetic data to identify the location and orientation of the modern Nubia‐Somalia plate boundary in the 300‐km‐wide zone between the southern Main Ethiopian Rift (MER) and Eastern Rift (ER) across the Mesozoic Anza rift in the Turkana Depression. This region exhibits lithospheric heterogeneity, 45 Ma‐Recent magmatism, and more than 1,500 m of base‐level elevation change, enabling the assessment of strain localization mechanisms. We relocate 1716 earthquakes using a new 1‐D velocity model. Using a new local magnitude scaling with station corrections, we find 1 ≤ ML ≤ 4.5, and a b‐value of 1.22 ± 0.06. We present 59 first motion and 3 full moment tensor inversions, and invert for opening directions. We use complementary geodetic displacement vectors and strain rates to describe the geodetic strain field. Our seismic and geodetic strain zones demonstrate that only a small part of the 300 km‐wide region is currently active; low elevation and high‐elevation regions are active, as are areas with and without Holocene magmatism. Variations in the active plate boundary's location, orientation and strain rate appear to correspond to lithospheric heterogeneities. In the MER‐ER linkage zone, a belt of seismically fast mantle lithosphere generally lacking Recent magmatism is coincident with diffuse crustal deformation, whereas seismically slow mantle lithosphere and Recent magmatism are characterized by localized crustal strain; lithospheric heterogeneity drives strain localization

Lithospheric low-velocity zones associated with a magmatic segment of the Tanzanian Rift, East Africa

International audienceRifting in a cratonic lithosphere is strongly controlled by several interacting processes including crust/mantle rheology, magmatism, inherited structure and stress regime. In order to better understand how these physical parameters interact, a 2 yr long seismological experiment has been carried out in the North Tanzanian Divergence (NTD), at the southern tip of the eastern magmatic branch of the East African rift, where the southward-propagating continental rift is at its earliest stage. We analyse teleseismic data from 38 broad-band stations ca. 25 km spaced and present here results from their receiver function (RF) analysis. The crustal thickness and Vp/Vs ratio are retrieved over a ca. 200 × 200 km2 area encompassing the South Kenya magmatic rift, the NTD and the Ngorongoro-Kilimanjaro transverse volcanic chain. Cratonic nature of the lithosphere is clearly evinced through thick (up to ca. 40 km) homogeneous crust beneath the rift shoulders. Where rifting is present, Moho rises up to 27 km depth and the crust is strongly layered with clear velocity contrasts in the RF signal. The Vp/Vs ratio reaches its highest values (ca. 1.9) beneath volcanic edifices location and thinner crust, advocating for melting within the crust. We also clearly identify two major low-velocity zones (LVZs) within the NTD, one in the lower crust and the second in the upper part of the mantle. The first one starts at 15–18 km depth and correlates well with recent tomographic models. This LVZ does not always coexist with high Vp/Vs ratio, pleading for a supplementary source of velocity decrease, such as temperature or composition. At a greater depth of ca. 60 km, a mid-lithospheric discontinuity roughly mimics the step-like and symmetrically outward-dipping geometry of the Moho but with a more slanting direction (NE–SW) compared to the NS rift. By comparison with synthetic RF, we estimate the associated velocity reduction to be 8–9 per cent. We relate this interface to melt ponding, possibly favouring here deformation process such as grain-boundary sliding (EAGBS) due to lithospheric strain. Its geometry might have been controlled by inherited lithospheric fabrics and heterogeneous upper mantle structure. We evidence that crustal and mantle magmatic processes represent first order mechanisms to ease and locate the deformation during the first stage of a cratonic lithospheric breakup

Fault-magma interactions during early continental rifting: Seismicity of the Magadi-Natron-Manyara basins, Africa

International audienceAlthough magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2-D studies. We analyze seismicity data from a 13 month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3-D velocity model reveal lower crustal earthquakes beneath the central basins and along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87±0.03. Focal mechanisms for 65 earthquakes, and 13 from a catalogue prior to our array reveal an along-axis stress rotation of ∼60° in the magmatically active zone. New and prior mechanisms show predominantly normal slip along steep nodal planes, with extension directions ∼N90°E north and south of an active volcanic chain consistent with geodetic data, and ∼N150°E in the volcanic chain. The stress rotation facilitates strain transfer from border fault systems, the locus of early-stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Results indicate that earthquakes are largely driven by stress state around inflating magma bodies