Jurassic back-arc and Cretaceous hot

Lithos 112 (2009) 163–187
Contents lists available at ScienceDirect
Lithos
j o u r n a l h o m e p a g e : w w w. e l s ev i e r. c o m / l o c a t e / l i t h o s
Jurassic back-arc and Cretaceous hot-spot series In the Armenian
ophiolites — Implications for the obduction process
Yann Rolland a,⁎, Ghazar Galoyan a,b, Delphine Bosch c, Marc Sosson a, Michel Corsini a,
Michel Fornari a, Chrystèle Verati a
a
b
c
Géosciences Azur, Université de Nice Sophia Antipolis, CNRS, IRD, Parc Valrose, 06108 Nice cedex 2, France
Institute of Geological Sciences, National Academy of Sciences of Armenia, 24a Baghramian avenue, Yerevan, 375019, Armenia
Géosciences Montpellier, CNRS UMR-5243, Université de Montpellier II, Place E. Bataillon, 34095 Montpellier Cedex 05, France
a r t i c l e
i n f o
Article history:
Received 2 July 2008
Accepted 16 February 2009
Available online 10 March 2009
Keywords:
Nd–Sr–Pb isotopes
Armenian ophiolite
Back-arc
Obduction
Oceanic plateau
Tethys
Lesser Caucasus
a b s t r a c t
The identification of a large OIB-type volcanic sequence on top of an obducted nappe in the Lesser Caucaus of
Armenia helps us explain the obduction processes in the Caucasus region that are related to dramatic change
in the global tectonics of the Tethyan region in the late Lower Cretaceous. The ophiolitic nappe preserves
three distinct magmatic series, obducted in a single tectonic slice over the South Armenian Block during the
Coniacian–Santonian (88–83 Ma), the same time as the Oman ophiolite. Similar geological, petrological,
geochemical and age features for various Armenian ophiolitic massifs (Sevan, Stepanavan, and Vedi) argue
for the presence of a single large obducted ophiolite unit. The ophiolite, shows evidence for a slow-spreading
oceanic environment in Lower to Middle Jurassic. Serpentinites, gabbros and plagiogranites were exhumed
by normal faults, and covered by radiolarites. Few pillow-lava flows have infilled the rift grabens.
The ophiolite lavas have hybrid geochemical composition intermediate between Arc and MORB signatures:
(La/Yb) N = 0.6–0.9; (Nb/Th) N = 0.17–0.57; (143Nd/144Nd) i = 0.51273–0.51291; (87Sr/86Sr) i = 0.70370–
0.70565; (207Pb/204Pb)i = 15.4587–15.5411; (208Pb/204Pb)i = 37.4053–38.2336; (206Pb/204Pb)i = 17.9195–
18.4594. These compositions suggest they were probably formed in a back-arc basin by melting of a shallow
asthenosphere source contaminated by a deeper mantle source modified by subducted slab-derived products.
87
Sr/86Sr ratios and petrological evidence show that these lavas have been intensely altered by mid-oceanic
hydrothermalism as well as by serpentinites, which are interpreted as exhumed mantle peridotites.
The gabbros have almost the same geochemical composition as related pillow-lavas: (La/Yb)N = 0.2–2.3;
(Nb/Th)N = 0.1–2.8; (143Nd/144Nd)i = 0.51264–0.51276; (87Sr/86Sr)i = 0.70386–0.70557; (207Pb/204Pb)i =
15.4888–15.5391; (208Pb/204Pb)i = 37.2729–37.8713; (206Pb/204Pb)i = 17.6296–17.9683. Plagiogranites show
major and trace element features similar to other Neo-Tethyan plagiogranites (La/Yb)N = 1.10–7.92; (Nb/
Th)N = 0.10–0.94; but display a less radiogenic Nd isotopic composition than basalts [(143Nd/144Nd)i =
0.51263] and more radiogenic (87Sr/86Sr)i ratios. This oceanic crust sequence is covered by variable thicknesses
of unaltered pillowed OIB alkaline lavas emplaced in marine conditions. 40Ar/39Ar dating of a single-grain
amphibole phenocryst provides a Lower Cretaceous age of 117.3 ± 0.9 Ma, which confirms a distinct formation
age of the OIB lavas. The geochemical composition of these alkaline lavas is similar to plateau-lavas [(La/
Yb)N = 6–14; (Nb/Th)N = 0.23–0.76; (143Nd/144Nd)i = 0.51262–0.51271; (87Sr/86Sr)i = 0.70338–0.70551;
(207Pb/204Pb)i = 15.5439–15.6158; (208Pb/204Pb)i = 38.3724–39.3623; (206Pb/204Pb)i = 18.4024–19.6744].
They have significantly more radiogenic lead isotopic compositions than ophiolitic rocks, and fit the
geochemical compositions of hot-spot derived lavas mixed with various proportions of oceanic mantle. In
addition, this oceanic + plateau sequence is covered by Upper Cretaceous calc-alkaline lavas: (La/Yb)N = 2.07–
2.31; (Nb/Th)N = 0.08–0.15; (144Nd/143Nd)i = 0.51271–0.51282; (87Sr/86Sr)i = 0.70452–0.70478), which
were likely formed in a supra-subduction zone environment. During the late Lower to early Upper Cretaceous
period, hot-spot related magmatism related to plateau events may have led to significant crustal thickening in
various zones of the Middle-eastern Neotethys. These processes have likely hindered subduction of some of the
hot and thickened oceanic crust segments, and allowed them to be obducted over small continental blocks such
as the South Armenian Block.
© 2009 Elsevier B.V. All rights reserved.
⁎ Corresponding author. Tel.: +33 4 92 07 65 86.
E-mail address: [email protected] (Y. Rolland).
0024-4937/$ – see front matter © 2009 Elsevier B.V. All rights reserved.
doi:10.1016/j.lithos.2009.02.006
164
Y. Rolland et al. / Lithos 112 (2009) 163–187
1. Introduction
The role of Oceanic Plateaus in the obduction processes of oceanic
crust has still not been clearly established. We understand that their
larger crustal thickness and buoyancy as compared to ‘standard’
oceanic crust does not allow them to subduct, in particular when they
reach subduction zones soon after their formation (e.g., Ben-Avraham
et al., 1981; Cloos, 1993; Abbot and Mooney, 1995; Abbot et al., 1997;
Kerr and Mahoney, 2007). However, the reasons for oceanic crust
obduction onto continental margins are still debated: (i) is obduction
driven by subduction of continental crust? Or (ii) does it result from
the intrinsic nature of the oceanic crust? In the first case, ophiolites are
obducted due to the mechanical coupling of continental crust with the
dense subducting slab (e.g., O'Brien et al., 2001; Guillot et al., 2003).
Continental subduction may be facilitated by the thinned margins of
continental domains following earlier phases of divergence rifting
that precede oceanic crust emplacement (Guillot and Allemand,
2002). In the second case, a lower density of oceanic lithosphere
might result from intra-oceanic hot-spot and magmatic arc events,
which will lead to crustal thickening and a decrease in lithosphere
density (Cloos, 1993; Abbot and Mooney, 1995). The emplacement of
oceanic plateaus has a great influence either on the slab dip, but also
on the cessation of subduction and on the onset of obduction as is
proposed for the Ontong–Java plateau (Petterson et al., 1997). The
ability of an oceanic plateau to resist subduction and eventually be
transported onto continental crust depends on both crustal thickness
and plateau age (Kerr and Mahoney, 2007). The older a plateau, the
cooler and thus the less buoyant it will be. Alternative hypotheses for
obduction involve rapid inversion of tectonic plate motions and rapid
continental convergence (e.g., Agard et al., 2007). Obduction is
ascribed to the presence of young oceanic crust in the hanging-wall
of the subduction zone, as a result of subduction initiation at the MidOceanic Ridge (e.g., Boudier et al., 1988; Nicolas, 1989); or to scalping
of oceanic lithosphere (e.g., Agard et al., 2007 and references therein).
The case of Armenian ophiolites (Lesser Caucasus) is peculiar as
recent investigations (Galoyan et al., 2007, 2009; Rolland et al., in
press) have shown the presence of slow-spreading ophiolites in
several locations. Further, the ophiolites were tectonically transported
above the South Armenian Bloc or SAB (Zakariadze et al., 1983).
Although some blueschists are locally found, these affect oceanic
Fig. 1. Tectonic map of the Middle East — Caucasus area, with main blocks and suture zones, after Avagyan et al. (2005), modified.
Y. Rolland et al. / Lithos 112 (2009) 163–187
crust-derived rocks which underwent intra-oceanic subduction and
exhumation within accretionary prisms (Rolland et al., 2009). In
contrast, the underthrusted Armenian continental crust appears not
to have been metamorphosed by any subduction event. Therefore, the
obduction of the Armenian ophiolites might be explained by the
intrinsic nature of the oceanic crust. However, the slow-spreading
nature of the ophiolites, and in particular the fact that exhumed
mantle forms a large part of the reconstructed ophiolite is rather in
agreement with a relatively dense oceanic lithosphere.
In this paper, we report new geochemical data, including major and
trace elements and Nd, Sr, Pb isotopes, on magmatic series from several
Armenian ophiolites (i.e. Stepanavan (NW Armenia), Sevan (N
Armenia), Vedi (central Armenia); Fig. 2). We identify three superposed
levels of lavas corresponding to three distinct environments: (1) backarc, (2) ‘OIB’-like and (3) arc. Moreover, we suggest that these ophiolite
165
windows correlate with each other and be part of a unique obducted
nappe. Tectonic transport of this nappe onto the SAB can be dated to the
Coniacian–Santonian (88–83 Ma; Sokolov, 1977; Sosson et al., in press).
Finally, the influence of oceanic plateau event in oceanic lithosphere
rheology and its role in the obduction process is discussed.
2. Geological setting
During the Mesozoic, the Southern Margin of the Eurasian
continent has been featured by closure of the Palaeo-Tethys and
opening of the Neo-Tethys Ocean (e. g.; Sengör and Yilmaz, 1981;
Tirrul et al., 1983; Ricou et al., 1985; Dercourt et al., 1986; Stampfli and
Borel, 2002, Fig. 1). Later on, subductions, obductions, micro-plate
accretions, ranging mostly from the Cretaceous to the Eocene, and
finally continent–continent collision have occurred between Eurasia
Fig. 2. Sketch geological map of Armenia, with location of the studied area: 1 — Stepanavan area; 2 — Sevan area; 3 — Vedi area.
166
Y. Rolland et al. / Lithos 112 (2009) 163–187
and Arabia. The study of Armenian ophiolites allows us to unravel part
of this complex history. The ophiolites are located in the northern part
of the Lesser Caucasus region (Fig. 1). The Lesser Caucasus lies south of
the Great Caucasus range, between the Black and Caspian seas. Here,
the ophiolitic belt separates the SAB from the active Eurasian margin.
The SAB is correlated westwards to the Taurides-Anatolide blocks,
which were separated from Gondwana in the Early Mesozoic (LowerMiddle Jurassic) and accreted to the Eurasian margin in the Late
Mesozoic periods (Upper Cretaceous). The Gondwanian nature of the
SAB is shown by the Proterozoic age of the basement crust and the age
and lithologies of the overlying sedimentary series (Aghamalyan,
2004; Sosson et al., in press). The active Eurasian margin is formed by
a thick volcanic arc sequence formed above an active margin, resting
on a Paleozoic (Caledonian to Hercynian) crystalline basement
(Adamia et al., 1981). The ophiolites are situated in three geographic
zones (Fig. 2):
(1) The Stepanavan ophiolite situated in NW Armenia.
(2) The Sevan ophiolite located in NNE Armenia.
(3) The Vedi zone, disposed in a more southerly position, in the
centre of Armenia.
The two first ophiolites are interpreted to correlate with each other
along the Sevan–Akera suture zone at the northern rim of the SAB, and
at the southern edge of the European active continental margin
(Knipper, 1975; Adamia et al., 1980). The Vedi ophiolite is diversely
interpreted as being an obducted sequence above the SAB (Knipper
and Sokolov, 1977; Zakariadze et al., 1983), or within a suture zone
correlating with Central Iran and Alborz ophiolites (Sokolov, 1977;
Adamia et al., 1981). It is generally believed that the different ophiolite
locations may represent suture zones, and thus feature several paleosubduction zones (Aslanyan and Satian, 1977; Knipper and Khain,
1980; Adamia et al., 1981; Aslanyan and Satian, 1982).
A companion paper written on the geology of the Sevan ophiolite
has already put up in details the lithologies and radiometric age of this
ophiolite (Galoyan et al., 2009). Main features are summarized below;
these include:
(i) A high level of fractional crystallisation in the series, with
cumulate olivine gabbros and two pyroxene gabbros overlain and
intruded by amphibole-bearing gabbros and more differentiated
melts (diorites to plagiogranites). These melts are maximally
differentiated and are generally emplaced in ductile extensive
shear zones cross-cutting the gabbros. This complete differentiation series suggests small percent partial melts and long-lived
cooling as is proposed to occur in slow spreading ophiolite settings
(Lagabrielle et al., 1984; Lagabrielle and Cannat, 1990). Absolute
radiometric datings indicate oceanic crust emplacement in the
Middle Jurassic, constrained at 165–160 Ma by zircon U–Pb age of
one tonalite (160 ± 4 Ma; Zakariadze et al., 1990) and by 40Ar/39Ar
amphibole age on gabbro (165.3 ± 1.7 Ma; Galoyan et al., 2009).
(ii) Rare pillow lavas are found, with compositions ranging from
tholeiitic basalts to andesites. The density of the feeding dyke
swarms is reduced, as rare dolerite dykes have been found
crosscutting the series. The slight calc-alkaline composition is also
evidenced by Nb–Ta negative anomalies, which agree with some
slab-derived contamination. These geochemical features support
slow spreading in a back-arc setting.
(iii) Peridotites are frequent and often exhumed as a result of intraoceanic extension. They are generally serpentinized, and witness
further hydrothermal alteration when exhumed at the contact with
marine water (‘listwenites’). The mineralogical nature of the mantlederived ultramafic rocks is still difficult to assess. The previous
petrographical investigations on the serpentinized ultramafics
suggest that protoliths were mantle-derived with various composi-
tions ranging from lherzolites to harzburgites and dunites (e.g.,
Melikyan et al.,1967; Harutyunyan,1967; Palandjyan,1971; Abovyan,
1981; Ghazaryan, 1987; Zakariadze et al., 1990). Undeformed
ultramafics have intrusive cross-cutting relationships and bear a
cumulative magmatic origin, shown by poikilitic texture of olivine
inclusions within large enstatite crystals (up to 10–15 mm;
Palandjyan, 1971). We have observed similar textures, together
with cumulative strata, contained in magmatic pods cross-cutting
serpentinites in the Stepanavan area (Galoyan et al., 2007). These
latter serpentinites are strongly deformed and altered. The ductile
character of deformation is in agreement with a mantle origin for
these rocks.
(iv) Radiolarites are found interlayered or disconformably overlying
the various lithologies described earlier. The fact that they overlie
gabbros, plagiogranites and serpentinites shows that these rocks were
uplifted and denuded by normal faults. Radiolarite datings undertaken in the different ophiolites all agree with oceanic accretion in the
Middle–Upper Jurassic (Danelian et al., 2007, 2008).
The ophiolitic sequences are weakly deformed with anchizonal
metamorphism. Only some outcrops show evidence of small shear
zones ascribed to the ophiolite obduction in the Coniacian–Santonian
(Sokolov, 1977; Zakariadze et al., 1983). HP metamorphism is
described in the Stepanavan region (Fig. 2), where blueschists outcrop
in small km2-size tectonic windows below the ophiolite. Timing of
metamorphism from radiometric 40Ar/39Ar phengite datings indicates
a HP metamorphic peak at ca 95 Ma, and MP–MT retrogression at 73–
71 Ma (Rolland et al., 2009).
The ophiolite series are locally overlain by (1) alkaline lavas, which
have a Lower Cretaceous age, though with very large error bars,
ranging from 120 to 95 ± 20 Ma, (Baghdasaryan et al., 1988; Satian and
Sarkisyan, 2006); and (2) Upper Cretaceous andesites and detrital
series (Dali valley; Stepanavan; Galoyan et al., 2007). The alkaline lavas
are alternatively interpreted as (1) intra-continental rifting (Satian
et al., 2005) in the Vedi area, and (2) plume-derived Ocean Island
magmatism above the oceanic crust before the obduction (Galoyan
et al., 2007, 2009). The calc-alkaline series are ascribed to intra-oceanic
arc emplacement above this oceanic crust sequence and implies the
presence of a subduction zone between the ophiolite and the SAB,
featured by the Stepanavan blueschists (Rolland et al., 2009). These
two magmatic sequences closely predate the ophiolite obduction onto
the SAB during the Coniacian–Santonian (Sokolov, 1977).
3. Analytical methods
Mineral compositions were determined by electron probe microanalysis (EMP). The analyses are presented in Figs. 4 and 6. They were
carried out using a Cameca Camebax SX100 electron microprobe at
15 kV and 1 nA beam current, at the Blaise Pascal University
(Clermont-Ferrand, France). Natural samples were used as standards.
For 40Ar/39Ar dating of the alkaline suite, fresh amphibole grains
were separated from the Vedi ophiolite unaltered sample AR-05-104.
Geochronology of amphiboles was performed by laser 40Ar/39Ar
dating. Results are presented in Table 1 and Fig. 7. Amphibole crystals
were separated under a binocular microscope. The samples were then
irradiated in the nuclear reactor at McMaster University in Hamilton
(Canada), in position 5c, along with Hb3gr hornblende neutron
fluence monitor, for which an age of 1072 Ma is adopted (Turner et al.,
1971). The total neutron flux density during irradiation was 9.0 × 1018
neutron cm− 2. The estimated error bar on the corresponding 40Ar⁎/
39
ArK ratio is ±0.2% (1σ) in the volume where the samples were set.
Three amphibole grains (~ 500 μm in diameter) were chosen for
analysis by the laser UV spectrometer in Géosciences Azur laboratory
at the Nice University. Analyses were done by step heating with a
Y. Rolland et al. / Lithos 112 (2009) 163–187
167
Table 1
Summary of amphibole 40Ar/39Ar dating results from the trachybasalt samples AR-05-104 and AR-05-70.
Step
Laser power (mW)
Atmospheric cont (%)
39
Ar (%)
37
ArCa/39ArK
40
Ar⁎/39ArK
Age (Ma ± 1σ)
Amphibole AR-05-104, J = 44.55, plateau age: 117.3 ± 0.9 Ma (92.4% 39Ar), isochron age: 117.5 ± 0.8 Ma (MSWD: 0.78)
1
400
99.99
0.26
8.84
2
500
95.85
0.88
9.90
3
550
92.07
0.46
2.24
4
650
24.82
6.00
4.15
5
718
8.79
6.40
5.02
6
750
4.17
22.09
5.46
7
800
0.00
9.24
5.54
8
1111
0.81
54.68
5.75
–
4.18
1.19
3.48
3.65
3.64
3.78
3.70
–
132.1
38.5
110.6
115.8
115.6
119.7
117.5
±
±
±
±
±
±
±
±
–
41.2
27.3
3.5
1.9
1.3
1.6
0.5
Amphibole AR-05-70 (1), J = 4.57, plateau age: –, isochron age: 107.8 ± 18 Ma (MSWD: 0.66)
1
340
99.99
4.44
2
460
96.12
25.09
3
560
87.34
52.41
4
640
70.78
12.68
5
640
58.83
2.15
6
1111
62.59
3.23
1.44
3.92
3.08
7.08
46.80
8.09
–
2.41
6.94
11.05
29.02
38.33
–
15
42
67
171
223
±
±
±
±
±
±
–
17
7
16
87
55
Amphibole AR-05-70 (2), J = 35.20, plateau age: –, isochron age: 114.5 ± 37 Ma (MSWD: 15)
1
380
144.64
6.23
2
500
113.84
17.22
3
772
92.47
16.80
4
1093
50.17
53.82
5
1190
123.34
0.71
6
4000
86.24
5.22
2.35
1.31
1.40
4.52
194.72
30.61
–
–
0.90
4.18
–
16.57
–
–
37
166
–
584
±
±
±
±
±
±
–
–
36
10
–
187
50 W CO2 Synrad 48–5 continuous laser beam. Measurement of
isotopic ratios was done with a VG3600 mass spectrometer, equipped
with a Daly detector system; see detailed procedures in Jourdan et al.
(2004). The typical blank values for extraction and purification of the
laser system are in the range 4.2–8.75, 1.2–3.9, and 2–6 cm3 STP for
masses 40, 39 and 36, respectively. The mass-discrimination was
monitored by analyzing air pipette volume at regular intervals. Decay
constants are those of Steiger and Jäger (1977). Uncertainties in
apparent ages in Table 1 are given at the 1σ level and do not include
the error on the 40Ar⁎/39Ark ratio of the monitor.
Thirty-seven samples of magmatic rocks from the Sevan, Stepanavan and Vedi ophiolites have been analyzed for major and trace
elements including Rare Earth Elements (REE; Table 2). Samples were
analyzed at the C.R.P.G. (Nancy, France). Analytical procedures and
analyses of standards can be found on the following website (http://
www.crpg.cnrs-nancy.fr/SARM).
For isotope measurements, powdered samples were weighed to
obtain approximately 100 to 200 ng of Sr, Nd and Pb. A leaching step with
6N HCl during 30 min at 65 °C was done before acid digestion. After
leaching, residues have been rinsed three times in purified milli-Q H2O.
Sr, Nd and Pb blanks for the total procedures were less than 50 pg, 15 pg
and 30 pg, respectively. Lead isotopes were measured by multi-collector
inductively-coupled plasma mass spectrometry (MC-ICP-MS; VG Plasma
54) at the Ecole Normale Supérieure in Lyon (ENSL). Details about isotope
chemical separations and analytical measurements including reproducibility, accuracy and standards, can be found in Bosch et al. (2008) and
references therein. The Nd and Sr isotopic data were measured on a
Finnigan MAT261 multicollector mass spectrometer at the Geochemical
Laboratory, Paul Sabatier University of Toulouse. 87Sr/86Sr was normalised to 8.3752, NBS standard was measured to 0.710250 (±15). 143Nd/
144
Nd ratio was normalised to a value of 146Nd/144Nd of 0.71219; measure
of Rennes standard was 0.511965 (±12).
4. Results
4.1. Field relationships
Synthetic logs are drawn on Fig. 3, showing the lithological
associations and the structural relationships in each of the three
studied ophiolitic zones.
4.1.1. Description of the ophiolitic units
In Stepanavan (Fig. 3A, B), ophiolite sections exhibit abundant
serpentinites, cross-cut by normal fault and shear zones in which
gabbro-norites, gabbros and plagiogranites are intrusive and
deformed (see Galoyan et al., 2007, 2009 for details and crosssections). Laterally, thick layers of pillow basalts are observed which
interlayer and are covered by radiolarites. On top of the ophiolite
section, a thin layer of alkaline lava flows is found. Above, these lavas
are eroded and unconformably covered by Upper Cretaceous
conglomerates and limestones, and calc-alkaline pillow basalts or
graywackes. The ophiolite sequence is thrusted over a blueschist facies
metamorphic sole.
In the Sevan area, sections are extremely variable laterally (Fig. 3C–E;
see Galoyan et al., 2009 for details and cross-sections). Pillow lavas are
rare, and serpentinites were frequently exhumed. Intense hydrothermal
alteration (‘listwenites’) has transformed the uppermost part of
exhumed serpentinites. The feeding doleritic dyke swarm is extremely
scarce. Large intrusive pods of amphibole-bearing gabbros and
plagiogranites are also exhumed and covered by radiolarites. Normal
faults are observed, and are interpreted as the cause of such lateral
variations, by vertical uplifting of footwall sections, and local infilling of
axial rift valleys, following the model of slow-spreading ophiolite (e.g.,
Lagabrielle et al., 1984; Lagabrielle and Cannat, 1990). Locally, thick
sequences of alkaline pillow lavas are observed. The ophiolite is locally
eroded, and unconformably covered by basal conglomerates and soils,
and an Upper Cretaceous section of limestones comprising graywackes
interlayers.
In the Vedi area, the ophiolite section is much thinner (Fig. 3F–H,
see Galoyan et al., 2009 for details and cross-sections). The basal
tectonic contact is exposed, with the top oriented to the south sense of
shear. At the base, the ophiolite rests on a serpentinite layer by a
tectonic contact. The ophiolite is intensely sheared above the basal
contact with boudins of tholeiitic basalts (Fig. 3H). Laterally, the
ophiolite consists mainly of gabbros (Fig. 3G) or serpentinites, which
suggests a similar lithology as in the Sevan and Stepanavan areas.
However, the different parts of the ophiolite are dismembered and
displaced from each other as a result of obduction deformation. Above
the ophiolites, layers of radiolarites are found below a very thick
section of alkaline pillow lavas (Fig. 3H). This section is of variable
thickness depending on the location. This may be due to lateral
168
Table 2
Representative whole-rock analyses of samples from the Sevan, Stepanavan and Vedi areas, major oxides are in wt.%, and trace elements and REE in ppm.
Groups
Sevan ophiolitic series
No.
Flaser
gabbro
Olivine
gabbro
Sevan Alkaline series
Olivine
gabbro
Gabbro
Gabbronorite
Hornblende
gabbro
Diorite
Diorite
Plagiogranite
Diabase
Trachyandesite
Basaltic
tracyandesite
Stepanavan ophiolitic series
Andesite
Basanite
Trachybasalt
Trachyandesite
Websterite
Hornblende
gabbro
Hornblende
gabbro
AR-03-25
AR-05-86
G150
AR-03-39
AR-03-24
AR-03-10
AR-04-218
AR-03-23
AR-03-19
AR-03-02
G154
AR-03-17
AR-03-34
G142
AR-05-80
AR-03-33
AR-04-03
AR-04-16
AR-04-45D
SiO2
Al2O3
Fe2O3
MnO
MgO
CaO
Na2O
K2O
TiO2
P2O5
LOI
Total
Mg#
Rb
Sr
Y
Zr
Nb
Ba
Hf
Ta
Pb
Th
U
V
Cr
Co
Ni
Cu
Zn
La
Ce
Pr
Nd
Sm
Eu
Gd
Tb
Dy
Ho
Er
Tm
Yb
Lu
Eu/Eu⁎
(La/Sm)N
(La/Yb)N
45.36
13.32
14.91
0.14
7.44
10.89
3.15
0.22
2.55
0.31
1.46
99.8
52.1
2.67
231.7
61.81
175.9
4.51
36.77
4.25
0.35
–
0.08
0.11
440.7
35.43
26.81
77.98
18.16
58.57
14.89
34.75
5.057
24.09
7.38
3.11
9.32
1.60
10.47
2.22
6.37
0.94
6.28
0.95
1.15
1.27
1.60
48.09
16.72
5.94
0.11
10.51
14.1
1.65
–
0.27
0.02
2.8
100.2
79.3
0.48
102
6.41
5.28
0.08
4.1
0.22
0.01
–
–
–
138
802
40.2
188
102
28.3
0.33
0.94
0.18
1.23
0.58
0.34
0.87
0.16
1.14
0.24
0.65
0.10
0.65
0.09
1.45
0.36
0.34
48.39
15.63
6.44
0.12
10.48
16.65
1.16
–
0.29
0.04
0.65
99.8
77.9
–
101
7.24
4.85
–
–
0.21
–
–
–
–
190
412
41.7
131
111
29.6
0.36
1.02
0.19
1.28
0.61
0.32
0.98
0.18
1.24
0.26
0.74
0.11
0.72
0.11
1.28
0.37
0.34
49.49
14.11
11.59
0.18
6.79
9.38
3.52
0.29
1.32
0.14
2.92
99.7
56.1
3.05
189.5
30.05
75.20
1.55
120.8
2.12
0.11
–
0.29
0.09
319.9
94.50
34.77
32.35
52.59
86.18
3.39
9.57
1.57
8.27
2.93
1.07
4.01
0.71
4.89
1.03
3.02
0.46
3.03
0.48
0.95
0.73
0.76
50.60
7.20
7.77
0.17
15.29
17.65
0.48
–
0.20
0.06
0.79
100.0
81.1
–
58.2
7.84
5.24
–
14.14
0.2
–
–
–
–
195.9
810.3
41.88
136.2
142.3
47.98
0.28
0.99
0.22
1.40
0.67
0.28
1.07
0.20
1.36
0.29
0.83
0.13
0.86
0.13
1.01
0.26
0.22
50.68
18.17
9.09
0.16
6.75
9.26
3.21
0.15
0.36
0.05
1.21
99.1
61.6
0.84
303.7
11.93
22.89
0.50
34.36
0.81
0.04
–
0.14
0.07
222.2
104.3
33.71
29.23
44.81
71.99
1.60
4.51
0.74
3.80
1.34
0.48
1.70
0.30
1.92
0.41
1.22
0.19
1.32
0.21
0.97
0.75
0.82
55.09
13.45
8.51
0.15
9.84
8.76
2.59
0.17
0.24
0.04
1.87
100.7
71.6
1.42
207.4
7.82
21.08
0.32
34.9
0.76
0.03
–
0.07
0.04
158.1
562.7
39.6
148.1
15.63
73.35
1.44
3.73
0.57
2.87
0.94
0.35
1.13
0.20
1.31
0.28
0.82
0.13
0.92
0.15
1.03
0.97
1.06
57.41
14.10
8.84
0.14
2.24
4.92
6.36
0.12
0.87
0.16
4.33
99.5
35.6
1.81
212.5
25.22
75.07
1.83
55.59
2.27
0.14
2.62
0.75
0.25
122.6
136.9
15.76
10.09
21.02
50.45
4.62
11.73
1.82
9.05
2.91
0.97
3.68
0.69
4.60
0.97
2.95
0.47
3.25
0.52
0.91
1.00
0.96
74.91
12.32
3.58
0.03
0.42
3.05
4.20
0.31
0.21
0.04
0.57
99.6
20.2
2.22
145.2
27.65
71.64
2.15
65.83
2.48
0.10
1.32
1.09
0.62
50.8
1464
5.52
37.19
6.91
9.98
5.38
12.66
1.78
8.36
2.68
0.67
3.50
0.64
4.33
0.94
2.94
0.46
3.30
0.53
0.67
1.26
1.10
46.02
16.29
8.39
0.13
7.73
10.68
3.53
0.37
1.26
0.17
4.61
99.2
66.6
13.97
630.9
23.93
127.4
2.95
285.3
2.91
0.24
1.51
1.01
0.29
179.3
277.0
38.26
54.43
58.20
66.89
7.07
18.32
2.71
12.68
3.53
1.36
3.95
0.67
4.21
0.84
2.45
0.36
2.39
0.37
1.12
1.26
2.00
53.70
14.09
11.35
0.15
4.52
3.64
6.07
–
1.36
0.12
4.85
99.9
46.2
–
28.76
27.82
81.43
1.77
16.09
2.36
0.14
1.61
0.55
0.35
334.9
–
29.48
8.95
102.8
67.20
4.25
11.14
1.67
8.76
2.93
1.14
4.04
0.71
4.78
1.01
2.99
0.46
3.16
0.51
1.01
0.91
0.91
54.27
15.16
12.36
0.19
3.74
4.38
6.63
–
1.33
0.15
1.78
100
39.5
–
49.82
29.88
85.60
1.44
19.01
2.38
0.12
–
0.37
0.27
321.8
–
24.33
5.14
15.61
80.0
3.91
10.93
1.81
9.35
3.26
1.15
4.18
0.75
4.99
1.06
3.11
0.47
3.22
0.51
0.95
0.75
0.82
55.48
14.13
12.45
0.13
4.07
5.49
3.96
0.61
1.17
0.11
2.18
99.8
41.6
5.05
102.5
27.49
54.39
1.69
16.65
1.63
0.12
–
0.33
0.14
347.4
251.7
27.09
16.09
5.23
19.68
2.69
7.17
1.15
6.17
2.32
0.79
3.36
0.62
4.34
0.95
2.82
0.44
2.98
0.47
0.86
0.73
0.61
40.63
14.40
11.70
0.29
4.15
11.40
3.64
1.42
2.06
0.48
10.01
100.2
43.4
31.89
147.0
25.58
153.4
40.63
166.7
3.57
2.99
4.98
4.13
1.43
257.4
33.86
38.26
34.29
61.86
100.1
32.41
64.63
7.64
29.89
6.02
1.97
5.58
0.82
4.69
0.90
2.48
0.35
2.33
0.36
1.04
3.39
9.39
43.80
17.58
9.46
0.11
6.70
4.95
4.09
2.24
1.68
0.44
8.89
99.9
60.5
44.83
330.8
17.37
131.4
17.95
299.0
2.86
1.19
4.63
4.06
0.94
271.7
21.68
31.44
28.76
54.48
91.55
29.46
59.03
6.92
27.11
5.21
1.63
4.45
0.62
3.42
0.61
1.64
0.22
1.44
0.22
1.03
3.56
13.83
51.57
14.34
6.06
0.11
0.77
9.78
6.34
0.56
1.98
1.08
6.65
99.3
21.7
7.93
341.5
52.91
411.2
49.22
168.7
9.03
3.80
4.45
5.15
4.61
286.1
162.0
24.06
13.83
18.34
92.74
48.54
107.1
13.38
56.55
12.98
4.14
12.43
1.84
10.27
1.87
4.80
0.64
4.05
0.60
1.0
2.35
8.09
53.24
1.03
6.02
0.15
23.18
16.52
0.12
–
0.05
0.03
0.54
100.9
90.0
–
11.79
1.05
–
–
3.71
–
–
–
–
–
135.1
2804
55.82
361.3
340.3
25.68
–
0.15
0.02
0.15
0.08
0.03
0.13
0.03
0.18
0.04
0.12
0.02
0.12
0.02
0.85
0.0
0.0
47.30
14.39
12.90
0.21
9.11
10.14
2.93
0.19
1.18
0.07
1.76
100.2
60.4
1.12
125.4
20.91
42.74
1.01
32.71
1.21
0.08
–
0.18
0.05
324.7
236.4
51.46
78.0
–
60.28
2.40
6.37
1.08
5.76
2.09
0.96
2.98
0.53
3.53
0.74
2.15
0.32
2.10
0.32
1.17
0.72
0.77
53.77
14.00
8.92
0.15
7.81
6.98
3.34
2.42
0.16
0.05
2.44
100.1
65.6
30.17
213.4
5.79
21.26
2.14
228.1
0.63
0.21
3.61
1.27
0.43
94.47
324.2
31.5
101.7
189.9
60.33
3.06
6.28
0.65
2.30
0.49
0.19
0.55
0.10
0.81
0.19
0.64
0.12
0.91
0.17
1.12
3.94
2.27
Y. Rolland et al. / Lithos 112 (2009) 163–187
Sample
Table 2 (continued )
Stepanavan ophiolitic series
Stepanavan alkaline series
Stepanavan calk-alkaline
series
Vedi ophiolitic series
Vedi Alkaline series
No.
Plagiogranite Basaltic trachy- Basaltic trachy- Basaltic trachy- Basaltic
Diabase
andesite
andesite
andesite
trachy-andesite
Olivine
basalt
Basaltic
trachyandesite
Basaltic Hornblende
trachygabrro
andesite
Diorite
Plagiogranite Basalt
Basaltic
andesite
Basalt
Trachybasalt
Basaltic
Trachydacite
trachyandesite
Sample
AR-04-44
AR-04-20
AR-04-30
AR-06-02
AR-03-53
AR-0405
AR-0432
AR-0440A
AR-0431
AR-05-113
AR-05110
AR-05-111
AR-05114
AR-05106
AR-0578
AR-05-104
AR-05-102
AR-04-75
SiO2
Al2O3
Fe2O3
MnO
MgO
CaO
Na2O
K2O
TiO2
P2O5
LOI
Total
Mg#
Rb
Sr
Y
Zr
Nb
Ba
Hf
Ta
Pb
Th
U
V
Cr
Co
Ni
Cu
Zn
La
Ce
Pr
Nd
Sm
Eu
Gd
Tb
Dy
Ho
Er
Tm
Yb
Lu
Eu/Eu⁎
(La/Sm)N
(La/Yb)N
75.35
12.20
2.71
0.03
0.77
2.05
5.03
–
0.11
0.02
1.07
99.3
38.1
0.58
91.04
1.23
6.48
0.35
20.58
0.15
–
–
0.02
0.01
30.74
421.5
7.71
24.86
189.5
23.13
2.42
3.90
0.42
1.58
0.29
0.35
0.24
0.04
0.22
0.05
0.15
0.03
0.21
0.04
4.05
5.25
7.92
51.53
14.69
14.81
0.23
4.15
4.86
5.74
0.18
1.62
0.13
1.89
99.8
37.7
1.4
61.01
35.91
86.0
1.62
19.58
2.54
0.13
2.29
0.43
0.12
459.8
99.23
38.85
22.76
64.58
130.6
4.23
11.08
1.88
9.88
3.44
1.28
4.69
0.87
5.93
1.28
3.79
0.58
3.89
0.61
0.98
0.77
0.73
48.55
13.29
8.67
0.15
6.86
10.49
4.74
0.24
1.08
0.11
6.01
100.2
63.1
7.94
95.7
26.52
68.43
1.9
21.73
1.83
0.15
1.93
0.19
0.09
305.4
316.7
42.98
109.1
132.8
81.13
2.53
7.37
1.31
7.04
2.55
0.99
3.49
0.65
4.34
0.93
2.75
0.42
2.84
0.44
1.02
0.63
0.60
45.37
14.27
13.52
0.32
6.22
4.09
1.53
5.37
3.12
1.31
5.11
100.2
49.8
45.3
157
40.6
268
52.8
608
5.97
3.97
7.04
5.39
1.34
172
–
26.3
6.5
20.4
177
50.3
96.1
12.1
53.3
11.2
3.87
10.5
1.50
8.11
1.49
3.75
0.51
3.18
0.45
1.09
2.83
10.68
48.54
15.01
12.65
0.27
4.25
5.33
3.93
2.69
2.64
1.08
3.16
99.5
42.3
33.17
322.6
44.41
294.4
57.95
578.3
6.51
4.20
2.54
5.98
1.46
94.35
25.77
16.75
–
9.48
137.1
50.59
107.0
12.97
53.32
11.26
4.08
10.41
1.52
8.69
1.57
4.16
0.57
3.64
0.56
1.15
2.83
9.38
50.19
13.91
13.73
0.24
3.27
5.85
5.11
0.42
3.39
0.67
2.94
99.7
34.0
7.62
198.8
51.24
373.5
42.33
156.6
8.01
3.24
2.24
4.65
1.20
201.6
–
31.93
–
14.94
152.7
40.02
85.12
10.85
45.27
10.35
3.39
10.3
1.63
9.50
1.80
4.93
0.70
4.51
0.69
1.0
2.43
5.99
49.15
18.53
10.19
0.16
5.25
8.25
4.36
0.52
0.86
0.14
3.12
100.5
52.7
9.61
520.3
16.0
44.4
2.14
133.9
1.25
0.17
7.22
0.72
0.19
241.5
21.05
29.51
15.55
12.6
150.9
4.93
11.46
1.74
8.35
2.36
0.94
2.60
0.44
2.80
0.57
1.62
0.24
1.60
0.25
1.15
1.31
2.07
49.79
15.80
8.82
0.15
3.54
9.12
3.54
1.24
1.07
0.20
7.27
100.6
48.4
18.12
303.8
24.4
99.18
2.29
239.1
2.69
0.18
3.42
1.46
0.67
279.1
31.91
29.05
22.38
188.1
86.53
7.87
18.51
2.68
12.58
3.47
1.13
3.89
0.65
4.09
0.85
2.46
0.38
2.54
0.39
0.94
1.43
2.09
52.20
17.05
9.28
0.16
3.59
6.56
4.61
1.00
0.94
0.18
5.36
100.9
44.5
18.45
282.3
24.19
95.3
3.32
213.6
2.61
0.26
5.66
1.67
0.58
263.7
73.83
27.36
18.95
170.9
100.0
8.69
18.05
2.53
11.65
3.15
1.04
3.57
0.61
3.96
0.83
2.46
0.38
2.54
0.40
0.95
1.74
2.31
45.09
21.24
4.32
0.07
7.88
14.29
2.12
0.17
0.16
–
5.08
100.4
79.7
2.41
492.7
3.70
4.30
0.08
131.6
0.16
–
–
–
–
89.5
785.2
30.98
130.7
92.4
24.06
0.231
0.71
0.14
0.78
0.36
0.24
0.53
0.10
0.64
0.14
0.38
0.06
0.35
0.06
1.64
0.40
0.44
58.57
16.03
6.66
0.11
5.18
7.08
3.94
0.47
0.33
0.04
2.08
100.5
62.7
4.12
254
10.1
42.3
0.49
57.5
1.3
0.04
1.1
0.42
0.12
185
123
22.8
42.4
18.3
50.7
2.09
4.84
0.72
3.74
1.20
0.42
1.46
0.25
1.65
0.34
1.01
0.16
1.14
0.18
0.97
1.10
1.24
70.45
14.69
4.44
0.08
1.06
3.87
4.47
0.2
0.43
0.09
1.04
100.8
34
1.11
161
13.1
86.5
0.59
33.5
2.35
0.05
1.12
0.48
0.13
44.5
9.1
6.6
6.4
–
39
2.91
6.15
0.81
4.48
1.45
0.69
1.81
0.32
2.11
0.46
1.35
0.21
1.52
0.25
1.31
1.26
1.29
47.5
16.17
8.76
0.15
8.46
8.57
3.99
0.72
0.93
0.09
4.78
100.1
67.6
3.51
134.6
20.78
54.9
0.69
141.7
1.48
0.07
–
0.15
0.06
211.3
416.7
42.44
195.5
14.05
65.52
1.96
6.25
1.09
5.9
2.15
0.90
3.0
0.53
3.48
0.74
2.08
0.32
2.15
0.33
1.08
0.57
0.62
48.3
15.03
10.14
0.16
5.97
8.01
3.96
0.18
1.2
0.12
7.14
100.2
56.0
4.14
110.1
26.91
73.98
2.46
16.2
2.04
0.20
–
0.22
0.10
238.3
324.7
49.08
122.6
81.78
94.28
3.07
8.51
1.47
7.97
2.88
1.12
4.02
0.70
4.67
0.97
2.78
0.42
2.83
0.45
1.0
0.67
0.73
44.58
12.52
9.36
0.12
2.63
15.53
3.83
–
2.35
0.33
9.17
100.4
37.8
0.63
153.5
24.94
160.8
23.22
1097
3.91
1.76
1.53
2.085
0.643
240.7
50.2
28.9
26.28
28.41
106.8
18.44
39.11
5.02
21.85
5.50
1.93
5.72
0.85
4.85
0.88
2.30
0.31
1.95
0.29
1.05
2.11
6.39
44.64
15.41
11.99
0.14
4.85
7.85
4.24
0.96
3.67
0.85
5.04
99.6
46.6
10.48
926
36.26
318.8
67.52
444.9
6.81
4.88
1.25
4.60
1.18
219.8
4.22
31.4
21.81
50.82
146.8
49.35
109.7
14.15
59.18
12.52
4.15
10.96
1.48
7.84
1.32
3.26
0.42
2.49
0.35
1.08
2.48
13.39
50.39
16.2
7.76
0.13
5.05
7.16
3.24
1.96
2.36
0.64
5.07
99.9
58.4
27.1
643
22.4
260
43.3
659
6.03
3.29
6.71
8.39
1.8
135
136
64.3
126
42.2
134
50.7
91.8
9.83
42.9
8.47
2.67
7.36
0.99
5.01
0.81
1.92
0.25
1.46
0.22
1.03
3.77
23.44
59.61
17.48
7.64
0.12
1.11
1.89
6.37
2.4
0.72
0.25
2.14
99.7
23.16
64.89
260.4
58.18
680.5
82.24
422
14.76
5.99
4.30
12.28
2.78
5.26
66.9
4.98
–
9.93
167.9
74.8
142.6
15.84
59.79
12.64
3.68
11.69
1.877
10.97
2.08
5.77
0.86
5.84
0.88
0.93
3.73
8.64
Y. Rolland et al. / Lithos 112 (2009) 163–187
Groups
169
170
Y. Rolland et al. / Lithos 112 (2009) 163–187
Fig. 3. Representative geological logs of the Stepanavan, Sevan and Vedi ophiolites.
Y. Rolland et al. / Lithos 112 (2009) 163–187
variations in the amount of erupted lavas, or may be explained by
tectonic scalping of the ophiolite upper part during obduction. These
alkaline lava flows are well exposed and preserved in the Vedi ophiolite.
They are made of amphibole-bearing basaltic pillows. The pillows are
larger (metre scale) than the ophiolites ones (several decimetre scale),
and interlayer with thin pink limestones. At the front of the obduction in
the Vedi zone, an olistolith formation exhibits conglomerates and slided
blocks in a muddy matrix (Fig. 3F). The olistostrom age is Coniacian–
Santonian (nanofossils, Carla Muller, com. pers.), and it connects
progressively below and above with Lower and Upper Coniacian reef
limestones, respectively. Therefore, the obduction age can be bracketed
to the Coniacian–Santonian, which agrees with former estimates
(Sokolov, 1977). Laterally, always in the Vedi zone, the upper part of
the ophiolite is made of kilometre scale slided blocks, mainly comprised
of alkaline pillow basalts and calc-schists. These blocks slide on a
greenish mudstone rock, probably originated from the ophiolite
alteration. The Upper Coniacian uncomformity is variably marked by
conglomerates, marls and reef limestones.
4.1.2. General features of the Armenian ophiolites: evidence for LOT
features
As emphasized in Galoyan et al. (2009) in the Sevan area, and by
Galoyan et al. (2007) in the Stepanavan area, the lithologies found in
all the exposed Armenian ophiolites are in good agreement with the
hypothesis of a slowly expanding spreading centre, as described for
the western Alps ophiolites and exposed earlier in Section 2 (Nicolas
and Jackson, 1972; Nicolas, 1989).
The similar lithological and age features found in the several
Armenian ophiolites suggest that they were part of the same oceanic
crust section. This has to be confirmed by the comparison of
geochemical data from each zone. The presence of three magmatic
series: ophiolite s.s. (tholeiitic), ‘OIB’ (alkaline) and arc (calc-alkaline)
in the same structural position (from bottom to top, respectively) has
been evidenced in the three zones. Isotopic geochemistry on the three
series will allow us to constrain the nature of sources and to identify
the magmatic processes that existed prior to ophiolite obduction.
4.2. Petrography and mineral chemistry
The field and microscopic analyses of Armenian ophiolite magmatic rocks show a continuous magmatic succession from ultramafic
cumulates (wherlites, websterites) to gabbros and plagiogranites,
cross-cutting intensely altered serpentinites in each of the three
studied ophiolites. All these lithologies were exhumed in the footwall
below normal faults and covered by pillow-basalts.
4.2.1. Serpentinites
The study of serpentinite mineralogy is difficult due to intense
serpentinization. However, EMP analysis of chromiferous spinel relicts
from Sevan Tsapatagh area (sample AR-05-80 in Galoyan et al., 2009;
Fig. 3) reflects still unaltered mineral compositions. EMP analyses
show a very narrow compositional range in Cr# (Cr / Cr + Al = 0.71–
0.73) and Mg# (Mg / Mg + Fe = 0.58–0.59) ratios (Fig. 4). These Cr#
compositions are more elevated than those of abyssal peridotites
(Brynzia and Wood, 1990) and agree with a fore-arc peridotites
composition (Parkinson and Pearce, 1998). However, Mg# values are
slightly lower than those of Parkinson and Pearce (1998), which is
ascribed to high partial melting in such context. However, we cannot
exclude any hydrothermal process, as the primary nature of chromites
is uncertain.
4.2.2. Ophiolite plutonic rocks
Wehrlites are found in the Stepanavan ophiolite (see Galoyan et al.,
2007). They have a poikilitic texture showing numerous clinopyroxene crystals with diopside composition (Wo45–47En48–50Fs2–4),
included in large olivine Fo87–88 (N60–65%) porphyric grains.
171
Fig. 4. Chemistry of Cr-spinel from Armenia with respect to the compositional fields of
Abyssal peridotites (1; Brynzia and Wood, 1990) and arc-related peridotites (Mariana
seamount peridotites (2) and dunites (3); Parkinson and Pearce, 1998).
Gabbros are the most abundant rocks in the crustal complex, and
are found in each ophiolite zone (Galoyan et al., 2007, 2009). Their
petrography evolves from cumulate-banded olivine gabbros in their
lower part towards more leucocratic plagioclase-rich gabbros in the
upper part (Abovyan, 1981; Ghazaryan, 1987, 1994). The cumulative
banded olivine gabbros and websterites are found locally, only in the
Sevan and Stepanavan areas, while more leucocratic gabbros are
widespread in the three zones.
Olivine gabbros found in Stepanavan and Sevan ophiolites (Galoyan et
al., 2007, 2009) are fresh, massive, and fine- to medium-grained (0.5 to
2 mm). They have cumulate, ophitic textures and consist of plagioclase
(~60–65%; An68–74, An80–89), olivine (~5–10%; Fo72–76), and clinopyroxene (~25–35%). Clinopyroxene is of augite (Wo39–44En45–48Fs11–13)
and diopside (Wo45En44Fs11) types. Some enstatite orthopyroxenes
(Wo2En75Fs23) are also found rimming olivine porphyrocrysts.
Websterites found in Stepanavan and Sevan ophiolites (Galoyan
et al., 2007, 2009) have a granular texture with large 2–8 mm
porphyrocrysts of orthopyroxene (30–70%), clinopyroxene (70–30%)
and olivine grains (0–35%; Fig. 5A). Orthopyroxenes are enstatite-rich
(Wo1–5En59–84Fs11–37) and clinopyroxenes are augites (Wo35–42En36–
40Fs15–19), olivine is relatively rich in forsterite (Fo84–88). Gabbronorites (from Stepanavan; Galoyan et al., 2007) have a gabbroic texture,
with plagioclase (10–60%, 1–3 mm), clino- and ortho-pyroxene.
Plagioclase is of bytownite type (An80–85), while orthopyroxenes (1–
4 mm) are enstatites (Wo2–5En59–61Fs34–37), and clinopyroxenes are
augites (Wo35–42En36–40Fs15–19).
Mesocratic to leucocratic gabbros of the upper section found in the
three ophiolites (Galoyan et al., 2007, 2009; Rolland et al., in press) are
massive, fine- to medium-grained and have gabbroic (or gabbro-ophitic),
xenomorphic granular texture (0.5–4 mm), with plagioclase (~40–65%;
An50–75, An72–93), clinopyroxene (8–45%; augite) and hornblende (0–
40%), without any olivine. Accessory minerals (1–10%) are apatite,
titanomagnetite, ilmenite and rarely quartz. The hornblende-rich gabbros
(Galoyan et al., 2009; Rolland et al., in press) have coarse granular
textures (Fig. 5B), with ~50–65% euhedral to subhedral plagioclase
(An54–58) and (~35–50%) anhedral to subhedral amphibole. Some brown
Ti-rich euhedral hornblende is presumed to be a primary mineral; while a
Ti-poor subhedral to xenomorphic green magnesio-hornblende (Leake et
al., 1997) is thought to be a secondary phase formed by hydrothermal alteration as it replaces generally the brown type. The augite
(Wo40–42En39–47Fs11–14), diopside (Wo45–48En40–44Fs8–15), and enstatite
(Wo2En57Fs41) relicts (5–10%) are found in the crystals of magnesiohornblendes that replace the pyroxenes. However, it is not related to
shear zones and fractures, and is thus a late magmatic mineral. In
leucocratic gabbros (Galoyan et al., 2009), the clinopyroxene (augite
172
Y. Rolland et al. / Lithos 112 (2009) 163–187
Wo40–41En33–35Fs18–19) content does not exceed 25%. Normal zoning is
observed in plagioclase (from An85 to An60), which is frequently
altered. Clinopyroxenes have alkaline to slightly tholeiitic compositions
(0.8 b Na+ Cab 0.9; Leterrier et al., 1982, Fig. 6). Pegmatitic gabbros
crosscutting the plutonic sequence (Vedi zone; Rolland et al., accepted)
are composed of plagioclase and hornblende, and are mainly altered
into chlorite, carbonate, sericite, albite, quartz, actinolite, etc.
Diorites occur as small intrusive bodies within the gabbro units in
Sevan and Vedi zones (Palandjyan, 1971; Abovyan, 1981; Ghazaryan,
1987, 1994). They have a porphyritic to subhedral granular (1–4 mm)
Y. Rolland et al. / Lithos 112 (2009) 163–187
173
Fig. 6. Chemical compositions of studied clinopyroxenes plotted in the Ti vs. (Na + Ca) diagram of Leterrier et al. (1982). Note that a majority of data plot in the Alkaline compositional
field, and a minority is in the Toleiitic part.
texture and have relatively similar hornblende contents (5–30%) as
gabbros. Plagioclase (~65–70%) is albite-rich (An34–38) and accessory
minerals (quartz, opaque oxides) are rare. Amphibole grains are
magnesio-hornblendes in composition, sometimes rimmed by actinolite fringes and epidote aggregates. In Sevan zone, diorites grade
into quartz-diorites (quartz 5–10%), laterally and upwards in the
series.
Plagiogranites are found in the three zones (Galoyan et al., 2007,
2009; Rolland et al., in press). They appear to be dioritic intrusives
most differentiated components, forming diffuse segregations or
discontinuous networks of veins. Plagiogranites have local coarse
pegmatic, or hypidiomorphic to xenomorph granular (0.5–4 mm)
textures. They are formed by 40–65% plagioclase (An15–30), 25–45%
quartz, minor biotite (b5%), ortho-amphibole (b5%; Stepanavan), Kfeldspar (0–10%, microcline; Fig. 5C) and accessory phases (titanomagnetite, hematite, sphene and apatite).
4.2.3. Ophiolite volcanic and subvolcanic rocks
Diabases are present in several locations (Sevan and Stepanavan
areas; Galoyan et al., 2007, 2009) as isolated dikes, crosscutting the
layered gabbros. They are generally altered (chlorite, epidote, carbonates) and have a subdoleritic texture composed of plagioclase (60–65%;
An65–75) and two clinopyroxenes (augite Wo41–44En44–47Fs11–13 and
diopside Wo45En37Fs18).
The volcanic rocks of the Armenian ophiolites we studied are
present as pillowed and massive lava flows and pillowed breccias. In
general, they show signs of hydrothermal alteration but relict igneous
textures are preserved. In the three locations basalts and basaltic
andesites are vesicular (1–5 mm, filled with carbonate-calcite, chlorite
and quartz) and largely aphyric (intersertal, spilitic, microdoleritic
and variolitic, up to 1.5–2 mm in diameter), composed mainly of
albitized plagioclase and/or plagioclase–clinopyroxene microlites, Timagnetite and hematite microlites, in a devitrified (calcite + chlorite)
groundmass (Fig. 5D).
4.2.4. Alkaline lavas
The alkaline basalts are found in the three zones on top of the
ophiolite section as large massive pillow-lavas or as diabase dykes, but
their relationships with the ophiolite (s.s.) pillow lavas remain unclear.
The first group of alkaline rocks displays large vesicles (0.5–3 mm), filled
with carbonates and rarely chlorites, and have both phyric and aphyric
(Fig. 5E). They have intersertal textures, with plagioclase megacrysts
(~5%; 0.5–2 mm), microliths and opaque minerals (3–10%), surrounded
by a calcite–chlorite mesostase. The second group (Vedi and Stepanavan
zones, e.g., samples AR-05-70 and AR-05-104 dated by 40Ar/39Ar) have
doleritic (Fig. 5F) to ophitic textures. They are mainly composed of
plagioclase (~40–55%; 1–3 mm), clinopyroxene (10–30%; 1–4 mm),
amphibole (~25%; 1–3 mm) and accessory Ti-magnetite (N5–10%),
apatite (~3%; prismatic, acicular, 0.5–1.5 mm) and rarely biotite. Apatites
are present in the plagioclase crystals and in the vitreous interstices,
which are filled by carbonates or carbonates-chlorites. The tabular
plagioclase laths show a transitional zoning with bytownite to labrador
(An72–60) or labrador to andesine (An55–32) compositions. Thin rims of
pure albite (Ab — 98%) are also present. Clinopyroxenes are generally
chloritized, but still preserve diopside compositions (Wo49En35Fs16).
The amphibole is a kaersutite (Leake et al., 1997), with zoning from
kaersutite to ferro-kaersutite from core to rim. Some samples show
abundant calcite-filled veins and pockets.
A few dacitic sills and dikes occur among the basaltic pillow lava
flows in the Vedi valley. As in the pillow basalts, plagioclase is the
main mineral phase and Fe-oxides are present (~5–10%; Fig. 5G).
Some 1–2 mm large unzoned plagioclase phenocrysts of oligoclase-
Fig. 5. Microphotograph of representative magmatic rock types from the Armenian ophiolite complex. Plutonic and volcanic ophiolite series: (A) subautomorph granular texture of a cumulate
banded websterite (sample AR-04-36, Stepanavan area, Cheqnagh valley; see Galoyan et al., 2007); (B) coarse-grained hornblende gabbro with normally zoned plagioclases (sample AR-05110, Vedi area, massif of Qarakert, see Rolland et al., in press-b); (C) xenomorph granular texture of a microcline (Mc) bearing plagioclase rich leucogranite (sample AR-05-109, in the same
massif, see Rolland et al., accepted); (D) aphyric, intersertal (spilitic) and variolitic basalt composed of mainly albitized plagioclase, Ti-magnetite and hematite microlites, in a devitrified
groundmass (sample AR-05-106, Vedi area, Khosrov valley, see Rolland et al., accepted). Alkaline series: (E) aphyric, intersertal basalt, totally devoid of phenocrysts, and composed of
carbonatized plagioclase microlites and opaque minerals (~5%) in a chlorite–carbonate groundmass (sample AR-05-80, Sevan area, Tsapatagh valley, see Galoyan et al., 2009); (F) doleritic
texture in a trachybasalt composed of plagioclase, chloritized clinopyroxene, kaersutite (Krs), Ti-magnetite and apatite (sample AR-05-104, Vedi area, see Rolland et al., accepted); (G) phyric
trachydacite with a hyalopilitic to cryptocrystalline texture (sample AR-04-75, Vedi valley, see Rolland et al., accepted). Calc-alkaline series: olivine-bearing, plagioclase phyric (15–40%) basalt
with a microcrystalline (plagioclase, quartz, opaque minerals) texture from pillow lavas suit (sample AR-04-32, Stepanavan area, Herher valley, see Galoyan et al., 2007), in which the olivine
phenocrysts are entirely pseudomorphosed to quartz and rims of iron oxides. From (A) to (C) under crossed nichols, and (D) to (H) under parallel nichols. Scale bar is for all photographs.
174
Y. Rolland et al. / Lithos 112 (2009) 163–187
andesine compositions are distributed in the fine-grained devitrified
groundmass made of albitic plagioclase, opaque microlites, and
carbonate-quartz-chlorite aggregates.
4.2.5. Calc-alkaline lavas of Stepanavan zone
They consist of large pillow-lavas of basaltic and basaltic andesitic
compositions with micro-cryptocrystalline (Fig. 5H) to intersertal textures
formed of large phenocrysts (2–7 mm) and microliths of andesineoligoclase plagioclase, and minor augite (Wo36–38En42–43Fs13–15) clinopyroxenes. These lavas overlie Upper Cretaceous limestones, unconformably lying on the ophiolite stricto sensu.
4.3.
40
Ar/39Ar dating
Complementary to previously published datings (Galoyan et al.,
2009) obtained from the ophiolite, which span the Middle Jurassic, we
provide here the first unambiguous dating of the alkaline suite. In the
Stepanavan and Sevan regions, the alkaline lavas have been deformed
and altered in the late collisional evolution, so that preserved and
unaltered amphibole-bearing lavas were only sampled in the Vedi
area.
Three analyses have been done on amphibole single grains from
two trachybasalt samples (AR-05-70 and AR-05-104) from the Vedi
ophiolite, which is described in Section 4.2 though only one is
considered fully successful. These datings are listed in Table 1 and the
successful one is presented in Fig. 7.
In the two datings for sample AR-05-70 (Table 1), the 39Ar content
was very low, so no plateau age can be calculated. However, from the
36
Ar/40Ar versus 39Ar/40Ar plots, it was possible to estimate isochron
ages, with large errors, of 108 ± 18 Ma (MSWD: 0.66) and 115 ± 37 Ma
(MSWD: 15). The very low 39Ar content is interpreted as a
consequence of the low K content of amphibole, which likely resulted
from pyroxene destabilisation in this sample.
In the dating of sample AR-05-104, a well-constrained plateau of
117.3 ± 0.9 Ma (2σ) was obtained, with 92% of released 39Ar (Fig. 7A).
The average 37ArCa/39ArK ratio is similar as the EPM value of the
amphibole from ~40 in low temperature steps, decreasing steadily to
~ 30 in high temperature steps (Fig. 7B). An isochron age of 117.5 ±
0.8 Ma (MSWD: 0.77) is obtained using the five steps of the plateau
age estimate (steps 3–7), with an initial 40Ar/36Ar ratio close to the
atmospheric value [(40Ar/36Ar)0 = 238 ± 4%; Fig. 7C]. Including the
step 1 of lower temperature we calculate a similar within-error
isochron age of 117.5 ± 0.8 Ma (2 σ).
The above age of 117.5 ± 0.8 Ma is more precise than the previous
ages determined by the whole-rock K–Ar method (Baghdasaryan
et al., 1988) which ranged between 114 and 97 Ma. We interpret
these younger K–Ar ages as resulting from alteration of the vitreous
matrix. These ~ 118 Ma Albian ages are also in agreement with age
estimates undertaken by Satian and Sarkisyan (2006) who provided
whole-rock Rb/Sr errorchron ages between 120 and 95 Ma. The
Alkaline sequence is thus undoubtedly younger than the ophiolite by
about 50 Ma, unlike some generally admitted views that all the
alkaline, calc-alkaline and tholeiitic series part of the obducted
sequence were formed in the same mid-oceanic context (Sokolov
1977; Knipper and Khain, 1980).
4.4. Major-trace-REE geochemistry
The geochemical analyses of the ophiolitic rocks from the Sevan
ophiolite are of relatively alkaline composition in comparison to
MORB. Major element data of pillow — lavas and plutonic rocks show
that they have predominantly basalt to trachybasalt compositions.
4.4.1. Major elements
Major element analysis of plutonic rocks ranges from gabbros to
granites (plagiogranites) with intermediate dioritic compositions
Fig. 7. 40Ar/39Ar amphibole dating results of trachybasalt sample AR-05-104 from the
Vedi alkaline suite.
(Fig. 8A). These magmatic rocks plot in a large domain comprised
between alkaline and tholeiitic tendencies of the TAS diagram (Le
Maitre et al., 1989). In the AFM diagram (Fig. 8B) most rocks lie close
to the limit between the tholeiitic and calc-alkaline fields.
1. Overall, the rocks of the ophiolitic suite are enriched in MgO and
more depleted in TiO2, K2O and P2O5 relative to the alkaline suite
(Figs. 8–10; Table 1). Compared to the plutonic rocks of the same
series, the volcanic rocks from the different areas plot in the same
compositional range (from basalts to andesites and trachyandesites) and are slightly Na2O richer.
2. The alkaline lavas from different zones plot in the same range,
varying compositionally from basanite-trachybasalt to basaltic
Y. Rolland et al. / Lithos 112 (2009) 163–187
175
Fig. 8. Plots of magmatic rocks (ophiolitic, alkaline and calc-alkaline series) in the (A) (Na2O + K2O) vs. SiO2 (Le Maitre et al., 1989) and (B) AFM (Irvine and Baragar, 1971) diagrams.
trachyandesite and trachyandesite, and are clearly in the calcalkaline/alkaline domain of the AFM diagram (Fig. 8A, B). One of
the most significant features of the alkaline lavas is their higher
TiO2, K2O and P2O5 contents.
3 The arc-type calk-alkaline lavas, have trachybasalt and basaltic
trachyandesite compositions in TAS diagram (Fig. 8A). They occupy
a transitional position between ophiolitic and alkaline domains in
Harker's diagram (Fig. 9), except lower TiO2 and higher Al2O3,
which depend on the abundance of plagioclase in such rocks.
Regarding the spread of compositional variations in major
elements within the series, it appears that only rough correlations
can be seen in the plots of SiO2 vs. other oxides (Fig. 9). Even the most
immobile elements during alteration processes, such as Al2O3, MgO
and TiO2 (e.g., Staudigel et al., 1996) do not show any clear correlations
(Fig. 9). In particular, Large Ion Lithophile Elements (LILE) such as Na
and K have scattered compositions, even in individual magmatic
suites. Such variations are ascribed to a combination of alteration and
magmatic processes. The occurrence of a long-lasting hydrothermal
event, ascribed to the slow-spreading oceanic environment is
indicated by scattered 40Ar/39Ar ages within individual gabbro
samples (Galoyan et al., 2009; Rolland et al., in press). Thus, we
ascribe the most important part of the elemental variability in the
ophiolitic suite (gabbros, diorites, plagiogranites and ophiolitic lavas)
to spilitization process in an oceanic environment while variations in
the alkaline and calc-alkaline volcanic rocks is ascribed to some
magmatic cause (variations in the source components, as is highlighted by the isotopic compositions, Section 4.5).
Thin section observations and previous studies of the Armenian
ophiolites (e.g., Palandjyan, 1971; Abovyan, 1981; Ghazaryan, 1994)
have shown that the whole ophiolitic sequence apart from the
alkaline lavas has been affected by oceanic low-temperature
hydrothermal alteration events. These processes induced modification of the whole-rock chemistry, as revealed by the increase of LOI
(Table 1).
4.4.2. Trace elements
High field strength elements (HFSE) are not mobilized during
alteration and their contents reflect, without ambiguity, those of their
parental magma (Staudigel et al., 1996). Contents in these trace
elements confirm the presence of three clearly distinct magmatic
suites, as defined in the previous section.
1. Basalts and gabbros of the ophiolite suite show strong enrichments
in LILE (Large Ion Lithophile Elements: Ba, Rb, K and Th), up to ten
times MORB values. They have negative anomalies in Nb–Ta and Ti
(Fig. 11A, B), which is generally indicative of volcanic island arc
environments (e.g. Taylor and McLennan 1985; Plank and Langmuir, 1998). However, trace element contents remain low relative
to volcanic arc lavas, which indicates a setting with little fractional
crystallization, in agreement with a back-arc setting (e.g., Galoyan
et al., 2009).
2. Overall, the concentrations of each element in the alkaline basalts
largely exceed the concentrations in the basalts from ophiolitic
series (Fig. 11C). Moreover, alkaline series basalts are characterized
by high abundances of LILE, high field strength elements (Nb, Ta, Zr
and Ti) and light rare-earth elements (LREE).
3 The calc-alkaline suite rocks show strong depletions in Nb and Ta,
relative to Th and La, and slight Ti negative anomaly (Fig. 11D). They
globally show slightly stronger enrichments in LREE and LILE than
the ophiolite suite rocks.
These differences in normalized element patterns support that
these basalts are not petrogenetically related and were most likely
derived from melts formed in different tectonic settings: (1) A backarc setting with slow-spreading rates, (2) Ocean-island within-plate
setting and (3) volcanic island arc.
Differences within the different suites can be related to magmatic
processes such as fractional crystallization and magma mixing or to
alteration processes. To analyse the importance of alteration, some
trace and two major elements are plotted versus Zr and Th (Fig. 10). Zr
is an incompatible element (for basaltic to andesitic lavas) wellknown to remain stable during alteration or weathering processes, so
it was used as reference to test the mobility of the other trace elements
(Fig. 10). The trace element composition of both ophiolitic series and
arc-type lavas plots in a restricted range of values for Zr (0–127 ppm),
while alkaline lavas with higher Zr are characterized by a large
compositional range (131–411 ppm; Table 2). This compositional
spread is ascribed to various levels of fractional crystallization with
one trachydacite sample having very high Zr content (681 ppm).
The concentration of zirconium normally increases in response to
magmatic processes such as fractional crystallisation, except for the
most differentiated lavas in which it has fractionated. Enrichment in Zr
is positively correlated with that of major elements as Ti, but there is
no clear correlation with SiO2, which are ascribed to a slight
176
Y. Rolland et al. / Lithos 112 (2009) 163–187
Fig. 9. Harker variation diagrams showing the compositions of the three (ophiolitic, alkaline and calc-alkaline) series.
fractionation of Si within each series. Compositional variations in SiO2
and TiO2 in alkaline lavas correlate well with variations in Zr (Fig. 10),
while traces and REE such as Nd and Th contents show slight positive
correlation with Zr contents, which are ascribed to fractionation of
these elements. Such process is also shown by the composition of
ultramafic cumulative plutonic rocks (websterite from Stepanavan
and gabbronorite from Sevan), which have lower trace element
contents than associated lavas due to their cumulative origin. In
contrast, Sr does not correlate well with Zr, and part of scattering of Sr
data may be due to its mobility due to alteration processes of
plagioclase. This is confirmed by a similar mobility of Ca, as shown in
the Ti vs. Ca diagram (Fig. 9).
4.4.3. REE geochemistry
In the chondrite-normalized rare earth element (REE) diagrams,
analysed ophiolite basalts and gabbros have flat and parallel REE
spectra in chondrite-normalized plots [(La/Yb)N = 0.6–0.9], showing
some slight depletions in LREE and a slight enrichment in MREE
(Fig. 11E, F). No extensive Eu anomalies were observed (Eu/
Eu⁎ = 0.95–1.15), which show that plagioclase has remained almost
stagnant, and is enriched in the final liquid. The concentration of REE
for volcanic rocks varies from 8 to 30 times chondrite and for gabbros
varies between 1 and 15 times chondrite, for exception a flaser gabbro
— 60 times (sample AR-03-25). These features are interpreted as a
result of extreme crystal fractionation involving plagioclase, clinopyroxene, orthopyroxene and, to a lesser extent, to olivine accumulation
(Pallister and Knight, 1981).
The websterite and gabbronorite have the lowest concentrations of
REE (0.1–0.9 and 1–5 times chondrite respectively) with patterns
characterized by depletion in LREE (Fig. 11F). One hornblende gabbro
(sample AR-04-45D from Stepanavan ophiolite) is characterized by
LREE enrichment ((La/Yb)N = 2.27) and some depletion in MREE (a
convex downward pattern) with smaller positive Eu anomalies (Eu/
Eu⁎ = 1.12).
Y. Rolland et al. / Lithos 112 (2009) 163–187
The diorites REE patterns (6–20 times chondrite) and plagiogranites are parallel to those of the gabbros, with smaller enrichment in
LREE ((La/Yb)N = 1.1). The most differentiated plagiogranite (sample
AR-04-44 from Stepanavan) is characterized by more depletion in the
middle to heavy REE compared to other plagiogranites, and strongly
positive Eu anomalies (Eu/Eu⁎ = 4.05) ascribed to high plagioclase
contents. Such features indicate a cumulative effect of plagioclase.
In contrast, chondrite-normalized REE patterns of alkaline lavas
(Fig. 11G) show huge LREE enrichments and HREE depletions [(La/
Yb)N = 6–14], being representative of intraplate continental basalts, as
compared to ophiolite lavas. Meanwhile, no extensive Eu anomalies are
observed (Eu/Eu⁎ = 0.95–1.15). One trachydacite sample (AR-04-75) is
featured by significant enrichments in trace elements, which is
explained by a high degree of fractional crystallization, as its REE
pattern is parallel to those of the basanite–trachyandesite series.
177
Chondrite-normalized REE patterns of calc-alkaline lavas are
strongly parallel and form a narrow domain (Fig. 11H). They have
similar HREE contents as volcanics of previous series with significantly
more depleted LREE contents than alkaline series rocks [(La/Yb)N =
2.1–2.3].
These differences in trace elements contents between the three
studied series further support that these basalts are petrogenetically
unrelated and, most likely derived from melts formed in different
tectonic settings.
4.5. Nd, Sr, Pb isotope geochemistry
4.5.1. Ophiolite series
Initial ɛNdi values of the ophiolitic lavas from the different studied
zones range from + 5.9 to +9.5 (Table 3) intermediate between
Fig. 10. Major and trace elements vs. Zr and Th diagrams. Major and trace elements are chosen to investigate the effects of alteration on the Sr, Nd and Pb isotopic systems;
explanations in the text.
178
Y. Rolland et al. / Lithos 112 (2009) 163–187
The initial Pb isotopic ratios in ophiolitic rocks range from 37.273 to
38.234 for 208Pb/204Pb, from 15.459 to 15.541 for 207Pb/204Pb and
from 17.630 to 18.459 for 206Pb/204Pb. In the Pb–Pb isotope diagrams
both ophiolitic volcanic and plutonic rocks plot on or close the MORB
domain (Fig. 12B).
4.5.2. Alkaline series
Overall, alkaline lavas show, in comparison to ophiolite rocks,
lower initial ɛNdi values ranging from +2.1 to + 4.0 but have a similar
range of (87Sr/86Sr)i ratios. In (143Nd/144Nd)i vs. (87Sr/86Sr)i diagram
(Fig. 12A) most of these rocks are located in the OIB field. In the
(208Pb/204Pb)i vs. (206Pb/204Pb)i, and (207Pb/204Pb)i vs. (206Pb/204Pb)i
diagrams these samples overlap various specific OIB provinces such
as Kerguelen, Samoa and Society and Marquises (Fig. 12B).
4.5.3. Calc-alkaline series
For the calc-alkaline lavas from Stepanavan, the initial ɛNdi values
and the (87Sr/86Sr)i ratios range from +3.8 to + 5.9 and from 0.7045
to 0.7048, respectively. Thus ɛNdi ratios appear to be intermediate
between those of ophiolitic and alkaline rocks, and Sr isotopic ratios
plot in the same range as these two series (Fig. 12A).
5. Discussion
Fig. 10 (continued).
typical MORBs and OIBs values and indicate a source region that
experienced long-term depletion in LREE (Fig. 11E). The initial (87Sr/
86
Sr)i ratios range from 0.7037 to 0.70565 for these lavas, which are
significantly too high for typical tholeiitic MORB lavas and thus, not
considered as the primary magmatic Sr isotope signature of these
rocks. The gabbros have ɛNdi (+4.3 to + 6.5) and initial Sr isotopic
values (0.70386 to 0.70557) in the same range as volcanic rocks.
In the Nd–Sr isotope diagram both ophiolitic volcanic and
plutonic rocks exhibit a significant increase of the initial Sr ratios
relatively to MORB (Fig. 12A). Moreover, this increase of (87Sr/86Sr) i
ratios positively correlates with Sr, Ba, Rb and K2O contents. This
shift towards 87Sr/86Sr radiogenic ratios is commonly attributed to
exchange between rocks and seawater during oceanic crust
hydrothermal alteration (e.g. McCulloch et al., 1981; Kawahata
et al., 2001).
Ophiolites of the Armenian Lesser Caucasus region are generally
separated into three distinct zones: (1) The Sevan-Akera zone in the
North (Knipper, 1975; Adamia et al., 1980), (2) The Zangezur zone in
the centre (Aslanyan and Satian, 1977; Knipper and Khain, 1980;
Adamia et al., 1981) and (3) The Vedi zone in the south (Knipper and
Sokolov, 1977; Zakariadze et al., 1983).
Due to the importance of Cenozoic volcanism spread over most of
Armenia (Fig. 2), it is still difficult to conclude only from geological
mapping whether the different ophiolites correlate with each other,
or if they represent various suture zones delimitating several
continental micro-blocks. For this reason, we have undertaken field
investigations in various ophiolites: Stepanavan (Galoyan et al.,
2007) and Sevan (Galoyan et al., 2009), along the northern rim of
Armenia; and the Vedi ophiolite, in the centre of Armenia. As
emphasized in the following discussion, the use of Nd, Sr and Pb
isotopes in complement to conventional major and trace element
data allow us to correlate the ophiolites with each other. These
ophiolites show some similarities and differences in their structure
and lithological successions, but these features remain compatible
with a single oceanic domain origin. This domain opened in the
Lower-Middle Jurassic and has undergone several phases of magmatic emplacement, for which we find evidence in each of the
geographic zones investigated. These correlations provide insight
into the evolution of the Tethyan domain, and in particular allow us
to propose a geodynamical model for the obduction of the ophiolite
over the Armenian block.
In the following discussion, we will evaluate the following points:
1. Petrographically and geochemically, the Armenian ophiolites are
similar to island-arc tholeiites. Such geochemical features are
typical for oceanic crust, formed in a back-arc setting with melting
of a shallow asthenospheric source contaminated by slab-derived
fluids (Saunders and Tarney, 1984). Such a hypothesis has already
been proposed for ophiolitic gabbros from Turkey (Kocak et al.,
2005), but has to be evaluated considering isotopic compositions
and partial melting constraints.
2. Alkaline lavas of variable thickness have covered this ophiolitic
sequence. Their origin has to be considered. (i) Do they also derive
from the same ophiolitic series? (ii) Did they form in an island-arc
setting or (iii) in an oceanic island/plateau environment? The
source of alkaline lavas will be discussed below regarding the Sr–
Nd–Pb isotopic data. The occurrence of alkaline magmatism prior
Y. Rolland et al. / Lithos 112 (2009) 163–187
179
Fig. 11. Trace and REE plots of the three studied magmatic suites. The multi-element spider diagrams are normalized to the N-MORB values of Sun and McDonough (1989), and REE
plots are normalized to the Chondrite values of Evensen et al. (1978). Patterns for the studied magmatic rocks: ophiolitic volcanic (A, E) and plutonic (B, F) series; OIB type alkaline
series (C, G), and arc type calk-alkaline series (D, H).
180
Table 3
Sr, Nd and Pb isotopic analyses of samples from ophiolitic complexes of Armenia.
Sevan
(206Pb/204Pb)
2σ
(238U/204Pb)
(206Pb/204Pb)i
(207Pb/204Pb)
2σ
(235U/204Pb)
(207Pb/204Pb)i
(208Pb/204Pb)
2σ
(232Th/204Pb)
(208Pb/204Pb)i
143
Nd/144Nd
2σ
(147Sm/144Nd)
(143Nd/144Nd)i
ɛNdi
(87Sr/86Sr)
2σ
(87Rb/86Sr)
(87Sr/86Sr)i
ɛSri
Stepanavan
165
165
AR-03-25
18.2577
0.00049
11.1328
17.9683
15.5033
0.00055
0.0818
15.4888
37.9444
0.0018
8.8411
37.8713
0.512957
0.000011
0.18522
0.51276
6.47
0.704433
0.000011
0.1326
0.70412
− 2.61
165
Vedi
165
165
165
117
117
117
165
117
117
95
95
165
165
165
117
117
117
AR-03-24 AR-03-10
AR-04-218
AR-03-02
G154
G142
AR-05-80
AR-04-44
AR-04-30
AR-03-53
AR-04-05
AR-04-32
AR-04-32
AR-05-113
AR-05-114
AR-05-106
AR-05-104
AR-05-78
AR-04-75
17.9321
0.00084
10.2463
17.6658
15.5344
0.00084
0.0753
15.5211
37.8028
0.0025
8.4704
37.7328
0.512957
0
0.29064
0.51264
4.25
0.70401
0.00001
0.0058
0.70399
−4.39
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
n.d.
0.512911
0.000044
0.21337
0.51268
4.98
0.703875
0.000053
0.008
0.70386
− 6.39
0.512848
0.000039
1.80575
0.5109
− 29.80
0.703134
0.000048
1.2024
0.70031
− 56.69
0.512929
0.000039
0.16855
0.51275
6.27
0.704589
0.000048
0.064
0.70444
1.88
0.51275
3.4E− 05
0.11069
0.51263
4.04
0.70471
4.6E− 05
0.0184
0.70466
5.12
0.51315
0.00003
0.21883
0.51291
9.53
0.706076
0.000039
0.3765
0.70519
12.59
0.51282
3.7E− 05
0.17125
0.51271
3.8
0.70485
4.8E− 05
0.0534
0.70478
5.54
0.5129
4E−05
0.1668
0.5128
5.87
0.7048
4E−05
0.1725
0.7045
1.91
17.8966
0.00057
10.274
17.6296
15.5524
0.00065
0.0755
15.5391
37.7996
0.0024
63.7004
37.2729
0.512989
0.000042
0.22032
0.51275
6.35
0.705745
0.000013
0.0754
0.70557
17.92
18.0686
0.00065
5.7363
17.9195
15.4661
0.00064
0.0421
15.4587
37.9304
0.0017
63.5101
37.4053
0.512965
0.000007
0.22032
0.51273
5.88
0.70509
0.00001
0.0754
0.70491
8.62
18.4101
0.00051
10.6564
18.1073
15.497
0.00049
0.0783
15.4897
38.1116
0.0015
64.1466
37.5546
0.513007
0.00006
0.21848
0.51277
6.74
0.70591
0.00001
0.1087
0.70565
19.15
20.6925
0.00076
62.4439
19.6744
15.6653
0.0005
0.4587
15.6158
40.6747
0.0016
252.2401
39.3623
0.512709
0.00008
0.1279
0.51262
2.3
0.70453
0.000009
0.0327
0.70448
1.47
19.2867
0.00047
26.8011
18.8497
15.5651
0.00055
0.1968
15.5439
39.2226
0.0021
89.804
38.7553
0.512813
0.000011
0.15218
0.51271
4.01
0.704272
0.00001
0.6009
0.70338
−14.11
n.d.
n.d.
19.2791
0.00051
34.3464
18.719
15.5926
0.00035
0.2523
15.5654
39.6069
0.0013
138.024
38.8887
0.51275
6E− 06
0.13822
0.51265
2.86
0.70568
9E− 06
0.1109
0.70551
16.15
n.d.
n.d.
19.2349
0.00068
36.8567
18.6339
15.5901
0.00061
0.2707
15.5609
39.5224
0.002
155.9984
38.7107
0.512696
0.00001
0.12767
0.51261
2.05
0.705421
0.000011
0.2974
0.70498
8.57
n.d.
n.d.
18.61
0.00092
12.7305
18.4024
15.5803
0.00123
0.0935
15.5702
38.6686
0.0034
56.9259
38.3724
0.512773
0.000008
0.11619
0.51269
3.7
0.704072
0.000011
0.392
0.70349
− 12.56
n.d.
n.d.
19.1457
0.00046
18.2699
18.8478
15.5681
0.00061
0.1342
15.5536
39.1236
0.002
54.4859
38.8401
0.512754
0.000008
0.12176
0.51267
3.26
0.70636
0.0001
0.6276
0.70543
14.97
n.d.
n.d.
18.8172
0.00062
13.7665
18.4594
15.559
0.00061
0.1011
15.5411
38.4167
0.0025
22.1422
38.2336
0.512966
0.000038
0.20223
0.51275
6.28
0.705061
0.000012
0.5785
0.7037
−8.55
n.d.
n.d.
0.512765
0.000037
0.12781
0.51268
3.39
0.705331
0.00005
0.7208
0.70427
− 1.59
Notes: isotopic data (2σ error) are corrected for in situ decay assuming a mean age of 165 Ma of the ophiolite (Galoyan et al., 2009), 117 Ma for the alkaline series (this paper) and 95 Ma for the calc-alkaline series from palaeontogical and Ar-Ar
dating (Rolland et al., in press).
i: initial ratios calculated at 165, 104 and 95 Ma respectively. εNdi calculated with actual (143Nd/144Nd)CHUR = 0.512638 and (147Sm/144Nd)CHUR = 0.1967 (Wasserburg et al., 1981). εSri calculated with actual (87Sr/86Sr)CHUR = 0.7045 and (87Rb/86Sr)
CHUR = 0.0827 (Wasserburg et al., 1981).
Pb isotopic ratios measured with external precision of ca. 250-300 ppm for the 206, 207, 208Pb/204Pb ratios.
Y. Rolland et al. / Lithos 112 (2009) 163–187
Locality
age (Ma)
sample
Y. Rolland et al. / Lithos 112 (2009) 163–187
Fig. 12. Plots of the magmatic rocks in (A) Sr–Nd; (B)
Northern Hemisphere Reference Line (Hart, 1984).
206
Pb/204Pb vs.
208
Pb/204Pb and
to obduction in the late Lower Cretaceous may be of importance for
the obduction model of the ophiolite crustal sequence.
3. Finally, the calk-alkaline lavas are Upper Cretaceous in age. These
volcanic arc-related series likely formed during closure of the NeoTethys ocean. Their geochemical features will be considered to
evaluate this hypothesis.
5.1. Significance of Armenian ophiolites: MOR or back-arc setting?
Armenian ophiolitic series are shown to be of slight alkaline to
tholeiitic character, ranging from basalts to basaltic andesites and basaltic
trachyandesites. Spider diagrams show clear Nb–Ta negative anomalies
(Fig. 11A, B), LILE enrichments and flat to slightly LREE-enriched spectra.
Their isotopic compositions are significantly more radiogenic in 87Sr/86Sr
and slightly less radiogenic in Nd isotopes than typical MORB compositions (Fig. 12A). These observations do not support a geochemical
“normal” ophiolitic crust and are more probably in agreement with
typical volcanic arc settings, in which enrichments in LILE, LREE result
from slab fluids/melts contaminations (Pearce et al., 1984).
5.1.1. Source components
The isotopic compositions of the ophiolitic magmatic rocks lie at
the limit of the MORB domain and overlap the OIB field (Fig. 12A, B).
Nevertheless, their flat REE spectra together with their “enriched”
isotopic character suggest partial melting from a spinel-bearing
mantle with small percent partial melts similar as to those known
for slow spreading ridges (Lagabrielle, 1987). This hypothesis will be
tested by a non-modal batch partial melting model in the following
section. Further, the measured Nb–Ta negative anomalies combined
with slight LILE enrichments are indicative of a volcanic arc setting.
Finally, the emplacement depth of pillow lava flows was clearly
abyssal as shown by the deposition of radiolarite interlayers. There-
181
207
Pb/204Pb isotopic diagrams, with the fields of MORB and of various OIB contexts. NHRL,
fore, it is likely that the ophiolitic rocks were produced by melting of a
depleted mantle source contaminated by hydrothermal slab-derived
fluids in a back-arc basin environment.
To estimate the level of slab-derived contamination in the
formation of the ophiolitic rocks, mixing curves have been drawn on
Fig. 13 between different components: depleted mantle pole (MORB),
Enriched Mantle 1 and Enriched Mantle 2 [EM1 and EM2, respectively; Zindler and Hart (1986), Salters and White (1998) and Hanan
et al. (2000)]. This isotopic modelling suggests that basaltic ophiolite
lava composition results from contamination of a typical MORB by a
mixed source composed of 1–4% EM2 and 2–5% EM1. Such degrees of
contamination appear to be relatively high in a back-arc setting, and
suggest the participation of subducted slab sediments in the source
(EM1) and a possible fertile E-MORB type source (EM2).
5.1.2. Partial melting estimates
Fig. 14 shows calculated REE composition for melts produced by a
depleted mantle type source (MORB) partial melting (Fig. 14A) and by
an enriched mantle source (EM2) partial melting (Fig. 14B). The REE
patterns of the ophiolitic suite fit those of melts produced by the partial
non-modal batch melting from 4% to 10% of a spinel-bearing mantle,
which mineralogical compositions spinel 5%, olivine 55%, orthopyroxene 20% and clinopyroxene 20%, using the partition coefficient factors
of McKenzie and O'Nions (1991), Johnson (1994), and Nikogosian and
Sobolev (1997). Such composition is rather similar to that of a slightly
depleted mantle source (e.g., Juteau and Maury, 1997).
5.2. Origin of alkaline lavas: source components and geodynamic
significance?
Mineral chemistry and geochemistry of the alkaline volcanic series
of Sevan, Stepanavan and Vedi ophiolites is similar to that of OIBs. As
182
Y. Rolland et al. / Lithos 112 (2009) 163–187
Fig. 13. Isotopic diagrams showing mixing curves between the different mantle end-members. 143Nd/144Nd vs. 206Pb/204Pb (A) and 87Sr/86Sr (B) isotopic diagrams; and 206Pb/204Pb
vs. 87Sr/86Sr (C) and 207Pb/204Pb (B). Compositions of end-members used in the calculation of mixing curves (after Hart, 1984; Zindler and Hart, 1986; Sun and McDonough, 1989;
Eisele et al., 2002) are the following. HIMU: (87Sr/86Sr) = 0,703; (143Nd/144Nd) = 0,51285; (206Pb/204Pb) = 21,5; (207Pb/204Pb) = 15,82; (208Pb/204Pb) = 40; [Sr] = 120; [Nd]A = 6,5;
[Pb] = 0,4. EM2 pole: (87Sr/86Sr) = 0,71682; (143Nd/144Nd) = 0,51216; (206Pb/204Pb) = 18,99; (207Pb/204Pb) = 15,65; (208Pb/204Pb) = 39,5; [Sr] = 218; [Nd] = 34; [Pb] = 25. EM1
pole: (87Sr/86Sr) = 0,705; (143Nd/144Nd) = 0,5122; (206Pb/204Pb) = 16,8; (207Pb/204Pb) = 15,45; [Sr] = 513; [Nd] = 33; [Pb] = 3,5; DMM: (87Sr/86Sr) = 0,7022; (143Nd/144Nd) =
0,513075; (206Pb/204Pb) = 17,3; (207Pb/204Pb) = 15,4; (208Pb/204Pb) = 37,5; [Sr] = 11,3; [Nd] = 1,12; [Pb] = 0,0489. The mixing curves equations are from Faure (1986).
Fig. 14. Modelling effect of non-modal batch melting of several mantle sources, and comparison with ophiolite and alkaline REE spectra compositional domains obtained in this study.
Calculated composition (A) of the depleted peridotite ‘ophiolite’ source: olivine 55%, orthopyroxene 20%, clinopyroxene 20% and spinel 5%; and (B) of the enriched peridotite
‘alkaline’ source: olivine 54%, orthopyroxene 20%, clinopyroxene 20%, garnet 2% and spinel 4% (Salters and Stracke, 2004).
Y. Rolland et al. / Lithos 112 (2009) 163–187
Fig. 15. Geodynamic reconstitution of the Lesser Caucasus in the Middle Jurassic to Upper Cretaceous periods.
183
184
Y. Rolland et al. / Lithos 112 (2009) 163–187
shown in the Mineral Chemistry section, pyroxenes are slightly
alkaline. The alkaline lava samples (Fig. 14) show strong enrichments in incompatible elements (up to 100 times chondrite
values). Partial melting degree calculation suggests that they may
be derived from ~ 20% non-modal batch melting of an enriched
spinel-garnet-bearing mantle source characterized by La and Lu
concentrations 3.3 times and 1.8 times the chondrite mantle,
respectively. The isotopic composition of these lavas plots in the
field of OIBs in agreement with an enriched alkaline mantle source.
In the isotopic plots of Fig. 13, the obtained isotopic data plots
mainly on the DMM-EM2 mixing curve, with ~ 2–3% of EM2, ~ 5–
15% HIMU and almost no EM1 contamination which suggests that
subduction-derived contamination may not be envisaged. The
isotopic compositions of the studied Armenian alkaline series are
thus rather in agreement with an OIB-type source. In the Vedi area,
Satian et al. (2005) already pointed out the alkaline character of
the lava series, which they interpreted as volcanic series formed in
an intra-continental rift. These lavas were emplaced above, and
formed ~ 50 Ma after the ophiolites. Moreover, they are interstratified and overlain by shallow marine reef limestone. Thus, all
these features are in agreement with a plume event that occurred
in an intra-oceanic setting.
Such alkaline magmatism is widely documented in the MiddleEast region, along the Arabian and Indian platforms, in relationship
with the formation of the Neo-Tethys ocean (e.g. Lapierre et al., 2004).
Similar Cretaceous alkaline series are found above the Iranian
ophiolite (Ghazi and Hassanipak, 1999), and in Turkey (Norman,
1984; Tüysüz et al., 1995; Tankut et al., 1998; Okay, 2000). However, it
is still difficult to relate these alkaline events due to their geographical
and temporal distance and to paucity of radio-chronological and Sr,
Nd, Pb isotopic data.
5.3. Reconstruction of the ‘ophiolite’ history
From all the available geological data, we propose the following
model for the evolution of the Armenian Ophiolite (Fig. 15):
1. The SAB is of Gondwanian origin according to lithological
associations found in central and SE Armenia (Knipper and
Khain, 1980; Kazmin et al., 1987; Aghamalyan, 2004). Therefore,
it is likely that the Sevan oceanic basin opened in response to the Ndipping subduction of Neotethys to the south of Eurasia (Fig. 15—
stage 1). The continuation of Paleotethys in the area is a matter of
debate, as the westward continuation of the Cimmeride orogenic
system is not identified in the Lesser Caucasus (e.g., Sengör, 1984,
1990). Emplacement of the ophiolite occurred in the Lower-to
Middle Jurassic (Galoyan et al., 2009). The older age of the Vedi
ophiolite (178.7 ±2.6 Ma; Rolland et al., accepted), with respect to
that of Sevan (160–165 Ma, Zakariadze et al., 1990; Galoyan et al.,
2009) implies that it was at the southern rim of the back-arc
system. The structural setting of the ophiolite obduction indicates
clearly that oceanic crust of the back-arc basin was formed between
the SAB and the Eurasian active margin.
2. Emplacement of an Oceanic Island/Plateau above the back-arc oceanic
crust during the late Lower Cretaceous (40Ar/39Ar age of 117.5±0.8 Ma
in this paper; Fig. 15—Stage 2).
3. The calc-alkaline lavas disconformably overlie the ophiolite and
related alkaline series (Galoyan et al., 2007). These lavas have
similar geochemical features as volcanic arc series, including the
isotopic Sr–Nd composition. Their emplacement is bracketed in the
Upper Cretaceous, as for the high pressure metamorphism
constrained in the Stepanavan area (Meliksetyan et al., 1984 and
references therein), constrained at about 95–90 Ma (Rolland et al.,
2009). Therefore this magmatic event can be related to the
subduction of Neo-Tethys ocean prior to the obduction of the
Armenian Ophiolites onto the SAB.
4. After this, the SAB enters the subduction zone in the Turonian (95–
88 Ma), which triggers a “collision” with the thickened oceanic
crust. During this process, part of the volcanic arc has probably
been subducted below the obducted oceanic sequence and metamorphosed in the blueschist facies (Rolland et al., 2009). The
large variety of lithologies comprising metabasites, marls and
conglomerates in a pelitic matrix, within the Stepanavan blueschists, is in agreement with such a scenario.
5. The obduction of the ‘ophiolite’ section over the SAB is further
constrained by the Lower Coniacian frontal flysch sequences,
found below and in front of the Vedi obducted sequence. The calcalkaline series found above the Stepanavan ophiolite show that
a volcanic arc was active during this time above the obducted
sequence.
6. The end of the obduction is constrained by Upper Coniacian fauna
in sediments unconformably overlying the ophiolite. Blocking of
the subduction below the Eurasian margin may stop at 73–71 Ma,
as shown by Ar–Ar age of MT-LP metamorphism in the Stepanavan
blueschists and the general tectonic uplifting of the region,
witnessed by erosion and absence of sedimentary record during
the Upper Cretaceous–Paleocene (Rolland et al., 2009). This 73–
71 Ma event is thus interpreted as the insight of ‘collision’.
5.4. Implications of a plateau/OIB event on the ophiolite obduction?
The alkaline series show features of Plume-related magmas. What
is the significance of such plume magmatism, and what is its
consequence with the obduction and preservation of the Armenian
ophiolites?
Only 3% of the current oceanic floor is composed of plumerelated crust, of which oceanic plateaus are the largest part
(Petterson et al., 1997). In Armenia, this alkaline event seems to
be relatively large due to the presence of alkaline lavas over the
ophiolite in all the studied sections. This large size is thus in
agreement with an oceanic plateau event, in which large volumes of
lavas are erupted during volcanic emplacement in a small time
range.
Only several examples of obducted plateau series have been
claimed worldwide, amongst which the Wrangellia terrane of Alaska
and British Columbia (e.g. Richards et al., 1991), and Gorgona Island in
Columbia (Duncan and Hargreaves, 1984; Storey et al., 1991), but such
examples remain relatively uncommon. The paucity of obducted
plateau sequences may be explained by the fact that they are not easily
recognized in the geological records. However, their potential in the
blocking or reversion of polarity of subduction zones has been noted
in numerous cases (e.g. Petterson et al., 1997; Kerr et al., 2003; Kerr
and Mahoney, 2007). For instance large oceanic plateaus can cause the
reversal of subduction polarity, as did the Ontong Java Plateau
(Coleman and Kroenke, 1981). Cloos (1993) calculated that basaltdominated oceanic plateau crust must exceed 17 km thickness to
survive subduction, and about 30 km to cause any significant
‘collisional’-type deformation.
The Armenian ophiolites show evidence for the obduction of a
single oceanic crust sequence above the SAB, as similar geological,
petrological, geochemical and age features are found in the three
studied Armenian ophiolitic massifs (Sevan, Stepanavan, and Vedi).
The oceanic crust s.s. corresponds to a slow-spreading ophiolite
formed in the Lower-Middle Jurassic in a back-arc basin by 4–10%
melting of a shallow asthenosphere spinel-bearing source contaminated by subducted slab-derived products. Alkaline volcanic
series with OIB-type geochemical features are found above the
ophiolite sequence in each of the studied areas, which late Lower
Cretaceous age has been constrained above by the 40Ar/39Ar
method on amphibole at 117.3 ± 0.9 Ma. Therefore, this Armenian
‘plume event’ shortly predates the Coniacian–Santonian (88–
83 Ma) obduction of the Armenian ophiolitic sequence. Therefore,
Y. Rolland et al. / Lithos 112 (2009) 163–187
the thickened and hot oceanic plate had a low density when it
overrode the SAB continental margin, which suggests that the
plume event has likely played an important role in the obduction
process.
The original width and thickness of the Armenian plume related
series are difficult to assess. However, this volcanic series covered
each of the studied Armenian ophiolites, which suggests a N104 km2
surface regarding the initial ophiolite surface prior to horizontal
shortening, and N105–106 km2 if they can be correlated to similar
settings in Turkey and Iran. Further, these lavas have significantly
more radiogenic lead isotopic compositions than ophiolitic rocks,
and are related to ~ 20% melting of an enriched garnet-spinelbearing source. These high partial melting estimates are in the range
of plateau events melting estimates (14–26%, Hauff et al., 1997;
Révillon et al., 2000; Herzberg and O'Hara, 2002; Kerr et al., 2002).
The melting of a garnet-bearing source suggests the contribution of
a deep mantle source below a thick lithosphere sequence or the
melting of a relatively thick oceanic crust (N30 km, as is constrained
by the stability of garnet in metabasites; Rapp et al., 1991). Such
thickness has been estimated for the Ontong Java plateau (e.g.,
Gladczenko et al., 1997; Richardson et al., 2000; Miura et al., 2004).
Expected mantle lithosphere “root”, below such thick oceanic crust
could be N300 km (Richardson et al., 2000). Therefore, the
Armenian hot-spot related magmatism likely features the formation
of an Oceanic Plateau or a large oceanic island, with significant
crustal thickening. Further, the age of the Armenian plume-related
series is within range of a period of major oceanic plateau formation
in the late Lower Cretaceous; fitting precisely the age of formation of
one of the largest plateaus, the Ontong Java plateau (Tarnudo et al.,
1991). These alkaline series are also locally covered by an arcderived calc-alkaline volcanic sequence, which was likely formed in
a supra-subduction zone environment. Further evidence of this
subduction is provided by blueschists series dated at 95–90 Ma
(Rolland et al., 2009).
Therefore both oceanic Plateau and volcanic arc formations shortly
pre-dated the obduction, which occurred in the Coniacian–Santonian
(88–83 Ma; Sokolov, 1977). Crustal thickening related to plateau and
arc events are thought to have increased crustal buoyancy (e.g., Cloos,
1993; Abbot and Mooney, 1995; Abbot et al., 1997; Kerr and Mahoney,
2007). Such low buoyancy likely hindered subduction of the oceanic
crust and allowed it to be obducted over the SAB continental crust.
Such process is not unlikely in other obduction contexts of the PeriTethyan region, especially in the Caucasus–Middle East segment, but
the obducted ophiolite sections have to be analysed in detail to find if
whether alkaline series of similar age and geochemical signatures may
be present.
Acknowledgements
This work was supported by the Middle East Basins Evolution
project jointly supported by a consortium including oil companies
and the CNRS. Many thanks to the MEBE program coordinators E.
Barrier and M. Gaetani for their support and encouragements, and M.
F. Brunet for coordinating the project. Analytical data were acquired
with the help of the Geosciences Azur Laboratory, in which we thank
L. Vacher and J.P. Goudour for their involvement during data
acquisition. We also thank the support of the French Embassy at
Yerevan for the MAE PhD grant of G. Galoyan. The paper was
significantly enhanced by the work of two anonymous reviewers and
the editor N. Eby, while English language was improved by G. Nolet.
References
Abbot, D., Mooney, W., 1995. The structural and geochemical evolution of the
continental crust: support for the oceanic plateau model of continental growth.
Reviews of Geophysics 33, 231–242.
185
Abbot, D., Drury, R., Mooney, W., 1997. Continents as lithological icebergs: the
importance of buoyant lithospheric roots. Earth and Planetary Science Letters
149, 15–27.
Abovyan, S.B., 1981. The mafic–ultramafic complexes of the ophiolitic zones in
Armenian SSR. Izvestia Academy of Science Armenian SSR. 306 pp. (in Russian).
Adamia, S., Bergougnan, H., Fourquin, C., Haghipour, A., Lordkipanidze, M., Ozgül, N.,
Ricou, L., Zakariadze, G., 1980. The Alpine Middle East between the Aegean and the
Oman traverses. 26th International Geological Congress Paris C5, 122–136.
Adamia, S., Bergougnan, H., Fourquin, C., Haghipour, A., Lordkipanidze, M., Ozgül, N.,
Adamia, S., Chkhotua, T., Kekelia, M., Lordkipanidze, Shavishili, I., Zakariadze, G.,
1981. Tectonics of the Caucasus and the adjoining regions: implications for the
evolution of the Tethys ocean. Journal of Structural Geology 3, 437–447.
Agard, P., Jolivet, L., Vrielynck, B., Burov, E., Monié, P., 2007. Plate acceleration: the
obduction trigger? Earth and Planetary Science Letters 258, 428–441.
Aghamalyan, V.A., 2004. The Lesser Caucasus earth crust formation and evolution in the
collision zone of Palaeo-Tethys. 5th International Symposium on Eastern Mediterranean Geology, Proceedings, vol. 1, pp. 17–20. Thessaloniki, Greece.
Aslanyan, A.T., Satian, M.A., 1977. On the geological features of Transcaucasian ophiolitic
zones. Izvestia Academy of Science Armenian SSR 4–5, 13–26 (in Russian).
Aslanyan, A.T., Satian, M.A., 1982. Middle Cretaceous ophiolite zones of Transcaucasus
and tectonic reconstructions. Ofioliti 7, 131.
Avagyan, A., Sosson, M., Philip, M.H., Karakhanian, A., Rolland, Y., Melkonyan, R., Rebai,
S., Davtyan, V., 2005. Neogene to Quaternary stress field evolution in Lesser
Caucasus and adjacent regions using fault kinematics analysis and volcanic cluster
data. Geodinamica Acta 18, 401–416.
Baghdasaryan, G.P., Vardanyan, A.V., Satian, M.A., 1988. On the age of the volcanic rocks of
ophiolitic association. Izvestia Academy of Science Armenian SSR 6, 11–18 (in Russian).
Ben-Avraham, Z., Nur, A., Jones, D., Cox, A., 1981. Continental accretion: from oceanic
plateaus to allochtonous terranes. Science 213, 47–54.
Bosch, D., Blichert-Toft, J., Moynier, F., Nelson, B.K., Telouk, P., Gillot, P.Y., Albarède, F.,
2008. Pb, Hf and Nd isotope compositions of the two Réunion volcanoes, Indian
Ocean: a tale of two small-scale mantle “blobs”. Earth and Planetary Science Letters
265, 748–768.
Boudier, F., Ceuleneer, G., Nicolas, A., 1988. Shear zones, thrusts and related magmatism
in the Oman ophiolite: initiation of thrusting on an oceanic ridge. Tectonophysics
151, 275–296.
Brynzia, L.T., Wood, B.J., 1990. Oxygen thermobarometry of abyssal spinel peridotites;
the redox state and C–O–H volatile composition of the Earth's sub-oceanic upper
mantle. American Journal of Science 290, 1093–1116.
Cloos, M., 1993. Lithospheric buoyancy and collisional orogenesis: subduction of oceanic
plateaus, continental margins, island arcs, spreading ridges, and seamounts.
Geological Society of America Bulletin 105, 715–737.
Coleman, P.J., Kroenke, L.W., 1981. Subduction without volcanism in the Solomon
Islands Arc. Geo-Marine Letters 1, 129–134.
Danelian, T., Galoyan, G., Rolland, Y., Sosson, M., 2007. Palaeontological, (Radiolarian)
Late Jurassic age constraint for the Stepanavan ophiolite, Lesser Caucasus, Armenia.
Proceedings of the 11th International Congress, Athens, May 2007. Bulletin of the
Geological Society of Greece, vol. 37.
Danelian, T., Asatryan, G., Sosson, M., Person, A., Sahakyan, L., Galoyan, G., 2008. Discovery
of two distinct Middle Jurassic Radiolarian assemblages in the sedimentary cover of
the Vedi ophiolite, Lesser Caucasus, Armenia. Comptes Rendus Palevol 7, 327–334.
Dercourt, J., Zonenschain, L.P., Ricou, L.E., Kazmin, V.G., Le Pichon, X., Knipper, A.L.,
Grandjacquet, C., Sbortshikov, I.M., Geyssant, J., Lepvrier, C., Pechersky, D.H., Boulin,
J., Sibuet, J.C., Savostin, L.A., Sorokhtin, O., Westphal, M., Bazhenov, M.L., Lauer, J.P.,
Biju-Duval, B., 1986. Geological evolution of the tethys belt from the Atlantic to the
pamirs since the LIAS. Tectonophysics 123, 241–315.
Duncan, R.A., Hargreaves, R.B.,1984. Plate Tectonic Evolution of the Caribbean Region in the
Mantle Reference Frame. Geological Society of America Memoir 162, 81–93.
Eisele, J., Sharma, M., Galer, S.J.G., Blichert-Toft, J., Devey, C.W., Hofmann, A.W., 2002. The role
of sediment recycling in EM-1 inferred from Os, Pb, Hf, Nd, Sr isotope and trace element
systematics of the Pitcairn hotspot. Earth Planetary Science Letters 196, 197–212.
Evensen, N.M., Hamilton, P.J., O'Nios, R.K., 1978. Rare earth abundances in chondritic
meteorites. Geochimica et Cosmochimica Acta 42, 1199–1212.
Faure, G., 1986. Principles of isotope geology, 2nd edition. Wiley, New York.
Galoyan, G., Rolland, Y., Sosson, M., Corsini, M., Melkonyan, R., 2007. Evidence for
superposed MORB, oceanic plateau and volcanic arc series in the Lesser Caucasus,
Stepanavan, Armenia. Comptes Rendus Geosciences 339, 482–492.
Galoyan, G., Rolland, Y., Sosson, M., Corsini, M., Billo, S., Verati, C., Melkonyan, R., 2009.
Geochemistry and 40Ar/39Ar dating of Sevan Ophiolites, Lesser Caucasus, Armenia):
evidences for Jurassic Back-arc opening and hot spot event between the South
Armenian Block and Eurasia. Journal of Asian Earth Sciences 34, 135–153.
doi:10.1016/j.jseaes.2008.04.002.
Ghazaryan, H.A., 1987. Stratified gabbros of ophiolitic series of the south-eastern part of
the Sevan mountain range. Tipomorfizm and the Parageneses of the Minerals of
Armenian SSR. Izvestia Academy of Science Armenian SSR, pp. 122–139 (in Russian).
Ghazaryan, H.A., 1994. Pecularities of the geological structures and petrogenesis of
ophiolite gabbroids, on the examples of the Sevan and Vedi ophiolite zones of
Armenia. Izvestia Academy of Science Armenian SSR 3, 19–31 (in Russian).
Ghazi, A.M., Hassanipak, A.A., 1999. Geochemistry of subalkaline and alkaline extrusives
from the Kermanshah ophiolite, Zagros Suture Zone, Western Iran: implications for
Tethyan plate tectonics. Journal of Asian Earth Sciences 17, 319–332.
Gladczenko, T.P., Coffin, M.F., Eldholm, O., 1997. Crustal structure of the Ontong Java
Plateau: modeling of new gravity and existing seismic data. Journal of Geophysical
Research— Solid Earth 102, 22711–22729.
Guillot, S., Allemand, P., 2002. 2D thermal modelling of the early evolution of the
Himalayan belt. Journal of Geodynamics 34, 77–98.
186
Y. Rolland et al. / Lithos 112 (2009) 163–187
Guillot, S., Garzanti, E., Baratoux, D., Marquer, D., Mahéo, G., de Sigoyer, J., 2003.
Reconstructing the total shortening history of the NW Himalaya. Geochemistry
Geophysics Geosystems 4 (7), 1064. doi:10.1029/2002GC000484.
Hanan, B.D., Blichert-Toft, J., Kingsley, R., Schilling, J.G., 2000. Depleted Iceland mantle
plume geochemical signature: artifact of multi-component mixing? Geochemistry
Geophysics Geosystems 1, 1999GC000009.
Hart, S.R., 1984. A large-scale isotope anomaly in the Southern Hemisphere mantle.
Nature 309, 753–757.
Harutyunyan, G.S., 1967. The breakdown of the intrusions of the north-western part of the
Sevan ridge dependent on age. Izvestia Academy of Science Armenian SSR vol. 1–2,
42–52 (in Russian).
Hauff, F., Hoernle, K., Schmincke, H.U., Werner, R., 1997. A mid Cretaceous origin fort he
Galapagos hot-spot: volcanological, petrological and geochemical evidence from
Costa Rican oceanic crustal segments. Geologische Rundschau 86, 141–155.
Herzberg, C., O'Hara, M.J., 2002. Plume-associated ultramafic magmas of Phanerozoic
age. Journal of Petrology 43, 1857–1883.
Irvine, T.N., Baragar, W.R.A., 1971. A guide to the chemical classification of the common
volcanic rocks. Canadian Journal of Earth Sciences 8, 523–548.
Johnson, K.T.M., 1994. Experimental cpx/ and garnet/melt partitioning of REE and other
trace elements at high pressures; petrogenetic implications. Mineralogical
Magazine 58, 454–455.
Jourdan, F., Féraud, G., Bertrand, H., Kampunzu, A.B., Tshoso, G., Le Gall, B., Tiercelin, J.J.,
Capiez, P., 2004. The Karoo triple junction questioned: evidence for Jurassic and
Proterozoic 40Ar/39Ar ages and geochemistry of the giant Okavango dyke swarn,
Botswana. Earth and Planetary Science Letters 222, 989–2006.
Juteau, T., Maury, R., 1997. Géologie de la croûte océanique. Masson, Paris. 367 pp.
Kawahata, H., Nohara, M., Ishizuka, H., Hasebe, S., Chiba, H., 2001. Sr isotope
geochemistry and hydrothermal alteration of the Oman ophiolite. Journal of
Geophysical Research 106 (B6), 11083–11099.
Kazmin, V.G., Sbortshikov, I.M., Ricou, L.-E., Zonenshain, L.P., Boulin, J., Knipper, A.L., 1987.
Volcanic belt-indicators of the Mesozoic–Cenozoic active outskirts of Eurasia. In:
Monin, A.S., Zonenshain, L.P. (Eds.), History of the Tethys Ocean. Academy of Sciences
of the USSR, P.P. Shirshov Institute of Oceanology, Moscow, pp. 58–74 (in Russian).
Kerr, A.C., Mahoney, J.J., 2007. Oceanic plateaus: problematic plumes, potential
paradigms. Chemical Geology 241, 332–353.
Kerr, A.C., Tarney, J., Kempton, P.D., Spadea, P., Nivia, A., Marriner, G.F., Duncan, R.A.,
2002. Pervasive mantle plume head heterogeneity: evidence from the Late
Cretaceous Caribbean–Colombian oceanic plateau. Journal of Geophysical
Research–Solid Earth 107. doi:10.1029/2001JB000790.
Kerr, A.C., White, R.V., Thompson, P.M.E., Tarney, J., Saunders, A.D., 2003. No oceanic
plateau — no Caribbean Plate? The seminal role of an oceanic plateau. In: Car, C., Buffler,
R.T., Blickwede, J. (Eds.), The Gulf of Mexico and Caribbean Region: Hydrocarbon
habitats, Basin formation and Plate tectonics. AAPG Memoir, vol. 79, pp. 126–268.
Knipper, A.L., 1975. The oceanic crust in the alpine belt. Tr. GIN NAS USSR 267, 207
(in Russian).
Knipper, A.L., Khain, E.V., 1980. Structural position of ophiolites of the Caucasus. Ofioliti
Special Issue 2, 297–314.
Knipper, A.L., Sokolov, S.D., 1977. Vedi ophiolites (Armenia) autochthon or allochthon?
Geotektonics 10, 261–269.
Kocak, K., Isık, F., Arslan, M., Zedef, V., 2005. Petrological and source region
characteristics of ophiolitic hornblende gabbros from the Aksaray and Kayseri
regions, central Anatolian crystalline complex, Turkey. Journal of Asian Earth
Sciences 25, 883–891.
Lagabrielle, Y., 1987. Les ophiolites: marqueurs de l'histoire tectonique des domaines
océaniques. Thèse Doctorat d'Etat, Univ. Bretagne Occidentale, Brest, 350 p.
Lagabrielle, Y., Cannat, M., 1990. Alpine Jurassic ophiolites resemble the modern central
Atlantic basement. Geology 18, 319–322.
Lagabrielle, Y., Polino, R., Auzende, J.M., Blanchet, R., Caby, R., Fudral, S., Lemoine, M.,
Mevel, C., Ohnenstetter, M., Robert, D., Tricart, P., 1984. Les témoins d'une
tectonique intra-océanique dans le domaine téthysien: analyse des rapports entre
les ophiolites et leurs couvertures métasédimentaires dans la zone piémontaise des
Alpes franco-italiennes. Ofioliti 9, 67–88.
Lapierre, H., Samper, A., Bosch, D., Maury, R.C., Béchennec, F., Cotton, J., Demant, A.,
Brunet, P., Keller, F., Marcoux, J., 2004. The Tethyan plume: geochemical diversity of
Middle Permian basalts from the Oman rifted margin. Lithos 74, 167–198.
Leake, B.E., Woolley, A.R., Arps, C.E.S., Birch, W.D., Gilbert, M.C., Grice, J.D., Hawthorne, F.C.,
Kato, A., Kisch, H.J., Krivovichev, V.G., Linthout, K., Laird, J., Mandarino, J.A., Maresch, W.V.,
Nickel, E.H., Rock, N.M.S., Schumacher, J.C., Smith, D.C., Stephenson, N.C.N., Ungaretti, L.,
Whittaker, E.J.W., Youzhi, G., 1997. Nomenclature of amphiboles: report of the
subcommittee on amphiboles of the international mineralogical association, commission
on new minerals and mineral names. American Mineralogist 82, 1019–1037.
Le Maitre, R.W., Bateman, P., Dudek, A., Keller, J., Lameyre Le Bas, M.J., Sabine, P.A.,
Schmid, R., Sorensen, H., Streckeisen, A., Woolley, A.R., Zanettin, B., 1989. A
classification of igneous rocks and glossary of terms. Blackwell, Oxford.
Leterrier, J., Maury, R., Thonon, P., Girard, D., Marchal, M., 1982. Clinopyroxene
composition as a method of identification of the magmatic affinities of palaeovolcanic series. Earth and Planetary Science Letters 59, 139–154.
McCulloch, M.T., Gregory, R.T., Wasserburg, G.J., Taylor, J., 1981. Sm–Nd, Rb–Sr, and 18O/
16
O systematics in an oceanic crustal section: evidence from the Samail ophiolite.
Journal of Geophysical Research 86, 2721–2736.
McKenzie, D., O'Nions, R.K., 1991. Partial melt distribution from inversion of rare earth
element concentrations. Journal of Petrology 32, 1021–1091.
Meliksetyan, B.M., Baghdasaryan, G.P., Ghukasyan, R.Kh., 1984. Isotopic-geochemical
and geochronological investigations of eclogite–amphibolites associated with
ophiolites of Sevan–Amasian belt, Amasian Massif. Izvestia Academy of Science
Armenian SSR 1, 3–22 in Russian.
Melikyan, L.S., Palandjyan, S.A., Chibukhchyan, Z.H., Vardazaryan, J.S., 1967. To a
question about the geological position and the age of an ophiolitic series of the zone
Shirak–Sevan–Akera of the Lesser Caucasus. Izvestia Academy of Science Armenian
SSR 1–2, 21–41 in Russian.
Miura, S., Suyehiro, K., Shinohara, M., Takahashib, N., Araki, E., Taira, A., 2004.
Seismological structure and implications of collision between the Ontong Java
Plateau and Solomon Island Arc from ocean bottom seismometer airgun data.
Tectonophysics 389, 191–220.
Nicolas, A., 1989. Structures of ophiolites and dynamics of the oceanic lithosphere.
Kluwer Academy Publisher, Dordrecht, The Netherlands. 367 pp.
Nicolas, A., Jackson, E.D., 1972. Répartition en deux provinces des péridotites des
chaînes alpines longeant la Méditerranée: implications géotectoniques. Bulletin
Suisse de Minéralogie et Pétrologie 53, 385–401.
Nikogosian, I.K., Sobolev, A.V., 1997. Ion microprobe analysis of melt Inclusions in
olivine: experience in estimating the olivine-melt partition coefficients of trace
elements. Geochemistry of Interiors 35, 119–126.
Norman, T.N., 1984. The role of the Ankara Melange in the development of Anatolia,
Turkey. Geological Society Special Publications 17, 441–447.
O'Brien, P., Zotov, N., Law, R., Khan, A.M., Jan, M.Q., 2001. Coesite in Himalayan eclogite
and implication for models of India–Asia collision. Geology 29, 435–438.
Okay, A.I., 2000. Was the Late Triassic orogeny in Turkey caused by the collision of an
oceanic plateau? Geological Society Special Publications 173, 25–41.
Palandjyan, S.A., 1971. The petrology of ultrabasites and gabbroic rocks of the Sevan
mountain chain. Izvestia Academy of Science Armenian SSR. 201 pp. (in Russian).
Pallister, J.S., Knight, R.J., 1981. Rare earth element geochemistry of the Samail ophiolite
near Ibra, Oman. Journal of Geophysical Research 86, 2673–2697.
Parkinson, I.J., Pearce, J.A., 1998. Peridotites from the Izu–Bonin–Mariana forearc, ODP
Leg 125): evidence for mantle melting and melt–mantle interaction in a suprasubductive zone setting. Journal of Petrology 39, 1577–1618.
Pearce, J.A., Lippard, S.J., Roberts, S., 1984. Characteristics and tectonic significance of
supra-subduction zone ophiolites. In: Kokelaar, B.P., Howells, M.F. (Eds.), Marginal
Basin Geology. Geological Society Special Publications, vol. 15, pp. 77–94.
Petterson, M.G., Neal, C.R., Mahoney, J.J., Kroenke, L.W., Saunders, A.D., Babbs, T.L.,
Duncan, R.A., Tolia, D., McGrail, B., 1997. Structure and deformation of north and
central Malaita, Solomon Islands: tectonic implications for the Ontong Java Plateau–
Solomon arc collision, and for the fate of oceanic plateaus. Tectonophysics 283,
1–33.
Plank, T., Langmuir, C.H., 1998. The chemical composition of subducting sediment and
its consequences for the crust and mantle. Chemical Geology 145, 325–394.
Rapp, R.P., Watson, E.B., Miller, C.F., 1991. Partial melting of amphibolite/eclogite and the
origin of Archean trondhjemites and tonalites. Precambrian Research 94,
4619–4633.
Révillon, S., Arndt, N.T., Chauvel, C., Hallot, E., 2000. Geochemical study of ultramafic
volcanic and plutonic rocks from Gorgona Island Colombia: the plumbing system of
an oceanic plateau. Journal of Petrology 41, 1127–1153.
Richards, M.A., Jones, D.L., Duncan, R.A., DePaolo, D.J., 1991. A mantle plume initiation
model for the Wrangellia flood basalt and other oceanic plateus. Science 252,
263–267.
Richardson, W.P., Okal, E.A., VanderLee, S., 2000. Rayleigh-wave tomography of the
Ontong Java Plateau. Physics of the Earth and Planetary Interiors 118, 29–51.
Ricou, L.E., Zonenshain, L.P., Dercourt, J., Kazmin, V.G., Le Pichon, X., Knipper, A.L.,
Grandjacquet, C., Sborshchikov, I.M., Geyssant, J., Lepvrier, C., Pechersky, D.M.,
Boulin, J., Sibuet, J.C., Savostin, L.A., Sorokhtin, O., Westphal, M., Bazhenov, M.L.,
Lauer, J.P., Biju-Duval, B., 1985. Méthodes pour l'établissement de neuf cartes
paléogéographiques de l'Atlantique au Pamir depuis le Lias. Bulletin de la Société
Géologique de France 8, 625–635.
Rolland, Y., Billo, S., Corsini, M., Sosson, M., Galoyan, G., 2009. Blueschists of the
Amassia–Stepanavan Suture Zone, Armenia): linking Tethys subduction history
from E-Turkey to W-Iran. International Journal of Earth Sciences 98, 533–550.
Rolland, Y., Billo, S., Corsini, M., Sosson, M., Galoyan, G. in press. The Armenian Ophiolite:
insights for Jurassic Back-arc formation, Lower Cretaceous hot-spot magmatism, and
Upper Cretaceous obduction over the South Armenian Block. Geological Society
Special Publication.
Salters, V.J.M., White, W.M., 1998. Hf isotope constraints on mantle evolution. Chemical
Geology 145, 447–460.
Salters, V., Stracke, A., 2004. Composition of the depleted mantle. Geochemistry
Geophysics Geosystems 5. doi:10.1029/2003GC000597.
Saunders, A.D., Tarney, J., 1984. Geochemical characteristics of basaltic volcanism within
back-arc basins. In: Kokelaar, B.P., Howells, M.F. (Eds.), Marginal Basin Geology.
Geological Society Special Publications, vol. 16, pp. 59–76.
Satian, M.A., Sarkisyan, E.A., 2006. On lithodinamics of Mesozoic volcanogenesedimentary complex of the Sevan range, Sevan–Akera ophiolite zone, Armenia.
Izvestia Academy of Science Armenian SSR 2, 19–26 (in Russian).
Satian, M.A., Stepanyan, J.H., Sahakyan, L.H., Mnatsakanyan, A.Kh., Ghukasyan, R.Kh.,
2005. The Mesozoic lamprophyric diatremes of Vedi zone, Armenia. Izdatelstvo
‘Nairi’ Yerevan. 148 pp. (in Russian).
Sengör, A.M.C., 1984. The Cimmeride Orogenic System and the Tectonics of Eurasia.
Geological Society of America Special Paper, vol. 195, p. 82.
Sengör, A.M.C., 1990. A new model for the late Paleozoic–Mesozoic tectonic evolution of
Iran and implications for Oman. In: Robertson, A.H., Searle, M.P., Ries, A.C. (Eds.),
The Geology and tectonics of the Oman region. Geological Society Special
Publication, vol. 49, pp. 797–831.
Sengör, A.M.C., Yılmaz, Y., 1981. Tethyan evolution of Turkey: a plate tectonic approach.
Tectonophysics 75, 181–241.
Sokolov, S.D., 1977. The olistostroms and ophiolitic nappes of the Lesser Caucasus.
Izdatelstvo Nauka, Moscow. 92 pp. (in Russian).
Y. Rolland et al. / Lithos 112 (2009) 163–187
Sosson, M., Rolland, Y., Danelian, T., Muller, C., Melkonyan, R., Adamia, S., Kangarli, T.,
Avagyan, A., Galoyan, G., in press. Subductions, obduction and collision in the Lesser
Caucasus (Armenia, Azerbaijan, Georgia), new insights. Geological Society Special
Publication.
Stampfli, G.M., Borel, G.D., 2002. A plate tectonic model for the Palaeozoic and Mesozoic
constrained by dynamic plate boundaries and restored synthetic oceanic isochrons.
Earth and Planetary Science Letters 196, 17–33.
Staudigel, H., Plank, T., White, B., Schmincke, H.-U., 1996. Geochemical fluxes during
seafloor alteration of the basaltic upper oceanic crust: DSDP Sites 417 and 418. In:
Bebout, G.E., Scholl, D.W., Kirby, S.H., Platt, J.P. (Eds.), Subduction: Top to Bottom.
AGU Monograph, pp. 19–38.
Steiger, R.H., Jäger, E., 1977. Subcomission on geochronology: convention of the use of
decay constants in geo- and cosmochronology. Earth Planetary Science Letters 36,
359–362.
Storey, M., Mahoney, J.J., Kroenke, L.W., Saunders, A.D., 1991. Are oceanic plateaus sites
of komatiite formation? Geology 19, 376–379.
Sun, S.S., McDonough, W.F., 1989. Chemical and isotopic systematics of oceanic basalts:
implications for mantle composition and processes. In: Saudners, A.D., Norry, M.J.
(Eds.), Magmatism in Ocean Basins. Geological Society Special Publications, vol. 42,
pp. 313–345.
Tankut, A., Dilek, Y., Önen, P., 1998. Petrology and geochemistry of the Neo-Tethyan
volcanism as revealed in the Ankara melange, Turkey. Journal of Volcanology and
Geothermal Research 85, 265–284.
Tarnudo, J.A., Sliter, W.V., Kroenke, L.W., Leckie, M., Mahoney, J.J., Musgrave, R.J., Storey,
M., Winterer, E.L., 1991. Rapid formation of the Ontong Java Plateau by Aptian plume
volcanism. Science 254, 399–403.
187
Taylor, S.R., McLennan, S.M., 1985. The continental crust: its composition and evolution.
Blackwell. 312 pp.
Tirrul, R., Bell, R., Griffis, R.J., Camp, V.E., 1983. The Sistan Suture zone of eastern Iran.
Geological Society of America Bulletin 94, 134–150.
Turner, G., Huneke, J.C., Podose, F.A., Wasserburg, G.J., 1971. 40Ar/39Ar ages and cosmic
ray exposure ages of Apollo 14 samples. Earth and Planetary Science Letters 12,
15–19.
Tüysüz, O., Dellaloğlu, A.A., Terzioğlu, N., 1995. A magmatic belt within the Neo-Tethyan
suture zone and its role in the tectonic evolution of northern Turkey. Tectonophysics 243, 173–191.
Wasserburg, B.J., Jacobsen, S.B., De Paolo, D.J., McCulloch, M.T., Wen, T., 1981. Precise
determination of Sm/Nd ratios, Sm and Nd isotopic abundances in standard
solutions. Geochimica et Cosmochimica Acta 45, 2311–2323.
Zakariadze, G.S., Knipper, A.L., Sobolev, A.V., Tsameryan, O.P., Dimitriev, L.V.,
Vishnevskaya, V.S., Kolesov, G.M., 1983. The ophiolite volcanic series of the Lesser
Caucasus. Ofioliti 8, 439–466.
Zakariadze, G.S., Knipper, A.L., Bibikova, E.V., Silantiev, S.A., Zlobin, S.K., Gracheva, T.V.,
Makarov, S.A., Kolesov, T.M., 1990. The setting and age of the plutonic part of the NE
Sevan ophiolite complex. Academy of sciences of the USSR Geological series 3,
17–30 (in Russian).
Zindler, A., Hart, S.R., 1986. Chemical geodynamics. Annual Reviews in Earth and
Planetary Sciences 14, 493–571.