# Petrogenesis of Quaternary Shoshonitic Volcanism in NE Iran (Ardabil): Implication for Postcollisional Magmatism

Authors: Habib Shahbazi Shiran
Journal: Journal of Geological Research
 (2013)

---

## Abstract

Trachyandesites, trachytes, andesites, and pyrocalstic rocks, with shoshonitic signature, are the main Quaternary volcanic rocks in the Sabalan region (Ardabil). Plagiocalse, K-feldspar, biotite associated with clinopyroxene, and glass are the main constituents of these lavas. Plagioclases are andesine to labradorite while clinopyroxenes have augitic composition. The Sabalan volcanic rocks show enrichment in LREEs (relative to HREEs) and are characterized by enrichment in LILEs and depletion in HFSEs. Petrological observations, along with rare earth and trace elements geochemistry, suggest shoshonitic signature for Sabalan lavas. This signature highlights derivation from a subduction-related source. The Sabalan volcanic rocks are isotopically characterized by derivation from an enriched mantle source with a tendency to plot in the fields defined by island-arc basalts (IAB) and OIBs (inεNd versus 87Sr/86Sr diagram). The geochemical and isotopic characteristics of the Sabalan lavas suggest that their magma has been issued via low degree partial melting of a subduction-metasomatized continental lithospheric mantle. The formation of these lavas is related to slab steepening and breakoff in a postcollisional regime.


---

## Body

## 1. Introduction

Cenozoic magmatism is well known from the peri-Arabic region, north of the Bitlis-Zagros Suture Zone (Figure1), and ranges in age from Eocene to Plio-Quaternary. The timing of the eruptions and pooling the plutons mostly coincides with and postdates a series of continental collisional events in the northern Bitlis-Zagros Suture Zone region [1, 2]. The Eocene magmatism in the peri-Arabic region (both in the Urumieh-Dokhtar magmatic belt and south of the Anatolian block) is calc-alkalic to shoshonitic, which resulted from subduction of the Neotethyan Ocean beneath the central Iranian and Anatolian blocks. The Plio-Quaternary (and Quaternary) alkaline-ultrapotassic magmatism occupies most of the Turkish-Iranian high plateau (NW Iran and NE Turkey) and is mostly characterized by within-plate and/or subduction-related geochemical characteristics [3–5]. The geochemical variations of late Cenozoic lavas indicate a progression from calc-alkaline to alkaline compositions with time [2].Figure 1
Simplified tectonic map of the eastern Mediterranean-Persian Gulf region showing the active plate boundaries and post-collisional volcanic rocks in the Peri-Arabian region modified after [6]. Abbreviation: PSF = Pampak-Sevan Fault; EKP = Erzurum-Kars Plateau; EAF = East Anatolian Fault; IAESZ = Izmir-Ankara-Erzincan Suture Zone; EF = Ecemis Fault; NAF = North Anatolian Fault; DSFZ = Dead Sea Fault Zone.The Turkish-Iranian high plateau is bounded on the north by the Eastern Pontide arc and Lesser Caucasus magmatic belt and to the south by continental blocks including Bitlis, Puturg, Sanandaj, and Sirjan blocks. Plio-Quaternary and/or Quaternary volcanic cones and flows are scattered in this plateau (such as Mt. Ararat, Sahand, Sabalan, Nemrut, and others) and most of them indicate postcollisional shoshonitic-alkaline affinities [7–10]. The Sabalan volcano (Ardabil) is among the youngest volcanic calderas, in eastern border of the Turkish-Iranian Plateau, consisting of Miocene trachyandesites and Plio-Quaternary trachy-andesite, pumiceous andesites, and dacites with ultrapotassic-shoshonitic signatures. The aim of this study is to recognize the major, trace, REE and Sr-Nd isotopes characteristics of the Sabalan Quaternary lavas and to present a consistent tectonomagmatic model for the formation and evolution of the Sabalan volcano.


## 2. Regional Geology

The southwestern Sabalan region is characterized by Miocene magmatism with occurrence of calc-alkaline and high K calc-alkaline to shoshonitic lavas [11]. Moreover, late Miocene sodic alkaline lavas are common in south of the Sabalan region [12]. Eocene and Oligo-Miocene lavas are also distributed in NW regions of Ardabil, near Meshkin Shahr-Ahar volcanic belt. The Eocene lavas are mostly andesites, trachyandesites, trachybasalts and analcime-bearing tephrites with shoshonitic to alkaline geochemical signatures. Volcanic breccias, tuffs, and interbedded lavas, with columnar joints are common in the N-NE of the Sabalan region. Andesitic to analcime-bearing tephritic dikes crosscut the pyroclastic sequence, mainly at NE of Meshkin Shahr. Eocene granitic, monzonitic, and monzogabbroic intrusions are common in the Ahar-Meshkin Shahr magmatic belt. These plutonic rocks have shoshonitic geochemical affinity [13].The Sabalan volcano consists of both late Miocene (old) and Quaternary (young) trachyandesitic and trachytic to dacitic lavas with huge bodies of ignimbrites-density pyroclastic flows [14]. The Sabalan magmatic activity can be divided into two stages. Stage 1 is huge masses of ignimbrites and volcanic ashes associated with interbedded andesitic to trachyandesitic lavas, erupted in early stages. After several eruptional phases, during collapse of the central domain, a caldera was formed. Stage 2 is dominant with growth of dacitic to trachyandesitic domes in the central parts of the old caldera. The Quaternary Sabalan lavas were described as high K calc-alkaline to shoshonitic by Dostal and Zebri [15] with highly differentiated REE patterns in chondrite-normalized diagrams. Didon and Germain [16] also suggested high K calc-alkaline and shoshonitic signatures for the Sabalan volcanic rocks.Trachybasalts are minor in the Sabalan volcano and occur mostly as old, Miocene lava flows (Figure2). They are distributed in the southwestern parts of the Sabalan caldera. These rocks grade upward into trachyandesitic lavas and occasionally are interbedded with pyroclastic rocks. Trachyandesites-porphyritic andesites occur mostly in the southern to western parts of the Sabalan caldera and are associated with pumiceous andesites. Trachytes are characterized by the presence of amphibole and biotite crystals. These rocks usually are vesicular. Quaternary porphyritic trachyandesites are mostly distributed in the central parts of the caldera and are the main constituents of the Quaternary Sabalan lavas (Figure 2). These lavas (trachyandesites) are intebedded with ignimberites. Dacites are minor in the Sabalan region and occur mostly as aphyric lavas associated with trachyandesites.Figure 2
Simplified geological map of the Sabalan region with emphasis on the distribution of Sabalan Miocene to Quaternary and Paleogene volcanic rocks (modified after 1/100000 Geological Map of Meshkin Shahr, Geological Survey of Iran).Ignimbrites constitute the major components of the pyroclastic eruptions and contain compressed fiammes with andesitic composition. Pumice and scoria flows are common. Lahar deposits, showing the latest stages of the volcanic activity, are predominant in the Sabalan region.


## 3. Petrography of Sabalan Quaternary Volcanic Rocks

Here, we adopted a classification based on the occurrence type of the rocks in the field, petrographical studies and mineral modes. In this way, we subdivided the rock samples into six groups.(1)
Trachyandesites. These rocks are porphyritic with large (1-2 cm) K-feldspar phenocrysts with slight alteration into clay minerals. Biotite is another rock-forming mineral and occurs as euhedral to subhedral phenocrysts. Plagioclase shows zoning and sometimes alteration into clay minerals and sericite. The groundmass indicates alteration into clay and chlorite, pseudomorph phases after glass shards. Plagioclase and K-feldspar microlites occur in the groundmass. Amphibole in trachyandesites is rare while clinopyroxene in more mafic varieties is common. These rocks show hyaloporphyritic texture.(2)
Dacitic Lavas. K-feldspar, plagioclase, biotite, and quartz are predominant phases of dacitic lavas. These minerals occur as small crystals within groundmass and/or as microphenocrysts. K-feldspar and plagioclase show slight alteration into sericite and clay minerals. Biotite displays brown paleochroism and is fresh. These rocks are mostly aphyric and most samples show holohyaline texture.(3)
Trachytes. Trachytes have more K-feldspar than plagioclase and biotite and no quartz compared to dacites. Plagioclase phenocrysts show oscillatory to patchy zoning. K-feldspar shows slight alteration into clay minerals. Plagioclase microlites and glass shards are common in the groundmass. These rocks show intersertal to hyalo-microlitic porphyritic texture.(4)
Andesites. Andesites are characterized by higher amounts of plagioclase than K-feldspar and occurrence of clinopyroxene microphenocrysts. Plagioclases show normal to oscillatory zoning. Plagioclases are less altered into clay than K-feldspar. Clinopyroxenes are subhedral to anhedral and occasionally show twining. Plagioclase microlites are dominant in the groundmass of the rocks.(5)
Ignimbrites. These rocks are characterized by high amount of biotite associated with K-feldspar, plagioclase, and rarely quartz with flowed glass shards (eutaxitic texture). These rocks sometimes show layering, dark layers rich in biotite, and white layers rich in feldspar, due to the gravity flow of the crystal aggregates after explosion of the volcano. Oxidized andesitic fiammes are characteristic of some ignimbrite flows.(6)
Pumice and Scoria. Andesitic scoria flows with large and abundant vesicles are common. Pumices and compressed pumices with fiammes are other rocks that accompany the pyroclastic units. Andesitic scoria flows have large amounts of glass, with rare plagioclase and clinopyroxene microphenocrysts.


## 4. Geochemistry of Quaternary Lavas

### 4.1. Analytical Methods

For this study we selected ten fresh samples from Quaternary Sabalan volcanic rocks for major, trace, and REE elements analyzed at ALS, Canada, using ICP-AES for major elements and ICP-MS for trace and REE elements. The uncertainty (1 sigma) is ~2% for major elements and 5 to 10% for trace elements (depending on concentration). Mineral analyses have been performed at Paris VI University using Cameca SX 100. Accelerating voltage, beam current, and beam diameter for the analyses were 20 kV, 20 nA, and 3μm, respectively.The isotopic analyses were carried out at the Laboratório de Geologia Isotópica da Universidade de Aveiro (Portugal). For Sr-Nd isotopes, samples were dissolved with HF/HNO3 in Teflon Parr acid digestion bombs at 200°C. After evaporation of the final solution, the samples were dissolved with HCl (6 N) and dried down. The elements for analysis were purified using a conventional two-stage ion chromatography technique: (i) separation of Sr and REE elements in ion exchange column with AG8 50 W Bio-Rad cation exchange resin; (ii) purification of Nd from other lanthanide elements in columns with Ln resin (ElChroM Technologies) cation exchange resin. All reagents used in sample preparation were subboiling distilled, and water was produced by a Milli-Q Element (Millipore) apparatus. Sr was loaded on a single Ta filament with H3PO4, whereas Nd was loaded on a Ta side filament with HCl, in a triple filament arrangement. Both elements were determined using a Multi-Collector Thermal Ionisation Mass Spectrometer (TIMS) VG Sector 54. Data were obtained in dynamic mode with peak measurements at 1-2 V for 88Sr and 0.8–1.5 V for 144Nd. Sr and Nd isotopic ratios were corrected for mass fractionation relative to 88Sr/86Sr = 0.1194 and 146Nd/144Nd = 0.7219. During this study, the SRM-987 standard gave an average value of 87Sr/86Sr = 0.710263 (6) (N=20; 95% c.l.) and the JNdi-1 standard yielded 143Nd/144Nd = 0.512098 (2) (N=29; 95% c.l.) (143Nd/144Nd data are normalized to the La Jolla Nd standard).


### 4.2. Mineral Chemistry

In this section, we focus on mineral chemistry of Quaternary Sabalan trachyandesites, the predominant rocks in the region.


#### 4.2.1. Plagioclase

Plagioclases in trachyandesites are characterized by andesine to labradorite composition (37 to 65% An, Table1). Their K2O content ranges between 0.3 and 1.4 wt% (Table 1). Their variable An content indicates spots from various bands of zoning within mineral.Table 1
Composition of clinopyroxene, plagioclase, and glass shards in the Sabalan volcanic rocks.
(a)
Name
S2
S2
S5
S5
S5
S9
S9
S9
S9
S9
S4
S4
S4
Mineral
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
SiO2
53.16
52.94
51.33
52.40
52.56
50.13
52.32
51.67
52.81
50.47
50.64
52.66
52.16
TiO2
0.82
0.47
0.76
0.95
0.42
0.80
0.46
0.71
0.55
0.95
0.59
0.28
0.45
Al2O3
2.22
1.99
2.44
2.52
1.34
2.67
1.22
1.93
1.11
2.10
2.78
2.71
2.51
Cr2O3
0.14
0.00
0.31
0.00
0.16
0.00
0.00
0.02
0.00
0.00
0.02
0.72
0.18
FeO*
8.26
8.94
8.74
8.93
9.47
8.80
9.07
8.38
8.86
9.44
8.39
6.36
6.46
MnO
0.37
0.74
0.57
0.51
0.59
0.55
0.36
0.19
0.53
0.42
0.36
0.31
0.21
MgO
15.07
15.01
14.67
14.29
14.84
14.87
14.78
15.19
15.39
15.35
15.37
16.60
16.81
CaO
20.05
19.37
19.80
19.93
19.70
20.10
19.16
20.51
20.12
20.06
20.24
19.60
20.19
Na2O
0.46
0.46
0.59
0.45
0.49
0.37
0.34
0.46
0.42
0.68
0.32
0.38
0.40
K2O
0.00
0.01
0.00
0.00
0.00
0.02
0.01
0.05
0.03
0.05
0.05
0.01
0.07
NiO
0.07
0.19
0.00
0.03
0.00
0.00
0.00
0.01
0.18
0.03
0.06
0.08
0.00
BaO
0.28
0.00
0.00
0.00
0.00
0.00
0.21
0.00
0.00
0.00
0.00
0.43
0.14
F
0.00
0.00
0.00
0.25
0.08
0.00
0.12
0.00
0.16
0.12
0.00
0.04
0.00
Cl
0.00
0.00
0.02
0.05
0.09
0.00
0.00
0.00
0.00
0.05
0.03
0.00
0.00
Total
100.89
100.10
99.23
100.32
99.74
98.31
98.05
99.11
100.14
99.70
98.83
100.18
99.57
Mg#
0.76
0.75
0.75
0.74
0.74
0.75
0.74
0.76
0.76
0.74
0.77
0.82
0.82
Cr#
0.04
0.00
0.08
0.00
0.07
0.00
0.00
0.01
0.00
0.00
0.00
0.15
0.04
En
0.44
0.44
0.43
0.43
0.43
0.43
0.44
0.44
0.44
0.44
0.44
0.48
0.48
Fs
0.14
0.15
0.15
0.15
0.15
0.14
0.15
0.14
0.14
0.15
0.14
0.10
0.10
Wo
0.42
0.41
0.42
0.43
0.41
0.42
0.41
0.43
0.42
0.41
0.42
0.41
0.42
Si
1.95
1.96
1.93
1.94
1.97
1.90
1.98
1.94
1.96
1.90
1.91
1.94
1.93
Ti
0.02
0.01
0.02
0.03
0.01
0.02
0.01
0.02
0.02
0.03
0.02
0.01
0.01
Al
0.10
0.09
0.11
0.11
0.06
0.12
0.05
0.09
0.05
0.09
0.12
0.12
0.11
Cr
0.00
0.00
0.01
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.02
0.01
Fe*
0.25
0.28
0.27
0.28
0.30
0.28
0.29
0.26
0.28
0.30
0.26
0.20
0.20
Mn
0.01
0.02
0.02
0.02
0.02
0.02
0.01
0.01
0.02
0.01
0.01
0.01
0.01
Mg
0.83
0.83
0.82
0.79
0.83
0.84
0.83
0.85
0.85
0.86
0.86
0.91
0.93
Ca
0.79
0.77
0.80
0.79
0.79
0.82
0.78
0.82
0.80
0.81
0.82
0.77
0.80
Na
0.03
0.03
0.04
0.03
0.04
0.03
0.03
0.03
0.03
0.05
0.02
0.03
0.03
K
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
Ni
0.00
0.01
0.00
0.00
0.00
0.00
0.00
0.00
0.01
0.00
0.00
0.00
0.00
Total
3.99
4.00
4.02
3.99
4.01
4.03
3.99
4.02
4.01
4.05
4.03
4.00
4.02
(b)
Name
S2
S2
S2
S9
S9
S9
S9
S5
S5
S4
S4
S4
Mineral
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
SiO2
56.984
57.374
57.797
58.691
58.503
58.330
55.219
51.093
50.666
52.338
52.592
55.365
TiO2
0.162
0.067
0.080
0.080
0.117
0.052
0.073
0.073
0.022
0.058
0.000
0.022
Al2O3
26.210
26.159
25.617
25.602
24.411
23.408
26.422
29.910
29.617
28.542
28.959
27.208
Cr2O3
0.000
0.063
0.000
0.032
0.000
0.000
0.032
0.000
0.114
0.032
0.000
0.292
FeO
0.542
0.593
0.683
0.450
0.536
0.400
0.596
0.533
0.938
0.691
0.517
1.065
MnO
0.061
0.000
0.000
0.087
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
MgO
0.028
0.017
0.013
0.007
0.023
0.023
0.040
0.030
0.101
0.061
0.055
0.088
CaO
8.802
8.795
7.656
7.085
7.094
7.526
9.544
12.947
12.451
12.060
12.005
9.777
Na2O
5.583
5.826
6.191
6.558
6.683
6.862
5.769
4.245
3.722
4.808
4.873
5.364
K2O
0.942
1.006
1.203
1.083
1.341
1.443
0.918
0.261
0.255
0.425
0.313
0.597
NiO
0.000
0.041
0.041
0.000
0.055
0.000
0.000
0.136
0.000
0.055
0.000
0.150
BaO
0.000
0.000
0.143
0.215
0.358
0.000
0.071
0.071
0.215
0.071
0.144
0.071
F
0.000
0.000
0.000
0.205
0.041
0.123
0.124
0.168
0.000
0.000
0.000
0.083
Cl
0.000
0.016
0.005
0.000
0.026
0.032
0.016
0.000
0.068
0.000
0.011
0.032
Total
99.314
99.955
99.430
100.095
99.189
98.199
98.824
99.466
98.169
99.141
99.469
100.115
An%
46.560
45.482
40.596
37.384
36.972
37.737
47.760
62.762
64.896
58.091
57.653
50.181
Si
2.582
2.587
2.618
2.639
2.662
2.680
2.535
2.350
2.358
2.409
2.409
2.511
Ti
0.006
0.002
0.003
0.003
0.004
0.002
0.003
0.003
0.001
0.002
0.000
0.001
Al
1.400
1.390
1.368
1.357
1.309
1.267
1.430
1.621
1.624
1.548
1.563
1.454
Cr
0.000
0.002
0.000
0.001
0.000
0.000
0.001
0.000
0.004
0.001
0.000
0.010
Fe*
0.021
0.022
0.026
0.017
0.020
0.015
0.023
0.021
0.037
0.027
0.020
0.040
Mn
0.002
0.000
0.000
0.003
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
Mg
0.002
0.001
0.001
0.000
0.002
0.002
0.003
0.002
0.007
0.004
0.004
0.006
Ca
0.427
0.425
0.372
0.341
0.346
0.370
0.469
0.638
0.621
0.595
0.589
0.475
Na
0.491
0.509
0.544
0.572
0.590
0.611
0.514
0.379
0.336
0.429
0.433
0.472
K
0.054
0.058
0.070
0.062
0.078
0.085
0.054
0.015
0.015
0.025
0.018
0.035
Ni
0.000
0.001
0.001
0.000
0.002
0.000
0.000
0.005
0.000
0.002
0.000
0.005
Total
4.985
4.998
5.002
4.996
5.013
5.032
5.031
5.034
5.003
5.042
5.035
5.009
(c)
Name
S2
S2
S5
S5
S9
S9
S9
S4
S4
Mineral
Glass
Glass
Glass
Glass
Glass
Glass
Glass
Glass
Glass
SiO2
75.32
75.748
73.41
74.528
68.305
72.28
66.705
72.307
72.793
TiO2
0.829
0.827
0.807
0.756
0.212
0.952
0.33
1.113
1.146
Al2O3
12.863
12.567
14.018
14.063
17.432
13.495
18.55
14.152
14.029
Cr2O3
0.064
0.048
0.032
0.032
0.032
0.114
0
0
0.08
FeO
0.773
1.853
0.773
0.732
0.578
1.4
0.908
1.424
2.7
MnO
0
0
0.009
0.043
0.094
0
0
0
0.035
MgO
0.041
0.041
0.003
0.025
0.036
0.051
0.046
0.051
0.065
CaO
0.315
0.259
0.186
0.259
1.087
0.316
2.163
0.327
0.227
Na2O
2.599
2.205
2.536
2.99
5.346
3.317
6.42
2.062
2.092
K2O
6.125
7.171
7.93
8.169
6.145
5.993
3.645
7.377
5.853
NiO
0
0.027
0
0
0.055
0.041
0
0.055
0
BaO
0
0
0.071
0
0.644
0
0
0
0
F
0.21
0
0.042
0
0.041
0.376
0.285
0
0
Cl
0.203
0.117
0.096
0.074
0.058
0.165
0.085
0.159
0.042
Total
99.343
100.863
99.727
101.671
100.065
98.5
99.139
99.028
99.063
Mg
#: Mg number, Feo*: total iron.


#### 4.2.2. Clinopyroxene

Clinopyroxene of trachyandesites are characterized by augitic composition (Figure3(a)). Enstatite (En) end-member of these clinopyroxenes ranges from 43 to 48% while their ferrosillite (Fs) component is low, varying between 10 and 15% (Table 1). TiO2 and Al2O3 range between 0.28–0.95 and 1.1–2.8 wt%, respectively, while Cr2O3 and Na2O contents show lower values of 0.1–0.7 and 0.3–0.7 wt%, respectively (Table 1). Their Ti and Na contents are different from those in alkaline rocks, characterized by high concentrations of these elements (titano-augite or aegirine-augite). In Ti + Cr versus Ca diagram of [17], trachyandesite clinopyroxenes have a tendency to plot in volcanic arc field (Figure 3(b)).(a) Pyroxene composition of the Sabalan lavas in Wo-En-Fs diagram. (b) Composition of pyroxenes in a plot of Ti + Cr versus Ca [17], indicating their formation in an arc-related setting.
(a)
(b)


#### 4.2.3. Glass Shards

Glass shards in the groundmass of trachyandesites have minor traces of alteration and then were selected for EMP analysis. Their SiO2 content is variable, from 66.7 to 75.7 wt% (Table 1), higher than SiO2 content of whole rocks (57.6–65.4 wt%). The K2O content of these glasses is between 3.6 and 8.1 wt%, with low content of MgO (0.003–0.07 wt%) and TiO2 (0.2–1.1 wt%).


### 4.3. Whole Rock Chemistry

The Sabalan lavas have slight alteration, characterized by low content of LOI values (1.1–1.9 wt%, Table2), except sample S08-7 with 3.2 wt% LOI. The lavas are intermediate to acidic in composition, with SiO2 content ranging from 57.6 to 65.4 wt% (Table 2). In total alkalis versus SiO2 diagram of [18], the Sabalan lavas (bulk rock) show tendency to plot in trachyandesitic to trachydacitic fields, similar to other shoshonitic lavas from NW Iran (Figure 4). These lavas tend to plot in trachydacite and rhyolite domains based on the composition of the glass shards of the rock groundmass. K2O and TiO2 contents of the Sabalan lavas vary between 2.5–4.8 and 0.7–0.9 wt% (Table 2). In K2O against SiO2 diagram [19], the lavas show affinity to high-K calc-alkaline and shoshonitic rocks. This geochemical signature is common for other postcollisional rocks from NW Iran and NE Turkey [2] (Figure 5). In chondrite-normalized REE diagram (Figure 6), Sabalan lavas show a fractioned REE trend, with enrichment in LREEs relative to HREEs (La(n)/Yb(n) = 18.8–20.4). Some samples show slight depletion in Eu, indicating plagioclase fractionation.Table 2
Whole rock and Sr-Nd isotope composition of the Sabalan volcanic rocks.
Name
S08-3
S08-1
S08-2
S08-4
S08-5
S08-6
S08-7
S08-8
S08-9
S08-10
SiO2
61.2
58.8
65.4
60.4
64.9
57.6
63.5
64.4
64.6
65.2
Al2O3
16.8
16.8
15
16.5
15.2
16.5
14.6
14.9
14.95
15
FeOt
5.09
6.06
3.69
5.08
3.97
6.08
3.62
3.73
3.82
3.75
CaO
5.45
6.17
2.62
5.37
2.64
5.8
2.61
2.69
2.73
2.65
MgO
2.4
2.46
1.09
2.41
1.08
2.99
1.08
1.1
1.1
1.06
Na2O
3.57
4.01
4.11
3.52
4.25
3.97
4.08
4.06
4.03
4.02
K2O
2.86
2.83
4.68
2.81
4.6
2.49
4.17
4.58
4.69
4.78
TiO2
0.7
0.89
0.7
0.7
0.71
0.88
0.68
0.71
0.71
0.7
MnO
0.1
0.1
0.09
0.1
0.09
0.11
0.09
0.09
0.09
0.08
P2O5
0.27
0.42
0.25
0.26
0.25
0.41
0.24
0.25
0.28
0.25
LOI
1.54
1.68
1.32
1.49
1.88
1.09
3.24
1.36
1.59
1.45
Total
100
100.5
99.2
98.8
99.8
98.1
98.1
98
98.7
99.1
Ba
844
1500
829
809
1140
784
1015
748
744
775
Ce
78.5
98.1
127.5
78.3
127.5
90
124.5
125.5
123.5
125
Cr
60
90
450
40
140
90
80
30
230
20
Cs
2.07
1.63
5.07
2.01
4.98
1.5
4.92
4.81
4.78
4.89
Dy
3.15
3.84
4.06
3.18
4.2
3.55
4.04
3.99
4.06
4.06
Er
1.88
2.27
2.48
1.87
2.6
2.01
2.47
2.5
2.45
2.56
Eu
1.29
1.65
1.36
1.31
1.4
1.53
1.33
1.37
1.39
1.35
Gd
4.62
5.78
6.27
4.54
6.43
5.29
6.14
6.13
6.15
6.18
Hf
4.9
5.7
8.7
4.8
9
5.1
8.5
8.6
8.5
8.7
Ho
0.64
0.77
0.82
0.63
0.85
0.7
0.84
0.81
0.82
0.83
La
46.6
57.5
73.4
46.3
73.5
52
71.9
72
70.4
72
Lu
0.27
0.32
0.41
0.27
0.42
0.27
0.39
0.39
0.39
0.4
Nb
20.9
28.3
48.8
21
49.3
25.1
47.1
46.9
46.4
47.5
Nd
27.6
35.4
41.4
27.4
41.5
32.3
41
40.7
40.7
40.4
Ni
28
32
8
28
9
36
6
7
8
7
Pb
13
16
21
14
23
13
20
24
19
23
Pr
7.9
10.05
12.45
7.92
12.5
9.1
12.4
12.4
12.25
12.3
Rb
72
71.5
180.5
72.4
180
64.5
173.5
176.5
176
179
Sm
4.79
5.88
6.61
4.66
6.62
5.46
6.52
6.58
6.46
6.44
Sr
626
837
401
622
410
751
448
415
421
415
Ta
1.4
1.8
3.4
1.4
3.6
1.6
3.3
3.4
3.3
3.4
Tb
0.61
0.75
0.79
0.59
0.83
0.69
0.78
0.77
0.77
0.78
Th
9.28
10.1
26.3
9.2
26.7
9.08
25.7
25.5
25.1
25.8
Tm
0.26
0.31
0.38
0.26
0.38
0.28
0.36
0.36
0.36
0.37
U
2.75
2.67
7.11
2.72
7.11
2.37
6.87
6.89
6.81
6.98
V
100
121
59
99
58
118
58
58
60
57
Y
17.7
21.2
23.7
17.9
23.9
19.2
23.5
23.4
23.3
23.6
Yb
1.78
2.02
2.61
1.77
2.7
1.85
2.58
2.55
2.52
2.56
Zn
72
89
102
98
81
89
72
95
68
92
Zr
210
253
368
211
376
230
356
356
354
357
87Sr/86Sr
0.70511
0.70452
0.70443
143Nd/144Nd
0.51273
0.51269
0.51271
εNd
1.8
1.04
1.4Figure 4
Total alkali versus SiO2 classification diagram of the Sabalan volcanic rocks [18]. Data source of NW Iran lavas is from [2].Figure 5
K2O versus SiO2 diagram [19] of Sabalan volcanic rocks. Data source of NW Iran lavas is from [2].Chondrite-normalized REE and mantle-normalized trace elements patterns of the Sabalan volcanic rocks. Normalization data are from [20].The Sabalan lavas show conspicuous depletion in Nb, Ta, Ti, and P (e.g., Nb(n)/La(n) = 0.5-0.6) with positive anomalies in most LIL elements including Rb, Th, U, and Pb (e.g., Th(n)/La(n) = 1.4–2.9 and U(n)/La(n) = 1.5–3.2) (Figure 6). These rare earth and trace elements patterns are mostly consistent with shoshonitic signature of these rocks, indicating magmatism along an active continental margin.


### 4.4. Nd-Sr Isotopes

Trachyandesitic samples have uniform initial (t = 2 Ma) 87Sr/86Sr values (0.704-0.705) with εNd ranging from ~+1 to +1.8. In the initial 87Sr/86Sr versus εNd variation diagram (Figure 7), the Sabalan volcanic rocks are characterized by derivation from an enriched mantle source, and all samples show tendency to plot in a field defined by island-arc (IAB) and ocean-island (OIB) basalts. The samples mostly show affinity to plot near the Ghoshchi alkali gabbros (unpublished data), Suffi-abad I-type granites [21], and Eastern Pontides adakites [22] (Figure 7). The Sabalan shoshonitic lavas have higher εNd than Quaternary Salmas lavas (unpublished data), indicating mantle source of the Sabalan lavas slightly depleted compared to that of Salmas lavas.εNd versus 87Sr/86Sr diagrams for Sabalan lavas. Data from other Turkish-Iranian postcollisional intrusive and extrusive rocks are shown for comparison; right-hand diagram is modified after [23].


## 4.1. Analytical Methods

For this study we selected ten fresh samples from Quaternary Sabalan volcanic rocks for major, trace, and REE elements analyzed at ALS, Canada, using ICP-AES for major elements and ICP-MS for trace and REE elements. The uncertainty (1 sigma) is ~2% for major elements and 5 to 10% for trace elements (depending on concentration). Mineral analyses have been performed at Paris VI University using Cameca SX 100. Accelerating voltage, beam current, and beam diameter for the analyses were 20 kV, 20 nA, and 3μm, respectively.The isotopic analyses were carried out at the Laboratório de Geologia Isotópica da Universidade de Aveiro (Portugal). For Sr-Nd isotopes, samples were dissolved with HF/HNO3 in Teflon Parr acid digestion bombs at 200°C. After evaporation of the final solution, the samples were dissolved with HCl (6 N) and dried down. The elements for analysis were purified using a conventional two-stage ion chromatography technique: (i) separation of Sr and REE elements in ion exchange column with AG8 50 W Bio-Rad cation exchange resin; (ii) purification of Nd from other lanthanide elements in columns with Ln resin (ElChroM Technologies) cation exchange resin. All reagents used in sample preparation were subboiling distilled, and water was produced by a Milli-Q Element (Millipore) apparatus. Sr was loaded on a single Ta filament with H3PO4, whereas Nd was loaded on a Ta side filament with HCl, in a triple filament arrangement. Both elements were determined using a Multi-Collector Thermal Ionisation Mass Spectrometer (TIMS) VG Sector 54. Data were obtained in dynamic mode with peak measurements at 1-2 V for 88Sr and 0.8–1.5 V for 144Nd. Sr and Nd isotopic ratios were corrected for mass fractionation relative to 88Sr/86Sr = 0.1194 and 146Nd/144Nd = 0.7219. During this study, the SRM-987 standard gave an average value of 87Sr/86Sr = 0.710263 (6) (N=20; 95% c.l.) and the JNdi-1 standard yielded 143Nd/144Nd = 0.512098 (2) (N=29; 95% c.l.) (143Nd/144Nd data are normalized to the La Jolla Nd standard).


## 4.2. Mineral Chemistry

In this section, we focus on mineral chemistry of Quaternary Sabalan trachyandesites, the predominant rocks in the region.


### 4.2.1. Plagioclase

Plagioclases in trachyandesites are characterized by andesine to labradorite composition (37 to 65% An, Table1). Their K2O content ranges between 0.3 and 1.4 wt% (Table 1). Their variable An content indicates spots from various bands of zoning within mineral.Table 1
Composition of clinopyroxene, plagioclase, and glass shards in the Sabalan volcanic rocks.
(a)
Name
S2
S2
S5
S5
S5
S9
S9
S9
S9
S9
S4
S4
S4
Mineral
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
SiO2
53.16
52.94
51.33
52.40
52.56
50.13
52.32
51.67
52.81
50.47
50.64
52.66
52.16
TiO2
0.82
0.47
0.76
0.95
0.42
0.80
0.46
0.71
0.55
0.95
0.59
0.28
0.45
Al2O3
2.22
1.99
2.44
2.52
1.34
2.67
1.22
1.93
1.11
2.10
2.78
2.71
2.51
Cr2O3
0.14
0.00
0.31
0.00
0.16
0.00
0.00
0.02
0.00
0.00
0.02
0.72
0.18
FeO*
8.26
8.94
8.74
8.93
9.47
8.80
9.07
8.38
8.86
9.44
8.39
6.36
6.46
MnO
0.37
0.74
0.57
0.51
0.59
0.55
0.36
0.19
0.53
0.42
0.36
0.31
0.21
MgO
15.07
15.01
14.67
14.29
14.84
14.87
14.78
15.19
15.39
15.35
15.37
16.60
16.81
CaO
20.05
19.37
19.80
19.93
19.70
20.10
19.16
20.51
20.12
20.06
20.24
19.60
20.19
Na2O
0.46
0.46
0.59
0.45
0.49
0.37
0.34
0.46
0.42
0.68
0.32
0.38
0.40
K2O
0.00
0.01
0.00
0.00
0.00
0.02
0.01
0.05
0.03
0.05
0.05
0.01
0.07
NiO
0.07
0.19
0.00
0.03
0.00
0.00
0.00
0.01
0.18
0.03
0.06
0.08
0.00
BaO
0.28
0.00
0.00
0.00
0.00
0.00
0.21
0.00
0.00
0.00
0.00
0.43
0.14
F
0.00
0.00
0.00
0.25
0.08
0.00
0.12
0.00
0.16
0.12
0.00
0.04
0.00
Cl
0.00
0.00
0.02
0.05
0.09
0.00
0.00
0.00
0.00
0.05
0.03
0.00
0.00
Total
100.89
100.10
99.23
100.32
99.74
98.31
98.05
99.11
100.14
99.70
98.83
100.18
99.57
Mg#
0.76
0.75
0.75
0.74
0.74
0.75
0.74
0.76
0.76
0.74
0.77
0.82
0.82
Cr#
0.04
0.00
0.08
0.00
0.07
0.00
0.00
0.01
0.00
0.00
0.00
0.15
0.04
En
0.44
0.44
0.43
0.43
0.43
0.43
0.44
0.44
0.44
0.44
0.44
0.48
0.48
Fs
0.14
0.15
0.15
0.15
0.15
0.14
0.15
0.14
0.14
0.15
0.14
0.10
0.10
Wo
0.42
0.41
0.42
0.43
0.41
0.42
0.41
0.43
0.42
0.41
0.42
0.41
0.42
Si
1.95
1.96
1.93
1.94
1.97
1.90
1.98
1.94
1.96
1.90
1.91
1.94
1.93
Ti
0.02
0.01
0.02
0.03
0.01
0.02
0.01
0.02
0.02
0.03
0.02
0.01
0.01
Al
0.10
0.09
0.11
0.11
0.06
0.12
0.05
0.09
0.05
0.09
0.12
0.12
0.11
Cr
0.00
0.00
0.01
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.02
0.01
Fe*
0.25
0.28
0.27
0.28
0.30
0.28
0.29
0.26
0.28
0.30
0.26
0.20
0.20
Mn
0.01
0.02
0.02
0.02
0.02
0.02
0.01
0.01
0.02
0.01
0.01
0.01
0.01
Mg
0.83
0.83
0.82
0.79
0.83
0.84
0.83
0.85
0.85
0.86
0.86
0.91
0.93
Ca
0.79
0.77
0.80
0.79
0.79
0.82
0.78
0.82
0.80
0.81
0.82
0.77
0.80
Na
0.03
0.03
0.04
0.03
0.04
0.03
0.03
0.03
0.03
0.05
0.02
0.03
0.03
K
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
Ni
0.00
0.01
0.00
0.00
0.00
0.00
0.00
0.00
0.01
0.00
0.00
0.00
0.00
Total
3.99
4.00
4.02
3.99
4.01
4.03
3.99
4.02
4.01
4.05
4.03
4.00
4.02
(b)
Name
S2
S2
S2
S9
S9
S9
S9
S5
S5
S4
S4
S4
Mineral
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
SiO2
56.984
57.374
57.797
58.691
58.503
58.330
55.219
51.093
50.666
52.338
52.592
55.365
TiO2
0.162
0.067
0.080
0.080
0.117
0.052
0.073
0.073
0.022
0.058
0.000
0.022
Al2O3
26.210
26.159
25.617
25.602
24.411
23.408
26.422
29.910
29.617
28.542
28.959
27.208
Cr2O3
0.000
0.063
0.000
0.032
0.000
0.000
0.032
0.000
0.114
0.032
0.000
0.292
FeO
0.542
0.593
0.683
0.450
0.536
0.400
0.596
0.533
0.938
0.691
0.517
1.065
MnO
0.061
0.000
0.000
0.087
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
MgO
0.028
0.017
0.013
0.007
0.023
0.023
0.040
0.030
0.101
0.061
0.055
0.088
CaO
8.802
8.795
7.656
7.085
7.094
7.526
9.544
12.947
12.451
12.060
12.005
9.777
Na2O
5.583
5.826
6.191
6.558
6.683
6.862
5.769
4.245
3.722
4.808
4.873
5.364
K2O
0.942
1.006
1.203
1.083
1.341
1.443
0.918
0.261
0.255
0.425
0.313
0.597
NiO
0.000
0.041
0.041
0.000
0.055
0.000
0.000
0.136
0.000
0.055
0.000
0.150
BaO
0.000
0.000
0.143
0.215
0.358
0.000
0.071
0.071
0.215
0.071
0.144
0.071
F
0.000
0.000
0.000
0.205
0.041
0.123
0.124
0.168
0.000
0.000
0.000
0.083
Cl
0.000
0.016
0.005
0.000
0.026
0.032
0.016
0.000
0.068
0.000
0.011
0.032
Total
99.314
99.955
99.430
100.095
99.189
98.199
98.824
99.466
98.169
99.141
99.469
100.115
An%
46.560
45.482
40.596
37.384
36.972
37.737
47.760
62.762
64.896
58.091
57.653
50.181
Si
2.582
2.587
2.618
2.639
2.662
2.680
2.535
2.350
2.358
2.409
2.409
2.511
Ti
0.006
0.002
0.003
0.003
0.004
0.002
0.003
0.003
0.001
0.002
0.000
0.001
Al
1.400
1.390
1.368
1.357
1.309
1.267
1.430
1.621
1.624
1.548
1.563
1.454
Cr
0.000
0.002
0.000
0.001
0.000
0.000
0.001
0.000
0.004
0.001
0.000
0.010
Fe*
0.021
0.022
0.026
0.017
0.020
0.015
0.023
0.021
0.037
0.027
0.020
0.040
Mn
0.002
0.000
0.000
0.003
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
Mg
0.002
0.001
0.001
0.000
0.002
0.002
0.003
0.002
0.007
0.004
0.004
0.006
Ca
0.427
0.425
0.372
0.341
0.346
0.370
0.469
0.638
0.621
0.595
0.589
0.475
Na
0.491
0.509
0.544
0.572
0.590
0.611
0.514
0.379
0.336
0.429
0.433
0.472
K
0.054
0.058
0.070
0.062
0.078
0.085
0.054
0.015
0.015
0.025
0.018
0.035
Ni
0.000
0.001
0.001
0.000
0.002
0.000
0.000
0.005
0.000
0.002
0.000
0.005
Total
4.985
4.998
5.002
4.996
5.013
5.032
5.031
5.034
5.003
5.042
5.035
5.009
(c)
Name
S2
S2
S5
S5
S9
S9
S9
S4
S4
Mineral
Glass
Glass
Glass
Glass
Glass
Glass
Glass
Glass
Glass
SiO2
75.32
75.748
73.41
74.528
68.305
72.28
66.705
72.307
72.793
TiO2
0.829
0.827
0.807
0.756
0.212
0.952
0.33
1.113
1.146
Al2O3
12.863
12.567
14.018
14.063
17.432
13.495
18.55
14.152
14.029
Cr2O3
0.064
0.048
0.032
0.032
0.032
0.114
0
0
0.08
FeO
0.773
1.853
0.773
0.732
0.578
1.4
0.908
1.424
2.7
MnO
0
0
0.009
0.043
0.094
0
0
0
0.035
MgO
0.041
0.041
0.003
0.025
0.036
0.051
0.046
0.051
0.065
CaO
0.315
0.259
0.186
0.259
1.087
0.316
2.163
0.327
0.227
Na2O
2.599
2.205
2.536
2.99
5.346
3.317
6.42
2.062
2.092
K2O
6.125
7.171
7.93
8.169
6.145
5.993
3.645
7.377
5.853
NiO
0
0.027
0
0
0.055
0.041
0
0.055
0
BaO
0
0
0.071
0
0.644
0
0
0
0
F
0.21
0
0.042
0
0.041
0.376
0.285
0
0
Cl
0.203
0.117
0.096
0.074
0.058
0.165
0.085
0.159
0.042
Total
99.343
100.863
99.727
101.671
100.065
98.5
99.139
99.028
99.063
Mg
#: Mg number, Feo*: total iron.


### 4.2.2. Clinopyroxene

Clinopyroxene of trachyandesites are characterized by augitic composition (Figure3(a)). Enstatite (En) end-member of these clinopyroxenes ranges from 43 to 48% while their ferrosillite (Fs) component is low, varying between 10 and 15% (Table 1). TiO2 and Al2O3 range between 0.28–0.95 and 1.1–2.8 wt%, respectively, while Cr2O3 and Na2O contents show lower values of 0.1–0.7 and 0.3–0.7 wt%, respectively (Table 1). Their Ti and Na contents are different from those in alkaline rocks, characterized by high concentrations of these elements (titano-augite or aegirine-augite). In Ti + Cr versus Ca diagram of [17], trachyandesite clinopyroxenes have a tendency to plot in volcanic arc field (Figure 3(b)).(a) Pyroxene composition of the Sabalan lavas in Wo-En-Fs diagram. (b) Composition of pyroxenes in a plot of Ti + Cr versus Ca [17], indicating their formation in an arc-related setting.
(a)
(b)


### 4.2.3. Glass Shards

Glass shards in the groundmass of trachyandesites have minor traces of alteration and then were selected for EMP analysis. Their SiO2 content is variable, from 66.7 to 75.7 wt% (Table 1), higher than SiO2 content of whole rocks (57.6–65.4 wt%). The K2O content of these glasses is between 3.6 and 8.1 wt%, with low content of MgO (0.003–0.07 wt%) and TiO2 (0.2–1.1 wt%).


## 4.2.1. Plagioclase

Plagioclases in trachyandesites are characterized by andesine to labradorite composition (37 to 65% An, Table1). Their K2O content ranges between 0.3 and 1.4 wt% (Table 1). Their variable An content indicates spots from various bands of zoning within mineral.Table 1
Composition of clinopyroxene, plagioclase, and glass shards in the Sabalan volcanic rocks.
(a)
Name
S2
S2
S5
S5
S5
S9
S9
S9
S9
S9
S4
S4
S4
Mineral
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
Cpx
SiO2
53.16
52.94
51.33
52.40
52.56
50.13
52.32
51.67
52.81
50.47
50.64
52.66
52.16
TiO2
0.82
0.47
0.76
0.95
0.42
0.80
0.46
0.71
0.55
0.95
0.59
0.28
0.45
Al2O3
2.22
1.99
2.44
2.52
1.34
2.67
1.22
1.93
1.11
2.10
2.78
2.71
2.51
Cr2O3
0.14
0.00
0.31
0.00
0.16
0.00
0.00
0.02
0.00
0.00
0.02
0.72
0.18
FeO*
8.26
8.94
8.74
8.93
9.47
8.80
9.07
8.38
8.86
9.44
8.39
6.36
6.46
MnO
0.37
0.74
0.57
0.51
0.59
0.55
0.36
0.19
0.53
0.42
0.36
0.31
0.21
MgO
15.07
15.01
14.67
14.29
14.84
14.87
14.78
15.19
15.39
15.35
15.37
16.60
16.81
CaO
20.05
19.37
19.80
19.93
19.70
20.10
19.16
20.51
20.12
20.06
20.24
19.60
20.19
Na2O
0.46
0.46
0.59
0.45
0.49
0.37
0.34
0.46
0.42
0.68
0.32
0.38
0.40
K2O
0.00
0.01
0.00
0.00
0.00
0.02
0.01
0.05
0.03
0.05
0.05
0.01
0.07
NiO
0.07
0.19
0.00
0.03
0.00
0.00
0.00
0.01
0.18
0.03
0.06
0.08
0.00
BaO
0.28
0.00
0.00
0.00
0.00
0.00
0.21
0.00
0.00
0.00
0.00
0.43
0.14
F
0.00
0.00
0.00
0.25
0.08
0.00
0.12
0.00
0.16
0.12
0.00
0.04
0.00
Cl
0.00
0.00
0.02
0.05
0.09
0.00
0.00
0.00
0.00
0.05
0.03
0.00
0.00
Total
100.89
100.10
99.23
100.32
99.74
98.31
98.05
99.11
100.14
99.70
98.83
100.18
99.57
Mg#
0.76
0.75
0.75
0.74
0.74
0.75
0.74
0.76
0.76
0.74
0.77
0.82
0.82
Cr#
0.04
0.00
0.08
0.00
0.07
0.00
0.00
0.01
0.00
0.00
0.00
0.15
0.04
En
0.44
0.44
0.43
0.43
0.43
0.43
0.44
0.44
0.44
0.44
0.44
0.48
0.48
Fs
0.14
0.15
0.15
0.15
0.15
0.14
0.15
0.14
0.14
0.15
0.14
0.10
0.10
Wo
0.42
0.41
0.42
0.43
0.41
0.42
0.41
0.43
0.42
0.41
0.42
0.41
0.42
Si
1.95
1.96
1.93
1.94
1.97
1.90
1.98
1.94
1.96
1.90
1.91
1.94
1.93
Ti
0.02
0.01
0.02
0.03
0.01
0.02
0.01
0.02
0.02
0.03
0.02
0.01
0.01
Al
0.10
0.09
0.11
0.11
0.06
0.12
0.05
0.09
0.05
0.09
0.12
0.12
0.11
Cr
0.00
0.00
0.01
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.02
0.01
Fe*
0.25
0.28
0.27
0.28
0.30
0.28
0.29
0.26
0.28
0.30
0.26
0.20
0.20
Mn
0.01
0.02
0.02
0.02
0.02
0.02
0.01
0.01
0.02
0.01
0.01
0.01
0.01
Mg
0.83
0.83
0.82
0.79
0.83
0.84
0.83
0.85
0.85
0.86
0.86
0.91
0.93
Ca
0.79
0.77
0.80
0.79
0.79
0.82
0.78
0.82
0.80
0.81
0.82
0.77
0.80
Na
0.03
0.03
0.04
0.03
0.04
0.03
0.03
0.03
0.03
0.05
0.02
0.03
0.03
K
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
Ni
0.00
0.01
0.00
0.00
0.00
0.00
0.00
0.00
0.01
0.00
0.00
0.00
0.00
Total
3.99
4.00
4.02
3.99
4.01
4.03
3.99
4.02
4.01
4.05
4.03
4.00
4.02
(b)
Name
S2
S2
S2
S9
S9
S9
S9
S5
S5
S4
S4
S4
Mineral
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
Plag
SiO2
56.984
57.374
57.797
58.691
58.503
58.330
55.219
51.093
50.666
52.338
52.592
55.365
TiO2
0.162
0.067
0.080
0.080
0.117
0.052
0.073
0.073
0.022
0.058
0.000
0.022
Al2O3
26.210
26.159
25.617
25.602
24.411
23.408
26.422
29.910
29.617
28.542
28.959
27.208
Cr2O3
0.000
0.063
0.000
0.032
0.000
0.000
0.032
0.000
0.114
0.032
0.000
0.292
FeO
0.542
0.593
0.683
0.450
0.536
0.400
0.596
0.533
0.938
0.691
0.517
1.065
MnO
0.061
0.000
0.000
0.087
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
MgO
0.028
0.017
0.013
0.007
0.023
0.023
0.040
0.030
0.101
0.061
0.055
0.088
CaO
8.802
8.795
7.656
7.085
7.094
7.526
9.544
12.947
12.451
12.060
12.005
9.777
Na2O
5.583
5.826
6.191
6.558
6.683
6.862
5.769
4.245
3.722
4.808
4.873
5.364
K2O
0.942
1.006
1.203
1.083
1.341
1.443
0.918
0.261
0.255
0.425
0.313
0.597
NiO
0.000
0.041
0.041
0.000
0.055
0.000
0.000
0.136
0.000
0.055
0.000
0.150
BaO
0.000
0.000
0.143
0.215
0.358
0.000
0.071
0.071
0.215
0.071
0.144
0.071
F
0.000
0.000
0.000
0.205
0.041
0.123
0.124
0.168
0.000
0.000
0.000
0.083
Cl
0.000
0.016
0.005
0.000
0.026
0.032
0.016
0.000
0.068
0.000
0.011
0.032
Total
99.314
99.955
99.430
100.095
99.189
98.199
98.824
99.466
98.169
99.141
99.469
100.115
An%
46.560
45.482
40.596
37.384
36.972
37.737
47.760
62.762
64.896
58.091
57.653
50.181
Si
2.582
2.587
2.618
2.639
2.662
2.680
2.535
2.350
2.358
2.409
2.409
2.511
Ti
0.006
0.002
0.003
0.003
0.004
0.002
0.003
0.003
0.001
0.002
0.000
0.001
Al
1.400
1.390
1.368
1.357
1.309
1.267
1.430
1.621
1.624
1.548
1.563
1.454
Cr
0.000
0.002
0.000
0.001
0.000
0.000
0.001
0.000
0.004
0.001
0.000
0.010
Fe*
0.021
0.022
0.026
0.017
0.020
0.015
0.023
0.021
0.037
0.027
0.020
0.040
Mn
0.002
0.000
0.000
0.003
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
Mg
0.002
0.001
0.001
0.000
0.002
0.002
0.003
0.002
0.007
0.004
0.004
0.006
Ca
0.427
0.425
0.372
0.341
0.346
0.370
0.469
0.638
0.621
0.595
0.589
0.475
Na
0.491
0.509
0.544
0.572
0.590
0.611
0.514
0.379
0.336
0.429
0.433
0.472
K
0.054
0.058
0.070
0.062
0.078
0.085
0.054
0.015
0.015
0.025
0.018
0.035
Ni
0.000
0.001
0.001
0.000
0.002
0.000
0.000
0.005
0.000
0.002
0.000
0.005
Total
4.985
4.998
5.002
4.996
5.013
5.032
5.031
5.034
5.003
5.042
5.035
5.009
(c)
Name
S2
S2
S5
S5
S9
S9
S9
S4
S4
Mineral
Glass
Glass
Glass
Glass
Glass
Glass
Glass
Glass
Glass
SiO2
75.32
75.748
73.41
74.528
68.305
72.28
66.705
72.307
72.793
TiO2
0.829
0.827
0.807
0.756
0.212
0.952
0.33
1.113
1.146
Al2O3
12.863
12.567
14.018
14.063
17.432
13.495
18.55
14.152
14.029
Cr2O3
0.064
0.048
0.032
0.032
0.032
0.114
0
0
0.08
FeO
0.773
1.853
0.773
0.732
0.578
1.4
0.908
1.424
2.7
MnO
0
0
0.009
0.043
0.094
0
0
0
0.035
MgO
0.041
0.041
0.003
0.025
0.036
0.051
0.046
0.051
0.065
CaO
0.315
0.259
0.186
0.259
1.087
0.316
2.163
0.327
0.227
Na2O
2.599
2.205
2.536
2.99
5.346
3.317
6.42
2.062
2.092
K2O
6.125
7.171
7.93
8.169
6.145
5.993
3.645
7.377
5.853
NiO
0
0.027
0
0
0.055
0.041
0
0.055
0
BaO
0
0
0.071
0
0.644
0
0
0
0
F
0.21
0
0.042
0
0.041
0.376
0.285
0
0
Cl
0.203
0.117
0.096
0.074
0.058
0.165
0.085
0.159
0.042
Total
99.343
100.863
99.727
101.671
100.065
98.5
99.139
99.028
99.063
Mg
#: Mg number, Feo*: total iron.


## 4.2.2. Clinopyroxene

Clinopyroxene of trachyandesites are characterized by augitic composition (Figure3(a)). Enstatite (En) end-member of these clinopyroxenes ranges from 43 to 48% while their ferrosillite (Fs) component is low, varying between 10 and 15% (Table 1). TiO2 and Al2O3 range between 0.28–0.95 and 1.1–2.8 wt%, respectively, while Cr2O3 and Na2O contents show lower values of 0.1–0.7 and 0.3–0.7 wt%, respectively (Table 1). Their Ti and Na contents are different from those in alkaline rocks, characterized by high concentrations of these elements (titano-augite or aegirine-augite). In Ti + Cr versus Ca diagram of [17], trachyandesite clinopyroxenes have a tendency to plot in volcanic arc field (Figure 3(b)).(a) Pyroxene composition of the Sabalan lavas in Wo-En-Fs diagram. (b) Composition of pyroxenes in a plot of Ti + Cr versus Ca [17], indicating their formation in an arc-related setting.
(a)
(b)


## 4.2.3. Glass Shards

Glass shards in the groundmass of trachyandesites have minor traces of alteration and then were selected for EMP analysis. Their SiO2 content is variable, from 66.7 to 75.7 wt% (Table 1), higher than SiO2 content of whole rocks (57.6–65.4 wt%). The K2O content of these glasses is between 3.6 and 8.1 wt%, with low content of MgO (0.003–0.07 wt%) and TiO2 (0.2–1.1 wt%).


## 4.3. Whole Rock Chemistry

The Sabalan lavas have slight alteration, characterized by low content of LOI values (1.1–1.9 wt%, Table2), except sample S08-7 with 3.2 wt% LOI. The lavas are intermediate to acidic in composition, with SiO2 content ranging from 57.6 to 65.4 wt% (Table 2). In total alkalis versus SiO2 diagram of [18], the Sabalan lavas (bulk rock) show tendency to plot in trachyandesitic to trachydacitic fields, similar to other shoshonitic lavas from NW Iran (Figure 4). These lavas tend to plot in trachydacite and rhyolite domains based on the composition of the glass shards of the rock groundmass. K2O and TiO2 contents of the Sabalan lavas vary between 2.5–4.8 and 0.7–0.9 wt% (Table 2). In K2O against SiO2 diagram [19], the lavas show affinity to high-K calc-alkaline and shoshonitic rocks. This geochemical signature is common for other postcollisional rocks from NW Iran and NE Turkey [2] (Figure 5). In chondrite-normalized REE diagram (Figure 6), Sabalan lavas show a fractioned REE trend, with enrichment in LREEs relative to HREEs (La(n)/Yb(n) = 18.8–20.4). Some samples show slight depletion in Eu, indicating plagioclase fractionation.Table 2
Whole rock and Sr-Nd isotope composition of the Sabalan volcanic rocks.
Name
S08-3
S08-1
S08-2
S08-4
S08-5
S08-6
S08-7
S08-8
S08-9
S08-10
SiO2
61.2
58.8
65.4
60.4
64.9
57.6
63.5
64.4
64.6
65.2
Al2O3
16.8
16.8
15
16.5
15.2
16.5
14.6
14.9
14.95
15
FeOt
5.09
6.06
3.69
5.08
3.97
6.08
3.62
3.73
3.82
3.75
CaO
5.45
6.17
2.62
5.37
2.64
5.8
2.61
2.69
2.73
2.65
MgO
2.4
2.46
1.09
2.41
1.08
2.99
1.08
1.1
1.1
1.06
Na2O
3.57
4.01
4.11
3.52
4.25
3.97
4.08
4.06
4.03
4.02
K2O
2.86
2.83
4.68
2.81
4.6
2.49
4.17
4.58
4.69
4.78
TiO2
0.7
0.89
0.7
0.7
0.71
0.88
0.68
0.71
0.71
0.7
MnO
0.1
0.1
0.09
0.1
0.09
0.11
0.09
0.09
0.09
0.08
P2O5
0.27
0.42
0.25
0.26
0.25
0.41
0.24
0.25
0.28
0.25
LOI
1.54
1.68
1.32
1.49
1.88
1.09
3.24
1.36
1.59
1.45
Total
100
100.5
99.2
98.8
99.8
98.1
98.1
98
98.7
99.1
Ba
844
1500
829
809
1140
784
1015
748
744
775
Ce
78.5
98.1
127.5
78.3
127.5
90
124.5
125.5
123.5
125
Cr
60
90
450
40
140
90
80
30
230
20
Cs
2.07
1.63
5.07
2.01
4.98
1.5
4.92
4.81
4.78
4.89
Dy
3.15
3.84
4.06
3.18
4.2
3.55
4.04
3.99
4.06
4.06
Er
1.88
2.27
2.48
1.87
2.6
2.01
2.47
2.5
2.45
2.56
Eu
1.29
1.65
1.36
1.31
1.4
1.53
1.33
1.37
1.39
1.35
Gd
4.62
5.78
6.27
4.54
6.43
5.29
6.14
6.13
6.15
6.18
Hf
4.9
5.7
8.7
4.8
9
5.1
8.5
8.6
8.5
8.7
Ho
0.64
0.77
0.82
0.63
0.85
0.7
0.84
0.81
0.82
0.83
La
46.6
57.5
73.4
46.3
73.5
52
71.9
72
70.4
72
Lu
0.27
0.32
0.41
0.27
0.42
0.27
0.39
0.39
0.39
0.4
Nb
20.9
28.3
48.8
21
49.3
25.1
47.1
46.9
46.4
47.5
Nd
27.6
35.4
41.4
27.4
41.5
32.3
41
40.7
40.7
40.4
Ni
28
32
8
28
9
36
6
7
8
7
Pb
13
16
21
14
23
13
20
24
19
23
Pr
7.9
10.05
12.45
7.92
12.5
9.1
12.4
12.4
12.25
12.3
Rb
72
71.5
180.5
72.4
180
64.5
173.5
176.5
176
179
Sm
4.79
5.88
6.61
4.66
6.62
5.46
6.52
6.58
6.46
6.44
Sr
626
837
401
622
410
751
448
415
421
415
Ta
1.4
1.8
3.4
1.4
3.6
1.6
3.3
3.4
3.3
3.4
Tb
0.61
0.75
0.79
0.59
0.83
0.69
0.78
0.77
0.77
0.78
Th
9.28
10.1
26.3
9.2
26.7
9.08
25.7
25.5
25.1
25.8
Tm
0.26
0.31
0.38
0.26
0.38
0.28
0.36
0.36
0.36
0.37
U
2.75
2.67
7.11
2.72
7.11
2.37
6.87
6.89
6.81
6.98
V
100
121
59
99
58
118
58
58
60
57
Y
17.7
21.2
23.7
17.9
23.9
19.2
23.5
23.4
23.3
23.6
Yb
1.78
2.02
2.61
1.77
2.7
1.85
2.58
2.55
2.52
2.56
Zn
72
89
102
98
81
89
72
95
68
92
Zr
210
253
368
211
376
230
356
356
354
357
87Sr/86Sr
0.70511
0.70452
0.70443
143Nd/144Nd
0.51273
0.51269
0.51271
εNd
1.8
1.04
1.4Figure 4
Total alkali versus SiO2 classification diagram of the Sabalan volcanic rocks [18]. Data source of NW Iran lavas is from [2].Figure 5
K2O versus SiO2 diagram [19] of Sabalan volcanic rocks. Data source of NW Iran lavas is from [2].Chondrite-normalized REE and mantle-normalized trace elements patterns of the Sabalan volcanic rocks. Normalization data are from [20].The Sabalan lavas show conspicuous depletion in Nb, Ta, Ti, and P (e.g., Nb(n)/La(n) = 0.5-0.6) with positive anomalies in most LIL elements including Rb, Th, U, and Pb (e.g., Th(n)/La(n) = 1.4–2.9 and U(n)/La(n) = 1.5–3.2) (Figure 6). These rare earth and trace elements patterns are mostly consistent with shoshonitic signature of these rocks, indicating magmatism along an active continental margin.


## 4.4. Nd-Sr Isotopes

Trachyandesitic samples have uniform initial (t = 2 Ma) 87Sr/86Sr values (0.704-0.705) with εNd ranging from ~+1 to +1.8. In the initial 87Sr/86Sr versus εNd variation diagram (Figure 7), the Sabalan volcanic rocks are characterized by derivation from an enriched mantle source, and all samples show tendency to plot in a field defined by island-arc (IAB) and ocean-island (OIB) basalts. The samples mostly show affinity to plot near the Ghoshchi alkali gabbros (unpublished data), Suffi-abad I-type granites [21], and Eastern Pontides adakites [22] (Figure 7). The Sabalan shoshonitic lavas have higher εNd than Quaternary Salmas lavas (unpublished data), indicating mantle source of the Sabalan lavas slightly depleted compared to that of Salmas lavas.εNd versus 87Sr/86Sr diagrams for Sabalan lavas. Data from other Turkish-Iranian postcollisional intrusive and extrusive rocks are shown for comparison; right-hand diagram is modified after [23].


## 5. Discussion

### 5.1. Petrogenesis

The Sabalan volcanic rocks are characterized by a fractionated REE pattern, with high LREE/HREE ratio, resembling shoshonitic rocks from active continental margins and/or postcollisional magmatism after closure of oceanic basins. Depletion in high field strength elements and enrichment in large lithophile elements are conspicuous. Such shoshonitic volcanism is characteristic of the Turkish-Iranian high plateau [2, 24]. The glass shards in the groundmass of these rocks are characterized by more fractionated nature, with higher SiO2 and K2O contents, indicating quenching products of more-evolved melt after fractionation of mafic minerals. Trachyandesite clinopyroxenes have a tendency to plot at volcanic arc rocks field in Ti + Cr versus Ca diagram, indicating an arc-related environment for the formation of Sabalan volcanic rocks. These characteristics show that these clinopyroxenes are crystallized from an arc-related calc-alkaline magma. The Sabalan lavas are similar to Goshchi alkali gabbros both of them are characterized by derivation from an isotopically enriched mantle source (OIB mantle source). They are clearly different from crustal contaminated, postcollisional Islamic peninsula shoshonites and ultrapotassic rocks. When compared to other Iranian postcollisional Plio-Quaternary lavas, like the Salmas basalts, the Sabalan rocks show higher εNd values, clearly indicating a slight depleted mantle source for the Sabalan lavas and/or contamination of the Salmas lavas with crustal materials during their ascent to the surface. The interaction with crustal rocks and/or derivation from a subduction-contaminated source is also clear from Th/Yb versus Ta/Yb diagram [25] for the Sabalan lavas (Figure 8). In this diagram, the Sabalan lavas show higher Th/Yb ratio, far from within plate mantle trend, resembling those lavas that erupted above an active continental margin. Higher Th content could have also resulted from interaction of ascending lavas with crustal rocks, via assimilation-fractional crystallization (AFC) processes. The Sabalan lavas seem to show contamination with high Th (and high 87Sr/86Sr) crustal rocks during their ascent to the surface beneath the volcano. This could be inferred from rare, partially dissolved crustal xenoliths in the Sabalan lavas, which can be evidence for crust contamination. In order to evaluate the mantle-melting regime of the Sabalan lavas, we used the Sm/Yb versus La/Sm diagram [26]. The Sabalan lavas are relatively similar to Salmas Quaternary basalts, characterized by low degree of partial melting (<0.1%) of a mantle source containing clinopyroxene-garnet lherzolites (Figure 9). Because the Sabalan lavas are characterized by higher La/Sm ratio, it seems that their mantle source contains more modal garnet than model source in Figure 9 and/or contains other minor phases including phlogopite and/or apatite.Figure 8
Th/Yb versus Ta/Yb diagram [25] for the Sabalan Quaternary lavas.Figure 9
Sm/Yb versus La/Sm diagram for the Sabalan volcanic rocks with melting curves obtained using the nonmodal batch melting modelling, modified after [26].


### 5.2. Geodynamic Setting

Plio-Quaternary (to Quaternary) volcanism in the Peri-Arabian region (in Turkish-Iranian high plateau) postdates the subduction of Neotethyan lithosphere since ~13 Ma [2]. These volcanic rocks become more alkaline with time and towards the south [2, 4, 5, 27]. However, most of volcanic rocks, like the Salmas and/or Maku Quaternary lavas, still show a SSZ geochemical signature including enrichment in LILE and/or depletion in HFSE, despite the lack of subduction system at that time. Geochemical and geodynamic evidence also shows a subduction setting for the genesis of Sabalan lavas: fractionated REEs patterns (high LREEs/HREEs ratio), enrichment in Th, Ba, Rb, U, Pb, and depletion in Nb, Ta, and Ti. Such geochemical characteristics are common in subduction-related environments [25]. These observations and isotopic characteristics of the lavas suggest that their magmas were probably derived from partial melting of a subduction-metasomatized continental lithospheric mantle in the garnet/spinel lherzolite field [2]. Widespread volcanism across the Turkish-Iranian plateau throughout the late Cenozoic until Plio-Quaternary is strictly ascribed to postcollisional mantle upwelling during slab steeping and breakoff beneath a subduction-accretion complex [4, 28].


## 5.1. Petrogenesis

The Sabalan volcanic rocks are characterized by a fractionated REE pattern, with high LREE/HREE ratio, resembling shoshonitic rocks from active continental margins and/or postcollisional magmatism after closure of oceanic basins. Depletion in high field strength elements and enrichment in large lithophile elements are conspicuous. Such shoshonitic volcanism is characteristic of the Turkish-Iranian high plateau [2, 24]. The glass shards in the groundmass of these rocks are characterized by more fractionated nature, with higher SiO2 and K2O contents, indicating quenching products of more-evolved melt after fractionation of mafic minerals. Trachyandesite clinopyroxenes have a tendency to plot at volcanic arc rocks field in Ti + Cr versus Ca diagram, indicating an arc-related environment for the formation of Sabalan volcanic rocks. These characteristics show that these clinopyroxenes are crystallized from an arc-related calc-alkaline magma. The Sabalan lavas are similar to Goshchi alkali gabbros both of them are characterized by derivation from an isotopically enriched mantle source (OIB mantle source). They are clearly different from crustal contaminated, postcollisional Islamic peninsula shoshonites and ultrapotassic rocks. When compared to other Iranian postcollisional Plio-Quaternary lavas, like the Salmas basalts, the Sabalan rocks show higher εNd values, clearly indicating a slight depleted mantle source for the Sabalan lavas and/or contamination of the Salmas lavas with crustal materials during their ascent to the surface. The interaction with crustal rocks and/or derivation from a subduction-contaminated source is also clear from Th/Yb versus Ta/Yb diagram [25] for the Sabalan lavas (Figure 8). In this diagram, the Sabalan lavas show higher Th/Yb ratio, far from within plate mantle trend, resembling those lavas that erupted above an active continental margin. Higher Th content could have also resulted from interaction of ascending lavas with crustal rocks, via assimilation-fractional crystallization (AFC) processes. The Sabalan lavas seem to show contamination with high Th (and high 87Sr/86Sr) crustal rocks during their ascent to the surface beneath the volcano. This could be inferred from rare, partially dissolved crustal xenoliths in the Sabalan lavas, which can be evidence for crust contamination. In order to evaluate the mantle-melting regime of the Sabalan lavas, we used the Sm/Yb versus La/Sm diagram [26]. The Sabalan lavas are relatively similar to Salmas Quaternary basalts, characterized by low degree of partial melting (<0.1%) of a mantle source containing clinopyroxene-garnet lherzolites (Figure 9). Because the Sabalan lavas are characterized by higher La/Sm ratio, it seems that their mantle source contains more modal garnet than model source in Figure 9 and/or contains other minor phases including phlogopite and/or apatite.Figure 8
Th/Yb versus Ta/Yb diagram [25] for the Sabalan Quaternary lavas.Figure 9
Sm/Yb versus La/Sm diagram for the Sabalan volcanic rocks with melting curves obtained using the nonmodal batch melting modelling, modified after [26].


## 5.2. Geodynamic Setting

Plio-Quaternary (to Quaternary) volcanism in the Peri-Arabian region (in Turkish-Iranian high plateau) postdates the subduction of Neotethyan lithosphere since ~13 Ma [2]. These volcanic rocks become more alkaline with time and towards the south [2, 4, 5, 27]. However, most of volcanic rocks, like the Salmas and/or Maku Quaternary lavas, still show a SSZ geochemical signature including enrichment in LILE and/or depletion in HFSE, despite the lack of subduction system at that time. Geochemical and geodynamic evidence also shows a subduction setting for the genesis of Sabalan lavas: fractionated REEs patterns (high LREEs/HREEs ratio), enrichment in Th, Ba, Rb, U, Pb, and depletion in Nb, Ta, and Ti. Such geochemical characteristics are common in subduction-related environments [25]. These observations and isotopic characteristics of the lavas suggest that their magmas were probably derived from partial melting of a subduction-metasomatized continental lithospheric mantle in the garnet/spinel lherzolite field [2]. Widespread volcanism across the Turkish-Iranian plateau throughout the late Cenozoic until Plio-Quaternary is strictly ascribed to postcollisional mantle upwelling during slab steeping and breakoff beneath a subduction-accretion complex [4, 28].


## 6. Conclusions

The Sabalan Quaternary volcanism includes a sequence of trachyandesites, latites, trachytes with huge bodies of ignimbrites, and pyroclastic rocks. Plagioclase with andesine to labradorite and clinopyroxene with augitic composition are predominant phases in these volcanic rocks. Lavas show shoshonitic geochemical signature with a fractioned REEs trend, enrichment in LREEs relative to HREEs, and depletion in HFSEs. The Sabalan volcanic rocks are isotopically characterized by derivation from an enriched mantle source (with/without crustal influence). All samples show tendency to plot in a field defined by island-arc (IAB) and ocean-island (OIB) basalts, in terms of theirεNd value and 87Sr/86Sr ratio. The lavas are characterized by higher Th/Yb ratio, resembling those lavas that erupted above an active continental margin with geochemical evolution via assimilation-fractional crystallization (AFC) process. The geochemical and isotopic signatures of the Sabalan lavas suggest that their melts have been issued via low degree partial melting of a subduction-metasomatized continental lithospheric mantle during Neotethyan slab steeping and break-off beneath the Turkish-Iranian high plateau.



---

*Source: 735498-2013-12-07.xml*