April 2017 LIP of the Month

The Eastern rift branch of East Africa: an LIP-style continental rift exhibiting massive magma and volatile production

James D. Muirhead1, Cynthia J. Ebinger2, Simon A. Kattenhorn3, Hyunwoo Lee4, Tobias P. Fischer4, Steven Roecker5, Sara Jaye Oliva2

1Syracuse University, New York, USA

(James Muirhead jmuirhead@ymail.com; https://jamesdmuirhead.wordpress.com/)

2Tulane University, Louisiana, USA

3University of Alaska Anchorage, Alaska, USA

4University of New Mexico, New Mexico, USA

5Rensselaer Polytechnic Institute, New York, USA

The following is based on a series of recent publications that include, but are not limited to, Lee et al. (2016), Muirhead et al. (2016), Lee et al. (2017), Roecker et al. (2017), and Weinstein et al. (in review). For full details, please refer to these papers.

  1. Introduction

Continental rift zones are sites of lithospheric thinning and heating, commonly accompanied by magmatism and volatile transfer to the lithosphere and atmosphere. Some rifts evolve to seafloor spreading, with the heavily intruded and stretched passive margin increasing the surface area and volume of continental plates. Quantifying the flux of magma and volatiles during rifting is thus fundamental to our understanding of the growth of continental lithosphere and the deep carbon cycle.

Although some continental rifts lack the magma production rates of large igneous provinces (LIPs), the Eastern rift branch of the East African Rift System (EARS) fits under the category of an LIP (Ernst, 2014), and hence is consider here to be an “LIP-style continental rift”. The Kenya Rift sector alone has eruptive volumes in excess of 100,000 km3 (Guth, 2015) and has developed since ~25 Ma in an intra-plate continental setting (Ebinger et al., 2000), albeit within a nascent divergent margin. Although intruded magma volumes for the Kenya Rift are poorly constrained, they are likely significant even during earliest rift stages (first few million years), as illustrated by the multi-disciplinary CRAFTI (continental rifting in Africa: fluid-tectonic involvement) project described below.

Quantifying the fluxes of magma and volatiles during the earliest stages of rifting can help to address outstanding fundamental questions regarding how continental rifts initiate and evolve to seafloor spreading. Force balance analyses indicate that the break up of thick, strong lithosphere in East Africa requires more than tectonic forces (Bott et al., 1991), and thus continental rifting is inferred to be assisted by buoyancy forces and heat transfer from ascending magmatic intrusions (Buck, 2004; Bialas et al., 2010), and possibly from fluid-driven lithospheric weakening via metasomatism (e.g., Ebinger et al., 1997; Vauchez et al., 2005; Lindenfeld et al., 2012). However, to date no study has quantified the flux of magma and related volatiles into early-stage rifts in the Eastern rift branch, or demonstrated the importance of these rift components for strain localization processes during incipient rifting.

The CRAFTI project in the Natron-Magadi basin (Kenya-Tanzania border) has produced a complementary, multi-disciplinary dataset that allows us to probe the inner workings of continental rifts in the Eastern rift branch. Through this dataset we have constrained, for the first time, a model for the involvement magma and magmatic volatiles across a ~40 km-deep crustal section of incipient continental rifting. The results of geophysical studies constrain crustal structure and strain patterns, and mantle studies are underway (Plasman et al., 2017; Roecker et al., 2017). Below we summarize major findings from recently published geological, geophysical, and geochemical studies in the region. From these studies, we can make the following broad conclusions regarding early-stage rifting in the Eastern rift branch:

  • Early-stage rifting (<7 Ma) in the Eastern rift branch is accompanied by magma accumulation at the crust mantle boundary (Plasman et al., 2017). These magmas rise to ca. 12 km depth at the Natron and Magadi basins, with significant volumes emplaced in the lower crust beneath rifted basins, as well as below uplifted horsts between basins (Roecker et al., 2017; Weinstein et al., in review).
  • These magmas are generated in the sub-continental lithosphere, and their emplacement at crustal levels contributes to crustal weakening through heating and hydraulic fracturing of surrounding rocks (Lee et al., 2017). The magmatic intrusions also release massive volumes of magmatic CO2 along border and intra-rift faults at large distances from volcanoes (Lee et al., 2016).
  • The intrusion of massive volumes of magma and subsequent release of magmatic volatiles plays a critical role in strain localization and rift evolution processes in the Eastern rift branch of East Africa (Muirhead et al., 2016).
  1. Quantifying the flux of magmatic CO2 across rift-wide fault systems

To test if rift-wide faults transport magmatic volatiles to the surface, diffuse CO2 flux was measured in the Natron and Magadi basins with an accumulation chamber in fault zones (Lee et al., 2016) (Figs. 1 and 2). Gas samples were collected in pre-evacuated glass vials and analyzed for δ13C-CO2 (‰) and CO2 concentrations to constrain the source of the observed CO2 (i.e., biogenic or magmatic).

Carbon isotope compositions (δ13C-CO2 values range from -3.8 to -11.7‰) and elevated CO2 concentrations (mean = 929 ± 386 ppm) indicate a strong magmatic contribution to the observed CO2 (Fig. 2). High CO2 flux data in fault zones (mean values in Magadi of 36.6±11.0 g m-2 d-1) suggests that faults act as permeable pathways that facilitate the ascent of this CO2. CO2 flux data (Fig. 2) support annual outputs of 4.05±1.9 Mt yr-1 in the ~180 km-long Natron-Magadi basin (Lee et al., 2016). If this flux were applied to the entire Eastern rift branch, it would equate to an annual flux of 71 ± 33 Mt yr-1 of CO2, which could potentially increase the global flux from natural systems by 11% to 708 Mt yr-1 (estimates from Burton et al. (2013)). These results suggest the EARS is an important contributor to the deep carbon cycle, and reveal that significant volumes of magmatic volatiles are released through fault systems during early-stage rifting.

Figure 1: Study site of the CRAFT experiment in the Natron and Magadi basins near the Kenya-Tanzania border. (a) SRTM map of the Magadi–Natron basin showing active volcanic centers (Suswa, Gelai, and Oldoinyo Lengai). Also shown are diffuse CO2 measurement locations (green stars) and earthquake epicenters (red circles) from the CRAFTI seismic network (purple triangles; not all stations shown). (b) Location of the study area in the EARS. (c) Cross-sections in the center of the Magadi and Natron basins from A–A’ and B–B’ annotated in (a). Purple circles represent earthquake foci that are located 2 km from each section and yellow fill represents basin sediments. Figure from Lee et al. (2016).

Figure 2: CO2 flux and geochemistry from diffuse soil CO2 analyses in the Natron-Magadi basins from Lee et al. (2016). A: Contoured CO2 flux map (25 m grids) in a faulted graben in the Magadi study area. The map was generated using the geo-statistical program Surfer. Black lines represent faults. B: Carbon isotope (δ13C-CO2) vs the reciprocal of CO2 concentration from Lee et al. (2016). Mantle C (-6.5±2.5‰) and biogenic CO2 values (-20 to -25‰) values are from Sano and Marty (1995) and Lewicki and Brantley (2000), respectively. Arrows show mixing lines between the end-member values. Natron and Magadi basin CO2 data fit on a mixing line between air and mantle C, suggesting that the soil CO2 is sourced from crustal magmas.

  1. Spring systems

Thermal spring samples (T = 36.8–83.5 °C) in the Natron-Magadi basin were analyzed to investigate the sources of magmas and fluids in the rift, as well as their potential interactions with the surrounding crust (Lee et al., 2017) (Fig. 3). Helium isotope data support both mantle (10-66%) and crustal contributions to the spring fluids. 3He/4He ratios (R/Ra, Ra = 1.4 × 10-6), reaching as high as 4.0 Ra, are close to those observed for subcontinental lithospheric mantle-derived volatiles (6.1 ± 0.9 Ra), instead of MORB-type (8.0 Ra), which is consistent with the model suggested by Mana et al. (2015).

Production rates of 4He in the Natron and Magadi basins range from 3.74 × 1011 to 3.34 × 1014 atoms m2 s1, and are up to four orders of magnitude greater than mean production rates (4.18 × 1010 atoms m2 s1) from continental settings on Earth. These high 4He fluxes therefore cannot be attributed solely to alpha decay of U and Th. It is instead likely that the elevated production of crust-derived 4He is driven by heating and fracturing of country rock in response to magma intrusion at depth (Lee et al., 2017).

Figure 3: (a) Annotated digital elevation map of the Natron and Magadi basins, showing the distribution of springs (yellow stars) from Muirhead et al. (2016). Samples collected for geochemical analyses are indicated by green stars. Black lines are border faults. (b) Plot of R/Ra versus 4He/20Ne ratios from Lee et al. (2017) for samples from the Magadi and Natron basins as well as select locations along the Kenya rift. Mixing lines represent binary mixing relationships between air-saturated water and end member values, which have different fractions of crustal and mantle He (represented as % fmantle). Sub-continental lithospheric mantle is selected as the mantle source.

  1. Crustal structure of an early-stage rift basin

The 39-station CRAFTI seismic array located 3427 events in the study region with magnitudes ranging from 0.4 to 4.4 (Weinstein et al., in review) (Fig. 1). Moho depths, determined from receiver functions, transition from 41 km beneath the rift flank to 27 km below the rift valley (Plasman et al., 2017). Joint inversion of ambient noise, arrival time, and gravity data was performed by Roecker et al. (2017) to image Vp, Vs, and Vp/Vs variations within the crust, with resolution to about 25 km. Receiver function models indicate high Vp/Vs at the sometimes highly reflective crust-mantle boundary, suggesting sill intrusions, or magmatic underplate. Beneath the central Natron and Manyara basins and the Crater Highlands to the west, low S-wave speeds are observed at mid to lower crustal depths, and attributed to the presence of partial melt (Roecker et al., 2017). Seismicity occurs in these areas of magmatic underplate (20-30 km), suggesting that the intrusion process is ongoing (Weinstein et al., in review). Seismicity also occurs in pancake-shaped swarms in the mid-crust and along the boundaries of the large volume low-velocity zone beneath Oldoinyo Lengai and Gelai volcanoes (Weinstein et al., in review). Unusually low Vp/Vs ratios characterize the footwall of the Natron basin border fault, which is actively degassing large volumes of mantle-derived CO2 (Lee et al., 2016). The low Vp/Vs could be caused by CO2-filled rocks, and may be a diagnostic feature of active degassing (Roecker et al., 2017).

  1. Magmatic fluid release and strain localization

Fault systems in the Eastern rift branch are classically sub-divided into large (throws >1 km) border faults, which bound the half-graben basins, and smaller intra-rift faults, which occur along the basin floor in the rift center. As rift basins evolve to sea-floor spreading ridges, the locus of strain accommodation gradually migrates, or jumps, from border faults into intra-rift fault systems, which are underlain by zones of magmatism where dike intrusion accommodates a large proportion of the plate-boundary strain (Ebinger and Casey, 2001; Keir et al., 2006; Corti et al., 2009; Muirhead et al., 2016).

We investigated the processes controlling this transition from border-fault controlled rifting (asymmetric rift) to strain localization to the rift center (intra-rift faulting with accompanying dike intrusion at depth) through an analysis of the Natron and Magadi basins (Muirhead et al., 2016) (Fig. 4). Remote-sensing analyses (0.5 m aerial photography, and 30 m SRTM data) of intra-rift fault systems were employed to quantify the distribution of fault strain. 40Ar-39Ar dating of faulted lavas constrained the ages of the faults, and were then used to estimate time-averaged extension rates in the fault populations (Muirhead et al., 2016). Results revealed that intra-rift faulting in the Magadi basin accommodate 1.53-1.73 mm yr-1 of extension, whereas faults in the Natron basin accommodate 0.33-0.61 mm yr-1 of extension (Fig. 4). Accounting for extensional strains related to basin-scale flexure of border fault hanging walls (described in detail in Muirhead et al. (2016)), it is estimated that intra-rift faults in Magadi accommodate 1.34-1.60 mm yr-1 of the total 2 mmyr-1 of regional extension across the basin (estimates from Saria et al., 2014) (Fig. 5). In contrast, intra-rift faults in the Natron basin account for only 0.02-0.56 mm yr-1 of the 1.8 mmyr-1 of regional extension (Fig. 5). These strain analyses, however, cannot constrain what proportion of this extension is accommodated via dike intrusion at depth. For example, the 2007 dike intrusion below Gelai volcano at the southern end of the Natron basin is shown to have accommodated the majority of extensional strain during this discrete rifting event (Calais et al., 2008), with normal faulting accommodating extension in only the upper few kilometers of the crust. As static stress changes above the dike tips are shown to initiate faulting and graben formation in rifting environments (Rubin and Pollard, 1988), extensional strains measured at the surface in the Natron and Magadi basins are expected to mirror that accommodated at depth regardless of what proportion is accommodated through either dike intrusion or normal faulting (Muirhead et al., 2016).

In all, our fault system analyses suggest that strain in the 7 Ma Magadi basin has focused into the intra-rift fault system, which is underlain by a magma intrusion zone above the mid-crustal magma bodies imaged in seismic studies (Roecker et al., 2017; Weinstein et al., in review). In the 3 Ma Natron basin, time-averaged strain is primarily concentrated along the border fault (1.24-1.78 mm yr-1 extension rate) (Fig. 5). Also observed in the region is a clear link between areas of high strain and hydrothermal fluid release (Figs. 3 and 4). In the Natron basin, hot springs are observed along the active rift borders, whereas the focusing of strain into a system of intra-rift faults in Magadi has occurred with an accompanying transition to fluid release in the rift center. Geochemical analyses of spring fluids and diffusely degassing CO2 along these fault systems support significant magmatic contributions to these fluids (Figs. 2 and 3). Similar patterns are seen in active deformation patterns from discrete magmatic rifting episodes (e.g., Calais et al., 2008), and the inter-magmatic period (Weinstein et al., in review). Therefore, we concluded that the flow of magmatic fluids plays an important role in focusing upper crustal fault strain to the rift center in early-stage continental rift basins, thereby assisting in the strain localization processes that drive rift basin evolution and along-axis segmentation (i.e., Ebinger, 2005).


Figure 4: STRM image of the Natron and Magadi basins. Bold black line represents the border fault for each basin. Rift-normal distribution of intra-rift fault-related strain across two 60 km-long transects in the Magadi and Natron basins. (a) Annotated SRTM image showing the location of each transect and the polygons in which strain, and corresponding extension rate, was calculated. (b) Cumulative extension rate in the intra-rift fault population of the Magadi (X-X’). (c) Cumulative extension rate in the intra-rift fault population of the Natron (Y-Y’) basin. Extension rates in each polygon in (a) are summed from west to east, resulting in total time-averaged extension rates of 1.63±0.10 mm yr-1 (b) and 0.5±0.03 mm yr-1 (c). Figure from Muirhead et al. (2016).

Figure 5: Summary of strain data modified from Muirhead et al. (2016). (A) Distribution of fault strain in the 7 Ma Magadi basin. Black lines and gray lines represent border faults and intra-rift faults, respectively. Based on estimated extension rates, the greatest proportion of regional extension (i.e., modeled total strain) is accommodated in the intra-rift fault system (67%–80% of total strain). Analyses from Muirhead et al. (2016) give a mean border fault extension rate ranging from 0.40 to 0.66 mm yr–1. (B) Distribution of fault strain in the 3 Ma Natron basin. Based on estimated extension rates, the greatest proportion of regional extension (i.e., modeled total strain) is accommodated along the border fault (69%–99% of total strain). Analyses from Muirhead et al. (2016) give a mean border fault extension rate ranging from 1.24 to 1.78 mm yr–1.

  1. LIP-style continental rifting: continental breakup, crustal accretion and volatile exchange

Results from the CRAFTI project reveal the importance of magmatism and magmatic volatile release during continental rifting (Fig. 6). Indeed, LIP-style continental rifting, and its associated magmatism, is critical for the growth of continental lithosphere, volatile exchange between the lithosphere and atmosphere, and initiating continental breakup.

Figure 6: Conceptual summary figure illustrating rift processes at crustal depths in the early-stage (<7 Ma) Natron and Magadi basins in the Eastern rift branch of East Africa. Low S-wave speeds are observed at mid to lower crustal depths and attributed to the presence of partial melt (Roecker et al., 2017). Seismicity occurs in these areas of magmatic underplate (20-30 km), suggesting that the intrusion process is ongoing (Weinstein et al., in review). Discrete magma chambers developing in the lower to middle crust feed smaller dike-sill complexes at upper crustal levels (<12 km depth; Weinstein et al., in review). Large volumes of mantle-derived CO2 are emitted from zones exhibiting high time-averaged strain and unusually low Vp/Vs ratios (Lee et al., 2016; Muirhead et al., 2016; Roecker et al., 2017), which collectively suggest a feedback between fluid flow and strain localization in the upper crust of early-stage continental rift basins (Muirhead et al., 2016). In all, these data suggest that early-stage rifting in the Eastern rift branch is accompanied by volumetrically significant magmatism and magmatic volatiles.

The earliest stages of rifting in the Eastern rift branch are associated with volumetrically significant magmatic underplating beneath the faulted rift valleys and their uplifted flanks. Early in the basin’s formation, discrete magma chambers develop in the lower to middle crust, and they feed smaller dike-sill complex formation at upper crustal levels (<12 km depth; Calais et al., 2008; Muirhead et al., 2015; Weinstein et al., in review) (Fig. 1). In the Natron-Magadi basin, high-resolution tomographic models indicate that a single ~30 km-wide magma chamber feeds the active carbonatitic volcano (Oldoinyo Lengai) as well as stacked sills beneath a monogenetic cone complex (Roecker et al., 2017). The volume of the melt is largely unconstrained, as is the flux rate into the lower crust to sustain this persistent partial melt zone. Extruded volumes for the Kenya Rift are roughly estimated at 310,000 km3 (Guth, 2015). Given that ratios between the volumes of intruded and extruded products are shown to be greater than 3:1 for LIPs and volcanic rifted margins (Coffin and Eldholm, 1994), intruded volumes in the ~600 km-long Kenya Rift could easily exceed 300,000 km3. Assuming an average age of 15 Ma for basins along this part of the rift, these crude volume estimates suggest that magma flux rates could be on the order of 10s of km3/km/Myr, which is comparable to some volcanic arc systems (23-200 km3/km/Myr estimated for intraoceanic arcs; Jicha and Jagoutz, 2015).

Volatile production rates in the East African Rift are also high. Active volcanoes in EARS represents some Earth’s largest gas emitters, with Oldoinyo Lengai and Nyiragongo emitting 2.4 and 3.4 Mtyr-1 of CO2, respectively (Brantley and Koepenick, 1995; Sawyer et al., 2008). These persistently active volcanoes represent a commonly investigated source of volatile transfer to the atmosphere. A more hidden mechanism of volatile exchange occurs via CO2 ascent along deeply penetrating extensional fault systems, which facilitate degassing of volumetrically significant partial melt detected at middle to lower crustal depths (Lee et al., 2016; Roecker et al., 2017) (Fig. 1). This tectonic degassing represents a poorly constrained source of deeply-derived CO2 release in extensional rift systems globally. However, in the EARS it is clearly responsible for CO2 release on the order of 10s of Mtyr-1. These observations suggest that magma-rich continental rifts represent an important source of naturally occurring and potentially climate modulating CO2. Neglecting the short and rapid environmental changes that accompany plume-related magmatism during rift initiation (Burgess and Bowring, 2015), the massive and prolonged CO2 emissions associated with LIP-style continental rifting could contribute to prolonged greenhouse conditions like those of the Cretaceous (Kidder and Worsley, 2010; Ernst, and Youbi, 2017).

The forces driving rift initiation (e.g., gravitational potential energy, basal shear from mantle convection) are theoretically too small to rupture the thick, slow-moving continental lithosphere beneath East Africa (e.g., Bott et al.,1991; Buck, 2004; Stamps et al., 2014). Our work suggests that additional processes facilitate the earliest stages of continental rifting in East Africa: (1) the widespread percolation of mantle-derived fluids that metasomatize and weaken the mantle lithosphere (Lee et al., 2016; Muirhead et al., 2016; Weinstein et al., in review); (2) heat transfer from progressively shallow magma intrusions (Lee et al., 2017; Roecker et al., 2017); and (3) the buoyancy forces of an upwelling low-density mantle plume (e.g., Ebinger and Sleep, 1998; Stamps et al., 2014). The CRAFTI project also represents the first cross-disciplinary study constraining the importance of the above processes for localizing strain and driving rift basin evolution. Specifically, the high Vp/Vs and seismicity in the lower crust indicate that intrusion is ongoing, consistent with the ongoing and volumetrically large CO2 degassing (Lee et al., 2016). The Natron basin appears to be at a transition from border-fault controlled rifting to intrusion and shallow faulting above a zone of magma intrusion (Roecker et al., 2017; Weinstein et al., in press). Lastly, areas of long-term extensional fault strain coincide with zones of magmatic volatile release at the surface, suggesting a feedback between fluid flow and strain localization in the upper crust of early-stage continental rift basins (Muirhead et al., 2016).

Click to open/close ReferencesReferences

Bialas, R.W., Buck, W.R., and Qin, R., 2010, How much magma is required to rift a continent?: Earth and Planetary Science Letters, v. 292, p. 68–78, doi: 10.1016/j.epsl.2010.01.02.

Bott, M. H. P., Stein, S., Wortel, R., Kusznir, N., Fleitout, L., Zoback, M. L., and Cartwright, J., 1991, Sublithospheric Loading and Plate-Boundary Forces [and Discussion]: Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, v. 337, p. 83-93.

Brantley, S. L., and Koepenick, K. W., 1995, Measured carbon-dioxide emissions from Oldoinyo Lengai and the skewed distribution of passive volcanic fluxes: Geology, v. 23, p. 933-936.

Buck, W.R., Einarsson, P., and Brandsdottir, B., 2006, Tectonic stress and magma chamber size as controls on dike propagation: Constraints from the 1975–1984 Krafla rifting episode: Journal of Geophysical Research, v. 111, B12404, doi: 10.1029/2005jb003879.

Burgess, S.D. and Bowring, S.A., 2015, High-precision geochronology confirms voluminous magmatism before, during, and after Earth’s most severe extinction: Science Advances, v. 1, doi: 10.1126/sciadv.1500470.

Burton, M. R., Sawyer, G. M., and Granieri, D., 2013, Deep carbon emissions from volcanoes: Rev. Mineral. Geochem, v. 75, p. 323-354.

Calais, E., et al., 2008, Strain accommodation by slow slip and dyking in a youthful continental rift, East Africa: Nature, v. 456, p. 783–787, doi: 10.1038/nature07478.

Coffin, M.F., and Eldholm, O., 1994, Large igneous provinces: crustal structure, dimensions, and external consequences: Reviews of Geophysics, v. 32, p. 1-36.

Corti, G., 2009, Continental rift evolution: From rift initiation to incipient break-up in the Main Ethiopian Rift, East Africa: Earth-Science Reviews, v. 96, p. 1–53.

Ebinger, C.J., 2005, Continental break-up: The East African perspective: Astronomy & Geophysics, v. 46, 2390, doi: 10.1111/j.1468-4004.2005.46216.x.

Ebinger, C. J., and Sleep, N. H., 1998, Cenozoic magmatism throughout east Africa resulting from impact of a single plume: Nature, v. 395, p. 788-791, doi: 10.1038/27417.

Ebinger, C.J., and Casey, M., 2001, Continental breakup in magmatic provinces: An Ethiopian example: Geology, v. 29, p. 527–530.

Ebinger, C.J., Djomani, Y.P., Mbede, E., Foster, A., and Dawson, J.B., 1997, Rifting Archaean lithosphere: The Eyasi-Manyara-Natron rifts, East Africa: Journal of the Geological Society [London], v. 154, p. 947–960, doi: 10.1144/gsjgs.154.6.0947.

Ebinger, C.J., Yemane, T., Harding, D.J., Tesfaye, S., Kelley, S., and Rex, D.C., 2000, Rift deflection, migration, and propagation: Linkage of the Ethiopian and Eastern rifts, Africa: Geological Society of America Bulletin, v. 112, p. 163–176.

Ernst, R. E., 2014, Large Igneous Provinces: Cambridge University Press, 666 p.

Ernst, R.E. and Youbi, N., 2017, How Large Igneous Provinces affect global climate, sometimes cause mass extinctions, and represent natural markers in the geological record. Palaeogeography, Palaeoclimatology, Palaeoecology, (in press), doi: 10.1016/j.palaeo.2017.03.014

Guth, A.L., 2015, Volcanic volumes associated with the Kenya Rift: Recognition and correlation of preservation biases, in Wright, T.J., et al., eds., Magmatic Rifting and Active Volcanism: Geological Society of London Special Publication 420, p. 43–67, doi: 10.1144/sp420.3.

Jicha, B.R. and Jagoutz, O., 2015, Magma production rates for intraoceanic arcs: Elements, v. 11, p. 105-111.

Keir, D., Ebinger, C.J., Stuart, G.W., Daly, E., and Ayele, A., 2006, Strain accommodation by magmatism and faulting as rifting proceeds to breakup: Seismicity of the northern Ethiopian rift: Journal of Geophysical Research, v. 111, B05314, doi: 10.1029/2005jb003748.

Kidder, D.L. and Worsley, T.R., 2010, Phanerozoic large igneous provinces (LIPs), HEATT (haline euxinic acidic thermal transgression) episodes, and mass extinctions: Palaeogeography, Palaeoclimatology, Palaeoecology, v. 295, p. 162-191.

Lee, H., Muirhead, J.D., Fischer, T.P., Ebinger, C.J., Kattenhorn, S.A., Sharp, Z.D., and Kianji, G., 2016, Massive and prolonged deep carbon emissions associated with continental rifting: Nature Geoscience, v. 9, p. 145–149, doi: 10 .1038 /ngeo2622

Lee, H., Fischer, T.P., Muirhead, J.D., Ebinger, C.J., Kattenhorn, S.A., Sharp, Z.D., Kianji, G., Takahata, N. and Sano, Y., 2017, Incipient rifting accompanied by the release of subcontinental lithospheric mantle volatiles in the Magadi and Natron basin, East Africa: Journal of Volcanology and Geothermal Research, doi: 10.1016/j.jvolgeores.2017.03.017.

Lindenfeld, M., Rumpker, G., Link, K., Koehn, D., and Batte, A., 2012, Fluid-triggered earthquake swarms in the Rwenzori region, East African Rift—Evidence for rift initiation: Tectonophysics, v. 566, p. 95–104, doi: 10.1016/j.tecto .2012.07.010.

Mana, S., Furman, T., Turrin, B.D., Feigenson, M.D., and Swisher, C.C., 2015, Magmatic activity across the East African North Tanzanian Divergence Zone: Journal of the Geological Society [London], v. 172, p. 368–389, doi: 10.1144/jgs2014-072.

Muirhead, J.D., Kattenhorn, S.A., and Le Corvec, N., 2015, Varying styles of magmatic strain accommodation in the East African Rift: Geochemistry, Geophysics, Geosystems, v. 16, p. 2775–2795, doi: 10.1002/2015GC005918.

Muirhead, J.D., Kattenhorn, S.A., Lee, H., Mana, S., Turrin, B.D., Fischer, T.P., Kianji, G., Dindi, E., and Stamps, D.S., 2016, Evolution of upper crustal faulting assisted by magmatic volatile release during early-stage continental rift development in the East African Rift: Geosphere, v. 12, no. 6, p. 1–31, doi:10.1130

Plasman, M., Tiberi, C., Ebinger, C., Roecker, S., Gautier, S., Albaric, J., Peyrat, S., Perrot, S., Déverchère, J., LeGall., B., Tarits, P., Hautot, S., Ferdinand-Wambura, R., Muzuka, A., Mulibo, G., Mtelela, K., Msabi, M., Kianji, G., Rasendra, N., and Gama, R., 2017, Lithospheric structure of the North Tanzanian Divergence, East African rift, estimated from receiver functions: Geophysical Journal International, in press.

Roecker, S., Ebinger, C., Tiberi, C., Mulibo, G., Ferdinand-Wambura, R., Mtelela, K., Muzuka, A., Kianji, G., Gautier, S., Albaric, J., and Peyrat, S., 2017, 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: Geophysical Journal International, in press.

Rubin, A.M., and Pollard, D.D., 1988, Dike-induced faulting in rift zones of Iceland and Afar: Geology, v. 16, p. 413–417.

Saria, E., Calais, E., Stamps, D.S., Delvaux, D., and Hartnady, C.J.H., 2014, Present-day kinematics of the East African Rift: Journal of Geophysical Research, v. 119, p. 3584–3600, doi: 10.1002/2013JB010901.

Sawyer, G. M., Carn, S. A., Tsanev, V. I., Oppenheimer, C., and Burton, M., 2008, Investigation into magma degassing at Nyiragongo volcano, Democratic Republic of the Congo: Geochemistry, Geophysics, Geosystems, v. 9, p. Q02017, doi:10.1029/2007gc001829.

Stamps, D.S., Flesch, L.M., Calais, E., and Ghosh, A., 2014, Current kinematics and dynamics of Africa and the East African Rift System: Journal of Geophysical Research, v. 119, p. 5161–5186, doi: 10.1002/2013JB010717.

Vauchez, A., Dineur, F., and Rudnick, R., 2005, Microstructure, texture and seismic anisotropy of the lithospheric mantle above a mantle plume: insights from the Labait volcano xenoliths (Tanzania): Earth and Planetary Science Letters, v. 232, p. 295-314.

Weinstein, A., Ebinger, C. J., Oliva, S. J., Roecker, S., Tiberi, C., Aman, M., Lambert, C., Witkin, E., Albaric, J., Gautier, S., Muzuka, A., Mulibo, G., Kianji, G., Hadfield, R., Illsley-Kemp, F., Msabi, M., Ferdinand-Wambura, R., Perrot, S., Muirhead, J. D., Rodzianko, A., Fischer, T., in review, Volatile-assisted rifting in East Africa: Seismicity of the Magadi-Natron-Manyara basins, Africa: in review in Geochemistry, Geophysics, Geosystems.