Evidence of active magmatic rifting at the Ma’Alalta volcanic field (Afar, Ethiopia)

During continental rifting, strain and magmatism are believed to localize to narrow magmatic segments, while the rift margin is progressively abandoned. We integrate volcanological, geochemical, petrological and seismic data from the Ma’Alalta volcanic field (MVF) near the western margin of Afar, to show that the MVF is an active magmatic segment. Magmatism in MVF initiated with lava flows and large-volume, caldera-forming ignimbrites from a central edifice. However, the most recent magmatic activity shifted towards mafic lava fields, cinder cones and obsidian-rich silicic domes erupted from vents aligned NNW-SSE, defining a ~ 35-km-long magmatic segment. Along the same area, a NNW-SSE alignment of earthquakes was recorded by two local seismic networks (2005–2009 and 2011–2013). The geochemistry of the mafic rocks is similar to those of nearby axial volcanoes. Inferred magma storage depth from mineral geobarometry shows that a shallow, silicic chamber existed at ~ 5-km depth below the stratovolcano, while a stacked plumbing system with at least three magma storage levels between 9 and 24 km depth fed the recent basalts. We interpret the wide set of observations from the MVF as evidence that the area is an active magmatic segment, showing that localised axial extension can be heavily offset towards the rift margin.


Introduction
The Ma'Alalta volcanic field (MVF hereafter) is located in Afar (Ethiopia), ~ 50 km north of Dabbahu volcano (Fig. 1). MVF consists of a felsic stratovolcano with summit nested calderas surrounded by more recent mafic lava flows and scoria cones and felsic lava flows and domes. Due to its location at the edge of the Afar depression, the MVF was originally classified as a marginal unit (Barberi et al. 1974a;Wiart and Oppenheimer 2005). However, little is known about the petrology, geochemistry and tectonics of MVF and how they relate to the tectonics of continental breakup in Afar.
The Afar rift is a triple junction where the Main Ethiopian Rift, Red Sea and the Gulf of Aden rift arms meet. Volcanism in Afar initiated ~ 30 Ma ago and was initially dominated by flood basalts (Hofmann et al. 1997). Since at least ~ 1 Ma, active rifting and volcanism have focused to en-echelon axial magmatic segments (~ 70-km-long, ~ 20-km-wide) marking the current plate boundary of the rift (e.g. Dabbahu segment, Manda Hararo segment and Erta Ale segment) (Fig. 1a) (Oppenheimer and Francis 1998;Wright et al. 2012;Medynski et al. 2013). Holocene rift marginal volcanism has also produced large volumes of lavas, such as those erupted by the silicic stratovolcanoes of the Nabro-Mallahle alignment (marginal volcanism) or from basaltic shield volcanoes, such as Dabbayra (transverse volcanism) (Fig. 1a). Nevertheless, despite the peculiarity of the marginal volcanism, current knowledge about marginal volcanoes, their magma sources and plumbing system and the relationships with Editorial responsibility: N. Métrich nearby magmatic segments is partial because of their remote location and because scientific interest has focused mainly on axial volcanism.
We provide new major and trace element analyses and the first pressure and temperature estimates of the MVF rock samples. The rocks were collected during the French-Italian campaigns in 1967-1973(Varet 1978Barberi et al. 1977; http:// repos itori es. dst. unipi. it/ index. php/ afar-repos itori es). We integrated the volcanological, geochemical and petrological analyses with remote sensing-based mapping and local earthquake catalogues. We combined the diverse datasets in order to understand the characteristics of the MVF magmatic production, the relationship of the volcano with the other rift segments and the current magmatic and tectonic activity of Afar.

The Afar depression
The Afar depression is a triple junction between the Red Sea, Gulf of Aden and Main Ethiopian Rifts, created by the divergent motions of the Nubian, Arabian and Somalian plates. Afar is flanked by the Ethiopian Plateau to the west, by the Somalian Plateau to the south-east and the Danakil Block and Ali-Sabieh Block to the east. Initiation of volcanic activity in Ethiopia is thought to be caused by a mantle plume that produced the Ethiopian flood basalts (Marty and Gezahegn 1996;Pik et al. 1999;Kieffer et al. 2004) with the majority of the volume of basaltic magma erupted over only ~ 1 Myr at about 30 Ma (Zumbo et al. 1995;Hofmann et al. 1997). Starting from the Oligocene, magmatic activity occurred inside the newly forming Afar rift floor and progressively migrated towards the current magmatic segments through the Miocene (Wolfenden et al. 2005;Stab et al. 2016).
The Afar rift floor is mainly covered by the Stratoid Series (4-1 Ma), followed by the Gulf Basalt formation (1.1-0.6 Ma) (Barberi and Santacroce 1980;Kidane et al. 2003;Stab et al. 2016), while magmatism and extension are currently mainly localised to axial magmatic segments within the rift (Hayward and Ebinger, 1996;Wright et al. 2006; Fig. 1). In central and northern Afar, these magmatic segments (e.g. Dabbahu, Manda Hararo and Erta Ale) are oriented NNW-SSE, which is the orientation of the Red Sea (Fig. 1). The magmatic segments consist of fissural eruptions Fig. 1 a Regional tectonic map showing the three diverging plates (Nubia, Arabia and Somalia) and the main structures. In red the magmatic segments and in green marginal volcanoes (triangles) and transverse structures (segments). Bulletin of Volcanology (2021) 83: 38 associated with extensional fractures sometimes related to shield volcanoes (Oppenheimer and Francis 1998;Wright et al. 2006;Ferguson et al. 2010;Field et al. 2012Field et al. , 2013Medynski et al. 2015Medynski et al. , 2016Hagos et al. 2016). Volcanism close to the margins of the Afar floor is also reported. Barberi et al. (1970Barberi et al. ( , 1974b described transverse alignments of volcanoes (e.g. Dabbayra, Assab and Dubbi) situated both east and west of the magmatic segments in Afar (green segment in Fig. 1a) and oriented nearly orthogonally to the Red Sea trend. Furthermore, they also classified stratovolcanoes on either side of the rift, such as Nabro and Ma'Alalta, as marginal volcanoes (green triangles in Fig. 1a). These are large stratovolcanoes that produced caldera-forming explosive eruptions (Barberi et al. 1974a;Wiart and Oppenheimer 2005;Oppenheimer et al. 2019) and lava flows generally more alkaline with respect to the tholeiitic-transitional axial volcanism (Civetta et al. 1975;De Fino et al. 1978;Donovan et al. 2018). The most recent eruption of a marginal volcano was in 2011 at Nabro, which produced basaltic to basaltic trachyandesite tephra and lava (Hamlyn et al. 2018;Donovan et al. 2018). The MVF is located in the Northern sector of Afar, where the extension rate is ~ 20 mm/year in a ~ NE-SW direction (

The Ma'Alalta volcanic field
The MVF is surrounded by outcrops of Jurassic limestones and meta-sediment and meta-volcanic of the Neoproterozoic basement and is dominated by a 1745-m-high stratovolcano with nested calderas elongated NE-SW, approximately perpendicular to the Red Sea trend. The outer caldera is approximately 8 × 5 km, and the inner one is 5 × 2.5 km. At the base of the MVF sequence, Barberi et al. (1970) identified fissural basaltic lava flows, underlying three ignimbrite sheets alternated with pumice fall layers, and before the emplacement of the third ignimbrite unit, and basaltic scoriaceous deposits with associated small scoria cones. On top of this succession, they recognized obsidian flows, pyroclastic layers and porphyritic trachyte flows and domes that built the stratovolcano. The most recent products are felsic lava domes and obsidian flows, erupted along NNW-SSE fissures south of the stratovolcano. Barberi et al. (1974a) distinguished an earlier stratovolcano-building phase with felsic non-peralkaline lavas and a later phase that produced large-volume ignimbrites, intra-caldera domes and lava flows of peralkaline composition, particularly rich in sanidine phenocrysts (Barberi et al. 1970). The felsic domes and flows south of the stratovolcano are more peralkaline and silica-rich with respect to the stratovolcano products, and they are nearly aphyric (Barberi et al. 1970(Barberi et al. , 1974a. Wiart and Oppenheimer (2005) proposed a satellite-based map suggesting that the stratovolcano is flanked to the west and south-west by felsic domes and lava flows but also by recent mafic cinder cones and NNW-SSE fissural lava flows and to the east by more weathered and eroded lavas, named "Old Lava". K-Ar age determinations gave 0.55 ± 0.05 Ma for a lava sample at the base of the stratovolcano and 0.12 ± 0.05 Ma for a pantellerite obsidian south of the stratovolcano (Barberi et al. 1972).

Rock samples and mapping
For this study, we used forty rock samples from the Afar Repository of the University of Pisa (http:// repos itori es. dst. unipi. it/ index. php/ afar-repos itori es), collected between 1967 and 1973 in the MVF region (Fig. 2). The sampling sites have been identified using the original sampling map. The samples are mainly lavas with few pyroclastics (online resource Table 1) and cover the main phases of the MVF volcanic activity (Figs. 2 and 3). We examined the rock samples and thin sections and also investigated the spatial distribution of the erupted products using satellite images. We obtained a detailed map of the different lava flows of MVF by determining the type of erupted products through petrographic observations and reconstructing the spatial distribution of the products with Landsat satellite panchromatic images. We mapped all the volcanic and tectonic features, and by draping the satellite images onto the SRTM DEM (resolution of 1 arc-second), we interpreted the relative stratigraphy of the volcanic deposits. Based on the petrographic observations and the spatial distribution of the deposits, we selected lava samples for bulk rock and mineral chemical analysis.

Analytical methods for geochemistry and mineral chemistry
We selected 14 lava samples representative of the two phases of MVF and without signs of alteration in hand sample and thin section ( Fig. 3; Table 1). All the samples were crushed and milled in agate jars and analysed for major elements by means of XRF using the spectrometer ARL 9400 XP of the University of Pisa. The Loss on Ignition (LOI) was determined at 950 °C. Three of these samples were previously analysed for major elements (Barberi et al. 1974a, b;Barberi et al. 1970;Tazieff et al. 1969) and are listed in online resource Table 2 but have not been inserted in the figures (the other literature data are present in both online resource Table 2 and figures). Trace element analysis was performed on 13 samples (Table 1) using inductively coupled plasma mass spectrometry (ICP-MS) with a PerkinElmer NexION 300X of the University of Pisa. About 50-60 mg of each powdered sample was dissolved in a mixture of HF and HNO 3 on a hot plate at ~ 170 °C in screw-top perfluoroalkoxy (PFA) vessel, following the procedure used at the ICP-MS laboratory for igneous rocks. Ultrapure Milli-Q water was used for diluting the solution. A procedural blank and 3 references samples underwent the same procedure. The geochemical references materials used are RGM-1 (Rhyolite Glass Mountain) (Govindaraju 1994) and WS-E (Whin Sill Dolerite) (Govindaraju 1994) to verify the accuracy of the analysis and BE-N (alkali basalt) (Govindaraju 1980) to calibrate the equipment (Table 1). A solution containing 103 Rh, 187 Re and 209 Bi has been added to the unknown samples, blank and reference samples to be used as the internal standard.
Polished thin sections of five samples were studied by scanning electron microscopy (BSE imaging) and EDS microanalysis with the FEI Quanta 450 ESEM FEG provided with a QUANTAX XFlash Detector 6/10 at the University of Pisa. The mineral composition of five samples, one trachyte of the first phase and one rhyolite and three basalts of the second phase (online resource Table 3), was analysed using the Electron Probe MicroAnalyzer (EPMA) 8200 Super Probe at the Department of Earth Sciences "Ardito Desio" of the University of Milan. Analytical conditions for minerals were 15 kV accelerating voltage, 5 nA beam current, 3 μm beam size and counting time of 20 s on peaks and 10 s in background.

Seismicity
The seismological analysis was carried out using a combined dataset of local seismicity acquired during 2005-2013 in Afar (Fig. 2). The data were collected using consecutive local seismic networks during -2006(Ebinger et al. 2008), 2007-2009(Illsley-Kemp et al. 2018). The latter seismic network was more sensitive to earthquakes in MVF, as the closest station was ~ 23 km from the stratovolcano, while the minimum distance was ~ 50 km during the earlier networks. We used the published earthquake catalogues from the above studies and isolated all the earthquakes that occurred near the MVF between longitude 39. 90E-40.35E and latitude 12.70 N-13.20 N, giving a total of 229 events (Fig. 2

Geological map of MVF and relative chronology of the volcanic units
We identify two different phases of activity based on the relative chronology of the volcanic units, the petrographic analysis and the distribution of the erupted products in the satellite images (Fig. 3). The first phase started at 0.55 ± 0.05 Ma (Barberi et al. 1972) with the formation   Sample  of the stratovolcano and was later affected by a calderaforming stage, producing trachytic-rhyolitic lava flows and ignimbrites (Fig. 3, online resource Fig. 1). This activity was followed by the second phase (0.12 ± 0.05 Ma, Barberi et al. 1972) producing multi-vent mafic and felsic eruptions. This phase produced a NNW-SSE-trend of fissural basaltic-trachybasaltic lava flows and cinder cones and rhyolitic-trachytic lavas and domes (Fig. 3, online resource Fig. 1).

The stratovolcano and caldera-forming phase (first phase)
Based on the available samples of the first phase, two stages can be distinguished: a mainly effusive stage, consisting of metaluminous to peralkaline trachytic-rhyolitic rocks building the central edifice, called the stratovolcano felsic flows (SFF-1), and a caldera-forming stage, consisting of metaluminous to peralkaline trachytic pyroclastic flow deposits, called calderas felsic ignimbrites (CFI-2) which include the final silicic activity of the calderas ( Fig. 3; Table 1; online resource Table 1). We observed an upward variation in size, type and abundance of phenocrysts from the base to the top of the stratovolcano felsic lava flows. The seven samples from the top have up to 40 vol% of phenocrysts consisting of alkali feldspar reaching 1 cm in size and minor plagioclase, fayalite, slightly pleochroic clinopyroxene and, in one sample, arfvedsonite ( Fig. 4; online resource Fig. 2a; online resource Table 1). Plagioclase is mainly andesine and minor oligoclase (Ab 51-68 An 26-47 ); alkali feldspar is mainly anorthoclase (Ab 63-72 Or 15-36 ); clinopyroxene is augite (Wo 38-44 Fs 23-57 ) (online resource Fig. 3a; online resource Table 3). Apatite and Fe-Ti oxides are the main accessory phases. The groundmass is holocrystalline or nearly so, and it is composed of the same crystals as the phenocrysts except for fayalite. Two samples from the top rim of the inner caldera display phenocrysts of quartz in addition to alkali feldspar and minor pleochroic clinopyroxene, arfvedsonite and aenigmatite (online resource Fig. 2b; online resource Table 1). The three samples from the base of the stratovolcano have lower phenocryst abundance (up to 20 vol%) and smaller phenocryst in size compared to samples from the top (Fig. 4; online resource Fig. 2c; online resource Table 1). Also, no minerals indicative of peralkalinity are observed.

The multi-vent mafic and felsic phase (second phase)
During the second phase, fissural eruptions occurred along a ~ 35-km-long NNW-SSE trending zone. We identified three groups of products-old mafic lava flows (OMF-3), felsic lava flows and domes (FFD-4) and recent mafic lava flows (RMF-4)-erupted by fissural vents with numerous cinder cones forming monogenetic and polygenetic small volume volcanic edifices distributed along the NNW-SSE trending zone (Fig. 3).
The OMF-3 and RMF-4 groups (Fig. 3) are basalt and trachybasalt lava flows and cinder cones (online resource Fig. 1; Table 1; online resource Table 1) that differ in their morphological features and degree of weathering ( Fig. 5a,b). The majority of OMF-3 lavas and poorly preserved cinder cones (Fig. 5a) are east of the volcano, while the RMF-4, associated to relatively well-preserved cinder cones (Fig. 5b), are located along the northern sector of the NNW-SSE trending zone and at its southern tip. All the mafic products show olivine, clinopyroxene and plagioclase phenocrysts, with large variations in their relative abundance ( Fig. 4; online resource Fig. 5; online resource Table 1). The groundmass is from hypohyaline to hypocrystalline, with clinopyroxene, plagioclase and opaques. In the mafic samples, plagioclase is labradorite to bytownite (Ab 17-44 An 54-82 ), olivine is Fo 77-84 (the tephroite component is Tp 0.4-0.2 ), and clinopyroxene is augite with minor diopside (Wo 39-45 Fs 9-16 ) (online resource Table 3; online resource Fig. 3c). Presence of reversed zoning in clinopyroxene, plagioclase and olivine crystals has been observed. In one case (sample G8), the large variation of clinopyroxene core composition in the augite field towards the hedenbergite endmember (Wo 41-48 Fs 9-32 ) compared to the more homogeneous rims (Wo 43-47 Fs 9-13 ) indicates mixing processes with evolved and peralkaline lavas. Sample D170 shows common reverse zoning in olivine, again an evidence of mixing of magmas with different degrees of chemical evolution.

Relative chronology of the volcanic units
The oldest products of MVF are the SFF-1, dated at 0.55 ± 0.05 Ma (Barberi et al. 1972), followed by the CFI-2 (Fig. 3b). We can identify at least two major explosive events that produced ignimbrite deposits, based on the presence of two calderas.
Our interpretation of the relationship between the first phase and the OMF-3 differs from Wiart and Oppenheimer (2005). We suggest that the neat and relatively straight contact (N-S) on the eastern side of the stratovolcano (Fig. 5c) and the outcropping of the OMF-3 products on the western flank of the stratovolcano (online resource Fig. 6) indicate that they post-date the stratovolcano. Otherwise, the pyroclastic flows responsible for the ignimbrites would be easily emplaced on the OMF-3 east of the stratovolcano due to their proximity to the edifice and the absence of a topographic boundary. This would lead to an irregular contact due to the overlapping of the ignimbrites and a fanlike geometry of the felsic flows (instead of an abrupt stop at the contact with the OMF-3; Fig. 5c; online resource Fig. 6). Regarding the western flank of the stratovolcano, it is unlikely that the OMF-3, if emplaced before it, would result in topographic continuity (similar slope and altitude; online resource Fig. 6) with the edifice and not covered by the stratovolcano products. Moreover, an outcropping portion of the ignimbrite deposit on the eastern flank of the stratovolcano is completely surrounded by OMF-3 mafic lava flows. These last observations indicate that the products of the first phase were deposited before the OMF-3. Based on the degree of erosion of lava flows and cinder cones, we classify OMF-3 as older than RMF-4 ( Fig. 5a,b; online resource Fig. 6). The FFD-4 south of the stratovolcano is dated by Barberi et al. (1972) at 0.12 ± 0.05 Ma. At the small scale, the well-preserved fan-like geometry of the northernmost lava flow of FFD-4 (preserved also in the portion where FFD-4 and RMF-4 are not in contact) indicates that FFD-4 precedes RMF-4 but, at the same time, the undulated  Table 1 front at the contact with the basalts (online resource Fig. 6) would suggest the opposite. Due to these uncertain observations of the only contact between RMF-4 and FFD-4 and the heterogeneous colour and surface morphology of both RMF-4 and FFD-4 (darker and lighter portions in Fig. 3), we suggest that they were likely emplaced during roughly the Fig. 5 Pairs of Landsat 5 satellite images and the relative geological map. SFF-1, stratovolcano felsic flows; CFI-2, calderas felsic ignimbrites; OMF-3, old mafic flows; RFM-4, recent mafic flows; FFD-4, felsic flows and domes; BAS basement. Locations are shown in online resource Fig. 6. World Imagery Source as in Fig. 2 Bulletin of Volcanology (2021) 83: 38 same time period. Despite the lack of certainty about their age relationship, as a whole, they are the youngest products of the MVF eruptive history as they are emplaced after the stratovolcano activity, which also influences the flow directions of the upper products (Fig. 5d, online resource Fig. 6).

Major and trace element geochemistry
In the total-alkali vs silica diagram (TAS; online resource Fig. 1), the analysed MVF rocks fall in the basalt, trachybasalt, trachyte and rhyolite fields and define a compositional gap (Daly Gap), from 51 to 64 wt% SiO 2 .
The SFF-1 rocks are mainly trachytes with few rhyolites (online resource Fig. 1). The samples from the top of the edifice are prevalently peralkaline (comenditic trachyte and pantellerite), while the rocks at the base are metaluminous trachytes, although very close to peralkalinity, and a pantellerite (Table 1; online resource Table 2). The variable degree of chemical evolution of the samples defines a differentiation trend of decreasing Al 2 O 3 , CaO, TiO 2 and P 2 O 5 , corresponding to fractionation of feldspar, pyroxene, Ti-Fe oxides and apatite, with nearly constant alkalies, while FeO t increases in the most evolved rhyolitic samples (Fig. 6a). Accordingly, the elements Ba, Sr and Eu, strongly compatible in feldspars, reach the maximum depletion in the most evolved sample (Figs. 6b and 7).
The FFD-4 products are prevalently rhyolitic (online resource Fig. 1) and are more silica-and alkali-rich and more peralkaline (comendite and comenditic trachyte) with respect to the SFF-1 samples taken as a whole (Table 1; Fig. 6; online resource Table 2). The trace element patterns of FFD-4 samples are similar to each other and have greater negative anomalies in Sr, Eu, Ba, Ti and P compared to the stratovolcano (Fig. 7), which may indicate a higher degree of differentiation by fractionation of feldspar, Ti-Fe oxides and apatite. The G3 comenditic trachyte (the least evolved sample of the FFD-4) is the least enriched in REE and has the smallest negative Eu and Ba anomalies (Fig. 7) and shows the same ratios of highly incompatible elements (e.g. Zr/Nb; online resource Fig. 8).
The OMF-3 and RMF-4 products are basalts and trachybasalts (online resource Fig. 1; Table 1; online resource  Table 2). All the samples have relatively low Mg# (< 60) and are all olivine-hypersthene normative (online resource Table 4). They all have no Eu anomalies (Fig. 7). Fractionation of olivine, clinopyroxene and Cr-spinel is active in the mafic reservoirs, as shown by the range in MgO and CaO (Fig. 6a) and Ni and Cr (online resource Fig. 7).
Overall, the whole distribution of the compatible elements indicates differentiation by crystal fractionation of olivine, pyroxene, feldspar, Cr-spinel, Ti-Fe oxides and apatite starting from transitional basalts, as already proposed for Dabbahu volcano (Barberi et al. 1975;Field et al. 2013; Fig. 6a; online resource Fig. 7). The distribution of the incompatible elements is also similar to the Dabbahu volcano samples (online resource Fig. 7, 8; Hutchison et al. 2018) with the exception of Y and Zr, slightly depleted in the MVF products for the same Th content. Moreover, the distribution of the MVF incompatible elements (Nb, Zr, Ta, Y, Vs and Th in online resource Figs. 7 and 8) shows good positive correlations going through zero (e.g. Th/ Rb = 0.12 ± 0.02, Nb/Zr = 0.20 ± 0.02) suggesting a genetic link between mafic and the whole felsic products. Due to the small dataset and to the presence of the Daly Gap, it is however difficult to give a certain interpretation regarding the correlation between mafic and felsic products.
The comparison between the felsic rocks of SFF-1 and FFD-4 indicates different extent of fractional crystallization of feldspar, as suggested by the negative anomalies of Eu, Sr and Ba and of apatite and Ti-oxide as indicated by the negative anomaly of P and Ti, respectively (Fig. 7). Moreover, the Sr and Ba are remarkably lower in FFD-4 for the same Th (or Nb and Zr) content with respect to SFF-1 (Fig. 6b). Since accumulation of feldspar in SFF-1 can be excluded because it is not confirmed by Al 2 O 3 vs. SiO 2 (Table 1; Fig. 6a), a major fractionation of feldspar for FFD-4 at the same evolution degree of SFF-1 is hypothesized. The trace element differences indicate that the two felsic rock suites were originated by different differentiation paths.

Geothermobarometry
Our geothermobarometric calculations (Table 2) based on clinopyroxene of the SFF-1 trachyte indicate crystallization pressures on the order of 1.34-1.49 ± 1.15 kbar. The crystallization temperature was 914-925 ± 24 °C for clinopyroxene, 902-954 ± 36 °C for plagioclase and 857 ± 23 °C for alkali feldspar. A higher crystallization pressure is obtained for the RMF-4 mafic magma of the lava flow south of the stratovolcano (3.8-4.1 ± 1.4 kbar for clinopyroxene and 1.1-4.0 ± 3.8 kbar for plagioclase). The crystallization temperature was 1146-1176 ± 45 °C for clinopyroxene, 1091-1162 ± 36 °C for plagioclase and 1140-1178 ± 43 °C for olivine (only two analyses). The crystallization pressure obtained for the RMF-4 mafic magma of the lava flow north-west of the stratovolcano is still higher, 6.3-6.6 ± 1.4 kbar for clinopyroxene and 4.7-7.8 ± 3.8 kbar for plagioclase. The crystallization temperature was 1161-1193 ± 45 °C for clinopyroxene, 1134-1210 ± 36 °C for plagioclase and 1147-1184 ± 43 °C for olivine. Finally, the crystallization pressure obtained for the RMF-4 magma associated to a scoria cone in close proximity to the FFD-4 felsic centres, just south of the stratovolcano, is 2.5-2.8 ± 1.4 kbar. This estimation came from analysis of one clinopyroxene core and four rims which were found in equilibrium with melt. This is supported by ten analyses of rims close to the equilibrium range (KD (Fe-Mg) cpx−liq = 0.31 and 0.32) giving crystallization pressures of 2.4-2.7 ± 1.4 kbar. The crystallization temperature for core and rim analyses in equilibrium is 1125-1155 ± 45 °C. We then derived the corresponding depths assuming a 270 bar/km geobaric gradient (Best 2013) and using clinopyroxene results because of the lower errors. The obtained crystallization pressures correspond to a5-km depth for the SFF-1 trachyte and to 9-10-km, 14-15-km and 23-24 km depths for the mafic products.  Barberi et al. 1975. b Ba and Sr trace element diagrams. The Dabbahu volcano whole rocks analyses are from Hutchison et al. (2018). SFF-1, stratovolcano felsic flows; FFD-4, felsic flows and domes; RMF-4-OMF-3, recent-old mafic flows

Structures and seismicity
En-echelon faults which are ~ 2-km long strike roughly parallel to the Red Sea trend and are clearly visible in the younger units at the southern tip of the MVF, and another NE-SW trend is defined by the elongation of the MVF calderas (Fig. 8). Furthermore, the distribution, morphology and elongation direction of the cinder cones are consistent with the two main trends (Fig. 8; Table 3). We identified three different groups of cinder cones based on their spatial distribution, and when the cone shapes were preserved, we also measured the direction of cone alignment using at least 3 cones closer than 400 m from each other, and the elongation direction using cones with aspect ratio > 1.4. The largest group of cinder cones is along the NNW-SSE trending zone west of the stratovolcano with up to 70 cones distributed along a zone striking 332° N, roughly parallel to the Red Sea trend. Also, the alignment of the cones within this group has a trend varying between 327° E and 351° E, and the cones show predominant elongation varying from 321° E to 17° E, both subparallel to the Red Sea trend (Fig. 8). Altogether, the products of the second phase west of the stratovolcano occur over a ~ 35-km-long, NNW-SSE trend. East of the stratovolcano, two smaller groups of cones have been identified associated only with the OMF-3. A group of 19 cones show a trend of 317° N, similar to the western cones, while the remaining 26 cones are aligned 49° N, parallel to the elongation direction of the calderas. Most of the cones east of the calderas are not elongated nor aligned (Fig. 8).
The majority of the seismicity in MVF occurs along a ~ 35-km-long, NNW-SSE trending region with the alignment of earthquakes corresponding to the largest group of cinder cones and to the recent flows (FFD-4 and RMF-4) (Fig. 8). The magnitude varies between 1.06 and 3.37 and Fig. 7 Chondrite-normalized REE patterns and primitive mantle-normalized trace element spider diagrams. SFF-1, stratovolcano felsic flows; FFD-4, felsic flows and domes; RMF-4-OMF-3, recent-old mafic flows. Normalizing values after McDonough and Sun (1995) gives a mean of 2.22. Most earthquakes in the alignment are observed at the southern end with good correlation to where the most recent RMF-4 was erupted and where faults cutting these flows and surficial sediments occur. A smaller group of earthquakes is located in eastern parts of the MVF in the area between the OMF-3 and the basement, indicating that this region is not the main locus of seismicity at present. No significant seismicity occurs at the MVF calderas.

Evolution of the MVF: from caldera activity to rift segment development
Two main phases of activity were identified during the evolution of MVF, based on the integration of volcanological, geochemical, petrological, structural and seismological evidence.
The first phase began at 0.55 Ma (Barberi et al. 1972) with explosive and effusive activity from a central vent  that built the stratovolcano. The two caldera structures suggest that the volcano was fed by a shallow crustal magma chamber, at ~ 5-km depth according to geobarometry applied to the lavas of the stratovolcano summit. This phase is dominated by structures trending NE-SW as indicated by the caldera elongation. Furthermore, this trend has been observed also in the NE-SW alignment of OMF-3 cinder cones east of the stratovolcano. This cinder cone alignment suggests that a tectonic structure orthogonal to the current rift axis (transverse alignment) existed during the formation of the stratovolcano and controlled the emplacement of the OMF-3 cinder cones aligned NE-SW. The large ranges in the concentrations of incompatible trace elements of SFF-1 indicate efficient differentiation processes in the crustal chambers leading from trachytes to peralkaline rhyolites. An estimate of the extent of differentiation made using Zr as fractionation index due to its high incompatibility in peralkaline rocks indicates formation of the peralkaline rhyolite by 60% fractional crystallization starting from the trachytic melts. No mafic activity has been observed for the first phase meaning that the basic magmas could have remained at depth, unerupted, as suggested for some MER volcanoes (Peccerillo et al. 2003;Hutchison et al. 2016), or could be mostly covered by later activity and thus difficult to identify and distinguish from the later mafic products. It is also possible that the compositions of the mafic products of the second phase represent the parental magma of the stratovolcano felsic products, erupted after the first phase and away from the stratovolcano. Since the stratovolcano is now heavily eroded and currently lacks any significant seismic activity, we suggest that it is now dormant or extinct. During the second phase, both mafic and felsic magmas were erupted. The mafic magma storage depths range from ~ 9 to ~ 24 km, and the presence of reversed zoning in clinopyroxene, plagioclase and olivine crystals (online resource Table 3; online resource Fig. 3c) suggest a stacked plumbing system where the mafic magmas undergo differentiation by fractional crystallization and magma mixing. In addition, connections between the mafic and felsic feeding systems of the second phase are revealed by evidence of interaction between basalts and felsic magmas, indicated by the presence of two populations of clinopyroxene in the mafic lavas erupted nearby the sialic centres (as a whole Wo 41-48 Fs 9-32 ). This suggests that the mafic magmas encountered mush-like felsic pockets during ascent and captured high Fs clinopyroxene without other clear evidence of magma mixing. The felsic products of the second phase (0.12 Ma from Barberi et al. 1972) are mostly aphyric and effusive, systematically more evolved (except for G33) and more peralkaline compared to the stratovolcano products and occupy a very narrow chemical range. It is possible that these aphyric and acidic lavas originated from highly fractionated melts rising from a relatively deep reservoir slowly enough to allow efficient gas/melt decoupling, preventing explosive activity (Gonnermann and Manga 2007); efficient crystal-melt separation in the magma chamber or adiabatic decompression-induced crystal resorption may explain the aphyric texture. This phase is characterized by the development of structures parallel to the NNW-SSE Red Sea trend and a shift in activity to the western side of the stratovolcano, where the youngest structures including RMF-4, FFD-4 and cinder cones occur. The youngest volcanism and structures correlate spatially with the majority of the recent seismicity, suggesting that this region is the current locus of strain and magmatism. We also interpret that the seismicity at the southern end of the MVF is caused by faulting on NNW-SSE planes along the rift axis, as short en-echelon faults are mapped across a ~ 10-km-wide axial zone (Fig. 8). En-echelon faults can also feed fissure eruptions during episodes of magma migration.
The marked changes in the types of erupted products, the different magma chamber depths, the distribution of the recent seismic activity and the structural arrangement between the two evolutionary phases of the MVF show that the recent activity of MVF differs from the first phase and it is now typical of an active magmatic rift segment.
A petrogenetic discussion of the magmatic evolution of MVF was not the focus of this work due to lack of mafic and intermediate products and of a complete sampling coverage of the area. We are aware that new studies on the stratovolcano, and in particular on the associated mafic products if outcropping, are needed to test our geochemical and petrological interpretation. Nevertheless, a process of fractional crystallization from transitional basalts to rhyolite, similar to those described for Dabbahu volcano (Barberi et al. 1975;Field et al. 2013), can be envisaged for the MVF products, with FFD-4 as a whole reaching a higher extent of differentiation with respect to SFF-1. We also show differences in Sr and Ba for the same degree of evolution (Zr or Th content) indicating differences in the fractionation process between the two felsic phases. We therefore speculate that the described chemical differences, in particular the different trace element ratios of felsic rocks between the two phases (Sr and Ba; Fig. 6b) are not only due to different degrees of fractional crystallization but also due to major differences in the evolutionary processes. The major differences in the evolutionary processes are possibly linked to different P-T-H 2 O conditions of magma storage or variable extent of crustal assimilation that combined with the differences regarding the erupted products and the geobarometric estimates could be related to differences in the plumbing system architecture.

The MVF magma source and comparison to other rift volcanoes
Taken as a whole, the MVF products show a similar distribution of major and trace elements to Dabbahu volcano, except for the G33 lavas of the stratovolcano summit, reaching very high incompatible element contents ( Fig. 6; online resource Figs. 7 and 8). In particular, the RMF-4 phase of MVF and the on-axis volcanism of Dabbahu volcano (Field et al. 2013) are both characterized by eruption of a significant volume of transitional mafic magmas. The recent mafic flows, just like at Dabbahu volcano, are associated with aphyric and obsidian-rich comendite erupted from several scattered vents and producing small volume domes or lava flows.
Lithosphere thinning is an important process influencing the depth of partial melting in the rising asthenosphere beneath rift systems. The across-and along-rift variation in the extent of stretching and thinning of the lithosphere leads to partial melting by decompression at different depths and consequently to different geochemical composition of the basaltic melts (Rooney 2010). We compared the chemical composition of the MVF mafic products (MgO wt% > 4) to similar mafic rocks of other Afar volcanoes by means of rare earth element ratios (Fig. 9). We use the La/Sm ratio to evaluate the degree of partial melting and/or source enrichment and the ratio Tb/Yb for the depth of the melting column (Rooney 2010;Wang et al. 2002). In fact, magmas generated at a depth sufficient to have residual garnet in the source show a marked fractionation in the middle and heavy rare earth elements (e.g. Tb/Yb, Dy/Yb; Wang et al. 2002). High pressure clinopyroxene of a spinel lherzolite source behaves somewhat like garnet, but with much lower effect (Rooney 2010).
We considered the marginal Nabro volcano (Donovan et al. 2018) and the transverse Assab range (Teklay et al. 2010), the on-axis Dabbahu volcano (Barberi et al. 1975;Hoare et al. 2020), the Badi volcano in the Dabbahu segment , the products from the Dabbahu and Manda Hararo segments Barrat et al. 2003) and the Erta Ale (Barrat et al. 1998) magmatic segment fissural products. As a whole, the mafic products of these volcanic centres range from 0.8 to 3.8 in La N /Sm N and from 1.0 to 2.2 in Tb N /Yb N (Fig. 9), crossing the Tb N /Yb N values at the transition between the spineland the garnet-bearing mantle sources derived from magma modelling in the Basin and Range province (Wang et al. 2002). Contributions from melting a garnet-bearing lower crust can be ruled out due to the low crustal thickness, Fig. 9 Plot of Tb/Yb and La/Sm normalized to chondritic values (McDonough and Sun 1995). Only samples with MgO greater than 4 wt% are plotted. The division between spinel and garnet-dominated melting is derived from magma modelling in the Basin and Range province (Wang et al. 2002). As a comparison with the MVF, volcanism of different locations of the Afar region is shown: the marginal Nabro volcano (Donovan et al. 2018) and the transverse Assab range (Teklay et al. 2010), the on-axis Dabbahu volcano (Barberi et al. 1975;Hoare et al. 2020), the Badi volcano in the Dabbahu segment , the products from the Dabbahu and Manda Hararo segments Barrat et al. 2003) and Erta Ale (Barrat et al. 1998) rift axis fissural products. See text for explanation ranging from 15 km beneath the Danakil Depression to 25-30 km beneath the Danakil microplate . The fractionation of garnet under water-saturated conditions observed in some arc settings (Alonso-Perez et al. 2009) can also be excluded given the lack of hydrous phases in Afar mafic volcanism. Therefore, we interpret the variation in Tb N /Yb N as the result of changes in the melting column depth due to variations in lithospheric thickness. The marginal volcanism (i.e. Nabro and Assab range) is related to a deeper, garnet-bearing source (Tb N /Yb N > 1.8), while the MVF is similar to the axial volcanism (i.e. Badi volcano, Dabbahu volcano, Dabbahu and Manda Hararo segments and Erta Ale) and is related to a shallower source (Tb N /Yb N < 1.8, spinel lherzolite in the Wang et al. 2002 model for Basin and Range, while still in the garnet stability field according to Ferguson et al. (2013)). This difference is in agreement with the more evolved rifting stage (and thinner lithosphere) of the axial magmatic segments and strengthens the analogy between MVF recent volcanism and axial volcanism.
Therefore, we suggest that the geochemistry of the available mafic samples of the MVF, consistent with the other observations of this work, may be taken as an indication that the recent phase of the MVF is more typical of the axial volcanism. Additional work specifically on the MVF mafic rocks is necessary to fully define the mantle sources and melting conditions. However, based on several different lines of evidence from volcanology, geochemistry, seismicity and petrology, we propose that MVF evolved from an initial felsic central stratovolcano to an active magmatic segment.
Strain and magmatism are generally thought to progressively localize in-rift to narrow magmatic segments and the rift margin being progressively abandoned during magma-rich rifting (Hayward and Ebinger 1996;Wolfenden et al. 2005). However, here we have presented evidence of a young axial segment close to the rift margin of Afar ( Fig. 1 and Fig. 3), suggesting that the process of strain localization may be more complex. Studies from Iceland suggest that rift zones can initiate as several subparallel overlapping rifts during incipient rifting. Some zones do not fully evolve like the Öraefajökull-Snaefell volcanic flank zone which initially developed on a preexisting cross-rift structure and then was abandoned (Hards et al. 2000). On the other hand, subparallel overlapping rifts can be active at the same time in response to the plate migration with respect to a stable plume, like the Eastern and Western volcanic zones (LaFemina et al. 2005). Using this analogy, the origin of MVF could be linked to a pre-existing lithospheric weakness guiding where strain localizes (e.g. Muirhead and Kattenhorn 2018) or to the interaction between the plate boundary and a hotspot.

Conclusions
In this paper, we integrate volcanological, geochemical, petrological and seismic data for the Ma'Alalta volcanic field, near the western margin of Afar, and we present a new geological map coupled with geochemical and petrological analyses. Growth of the MVF started at around 0.55 Ma with the building of a stratovolcano, later modified by caldera-forming explosive eruptions and by effusion of silicic intracaldera domes. The volcanic structures indicate a NE-SW (across-rift) trending tectonic regime in this first phase. In contrast, the recent volcanic activity of MVF shifted towards mafic lava fields, cinder cones and peralkaline obsidian-rich silicic domes erupted from several vents aligned NNW-SSE, parallel to the Red Sea rift. Inferred magma storage depths from mineral geobarometry suggest a shallow chamber at ~ 5 km below the silicic stratovolcano, while a stacked plumbing system with polybaric magma storage levels between ~ 9 and ~ 24 km depth fed the recent basalts. The geochemical and petrographic features indicate distinct plumbing systems and differentiation paths of felsic magmas between the early (stratovolcano) and recent volcanic phases. The recent NNW-SSE-trending mafic rocks geochemical features are similar to the nearby Dabbahu and Badi volcanoes and differ from marginal activity. The pattern of volcanism and recent seismicity defines a ~ 35-km-long magmatic segment. We interpret the wide set of observations as evidence that the activity at the MVF evolved from a felsic central stratovolcano to an active rift segment, showing that localised axial extension can be heavily offset towards the rift margin.
Data availability The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.
Code availability Not applicable.

Conflict of interest
The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.