# Stratigraphic framework of pulsed volcanism in the central Ethiopian Plateau

Dereje Ayalew, Faysel Sefa Abdu, Raphaël Pik, Bekele Abebe, Worash Getaneh, Balemwal Atnafu

► To cite this version:

Dereje Ayalew, Faysel Sefa Abdu, Raphaël Pik, Bekele Abebe, Worash Getaneh, et al. Stratigraphic framework of pulsed volcanism in the central Ethiopian Plateau. Bulletin of Volcanology, 2024, 86 (4), pp.article id.33. 10.1007/s00445-024-01726-8. insu-04783362

HAL Id: insu-04783362

https://insu.hal.science/insu-04783362v1

Submitted on 30 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Stratigraphic framework of pulsed volcanism in central Ethiopian plateau

Dereje Ayalew

dereayal@yahoo.com

Addis Ababa University College of Law and Governance Studies

Faysel Sefa

Addis Ababa University

Raphaël Pik

CRPG CNRS-UL: Centre de Recherches Petrographiques et Geochimiques

Bekele Abebe

Addis Ababa University

Worash Getaneh

Addis Ababa University

Balemwal Atnafu

Addis Ababa University

## Research Article

**Keywords:** Basalt, Ethiopia, formation, magmatic pulse, plateau, rhyolite, stratigraphy

**Posted Date:** July 6th, 2023

**DOI:** https://doi.org/10.21203/rs.3.rs-3090728/v1

**License:** © This work is licensed under a Creative Commons Attribution 4.0 International License.
Read Full License

**Version of Record:** A version of this preprint was published at Bulletin of Volcanology on March 11th, 2024. See the published version at https://doi.org/10.1007/s00445-024-01726-8.
**Abstract**

This paper presents new stratigraphic correlation and naming for younger basalt lavas (<24 Ma) and enclosed rhyolites near the intersection of the southern Red Sea rift and the Ethiopian rift valley, based upon three volcano-stratigraphic sections located over about 30 kilometres, and they bracket the timing of eruption of the upper sequence with existing geochronology. The correlations are based on the petrography (thin sections) of ca. 32 samples taken from these three sections, combined with outcrop description. The stratigraphy is proposed as 5 new formations, partially placed into an existing stratigraphic framework, based on the dominant unit-basalt or rhyolite in the package of sequences. These new stratigraphic columns also include non-volcanic units (interlayered sedimentary strata, paleosols and weathering surfaces), which are not usually described in the volcanic pile, showing that there was more complexity in this succession than previously known. This new stratigraphic framework will encourage more geochronology and geochemistry studies that are sorely needed.

**Introduction**

Continental flood basalt (CFB) provinces form the largest volcanic features on Earth, typically both form in extensional environments and predate extension (Richards et al. 1989; Courtillot et al. 1999). There is a co-dependence between flood basalts and continental rifting but timing is complex. For instance, the Keweenawan LIP (large igneous province) of Mid-Continent rift was erupted into a developing rift (Elling et al. 2020), while the Deccan traps likely occurred prior to rifting (Sprain et al. 2019). On the other hand the flood basalts in East Africa seem to have erupted, at least partially, during localized phases of extension (Rooney 2017). Most CFB provinces appear to be characterized by vast thickness of laterally extensive lava flows of relatively evolved tholeiitic basalts (Hooper and Hawkesworth 1993; Hooper et al. 2007; Rooney 2017). High resolution geochronology has revealed that the peak of magmatism is quite a bit shorter than even <1 Ma (Bryan et al. 2010; Self et al. 2014).

Silicic volcanics are an important component of many flood basalt provinces (Halder et al. 2021). They usually occur in the upper parts of the flood basalt sequence. There is a problem or question (e.g., timing, magma plumbing system and origin) surrounding silicic volcanism in flood basalt provinces.

Stratigraphic studies are critical for understanding flood basalt provinces. Such studies highlight petrographic and geochemical variation in stratigraphy related to the magmatic differentiation system (Krans et al. 2018; Davis et al, 2021; Streck et al. 2023). Hence stratigraphic study is needed to document the temporal and spatial changes in magma chemistry and mineralogy, which are related to magma evolution.

Current review papers divide the Ethiopian flood basalts, related to plume impingement, into two sequences: the Oligocene traps phase (31~28Ma; Rooney 2017) and the Eocene flood basalt province/LIP (Rooney 2017; Steiner et al. 2021). The study area (Fig. 1) is unrelated to the early Oligocene traps episode. Instead, it corresponds to younger lavas (<24 Ma), referred to as the early Miocene
resurgence phase (Rooney 2017), in central Ethiopian plateau. Stratigraphic studies are lacking in the area, meaning there was not a detailed stratigraphic section. Much of the volcanic succession in this area is preserved due to the relatively young age compared to other flood basalts. Such a thing makes these lavas very useful in that much or all of their history is recorded. This study employs field observation and petrography in the context of stratigraphy, near the intersection of the southern Red Sea rift and the Ethiopian rift valley, with the aim to describe/define the formations and their stratigraphic units to help interpret the overall thickness, timescales, eruption evolution from basalt to silicic volcanism and pulsed eruption. It fully addresses the current state of knowledge regarding pulses of magmatism in East Africa (Rooney 2017; Krans et al. 2018; Rooney 2020a-d; Steiner et al. 2021).

## Geologic background

Recent reviews of magmatism in East Africa (Rooney, 2017; Krans et al. 2018; 2020a; 2020b; 2020c; 2020d; Steiner et al. 2021) detail pulses of magmatism in East Africa and contextualized those pulses into geologic process active during eruption. These review papers show that the flood basalt event (~45-28 Ma) is different tectonically and magmatically from the younger lavas (<24 Ma). Magmatism of this age is more closely linked to extensional events than to the impingement of the mantle plume that generated the flood basalts (while there is a plume component in the magma source).

Current review papers divide the Ethiopian flood basalts into two sequences: the Oligocene Traps phase (31~28 Ma; Rooney 2017) and the Eocene flood basalt province/LIP (Rooney 2017; Steiner et al. 2021). The oldest flood basalt lavas are found in southern Ethiopia (Davidson and Rex 1980; Ebinger et al. 1993; George et al. 1998). The volume of flood lavas erupted in southern Ethiopia is significant (Ebinger et al. 1993; George and Rogers 2002; Steiner et al. 2021). There are also several thick silicic volcanics in the flood basalts of Ethiopia and northern Kenya (Ebinger et al. 1993; Steiner et al. 2021). Specifically, the regional Jimma volcanic unit, the Amaro tuff, and the Arba Minch tuff are some of the documented examples. The Oligocene traps phase was followed by volumetrically less significant lavas manifested as shield volcanoes or isolated flows, referred to as the early Miocene resurgence phase (Rooney 2017).

Rooney (2017; 2020a) and citations therein lay out that magma pulses in East Africa occur due to the interplay of (1) impingement of a mantle plume (the Oligocene flood basalts in northwest and central Ethiopia, and the Eocene flood basalts in southwestern Ethiopia, northern Kenya, and potentially in South Sudan); and (2) the interplay between extension, decompression melting, and elevated mantle temperatures beneath East Africa (the younger lavas <24 Ma).

The age of the lavas presented (20-10 Ma) would suggest that they are intimately associated with periods of extension. Though these younger lavas may take on the form of stratiform sheet flows, similar to the Oligocene traps phase from a morphological sense, their magmatic origin differs from the initial flood basalt pulses. Consequently, they are treated as a different sequence than those of the Oligocene flood basalts (Oligocene traps phase). Nevertheless, the main flood basalts nowhere show contact with the younger lavas in the area.
# Methods

Traverse routes, made along road cuttings, were selected with the purpose of identifying all mappable units and constructing the local and composite stratigraphy. Our method entailed detailed field logging and regular sampling through the succession at each selected locality, followed by a geologic mapping. Three stratigraphic sections were examined, measured and sampled. We choose traverses that are bracketed by existing geochronology, but the bottom section of the sequences is undated. The 3 stratigraphic successions examined identify as mixed and layered sections. Basalts were sepaared out from rhyolites based upon the interpretation of hiatuses (non-volcanic units) such as sediments, or soils, or weathering surfaces. These non-volcanic units definitively mean that volcanism was not active at that spot for a certain amount of time. The larger units found below and above each hiatus are proposed as formations. The thickness of the volcanic pile was determined by measuring elevation variations recorded with a handheld GPS.

We sampled all identified flows (32 unique samples). Individual flows were identified in the field based upon lava texture (e.g., aphyric vs. porphyritic and modal mineralogy), used as a classification system (Renith 2014; Krans et al. 2018). The detailed petrography is given for each sample (Supp. Material). The details of the petrography for larger correlating units together with field observation are provided in the description of each formation. Figure 4 illustrates detailed sections that describe the stratigraphic units and locations of the studied samples. This eventually indicates the units and horizons (base, middle, top) of the samples where they are from.

# Results

The volcanic succession for the central Ethiopian plateau can be broadly separated into lower fissure-fed basalts and rhyolites overlain by upper point-sourced shield basalts. We used the formal stratigraphic term "Formation" to classify the different unit types which are mappable at this scale (Fig. 2). All units are formally named using a geographic name, followed by a lithic designation which is more descriptive than the generic term "Formation" (Reidel et al. 2013). The geographic designators were derived from the type locality/area where a complete sequence of each unit is exposed. Five formations are identified from base to top: Zaro Basalt, Gawna Rhyolite, Sasit Basalt, Alaji-Molale Rhyolite and Tarmaber-Guassa Basalt (Figs. 3 and 4), which are partially placed into an existing stratigraphic framework (Zanettin and Justin-Visentin 1973; Justin-Visentin et al. 1974). The formation ascriptions given in Table 1 are taken from the results of traverse and areal/local field mapping. The stratigraphic scheme constructed from this work is summarized as follows.

## Zaro Basalt

The oldest in the stratigraphic sequence is here named the Zaro Basalt. The Zaro Basalt lies in contact with the underlying Mesozoic sandstone (Debre Libanos Formation, Asefa, 1991). The unconformity is described by an uneven (undulating) surface as indicated by frequent lava ponding. The Zaro Basalt
occupies the lowest elevations (<2000 m above mean sea level) in the westward draining river valleys in the study area (Fig. 2). It has a total thickness of about 400 m (Fig. 4).

The Zaro Basalt forms subdued relief and the flows are poorly exposed. It consists dominantly of dark, basalt lava flows with a thin intercalation of rhyolite ignimbrite towards the top. Basalt flow thicknesses commonly vary between 10 and 40 m. The flow top of the lavas is typically scoriaceous. Individual flows are commonly columnar (Fig. 5). The flows are intensely fractured, typically with alteration surfaces.

The basalt flows exhibit variable textural features including aphyric and porphyritic varieties. The porphyritic varieties contain phenocrysts (up to 10 vol. %) of plagioclase (up to 2.8 mm in size) and olivine (1.5 mm in length) set in a microcrystalline groundmass consisting of plagioclase, olivine and opaque minerals (Fig. 6, Supp. Material). Olivine is often altered to iddingsite. The Zaro Basalt flows are variably altered and contain amygdales of quartz, calcite, chlorite and zeolite that are replacing phenocrysts and groundmass.

There is an intercalated sedimentary layer, dominantly sandstone with variable grain size and degree of weathering, in the middle of the succession (Fig. 3). A similar intercalation within basalts has been observed in the Jema River valley, farther west of the study area (Abbate et al. 2014; Passey et al. 2021). The presence of inter-lava sediments in these two localities and many other places suggests a considerable volcanic quiescence in the middle of the eruption of the Zaro Basalt. This also means that some areas were at lower elevation than others.

## Gawna Rhyolite

The early Zaro Basalt is unconformably overlain by thick rhyolite lavas (Gawna Rhyolite) in the western side of the mapped area. This is the first reported occurrence of glassy rhyolite lava flows associated with the Ethiopian basalt pile. These silicic lavas appear to find no equivalent in the northern Ethiopian plateau. The Gawna Rhyolite forms prominent vertical cliffs (Fig. 7) and comprises a pile up to 170 m thick (Fig. 4).

The rhyolite lava sequence displays variable textural features (Fig. 8). In the lower and middle portions, the flows exhibit vertically oriented, regular colonnade-like structures (Fig. 8a). In contrast, towards the top the lava flows consist of randomly oriented joints forming entablature (Fig. 8b). The Gawna Rhyolite lavas exhibit flow banding and flow folding texture. Rhyolite flow thicknesses vary between 5 and 50 m. The rhyolite lavas frequently show perlitic cracks, implying hydration of obsidian.

The Gawna Rhyolite is dark in color and is typically porphyritic containing quartz and alkali feldspar phenocrysts set in a glassy matrix (Fig. 6, Supp. Material). Towards the top, the Gawna Rhyolite is intensely altered and contains brown opals, white agate and green jasper/chert. The glassy (vitrophyre) texture of the Gawna Rhyolite possibly relates to rapid loss of dissolved gas by viscous lava fluxing at a higher rate to a cooler surface air condition (Hefferan and O'Brien 2010). The Gawna Rhyolite extends
laterally for at least tens of kilometres, indicating that it was fluidal and mobile due to perhaps reduced viscosity as a result of compositional factor.

## Sasit Basalt

Whole rock Rb-Sr isochron on intercalated rhyolites from this formation gave an age of ca. 20.7 Ma (Ayalew 2011). The Sasit Basalt is composed of dense, dark basalt flows with numerous thin intercalations of rhyolite ignimbrite, especially towards the upper portion (Fig. 4). It is preserved over much of the plateau, making the Sasit Basalt the most voluminous of this region. The basalt succession reaches to 440 m thick (Fig. 4).

The Sasit Basalt comprises of at least nine basaltic flows that extend subhorizontally for at least tens or more kilometers (Fig. 9). The flows are between 15 and 50 m thick. Individual flows are very extensive, lithologically uniform, consistent with sheet flow geometry. The flow top of Sasit Basalt commonly consists of a vesicular to scoriaceous basalt or a ruby to brecciated basalt (Fig. 10). It typically occupies up to 50% of the entire thickness of a flow. Flow top breccias typically consists of angular, vesicular fragments of basaltic rubble that lie above a zone of nonfragmented, vuggy basalt displaying entablature-like cooling joints with fan-shaped base (Fig. 10). It is likely created by autobrecciation process. The flow bottom of the flows typically consists of sparsely vesicular, fine-grained basalt, suggesting lava flowage on relatively dry ground (Fig. 10). Although the great majority of the Sasit Basalt flows are sheet flow, there are a few compound flows with numerous small distinct lobes (Supp. Material). Individual lobes are affected by shpheroidal weathering.

The Sasit Basalt flows are variably porphyritic (up to 12 vol. %) with phenocrysts of dominantly plagioclase (up to 3 mm in size) with subordinate olivine (about 1.5 mm long) and pyroxene (<1 mm in length) set in a microcrystalline groundmass comprising of lath-shaped plagioclase microlites, olivine, clinopyroxene and opaque minerals (Fig. 6, Supp. Material). The phenocrysts often form glomerophyric aggregates. The basalts are locally altered and that an agate-bearing basaltic layer has been identified at an elevation of ca. 2585 m above base level.

It is important to note that the Sasit Basalt includes numerous thin (up to 20 m thick) intercalations of silicic tuff and ignimbrite (Fig. 4), but are not mappable at this scale. This unit is conspicuous for containing carbonized tree trunks/branches (Fig. 11), indicating wet environment during the deposition of this unit. Towards the top, there is a thin rhyolite level which bears precious opal (Ayalew et al. 2020). Inter-lava sediments (conglomerate, sandstone and mudstone) are common between basalt lava flows (Supp. Material). Their presence implies there was considerable volcanic quiescence between the eruptions of successive lava flows.

## Alaji-Molale Rhyolite

The Alaji-Molale Rhyolite is dated at ca. 15.4 Ma (K-Ar dating, Zanettin et al., 1974 and Rb-Sr isochron, Ayalew et al. 2002). It comprises of silicic (rhyolite) ignimbrite and rests directly on the Sasit Basalt. The
Alaji-MolaleRhyolite is widely distributed over the plateau, forming the flat terminal surface with typical pyramidal topography (Supp. Material). The thickness of the rhyolite succession attains about 120 m (Fig. 4) which pinches out westward away from the Afar rift margin. The Alaji-MolaleRhyolite serves as a marker horizon to separates the sheet flow sequence from the overlying shield-like pile.

The Alaji-MolaleRhyolite is dominated by pyroclastic deposits with a thin intercalation of basaltic lava flow. Individual flows are between 10 and 35 m thick. The ignimbrite sequence maintains a similar thickness over kilometers and includes both welded and unwelded units. The welded varieties show a well-developed eutaxitic texture, implying their emplacement at relatively high temperature. Individual prominent flows are identified over a distance of a few tens of kilometers which form typical landscape elements such as steep walls and pyramids (Supp. Material).

The Alaji-MolaleRhyolite consists of an alternance of alkaline to peralkaline large ignimbrites and tuffaceous levels. It is variably porphyritic, containing phenocrysts of alkali feldspar and quartz enclosed within a glassy matrix (Fig. 6, Supp. Material). The peralkaline varieties include comendites; pantellerites are uncommon (Ayalew et al., 2002). The interbedded basalt flow is lithologically indistinguishable from underlying Sasit Basalt.

## Tarmaber-Guassa Basalt

The Tarmaber-Guassa Basalt is dated at ca. 10.4 Ma (K-Ar technique, Chernet et al. 1998), terminating the central Ethiopian magmatic event. It overlies unconformably the Alaji-MolaleRhyolite ignimbrites, separated by about 30 m thick sediments with variegated conglomerates, sandstones and mudstones. The Tarmaber-Guassa Basalt crops out close to the rift margin in the more elevated places (their elevation can reach to 3315 m above mean sea level) of the mapped area. It locally fills ancient erosional channels. The thickness of this formation is highly variable, in places reaching up to 410 m (Fig. 4).

The Tarmaber-Guassa Basalt forms still-preserved edifices (shield volcanoes), made up of more localized, lenticular basalt flows with scoriaceous flow top. The flows are frequently separated by red paleosol. They dip 3-5° away from the summit. The lavas are commonly scoriaceous and that vesicles are often filled with amygdales. Pyroclastic layers are abundant in some sequences. Although basalt flows dominate the sequence, differentiated products are also frequent.

The Tarmaber-Guassa Basalt flows are distinctive by being porphyritic (phenocrysts comprise up to half, sometimes more, of the rock volume); plagioclase, pyroxene and olivine phenocrysts characterize the flows (Fig. 6, Supp. Material). Nevertheless, there are also aphyric varieties. The character of the basalts is decidedly alkaline. More evolved products include dominante plagioclase- or olivine-phyric hawaiite/mugearite and subordinate aphyric trachyte and tephrophonolite.

# Discussion

## Overall thickness of the stratiform sheet flows
The basal flow of the sheet flowsequence is exposed in Jema gorge along the western side of the mapped area (Fig. 2). In this locality, the contact between the Mesozoic sandstone and overlying basalt is found at an altitude of ca. 1620 m above mean sea level (Fig. 4). The top of the stratiform sheet flows, away from plateau-Afar escarpment to minimize the effect at the rift margin where rift flank uplift is superposed on broader dynamic uplift, resides at an elevation of ca. 2910 m above base level (Fig. 4). Considering these elevations ca. 1620 m as the base and ca. 2910 m as the top of the sheet flowsequence, we estimated the present, preserved thickness of the basalt and enclosed rhyolite pile to be 1290 m.

Estimates of the thickness of the volcanic sequence are always subject to a certain degree of uncertainty due to the effects of erosion. It is worth mentioning that the calculated thickness should be considered as a minimum estimate because the calculation did not take into account the original stack before erosion. Hence, the total thickness of the basalt and rhyolite pile can exceed 1290 m and a somewhat larger figure applies. It is apparent that the sheet flow sequence in the area is far thinner than in the Oligocene flood basalts in the north (~2100 m thick, Ayalew et al. 2021). Nevertheless, it built up to the same order of thickness as other major flood basalt provinces (e.g., the Columbia River Basalt thickens 1500 m, Reidel et al. 2013).

## Frequency of silicic volcanism

Silicic volcanics are an important component of many flood basalt provinces (Halder et al. 2021). Thick silicic volcanics exist both in the flood basalts of Ethiopia and northern Kenya, and subsequent younger lavas (Rooney 2017; Steiner et al. 2021). They occur in the volcanic stratigraphy as ash layers between basaltic lava flows or as individual eruptive units, comparable in volume to individual mafic lava units. Cooling units for these felsic ignimbrites and lavas extend laterally up to several tens of kilometers, indicating voluminous fissure eruptions.

Widespread silicic volcanics occur in the study area. They are prolifically intercalated with and superimposed on the basalt flows, indicating periodic silicic eruption. This is supported by the first time thick rhyolite lava flows (Gawna Rhyolite) reported from the basalt succession of Ethiopia (Figs. 2 and 4). Neglecting the interlayered rhyolites, the larger felsic eruptives form about 30% by volume of the total volcanic pile (Fig. 4). This yields their value in labeling and identifying particular formations (Figs. 3 and 6). The presence at certain horizons and localities of prominent silicic ignimbrites and lavas capping the basalt pile leads to various petrogenetic models of silicic magma generation.

The dataset presented (20-10 Ma) suggest that these younger lavas are related to extension (Rooney 2017). This allows for the shallow ponding of mafic magma, and subsequent silicic magma generation possibly by fractional crystallization and ultimately periodic silicic eruptions. This is supported by the presence of extensive rhyolites in the lava sequence that illustrates the time-progressive evolution of magma. This suggests that the frequency and volumetric importance of the silicic volcanism increased progressively in the younger lavas (<24 Ma).
## Episodicity of Ethiopian volcanism

The current state of knowledge regarding magmatism in East Africa (Rooney, 2017; Krans et al. 2018; Rooney 2020a; 2020b; 2020c; 2020d; Steiner et al. 2021) detail pulses of magmatism in East Africa where the main flood basalt events (~45-28 Ma) end and subsequent shield building takes place afterward followed by rift-related volcanism. The age of the lavas presented (20-10 Ma) would suggest that they are ultimately linked with periods of extension. Though these younger lavas may take on the form of stratiform sheet flows, their magmatic origin (petrologically and geochemically) differ from the initial flood basalt pulses (Rooney 2017; 2020a and citations therein). As a result, they are treated as a different sequence than those of the Oligocene flood basalts (Oligocene traps phase).

The 3 stratigraphic sections logged (Fig. 4) bracket the timing of eruption of the upper sequence with existing geochronology, but the bottom section remains undated. Whole rock Rb-Sr isochron on intercalated rhyolites within the Sasit Basalt yield an age of ca. 20.7 Ma (Ayalew 2011), permitting an earlier start of volcanism. These geochronologic dates suggest that volcanism began in the area most likely in early Miocene, flowing the massive outpouring of flood basalts and enclosed rhyolites (the Oligocene traps phase 31~28 Ma).

The volcanic succession in the area is characterized by significant basalt flows in pulses, each capped by silicic volcanics. Long intervals between successive flows within or between individual formations are attested to by the common development of inter-lava sediments, paleosols and alteration surfaces (Supp. Material), implying long time spans of eruption. This Miocene activity typically defines as a significant basaltic event followed by a period of quiescence or silicic activity, consistent with the breakdown and merging of regional units in Rooney (2020a). Yet, this long period, suggesting pulsed magma eruption, has resulted in the Ethiopian pile building up to the same order of thickness as other major flood basalt provinces (e.g., central Atlantic magmatic province, Marzoli et al. 1999).

## Conclusion

Stratigraphy study, together with any published ages, on rift-related lavas (which may take on the form of stratiform sheet flows) in central Ethiopian plateau indicates that important episodes of magmatism have been active in the area, characterized by slow widespread basalts in pulses, each capped by silicic volcanic (ignimbrites and lavas). Five formations, separated by non-volcanic units (e.g., sediments, paleosols and weathering surfaces) are identified and mapped. These formations, partially placed into an existing stratigraphic framework, are named from base to top: Zaro Basalt, Gawna Rhyolite, Sasit Basalt, Alaji-Molale Rhyolite and Tarmaber-Guassa Basalt.

## Declarations

### Acknowledgements
The authors would like to thank Addis Ababa University ("thematic research" grant scheme) and Amhara Mines Resource Development Promotion Agency (geologic mapping of five Woredas in north Shoa province project) for their financial support. We are indebted to anonymous journal reviewers and Associate Editor Loïc Vanderkluysen who made useful comments on the manuscript.

## References

1. Abbate E, Bruni P, Ferretti MP, Delmer C, Laurenzi MA, Hagos M, Bedri O, Rook L, Sagri M, Libsekal Y (2014) The East Africa Oligocene intertrappean beds: regional distribution, depositional environments and Afro/Arabian mammal dispersals. *J Afr Earth Sci* 99: 463-489.

2. Assefa G (1991) Lithostratigraphy and environment of deposition of the late Jurassic-early Cretaceous sequence of the central part of northwestern Plateau. *N Jb Geol Paläont Abh* 182: 255-284; Stuttgart.

3. Ayalew D (2011) The relations between felsic and mafic volcanic rocks in continental flood basalts of Ethiopia: implication for the thermal weakening of the crust. *Geol Soc Lon Spec Publ* 357: 253-264.

4. Ayalew D, Barbey P, Marty B, Reisberg L, Yirgu G, Pik R (2002) Source, genesis, and timing of giant ignimbrite deposits associated with Ethiopian continental flood basalts. *Geochim Cosmochim Acta* 66(8):1429-1448.

5. Ayalew D, Getaneh W, Pik R, Atnafu B, Zemelak A, Belay E (2021) Stratigraphic framework of the northeastern part of the Ethiopian flood basalt province. *Bull Volcanol* 83. 10.1007/s00445-021-01482-z.

6. Ayalew D, Pik R, Gibson S, Yirgu G, Ali S, Assefa D (2020) Pedogenic origin of Mezezo opal hosted in Ethiopian Miocene rhyolites. *Can Mineral* 58(2): 231-246.

7. Bryan SE, Ukstins Peate I, Peate DW, Self S, Jerram DA, Mawby MR, Marsh JS, Miller JA (2010) The largest volcanic eruptions on Earth. *Earth Sci Rev* 102(3-4): 207-229.

8. Chernet T, Hart WK, Aronson JL, Walter RC (1998) New age constraints on the timing of volcanism and tectonism in the northern main Ethiopian rift-southern Afar transition zone (Ethiopia). *J Volcanol Geotherm Res* 80: 267-280

9. Coulié E, Quidelleur X, Gillot P-Y, Courtillot V, Lefèvre J-C, Chiesa S (2003) Comparative K–Ar and Ar/Ar dating of Ethiopian and Yemenite Oligocene volcanism: implications for timing and duration of the Ethiopian traps. *Earth Planet Sci Lett* 206:477-492.

10. Courtillot V, Jaupart C, Manighetti I, Tapponnier P, Besse J (1999) On causal links between flood basalts and continental breakup. *Earth Planet Sci Lett* 166: 177-195.

11. Davidson A, Rex DC (1980) Age of volcanism and rifting in southwestern Ethiopia. *Nature* 283: 657-658.

12. Davis WR, Collins MA, Rooney TO, Brown EL, Stein CA, Stein S, Moucha R (2021) Geochemical, petrographic, and stratigraphic analyses of the Portage Lake volcanics of the Keweenawan CFBP:
implications for the evolution of main stage volcanism in continental flood basalt provinces. Geol Soc Spec Publ 518: 67-100.

13. Ebinger CJ, Yemane T, WoldeGabriel G, Aronson JL, Walter RC (1993) Late Eocene-Recent volcanism and faulting in the southern main Ethiopian rift. J Geol Soc 150: 99-108.

14. Elling RP, Stein S, Stein CA, Keller GR (2020) Tectonic implications of the gravity signatures of the Midcontinent rift and Grenville Front. Tectonophysics 778: 228-369.

15. George, R., Rogers, N., 2002. Plume dynamics beneath the African plate inferred from the geochemistry of the Tertiary basalts of southern Ethiopia. Contrib Mineral Petrol 144: 286-304.

16. George R, Rogers N, Kelly S (1998) Earliest magmatism in Ethiopia: evidence for two mantle plumes in one flood basalt province. Geology 26: 923-926.

17. Halder M, Paul D, SensarmaS (2021) Rhyolites in continental mafic large igneous provinces: petrology, geochemistry and petrogenesis. Geosci Front 12: 53-80.

18. Hefferan K, O'Brien J (2010) Earth materials. Wiley-Blackwell pp. 608.

19. Hofmann C, Courtillot V, Féraud G, Rochette P, Yirgu G, Ketefo E, Pik R (1997) Timing of the Ethiopian flood basalt event and implications for plume birth and global change. Nature 389: 838-841.

20. Hooper PR, Camp VE, Reidel SP, Ross ME (2007) The origin of the Columbia River flood basalt province: plume versus nonplume models. In: Foulger GR, Jurdy DM (eds) Plates, plumes, and planetary processes. Geol Soc Am Spec Pap 430: 635-668.

21. Hooper PR, Hawkesworth CJ (1993) Isotopic and geochemical constraints on the origin and evolution of the Columbia River Basalt. J Petrol 34: 1203-1246.

22. Justin-Visentin E, Nicoletti M, Tolomeo L, Zanettin B (1974) Miocene and Pliocene volcanic rocks of the Addis Ababa - Debra Berhan area (Ethiopia). Bull Volcanol 38: 237-253.

23. Kieffer B, Arndt N, Lapierre H, Bastien F, Bosch D, Pecher A, Yirgu G, Ayalew D, Weis D, Jerram DA, Keller F, Meugniot C (2004) Flood and shield basalts from Ethiopia: magmas from the African Superswell. J Petrol 45:793-834.

24. Krans SR, Rooney TO, Kappelman J, Yirgu G, Ayalew D (2018) From initiation to termination: a petrostratigraphic tour of the Ethiopian low-Ti flood basalt province. Contrib Mineral Petrol 173(5): 37.

25. Marzoli A, Renne PR, Piccirillo EM, Ernesto M, Bellieni G, De Min A (1999) Extensive 200 million-year-old continental flood basalts of the central Atlantic magmatic province. Science 284 (5414): 616-618.

26. Mohr P, Zanettin B (1988) The Ethiopian food basalt province. In: Macdougall, J.D. (ed) Continental flood basalts. Kluwer Academic Publishers, pp 63-110.

27. Nelson WR, Hanan BB, Graham DW, Shirey SB, Yirgu G, Ayalew D, Furman T (2019) Distinguishing plume and metasomatized lithospheric mantle contributions to post-flood basalt volcanism on the southeastern Ethiopian plateau. J Petrol 60(5):1063-1094.
28. Passey SR, McLean C, Ayalew D (2021) Controls on sediment distribution in a volcanically-affected basin: insights from the Ethiopian flood basalt province. Geol Soc Lon Spec Publ DOI: 10.1144/SP520-2021-70.

29. Reidel SP, Camp VE, Tolan TL, Martin BS (2013) The Columbia River flood basalt province: stratigraphy, areal extent, volume, and physical volcanology. Geol Soc Am Spec Pap 497: 1-43.

30. Renjith ML (2014) Micro-textures in plagioclase from 1994-1995 eruption, Barren Island volcano: evidence of dynamic magma plumbing system in the Andaman subduction zone. Geosci Front 5: 113-126.

31. Richards M, Duncan R, Courtillot V (1989) Flood basalts and hot-spot tracks: plume heads and tails. Science 246: 103-107.

32. Rooney TO (2017) The Cenozoic magmatism of East-Africa: part I- flood basalts and pulsed magmatism. Lithos 286: 264-301.

33. Rooney TO (2020a) The Cenozoic magmatism of East Africa: Part II – rifting of the mobile belt. Lithos 360-361: 105291.

34. Rooney TO (2020b) The Cenozoic magmatism of East Africa: Part III – rifting of the craton. Lithos 105390.

35. Rooney TO (2020c) The Cenozoic magmatism of East Africa: Part IV- the terminal stages of rifting preserved in the Northern East African rift system. Lithos 360-361: 105381.

36. Rooney TO (2020d) The Cenozoic magmatism of East Africa: part V– magma sources and processes in the East African Rift. Lithos 360: 105296.

37. Self S, Schmidt A, Mather TA (2014) Emplacement characteristics, time scales, and volcanic gas release rates of continental flood basalt eruptions on Earth. Geol Soc Am Spec Paper 505.

38. Sprain CJ, Renne PR, Vanderkluysen L, Pande K, Self S, Mittal T (2019) The eruptive tempo of Deccan volcanism in relation to the Cretaceous-Paleogene boundary. Science 363: 866-870.

39. Steiner RA, Rooney TO, Girard G, Rogers N, Ebinger C, Peterson L, Phillips R (2022) Initial Cenozoic magmatic activity in East Africa: new geochemical constraints on magma distribution within the Eocene continental flood basalt province. Geol Soc Spec Publ 518(1): 435-465.

40. Streck MJ, Swenton VM, McIntosh W, Ferns ML, Heizler M (2023) Columbia River Rhyolites: age-distribution patterns and their implications for arrival, location, and dispersion of continental flood basalt magmas in the Crust. Geosciences 13(2): 46; https://doi.org/10.3390/geosciences13020046

41. Ukstins IA, Renne PR, Wolfenden E, Baker J, Ayalew D, Menzies M (2002) Matching conjugate rifted margins: $^{40}$Ar/$^{39}$Ar chrono-stratigraphy of pre- and syn-rift bimodal flood volcanism in Ethiopia and Yemen. Earth Planet Sci Lett 198: 289-306.

42. Zanettin B, Gregnanin A, Justin-Visentin E, Nicoletti M, Petrucciani C, Piccirillo EM, Tolomeo L (1974) Migration of the Oligocene-Miocene ignimbritic volcanism in the central Ethiopian plateau. Neues Jahrb Geol PM 9: 567-574.
43. Zanettin B, Justin-Visentin E (1973) Serie di vulcanite etiopiche I: la serie dell'altipiano etiopico centro-orientale. Boll Soc Geol It 92:313-327.

44. Zanettin B, Justin-Visentin E Piccirillo EM (1978) Volcanic succession, tectonics and magmatology in central Ethiopia. Atti Mem Accad Patavina Sc 90: 5-19.

Table 1

<table><caption>Table 1 Simplified stratigraphic columns for late Oligocene-Miocene volcanic successions on the central Ethiopia plateau.</caption><thead><tr><th>Locality</th><th>Profile coordinate</th><th>Elevation (m, top)</th><th>Formation</th><th>Dominant lithology</th><th>Thickness (m)</th><th>Age (Ma)</th></tr></thead><tbody><tr><td rowspan="4">Seladingay</td><td rowspan="4">9°57'17''N<br>39°37'49''E</td><td rowspan="4">3325</td><td>Tarmaber</td><td>basalt</td><td>410</td><td>10.4</td></tr><tr><td>Seladingay</td><td>rhy. ign.</td><td>120</td><td>15.4</td></tr><tr><td>Sasit</td><td>basalt</td><td>400</td><td>20.7</td></tr><tr><td>Gawna</td><td>rhy. lava</td><td>200</td><td></td></tr><tr><td rowspan="4">Sasit</td><td rowspan="4">9°55'29''N<br>39°29'44''E</td><td rowspan="4">2780</td><td>Seladingay</td><td>rhy. ign.</td><td>200</td><td></td></tr><tr><td>Sasit</td><td>basalt</td><td>460</td><td></td></tr><tr><td>Gawna</td><td>rhy. lava</td><td>300</td><td></td></tr><tr><td>Zaro</td><td>basalt</td><td>180</td><td></td></tr><tr><td rowspan="3">Gawna</td><td rowspan="3">9°55'41''N<br>39°22'39''E</td><td rowspan="3">2680</td><td>Sasit</td><td>basalt</td><td>460</td><td></td></tr><tr><td>Gawna</td><td>rhy. lava</td><td>300</td><td></td></tr><tr><td>Zaro</td><td>basalt</td><td>400</td><td></td></tr></tbody></table>

Abbreviations: rhy. Rhyolite, ign. Ignimbrite.

Age sources: Chernet et al. (1998), Coulié et al. (2003), Ayalew (2011).

Figures
<figure><img src="image_2.png" /><figcaption>Figure 1</figcaption></figure>

Topographic relief map showing the extent of flood basalts in Ethiopia and Yemen (modified from Rooney 2017). Numbers indicate ages of the earliest flood basalt in different parts of the Ethiopian plateau, highlighting the spatial diachroneity of magmatism. Colors represent different volcanism as shown in the legend. Inset rectangle (red) denotes location of the studied area which is part of the central Ethiopian plateau. Age source: Davidson and Rex, 1980; Ebinger et al. 1993; Hofmann et al. 1997; George et al. 1998; Ukstins et al. 2002; Coulie et al. 2003; Kieffer et al. 2004; Nelson et al. 2019.
<figure><img src="image_3.png" /><figcaption>Figure 2</figcaption></figure>

Geologic map of the investigated area. Google Earth imagery as a base-map has been used to construct the geologic map. Stratigraphic sections and traverse routes are also contributed to assemble the geologic map. Locations of the type localities, stratigraphic successions and field photos are shown.

<figure><img src="image_4.png" /><figcaption>Figure 3</figcaption></figure>

View of the lower three formations along tributaries of Jema River.
<figure><img src="image_5.png" /><figcaption>Figure 4</figcaption></figure>

Three measured stratigraphic logs on the study area.

<figure><img src="image_6.png" /></figure>
<caption>Figure 5</caption>

View of Zaro Basalt lying in contact with Mesozoic sandstone taken from Gawna Section. Note that the Zaro Basalt flows exhibit columnar cooling joints around the contact.

<figure><img src="image_7.png" /><figcaption>Figure 6</figcaption></figure>

Composite stratigraphy of the volcanic succession on the study area with representative petrography from each formation.
<figure><img src="image_8.png" /><figcaption>Figure 7</figcaption></figure>

View of Gawna Rhyolite lava flows forming steep cliff taken from Gawna section.
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8**</figcaption></figure>

Photo of contraction/cooling joints showing colonnade in the lower (a) and entablature in the upper (b) portions of the Gawna Rhyolite lavas. Note person for scale representing 1.6 m in length/height.
<figure><img src="image_11.png" /><figcaption>Figure 9</figcaption></figure>

View of Sasit Basalt comprising of numerous thin, laterally extensive sheet flows.
<figure><img src="image_12.png" /><figcaption>**Figure 10**</figcaption></figure>

Photo of lava flow geometry from Sasit Basalt illustrating flow top breccias above entablature-like joints and a sparsely vesicular to scoriaceous, fine-grained basalt flow top of an underlying flow. Note sub-horizontal segmentation patterns, reflecting inward progressive propagation of joint faces. Note person for scale representing 1.6 m in length/height.
<figure><img src="image_13.png" /><figcaption>**Figure 11**</figcaption></figure>

Fossilized woods (tree logs) within interlayered rhyolite from the lower part of the Sasit Basalt along Seladingay section. Note persons for scale, representing 1.82 m in length/height.

**Supplementary Files**

This is a list of supplementary files associated with this preprint. Click to download.

* Supp.Material1.docx