Skip Navigation


Journal of Petrology Advance Access originally published online on March 7, 2006
Journal of Petrology 2006 47(5):1017-1050; doi:10.1093/petrology/egl003
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Appendix 1
Right arrow All Versions of this Article:
47/5/1017    most recent
egl003v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (5)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by WEINSTEIN, Y.
Right arrow Articles by STEIN, M.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2006. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oxfordjournals.org

The Role of Lithospheric Mantle Heterogeneity in the Generation of Plio-Pleistocene Alkali Basaltic Suites from NW Harrat Ash Shaam (Israel)

Y. WEINSTEIN1,2,*, O. NAVON3, R. ALTHERR4 and M. STEIN3,{dagger}

1 DEPARTMENT OF GEOGRAPHY, BAR-ILAN UNIVERSITY, RAMAT-GAN 52900, ISRAEL
2 DEPARTMENT OF GEOCHEMISTRY, CENTER OF EARTH SCIENCE, UNIVERSITY OF GÖTTINGEN, 37077 GÖTTINGEN, GERMANY
3 INSTITUTE OF EARTH SCIENCES, THE HEBREW UNIVERSITY OF JERUSALEM, JERUSALEM 91904, ISRAEL
4 INSTITUTE OF MINERALOGY, UNIVERSITY OF HEIDELBERG, D-69120 HEIDELERG, GERMANY

RECEIVED AUGUST 9, 2004; ACCEPTED JANUARY 11, 2006


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Plio-Pleistocene volcanism in the Golan and Galilee (northeastern Israel) shows systematic variability with time and location: alkali basalts were erupted in the south during the Early Pliocene, whereas enriched basanitic lavas erupted in the north during the Late Pliocene (Galilee) and Pleistocene (Golan). The basalts show positive correlations in plots of ratios of highly to moderately incompatible elements versus the concentration of the highly incompatible element (e.g. Nb/Zr vs Nb, La/Sm vs La) and in diagrams of REE/HFSE (rare earth elements/high field strength elements) vs REE concentration (e.g. La/Nb vs La). Some of these correlations are not linear but upward convex. 87Sr/86Sr ratios vary between 0·7031 and 0·7034 and correlate negatively with incompatible element concentrations and positively with Rb/Sr ratios. We interpret these observations as an indication that the main control on magma composition is binary mixing of melts derived from two end-member mantle source components. Based on the high Sr/Ba ratios and negative Rb anomalies in primitive mantle normalized trace element diagrams and the moderate slopes of MREE–HREE (middle REE–heavy REE) in chondrite-normalized diagrams, we suggest that the source for the alkali basaltic end-member was a garnet-bearing amphibole peridotite that had experienced partial dehydration. The very high incompatible element concentrations, low K content, very low Rb contents and steep MREE–HREE patterns in the basanites are attributed to derivation from amphibole- and garnet-bearing pyroxenite veins. It is suggested that the veins were produced via partial melting of amphibole peridotites, followed by complete solidification and dehydration that effectively removed Rb and K. The requirement for the presence of amphibole limits both sources to lithospheric depths. The spatial geochemical variability of the basalts indicates that the lithosphere beneath the region is heterogeneous, composed of vein-rich and vein-poor domains. The relatively uniform 143Nd/144Nd ({varepsilon}Nd = 4·0–5·2) suggests that the two mantle sources were formed by dehydration and partial melting of an originally isotopically uniform reservoir, probably as a result of a Paleozoic thermal event.

KEY WORDS: basanites; lithospheric heterogeneity; magma mixing; amphibole peridotite; pyroxenites


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Alkali basalts form a small part of the global basaltic inventory. However, the association with continental rift environments and the rich assemblages of mantle xenoliths they carry make them an important source of information on the composition of the uppermost mantle beneath continental areas and on the mechanisms of basaltic magma generation. A longstanding discussion concerns the question of whether continental alkali basalts are derived from asthenospheric or lithospheric mantle sources or a mixture of both (e.g. Altherr et al., 1990Go; Wilson & Downes, 1991Go; Stein & Hofmann, 1992Go; Stein et al., 1997Go; Shaw et al., 2003Go). Another discussion concerns the mechanism of melting, either a passive, extension-related scenario (e.g. Chapin, 1979Go; Bohannon et al., 1989Go; Stein & Hofmann, 1992Go; Bradshaw et al., 1993Go; Shaw et al., 2003Go), or an active scenario, which involves a thermal anomaly, usually a mantle plume (e.g. Gass, 1970Go; White & McKenzie, 1989Go; Camp & Roobol, 1992Go; Wilson & Patterson, 2001Go). A third possibility, which involves melting of the lithosphere by heat conduction from a sub-lithospheric source, was discussed in the case of small-volume continental flood basalts (Turner et al., 1996Go; following Gallagher & Hawkesworth, 1992Go), but is also relevant to the generation of alkali basalts, as discussed below.

A third question, which is closely related to the above two, concerns the involvement of fluids in the generation of alkali basalts. Hydrous fluids are known to be important in the subduction environment, where they control the degree of melting (e.g. Stolper & Newman, 1994Go). There is growing evidence that hydrous minerals are also a common constituent in the source rocks of continental alkaline lavas (e.g. Francis & Ludden, 1995Go), but the role of ‘wet melting’ in the petrogenesis of alkali basalts has still not been well established. For instance, the origin of both East African and western Arabian basalts is mainly discussed in terms of melting a dry peridotitic source (e.g. White & McKenzie, 1989Go; Camp & Roobol, 1992Go). Stein et al. (1997)Go were the first to consider the possibility that the presence of hydrous minerals (i.e. amphibole) in the lithospheric mantle, metasomatic remnants of Late Proterozoic subduction, may actually be the main control on partial melting beneath Arabia. Similarly, Späth et al. (2001)Go considered dehydration of a metasomatized lithosphere as the cause for melting under the Chyulu Hills, southern Kenya.

In this study we report the temporal and spatial variations in bulk-rock chemical and Nd–Sr isotopic composition of Plio-Pleistocene alkali basalts from the northwestern part of Harrat ash Shaam, one of the largest volcanic provinces in the Arabian plate (Fig. 1). We show that most of the variability in the trace element composition of these basalts cannot be accommodated by simple models of variable degree melting of a homogeneous mantle source. Instead, we attribute this variability to mixing of melts derived from partial melting of pyroxenite veins and their hydrated host peridotites. We also suggest that the veins were produced during earlier magmatic processes within the lithosphere. Finally, we discuss the implication of the model for the melting regime beneath northern Arabia.


Figure 1
View larger version (50K):
[in this window]
[in a new window]
 
Fig. 1. Map of the volcanism in northern Israel (coordinates are the new Israel grid). The ages of the volcanic rocks are based on K–Ar ages mainly determined by Heimann (1990)Go, Shaliv (1991)Go, Mor (1993)Go, Heimann et al. (1996)Go and during the current study (the central Golan, Table 1). Inset map shows the distribution of Cenozoic magmatism in the western Arabian peninsula.

 
Geological background
This study focuses on basalts that erupted during the Plio-Pleistocene in the Golan and the Galilee areas of northern Israel, which are located at the northwestern part of Harrat ash Shaam, the largest of about 20 Cenozoic volcanic fields in the western Arabian Peninsula (Fig. 1). The volcanism in Arabia is Oligocene to Holocene in age (31 Ma to recent; Menzies et al., 1997Go; although ages as old as Middle Eocene have also been reported, e.g. Pallister, 1987Go), and is associated with rifting along the Red Sea. Volcanism in the Harrat ash Shaam commenced at about 25 Ma (Barberi et al., 1980Go; Ilani et al., 2001Go) and continued to recent times (Dubertret, 1955Go). At its northwestern edge, the Ash Shaam volcanic field is bisected by the Dead Sea Transform fault (Fig. 1), which separates the volcanism of the Galilee to the west from that of the Golan to the east. Volcanism in this area occurred in two main periods: Middle Miocene (17–9 Ma) and Plio-Pleistocene (5·5 Ma to recent, Mor, 1993Go; Shaliv, 1991Go). The center of the Middle Miocene volcanism is located in the Yizre'el Valley and southeastern Galilee (Fig. 1, Shaliv, 1991Go; Weinstein, 2000Go). The Plio-Pleistocene volcanism includes basaltic flows and cinder cones in the Golan and the eastern Galilee as well as basalts along the Dead Sea Transform. The latter are mostly known from subsurface drilling and are intercalated with lacustrine sediments (Mittlefehldt & Slager, 1986Go). This study focuses on the Plio-Pleistocene basalts from the Golan and Galilee.

Overall in the study area, there is a trend of decreasing age from south to north (e.g. Mor, 1993Go; Heimann et al., 1996Go; Fig. 1). West of the Dead Sea Transform, the focus of the Middle Miocene volcanism was in the Yizre'el Valley, locally reaching a thickness of more than 800 m. Early Pliocene activity was concentrated in eastern Lower Galilee (west of the Sea of Galilee), and Late Pliocene lavas erupted further north, in the Upper Galilee (Fig. 1). The thickness of the Galilee Pliocene volcanic sections usually does not exceed 100 m. East of the Dead Sea Transform, the southern Golan plateau is covered by a 20–200 m section of Early Pliocene basalts (Fig. 1). Further north, in the Yahudiya plateau of the central western Golan, there is an ~150 m thick Late Pliocene basaltic section. This section is probably underlain in the subsurface by Early Pliocene basalts, as indicated by a 3·3 Ma basaltic flow in the bottom of one of the local canyons (sample Wy-292, Table 1) and from seismic reflection profiles (Dafny et al., 2003Go). Pleistocene basalts are found only in the eastern and northern Golan (Fig. 1). In most of the eastern and northern Golan, the volcanic section is relatively thick (up to 700 m in the Bashanit area; Dafny et al., 2003Go; Fig. 1), but the age of the subsurface volcanic rocks is not known (there is one small exposure of Early Pliocene basalts in the northernmost part of the Golan; Heimann et al., 1996Go). A similar decreasing age pattern is observed along the Dead Sea Transform. Basalts beneath and around the Sea of Galilee are of Early Pliocene age (Heimann, 1990Go; Heimann et al., 1996Go; Hurwitz et al., 2002Go), whereas Late Pliocene and Early Pleistocene basalts are found only further north, in northern Korazim and in the subsurface and around the Hula Valley (Fig. 1, Heimann, 1990Go; Mor, 1993Go).


View this table:
[in this window]
[in a new window]
 
Table 1: Geochemistry of representative Plio-Pleistocene Golan–Galilee basalts

 
The main active tectonic features in the region are the Red Sea rift and the Dead Sea Transform fault (inset to Fig. 1). The Red Sea trends N25W. Rifting started at 26 Ma, accompanied (or pre-dated) by uplift and volcanism on the Arabian side (e.g. Omar & Steckler, 1995Go; Menzies et al., 1997Go). Seafloor spreading commenced at about 5 Ma along the southern part of the axial trough. The Dead Sea Transform fault separates the Arabian plate, which includes the Golan, from the Sinai sub-plate including the Yizre'el Valley and the Galilee. The fault trends north–south and the total amount of left-lateral slip along it is about 105 km (Quennell, 1958Go; Freund et al., 1970Go; Freund, 1978Go; Garfunkel et al., 1981Go). It is commonly accepted that this displacement has accumulated during the past 20 Myr (Bartov et al., 1980Go) and that an offset of 60 km had already occurred by the end of the Miocene (5 Ma), before the Early Pliocene volcanism started (Freund, 1978Go; Garfunkel et al., 1981Go). A series of pull-apart basins (e.g. the Dead Sea, the Sea of Galilee and the Hula Valley, Fig. 1) was developed along the Dead Sea Transform as a result of sideways stepping of the major fault planes (see Garfunkel, 1981Go; Garfunkel & Ben Avraham, 2001Go). Pliocene and Plio-Pleistocene basalts are buried in the Sea of Galilee and the Hula Valley basins, respectively, under hundreds to thousands of meters of fluvio-lacustrine sediments.


    SAMPLING AND ANALYTICAL METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Ninety-seven samples were collected from volcanic outcrops in the Golan and the Galilee, representing the entire Plio-Pleistocene sequence on both sides of the Dead Sea Transform. About 100 g of fresh rock pieces were crushed and ground to fine powder (<200 mesh) for chemical and isotopic analyses, and thin sections were prepared for petrographic studies.

Concentrations of major elements were determined by X-ray fluorescence (XRF) using a wavelength Philips PW-1404 system at Guttenberg University, Mainz, using glass discs. Calibration curves were established using international standards [recommended values from Govindaraju (1989)Go and Jochum et al. (1990)Go]. Fe was determined as Fe2O3(total); Fe2+ was assumed to be 90% of the total iron. Accuracy was within 1–2% for most elements (accuracies of Na and Mn are 4% and 5%, respectively). Reproducibility was measured by running duplicates within run and in separate runs and was found to be better than 2% for most elements and 4% for Na. In 12 samples, concentrations of major elements were determined by inductively coupled plasma atomic emission spectrometry (ICP-AES; Weinstein et al., 1994Go), and three of these samples were remeasured by XRF. The two methods agree within a few per cent for most of the elements, except for P concentrations, which were consistently higher (by 6–9%) in the ICP-AES measurements. For these three samples, we use the XRF data. Loss on ignition (LOI) was determined by heating 1 g of a sample to about 1100°C, and was corrected for oxidation of ferrous iron during heating. Calculated LOI values are considered as maximum values, as many basalts are relatively oxidized.

Concentrations of trace elements (V, Cr, Co, Ni, Cu, Zn, Ga, Rb, Sr, Y, Zr, Nb, Ba) were also determined by XRF using powder pellets. Accuracy is usually better than 10%relative. Ba concentrations measured on one of the standards (BHVO-1) were very different from the recommended values (86 and 139 ppm, respectively), whereas six other standards yielded accuracies better than 10%. This suggests that Ba data should be treated cautiously. Precision of trace elements was also better than 10% relative, except for Ba, for which precision was 16%.

Concentrations of Sc, Cs, Hf, Ta, Pb, Th, U and 14 rare earth elements (REE) in 40 samples were determined by inductively coupled plasma mass spectrometry (ICP-MS) using a VG Plasma Quad II+ system (University of Göttingen) with a peak jump mode [following the procedure described by Gao et al. (1999)Go]. Blanks are negligible for most of the elements, except for Gd and Dy (in some of the runs), which are, therefore, not used below for quantitative modeling. Accuracies are better than 10% for most REE and U (for Tb and Tm accuracy is 20%), whereas those for Ta and Hf are 15% (low stability of HFSE in the solution) and those for Pb and Th are 25% and 30%, respectively (a strong memory effect for Th). Therefore, both the Pb and Th data should be interpreted cautiously. Sc was consistently higher (by about 20%) than the recommended values for international standards as a result of interference with HCO2. Cs was lower by 20%. Therefore, both Sc and Cs measured by ICP-MS are not reported in this study.

The concentrations of Sc, Cs, Hf, Ta, Th, U and 11 REE in 15 samples were determined by instrumental neutron activation analysis (INAA) in the University of Karlsruhe, following the procedure described by Kramar & Puchelt (1982)Go. Accuracies are between 0 and 15%. Comparative analyses of six samples (three of them from the current study) show that there is a good agreement between INAA and ICP-MS measurements for La, Ce, Nd, Sm, Eu, Gd, Yb, Lu, Th and U (<20% difference), whereas Tb, Ho, Ta and Hf are in significant disagreement [for the last two this is probably because of the problem of high field strength element (HFSE) stabilization in solutions during ICP-MS analyses]. For the three comparative samples, the reported data (Table 1) are those from the INAA (except for the Pb and Pr data, which are from the ICP-MS measurements).

Rock powder (100 mg) was leached for about 40 min in warm 0·5N HCl to remove secondary Sr. Leaching experiments on some of the samples showed that the main difference in the Sr isotope ratios of the treated basalts is between no leaching and 0·5N HCl, whereas the effect of changing to stronger acid (e.g. 6N HCl, Table 1) was within the analytical error of the laboratory (an average difference of 2·3 x 10–5). After washing and drying, the samples were digested by 3 ml of 24N HF + 0·5 ml of 14N HNO3 for at least 8 h at 120°C, and evaporated at 140°C (some of the samples were heated for an additional 2 h at 180°C to break the fluorides). The evaporates were then dissolved in diluted HNO3 (120°C), checked for complete dissolution, evaporated and redissolved by 6N HCl (70–80°C). After evaporation to dryness, the residue was taken up by 3·2 ml of 2·6N HCl and centrifuged in Teflon tubes. Solutions (2 ml) were loaded on columns consisting of ~15 ml cation exchange (preconditioned by 20 ml of 2·6N HCl). The columns were then washed by additional 42 ml of 2·6N HCl (eluting most of the elements except for Sr and the REE), and then most of the Sr was eluted from the column by another 6 ml of 2·6N HCl. After 10 ml wash (14 ml in the last few samples) by 6N HCl, the REE were eluted by an additional 6 ml of 6N HCl. After drying, both Sr and REE aliquots were treated with HNO3 to destroy residual resin and dried again. REE aliquots were taken by 0·5 ml of 0·18N HCl and loaded on columns consisting of 0·5 g orthophosphoric acid (1·5N) mixed with 5 g Teflon beads (after conditioning the columns with 0·18N HCl). The columns were washed with an additional 7·5 ml of 0·18N HCl, and the Nd was eluted with an additional 6 ml and dried. All chemical procedures were carried out in a clean environment, and all the acids used were double distilled. Blanks were negligible. Isotopic analyses were carried out on a Finnigan MAT 262 mass spectrometer at the Geochemistry Institute, University of Göttingen, equipped with nine Faraday cups and using a multicollector static mode for all elements. For Sr isotope analysis, about 4 ng were loaded with phosphoric acid onto Ta single filaments (except for the standard NBS-987, which was loaded onto Re filaments). Sr analyses were corrected for mass fractionation by normalizing to 88Sr/86Sr = 8·375209. The NBS-987 standard was routinely measured, giving a mean 87Sr/86Sr value of 0·710229 and a reproducibility (2{sigma}) of 18 x 10–6 (n = 13). The mean of all NBS-987 measured in the laboratory during the period September 1995 to July 1996 was 0·710244 (2{sigma} = 32 x 10–6). For Nd isotope analysis, about 1 ng was loaded onto Re double filaments. Isotope ratios were normalized to 146Nd/144Nd = 0·7219. The measured mean 143Nd/144Nd for the La Jolla Nd standard was 0·511848 with a reproducibility of 16 x 10–6 (2{sigma}, n = 11), whereas the mean for the period September 1995 to July 1996 was 0·511851 (2{sigma} = 18 x 10–6). All isotopic ratios measured (Sr, Nd) are considered as initial ratios because of the relatively young ages (0–5 Ma) of the basalts studied and their relatively low Rb/Sr ratios.

K–Ar ages were determined by A. Heimann (Geological Survey of Israel). Details of the analytical procedures have been given by Heimann et al. (1996)Go.


    PETROGRAPHY AND CHEMISTRY OF THE GOLAN–GALILEE BASALTS
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
General description and petrography of the Golan–Galilee basalts
The volcanic rocks of the Galilee and the Golan are mainly composed of basaltic flows. Scoriaceous rocks are mostly restricted to the young cinder cones of the northern Golan (Pleistocene) and of the Upper Galilee (Late Pliocene), although a few scoria exposures are also found amongst the Early Pliocene rocks of the southern Golan. Tuff deposits are very rare, and are not included in this study.

We subdivide the Plio-Pleistocene Golan–Galilee volcanism into six units according to age and location (Fig. 1). These are, in order: Early Pliocene basalts from the southeastern (Lower) Galilee, Late Pliocene basalts from the northeastern (Upper) Galilee, Early Pliocene basalts from the southern Golan, Late Pliocene basalts from the western Golan, Early Pleistocene basalts from the eastern and northern Golan, and Late Pleistocene basalts and scoria also from the eastern and northern Golan. In most figures, we refer to both Golan Pleistocene units as one group.

Most Golan–Galilee samples typically show porphyritic to glomeroporphyritic textures, with the phenocrysts usually dominated by olivine (Fo78 to Fo85, 5–20% of the rock; Table 1), which is often partly or completely replaced by iddingsite. Clinopyroxene phenocrysts occur in smaller proportions (<5%), and they are rare in the Early and Late Pliocene basalts. Early Pleistocene basalts from the Golan are the only group that shows similar proportions of olivine and clinopyroxene phenocrysts (2–3%). Plagioclase phenocrysts (labradorite, usually <1%) are observed in about 20% of the samples, mainly those of Early Pliocene age. The groundmass of all Early Pliocene and of Golan Late Pliocene basalts is mostly medium grained, whereas that of Galilee Late Pliocene and of Golan Pleistocene basalts is usually fine grained, with a pilotaxitic texture dominated by plagioclase laths and prismatic crystals of clinopyroxene. Intergranular nepheline is also common. The groundmass of the scoriaceous rocks usually shows a hyalopilitic texture, although the glass is usually replaced. Calcite amygdules are common in the Pliocene basalts.

Mantle xenoliths are common in the Golan and Galilee basalts, though usually restricted to the basanitic units. Spinel peridotites are found both in the Golan Late Pleistocene and in the Upper Galilee Late Pliocene basanites. Pyroxenites are found just in one Pleistocene maar from the northern Golan (Birket Ram, Fig. 1). The xenolith assemblage of this maar is especially large and variable, containing peridotites, dry and hydrous pyroxenites, mafic granulites and megacrysts of kaersutite and microcline (Mittlefehldt, 1984Go; Stein, 1987Go). Mafic granulites (as well as anorthoclase megacrysts) were also found in an Early Pliocene Lower Galilee site (Gazit, 2005Go).

Major and trace element geochemistry of the basalts
The bulk-rock major and trace element compositions of representative Golan–Galilee basalts are given in Table 1. The complete dataset is available for downloading from the Journal of Petrology website at http://www.petrology.oupjournals.org. The samples mainly plot in the basanite, hawaiite and alkali basalt fields of the total alkalis vs silica (TAS) diagram (Fig. 2), and show a trend of increasing total alkali concentration (Na2O + K2O) with decreasing silica. SiO2 concentrations range between 43% and 51% and total alkalis between 3·5% and 7·5%. The spectrum of compositions (alkali basalt–hawaiite–basanite) is similar to that of the Neogene–Quaternary basalts from northeastern Jordan (another part of the Harrat ash Shaam volcanic field; Shaw et al., 2003Go), although the Golan–Galilee basalts are slightly higher in SiO2 (Fig. 2). Miocene basalts from the Yizre'el Valley and the southeastern Galilee (Fig. 1) are somewhat different. These extend to lower SiO2 concentrations than the Plio-Pleistocene basalts (including nephelinites with 38–40% SiO2, not shown in Fig. 2; Weinstein, 2000Go) but show no systematic increase of alkalis with decreasing SiO2 (Fig. 2).


Figure 2
View larger version (23K):
[in this window]
[in a new window]
 
Fig. 2. TAS (total alkalis vs SiO2; Le Bas et al., 1986Go) diagram showing the classification of the Golan and Galilee Plio-Pleistocene basalts. Also shown are fields for Middle Miocene basalts from the Yizre'el Valley and southeastern Galilee (Weinstein, 2000Go; excluding two nephelinites with 38 and 40% SiO2; dashed line) and Miocene to Pleistocene basalts from northeastern Jordan (Shaw et al., 2003Go; dotted line). All data are plotted on a volatile-free basis.

 
Most Golan–Galilee samples have MgO >7% (up to 11·6%) and Mg numbers from 0·55 to 0·66 (Fig. 3a). The hawaiites typically have lower MgO concentrations, implying a relatively higher degree of differentiation. Alumina (12–15% Al2O3), and to a lesser extent Na2O (2·6–5·1%), correlate negatively with MgO (Fig. 3b and c). All other major elements (e.g. CaO, Fig. 3d) do not show any correlation with MgO (either general or within-suite). Nickel (3–326 ppm), Cr (5–423 ppm) and, to a lesser extent, Co (22–68 ppm) correlate positively with MgO concentration (Fig. 3g and h), reflecting olivine and clinopyroxene fractionation. On the other hand, incompatible elements, as well as TiO2 and P2O5, do not correlate with MgO (Figs 3e, f and 4a–d), implying that the variability in their concentrations is not controlled by fractional crystallization. However, they correlate positively with each other and negatively with SiO2 (Fig. 4e–h).


Figure 3
View larger version (21K):
[in this window]
[in a new window]
 
Fig. 3. (a) Mg numbers, (b)–(f) major elements and (g, h) compatible trace elements vs MgO (wt %) in the Golan–Galilee basalts. Symbols are the same as in Fig. 2. Lines show simulations of batch differentiation of (1) olivine, (2) clinopyroxene and (3) a mixture of 60% olivine and 40% clinopyroxene.

 

Figure 4
View larger version (19K):
[in this window]
[in a new window]
 
Fig. 4. Concentration of incompatible trace elements vs wt % MgO and wt % SiO2 in the Golan–Galilee basalts. Symbols are as in Fig. 2.

 
Primitive mantle normalized incompatible trace element patterns (Fig. 5a–c) are similar to those of ocean island basalts (OIB, Fig. 5c), with strong enrichment in most highly incompatible elements. This enrichment is stronger in the basanites than in the alkali basalts. The REE show a continuous decrease of normalized concentrations from La to Yb (Fig. 6; except for a slight positive Eu anomaly in some of the samples), with (La/Yb)CN ratios of 5–20 in the alkali basalts and 18–40 in the basanites (CN indicates chondrite normalized). Concentrations of heavy REE (HREE) in all Golan–Galilee basalts are lower than in typical mid-ocean ridge basalt (MORB) and OIB. Slopes of middle REE (MREE) to HREE are similar to those of OIB in the alkali basalts (Tb/Yb(CN) ratios of 1·6–3·4), but higher in the basanites (2·8–4·2). All Golan–Galilee basalts show prominent anomalies in primitive mantle normalized trace element diagrams (Fig. 5). These include positive anomalies in Ba, Nb, Sr, P and Zr, and negative anomalies in Rb. The positive P anomalies and negative Rb anomalies are more prominent in the basanites than in the alkali basalts (Fig. 5a and b), and the Ba and Nb anomalies do not always exist. In addition, the basanites show negative anomalies of K and Pb (Fig. 5a and b). In the alkali basalts and in few basanites, there are also negative Th and U anomalies. Basalts from northeastern Jordan show similar Nb, Sr, P, Rb and K anomalies, but they do not show positive Ba and P anomalies (Shaw et al., 2003Go). In general, the ratios of the REE to their neighboring HFSE in Fig. 5 (e.g. La/Nb, Sm/Zr and Eu/Ti) increase from the alkali basalts to the basanites (Fig. 7a and b). Similarly, P/Zr ratios also increase from alkali basalts to the basanites (Fig. 7c). Basalts from northeastern Jordan show similar trends of increasing La/Nb from alkali basalts to basanites (Fig. 7a; Shaw et al., 2003Go), although their overall La (and Nb) concentrations are lower than in the Golan–Galilee. On the other hand, in volcanic rocks from the Dead Sea area (Shaw, 2003Go) La/Nb ratios decrease from alkali basalts to basanites (Fig. 7a).


Figure 5
View larger version (27K):
[in this window]
[in a new window]
 
Fig. 5. Primitive mantle normalized trace element diagram showing the composition of representative basalts from (a) Golan, (b) Galilee and (c) averages and end-member compositions (Wy-16 and Wy-337; see text). Symbols are as in Fig. 2, but Pleistocene samples are subdivided into Early ({diamond}) and Late ({diamondsuit}). Dashed lines are for alkali basalts and continuous lines for basanites. Primitive Mantle values are from McDonough & Sun [1995Go; for Pb, we use the value of 0·071 ppm, recommended by Sun & McDonough (1989)Go]. The OIB composition is from Sun & McDonough (1989)Go.

 

Figure 6
View larger version (18K):
[in this window]
[in a new window]
 
Fig. 6. Chondrite-normalized REE patterns for Golan–Galilee alkali basalts and basanites. Symbols and lines are as in Fig. 5. Chondrite concentrations are from McDonough & Sun (1995)Go.

 

Figure 7
View larger version (21K):
[in this window]
[in a new window]
 
Fig. 7. (a), (b) REE/HFSE vs REE concentration; (c) P/Zr vs wt % P2O5; (d) and (e) ratios of highly incompatible to less incompatible elements vs concentration of the highly incompatible element (H/I vs H; see text). Symbols are as in Fig. 2. The field of alkali basalts is encircled by a dashed line and that of basanites by a continuous line. Fields for Miocene basalts from the Yizre'el Valley and southeastern Galilee (Weinstein, 2000Go) and basalts from northeastern Jordan and the Dead Sea region (Shaw, 2003Go) are shown in (a) and (d) (a few outliers were excluded in the construction of these fields).

 
A central observation in this study is the patterns defined by the Golan–Galilee basalts in plots of some incompatible element ratios vs concentration (e.g. La/Nb vs La, Nb/Zr vs Nb and La/Sm vs La, Fig. 7a–e). In these, the Pliocene alkali basalts define steeper trends than the Pleistocene Golan basanites. The general trends defined by the Golan samples in these diagrams are convex-upward, which we interpret below as mixing curves.

Spatial and temporal variability
In both the Golan and Galilee, basalt chemistry appears to change with both time and geographic location, from alkali basalts in the south to mainly basanites in the north. In the Golan, Early Pliocene lavas from the south are mainly alkali basalts, with a few hawaiites (Fig. 2). In the western Golan, the Late Pliocene basalts are also alkali basalts, but with lower silica contents (almost all samples contain normative nepheline compared with only 60% of the Early Pliocene samples). Pleistocene lavas from the northern and eastern Golan are mainly basanites and hawaiites (Fig. 2). In the Galilee, Early Pliocene lavas from the south (Lower Galilee) are mainly alkali basalts and hawaiites, whereas Late Pliocene lavas from the north (Upper Galilee) are all basanitic (Fig. 2).

The spatial–temporal geochemical variability is also reflected in the concentrations of incompatible elements and in H/I ratios (where H indicates a hygromagmatophile, highly incompatible element, and I indicates a moderate incompatible element, e.g. Nb/Zr and La/Sm). Basanites (Golan Pleistocene and Galilee Late Pliocene) have higher incompatible element concentrations and ratios than the alkali basalts (Early Pliocene and most of the Golan Late Pliocene samples; see Figs 5a–c and 7d, e). Also, most of the basanites contain >2·5% TiO2 and >0·9% P2O5 (averages 2·8% and 1·2%, respectively) whereas the alkali basalts have lower concentrations (averages of 2·3% and 0·5%, respectively; Table 1).

Ratios of light REE (LREE) and MREE to HFSE of similar compatibility (e.g. La/Nb, Sm/Zr and Eu/Ti) also increase in the Golan, from the southern Early Pliocene alkali basalts to the northern and eastern Pleistocene basanites (Table 2, Fig. 7a and b). The behavior of P is similar to that of the REE; P/Zr increases from 7–15 in the southern Golan to 15–30 in the northern Golan (Fig. 7c). There is no systematic variability (between suites or within an individual suite) in MgO and compatible element concentrations, except for the Pleistocene basanites from the northern Golan, for which Weinstein et al. (1994)Go observed that Late Pleistocene basalts have higher MgO and Ni concentrations than Early Pleistocene basalts.


View this table:
[in this window]
[in a new window]
 
Table 2: Key element ratios in samples Wy-337 and Wy-16 and averages in the various units

 
An interesting spatial–temporal pattern is shown in the trace element and P concentrations of basalts from the central Golan (Fig. 8). More than half of the samples were taken from flows along a 150 m high cliff section in the Zavitan river in the Yahudiya plateau (western Golan, Fig. 1), spanning a period from the Early to Late Pliocene (3·3–1·8 Ma, Table 1 and Fig. 8). The rest of the samples were collected along a line drawn from the Zavitan river to the Bashanit area in the eastern Golan (A–A' in Fig. 1). The samples along this line are all Pleistocene in age (1·3–0·4 Ma, Table 1 and Fig. 8), but there is no clear trend of decreasing age toward the NE (Fig. 8), even though the easternmost sample is also the youngest. In general, the concentration of P and incompatible elements in the Pleistocene samples is higher than in the Pliocene ones. However, whereas during the Late Pliocene there is a clear temporal trend of increasing concentrations with time, in the Pleistocene samples the main change is with location (towards the NE) and not with time (e.g. concentrations of all three elements increase from the 0·7 Ma sample to the 1·3 Ma one), even though the youngest sample also has the highest concentrations. We note that the location of the vents that fed the Late Pliocene basalts of the Yahudiya plateau is unknown. They could either be further to the east (now covered by Pleistocene basalts, Fig. 1) or unexposed local fissures in the western Golan. Thus, the control on the variation in basalt composition in the central Golan is both temporal and spatial.


Figure 8
View larger version (12K):
[in this window]
[in a new window]
 
Fig. 8. Concentrations of (a) Nb, La and (b) P2O5 in a Plio-Pleistocene section from the central Golan. The section is a combination of (1) a Late Pliocene section (3·3–1·7 Ma, Table 1) taken along the canyon wall of the Zavitan River in the Yahudiya plateau and (2) a Pleistocene section of surface basalts along a line (A–A', Fig. 1) from the Yahudiya plateau to the NE. The numbers inside the figure are K–Ar ages of some of the samples (in Ma).

 
Unlike the Golan–Galilee, Neogene–Quaternary basalts from northeastern Jordan do not show any temporal or spatial patterns (Shaw et al., 2003Go).

Sr and Nd isotope geochemistry
The Sr–Nd isotopic compositions of the Golan–Galilee alkali basalts and basanites are given in Table 1 and plotted in Fig. 9a and b. Nd isotope ratios are remarkably uniform in the range {varepsilon}Nd = 4·0–5·2, whereas 87Sr/86Sr varies between 0·7031 and 0·7034, a composition close to that of PREMA (PREvalent MAntle, Wörner et al., 1986Go). A similar range of Sr–Nd isotope compositions has been reported for Miocene basalts from the Yizre'el Valley, southeastern Galilee and the southern Golan ({varepsilon}Nd 4·1–5·4 and 87Sr/86Sr 0·7031–0·7034, Weinstein, 2000Go; Stein & Hofmann, 1992Go). The isotopic compositions of the Golan–Galilee basalts fall within the field of basalts from Saudi Arabia (Altherr et al., 1990Go, Fig. 9a), although the latter show larger ranges in both 87Sr/86Sr and {varepsilon}Nd (0·7028–0·7035 and 3·6–7·6, respectively), extending towards the Red Sea MORB field (e.g. Volker et al., 1993Go). Basalts from northeastern Jordan (Shaw et al., 2003Go) show a significantly wider range of Sr isotope compositions (0·70305–0·70415), whereas their Nd isotope ratios ({varepsilon}Nd = 4·1–6·4) are relatively similar to those of the Golan–Galilee basalts. Interestingly, Sr and Nd isotopes are strongly negatively correlated in the Jordanian basalts (Shaw et al., 2003Go), probably as a result of contamination with an upper crust component. This trend is not observed in the Golan–Galilee and is just slightly shown by the Saudi Arabian basalts (Altherr et al., 1990Go).


Figure 9
View larger version (27K):
[in this window]
[in a new window]
 
Fig. 9. Sr and Nd isotope ratios in the Golan basalts (including samples from Stein & Hofmann, 1992Go), Miocene samples from Weinstein (2000)Go ({circ}, Yizre'el Valley and southeastern Galilee; •, one southern Golan sample) and Mezozoic samples ({diamond}, Stein & Hofmann, 1992Go). Other symbols are as in Fig. 2. Also shown are: (a) fields of compositions of Red Sea MORB (Hegner & Pallister, 1989Go; Altherr et al., 1990Go; Volker et al., 1993Go; F. Volker, unpublished data, 1994); Saudi Arabian basalts (Altherr et al., 1990Go); northeastern Jordan (Harrat ash Shamm) basalts (Shaw et al., 2003Go); Yemen Quaternary basalts (Baker et al., 1997Go); lithospheric mantle xenoliths from Israel and Arabia (Brueckner et al., 1988Go; Henjes-Kunst et al., 1990Go; Stein et al., 1993Go) and the compositions of PREMA (Wörner et al., 1986Go) and the Afar Plume (Vidal et al., 1991Go; Schilling et al., 1992Go; Deniel et al., 1994Go); (b) fields of OIB, Arabian lower and upper crust (Duyverman et al., 1982Go; Cohen et al., 1984Go; Hegner & Pallister, 1989Go; Henjes-Kunst et al., 1990Go; McGuire & Stern, 1993Go; Stein et al., 1993Go; O. Gazit, unpublished data, 2005) and East African lower crust (Cohen et al., 1984Go). BE, Bulk Earth.

 
The Sr isotope ratios of the Golan–Galilee basalts are negatively correlated with incompatible element concentrations (e.g. Nb and Sr, Fig. 10a and b) and positively correlated with Rb/Sr ratios (Fig. 10c) and 1/Sr (Fig. 10d).


Figure 10
View larger version (15K):
[in this window]
[in a new window]
 
Fig. 10. 87Sr/86Sr ratios vs (a) Nb, (b) Sr, (c) 87Rb/86Sr and (d) 1/Sr in the Golan–Galilee basalts (symbols as in Fig. 2), including data from Stein & Hofmann (1992)Go, and in Miocene basalts (Weinstein, 2000Go; symbols as in Fig. 9) and two Late Pliocene to Pleistocene samples from the Korazim area north of the Sea of Galilee ({blacktriangleup}, Y. Weinstein, unpublished data, 1998). Also shown in (c) are the compositions of calculated APR and PYV sources ({blacksquare}; Rb/Sr ratios are from Tables 4 and 5, and the 87Sr/86Sr values are those of Wy-337 and Wy-52, respectively; see text).

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
In the following discussion we focus on the nature and the causes of the compositional variability observed in the Golan–Galilee basalts. Stein & Goldstein (1996)Go and Stein et al. (1997)Go have proposed that all Phanerozoic basalts from Arabia have a lithospheric rather than an asthenospheric mantle source. Following this approach, we show below that the compositional variability of the Golan–Galilee basalts reflects the heterogeneous nature of the underlying lithospheric mantle.

Alteration and crustal contamination processes
Alteration probably did not have a significant effect on the chemistry of the Golan–Galilee basalts. Almost all the samples are fresh, showing no petrographic evidence for alteration except for the presence of iddingsite that partially or completely replaces olivine. However, the scatter shown by the Galilee basanites in plots of Rb (and to a lesser extent K) versus SiO2 (Fig. 4h) or versus other incompatible elements (not shown) implies that this group was subjected to some secondary process that affected the concentration of the more mobile elements. As the Golan basalts show good correlations in the above plots (Fig. 4h) and do not show any other signs of alteration, we chose to focus on this group in the petrogenetic discussion below.

The degree of crustal contamination of the Golan–Galilee basalts, as well as of other Phanerozoic basalts from Israel (Stein & Hofmann, 1992Go), is probably fairly small on the basis of their Sr and Nd isotope compositions (Fig. 9b). Upper crustal contamination should result in a significantly larger range of variability and higher Sr isotope ratios. A contribution from the lower crust cannot be absolutely ruled out, but should probably result in some variability in Nd isotope composition (Fig. 9b), unlike the uniform values observed in the Golan–Galilee basalts (Fig. 8a). Similarly, Altherr et al. (1990)Go showed that most Saudi Arabian basalts carry near pristine mantle isotopic signatures and ruled out crustal contamination. We note, however, that some Arabian basalts, e.g. from Yemen and especially those from nearby northeastern Jordan, do show a crustal Nd–Sr isotopic signature (Baker et al., 1997Go; Shaw et al., 2003Go; Fig. 9a).

Differentiation of the primary magmas
The Golan–Galilee basalts are clearly not primary in composition but are, nevertheless, relatively primitive. Their Mg numbers are between 0·66 and 0·37 (Table 1), MgO contents are usually <12 wt %, and Ni and Cr are <350 ppm and <500 ppm, respectively. This indicates that all the basalts have experienced a certain degree of crystal fractionation. However, most of the incompatible elements show no correlation with MgO concentration (e.g. Fig. 4a–d). This indicates that the wide concentration range of incompatible elements (e.g. a factor of five in Nb and nine in La concentrations) in the Golan–Galilee basalts was not generated by differentiation processes. Rather, the variability in incompatible element concentrations could have originated from mantle source heterogeneity or from variable degree partial melting. The negative correlation of incompatible elements (e.g. Nb, La, Ba, Rb) with SiO2 (Fig. 4e–h) also supports this conclusion.

Because all of the Golan–Galilee magmas underwent differentiation, their original (primary) incompatible element composition must have been modified to some extent. Based on the good correlation of Ni and MgO (Fig. 3g), we assume that most primitive magmas followed more or less similar differentiation paths. Therefore, to better discuss melting regimes and mantle source compositions, we have corrected the Golan–Galilee basalt compositions to a reference value of 10% MgO by adding various minerals to the bulk-rock composition, following the approach of Klein & Langmuir (1987)Go. This correction procedure should restore the magma composition close enough to that of the primary magmas (e.g. Takahashi & Kushiro, 1983Go; Falloon & Green, 1988Go). About 10% of the samples have MgO >10% (Table 1, Fig. 3) and required no correction.

The dominance of olivine and clinopyroxene in the phenocryst population and the good correlation of Ni, Cr and Mg numbers with MgO concentrations (Fig. 3a, g and h) suggest that these minerals were the main fractionating phases. Plagioclase fractionation may be ruled out based on the negative correlation between Al2O3 and MgO (Fig. 3b), the high Sr concentrations, the rarity of Eu negative anomalies (Fig. 6) and the minor presence of plagioclase phenocrysts in the basalts. Amphibole is also ruled out because of the absence of amphibole phenocrysts. As indicated by the fractionation simulations shown in Fig. 3a, b, d and e, clinopyroxene cannot be the main fractionating phase. On the other hand, olivine alone cannot account for the trend observed in the CaO vs MgO diagram (Fig. 3d) and probably also not for that in the Al2O3 vs MgO diagram (Fig. 3b, too steep). Altogether, the data trends are best explained by a fractionating mineral assemblage with an olivine:clinopyroxene ratio of about 60:40, and we use these proportions to correct to 10 wt % MgO. In a few samples, where large deviations from the 60:40 curve were noticed, the proportions were changed accordingly. As the distribution coefficients of both clinopyroxene and olivine are <1 for all incompatible elements, the effect of the correction procedure on the incompatible elements is mainly one of dilution, by up to 17–18% for the highly incompatible elements (e.g. La) and less for others.

The 60:40 ol:cpx ratio assumed in the fractionating assemblage is very different from the usually observed strong dominance of olivine in the phenocryst population of basalts (Table 1). The phenocryst assemblage probably reflects a later, shallower, crystallization stage in which olivine is expected to dominate (e.g. Takahashi & Kushiro, 1983Go; Sack et al., 1987Go; Weinstein et al., 1994Go).

Melting of a heterogeneous lithospheric mantle
On the basis of the above discussion, we consider that the compositional spectrum observed in the Golan–Galilee basalts, especially that of incompatible elements, was not produced by crystal fractionation processes, but reflects mantle source heterogeneity and/or partial melting processes.

Considering the young age of the Golan–Galilee basalts and their low Rb/Sr ratios, the variability of Sr isotope ratios (Fig. 9a) precludes the derivation of the parent magmas from a homogeneous mantle source. Furthermore, the negative correlation of Sr isotope ratios with incompatible element concentrations (e.g. Nb and Sr, Fig. 10a and b) and the positive correlation with 87Rb/86Sr and with 1/Sr (Fig. 10c and d) suggest that the spectrum of melt compositions was produced by binary mixing. Mixing is further supported by the upward-convex patterns in some of the H/I vs H diagrams (e.g. Fig. 11a–d) that differ from the linear trends expected in a case of variable degrees of melting (see Hofmann & Feigenson, 1983Go). The linear trends observed in diagrams of MREE/HREE vs LREE/HREE (e.g. Tb/Yb vs La/Yb, Fig. 11f) also differ from those expected in variable degree of partial melting.


Figure 11
View larger version (21K):
[in this window]
[in a new window]
 
Fig. 11. (a)–(e) Incompatible element ratios vs concentration in the Golan–Galilee basalts. All data are corrected to 10% MgO (see text). Symbols are as in Fig. 2. The fields of alkali basalts and basanites are as in Fig. 7. Also shown are mixing lines (continuous; not numbered) and simulations (dashed curves) of non-modal batch partial melting (Shaw, 1970Go) in the range 0·1–25% for five source mineralogies: 1, spinel lherzolite (ol:opx:cpx:sp ratios of 60:20:15:5 and –15:50:50:15 in the rock and in the melt, respectively); 2, garnet lherzolite (ol:opx:cpx:gt ratios of 60:20:10:10 and 0:30:35:35 for the rock and melt, respectively); 3, amphibole lherzolite (ol:opx:cpx:am ratios of 55:20:5:20 and –15:35:15:65, respectively); 4, phlogopite–garnet lherzolite (ol:opx:cpx:gt:phl 65:15:10:7:3 and 5:15:60:10:10, respectively); 5, amphibole–garnet clinopyroxenite (cpx:sp:gt:am 75:10:10:5 and 25:20:30:25, respectively). The source composition is arbitrarily assumed to be that of the Primitive Mantle (‘PM’, McDonough & Sun, 1995Go), except for Nb, which is assumed to be 2 x PM. Partition coefficients are listed in Table 3. (f) Th/Yb (MREE/HREE) vs La/Yb (LREE/HREE) in the basalts, including mixing lines and batch melting simulations with mineralogies similar to those in (a)–(e) and arbitrarily assuming source La and Tb are 2 x PM, while Yb is 1 x PM. (See text for further explanation.)

 
The incompatible element characteristics of the Golan–Galilee basalts (Figs 7 and 11a–e) could alternatively be considered as a combination of two linear arrays, each of which originates from variable degree melting of a distinct mantle source. However, in this scenario we should expect each group to cluster around the 87Sr/86Sr ratio of its mantle source. The continuous variation of Sr isotope ratios and their correlation with incompatible element concentrations (Fig. 10a and b) suggests that this is not the case. In addition, it is hard to fit the trace element data with a partial melting model involving two distinct mineralogical assemblages. Source composition and mineralogy that fit a given set of element ratios and concentrations cannot produce the trends of other sets of elements. For example, in the diagram of Nb/Zr vs Nb (Fig. 11a), the relatively moderate trend defined by the more enriched samples (mainly the basanites) is similar (parallel) to the slope of a spinel peridotite melting curve (curve 1 in Fig. 11), whereas the steep slope of the less enriched samples (<40 ppm Nb) is similar to that of a pyroxenite curve (curve 5). On the other hand, in the diagram of La/Nb vs La (Fig. 11b), the trend of the enriched samples is similar to that of an amphibole lherzolite (curve 3) or of a phlogopite peridotite melting curve (curve 4), whereas the steeper slope of the alkali basalts could not be produced by any source mineralogy. The opposite is true for La/Sm vs La (Fig. 11c), where the relatively flat slope defined by the enriched samples does not fit any of the melting curves. Moreover, the positive correlation shown by both basanites and alkali basalts in Sr/Nb vs Sr (Fig. 11e) does not match any of the melting simulations. Therefore, we consider that the incompatible element characteristics of the Golan–Galilee basalts are mainly controlled by binary mixing of partial melts from discrete mantle sources, especially in the case of the Golan basalts. Mixing of discrete mantle melts, as opposed to mixing of mantle components (e.g. variable degree of source enrichment), is our preferred model because we anticipate that in the latter the correlation of 87Sr/86Sr with incompatible element concentrations (Fig. 10a and b) would be erased during partial melting processes. Nevertheless, the rarity of reversely zoned phenocrysts in the Golan and Galilee basalts (Weinstein et al., 1994Go) implies that mixing occurred early in the ‘life’ of the magma, close to the source area, before the phenocrysts started crystallizing. Mixing as the main control on basalt chemistry was also identified by Shaw et al. (2003)Go in the neighboring northeastern Jordanian basalts.

We therefore conclude that the compositions of the Golan–Galilee basalts reflect mixing of melts derived from two distinct mantle sources. One end-member is the main source for the alkali basalts, and is characterized by relatively low concentrations of incompatible elements, low H/I ratios (e.g. Nb/Zr) and low REE/HFSE ratios (e.g. La/Nb). The second end-member is enriched in incompatible elements and in the above ratios, and is the main source for the basanites. The compositional range observed in the basanites (e.g. Fig. 11a–c and f) suggests that the enriched source is heterogeneous; namely, it includes a spectrum of enriched compositions.

It is interesting to note that alkali basalts and basanites from the Dead Sea area of Jordan (Shaw, 2003Go) show a very clear upward-convex trend on a diagram of Nb/Zr vs Nb (Fig. 7d). Miocene basalts from the Yizre'el Valley and the southeastern Lower Galilee (Weinstein, 2000Go) also show an upward-convex pattern on the same diagram (Fig. 7d). On the other hand, basalts from northeastern Jordan, which were also interpreted in terms of binary mixing by Shaw et al. (2003)Go, do show a positive correlation on this diagram but not an upward-convex pattern (Fig. 7d). In all three datasets, an upward-convex pattern is not identified on the La/Nb vs La diagram (Fig. 7a), probably because of the relatively small range of La concentrations in these datasets. However, whereas the Israeli Miocene and the northeastern Jordanian basalts do show a positive correlation (similar to the Golan–Galilee basalts), the Dead Sea samples show a clear negative correlation of La/Nb ratios with La concentrations (Fig. 7a). If the Nb/Zr vs Nb trends (Fig. 7d) are taken as evidence for binary mixing in petrogenesis of the Dead Sea basalts, then the enriched end-member in this case must have had a lower La/Nb ratio than the alkali basaltic end-member.

Unlike Sr isotopes, Nd isotope compositions are uniform in all the Golan–Galilee basalts ({varepsilon}Nd = 4·6 ± 0·6), suggesting a genetic relationship between the two mantle sources. Based on this observation, we suggest that the two mantle sources may have been derived from a common reservoir, which separated into the above two sources. The following discussion focuses on the characterization of the two sources and on their evolution.

Source of the alkali basalts
The melt derived from the alkali basaltic source can be closely characterized by sample Wy-337 from the southern Golan, as it has the lowest concentrations of highly incompatible elements. On a primitive mantle normalized trace element diagram (Fig. 5c), Wy-337 displays prominent peaks at Ba, Nb, Sr, Ti and Zr, and negative anomalies at Rb and U. Because there is no major mantle mineral that can significantly fractionate Rb from Ba, we attribute the low Rb/Ba (Fig. 5a, Table 2) to preferential removal of Rb by hydrous fluids during the petrogenesis of the basalts or the evolution of their mantle sources. Equilibration of the source with fluids may also cause the negative U anomaly (high D(fluid–rock), Keppler, 1996Go) and the positive Nb anomaly, which as a HFSE has a tendency to be strongly retained in the rock (Keppler, 1996Go; see Table 3) compared with other elements (including the REE).


View this table:
[in this window]
[in a new window]
 
Table 3: Mineral–melt and fluid–melt partition coefficients used in this study

 
The relatively high Sr/Ba and the higher than PM Nb/Ba of Wy-337 (0·71 and 1·10, respectively, PM-normalized ratios) suggest that Ba may have been buffered by a residual hydrous phase. This phase was probably an amphibole rather than a phlogopite, as buffering by phlogopite would cause extremely low Ba concentrations in the partial melts (DBa(phl-melt) >1; Adam et al., 1993Go; LaTourrette et al., 1995Go), and the Nb/Ba ratio would, therefore, be much higher than that of the primitive mantle. The presence of phlogopite is also not favored by the relatively low K/Nb ratio (<1, PM ratios) in the alkali basalts (Fig. 5a–c).

It is, therefore, suggested that the source for the alkali basaltic end-member melt was an amphibole-bearing peridotite (designated as APR) that underwent partial dehydration, which caused the preferred removal of Rb and U. The presence of amphibole limits the depth of the source to <90 km, as 30 kbar is the maximum pressure for amphibole stability (Green, 1973Go; Brey et al., 1983Go; Wallace & Green, 1988Go). The low Sm/Nd (0·73 PM, Table 2) and the relatively high Tb/Yb (1·7 PM, Table 2 and Fig. 5c) of Wy-337 suggest buffering by garnet, which implies that the APR source was located deeper than 60 km (~20 kbar, Mysen & Boettcher, 1975Go). Combining the two constraints, the APR source appears to have been located at 60–90 km depth.

A source depth shallower than 90 km implies a lithospheric rather than a deeper mantle origin, as the low heat flow beneath Israel (~40 mW/m2, Eckstein, 1979Go) predicts a lithosphere thickness in excess of 100 km. Unlike McGuire (1988)Go, Nasir (1992)Go and Shaw et al. (2003)Go, we consider that the elevated temperatures calculated for peridotitic xenoliths (e.g. 1015°C at 15 kbar in western Arabia, McGuire, 1988Go; McGuire & Bohannon, 1989Go; Nasir, 1992Go; Stein et al., 1993Go) hosted in Arabian Cenozoic basalts represent local heating by magmas (although not the host magmas) and not regional elevated geotherms (see Stein et al., 1993Go).

A lithospheric mantle source for the Israeli and Arabian basalts was previously suggested by Stein & Goldstein (1996)Go, who showed that all juvenile Late Proterozoic magmas from the Arabian–Nubian Shield (ANS) and Phanerozoic alkali basalts from the Arabian plate lie on the same trends in Nd, Sr and Pb isotope evolution diagrams. These trends differ significantly from the growth curves of asthenospheric MORB-source mantle. Moreover, isotope ratios of Israeli Mezozoic basalts (Triassic to Cretaceous) are very similar to those of the Plio-Pleistocene Golan–Galilee ones (corrected initial 87Sr/86Sr ratios of 0·7029–0·7031 and {varepsilon}Nd of 4·5–5·4, Fig. 9a), probably implying a closer genetic relationship between these two groups than between the coeval Golan–Galilee and Arabian basalts. Stein & Goldstein (1996)Go concluded that the sources of these magmas have remained associated with the ANS since its formation during the Late Proterozoic Pan-African orogeny and that these sources are, therefore, located in the lithospheric mantle.

Following Stein et al. (1997)Go, it is suggested that the incompatible element characteristics of the APR source were acquired during Late Proterozoic Pan-African events as a consequence of fluxing of the newly formed mantle wedge by slab-derived fluids. In the slab, the fluids equilibrated with eclogite and peridotite. In the wedge, a hydration front propagated upwards, leaving behind amphibole peridotite. Trace elements carried by the fluids partitioned between the fluid and the hydrated peridotite. This chromatographic effect (Navon & Stolper, 1987Go) led to the stripping of Nb and other trace elements with high D(amphibole peridotite–fluid) from the migrating fluid and their preferential retention at deeper wedge levels. Lanthanum and other trace elements with lower D(peridotite–fluid) remained in the fluid until shallower wedge levels. As a result, Nb-enriched zones were produced in deeper parts of the wedge, whereas Nb-depleted zones were generated at shallower levels or away from the fluid migration channels. The low La/Nb and Th/Nb ratios of the alkali basaltic melts imply derivation from a Nb-enriched zone.

In Table 4 and Fig. 12 we present a possible model for the petrogenesis of the alkali basaltic Golan end-member (Wy-337). It includes: (1) equilibration of hydrous fluids with a subducting slab (an arbitrary eclogite/peridotite ratio of 1:1 and a fluid/rock ratio of 0·005); (2) infiltration of these fluids and hydration of the mantle wedge, producing amphibole peridotites with Nb-enriched and Nb-depleted zones; (3) partial dehydration of the amphibole peridotite (fluid/rock ratio of 0·005), which significantly reduced, although it did not eliminate, the content of amphibole in the source rock and which caused the large decrease in the Rb content of the source rock; finally, (4) partial melting (~13%), which produced the alkali basaltic end-member melt. Further details of the model can be found in the notes to Table 4.


Figure 12
View larger version (11K):
[in this window]
[in a new window]
 
Fig. 12. Primitive Mantle-normalized trace element diagram of the simulated alkali basalt end-member melt, compared with the composition of the alkali basalt sample Wy-337 from the southern Golan ({square}). Also shown is the composition of the amphibole peridotite source after partial dehydration (APR). (See Table 4 and text for details.)

 

View this table:
[in this window]
[in a new window]
 
Table 4: Forward batch melting simulation of production of the alkali basaltic end-member melt (Wy-337)

 
In summary, the alkali basaltic end-member magmas originated by partial melting of a gt–amph-peridotite (APR). This source composition was formed by Late Proterozoic interaction of peridotite with slab-derived fluids, followed by partial dehydration.

Source of the basanites
For the characterization of the partial melt of the enriched end-member, we used the Golan basanite Wy-16, which is one of the samples with the highest concentrations of highly incompatible elements (Table 1) and the steepest trace element pattern (Fig. 5a and c). However, as mentioned above, the enriched end-member does not have a uniform composition. Therefore, sample Wy-16 does not necessarily represent the exact mixing end-member for all the Golan magmas.

Compared with the alkali basalts, the enriched end-member (Wy-16) has (Table 2):

(1) a stronger negative Rb anomaly (Rb/Ba = 0·35, Rb/Sr = 0·47, PM ratios), which following the above reasoning, implies an intense leaching by fluids;

(2) a negative K anomaly;

(3) a steep LREE–MREE slope (Sm/Nd = 0·5), which suggests either a dominance of clinopyroxene or a major role for garnet in the source;

(4) a very steep MREE–HREE slope (Tb/Yb = 4·2), which implies either a high garnet content in the source or a very low degree of melting;

(5) REE/HFSE ratios higher than those of the alkali basalts (La/Nb approaching unity, compared with 0·5 in Wy-337), which may point to the involvement of a Nb-depleted source;

(6) a positive P anomaly (high P/Nd ratio), which could reflect the dominance of clinopyroxene (phosphorus has a low D(Cpx–melt) compared with its neighboring REE elements, Liotard et al., 1988Go; see Table 3) or the involvement of apatite.

It is evident that the source for the basanitic melts has to be significantly more enriched than the APR source. In fact, an unreasonably low degree of partial melting (<0·1%) would be needed to produce the Wy-16 La concentration of 134 PM from the APR source. Also, an APR-like source could not produce the very low Rb/Ba of the basanites (Table 2) and the negative K anomaly. To achieve these values, fluid/rock ratios significantly greater than in the case of the alkali basalt source would be required. For example, complete dehydration of an APR source (30% amphibole carrying 1–2 wt % H2O) would produce a fluid/rock ratio of <0·6%. An external source for the fluids (except for the slab-derived fluids) is rejected, as western Arabia area has been far from any subduction-related plate margin since the Late Proterozoic. Late-stage Rb and K loss during degassing on eruption or secondary alteration is also rejected based on the absence of any correlation between low Rb or K and petrographic evidence for alteration, and the correlation of 87Sr/86Sr with Rb/Sr (Fig. 10c), which suggests a mantle source effect. Thus, we conclude that the enriched source was probably rich in clinopyroxene, contained garnet, was highly enriched in incompatible elements (compared with the APR source) and had interacted with a relatively large volume of fluid.

Following the above assumption about a genetic relationship between the two mantle sources, we suggest that the enriched chemical signature observed in Wy-16 and other basanites was generated through a two-stage melting process (Table 5), similar to the process suggested by Francis (1991)Go and more recently by Thompson et al. (2005)Go, as follows. The first stage included the formation of the enriched source by partial melting of an original amphibole peridotite source. During this early melting event, water and other fluids concentrated in the melt. Assuming 1–2% H2O in the amphibole, 9% melting (Table 5, row 9) while dehydrating 18% of the peridotitic amphibole (see Table 5), we obtain a melt with 2–4% H2O. This melt then segregated and completely crystallized, forming pyroxenite veins (designated PYV) within the mantle peridotite (Fig. 13). The crystallization depth could be either shallower than or similar to that of the parental peridotite. Upon crystallization, the melt lost most of its fluids (Table 5, row 11), which preferentially carried Rb and K out of the veins. More recently, during the Plio-Pleistocene, the veins were partially remelted (12%, similar to the APR source) to yield the basanitic end-member magmas (Table 5, row 15; Fig. 13). Upon melting, the magmas equilibrated with neighboring peridotites (Table 5, rows 16 and 17), acquiring close to primary Mg numbers with very little effect on their final trace element composition.


Figure 13
View larger version (12K):
[in this window]
[in a new window]
 
Fig. 13. Primitive Mantle-normalized trace element diagram of the simulated basanitic end-member melt, compared with the composition of the basanite sample Wy-16 ({diamondsuit}). Also shown are the compositions of the amphibole peridotite (no dehydration) and of the pyroxenitic vein source (PYV) after dehydration. (See Table 5 and text for details.)

 

View this table:
[in this window]
[in a new window]
 
Table 5: Forward simulation of two-stage melting to produce the basanitic end-member melt (Wy-16)

 
As a fully crystallized H2O-rich basaltic magma, the veins probably contained some aluminum-bearing phases, e.g. amphibole, spinel and garnet. This is supported by the presence of amphibole pyroxenite xenoliths with variable amounts of amphibole in tuffs from the Golan and spinel- and amphibole-bearing garnet clinopyroxenite xenoliths elsewhere in Israel (e.g. KM-140 with 75% clinopyroxene, 10% garnet, 10% spinel and 5% amphibole; Stein, 1987Go). Nevertheless, the amphibole content of the pyroxenite source has to be limited to a few per cent, as a larger fraction would strongly buffer Rb and Ba during dehydration, which would result in higher concentrations of these elements in the final melt. A close fit to Wy-16 composition is obtained using a non-dehydrated amphibole peridotite in the first stage (Table 5, row b) and a pyroxenite in the second stage, composed of 74% clinopyroxene, 10% orthopyroxene, 10% spinel, 5% garnet and 1% amphibole (Table 5, row d).

The calculated vein composition is similar to that of amphibole clinopyroxenite xenoliths found in Pleistocene tuffs from Birket Ram (Golan Heights, Stein, 1987Go; Stein et al., 1993Go), which were described by Mittlefehldt (1984)Go as conduit-filling melts. These xenoliths have Rb/Sr PM-normalized ratios of 0·10–0·17, in close agreement with the ratio of 0·23 for the calculated veins (although the Rb and Sr concentrations in some of these xenoliths are higher than the calculated values by up to a factor of two), and have an 87Sr/86Sr isotope ratio of 0·7031, similar to those of the basanites. The Nd isotope ratios of the Birket Ram xenoliths are also very similar to the isotope ratios of the basanites ({varepsilon}Nd = 5·1 ± 0·6 and 4·4 ± 0·4, respectively, Table 1). Pyroxenites from other localities in Israel (Mt. Carmel and Makhtesh Ramon; Stein, 1987Go) and northeastern Jordan (J. Baker, personal communication, 2005) are less similar to our calculated vein composition and are usually considered to be cumulates from basaltic magmas (Bonen, 1980Go; Stein, 1987Go).

The nature of pyroxenite melt compositions has been discussed in a number of studies since the 1960 s (e.g. O'Hara, 1968Go). Kogiso et al. (2004)Go, in a recent compilation, concluded that at pressures higher than 2 GPa, the compositions of partial melts of pyroxenites are controlled by the existence of a garnet–pyroxene thermal divide. Silica-excess pyroxenites produce silica-rich melts, whereas melts derived from silica-deficient pyroxenites are silica poor. Kogiso et al. (2004)Go also showed that with increasing degrees of melting, partial melts do not cross the divide. We assume that the composition of the early amphibole peridotitic melt (first stage) was alkali basaltic. Therefore, the resulting pyroxenite veins (PYV) were probably silica deficient, and the second generation melts derived from these veins would, therefore, be silica poor, similar to the basanitic melts. Moreover, according to Kogiso et al. (2004)Go, partial melts of silica-deficient pyroxenites tend to have lower silica contents than peridotitic melts, in agreement with the relatively low silica content in the Golan–Galilee basanites.

The two-stage melting scenario explains well the high concentrations of REE in the basanites, while avoiding the need for very low degrees of partial melting, and actually produces both end-member melts by similar degrees of melting (12–13%). The lower U concentration predicted by our model (Fig. 13) is probably the result of our use of the extremely high D(fluid–melt) values of Keppler (1996)Go, in which the fluid was a 5 M HCl solution. Values of D(fluid–melt) in Cl-free solutions are usually significantly lower, especially for U (0·36 and 0·003, respectively; Keppler, 1996Go). Using an intermediate value, between those for 5 M HCl and Cl-free solutions, would produce a concentration similar to those in the basanites. Nb concentration in the calculated melt is overestimated (La/Nb in Wy-16 equals 1·0 compared with the calculated value of 0·5; Table 5 and Fig. 13). The lower Nb in the basanites may point to a contribution from a Nb-depleted zone (Stein et al., 1997Go) during the early melting stage. Contribution from such a depleted zone is probably also the reason for the relatively low concentrations of Zr and TiO2 in the basanites (Fig. 13). It should be remembered that the enriched end-member is not homogeneous, and some samples (e.g. Wy-292) define a trend toward an enriched end-member with lower La/Nb, which reflects a lower contribution from the Nb-depleted source (Fig. 11b). It is interesting to note that low Nb basalts, which probably originated from Nb-depleted zone, have also been identified in the Dead Sea area (Shaw, 2003Go; Fig. 7a). An alternative approach to explaining the high La/Nb in the basanites (and the high P concentrations in the basanitic magmas; Table 1) may be the dissolution of apatite (LREE rich and HFSE poor) during the formation of the vein melts.

The presence of some amphibole in the veins, and the suggestion that they were derived from partial melting of amphibole peridotites, places the PYV source within the lithosphere, in agreement with our earlier comment about the uniform 143Nd/144Nd ratios and with the suggestion of Stein & Goldstein (1996)Go that all Phanerozoic Arabian basalts were derived from lithospheric mantle sources. This contrasts with the model proposed by Shaw et al. (2003)Go for the northeastern Jordanian basalts, in which they hypothesized that the enriched magmas (basanites) were derived from an asthenospheric OIB-type source. The main disagreement between the two models concerns lithospheric thickness, and whether or not garnet-bearing mantle rocks are stable at lithospheric pressures. As stated above, in disagreement with Shaw et al. (2003)Go and McGuire (1988)Go, we think that the Arabian–Nubian lithosphere in our study area is thicker than 100 km.

Timing of formation of the two sources
The model developed for the evolution of the sources of the Golan basanites and alkali basalts includes two main stages:

(1) Late Proterozoic subduction-related metasomatism of the Arabian–Nubian lithospheric mantle, which resulted in the formation of amphibole peridotites;

(2) lithospheric differentiation processes, including dehydration and low-degree partial melting of these amphibole peridotites that resulted in the formation of the APR and the PYV sources, respectively.

These events could be related to one of a number of thermal or extension events that the Arabian–Nubian lithosphere has experienced since its formation; these could have occurred either at the end of the Late Proterozoic or at any time during the Phanerozoic. The correlation of 87Sr/86Sr with Rb/Sr ratios in the Golan–Galilee basalts (Fig. 10c) cannot be used as an isochron to date these events, as the linear trend observed is interpreted here as a mixing line. However, if the relatively uniform initial Nd isotope ratios (Fig. 9a) of both basanites and alkali basalts are taken as evidence for derivation of their sources from a single mantle protolith, and assuming that both events (dehydration and melting) occurred simultaneously, then the different Sr isotopes compositions of the basanites and alkali basalts could reflect differences in the Rb/Sr ratios of their respective sources acquired during the dehydration and melting processes. Using a two-point model-age calculation for the two sources (PYV and the APR; Tables 4 and 5, Fig. 10c), and using the respective isotopic compositions of Wy-337 and Wy-52 (there are no isotopic data for Wy-16), yields a differentiation age of 313 Ma. Forcing the model source composition to produce melt Rb and Sr concentrations identical to those of Wy-16 and Wy-337 yields older ages (517 Ma). The general slope of the data in Fig. 10c defines an intermediate slope, which, following the same reasoning, yields an age of about 400 Ma. Therefore, we suggest a scenario in which amphibole–garnet peridotites, which were produced in the Late Proterozoic, underwent dehydration and/or partial melting during the Paleozoic, resulting in the formation of the PYV and APR sources. This differentiation event did not produce a large difference in the Sm/Nd ratios of the APR and the PYV sources (0·78 and 0·66, PM-normalized ratios; Tables 4 and 5, respectively), explaining the relative uniformity in the 143Nd/144Nd ratios of the basalts.

Although there is no evidence of magmatic activity in the region during Late Paleozoic times, this was nevertheless a period of major uplift and regional erosion (Garfunkel, 1998Go). A major thermal event in the Israel–Sinai area during the Late Devonian is also implied by zircon fission-track ages (Kohn et al., 1992Go), K–Ar ages and Rb–Sr isochrons for Cambrian illites (Segev & Steinitz, 1986Go), and a Nd isotope isochron age of 400 Ma calculated for a suite of pyroxenite xenoliths from northern Israel (Stein, 1987Go).

Lithospheric mantle heterogeneity
The main message of the above discussion is that the spectrum of compositions in the Golan–Galilee basanites and alkali basalts reflects mixing of melts derived from garnet pyroxenite veins and amphibole–garnet peridotites, respectively. We suggest that the proportion of PYV melt in the mixed product is dependent on the relative proportion of veins in the source; the greater the proportion of veins, the larger the contribution of their melt to the mixed magma and the more enriched the magma is. We, therefore, envisage that the lithospheric mantle beneath northern Israel is composed of several domains that are no more than tens of kilometers in size. The lithospheric domains that supplied the magmas to the southern Golan, as well as to the Lower Galilee, are vein poor, whereas those supplying magmas to the Upper Galilee and to the northern and eastern Golan are vein rich. The Yizre'el Valley, south of the Lower Galilee, was also supplied by magmas from a vein-rich lithospheric domain, albeit at an earlier time (Middle Miocene, Weinstein, 2000Go). According to our model, the change in magma composition during the Plio-Pleistocene reflects the migration of melting from vein-poor to vein-rich lithospheric domains.

Most of the lavas from the Lower Galilee and southern Golan are compositionally indistinguishable (e.g. Figs 2, 5c and 7). We, therefore, suggest that the southern and western Golan and the Lower Galilee received their magmas from a common vein-poor mantle domain. A single lithospheric domain underlying both the southern Golan and the Lower Galilee is unlikely, as the two regions lie on opposite sides of the Dead Sea Transform, and restoring the 105 km left-lateral displacement along the transform (Freund, 1978Go) places the southern Golan and its underlying lithospheric mantle at the latitude of Jerusalem, clearly disconnected from the Lower Galilee. We, therefore, propose that both areas were fed by magmas that originated from a vein-poor lithospheric mantle domain in the sub-Golan lithosphere. We prefer the sub-Golan over the Lower Galilee as the source of the magmas for the following reasons. First, in the Golan the APR magmatism continued longer (Early to Late Pliocene) than in the Lower Galilee (Early Pliocene only). Second, the volume of APR magmatism is larger in the Golan, including both the Early Pliocene in the south and the Late Pliocene in the central Golan. This volume may be even larger if the basalts to the south (northern Jordan) and to the east (Syria) are included. Third, the southern Golan lavas are more uniform alkali basalts, whereas those of the Galilee also include more enriched samples (e.g. hawaiites and even one basanite, Fig. 2). It could be that the lithospheric domains to the west of the DST are mainly vein rich, manifested by the Upper Galilee basanites and by the Miocene basanites and nephelinites from the Yizre'el Valley (Weinstein, 2000Go), and that some of these magmas mixed with sub-Golan-derived APR magmas during the Early Pliocene.

In contrast to the proposed shared APR lithospheric domain, vein-rich domains should be envisaged as isolated pockets, with different compositions under Galilee and under the northern Golan, as suggested by the relatively wide compositional spectrum in diagrams such as La/Nb and La/Sm vs La (see Fig. 11). This is in accordance with the 105 km displacement along the DST, as restoring for the displacement will place these two areas far away from each other.

Cause of Plio-Pleistocene partial melting
The Levant area, including Israel, has experienced several volcanic episodes during the Phanerozoic. These include post-orogenic magmatic activity at the end of the Pan-African period, during the Late Proterozoic to Early Cambrian, and intra-plate basaltic volcanism during the Triassic, the Late Jurassic–Early Cretaceous, the Late Cretaceous and the Neogene–Quaternary (e.g. Garfunkel, 1989Go). The hydrous nature of the Arabian–Nubian lithospheric mantle, proposed here and elsewhere (e.g. Stein et al., 1997Go; Baker et al., 1997), is consistent with the common occurrence of hydrous mantle xenoliths in Arabian basalts (e.g. Mittlefehldt, 1984Go; Stein et al., 1993Go; Chazot et al., 1996Go; Downes et al., 2004Go; J. Baker, personal communication, 2005) and provides a reasonable explanation for these periodic partial melting events, as the wet peridotite solidus is about 300–400°C lower than the dry solidus in the pressure range 20–30 kbar (Fig. 14 and references therein). It is clear from Fig. 14 that under the PT conditions of the Arabian lithosphere (Stein et al., 1993Go), a volatile-saturated lithosphere should be partially molten at pressures in excess of about 25 kbar. Although we cannot disregard the possibility of free volatiles in the source, it is more reasonable to assume fluid-absent conditions (e.g. the H2O-undersaturated or the CO2–H2O-undersaturated solidus curves in Fig. 14). Under such conditions no melting is expected at pressures lower than about 30 kbar. However, lithospheric rocks are very close to their solidus.


Figure 14
View larger version (26K):
[in this window]
[in a new window]
 
Fig. 14. Geotherms [dashed curves; Falloon & Green (1990)Go for oceanic and shield geotherms and Stein et al. (1993)Go for that of the Arabian–Nubian massif] and solidus curves (solid lines) of dry and H2O–CO2-saturated and -undersaturated peridotites (Green, 1973Go; Olafsson & Eggler, 1983Go; Falloon & Green, 1989Go, 1990Go).

 
The solidus curves illustrated in Fig. 14 are for peridotites. Wet melting experiments are limited to peridotites and there are almost no data for melting of amphibole pyroxenites (the ‘veins’). For dry pyroxenites, experimental data show that their solidi are either lower by 50–200°C (Green & Ringwood, 1967Go; Chapman, 1975Go; Thompson, 1975Go) or similar to those of dry peridotites (Ito & Kennedy, 1973Go; see compilation by Hirschmann & Stolper, 1996Go). Kogiso et al. (2004)Go have recently demonstrated that pyroxenite solidi are actually higher than those of peridotites with the same alkali contents or Mg-number. However, they noted that, as most pyroxenites have higher alkali contents than peridotites, the actual solidus temperatures of mantle pyroxenites would probably be lower than those of peridotites. Experimental data for phlogopite pyroxenites show that at 20–30 kbar the solidus temperatures are lower than those of phlogopite peridotites (Lloyd et al., 1985Go; Menner & Dunn, 1995Go). Therefore, it is probable that the solidus of amphibole pyroxenite is also somewhat lower than that of amphibole peridotite. Nevertheless, in our study we assume similar solidus temperatures for amphibole peridotite and amphibole pyroxenite.

In summary, under the geothermal conditions of the Arabian–Nubian lithosphere, the sources for the Neogene–Quaternary magmas—both peridotites and pyroxenites—were probably close to their solidus, but still within the sub-solidus field. Two mechanisms have been proposed as possible causes for deviation from these conditions to produce partial melting:

(1) decompression as a result of extension (e.g. Bohannon et al., 1989Go; Stein et al., 1993Go), associated with rifting along the Red Sea;

(2) thermal perturbation (e.g. Garfunkel, 1989Go; Camp & Roobol, 1992Go), possibly in connection with mantle plume activity under Afar, further to the south.

Unlike the Middle Miocene period, when lithospheric extension resulted in the formation of large sedimentary basins throughout the Levant, including the Lower Galilee–Yizre'el Valley (e.g. Shaliv, 1991Go) and the southern Golan basins, there is no evidence for significant extension in the region during Plio-Pleistocene times, apart from localized pull-apart basins along the DST. Also, the solidus curves of H2O- or CO2–H2O-undersaturated peridotites are very steep in the pressure range 20–30 kbar (Fig. 14; see also compilation by McKenzie & O'Nions, 1995Go), throwing doubt on the efficiency of decompression as a major cause for melting in this case. Therefore, we prefer the thermal perturbation scenario. Thermal perturbation means relatively high temperatures in the sub-lithospheric mantle. High temperatures may be caused by a (small) hotspot (Garfunkel, 1989Go) or by a hot asthenospheric lobe, channeled laterally from the Afar area in the south (Camp & Roobol, 1992Go), similar to the hot asthenospheric mantle ‘fingers’ that have been imaged by seismic tomography beneath the volcanic fields of central and western Europe (Granet et al., 1995Go; Wilson & Patterson, 2001Go) and postulated to occur beneath the Cenozoic fields of central and northern Africa (Wilson & Guiraud, 1998Go). Whatever the exact scenario, we suggest that lithospheric melting was caused by heat transferred by conduction from thermally anomalous zone within the sub-lithospheric mantle.

We wish to emphasize that, in our model, partial melting was not triggered or immediately preceded by metasomatic events, as suggested by Bailey (1982)Go or Wedepohl et al. (1994)Go. Instead we propose that the lithosphere was pre-conditioned by metasomatic ‘events’ that ended hundreds of millions of years earlier (Late Proterozoic). There is no need for introduction of extra-lithospheric fluids and incompatible elements after the stabilization of the lithosphere, as is often suggested for intra-continental alkali magmatism (e.g. Menzies et al., 1987Go; Menzies & Halliday, 1988Go; Wilson & Guiraud, 1992Go) and for kimberlites (Roden & Murthy, 1985Go). This also means that the recent (Neogene) metasomatism identified by Henjes-Kunst et al. (1990)Go in xenoliths from some Arabian volcanic fields was insignificant in affecting the geochemistry of the mantle sources of Golan–Galilee basalts.

What is the reason for the south–north migration of volcanic activity in the Golan and Galilee? One explanation could be that this is the result of the northward propagation of a hot asthenospheric mantle lobe (e.g. Camp & Roobol, 1992Go). Alternatively, it could be a mantle source effect, reflecting the location of the vein-rich domains at somewhat shallower depths than the vein-poor domains in the lithospheric mantle, as a result of the early partial melting that probably left a relatively refractory zone at the base of the lithosphere under the vein-rich areas. Therefore, even if sub-lithospheric heating was equally distributed beneath the regional lithosphere, the vein-rich domains, being further away from the heat source, reacted later than the vein-poor domains. The first explanation (the hot asthenosphere propagation scenario) is supported by the fact that vein-rich melting already occurred during the Miocene under the Yizre'el Valley (Weinstein, 2000Go), to the south of the main Plio-Pleistocene activity. However, the Miocene volcanism could alternatively be the result of the Miocene extension events, which were concentrated in this area (Shaliv, 1991Go). On the other hand, a scenario of a propagation of a hot asthenosphere lobe cannot be reconciled with the observations from the central Golan, where concentrations of incompatible elements and P2O5 increase with time (Pliocene to Pleistocene), although with only a small change in location (about 10 km or less, dependent on where the Pliocene basalts erupted; Fig. 8). Similarly, during the Miocene, most basanites were erupted during the later stages of the period (Weinstein, 2000Go), as would be expected in the case of shallower vein-rich domains. Therefore, we favor the shallower veins explanation, although the asthenosphere lobe propagation scenario cannot be ruled out.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
(1) The composition of Plio-Pleistocene basalts from the Golan and Galilee, Israel (northwestern Arabia), varies with time and location. Alkali basalts erupted in the south during the Early Pliocene, whereas during the Late Pliocene (west of the Dead Sea Transform) and Pleistocene (east of the transform) magmatic activity was concentrated in the north and was mainly basanitic.

(2) 87Sr/86Sr ratios vary between 0·7031 and 0·7034, and correlate negatively with incompatible element concentrations and positively with Rb/Sr ratios. This, together with the upward-convex patterns in plots of incompatible element ratios vs concentration (e.g. Nb/Zr vs Nb, Nd/Sm vs Nd and La/Nb vs La) and other evidence, suggests that the main control on magma composition is binary mixing of melts derived from two end-member sources.

(3) The alkali basalt end-member melt is characterized by relatively high Sr/Ba and Nb/Ba and by low Sm/Nd and relatively high Tb/Yb. Therefore, its source is identified as an amphibole- and garnet-bearing peridotite (APR), limiting its depth to 60–90 km and implying a lithospheric origin for these magmas. The negative anomalies of Rb in the alkali basalts are attributed to a partial dehydration process that shaped the final composition of the APR source.

(4) The basanitic end-member melt is characterized by very high incompatible element concentrations, negative anomalies of K and Rb (the latter stronger than in the alkali basalts), lower Sm/Nd than the alkali basalts and very high Tb/Yb. We propose that these characteristics were produced by a two-stage melting process and that the enriched source of the basanites is pyroxenite veins, hosted by the amphibole peridotites.

(5) The formation of the enriched source of the basanites involves the following stages: (a) early melting of lithospheric amphibole peridotite; (b) solidification of the partial melts to form pyroxenite veins (PYV) within the amphibole peridotites; (c) release of hydrous fluids during solidification of the veins and removal of Rb and K from the veins; (d) Plio-Pleistocene partial melting of the veins (plus reaction with the peridotitic wall-rock) to produce the basanitic end-member.

(6) Based on the above we envisage that the lithosphere beneath northeastern Israel is composed of domains (tens of kilometers scale) that differ in their content of pyroxenite veins.

(7) The uniform 143Nd/144Nd ratios may suggest that, in the past, the two sources shared a common source that differentiated by dehydration and partial melting, probably as a result of a Paleozoic thermal event.

(8) Melting during the Plio-Pleistocene probably occurred as a result of heating of the base of the lithosphere by sub-lithospheric, anomalously hot, mantle.

(9) This thermal event migrated from south to north, affecting first (Early Pliocene) vein-poor lithospheric domains and subsequently (Late Pliocene in the Galilee and Pleistocene in the Golan) shallower vein-rich domains.


    ACKNOWLEDGEMENTS
 
This work could not have been accomplished without the help and advice of G. Wörner, who kindly allowed us to use the analytical facilities of the Institute of Geochemistry at the University of Göttingen. We are grateful to the staff of the Geological Survey of Israel, including Y. Mizrahi, S. Ashkenazi, M. Arnon, E. Ram and others who assisted in the field work. The staff of the Institute of Mineralogy at the University of Heidelberg and the Institute of Mineralogy at the Gutenberg University, Mainz, helped in the analytical work. We wish to express special thanks to K. Simon and A. Eisenhauer for their support and technical help, and to A. Heimann for the K–Ar dating. G. Steinitz, Z. Garfunkel, D. Mor, G. Shaliv, A. Rocholl and K. Munker assisted with information and fruitful discussions. Reviews by H. Downes, A. Le Roex and J. Baker were of great help. This study was supported by the Israel Science Foundation (ISF Grant 50399 to M. Stein and Z. Garfunkel).


    FOOTNOTES
 
{dagger} Present address: Geological Survey of Israel, Jerusalem, Israel. Back


* Corresponding author. Telephone: 972-54-520-4873. Fax: 972-3-534-4430. E-mail: weinsty{at}mail.biu.ac.il


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 SAMPLING AND ANALYTICAL METHODS
 PETROGRAPHY AND CHEMISTRY OF...
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Adam, J., Green, T. H. & Sie, S. H. (1993). Proton microprobe determined partitioning of Rb, Sr, Ba, Y, Zr, Nb and Ta between experimentally produced amphiboles and silicate melts with variable F content. Chemical Geology 109, 29–49.[CrossRef][Web of Science]

Altherr, R., Henjes-Kunst, F. & Baumann, A. (1990). Asthenosphere versus lithosphere as possible sources for basaltic magmas erupted during formation of the Red Sea: constraints from Sr, Pb and Nd isotopes. Earth and Planetary Science Letters 96, 269–286.[CrossRef][Web of Science]

Bailey, D. K. (1982). Mantle metasomatism—continuing chemical change within the earth. Nature 296, 525–530.[CrossRef]

Baker, J. A., Menzies, M. A., Thirlwall, M. F. & Macpherson, C. J. (1997). Petrogenesis of Quaternary intraplate volcanism, Sana'a, Yemen: implications for plume–lithosphere interaction and polybaric melt hybridization. Journal of Petrology 38, 1359–1390.[CrossRef][Web of Science]

Barberi, F., Capaldi, G., Gasperini, P., Marinelli, G., Santacroce, R., Scandone, R., Treuil, M. & Varet, J. (1980). Recent basaltic volcanism of Jordan and its implications on the geodynamic history of the Dead Sea shear zone. In: Carrelli, A. (ed.) Geodynamic Evolution of the Afro-Arabian Rift System. Atti die Convegni Lincei 47, 667–683.

Bartov, Y., Steinitz, G., Eyal, M. & Eyal, Y. (1980). Sinistral movement along the Gulf of Aqaba (Elat)—its age and relation to the opening of the Red Sea. Nature 285, 220–221.[CrossRef]

Beattie, P. (1993). Uranium–thorium disequilibria and partitioning on melting of garnet peridotite. Nature 363, 63–65.[CrossRef]

Beattie, P. (1994). What determines the values of mineral–melt partition coefficients? Mineralogical Magazine 58A, 63–64.

Bohannon, R. G., Naeser, C. W., Schmidt, D. L. & Zimmermann, R. A. (1989). The timing of uplift, volcanism and rifting peripheral to the Red Sea: a case for passive rifting? Journal of Geophysical Research 94B, 1683–1701.

Bonen, D. (1980). The Mezozoic basalts of Israel. Ph.D. thesis, The Hebrew University, Jerusalem (in Hebrew).

Bradshaw, T. K., Hawkesworth, C. J. & Gallagher, K. (1993). Basaltic volcanism in the southern Basin and Range: no role for a mantle plume. Earth and Planetary Science Letters 116, 45–62.[CrossRef][Web of Science]

Brey, G., Brice, W. R., Ellis, D. J., Green, D. H., Harris, K. L. & Ryabchikov, I. D. (1983). Pyroxene–carbonate reactions in the upper mantle. Earth and Planetary Science Letters 62, 63–74.[CrossRef][Web of Science]

Brueckner, H. K., Zindler, A., Seyler, M. & Bonatti, E. (1988). Zabargad and the isotopic evolution of the sub-Red-Sea mantle and crust. Tectonophysics 150(1–2), 163–176.[CrossRef][Web of Science]

Camp, V. E. & Roobol, M. J. (1992). Upwelling asthenosphere beneath western Arabia and its regional implications. Journal of Geophysical Research 97B, 15255–15271.[CrossRef]

Chapin, C. E. (1979). Evolution of the Rio Grande rift—a summary. In: Rieker, R. E. (ed.) Rio Grande Rift: Tectonics and Magmatism. Washington, DC: American Geophysical Union, pp. 1–5.

Chapman, N. A. (1975). An experimental study of spinel clinopyroxenite xenoliths from the Duncansby Ness vent, Caithness, Scotland. Contributions to Mineralogy and Petrology 51, 223–230.[CrossRef][Web of Science]

Chazot, G., Menzies, M. & Harte, B. (1996). Determination of partition coefficients between apatite, cpx, amphibole and melt in natural spinel lherzolites from Yemen: implications for wet melting of the lithosphere mantle. Geochimica et Cosmochimica Acta 60, 423–437.[CrossRef][Web of Science]

Cohen, R. S., O'Nions, K. R. & Dawson, J. B. (1984). Isotope geochemistry of xenoliths from East Africa: implications for development of mantle reservoirs and their interaction. Earth and Planetary Science Letters 68, 209–220.[CrossRef][Web of Science]

Dafny, E., Gvirtzman, H., Burg, A. & Fleischer, L. (2003). The hydrogeology of the Golan basalt aquifer, Israel. Israel Journal of Earth Sciences 52, 139–153.[CrossRef]

Dalpé, C. & Baker, D. R. (1994). Partition coefficients for REE between calcic amphibole and Ti-rich basanitic glass at 1·5 GPa, 1100°C. Goldschmidt Conference, Edinburgh. Mineralogical Magazine 58A, 207–208.

Deniel, C., Vidal, P., Coulon, C., Vellutini, P.-J. & Piguet, P. (1994). Temporal evolution of mantle sources during continental rifting: the volcanism of Djibouti (Afar). Journal of Geophysical Research 99B, 2853–2869.[CrossRef]

Downes, H., Beard, A. & Hinton, R. (2004). Natural experimental charges: an ion-microprobe study of trace element distribution coefficients in glass-rich hornblendite and clinopyroxenite xenoliths. Lithos 75, 1–17.[CrossRef][Web of Science]

Dubertret, L. (1955). The ossiferous seams of Khirbet El Umbachi and of Hebariye (Safa). Annales Archéologiques de Syrie IV–V, 59–75.

Duyverman, H. J., Harris, N. B. W. & Hawkesworth, C. J. (1982). Crustal accretion in the Pan-African: Nd and Sr isotope evidence from the Arabian shield. Earth and Planetary Science Letters 59, 315–326.[CrossRef][Web of Science]

Eckstein, Y. (1979). Review of heat flow data from the eastern Mediterranean region. Pure and Applied Geophysics 117, 150–159.[CrossRef][Web of Science]

Falloon, T. J. & Green, D. H. (1988). Anhydrous partial melting of peridotite from 8 to 35 kb and the petrogenesis of MORB. Journal of Petrology, Special Lithosphere Issue, 379–414.

Falloon, T. J. & Green, D. H. (1989). The solidus of carbonated, fertile peridotite. Earth and Planetary Science Letters 94, 364–370.[CrossRef][Web of Science]

Falloon, T. J. & Green, D. H. (1990). Solidus of carbonated fertile peridotite under fluid saturated conditions. Geology 18, 195–199.[Abstract/Free Full Text]

Francis, D. (1991). Some implications of xenolith glasses for the mantle sources of alkaline mafic magmas. Contributions to Mineralogy and Petrology 108, 175–180.[CrossRef][Web of Science]

Francis, D. & Ludden, J. (1995). The signature of amphibole in mafic alkaline lavas, a study in the northern Canadian Cordillera. Journal of Petrology 36, 1171–1191.[Abstract/Free Full Text]

Freund, R. (1978). The concept of a sinistral megashear. In: Serruya, C. (ed.) Lake Kinneret. Monographiae Biologicae 32, 27–31.

Freund, R., Garfunkel, Z., Zak, I., Goldberg, M., Weissbrod, T. & Derin, B. (1970). The shear along the Dead Sea rift. Philosophical Transactions of the Royal Society of London, Series A 267, 107–130.[CrossRef]

Gallagher, K. & Hawkesworth, C. J. (1992). Dehydration melting and the generation of continental flood basalts. Nature 358, 57–59.[CrossRef]

Gao, S., Ling, W., Qiu, Y., Lian, Z., Hartmann, G. & Simon, K. (1999). Contrasting geochemical and Sm–Nd isotopic compositions of Archean metasediments from the Kongling high-grade terrain of the Yangtze craton: evidence for cratonic evolution and redistribution of REE during crustal anatexis. Geochimica et Cosmochimica Acta 63, 2071–2088.[CrossRef][Web of Science]

Garfunkel, Z. (1981). Internal structure of the Dead Sea leaky transform (rift) in relation to plate kinematics. Tectonophysics 80, 81–108.[CrossRef][Web of Science]

Garfunkel, Z. (1989). Tectonic setting of Phanerozoic magmatism in Israel. Israel Journal of Earth Sciences 38, 51–74.

Garfunkel, Z. (1998). Constraints on the origin and history of the Eastern Mediterranean basin. Teconophysics 298: 5–35.

Garfunkel, Z. & Ben Avraham, Z. (2001). Basins along the Dead Sea Transform. In: Ziegler, P. A., Cavazza, W., Robertson, A. H. F. & Crasquin-Soleau, S. (eds) Peri-Tethys Memoir 6: Peri-Tethyan Rifts/Wrench Basins and Passive Margins. Mémoires du Muséum National d'Histoire Naturelle 186, 607–627.

Garfunkel, Z., Zak, I. & Freund, R. (1981). Active faulting in the Dead Sea rift. Tectonophysics 80, 1–26.[CrossRef][Web of Science]

Gass, I. G. (1970). Tectonic and magmatic evolution of the Afro-Arabian dome. In: Clifford, T. N. & Gass, I. G. (eds) African Magmatism and Tectonics. Edinburgh: Oliver & Boyd, pp. 285–300.

Gazit, O. (2005). Mafic granulite xenoliths from Qarnei Hittin and the evolution of the lower crust beneath northern Israel. Geological Survey of Israel Report GSI/11/05, 67 pp. (in Hebrew, with English abstract).

Govindaraju, K. (1989). Compilation of working values and sample description for 272 geostandards. Geostandard Newsletters, Special Issue 13, 1–113.

Granet, M., Wilson, M. & Achauer, U.(1995). Imaging a mantle plume beneath the French Massif Central. Earth and Planetary Science Letters 136, 281–296.[CrossRef][Web of Science]

Green, D. H. (1973). Experimental melting studies on a model upper mantle composition at high pressure under water-saturated and water-undersaturated conditions. Earth and Planetary Science Letters 19, 37–53.[CrossRef][Web of Science]

Green, D. H. & Ringwood, A. E. (1967). The genesis of basaltic magmas. Contributions to Mineralogy and Petrology 15, 103–190.[CrossRef]

Green, T. H. (1994). Experimental studies of trace element partitioning applicable to igneous petrogenesis—Sedona 16 years later. Chemical Geology 117, 1–36.[CrossRef][Web of Science]

Green, T. H., Blundy, J. D., Adam, J. & Yaxley, J.M. (2000). SIMS determination of trace element partition coefficients between garnet, clinopyroxene and hydrous basaltic liquids at 2–7·5 GPa and 1080–1200°C. Lithos 53, 165–187.[CrossRef][Web of Science]

Hart, S. R. & Dunn, T. (1993). Experimental cpx–melt partitioning of 24 trace elements. Contributions to Mineralogy and Petrology 113, 1–8.[CrossRef][Web of Science]

Hauri, E. H., Wagner, T. P. & Grove, T. L. (1994). Experimental and natural partitioning of Th, U, Pb and other trace elements between garnet, cpx and basaltic melts. Chemical Geology 117, 149–166.[CrossRef][Web of Science]

Hegner, E. & Pallister, J. S. (1989). Pb, Sr, Nd isotopic characteristics of Tertiary Red Sea rift volcanics from the central Saudi Arabian coastal plain. Journal of Geophysical Research 94B, 7749–7755.

Heimann, A. (1990). The development of the Dead Sea rift and its margins in northern Israel during the Pliocene and the Pleistocene. Geological Survey of Israel Report GSI/28/90, 82 pp. (in Hebrew).

Heimann, A., Steinitz, G., Mor, D. & Shaliv, G. (1996). The Cover Basalt Formation, its age and its regional and tectonic setting: implications from K–Ar and 40Ar/39Ar geochronology. Israel Journal of Earth Sciences 45, 55–71.

Henjes-Kunst, F., Altherr, R. & Baumann, A. (1990). Evolution and composition of the lithospheric mantle underneath the western Arabian peninsula: constraints from Sr–Nd isotope systematics of mantle xenoliths. Contributions to Mineralory and Petrology 105, 460–472.

Hirschmann, M. M. & Stolper, E. M. (1996). A possible role for garnet pyroxenite in the origin of the ‘garnet signature’ in MORB. Contributions to Mineralogy and Petrology 124, 185–208.[CrossRef][Web of Science]

Hofmann, A. W. & Feigenson, M. D. (1983). Case studies on the origin of basalt: 1. Theory and reassessment of Grenada basalts. Contributions to Mineralogy and Petrology 84, 382–389.[CrossRef][Web of Science]

Horn, I., Foley, S. F., Jackson, S. E. & Jenner, G. A. (1994). Experimentally determined partitioning of HFS- and selected transition elements between spinel and basaltic melt. Chemical Geology 117, 193–218.[CrossRef][Web of Science]

Hurwitz, S., Garfunkel, Z., Ben-Gai, Y., Reznikov, M., Rotstein, Y. & Gvirtzman, H. (2002). The tectonic framework of a complex pull-apart basin: seismic reflection observations in the Sea of Galilee, Dead Sea transform. Tectonophysics 359, 289–306.[CrossRef][Web of Science]

Ilani, S., Harlvan, Y., Tarawneh, K., Rabba, I., Weinsberger, R., Ibrahim, K., Peltz, S. & Steinitz, G. (2001). New K–Ar ages of basalts from the Harrat Ash Shaam volcanic field in Jordan: implications for the span and duration of the upper-mantle upwelling beneath the western Arabian plate. Geology 29, 171–174.[Abstract/Free Full Text]

Ito, K. & Kennedy, G. C. (1973). The composition of liquids formed by partial melting of eclogites at high temperatures and pressures. Journal of Geology 82, 383–392.

Jochum, K. P., Seufert, H. M. & Thirlwall, M. F. (1990). High sensitivity Nb analysis by spark-source mass spectrometry (SSMS) and calibration of XRF Nb and Zr. Chemical Geology 81, 1–16.[CrossRef][Web of Science]

Kelemen, P. B., Shimizu, N. & Dunn, T. (1993). Relative depletion of niobium in some arc magmas and the continental crust: partitioning of K, Nb, La and Ce during melt/rock reaction in the upper mantle. Earth and Planetary Science Letters 120, 111–134.[CrossRef][Web of Science]

Kennedy, A. K., Lofgren, G. E. & Wasserburg, G. J. (1993). An experimental study of trace element partitioning between olivine, orthopyroxene and melt in chondrules: equilibrium values and kinetic effects. Earth and Planetary Science Letters 115, 177–195.[CrossRef][Web of Science]

Keppler, H. (1996). Constraints from partitioning experiments on the composition of subduction-zone fluids. Nature 380, 237–240.[CrossRef]

Klein, E. M. & Langmuir, C. H. (1987). Global correlations of ocean ridge basalt chemistry with axial depth and crustal thickness. Journal of Geophysical Research 92(B8), 8089–8115.

Kogiso, T., Hirschmann, M. M. & Pertermann, M. (2004). High-pressure partial melting of mafic lithologies in the mantle. Journal of Petrology 45, 2407–2422.[Abstract/Free Full Text]

Kohn, B. P., Eyal, M. & Feinstein, S. (1992). A major Late Devonian–Early Carboniferous (Hercynian) thermotectonic event at the NW margin of the Arabian–Nubian shield: evidence from zircon track dating. Tectonics 11, 1018–1027.[Web of Science]

Kramar, U. & Puchelt, H. (1982). Reproducibility tests for INAA determinations with AGV-1, BCR-1 and GSP-1 and new data for 17 geochemical reference materials. Geostandards Newsletter 6, 221–227.

LaTourrette, T., Hervig, R. L. & Holloway, J. R. (1995). Trace element partitioning between amphibole, phlogopite and basanitic melt. Earth and Planetary Science Letters 135, 13–30.[CrossRef][Web of Science]

Le Bas, M. J., Le Maitre, R. W., Streckeisen, A. & Zanettin, B. (1986). A chemical classification of volcanic rocks based on the total alkali–silica diagram. Journal of Petrology 27, 745–750.[Abstract/Free Full Text]

Liotard, J. M., Briot, D. & Boivin, P. (1988). Petrological and geochemical relationships between pyroxene megacrysts and associated alkali-basalts from Massif Central (France). Contributions to Mineralogy and Petrology 98, 81–90.[CrossRef][Web of Science]

Lloyd, F. E., Arima, M. & Edgar, A. D. (1985). Partial melting of a phlogopite-clinopyroxenite nodule from southwest Uganda: an experimental study bearing on the origin of highly potassic continental rift volcanics. Contributions to Mineralogy and Petrology 91, 321–329.[CrossRef][Web of Science]

McDonough, W. F. & Sun, S.-S. (1995). The composition of the Earth. Chemical Geology 120, 223–253.[CrossRef][Web of Science]

McGuire, A. V. (1988). Petrology of mantle xenoliths from Harrat al Kishb: the mantle beneath western Saudi Arabia. Journal of Petrology 29, 73–92.[Abstract/Free Full Text]

McGuire, A. V. & Bohannon, R. G. (1989). Timing of mantle upwelling: evidence for a passive origin of the Red Sea rift. Journal of Geophysical Research 94B, 1677–1682.

McGuire, A. V. & Stern, R. J. (1993). Granulite xenoliths from western Saudi Arabia: the lower crust of the Late Precambrian Arabian–Nubian shield. Contributions to Mineralogy and Petrology 114, 395–408.[CrossRef][Web of Science]

McKenzie, D. & O'Nions, R. K. (1995). The source regions of OIB. Journal of Petrology 36, 139–159.

Menner, A. V. & Dunn, T. (1995). Amphibole and phlogopite stability in an initially amphibole-bearing spinel peridotite under water-undersaturated conditions: 1 to 2·5 GPa. EOS Transactions, American Geophysical Union, Fall Meeting, 697.

Menzies, M. A. & Halliday, A. (1988). Lithosphere mantle domains beneath the Archean and Proterozoic crust of Scotland. Journal of Petrology, Special Lithosphere Issue, 275–302.

Menzies, M. A., Rogers, N., Tindle, A. & Hawkesworth, C. J. (1987). Metasomatic and enrichment processes in lithospheric peridotites, an effect of asthenosphere–lithosphere interaction. In: Menzies, M. A. & Hawkesworth, C. J. (eds) Mantle Metasomatism. London: Academic Press, pp. 313–361,

Menzies, M. A., Baker, J., Chazot, G. & Al-Kadasi, M. (1997). Evolution of the Red Sea volcanic margin, western Yemen, In: Mahoney, J. & Coffin, M. (eds) Large Igneous Provinces: Continental, Oceanic and Planetary Flood Volcanism. Geophysical Monograph, American Geophysical Union 100, 29–43.

Mittlefehldt, D. W. (1984). Genesis of cpx–amphibole xenoliths from Birket Ram: trace element and petrologic constraints. Contributions to Mineralogy and Petrology 88, 280–287.[CrossRef][Web of Science]

Mittlefehldt, D. W. & Slager, Y. (1986). Petrology of the basalts and gabbros from the Zemah-1 drillhole, Jordan Rift Valley. Israel Journal of Earth Science 35, 10–22.

Mor, D. (1993). A time-table for the Levant volcanic province, according to K–Ar dating in the Golan Heights, Israel. Journal of African Earth Science 16(3), 223–234.

Mysen, B. O. & Boettcher, A. L. (1975). Melting of a hydrous mantle: I. Phase relations of natural peridotite at high pressures and temperatures with controlled activities of water, carbon dioxide and hydrogen. Journal of Petrology 16, 520–548.[Abstract/Free Full Text]

Nasir, S. (1992). The lithosphere beneath the northwestern part of the Arabian plate (Jordan): evidence from xenoliths and geophysics. Tectonophysics 201, 357–370.[CrossRef][Web of Science]

Navon, O. & Stolper, E. (1987). Geochemical consequences of melt percolation: the upper mantle as a chromatographic column. Journal of Geology 95, 285–307.[Web of Science]

Nielsen, R. L., Gallahan, W. E. & Newberger, F. (1992). Experimentally determined mineral–melt partition coefficients for Sc, Y and REE for olivine, orthopyroxene, pigeonite, magnetite and ilmenite. Contributions to Mineralogy and Petrology 110, 488–499.[CrossRef][Web of Science]

O'Hara, M. J. (1968). The bearing of phase equilibria studies in synthetic and natural systems on the origin and evolution of basic and ultrabasic rocks. Earth-Science Reviews 4, 69–134.

Olafsson, M. & Eggler, D. H. (1983). Phase relations of amphibole, amphibole–carbonate and phlogopite–carbonate peridotite: petrologic constraints on the asthenosphere. Earth and Planetary Science Letters 64, 305–315.[CrossRef][Web of Science]

Omar, G. I. & Steckler, M. S. (1995). Fission track evidence on the initial rifting of the Red Sea: two pulses, no propagation. Science 270, 1341–1344.[Abstract/Free Full Text]

Pallister, J. S. (1987). Magmatic history of Red Sea rifting: perspective from the central Saudi Arabian coastal plain. Geological Society of America Bulletin 98, 400–417.[Abstract/Free Full Text]

Quennell, A. M. (1958). The structural and geomorphic evolution of the Dead Sea Rift. Quarterly Journal of the Geological Society of London 114, 1–24.

Roden, M. F. & Murthy, V. R. (1985). Mantle metasomatism. Annual Review of Earth and Planetary Sciences 13, 269–296.[CrossRef][Web of Science]

Sack, R. O., Walker, D. & Carmichael, I. S. E. (1987). Experimental petrology of alkalic lavas: constraints on cotectics of multiple saturation in natural basic liquids. Contributions to Mineralogy and Petrology 96, 1–23.[CrossRef][Web of Science]

Schilling, J.-G., Kingsley, R. H., Hanan, B. B. & McCully, B. L. (1992). Nd–Sr–Pb isotopic variations along the Gulf of Aden: evidence for the Afar mantle plume–continental lithosphere interaction. Journal of Geophysical Research 97B, 10927–10966.[CrossRef]

Segev, A. & Steinitz, G. (1986). Dating of epigenetic manganese nodules and reset illites in the Cambrian Timna Formation, southern Israel. Abstracts of the ICOG 6, Cambridge. Terra Cognita 6, 112–113.

Shaliv, G. (1991). Stages in the tectonic and volcanic history of the Neogene basin in the Lower Galilee and the Valleys. Geological Survey of Israel Report GSI/11/91, 94 pp. (in Hebrew).

Shaw, D. M. (1970). Trace element fractionation during anatexis. Geochimica et Cosmochimica Acta 34, 237–243.[CrossRef][Web of Science]

Shaw, J. (2003). Geochemistry of Cenozoic volcanism and Arabian lithospheric mantle in Jordan. Ph.D. thesis, Royal Holloway University of London, 268 pp.

Shaw, J. A., Baker, J. A., Menzies, M. A., Thirlwall, M. F. & Ibrahim, K. M. (2003). Petrogenesis of the largest intraplate volcanic field on the Arabian plate (Jordan): a mixed lithosphere–asthenosphere source activated by lithospheric extension. Journal of Petrology 44, 1657–1679.[Abstract/Free Full Text]

Späth, A., Le Roex, A. P. & Opiyo-Akech, N. (2001). Plume–lithosphere interaction and the origin of continental rift-related alkaline volcanism—the Chyulu Hills volcanic province, southern Kenya. Journal of Petrology 42, 765–787.[Abstract/Free Full Text]

Stein, M. (1987). Ultra-basic xenoliths in Cenozoic and Mesozoic volcanic rocks in Israel. Ph.D. thesis, The Hebrew University, Jerusalem, 117 pp. (in Hebrew).

Stein, M. & Goldstein, S. L. (1996). From plume head to continental lithosphere in the Arabian–Nubian shield. Nature 382, 773–778.[CrossRef]

Stein, M. & Hofmann, A. W. (1992). Fossil plume head beneath the Arabian lithosphere? Earth and Planetary Science Letters 114, 193–209.[CrossRef][Web of Science]

Stein, M., Garfunkel, Z. & Jagoutz, E. (1993). Chronothermometry of peridotitic and pyroxenitic xenoliths: implications for the thermal evolution of the Arabian lithosphere. Geochimica et Cosmochimica Acta 57, 1325–1337.[CrossRef][Web of Science]

Stein, M., Navon, O. & Kessel, R. (1997). Chromatographic metasomatism of the Arabian–Nubian lithosphere. Earth and Planetary Science Letters 152, 75–91.[CrossRef][Web of Science]

Stolper, E. & Newman, S. (1994). The role of water in the petrogenesis of Mariana trough magmas. Earth and Planetary Science Letters 121, 293–325.[CrossRef][Web of Science]

Stosch, H.-G. (1982). REE partitioning between minerals from anhydrous spinel peridotite xenoliths. Geochimica et Cosmochimica Acta 46, 793–811.[CrossRef][Web of Science]

Sun, S.-S. & McDonough, W. F. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In: Saunders, A. D. & Norry, M. J. (eds) Magmatism in the Ocean Basins. Geological Society, London, Special Publications 42, 313–345.

Takahashi, E. & Kushiro, I. (1983). Melting of a dry peridotite at high pressures and basalt magma genesis. American Mineralogist 68, 859–879.[Abstract]

Thompson, R. N. (1975). Is upper mantle P contained in sodic gt? Earth and Planetary Science Letters 26, 417–424.[CrossRef][Web of Science]

Thompson, R. N., Ottley, C. J., Smith, P. M., Pearson, D. G., Dickin, A. P., Morrison, M.A., Leat, P. T. & Gibson, S. A. (2005). Source of the Quaternary alkalic basalts, picrites and basanites of the Potrillo volcanic field, New Mexico, USA: lithosphere or convecting mantle? Journal of Petrology 46(8), 1603–1643.[Abstract/Free Full Text]

Turner, S., Hawkesworth, C., Gallagher, K., Stewart, K., Peate, D. & Mantovani, M. (1996). Mantle plumes, flood basalts and thermal models for melt generation beneath continents: assessment of a conductive heating model and application to the Paraná. Journal of Geophysical Research 101B, 11503–11518.[CrossRef]

Van Westrenen, W., Blundy, J. & Wood, B. (1999). Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. American Mineralogist 84, 838–847.[Abstract]

Vidal, P., Deniel, C., Vellutini, P.-J., Piguet, P., Coulon, C., Vincent, J. & Audin, J. (1991). Changes of mantle sources in the course of a rift evolution: the Afar case. Geophysical Research Letters 18, 1913–1916.[Web of Science]

Volker, F., McCulloch, M. T. & Altherr, R. (1993). Submarine basalts from the Red Sea: new Pb, Sr and Nd isotopic data. Geophysical Research Letters 20, 927–930.[Web of Science]

Wallace, M. E. & Green, D. H. (1988). An experimental determination of primary carbonatite magma composition. Nature 335, 343–346.[CrossRef]

Wedepohl, K. H., Gohn, E. & Hartmann, G. (1994). Cenozoic alkali basaltic magmas of western Germany and their products of differentiation. Contributions to Mineralogy and Petrology 115, 253–278.[CrossRef][Web of Science]

Weinstein, Y. (2000). Spatial and temporal geochemical variability in basin-related volcanism, northern Israel. Journal of African Earth Sciences 30, 865–885.[CrossRef]

Weinstein, Y., Navon, O. & Lang, B. (1994). Fractionation of Pleistocene alkali-basalts from the northern Golan Heights, Israel. Israel Journal of Earth Sciences 43, 63–79.

White, R. S. & McKenzie, D. (1989). Magmatism at rift zones: the generation of volcanic continental margins and flood basalts. Journal of Geophysical Research 94, 7685–7730.

Wilson, M. & Downes, H. (1991). Tertiary–Quaternary extension related alkaline magmatism in western and central Europe. Journal of Petrology 32(4), 811–849.[Abstract/Free Full Text]

Wilson, M. & Guiraud, R. (1992). Magmatism and rifting in western and central Africa, from late Jurassic to recent times. Tectonophysics 213, 203–225.[CrossRef][Web of Science]

Wilson, M. & Guiraud, R. (1998). Late Permian to recent magmatic activity on the Afrian–Arabian margin of Tethys. In: MacGregor, D. S., Moody, R. T. J. & Clark-Lowes, D. D. (eds) Petroleum Geology of North Africa. Geological Society, London, Special Publications 132, 231–263.

Wilson, M. & Patterson, R. (2001). Intra-plate magmatism related to hot fingers in the upper mantle: evidence from the Tertiary–Quaternary volcanic province of western and central Europe. In: Ernst, R. & Buchan, K. (eds) Mantle Plumes: their Identification through Time. Geological Society of America, Special Papers 352, 37–58.

Wörner, G., Zindler, A., Staudigel, H. & Schmincke, H. U. (1986). Sr, Nd, and Pb isotope geochemistry of Tertiary and Quaternary alkaline volcanics from West Germany. Earth and Planetary Science Letters 79, 107–119.[CrossRef][Web of Science]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Eur J MineralHome page
F. Lucassen, G. Franz, R. L. Romer, and P. Dulski
Late Cenozoic xenoliths as a guide to the chemical - isotopic composition and thermal state of the upper mantle under northeast Africa
European Journal of Mineralogy, December 1, 2008; 20(6): 1079 - 1096.
[Abstract] [Full Text] [PDF]


Home page
Geological MagazineHome page
N. A. LEASE and A.-F. M. ABDEL-RAHMAN
The Euphrates volcanic field, northeastern Syria: petrogenesis of Cenozoic basanites and alkali basalts
Geological Magazine, September 1, 2008; 145(5): 685 - 701.
[Abstract] [Full Text] [PDF]


Home page
Geological Society, London, Special PublicationsHome page
M. Ismail, G. Delpech, J.-Y. Cottin, M. Gregoire, B. N. Moine, and A. Bilal
Petrological and geochemical constraints on the composition of the lithospheric mantle beneath the Syrian rift, northern part of the Arabian plate
Geological Society, London, Special Publications, January 1, 2008; 293(1): 223 - 251.
[Abstract] [Full Text] [PDF]


Home page
J PetrologyHome page
M.-S. Krienitz, K. M. Haase, K. Mezger, and M. A. Shaikh-Mashail
Magma Genesis and Mantle Dynamics at the Harrat Ash Shamah Volcanic Field (Southern Syria)
J. Petrology, August 1, 2007; 48(8): 1513 - 1542.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Appendix 1
Right arrow All Versions of this Article:
47/5/1017    most recent
egl003v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (5)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by WEINSTEIN, Y.
Right arrow Articles by STEIN, M.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?