1. Introduction
The Euphrates volcanic field is located in northeastern Syria near its borders with Iraq and Turkey. The volcanic field occurs at a tectonically active region and represents the northern segment of the Middle East Cenozoic volcanic province (Fig. 1). This province consists of large, discontinuous volcanic exposures extending from Sinai, Jordan and Israel, to Lebanon, Syria and Turkey (Dubertret, Reference Dubertret1955; Mouty et al. Reference Mouty, Delaloye, Fontignie, Piskin and Wagner1992; Shaw et al. Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003; Abdel-Rahman & Nassar, Reference Abdel-Rahman and Nassar2004; Weinstein et al. Reference Weinstein, Navon, Altherr and Stein2006; Lustrino & Sharkov, Reference Lustrino and Sharkov2006). Cenozoic volcanism has produced some other extensive volcanic provinces, although of very different styles, further south in east Africa (mostly Kenya and Ethiopia), in Yemen and in Arabia (Mohr, Reference Mohr1983; White & McKenzie, Reference White and McKenzie1989; Camp & Roobol, Reference Camp and Roobol1992; Camp, Roobol & Hooper, Reference Camp, Roobol and Hooper1992; Baker, Thirlwall & Menzies, Reference Baker, Thirlwall and Menzies1996; Volker et al. Reference Volker, Altherr, Jochum and McCulloch1997; Bertrand et al. Reference Bertrand, Chazot, Blichert-Toft and Thoral2003; Furman et al. Reference Furman, Bryce, Karson and Iotti2004, Reference Furman, Kaleta, Bryce and Hanan2006).
The Quaternary mafic lavas of northern Kenya were interpreted by Furman et al. (Reference Furman, Bryce, Karson and Iotti2004) as being derived from mixtures of plume and ambient mantle sources. Furman et al. (Reference Furman, Kaleta, Bryce and Hanan2006) demonstrated that the Tertiary mafic lavas of Kenya reflect isotopically distinct plume heads beneath Tanzania and Afar that are derived from the chemically heterogeneous South African superplume. In Ethiopia, the Cenozoic lava flows were interpreted to have been developed in association with the Afar plume (Barberi et al. Reference Barberi, Ferrara, Santacroce, Treuil and Varet1975; Mohr, Reference Mohr1983; Hart et al. Reference Hart, Wolde, Walter and Mertzman1989; Stewart & Rogers, Reference Stewart and Rogers1996; Barrat et al. Reference Barrat, Fourcade, Jahn, Cheminèe and Capdevila1998; Pik et al. Reference Pik, Deniel, Coulon, Yirgu and Marty1999; George & Rogers, Reference George and Rogers2002), and are clearly related to the East African rift system. According to Ebinger & Sleep (Reference Ebinger and Sleep1998), Cenozoic magmatism throughout east Africa resulted from impact of a single plume.
Baker, Thirlwall & Menzies (Reference Baker, Thirlwall and Menzies1996) demonstrated that the Oligocene flood basalts of western Yemen are also Afar plume products like those of similar age in Ethiopia. The Plio-Quaternary mafic lavas of western Yemen were interpreted to have been produced in response to small amounts of lithospheric extension that were metasomatized and hydrated by the Afar plume (Baker et al. Reference Baker, Menzies, Thirlwall and Macpherson1997, Reference Baker, Chazot, Menzies and Thirlwall1998). Several comprehensive studies were carried out on the Cenozoic continental flood basalts of Saudi Arabia (e.g. Camp & Roobol, Reference Camp and Roobol1989, Reference Camp and Roobol1992; Camp, Roobol & Hooper, Reference Camp, Roobol and Hooper1992). These were interpreted to have occurred in two distinct phases. Phase one (30 to 20 Ma) of tholeiitic to transitional lavas was related to passive-mantle upwelling during extension of the Red Sea basin, and phase two (12 Ma to Recent) that produced transitional to strongly alkalic lavas was related to active-mantle upwelling facilitated by minor continental extension (Camp & Roobol, Reference Camp and Roobol1992). According to White & McKenzie (Reference White and McKenzie1989) and Camp & Roobol (Reference Camp and Roobol1992), the source of upwelling is either a mantle plume centrally located beneath the West Arabian Swell, or an elongated and extended lobe of hot asthenosphere emanating from the Ethiopian mantle plume.
In sharp contrast, the Afar plume of Ethiopia is thought not to have been channelled northwestwards beneath the Arabian plate, and it has played no role in producing the Cenozoic volcanic field of Jordan (Shaw et al. Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003). These Jordanian alkali basalts and basanites were formed in response to a phase of lithospheric extension. Miocene basaltic magmatism in Sinai is interpreted by Baldridge et al. (Reference Baldridge, Eyal, Bartov, Steinitz and Eyal1991) to have been related to the opening of the Red Sea. Garfunkel (Reference Garfunkel1989) described the ‘Cover’ Cenozoic basalts, which are located in southeastern Galilee and in southern Golan heights, as mostly alkaline basalts of intraplate character and ocean-island basalt affinities and proposed that the Jordan Valley depression could have acted as a route for magma ascent. According to Weinstein et al. (Reference Weinstein, Navon, Altherr and Stein2006), basanites and alkali basalts of the southern Golan and Galilee region have two mantle sources (a dehydrated garnet-bearing amphibole peridotite, and amphibole-garnet-bearing pyroxinite veins) occurring at lithospheric depths, and were formed by dehydration and partial melting of an originally isotopically uniform reservoir. Abdel-Rahman & Nassar (Reference Abdel-Rahman and Nassar2004) interpreted the Pliocene alkali basalts of northern Lebanon to have erupted during an episode of localized extension, which occurred particularly at the junction between the Dead Sea–Ghab transform (Fig. 2a, b) and the Yammouneh restraining bend (that is, in a transtensional tectonic regime).
In their general study that included Cenozoic volcanism in Syria, Mouty et al. (Reference Mouty, Delaloye, Fontignie, Piskin and Wagner1992) presented preliminary geochemical data and produced whole-rock K–Ar ages ranging from 23.5 to 1.5 Ma. They demonstrated that there was a large gap (between 16 and 8 Ma) in Cenozoic volcanic activity, and interpreted this to correspond with the interval between the two stages of spreading of the Red Sea–Dead Sea rift system. Mouty et al. (Reference Mouty, Delaloye, Fontignie, Piskin and Wagner1992) further proposed that the rocks are mostly within-plate transitional to alkaline basalts, and that magmatism was associated with left-lateral movements along the Dead Sea rift. In their study of the Neogene alkaline mafic lavas of western Syria, Lustrino & Sharkov (Reference Lustrino and Sharkov2006) indicated that much of the volcanic activity in that region is concentrated near the Dead Sea transform fault system. They concluded that the origin of the lavas of western Syria reflects a strong lithospheric control on the loci of partial melting which developed as a result of mantle decompression formed as a consequence of the transition from strongly compressive to transtensional stresses, and that there is no evidence to relate volcanic activity to northwestward channelling of the Afar plume or to any other mantle plume.
Near the Bitlis suture zone, a large crustal segment (the Karacadağ volcanic plateau of southern Turkey; Fig. 2b), occurring to the north of the investigated volcanic field, is covered largely by Late Cenozoic lavas. K/Ar ages reported in Pearce et al. (Reference Pearce, Bender, De Long, Kidd, Low, Güner, Saroğlu, Yilmaz, Moorbath and Mitchell1990), and Notsu et al. (Reference Notsu, Fujitani, Ui, Matsuda and Ercan1995) on these transitional to alkaline mafic lavas range from Upper Pliocene to Recent. Pearce et al. (Reference Pearce, Bender, De Long, Kidd, Low, Güner, Saroğlu, Yilmaz, Moorbath and Mitchell1990) proposed that the lavas of southern Turkey were derived from the mantle lithosphere of the Arabian continent which had been enriched by small volumes of asthenospheric melts.
No modern petrological–geochemical data on the northeastern Syrian Cenozoic basalts are available in the published international literature, and their magma source and petrogenesis are far from certain. The aim of our study is as follows: (1) to provide petrological, chemical and isotopic data on the Euphrates volcanic field of Syria, (2) to characterize its mantle source and propose a petrogenetic model for its magma generation, and (3) to evaluate the tectonic regime of magma emplacement in the context of the geodynamic framework and nearby plate boundaries.
2. Geological setting
The Euphrates volcanic field of northeastern Syria is located between latitudes 35°00′ and 36°30′N and longitudes 39°30′ and 41°00′E (Fig. 2c). It occurs on the northern flank of the Arabian plate, and is situated near major plate boundaries and tectonically active regions (Fig. 2a, b). The Bitlis suture, the left-lateral East Anatolian fault, and the Zagros fold belt occur to the north and east of the investigated volcanic field. The left-lateral Dead Sea–Ghab transform fault system, located further to the west of the volcanic field, has been active since the end of the Miocene (Lyberis et al. Reference Lyberis, Yurur, Chorowicz, Kasapoglu and Gundogdu1992), and permits the northward movement of the Arabian plate and its collision with Eurasia along the Bitlis suture, occurring just to the north of the Euphrates volcanic field (Fig. 2a, b).
The Euphrates graben/fault system was developed by extension/transtension following a period of extensive uplift and erosion during the Jurassic and Early Cretaceous (A. Koopman, unpub. Syrian Petroleum Company Report, 1986; Sawaf et al. Reference Sawaf, Al-Saad, Gebran, Barazangi, Best and Chaimov1993). Litak et al. (Reference Litak, Barazangi, Beauchamp, Seber, Brew, Sawaf and Al-Youssef1997) interpreted the reactivation of deep-seated fractures beneath the Euphrates graben to have resulted as a consequence of Neogene collision of the Arabian plate with Eurasia along the Bitlis suture. The presence of such an extensional graben near this collisional zone could also be interpreted to reflect the ‘escape’ of Turkey along the Anatolian fault system.
According to Giannerini et al. (Reference Giannérini, Campredon, Féraud and Zakhem1988), a change in plate motion occurred during Plio-Quaternary times, when the Arabian plate rotated to the northeast. This resulted in a transversal tension along the Levantine fault, with the Palmyrian fold belt behaving like a sinistral, intracontinental-transform fault. It should be noted that the Euphrates volcanic field occurs near the intersection between the Euphrates graben/fault system and the northeastern extension of the Palmyride fold zone. The various volcanic exposures lie in direct continuation of the Palmyride structural zone (Fig. 2b).
Lava flows of the Euphrates volcanic field overlie undeformed Neogene and Quaternary deposits; few old sedimentary sections are exposed in the region. Regional mapping of Dubertret (Reference Dubertret1955) showed the Neogene and Quaternary volcanic rocks of Syria and Lebanon to be fissure-type eruptions. He stated that lava ascended along re-activated as well as Recent fault planes, forming voluminous flows. Ponikarov (Reference Ponikarov1967) also studied the general geology of Syria, which included the area under investigation.
The Euphrates volcanic field includes the Ein Abou-Jema locality, which covers an area of about 30 km2 of mafic–ultramafic lava flows, Himmat el-Jazira (400 km2), Feidat el-Mieza (680 km2) and Kaukab (100 km2). These consist of lava flows and minor pyroclastic flows of Middle and Upper Quaternary age, overlying Neogene sediments. The lavas are mainly composed of phyric varieties (olivine basalts), forming fine-grained massive rocks with variable degrees of vesicularity. The flows form plateaux typical of shield-type volcanoes, and the lower unit in some of the volcanic piles is topped by craters with lava flows and associated pyroclastic materials, reflecting mixed effusive and explosive activity.
More specifically, piles of lava flows at the Kaukab locality are about 10 m thick at the centre, diminishing to 3 m in thickness at the edge of the plateau. The thickness of other piles of lava flows ranges from 7 m to 10 m at the localities of Feidat El-Meiza and Ein Abou Jema, to almost 20 m at Himmat El-Jazira. Based on their position on top of Pleistocene (Lower Quaternary) terraces, Ponikarov (Reference Ponikarov1967) suggested that these lava flows are Upper Quaternary in age. More recent volcanic activity is indicated by the superposition of lava and scoria cones on these flows.
3. Methods of analysis
Whole rock major- and trace-element analyses were conducted at the Activation Laboratories (ACTLABS) in Ancaster, Ontario, Canada. The samples were run for major oxides on a combination of simultaneous and sequential Thermo Jarrell-Ash Enviro II ICP. Calibration was achieved using a combination of CANMET and USGS reference materials. Loss on ignition (LOI) was determined by heating powdered samples for 50 minutes.
For trace and rare earth element (REE) analyses, solutions of the fused samples prepared as above are spiked with three internal standards to cover the entire mass range. These are further diluted and are introduced into a Perkin Elmer-SCIEX Elan 6000 ICP/MS using a proprietary sample introduction methodology. The fusion process ensures total metals (particularly for elements like REE in resistate phases; this may not be the case for acid digestions, particularly for heavy REE and other elements, contained in refractory minerals like zircon, sphene or monazite). The primitive mantle values used for normalization are those of Sun & McDonough (Reference Sun, McDonough, Saunders and Norry1989).
The Sr and Nd isotopic analyses were conducted at McMaster University, Ontario, Canada. Before commencing the chemical procedures, sample powders were leached in 6M HCl for several hours in order to remove any traces of secondary Sr. Sr and Nd were separated in cation exchange columns; the separation was followed by purification in a second separation column according to the method of Richard, Shimizu & Allègre (Reference Richard, Shimizu and Allègre1976). The solution of Nd was mounted on double filaments and measured on a VG Isolab-54 mass spectrometer with multiple collectors in dynamic mode. During the course of the study, the La Jolla Nd-standard gave 143Nd/144Nd = 0.51185. The average precision of repetitions was ± 0.000015 (2σ). The solution of Sr was mounted on a single filament and measured by a VG Isomass–354 with multiple collectors. The Sr isotopes were normalized for fractionation to 86Sr/88Sr = 0.1194, related to a value of 0.710241 for the NBS-987 standard. The average precision of repetitions was ± 0.00002 (2σ).
4. Petrography
Rocks of the Euphrates volcanic field are predominantly phyric, with olivine, clinopyroxene and/or plagioclase phenocrysts embedded in a microcrystalline to cryptocrystalline groundmass. The groundmass consists primarily of plagioclase laths, small grains of clinopyroxene, olivine and opaque phases, in addition to minor amounts of alteration products (mainly iddingsite). These mafic volcanic rocks exhibit a variety of textures, including intersertal, porphyritic, glomeroporphyritic, ophitic, sub-ophitic, pilotaxitic, vesicular and rarely aphyric.
Plagioclase makes up about 45 to 55 vol.% of the whole rock. It forms crystals that range in size from tiny microlites to large prisms, occurring mostly as groundmass material (0.05 to 0.5 mm long microlaths), and rarely as microphenocrysts or phenocrysts (0.6 to 1.2 mm long) which may contain micro-inclusions of apatite needles. Clinopyroxene forms anhedral to subhedral crystals, which often show sector zoning. It varies in colour from pale beige to lilac, forming about 30–35 vol.% of the whole rock. Clinopyroxene is also abundant both as a phenocryst and as a groundmass phase. The proportions of clinopyroxene existing as phenocryst versus groundmass phases are highly variable. Clinopyroxene is rarely twinned, and occasionally contains inclusions of plagioclase microlaths and opaque phases. Olivine crystals occur as subhedral to euhedral, large, equant and occasionally skeletal phenocrysts, and in the form of minute anhedral groundmass grains; olivine forms 15–20 vol.% of the whole rock. Corroded and embayed olivine phenocrysts are observed in several samples. The olivine phenocrysts are largely fresh, but in some rocks the olivine is partially altered to iddingsite. Olivine crystals generally contain inclusions of opaque phases. The latter also occur as small cubes or acicular needles within the groundmass.
5. Geochemistry
5.a. Major, trace and rare-earth element geochemistry
Major and trace element data for 17 representative samples of the Euphrates volcanic field of northeastern Syria are presented in Table 1. Rocks of this volcanic field are largely ultrabasic, as they have very low silica values (38.2 to 45.5 wt% SiO2). They exhibit narrow major element compositional ranges: 9.5–12.9 wt% Al2O3, 8.7–13.0 wt% MgO, 13.0–15.5 wt% Fe2O3 (as total iron), 9.5–11.7 wt% CaO and 2.5–3.5 wt% TiO2 (Table 1). Based on the total alkali versus silica diagram of Le Bas & Streckeisen (Reference Le Bas and Streckeisen1991; not shown), the lavas are represented largely by nephelinites, basanites and alkali basalts. The Mg numbers (= molar Mg/(Mg+Fe2+), assuming a Fe3+/Fe2+ ratio of 0.2) range from 0.61 to 0.68 (with an average of 0.65). These values indicate that the magma may have experienced some degree of olivine and clinopyroxene fractionation.
Fe2O3* is total iron presented as Fe2O3, and Mg no. = (molar Mg/(Mg+Fe2+)) assuming Fe3+/Fe2+ = 0.20.
The Euphrates volcanic rocks also exhibit a relatively narrow trace element compositional range: Cr = 270–392 ppm, V = 203–289 ppm, Sr = 614–1140 ppm, Ba = 175–278 ppm, and Rb = 8–19 ppm, and are generally enriched in the high field strength elements (HFS) such as Zr (133–276 ppm), Y (17–28 ppm) and Nb (25–71 ppm; Table 1). Rocks of the investigated volcanic field exhibit elemental ratios, such as La/Nb (0.77), Zr/Nb (4.6), and Rb/Nb (0.30), similar to the average of some OIB lavas (La/Nb = 0.72, Zr/Nb = 4.1, and Rb/Nb = 0.37; Weaver, Reference Weaver1991).
Variation diagrams of major elements versus Zr (Fig. 3) indicate that P2O5, TiO2, total iron (as Fe2O3) and alkalis increase with increasing Zr, whereas CaO remains somewhat constant. Trace elements exhibit well-defined correlations with Zr; Nb, Y and Sr show a gradual increase with increasing Zr (Fig. 3). Since Zr and Y are incompatible in the main fractionating phases of basaltic magmas (olivine, pyroxene and plagioclase), the Zr/Y ratio is not normally affected by moderate amounts of fractional crystallization. However, as Zr is more incompatible in mantle phases than Y, the Zr/Y ratio tends to be higher when the degree of melting is small; variation of Zr/Y with Zr can thus be used to illustrate petrogenetic processes such as partial melting. The positive correlation of Zr/Y with Zr obtained for the Euphrates lavas (not shown) suggests that partial melting processes have played a significant role in producing the range of magma compositions observed (see Section 6). These trace elements, combined with the geological setting, confirm the within-plate character of the investigated Euphrates lavas.
The rare earth element compositions of seventeen samples from the Euphrates volcanic field are given in Table 2. The rocks are generally enriched in REE relative to primitive mantle, with the sum of REE ranging from 113 to 344 ppm. Primitive mantle-normalized REE patterns are illustrated in Figure 4a. Overall, the REE patterns are parallel to subparallel and strongly fractionated ((La/Yb)N = 19.6). The two nephelinite samples that possess the lowest SiO2 contents (38.2 and 38.6 wt% SiO2) have significantly higher LREE concentrations than the other samples. This is, most likely, the result of differences in the degrees of partial melting of the mantle source rock (see Section 6).
The primitive mantle-normalized incompatible element patterns of the Euphrates volcanic field (Fig. 4b) indicate that the rocks are enriched in the large ion lithophile (LIL) elements (K, Rb, Ba, Th and the LREE) and the incompatible HFS elements Nb, Zr and Ti relative to primitive mantle normalizing abundances. A comparison with trace element profiles of the St Helena alkali basalts (which are typical of HIMU-OIB: Chaffey, Cliff & Wilson, Reference Chaffey, Cliff, Wilson, Saunders and Norry1989) show that the Euphrates lavas are similar in trace-element composition (Fig. 4b). Thus, the overall chemical traits of the investigated lavas reflect strong affinities to OIB. According to Pearce et al. (Reference Pearce, Bender, De Long, Kidd, Low, Güner, Saroğlu, Yilmaz, Moorbath and Mitchell1990), melting of ‘plume’ asthenosphere, or lithosphere enriched by small volume melts from the asthenosphere produces ‘humped’ patterns with Ba, Th and Nb the most enriched elements and the degree of enrichment directly relating to the degree of incompatibility. The marked depletion of Y and the heavy REE in the Euphrates lavas may indicate that garnet was a residual phase during the partial melting event. Thus, the Euphrates lavas can be interpreted in terms of a garnet-bearing mantle source (see Section 6).
5.b. Nd and Sr isotopes
The 143Nd/144Nd ratios of representative samples from the Euphrates volcanic field range from 0.512868 ± 1 to 0.512940 ± 1 (εNd = 4.5 to 5.9), whereas the 87Sr/86Sr ratios range from 0.703090 ± 2 to 0.703520 ± 2 (Table 3). The isotopic compositions are plotted in Figure 5, along with data from other alkali basaltic suites. Also plotted in this diagram (Fig. 5) are the compositions of various mantle reservoirs (EMI, EMII, HIMU and N-MORB), taken from Hart (Reference Hart1988); EMI (‘enriched mantle 1’) and EMII (‘enriched mantle 2’) types of OIB may represent the addition of small amounts of subducted sediments; pelagic in the case of EMI and terrigenous in the case of EMII (Weaver, Reference Weaver1991; Hofmann, Reference Hofmann1997). HIMU (‘high μ’) refers to a high 238U/204Pb (μ) mantle end-member, and has the lowest 87Sr/86Sr of any OIB (Hofmann, Reference Hofmann1997), which is thought to be derived from subducted basaltic oceanic crust. Examples of HIMU-OIB are St Helena, Bouvet, Ascension, Austral Islands, Balleny Islands and the Azores; typical EMI-OIB are Tristan da Cunha, Gough, Kerguelen and Pitcairn; EMII-OIB: Society Islands, Samoa, Tutuila and Upolu (Weaver, Reference Weaver1991; Hofmann, Reference Hofmann1997).
Rocks of the Euphrates volcanic field are, as shown in Figure 5, isotopically similar to some Golan–Galilee alkali basalts (Weinstein et al. Reference Weinstein, Navon, Altherr and Stein2006), Jordanian alkali basalts (Shaw et al. Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003), Lebanese Cenozoic alkali basalts (Abdel-Rahman & Nassar, Reference Abdel-Rahman and Nassar2004), and to the Mesozoic alkali basalts of the Middle East (Laws & Wilson, Reference Laws and Wilson1997; Abdel-Rahman, Reference Abdel-Rahman2002). The investigated lavas also exhibit isotopic compositions similar to the St Helena alkali basalts (Staudigel et al. Reference Staudigel, Zindler, Hart, Leslie and Clague1984; Hofmann, Reference Hofmann1997) as they exhibit relatively high initial 143Nd/144Nd and low initial 87Sr/86Sr isotopic ratios, and are thus distinct from EMI-OIB and EMII-OIB (e.g. Weaver, Reference Weaver1991; Wilson, Reference Wilson1993). It should be noted that the St Helena alkali basalts and the Lebanese lavas were interpreted to represent small-degree partial melts of upper mantle material (Weaver, Reference Weaver1991; Abdel-Rahman & Nassar, Reference Abdel-Rahman and Nassar2004).
6. Discussion
6.a. Role of crustal contamination
On the variation diagram of Ta/Yb versus Th/Yb (after Pearce, Reference Pearce, Hawkesworth and Norry1983), the Euphrates volcanic rocks plot along the diagonal trend of mantle spectrum, and away from the trend of crustal contamination (Fig. 6). This diagram shows also that the examined rocks overlap with some Golan–Galilee alkali basalts (taken from Weinstein et al. Reference Weinstein, Navon, Altherr and Stein2006), but the Lebanese Cenozoic alkali basalts (Abdel-Rahman & Nassar, Reference Abdel-Rahman and Nassar2004) plot slightly off this trend. The degree of contamination that silicate melts may have experienced can be assessed via the use of certain chemical parameters. For example, basaltic lavas affected by crustal contamination exhibit K/P ratios > 7, La/Ta > 22, and La/Nb > 1.5 (e.g. Hart et al. Reference Hart, Wolde, Walter and Mertzman1989). The low values of such elemental ratios in rocks of the Euphrates volcanic field (K/P 3.4, La/Ta 13, La/Nb 0.77, on average), along with their Sr–Nd isotopic composition, and their low silica values (43.1 wt% SiO2, on average), all suggest that the magma was subjected to minimal crustal contamination. Magma ascent may have been rapid enough from the site of partial melting to the surface to escape contamination. As pointed out by Smith et al. (Reference Smith, Sánchez, Walker and Wang1999), the Nd content of most lower crustal xenoliths is too low (usually <10 ppm) to significantly change Nd-isotopic values without adding 70% to 85% lower crustal material. Such a large amount of contamination by crustal material is thermodynamically difficult because a considerable amount of heat is required to assimilate crustal rocks, and the magma would then cool quickly and perhaps ‘freeze’ in place. Moreover, this would have resulted in the presence of some lower crustal xenoliths within the lava flows, but the Euphrates lava flows contain no such crustal xenoliths.
Based on geological and geophysical data, Sawaf et al. (Reference Sawaf, Al-Saad, Gebran, Barazangi, Best and Chaimov1993) showed that a slight crustal thinning is present, going from 37 km south and north of the Euphrates graben to 35 km beneath it. The occurrence of the Euphrates volcanic field within such a tensional stress field is consistent with a lack of significant crustal contamination in these lavas; the lower 87Sr/86Sr ratios (0.70309 to 0.70352) of the Euphrates lavas could represent mantle source values. Thus, the geochemical and field characteristics, along with the nature of tectonism (extensional), suggest that the role of crustal contamination during magma evolution has been minimal.
6.b. Source characteristics and possible depth of melting
In terms of their Zr, Nb and Y compositions (Fig. 7a), the Euphrates lavas plot just above the plume-related mid-ocean ridge basalt (P-MORB), as they exhibit relatively higher concentrations of Zr, but lower concentrations of Y than transitional-, or normal-MORB (T-MORB or N-MORB: Menzies & Kyle, Reference Menzies, Kyle and Menzies1990; Melluso et al. Reference Melluso, Beccaluva, Brotzu, Gregnanin, Gupta, Morbidelli and Traversa1995). The examined rocks overlap with those representing the Golan–Galilee and Jordanian alkali basalts of Weinstein et al. (Reference Weinstein, Navon, Altherr and Stein2006) and Shaw et al. (Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003), respectively, but exhibit higher Zr/Y ratios than the Lebanese Cenozoic alkali basalts of Abdel-Rahman & Nassar (Reference Abdel-Rahman and Nassar2004).
The Euphrates lavas also exhibit elemental ratios ((Zr/Nb) = 4.6, (La/Nb) = 0.77, (Ba/Th) = 61, and (Rb/Nb) = 0.30, on average), similar to those characteristic of OIB (Weaver, Reference Weaver1991). Niobium, along with some other HFS elements, is found to be highly variable in lithospheric mantle melts. Thus, variations in the La/Nb ratios have been interpreted by some authors to reflect the style of metasomatic enrichment (small fraction convection mantle melt, or subduction-related metasomatism: Gibson et al. Reference Gibson, Thompson, Dickin and Leonardos1995). Considering the within-plate character of the Euphrates volcanic field, subduction-related metasomatism was an unlikely process in the region during Cenozoic times. Bradshaw & Smith (Reference Bradshaw and Smith1994) and Smith et al. (Reference Smith, Sánchez, Walker and Wang1999) have suggested that, since HFS elements (such as Nb) are depleted in the lithospheric mantle relative to the light REE (e.g. La), high Nb/La ratios (approximately > 1) indicate an OIB-like asthenospheric mantle source for basaltic magmas, and lower ratios (approximately < 0.5) indicate a lithospheric mantle source. The Nb/La and La/Yb ratios (averages of 1.3 and 27, respectively) are consistent with an asthenospheric mantle (OIB-like) source (Fig. 7b): the Euphrates lavas plot in the field of asthenospheric mantle and near average OIB. Similarly, the Golan–Galilee, the Jordanian and the Lebanese Cenozoic alkali basalts (Weinstein et al. Reference Weinstein, Navon, Altherr and Stein2006; Shaw et al. Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003; Abdel-Rahman & Nassar, Reference Abdel-Rahman and Nassar2004, respectively), plot in the field of asthenospheric mantle. Thus, trace element and isotopic data suggest that the investigated lavas have chemical characteristics similar to OIB (Fig. 7a), and were possibly derived from an asthenospheric mantle source. The relatively uniform 143Nd/144Nd ratios (εNd = 4.5 to 5.9) for the Euphrates volcanic field suggest that its magma was developed by partial melting of an isotopically uniform mantle source.
The presence of garnet as a residual phase in the source region of the lavas can be inferred from: (1) fractionation of the heavy rare-earth elements (HREE; McKenzie & O'Nions, Reference McKenzie and O'Nions1991; cf. Fig. 4a), and (2) the greater than chondritic Dy/Yb ratio (1.57); the Euphrates lavas have an average Dy/Yb ratio of 3.64. These rocks also have an average (Tb/Yb)N ratio of 3.5, which is slightly higher than that of the alkali basalts of Hawaii ((Tb/Yb)N range from 1.9 to 2.5); the Hawaiian basalts are considered to have been derived from a garnet-bearing mantle source (McKenzie & O'Nions, Reference McKenzie and O'Nions1991; Frey et al. Reference Frey, Garcia, Wise, Kennedy, Gurriet and Albarede1991, Reference Frey, Clague, Mahoney and Sinton2000). Further south, the source of the alkali basalts and basanites of the Jordanian Miocene to Recent lavas were related to melts of deep garnet-bearing lithospheric mantle, mixed with asthenospheric mantle melts (Shaw et al. Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003). These basaltic melts were formed in response to a phase of lithospheric extension. It should be noted that Shaw et al. (Reference Shaw, Baker, Menzies, Thirlwall and Ibrahim2003) concluded that the Afar plume of Ethiopia has not been channelled northwestwards beneath the Arabian plate, and played no role in producing the Arabian or Jordanian volcanic fields.
Using the study of Ellam (Reference Ellam1992) on the estimation of depths of melt segregation, curves corresponding to the Ce, Sm and Yb concentrations of the Euphrates lavas (Fig. 7c) indicate a melt segregation depth ranging from about 94 to 108 km (that is, within the garnet lherzolite zone), consistent with an asthenospheric mantle source as indicated in Figure 7b. It should be noted that the transition from garnet to spinel peridotite takes place between a depth of about 60 to 80 km for normal mantle and about 80 to 100 km within hot mantle plumes (McKenzie & O'Nions, Reference McKenzie and O'Nions1991; Lassiter, DePaolo & Mahoney, Reference Lassiter, DePaolo and Mahoney1995).
6.c. Petrogenesis of the Euphrates lavas
It is well established that the geochemical signatures of basaltic magmas rely on three principal factors: (1) the composition of the mantle source including the presence or absence of a volatile phase (H2O + CO2), (2) the percentage of fusion and (3) the pressure conditions (depth) of the system subject to fusion. Thus, the nature of the mantle source material, whether it is dominated by recycled oceanic or continental crust, or by recycled sedimentary components, and the processes associated with melting and melt migration determine the composition of lavas. A number of studies, including those of White (Reference White1985), Allègre et al. (Reference Allègre, Hamelin, Provost and Dupré1987), Hart (Reference Hart1988), Weaver (Reference Weaver1991), Gibson et al. (Reference Gibson, Thompson, Weska, Dickin and Leonardos1997), Abdel-Rahman & Kumarapeli (Reference Abdel-Rahman and Kumarapeli1999), Frey et al. (Reference Frey, Clague, Mahoney and Sinton2000), Gibb & Henderson (Reference Gibb and Henderson2006), Furman et al. (Reference Furman, Kaleta, Bryce and Hanan2006) and Weinstein et al. (Reference Weinstein, Navon, Altherr and Stein2006), reflect that mafic alkaline lavas are diverse geochemically and form by partial melting of a range of mantle sources. In order to assess the role of petrogenetic processes such as partial melting, modelling was performed (using the batch melting equations of Shaw, Reference Shaw1970 and Allègre & Minister, Reference Allègre and Minster1978), and the calculations were done using two model source compositions: (1) primitive mantle taken from Sun & McDonough (Reference Sun, McDonough, Saunders and Norry1989) and (2) a mixed (50% primitive–50% depleted) mantle source of McKenzie & O'Nions (Reference McKenzie and O'Nions1991). Spinel, garnet and clinopyroxene were assumed to decrease in abundance linearly with increasing degrees of partial melting, as they are typically consumed at less than 25% partial melting (McKenzie & O'Nions, Reference McKenzie and O'Nions1991; Lassiter, DePaolo & Mahoney, Reference Lassiter, DePaolo and Mahoney1995).
Model and melting proportions are given in Tables 4 and 5 and are in line with those used in other partial melting calculations (e.g. Hanson, Reference Hanson1980; McKenzie & O'Nions, Reference McKenzie and O'Nions1991; Witt-Eickschen & Kramm, Reference Witt-Eickschen and Kramm1997; Abdel-Rahman, Reference Abdel-Rahman2002). Modelling was performed using three different mantle mineral assemblages: spinel lherzolite, garnet lherzolite and spinel–garnet lherzolite, for both a primitive and a mixed source composition. The partition coefficients used are from McKenzie & O'Nions (Reference McKenzie and O'Nions1991). Partial melting calculations were performed for 1%, 2% and 3% partial melting. The results of the modelling (Tables 4, 5; Fig. 8) show that melting of a spinel-bearing source overestimates the HREE, and melting of a mixed source yields much lower LREE concentrations. Thus, neither depleted nor mixed primitive/depleted mantle material represents the mantle source for the Euphrates lavas and garnet is a required phase, but possibly with some minor spinel. The REE pattern of the calculated liquid produced by 1% batch partial melting of garnet lherzolite (of a primitive mantle composition) produces a reasonable fit, as it closely matches that of the average Euphrates volcanic rocks (Fig. 8).
See notes for Table 5.
Calculated melts produced by 1%, 2%, and 3% batch partial melting are no. 1 to no. 18, and no. 19 is the measured, average concentration of the Euphrates lavas. The starting mode, melt mode and mantle source type used to produce each of the calculated melts are as follows:
Melt no. 1–3; starting mode a, melt mode a, mixed source, for 1, 2, and 3% melting, respectively,
no. 4–6; starting mode b, melt mode b, mixed source, for 1, 2, and 3% melting, respectively,
no. 7–9; starting mode c, melt mode c, mixed source, for 1, 2, and 3% melting, respectively,
Melt no. 10–12; starting mode a, melt mode a, primitive source, for 1, 2, and 3% melting, respectively,
no. 13–15; starting mode b, melt mode b, primitive source, for 1, 2, and 3% melting, respectively,
no. 16–18; starting mode c, melt mode c, primitive source, for 3, 2, and 1% melting, respectively.
The composition of the calculated melt no. 18 (produced by 1% melting of garnet lherzolite of a primitive mantle source) closely matches that of the measured average composition of the Euphrates lavas (no. 19). See text for details.
Another modelling scenario has also been considered. Similar to the model suggested by Pearce et al. (Reference Pearce, Bender, De Long, Kidd, Low, Güner, Saroğlu, Yilmaz, Moorbath and Mitchell1990) for the Karacadağ volcanic plateau of southern Turkey, this modelling scenario is based on the assumption that the source is a mantle lithosphere, which was enriched in incompatible trace elements by metasomatism. The metasomatizing fluids could derive from either upwelling of a very deep primitive source such as garnet lherzolite, or by lateral migration from a MORB-like depleted source. The percolation of such fluids through the Euphrates lithospheric mantle source may have led to its enrichment in trace elements. With no geophysical or structural evidence of upwelling, the metasomatizing fluids may have been generated from a similarly depleted source by a very limited percentage of partial melting, which would be needed to reproduce the observed enrichment in incompatible elements. Thus, for this modelling attempt, the calculations were performed for 0.1% to 10% partial melting, using a variety of mantle source materials with variable mineralogies, and a number of different scenarios involving metasomatized sources. For the chemical compositions of mantle sources and the partition coefficients of the various elements, we have used values of McKenzie & O'Nions (Reference McKenzie and O'Nions1991). Model parameters, melting proportions and results are given in Table 6 and Figure 9. Within the various modelling calculations made, a match was produced as a result of a small degree of partial melting (2% melting; Fig. 9) of a garnet-bearing depleted source (garnet peridotite) enriched with 5% addition of metasomatizing fluids derived from a similar depleted source by 0.5% melting.
Source rocks Gt–Lhr, Gt–Sp–Lhr, and Sp–Lhr represent garnet lherzolite, garnet–spinel lherzolite and spinel lherzolite, respectively. DS1 shows the chemical composition of the garnet-bearing depleted source, MDS shows the chemical composition of fluids produced by 0.5% fusion of the depleted source, and DM5 shows the chemical composition of the garnet-bearing depleted mantle source with 5% addition of metasomatizing fluids derived from a similar depleted source by 0.5% fusion. The calculated melts listed here and shown in Figure 9 are produced by 10% to 1% partial melting of this metasomatized garnet-bearing depleted source. The average composition of the Euphrates lavas (Euphr.) closely resembles that of the calculated liquid produced by 2% partial melting of the metasomatized garnet-bearing depleted mantle source (see text for details).
However, the overall chemical and petrological characteristics are more consistent with the generation of the Euphrates magma by a small degree of partial melting (F = 1%) of a fertile, garnet-lherzolite, asthenospheric mantle source (cf. Fig. 8), possibly containing a minor spinel component.
6.d. The geodynamic framework
The Euphrates volcanic field is located near the Bitlis collision suture, which represents a major plate boundary (cf. Fig. 2). This tectonic feature was created by, and is continuously being reshaped as a result of, movements of both the African and the Arabian plates with respect to the Eurasian plate. The lithospheric mantle of the northern border of the Arabian plate was deformed during and after the collision which led to the formation of the Bitlis suture (current tectonic regime). This could have provoked the development of shear zones; the equivalent at great depth of detached or transform faults observed in the underlying continental base. According to Spera (Reference Spera, Menzies and Hawkesworth1987), the rise of magmas by such a network of deformation and fractures could explain the extrusion to the surface (the eruption) of alkaline magmas. Chemical characteristics of the Plio-Quaternary Euphrates lavas reflect their within-plate nature, consistent with magma generation in an extensional tectonic regime.
As indicated above, the Euphrates graben/fault system was developed by extension/transtension following a period of extensive uplift and erosion during the Jurassic and Early Cretaceous (Brew et al. Reference Brew, Barazangi, Al-Maleh and Sawaf2001). This graben is considered to represent manifestation of re-activated pre-Cenozoic structures that responded to forces acting along nearby Arabian plate boundaries (Cenozoic collision along the Bitlis suture: Sawaf et al. Reference Sawaf, Al-Saad, Gebran, Barazangi, Best and Chaimov1993). The graben becomes more pronounced in the complex region where it intersects the NE-trending Palmyride mountain belt northwest of Der Azzour. Based on seismic reflection data, Sawaf et al. (Reference Sawaf, Al-Saad, Gebran, Barazangi, Best and Chaimov1993) showed that Quaternary volcanism produced basaltic flows (including the Euphrates volcanic field) above some of the deeper penetrating faults (Palmyride trend) and along a northwestern (Euphrates) fault system where it cuts the western side of the Derro high, also indicating a deeply penetrating fault.
Sawaf et al. (Reference Sawaf, Al-Saad, Gebran, Barazangi, Best and Chaimov1993) also demonstrated that a slight crustal thinning, from 37 km south and north of the Euphrates depression to 35 km beneath it, can account for the subtle change in the long-wavelength Bouguer gravity anomaly over it. The geophysical study of Brew et al. (Reference Brew, Litak, Seber, Barazangi, Al-Iman and Sawaf1997) showed that basement depth on the north side of the Euphrates is about 6 km, while south of it the basement depth is at least 8.5 km, and those authors interpreted this difference in basement depth on adjacent sides of the graben system as an indication that the Euphrates system is a suture/shear zone. This suture may represent a zone of weakness between distinct crustal blocks, possibly inherited from Late Proterozoic accretion of the Arabian plate. Thus, the nature of the Euphrates volcanic rocks, combined with the regional geological context and the overall tectonic framework, suggest that volcanism has occurred in a localized tensional regime.
In their study of volcanism in East Anatolia, Pearce et al. (Reference Pearce, Bender, De Long, Kidd, Low, Güner, Saroğlu, Yilmaz, Moorbath and Mitchell1990) indicated that perturbation of the thickened lithosphere by delamination of the thermal boundary layer, perhaps coupled with local stretching associated with pull-apart basins on strike-slip fault systems, is sufficient to generate melt. Litak et al. (Reference Litak, Barazangi, Beauchamp, Seber, Brew, Sawaf and Al-Youssef1997) demonstrated that Neogene continental collision along the northern and eastern Arabian plate boundaries caused minor reactivation of the Euphrates fault system, in harmony with significant inversion and the main phase of uplift of the nearby Palmyride and Sinjar mountains. In their study of the Neogene volcanic activity of western Syria, Lustrino & Sharkov (Reference Lustrino and Sharkov2006) proposed that the transition from strongly compressive to transtensional stresses may have allowed for partial melting as a consequence of mantle decompression. In the Euphrates region, heat necessary for magma generation may have been the result of adiabatic decompression and melting. Reactivation of deep-seated fractures beneath the Euphrates graben was interpreted to have resulted as a consequence of Neogene collision of the Arabian plate with Eurasia along the Bitlis suture (Litak et al. Reference Litak, Barazangi, Beauchamp, Seber, Brew, Sawaf and Al-Youssef1997). The lavas may have penetrated the crust along such reactivated fractures beneath the Euphrates basin.
7. Conclusions
(1) The Plio-Quaternary Euphrates volcanic field of northeastern Syria occurs along the northern extension of the Palmyride zone at its intersection with the Euphrates graben. This volcanic field represents the northern segment of the Cenozoic volcanic province of the Middle East and is located near the collision (Bitlis) suture of the Arabian and the Eurasian plates. The lavas form relatively thick discontinuous successions. They are mostly phyric, consisting of about 45–55 vol.% plagioclase, 30–35% clinopyroxene, 15–20% olivine and opaque phases.
(2) Rocks of the Euphrates volcanic field are mostly basanites and alkali basalts, have a narrow range of major element compositions (SiO2, 38.2–45.5 wt%; MgO, 8.7–13.0 wt%), are alkaline in nature, and are enriched in Ti (2.5–3.5 wt% TiO2), Zr (133–276 ppm), Nb (25–71 ppm) and Y (17–28 ppm). The primitive mantle-normalized patterns are strongly fractionated ((La/Yb)N = 19.6) and conformable. The 143Nd/144Nd isotopic composition of the Euphrates lavas ranges from 0.512868 to 0.512940 and 87Sr/86Sr from 0.70309 to 0.70352. The overall chemical and isotopic compositions of these lavas reflect strong affinities to OIB. The relatively uniform 143Nd/144Nd (εNd = 4.5 to 5.9) for the Euphrates volcanic field suggests that its magma was developed by partial melting of an isotopically uniform mantle source.
(3) Elemental ratios such as K/P (3.4, on average), La/Ta (13) and La/Nb (0.77), and the very low silica values (43.1 wt% SiO2), suggest that crustal contamination did not play a significant role during magma evolution; the magmas probably experienced very rapid ascent.
(4) Petrogenetic modelling suggests that the magma could have been produced as a result of a small degree of partial melting of either (1) a garnet-bearing depleted source enriched with a small addition of metasomatizing fluids, or (2) a garnet-bearing fertile source. The overall chemical and petrological characteristics are more consistent with the generation of the Euphrates magma by a small degree of partial melting (F = 1%) of a primitive, garnet-lherzolite mantle source, possibly containing a minor spinel component. The Euphrates volcanic field displays the geochemical characteristics of within-plate lavas. Heat necessary for magma generation may have been the result of adiabatic decompression and melting. The Neogene collision of the Arabian plate with Eurasia along the Bitlis suture resulted in reactivation (beneath the Euphrates basin) of deep-seated fractures, along which lavas may have penetrated the crust.
Acknowledgements
We would like to note that the collaboration of Roger Laurent (Laval University, Québec) was essential in this study. We thank Alan Dickin at MacMaster University, Ontario, for isotope analysis work. We appreciate the help of A. A. Turkmani, Ministry of Petroleum and Mineral Resources, Syria, during the field work, as well as that from Muawia Barazangi and Graham Brew at Cornell University and Michel Delaloye of the University of Geneva. We acknowledge the technical support provided by M. Ijreiss of the American University of Beirut (Lebanon). This work was part of a larger research project made possible by grants from the Fonds pour la formation de chercheurs et l'aide à la recherche (Québec) and the Conseil de recherches en sciences humaines du Canada. Valuable comments made by the two referees (J. Baker and T. Furman) improved our contribution.