Petrology and petrogenesis of a tuff layer in the Upper Red Formation, north of Tabriz (NW Iran): shoshonitic magmatism with high Ba and Sr affinity

Document Type : Original Article

Authors

1 Assistant Professor, Faculty of Mining Engineering, Sahand University of Technology, Sahand New Town, Iran

2 Assistant Professor, Department of Geology, Payame Noor University, P.O. Box 19395-4697 Tehran, Iran

Abstract

Introduction
The studied volcanic rocks as a dacitic tuff layer intercalated with the Upper Red Formation (URF, Late Miocene) and are located in the vicinity of the North Tabriz fault. During the Neogene, the red highlands north of Tabriz fault (Eynali) were a different basin from its southern part, namely the Sahand Volcanic complex. Based on the studies carried out on the volcanic rocks and pyroclastics of Sahand Volcano, the volcanic centers of Sahand have been active intermittently from the Late Miocene to the Late Pleistocene (Ghauori, 2002; Ghalamghash et al., 2019). The Upper Red Formation consists of red conglomerate alternated with sandstone, shale and marl and is associated with evaporite units (Asadian, 1993). These sediments have been deposited following uplift in a back-arc basin and within the Neotethys volcanic arc in Central Iran (Shahabpour, 2007). The aim of the present study is to investigate the relation between the tuff layer of URF and the first volcanic manifestations of Sahand volcano.
Geological Background
Subduction of the Neotethys under the central Iranian plate, followed by the collision of the Iranian and the Arabian plates (continental-continental collision), is responsible for the development of four structural zones in Iran. These structural zones with northwest-southeast trend include Zagros-Folded-Thrust belt, Sanandaj-Sirjan metamorphic and magmatic zone and Urmia-Dokhter magmatic arc (Alavi, 1994; Mohajjel et al., 2003). Omrani et al. (2008) have divided the volcanic rocks of Urmia-Dokhtar magmatic arc (including the studied area) into two categories: Eocene and Miocene to Plio-Quaternary. Eocene volcanic rocks consist of andesite, tuff and intermediate pyroclastics with small amounts of basalt, andesite and rhyolite. Miocene to Plio-Quaternary volcanic rocks are composed of andesite to dacititc rocks with Late Miocene to Pliocene age, which are followed by mafic volcanic rocks (Jahangiri, 2007; Omrani et al., 2008). Dacitic domes belonging to Late Miocene in the north of Tabriz fault, with adakitic composition, intruded the Upper Red Formation or Eocene volcanic units (Jahangiri, 2007).
Analytical Methods
Due to the lack of textural and mineralogical diversity of the studied rocks, four fresh samples were sent to the laboratory of the SGS Company located in Toronto, Canada, for analysis of major, trace and rare earth elements with ICP-MS. In order to determine the chemical composition of the rock-forming minerals, a sample of the studied rocks, after preparing a thin-polished section, was analyzed with an electron microprobe (CAMECA SX100) device at the Mineral Processing Research Center of Iran. The analytical conditions for voltage, beam current and beam diameter were set to 15 kV, 20 nA and 5μm, respectively.
Discussion

Petrography

The studied rocks are dominated by the presence of quartz, plagioclase, alkali feldspar and biotite as phenocrysts with a glassy groundmass (Hyaloporphyry) Apatite is rare and calcite and iron oxides form the secondary minerals. Quartz as anhedral to subhedral with embayed texture accounts for about 15% of phenocrysts. Plagioclase is subhedral and forms for about 25% of the rock phenocrysts.

Mineral Chemistry

Plagioclase and potassium feldspar are Ab70An25Or5 and Ab32An1Or67, in composition respectively. Thermometry of the feldspars based on the Ab-An-Or diagram (Fuhrman and Lindsley, 1988; Nekvasil, 1992) shows that they are of relatively low temperature type (~700 ºC). The composition of micas varies from biotite to phlogopite in diagram of Fe/Fe+Mg vs. total Al and are classified as primary and re-equilibrated primary biotites on [(Fe*+Mn)-10*TiO2-MgO] diagram. The studied biotites belong to calc-alkaline orogenic suites originated from a crust-mantle mixed source.

Whole-Rock Geochemistry

The studied rocks have a distinct enrichment of LILE (i.e. Rb, Ba, Th, U, K) and LREE compared to HFSE (i.e. Ta, Nb, Ti, Zr, Hf, Y) and HREE. The rocks have high amounts of Sr (400-540 ppm) and Ba (930-1130 ppm) as well. The studied tuff indicates the features of metaluminous and high-K calc-alkaline to shoshonite magmatic suites, and has the characteristics of rare elements indicative of arc type magmatism. The Nb/Ta ratio in the studied samples varies from 14.7 to 15.8, which is higher than the predicted values for the continental (Taylor and Mclennan, 1985), but it is similar to arc volcanic rocks (Stolz et al., 1996). The above features in combination with the negative anomaly of Nb, Ta and Ti and the high ratios of Ba/La, Ba/Zr and Ba/Nb >30 (Gill, 1981) point to their similarity with magmas related to subduction.
Discussion and Conclusion
The lack of geological evidence in the region, indicating the existence of active subduction at the time of formation of the rocks under study; Thus, the observed geochemical features seem to be related to the origin rather than a tectonic origin. The enrichment of the studied rocks with some elements (i.e. Ba, Sr and Rb) requires extensive crystallization, crust contamination or very small partial melting. The studied rocks show non-adakitic characteristics, and therefore their genesis may be different from the types of adakitic rocks of Sahand Volcanic Complex. The bedrock of Sahand volcano is composed of Paleozoic-Mesozoic sedimentary deposits, Eocene volcanic rocks, lower Miocene deposits (Qom Formation) and Upper Red Formation (including marls, and red sandstones and gypsum belonging to the middle to late Miocene) (Abbassi et al., 2021). On the other hand, the oldest activity of Sahand volcano is attributed to the Late Miocene (Old Sahand in the division of Ghalamghash et al. (2019) with an age of ~8 Ma; and the thick pyroclastic sequence on the western slope of the volcano named Ghermeziqul Formation in the division of Moine Vaziri and Amine Sobhani (1977) with an age of 9-12 Ma). Therefore, considering the stratigraphic position of the tuff layer and its geochemical similarities with the non-adakitic eruptions of Sahand, it is likely the tuff layer was originated as the result of the first explosive activity of Sahand at the same time with the formation of Upper Red sediments (Late Miocene).

Keywords

Main Subjects


مقدمه

منطقة آذربایجان بخشی از پهنة ماگمایی البرز باختری-آذربایجان یا کمان ماگمایی تبریز-بزمان، در ‌شمال‌باختری ایران و در مقیاس بزرگ‌تر در پهنة برخوردی میان صفحة عربستان و اوراسیاست. پس از برخورد دو صفحة عربستان و اوراسیا و رخداد ستبرشدگی پوسته‌ای و بالاآمدگی، دور جدیدی از فعالیت‌های آتشفشانی در این منطقه رخ داده است که احتمالاً سنی نزدیک به 8/2 تا 12 میلیون سال پیش داشته است (Keskin et al., 1998).

بر پایة بررسی‌ سنگ‌های آتشفشانی و آذرآواری آتشفشان سهند و با توجه به داده‌های گاه‌شماری و روابط چینه‌شناسی حاکم بر سنگ‌های آتشفشانی سهند، مراکز آتشفشانی سهند از میوسن پسین تا پلیستوسن پایانی به‌طور متناوب در فعالیت بوده‌اند (Ghauori, 2002). در مراحل آغازین فعالیت آتشفشانی در سهند، گدازه‌های با ترکیب مختلف فوران کرده‌اند و در پایان میوسن، انتشار گستردة پونس و خاکستر آتشفشانی رخ داده است. در جدیدترین فعالیت آتشفشانی، پس از مراحل آرامش و خاموشی، گدازه‌های اسیدی مانند داسیت و ریوداسیت فوران کرده‌اند که همسن با گنبد‌های نیمه‌آتشفشانی دیگر در منطقة آذربایجان بوده‌اند و به مرحلة پایانی فعالیت آتشفشانی وابسته هستند (Ghauori, 2002). به باور جهانگیری و اشرفی (Jahangiri and Ashrafi, 2006)، وجود لایه‌های رسی همراه با مواد عارضه‌ای در بخش زیرین قدیمی‌ترین سکانس‌های آذرآواری سهند، نشانه‌ای از فعالیت فراتوماگمایی این آتشفشان در آغاز فعالیت آن به‌شمار می‌رود.

داده‌های سن‌سنجی به روش پتاسیم-آرگونِ معین وزیری و امین سبحانی (Moine Vaziri and Amine Sobhani, 1977) نشان‌دهندة فعالیت آتشفشان سهند از 12 میلیون سال پیش تا 400-140 هزار سال پیش هستند. قلمقاش و همکاران (Ghalamghash et al., 2019) بر پایة ویژگی‌های صحرایی، سنگ‌نگاری و زمین‌شیمیایی، آتشفشان سهند را به سه مرحلة سهند قدیمی، سهند میانی و سهند جوان دسته‌بندی کرده‌اند. محصولات فعالیت سهند در مراحل قدیمی و میانی شامل مواد آذرآواری و گدازه و سهند جوان ‌به‌صورت گنبدهایی با ترکیب حد واسط تا اسیدی بوده است. ترکیب سنگ‌شناسی سنگ‌های آتشفشانی سهند از آندزیت بازالتی تا آندزیت، داسیت و لاتیت متغیر است و سنگ‌ها به سری ماگمایی آهکی-قلیایی، آهکی-قلیایی پتاسیم بالا، آداکیتی و شوشونیتی (قلیایی پتاسیک) گرایش دارند (Moine Vaziri, 1999; Pirmohammadi et al., 2011; Chaharlang et al., 2016; Lechmann et al., 2018; Ghalamghash et al., 2019). به باور معین وزیری و امین سبحانی (Moine Vaziri and Amine Sobhani, 1977)، برخلاف مراحل آغازین، در پایان فعالیت آتشفشان سهند ماگمای غنی از K، Ba و Sr تولید شده است و ماگما‌های سازندة آتشفشان خاستگاه‌های متفاوتی داشته‌اند. با این وجود، بررسی‌های جدید نشان می‌دهد ماگما‌های غنی از K (شوشونیتی)، Ba و Sr در مراحل آغازین فعالیت سهند (سهند قدیمی) و تقریباً با فوران‌های آتشفشان سارای نیز همزمان تولید شده است (Lechmann et al., 2018; Ghalamghash et al., 2019).

گاه سنگ‌های سری شوشونیتی با قلیایی پتاسیک مترادف گرفته شده است (Muller and Groves, 1997). با وجود این، سنگ‌های غنی از پتاسیم مجموعة‏ شوشونیتی معمولاً با اشباع از سیلیس‌بودن شناخته می‏شوند؛ بدین‌گونه‌که معمولاً هیپرستن نورماتیو دارند؛ اما اندکی نفلین یا کوارتز نورماتیو دارند یا کلاً ندارند (Hyndman, 1985). مقدار کل آلکالی (Na2O+K2O) در شوشونیت‌ها، از 5 درصدوزنی بیشتر است. همچنین، نسبت K2O به Na2O آنها در سنگ‌های با سیلیس 55 درصدوزنی از یک بزرگ‌تر است و در سنگ‌های با درصد سیلیس 50 درصدوزنی بیشتر از 6/0 است. محتوی TiO2 آنها از 3/1درصدوزنی کمتر است (Morrison, 1980). سنگ‌های غنی از پتاسیم واقعی، بسیار تحت‏اشباع از سیلیس هستند و لوسیت نورماتیو چشمگیری دارند. چنین سنگ‌هایی ممکن است از لحاظ خاستگاه نسبت به سنگ‌های تقریباً اشباع از سیلیس تفاوت داشته باشند و شاید در یک گروه جداگانه باشند (Ferguson and Cundari, 1975; Morrison, 1980).

الگوهای سنگ‌زایی پیدایش ماگماهای شوشونیتی بسیار گسترده هستند. این الگو‌ها توصیف‌کنندة ترکیب متغیر و محیط زمین‌ساختی متفاوت این نوع ماگما‌ها هستند. ماگماهای شوشونیتی نشان‌دهندة مراحل پایانی میان فعالیت ماگمایی کلسیمی- قلیایی و قلیایی در طول تکامل محیط پس از برخورد و تأخیری هستند (Liegeois et al., 1998). بیشتر ماگماهای پتاسیک و التراپتاسیک در یکی از این سه محیط یافت می‌شوند:

1) در بالای پهنة فرورانش، در بخش ژرف این پهنه (Morrison, 1980

2) در محیط‌های کوهزایی پس از برخورد مانند تبت (Wang et al., 2006) و آلپ (Venturelli et al., 1984

3) به‌طور نادر در محیط‌های کششی دورن‌قاره‌ای مانند بخش جنوب‌باختری اوگاندا و کافت خاور آفریقا (Thompson, 1985).

منطقة مورد بررسی در نزدیکی گسل شمال تبریز و در ارتفاعات عون‌ابن‌علی (عینالی) مشرف به شهر تبریز جای دارد و سنگ‌های آتشفشانی آن به شکل لایة توف داسیتی درون سازند قرمز ‌بالایی یا URF [1](میوسن پسین) یافت می‌شوند که پهنه‌ای متفاوت از بخش جنوبی گسل تبریز یعنی مجموعة آتشفشانی سهند در طول نئوژن است (شکل 1). با توجه به نبود مطالعات پترولوژی روی این سنگ‌ها و بررسی این موضوع که آیا این ماگماتیسم در ارتباط با فعالیت آتشفشان سهند بوده است یا نه پژوهش کنونی انجام شد.

روش انجام پژوهش

پس از بررسی‌های صحرایی و مشخص‌شدن روابط سنی و صحرایی سنگ‌های گوناگون منطقه، شمار 20 نمونه از سنگ‌های آتشفشانی منطقه گردآوری شد. پس از تهیة مقاطع نازک از آنها، ‌بررسی میکروسکوپی سنگ‌ها و تعیین ویژگی‌های بافتی و کانی‌شناسی آغاز شد. ‌برای تعیین دقیق ترکیب شیمیایی کانی‌های سازنده، نمونه‌ای از سنگ‌ها پس از تهیة مقطع نازک- صیقلی با دستگاه ریزکاو الکترونی مدل CAMECA SX100 در مرکز تحقیقات فرآوری مواد معدنی ایران تجزیة شیمیایی نقطه‌ای شد (جدول‌های 1 و 2). تجزیة روی کانی‌های یادشده در شرایط ولتاژ 15کیلوولت، شدت جریان20 نانوآمپر و قطر پرتو 5 میکرون انجام گرفت. برای جدایش مقادیر آهن دو و سه ظرفیتی بیوتیت از برنامة MICA+ (Yavuz, 2003) استفاده شد. با توجه به نبود تنوع بافتی و کانی‌شناسی سنگ‌های منطقه، شمار 4 نمونه سالم از سنگ‌های منطقه ‌برای تجزیة عنصرهای اصلی، فرعی و ‌خاکی کمیاب با دستگاه ICP-MS، به آزمایشگاه شرکت SGS واقع در شهر تورنتو کشور کانادا فرستاده شدند. داده‌های به‌دست‌آمده از تجزیة نمونه‌ها در جدول 3 آورده شده‌اند. افزون‌بر داده‌های این ‌بررسی، ‌برای مقایسه و داشتن دید کلی از ‌ویژگی‌های زمین‌شیمیایی ماگماتیسم‌های مشابه در منطقه به‌ویژه مجموعة آتشفشانی سهند از داده‌های بررسی‌های پیشین در اینباره کمک گرفته شد. در مرحلة پایانی این پژوهش، داده‌های به‌دست‌آمده از بررسی‌های صحرایی، سنگ‌نگاری و زمین‌شیمی با یکدیگر تلفیق و تجزیه و تحلیل شدند. در پردازش داده‌های سنگ کل نرم‌افزار GCDkit 4.1 (Janousek et al., 2006) به‌کار برده شد.

شکل 1. نقشة زمین‌شناسی مجموعة آتشفشانی سهند در جنوب گسل تبریز، جایگاه منطقة بررسی‌شده و سازند قرمز ‌بالایی به سن میوسن پسین در شمال شهر تبریز (با تغییرات و ساده‌شده از مجموعه نقشه‌های 100000/1 زمین‌شناسی در محدودةآتشفشان سهند: تبریز (Asadian, 1993)، بستان‌آباد (Behrouzi et al., 1997)، اسکو (Khodabandeh and Amini Fazl, 1995)، آذرشهر (Ghadirzadeh, 2002)، قره‌آغاج (Rezaei and Fonoudi, 2009)، مراغه (Alavi and Shahrabi, 1975) و عجب‌شیر (Soltani Sisi, 2005). نقشة ایندکس، جایگاه منطقه مورد بررسی در ‌شمال‌باختری ایران را نشان می‌دهد.

Figure 1. Geological map of the Sahand Volcanic Complex in the south of the Tabriz fault along with the location of the studied area and the Upper Red Formation of Late Miocene age in the north of Tabriz (modified and simplified from 1:100000 geological quadrangle maps of Sahand Volcano region: Tabriz (Asadian, 1993), Bostan Abad (Behrouzi et al., 1997), Osku (Khodabandeh and Amini Fazl, 1995), Azar Shar (Ghadirzadeh, 2002), Ghareh Aghaj (Rezaei and Fonoudi, 2009), Maragheh (Alavi and Shahrabi,1975), and Ajab Shir (Soltani Sisi, 2005). The index map shows the location of the studied area in the northwest of Iran.

جدول 1. داده‌های تجزیه ریزکاوالکترونی کانی پلاژیوکلاز (بر پایة درصد وزنی) به‌همراه فرمول ساختاری به‌دست آمده برپایه 8 اتم اکسیژن (سازنده‌های پایانی: Ab: آلبیت؛ Or: ارتوکلاز و An: آنورتیت).

Table 1. Electron microprobe analysis of plagioclase mineral (in wt.%) along with the structural formula based on 8 oxygen atoms (the end members: Ab: albite, Or: orthoclase and An: anorthite).

Sample

kfs41

kfs42

kfs43

plg33

plg34

plg35

plg36

plg37

plg38

SiO2

61.81

63.93

63.59

57.97

58.84

57.42

58.35

59.60

59.47

TiO2

3.50

0.00

0.00

0.00

0.44

3.16

0.00

1.16

0.00

Al2O3

18.35

18.95

19.57

25.35

24.78

24.78

24.74

23.75

23.88

FeO

0.08

0.06

0.12

0.18

0.23

0.29

0.18

0.18

0.18

MgO

0.00

0.00

0.01

0.01

0.01

0.01

0.02

0.02

0.01

CaO

0.23

0.33

0.14

5.47

5.14

6.00

5.76

4.53

5.00

Na2O

3.63

3.75

3.61

8.62

8.70

7.81

8.81

8.89

8.79

K2O

11.60

11.83

11.92

0.95

0.89

0.75

1.07

1.08

0.95

Total

99.20

98.85

98.96

98.55

99.03

100.22

98.93

99.21

98.28

Si

5.93

5.93

5.87

5.28

5.35

5.30

5.34

5.45

5.43

Al

2.07

2.07

2.13

2.72

2.65

2.70

2.66

2.56

2.57

Ti

0.25

0.00

0.00

0.00

0.03

0.22

0.00

0.08

0.00

Fe2+

0.01

0.01

0.01

0.01

0.02

0.02

0.01

0.01

0.01

Mg

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Ca

0.02

0.03

0.01

0.53

0.50

0.59

0.56

0.44

0.49

Na

0.68

0.67

0.65

1.52

1.53

1.40

1.56

1.58

1.56

K

1.42

1.40

1.40

0.11

0.10

0.09

0.13

0.13

0.11

Cations

10.38

10.11

10.07

10.18

10.19

10.32

10.27

10.24

10.17

X

8.25

8.00

8.00

8.00

8.03

8.22

8.00

8.08

8.00

Z

2.12

2.11

2.07

2.18

2.16

2.10

2.27

2.16

2.17

Ab

31.90

32.00

31.30

70.30

71.70

67.20

69.40

73.50

72.20

An

1.10

1.60

0.70

24.70

23.40

28.50

25.10

20.70

22.70

Or

67.00

66.40

68.00

5.10

4.80

4.20

5.60

5.90

5.10

 جدول 2. داده‌های تجزیه ریزکاوالکترونی از کانی بیوتیت به‌همراه فرمول ساختاری بر پایه 22 اتم اکسیژن.

Table 2. Electron microprobe analysis of biotite mineral along with the structural formula based on 22 oxygen atoms.

Sample

bt30

bt31

bt32

bt39

bt40

bt44

bt45

SiO2

37.65

38.38

38.51

38.84

38.27

38.55

39.16

TiO2

4.51

5.16

3.07

2.58

3.50

4.04

3.24

Al2O3

13.07

13.65

12.97

14.71

15.77

14.76

14.87

Cr2O3

0.00

0.00

0.00

0.00

0.03

0.00

0.00

V2O3

0.20

0.21

0.23

0.10

0.19

0.23

0.25

FeO

13.86

12.78

13.96

14.04

14.45

13.97

13.71

MnO

0.23

0.25

0.25

0.23

0.22

0.26

0.26

MgO

16.09

17.13

16.53

13.00

13.82

13.41

13.88

CaO

0.09

0.00

0.00

1.60

0.05

0.00

0.00

Na2O

0.61

0.63

0.57

0.57

0.61

0.65

0.67

K2O

9.08

8.37

9.29

9.47

9.05

9.34

9.33

Total

95.39

96.56

95.38

95.14

95.96

95.21

95.37

Si

5.72

5.67

5.82

6.08

5.80

5.94

5.99

AlIV

2.28

2.33

2.18

1.92

2.20

2.06

2.01

AlVI

0.05

0.05

0.13

0.79

0.62

0.62

0.67

Ti

0.52

0.57

0.35

0.30

0.40

0.47

0.37

Fe3+IV

0.08

0.06

0.00

0.00

0.00

0.00

0.00

Fe3+VI

0.73

1.07

0.47

1.08

0.08

0.04

0.25

Fe2+

0.95

0.45

1.30

0.76

1.75

1.76

1.51

Mn

0.03

0.03

0.03

0.03

0.03

0.03

0.03

Mg

3.64

3.77

3.73

3.03

3.12

3.08

3.17

Ca

0.02

0.00

0.00

0.27

0.01

0.00

0.00

Na

0.18

0.18

0.17

0.17

0.18

0.19

0.20

K

1.76

1.58

1.79

1.89

1.75

1.84

1.82

Cations

15.96

15.76

15.96

16.33

15.94

16.03

16.02

Mg#

0.67

0.70

0.68

0.62

0.63

0.63

0.64

 جدول 3. داده‌های به‌دست‌آمده از تجزیة شیمیایی سنگ کل نمونه‌های معرف توف داسیتی.

Table 3. The whole rock chemical data of the representative samples of dacitic tuff.

Sample No.

TM-3

TM-4

TM-9

TM-12

SiO2

64.1

64.1

63.7

59.1

TiO2

0.37

0.36

0.4

0.37

Al2O3

14.4

14.6

14.8

14.2

Fe2O3t

3.34

3.46

3.48

2.75

MnO

0.06

0.07

0.08

0.04

MgO

0.89

1.01

1.07

2.31

CaO

2.49

3.19

2.87

2.91

Na2O

3.3

3.7

3.4

3.2

K2O

4.96

5.01

4.95

1.28

P2O5

0.19

0.17

0.19

0.17

LOI

1.77

2.5

2.09

8.69

Sum

95.9

98.3

97

95.1

Ba

1090

1130

1040

930

Mg#

35

37

38

62

K2O/Na2O

1.50

1.35

1.45

0.4

Sr

540

430

460

400

Ce

98.5

93.7

100

92.6

Cs

1.8

1.6

1.8

2.3

Dy

3.01

3.37

3.37

2.99

Er

1.66

1.79

1.98

1.59

Eu

1.06

0.92

1.05

1

Ga

17

17

18

17

Gd

3.81

3.93

4.09

3.78

Hf

5

5

5

4

Ho

0.6

0.66

0.7

0.58

La

60

54.8

59.5

54.7

Lu

0.29

0.29

0.32

0.23

Nb

25

25

27

24

Nd

31.6

30.4

32.5

29.9

Pr

9.91

9.43

9.97

9.29

Rb

123

124

128

35.4

Sm

5.3

5.1

5.5

4.9

Ta

1.7

1.7

1.7

1.6

Tb

0.55

0.59

0.65

0.55

Th

39.9

38.9

41.2

36.6

Tm

0.25

0.28

0.29

0.24

U

6.6

6.82

7.1

5.77

V

54

44

48

58

Y

17.2

19.9

20.5

16.3

Yb

1.8

1.8

2.1

1.6

Zr

172

172

179

137

Sr/Y

31.39

21.61

22.44

24.53

Eu/Eu*

0.72

0.63

0.68

0.71

LaN/YbN

22.47

20.53

19.10

23.05

زمین‌شناسی

پهنة نئوژن تبریز در ‌شمال‌باختری ایران بخشی از یک سیستم گسلی راستالغز در میوسن پایانی میان صفحه‌های برخوردیِ عربی و اوراسیا است (Axen et al., 2001 Allen et al., 2004;). امروزه پهنة تبریز یک پهنة میان کوهستانی به‌شمار می‌رود که از شمال و شمال‌خاوری با رسوبات آواری خشکی سازند قرمز ‌بالایی (کوه عون‌ابن‌علی یا عینالی)، به سن میوسن میانی تا پسین و در جنوب با آتشفشان سهند (12میلیون تا 140 هزارسال پیش) (Gansser, 1955; Fisher, 1968; Moine Vaziri and Amine Sobhani, 1977; Moine Vaziri, 1999) فراگرفته شده است. از سمت باختر، این پهنه به زمین‌های پست دریاچة ارومیه می‌رسد.

آتشفشان سهند از آتشفشان‌های جوان در بخش شمالی کمربند ماگمایی ارومیه-دختر است که از میوسن پایانی تا پلیوسن ‌به‌صورت متناوب فعالیت‌های انفجاری و خروج گدازه داشته است. گدازه‌های آتشفشانی و نیمه‌آتشفشانی اسیدی و تا اندازه‌ای حد واسط در این منطقه گسترش چشمگیری دارند. به‌گونه‌ای‌که می‌توان کامل‌ترین پیوندها با فعالیت آتشفشانی پس از برخورد از میوسن تا پلیوسن را در این منطقه دید (Moine Vaziri and Amine Sobhani, 1977).

فرورانش پوستة اقیانوسی نئوتتیس به زیر ایران مرکزی و به‌دنبال آن برخورد ورقه ایران و ورقه عربی (برخورد قاره-قاره‌ای) گسترش چهار پهنه ساختاری در ایران را در پی داشته است. این پهنه‌های ساختاری با روند ‌شمال‌باختری-جنوب‌خاوری شامل پهنة چین‌خوردة زاگرس (زاگرس مرتفع و کمربند چین‌خوردة زاگرس)، پهنة‌ دگرگونی سنندج-سیرجان و کمان ماگمایی ارومیه-دختر هستند ( Alavi, 1994; Mohajjel et al., 2003). کمان ماگمایی ارومیه-دختر از سنگ‌های آذرین درونی و بیرونی با بازة سنی ائوسن-کواترنری و راستای NW-SE ساخته شده است و بر پایة بررسی‌های زمین‌شیمیایی بیشتر شامل سنگ‌های آهکی-قلیایی پهنة فرورانش است (Jahangiri, 2007). این پهنه بخشی از فعالیت‌های آتشفشانی- درونی ترسیری در ایران مرکزی به‌شمار می‌رود. فعالیت‌های آتشفشانی این کمربند از کرتاسة پسین با فوران‌های آهکی-قلیایی آغاز شده است و در ائوسن (زمان اوج فعالیت) با گدازه‌های آهکی- قلیایی ادامه یافته است و در پایان ترشیری، گدازه‌های غنی از پتاسیم پدید آمده‌اند (Berberian and King, 1981). عمرانی و همکاران (Omrani et al., 2008) سنگ‌های آتشفشانی این کمان ماگمایی را به دو گروهِ سنگ‌های آتشفشانی ائوسن و سنگ‌های آتشفشانی میوسن تا پلیوکواترنری دسته‌بندی کرده‌اند. سنگ‌های آتشفشانی ائوسن متشکل از آندزیت، توف و آذرآواری‌های حد واسط با مقادیر ناچیزی بازالت، آندزیت و ریولیت هستند. سنگ‌های آتشفشانی میوسن تا پلیوکواترنری دربردارندة سنگ‌های آندزیتی تا داسیتی با سن میوسن پسین تا پلیوسن هستند که با سنگ‌های آتشفشانی مافیک دنبال می‌شوند (Jahangiri, 2007; Omrani et al., 2008). دم‌های داسیتی با سن میوسن پسین در شمال گسل تبریز که ترکیب آداکیتی دارند درون سازند قرمز ‌بالایی و یا واحد‌های آتشفشانی ائوسن تزریق شده‌اند (Jahangiri, 2007).

لایة توف داسیتیِ بررسی‌شده در شمال شهر تبریز و در کوه عون‌ابن‌علی (عینالی) رخنمون دارد. سن‌سنجی این توف‌ها به روش اثر شکافت آپاتیت، سن 1/1± 2/11 میلیون سال پیش (تورتونین) را نشان می‌دهد که با پیدایش سازند قرمز ‌بالایی در این منطقه همزمان است (Reichenbacher et al., 2011). سازند قرمز ‌بالایی در بخش بزرگی از مرکز و ‌شمال‌باختری ایران گسترش دارد و بخش بالایی سازند قم دانسته می‌شود. این سازند تناوبی از کنگلومرا، ماسه‌سنگ، شیل و مارن سرخ رنگ است و واحد‌های تبخیری آن را همراهی می‌کنند (Asadian, 1993). این رسوبات به‌دنبال بالاآمدگی در پهنة پشت کمان و درون کمانی آتشفشانی نئوتتیس در ایران مرکزی نهشته شده اند (Shahabpour, 2007).

این سنگ‌ها در صحرا به شکل لایه‌ای و روشن رنگ با ضخامت متغیر (نزدیک به 5 تا 50 متر) هستند که موازی و درون لایه‌هایی از مارن و ماسه‌سنگ‌های سرخ رنگ (سازند قرمز ‌بالایی) گسترش یافته‌اند (شکل‌های 2-D و 2-G). از دیدگاه گسترش جانبی، گسترش آن ممتد نیست و گاهی با گسل‌های فرعی جابجا شده است. مرزبالایی لایة یادشده ‌به‌صورت گسله با سازند قرمز ‌بالایی (شکل 2-A) یا پوشیده شده با آن بدون سطح فرسایشی دیده می‌شود. مرز زیرین آن در همبری مستقیم با سازند قرمز بالایی است و تغییر تدریجی میان آنها دیده نمی‌شود.

بخش بالاییِ لایة توفی کمابیش سخت و سنگی است (توف بلورین) و با رنگ ظاهری خاکستری مایل به سرخ و رنگ خاکستری بسیار روشن در سطح تازه دیده می‌شود (شکل‌های 2-E، 2-F و2-I)؛ اما بخش زیرین به رنگ خاکستری روشن است و خاکستر ریز تا درشت لاپیلی‌های هسته‌دار (کوچک‌تر از 10 سانتیمتر) نیز در آن یافت می‌شوند (شکل 2-H). انکلاو‌های مافیک در اندازه‌های کمتر از 5 سانتیمتر معمولاً درون لایة توفی دیده می‌شوند (شکل 2-J). سطح تازه این سنگ‌ها در نمونة دستی رنگ خاکستری روشن دارد و بلور‌های درشت میکا (بیوتیت) در آن را می‌توان با چشم غیر مسلح به‌آسانی شناسایی کرد (شکل 3-A). برخی نمونه‌های این سنگ‌ها دچار دگرسانی شدیدی شده‌اند و کانی‌های رسی، کلسیت و کوارتز در آنها پدید آمده است.

سنگ‌نگاری

همان‌گونه‌که در بخش پیشین گفته شد، لایة توفی از دو بخش شامل بخش سست و غیرسنگی، یعنی لاپیلی توف (دارای خاکستر ریز و لاپیلی‌های هسته‌دار) و بخش سنگی، یعنی توف بلورین داسیتی ساخته شده است. ازاین‌رو، در این بخش تنها به سنگ‌نگاری بخش سنگی آن پرداخته می‌شود. بافت این سنگ‌ها پورفیری با زمینة‌ شیشه‌ای (ویتروفیری) است و کوارتز، پلاژیوکلاز، آلکالی‌فلدسپار و بیوتیت از درشت بلور‌های آن به‌شمار می‌روند (شکل 3-B).

شکل 2. ویژگی‌های صحرایی لایة توفیِ ‌بررسی‌شده. A) ترانشة کنار جاده‌ای با رخنمونی از لایة توفی، ماسه‌سنگ و مارن‌های سرخ رنگ سازند قرمز بالایی (URF) و کنگلومرای پلیوسن که با همبری گسلی شناخته می‌شوند؛ B) تصویر لایة توفی از نزدیک؛ C) همبری سازند قرمز بالایی و کنگلومرای پلیوسن از نمای نزدیک؛ D) رخنمون لایة روشن رنگ توفی در میان سازند قرمز بالایی در یکی از ورودی‌های تفرجگاه کوه عینالی؛ E، F) لایة توفی که بخش بالایی آن با سطح هوازده مایل به قرمز و بخش زیرین آن با رنگ خاکستری روشن دیده می‌شود؛ G) لایة توفی روشن رنگ موازی لایه‌بندی و درون سازند قرمز بالایی؛ H) لاپیلی‌های هسته‌دار در لایة توفی؛ I) سطح هوازده و تازه لایة توفی؛ J) انکلاو‌های مافیک (کوچک‌تر از 5 سانتیمتر) درون سنگ‌های توفی.

Figure 2. The field characteristics of the studied tuff layer. A) Roadside trench with an outcrop of the tuff layer along with red sandstone and marl of the Upper Red Formation (URF) and Pliocene conglomerate, which are characterized by a fault contact; B) A close-up photo of the tuff layer; C) The contact of URF and the Pliocene conglomerate from a close-up view; D) The outcrop of a light color tuff as intercalation in URF, Eynali Mountain resort; E, F) The tuff layer whose the upper part is characterized by reddish weathered surface and its lower part by light gray color; G) The light-colored tuff layer parallel to the bedding within the URF; H) Armored lapilli in the tuff layer; I) The weathered and fresh surface of the tuff; J) Mafic enclaves (smaller than 5 cm) within the tuff.

شکل 3. تصویر نمونة توف بلورین مورد ‌بررسی در نمونة دستی و در زیر میکروسکوپ. A) نمونة دستی خاکستری روشن توف بلورین داسیتی با بلور‌های سیاه رنگ و ورقه‌ای بیوتیت؛ B) تصویر میکروسکوپی از توف بلورین داسیتی با بافت ویتروفیری و فنوکریست‌های شکل‌دار بیوتیت و بی‌شکل و یا شکستة پلاژیوکلاز (در نور پلاریزه متقاطع).

Figure 3. The hand-specimen and photomicrograph of the crystal tuff. A) Light gray hand-specimen of the dacitic crystal tuff with black and flake biotite crystals; B) The microscopic image of the dacitic crystal tuff with vitrophyric texture where phenocrysts of euhedral biotite and anhedral/crushed plagioclase are observed (in cross-polarized light).

آپاتیت کانی کمیاب و کلسیت و اکسید‌های آهن از کانی‌های ثانویة این سنگ‌ها هستند. کوارتز ‌به‌صورت بی‌شکل تا نیمه‌شکل‌دار گاهی بافت خلیجی دارد و نزدیک به 15 درصد فنوکریست‌ها را در بر می‌گیرد. پلاژیوکلاز ‌به‌صورت نیمه‌شکل‌دار و شکسته است و نزدیک به 25 درصد فنوکریست‌ها را در بر می‌گیرد. بر پایة تجزیة ریزکاو الکترونی ترکیب آن الیگوکلاز شناسایی شد (شکل 4-A). آلکالی‌فلدسپار ‌به‌صورت ‌نیمه‌شکل‌دار است و نزدیک به 10 درصد حجمی فنوکریست‌ها را تشکیل می‌دهد و بر پایة تجزیة ریزکاو الکترونی، ترکیب آن سانیدین شناسایی شد (شکل 4-A). بیوتیت ‌به‌صورت شکل‌دار تا ‌نیمه‌شکل‌دار است و نزدیک به 15 درصدحجمی فنوکریست‌ها را در بر می‌گیرد. بر پایة تجزیه میکروپروب ترکیب آن بیوتیت منیزیم‌دار تا فلوگوپیت به‌دست آمد (شکل 5).

شیمی‌کانی در ارتباط با شرایط فیزیکوشیمیایی ماگما

ترکیب پلاژیوکلازهای توف داسیتی ‌به‌صورت Ab70An25Or5 و ترکیب پتاسیم‌فلدسپارهای تجزیه‌شده ‌به‌صورت Ab32An1Or67 است (جدول 1). دماسنجی فلدسپار‌ها بر پایة نمودار Ab-An-Or ( Fuhrman and Lindsley, 1988; Nekvasil, 1992) نشان می‌دهد این کانی‌ها از نوع کمابیش کم دما (نزدیک به 700 درجة سانتیگراد) هستند (شکل 4-B). جدول 2 داده‌های تجزیة ریزکاوالکترونی میکا در نمونة توف داسیتی را نشان می‌دهد که در آن SiO2 برابر با 68/37 تا 16/39 درصدوزنی، Al2O3 برابر با 04/13 تا 87/14 درصدوزنی، MgO برابر با 13 تا 13/17 درصدوزنی و K2O برابر با 37/8 تا 47/9 درصدوزنی است. ترکیب میکا‌ها در نمودار Fe/Fe+Mg در برابر Al کل در گسترة میان بیوتیت و فلوگوپیت جای می‌گیرد (شکل 5) و در نمودار سه‌تایی (Al+Fe3++Ti)-Mg-(Fe2++Mn)، بیشتر نمونه‌ها در محدودة بیوتیت منیزیم‌دار هستند (شکل 6). همچنین، در نمودار سه‌تایی (Fe*+Mn)-10*TiO2-MgO) در مرز میان بیوتیت‌های اولیه و بیوتیت‌های بازتبلوریافته جای می‌گیرند (شکل 7).

شکل 4. A) ترکیب شیمیایی فلدسپار‌ها در نمودار رده‌بندی Ab-An-Or (Deer et al., 1992B) دماسنجی ماگمای سازندة توف داسیتی بر پایة ترکیب فلدسپارها.

Figure 4. A) Chemical composition of feldspars in the Ab-An-Or classification diagram (Deer et al., 1992); B) Thermometry of the magma forming dacitic tuff based on the composition of feldspars.

شکل 5. نامگذاری بیوتیت‌ها در نمودار Al در برابر Fe/(Fe+Mg) (Deer et al., 1992; Rieder et al., 1998) (نمودار B پایینی بخش هاشورخورده در نمودار A را نشان می‌دهد).

Figure 5. The classification of biotite in Al versus Fe/(Fe+Mg) diagram (Deer et al., 1992; Rieder et al., 1998) (Plot B shows the hatched area in plot A).

شکل 6. ترکیب بیوتیت‌ها در نمودار رده‌بندی Mg-(AlVI+Fe3++Ti)-(Fe2++Mn) (Foster, 1960).

Figure 6. Composition of biotites in the classification diagram of Mg-(AlVI+Fe3++Ti)-(Fe2++Mn) (Foster, 1960).

شکل 7. تعیین اولیه یا ثانویه‌بودن ترکیب بیوتیت‌ها در نمودار (FeO*+MnO)-10*TiO2-MgO (Nachit et al., 2005).

Figure 7. The primary and secondary classification of biotites in (FeO*+MnO)-10*TiO2-MgO diagram (Nachit et al., 2005).

برای تعیین سری ماگمایی و محیط تکتونوماگمایی ماگمای سازندة بیوتیت‌ها از نمودارهای پیشنهادیِ عبدالرحمان (Abdel-Rahman, 1994) استفاده شد (شکل 8). بر پایة این نمودارها، بیوتیت‌ها به سری ماگمایی کالک‌آلکالن کوهزایی متعلق هستند. گفتنی است در این نمودار‌ها محدوده‌ای برای سری‌های آلکالن کوهزایی لحاظ نشده است و چنین سری‌هایی با محدودة کالک‌آلکالن کوهزایی همپوشانی نشان می‌دهند (Ashrafi et al., 2009).

شکل 8 . بررسی سرشت ماگما در: A، B، C و E) نمودار اکسیدهای FeO*، MgO و Al2O3 بر پایة ترکیب بیوتیت (Abdel-Rahman, 1994D) نمودار Mg در برابر Al (Nachit et al., 1985).

Figure 8. Determining the magma nature in: A, B, C, E) FeO*, MgO and Al2O3 oxides diagrams based on biotite composition ( (Abdel-Rahman (1994); D) Mg versus Al diagram (Nachit et al., 1985).

همچنین، برای بررسی نقش پوسته و گوشته در پیدایش ماگمای سازندة بیوتیت‌ها نمودار MgO در برابر TFeO/(TFeO+MgO) به‌کار برده شد. بر پایة این نمودار، بیوتیت‌ها در محدودة منبع پوسته‌ای و گوشته‌ای مخلوط‌شده جای می‌گیرند (شکل 9).

شکل 9. ترکیب بیوتیت در نمودار MgO در برابر TFeO/(TFeO+MgO) (Zhou, 1986) (TFeO: آهن کل ‌به‌صورت FeO).

Figure 9. The biotite composition in the TFeO/(TFeO+MgO) versus MgO diagram (Zhou, 1986) (TFeO represents total iron as FeO).

بیوتیت یکی از کانی‌هایی است که میزبان تیتانیم است و غلظت آن با دما کنترل می‌شود. هنری وگوییدوتی (Henry and Guidotti, 2002) و هنری و همکاران (Henry et al., 2005) معادلة زیر را برای به‌دست‌آوردن دمای پیدایش بیوتیت پیشنهاد کرده‌اند.

که در آن 3594/2 -a=، 10-9 x 6482/4 b=، 7283/1 - c= و XMg=Mg/(Mg+Fe) است. بر پایة این معادله، دمای 700 تا800 درجة سانتیگراد برای پیدایش بیوتیت‌ها برآورد شد (شکل 10). ‌برای بررسی فوگاسیتة اکسیژن در ماگما، نمودار سه‌تایی Fe2+-Fe3+-Mg به‌کار برده شد. همان‌گونه‌که در شکل 11 دیده می‌شود، نمونه‌ها گستره‌ای از فوگاسیتة اکسیژن کم تا بالا را نشان می‌دهند.

شکل 10. دماسنجی بر پایة ترکیب بیوتیت با کمک نمودار Mg/(Mg+Fe) در برابر Ti (Henry et al., 2005).

Figure 10. Biotite thermometry in Mg/(Mg+Fe) versus Ti diagram (Henry et al., 2005).

 شکل 11. ارزیابی فوگاسیتة اکسیژن در ماگما بر پایة ترکیب بیوتیت در نمودار سه‌تایی Fe3+–Fe2+–Mg (Wones and Eugster, 1965).

Figure 11. Magma fugacity based on the biotite composition in Fe3+-Fe2+-Mg ternary diagram (Wones and Eugster, 1965).

شیمی سنگ کل

دامنة تغییرات مقادیر SiO2 نمونه‌های مورد ‌بررسی از 1/59 تا1/64 درصدوزنی در نوسان است. با توجه به اینکه این سنگ‌ها کم و بیش دارای نشانه‌هایی از دگرسانی بودند، برای بررسی‌های زمین‌شیمیایی آنها نمودار‌هایی بر پایة عنصرهای کم تحرک به‌کار برده شد. میزان LOI بیشتر از 5/2 درصدوزنی نمونه‌ها را می‌توان نشانة دگرسان‌بودن آنها دانست. مقدار LOI برای نمونه‌های TM3، TM4، TM9 و TM12 به‌ترتیب برابر با 77/1، 5/2، 09/2 و 69/8 درصدوزنی است و همان‌گونه‌که دربارة نمونة TM12 دیده می‌شود، مقدار LOI بسیار بالاست و نشان می‌دهد این نمونه بسیار دچار هوازدگی و دگرسانی شده است (جدول 3). در زیر میکروسکوپ نیز در این نمونه مقدارهای بالایی از کانی‌های رسی، کلسیت و هیدرواکسیدهای آهن دیده شدند.

برای نامگذاری این سنگ‌ها، نمودار SiO2 در برابر Na2O+K2O به‌کار برده شد. در این نمودار، نمونه‌ها در محدودة داسیت تا تراکی‌داسیت جای می‌گیرند (شکل 12-A). همچنین، در نمودار Nb/Y در برابر Zr/TiO2 نمونه‌ها در محدودة تراکی‌آندزیت جای گرفته‌اند (شکل 12-B).

برای تعیین سری ماگمایی نمونه‌ها نیز نمودارهای Ta/Yb در برابر Ce/Tb و Ta/Yb در برابر Th/Yb به‌کار برده شدند. همان‌گونه‌‌که دیده می‌شود، نمونه‌ها در محدودة سری شوشونیتی جای گرفته‌اند (شکل‌های12- C و 12-D).

سنگ‌های ‌بررسی‌شده می‌توانند بر پایة ضریب آلکالینیته (AI=Al-(K+Na); molecular)، ضریب آهک-آلکالی تعدیل‌شده (MALI= Na2O+K2O-CaO)، ضریب اشباع از سیلیس-فلدسپاتویید (FSSI= Q-[Lc+2(Ne+Kp)]/100; normative)، ضریب آهن (Fe*= FeO+0.9Fe2O3/(FeO+0.9Fe2O3+MgO)) و ضریب اشباع از آلومینیم (ASI= Al/(Ca-1.67P+Na+K); molecular) رده‌بندی شوند (Frost and Frost, 2008).

شکل 12. نامگذاری و تعیین سری ماگمایی سنگ‌های بررسی‌شده به‌همراه مقایسة آنها با سنگ‌های آتشفشانی سهند در: A) نمودار رده‌بندی SiO2 در برابر Na2O+K2O (Le Bas et al., 1986B) نمودار رده‌بندی بر پایة عنصرهای کم‌تحرک (Winchester and Floyd, 1977C) نمودار تعیین سری ماگمایی بر پایة Ta/Yb در برابر Ce/Tb؛ D) نمودار تعیین سری ماگمایی بر پایة Ta/Yb در برابر Th/Yb (Muller and Groves, 1997).

Figure 12. Classification and discrimination of the magmatic series of the studied rocks along with their comparison with Sahand volcanic rocks in: A) SiO2 versus Na2O+K2O classification diagram (Le Bas et al., 1986); B) Immobile elements classification diagram (Winchester and Floyd, 1977); C) Magmatic series in Ta/Yb versus Ce/Tb discrimination diagram; D) Magmatic series based on Ta/Yb versus Th/Yb discrimination diagram (Muller and Groves, 1997).

ضریب آهن برای تعیین غنی‌شدگی آهن سنگ‌های فلدسپاری در هنگام جدایش بلورین مذاب به‌کار برده می‌شود. سنگ‌های منیزیم‌دار در شرایط اکسیدان تحول می‌یابند و غنی‌شدگی اندکی از آهن را نشان می‌دهند؛ اما سنگ‌های آهن‌دار در فوگاسیتة اکسیژن کم دچار جدایش بلورین گسترده‌ای می‌شوند. بیشتر سنگ‌های آلکالن دنیا ویژگی آهن‌دار نشان می‌دهند. با وجود این، برخی ماگماهای حاشیه‌های همگرا (بخش‌های با کشش‌های محلی درون محیط زمین‌ساختی در کل همگرا) از منیزیم‌دار تا آهن‌دار در تغییر هستند (Frost and Frost, 2008). گمان می‌رود سنگ‌های منیزیم‌دار، شامل سنگ‌های بررسی‌شده، از طریق جدابش بلورین مذاب در شرایط اکسیدان و کمابیش آب‌دار پدید آمده باشند (شکل 13-A) (Frost and Frost, 2008). بر پایة ضریب MALI، نمونه‌ها در محدودة آلکالی-کلسیک هستند؛ اما بیشتر نمونه‌های مورد مقایسه (سنگ‌های آتشفشانی سهند) از روند کالک-آلکالی پیروی می‌کنند (شکل 13-B). از لحاظ ضریب اشباع از آلومینیم، بیشتر نمونه‌ها ویژگی متاآلومینوس نشان می‌دهند (شکل 13-C). همان‌گونه‌که در شکل 13-D دیده‌می‌شود همة نمونه‌ها به‌همراه بیشتر نمونه‌های مورد مقایسه، ضریب آلکالینیته و ضریب اشباع از سیلیس-فلدسپاتویید مثبت (FSSI>0) دارند.

در نمودار عنکبوتی عنصرهای ناسازگار نمونه‌ها به ترکیب گوشتة اولیه (Wood et al., 1979) بهنجار شده‌اند. این نمودار نشان می‌دهد نمونه‌ها غنی‌شدگی مشخصی از عنصرهای با شعاع یونی بزرگ (LILE) (مانند: روبیدیم، باریم، توریم، اورانیم، پتاسیم) و عنصرهای خاکی کمیاب سبک نسبت به عنصرهای با میدان پایداری قوی (HFSE) (مانند: تانتالیم، نیوبیم، تیتانیم، زیرکونیم، هافنیم، ایتریم) و عنصرهای خاکی کمیاب سنگین دارند (شکل 14-A). نمونه‌ها مقدار استرانسیم (400 تا 540 پی پی ام) و باریم (930 تا 1130 پی پی ام) بالایی دارند. همچنین، نمونه‌های بررسی‌شده از عنصرهایی مانند Ta، Nb، Ti ناهنجاری منفی نشان می‌دهند (شکل 5). ناهنجاری منفی تانتالیم و نیوبیم همانند سنگ‌های وابسته به پهنه‌های فرورانش ‌حاشیة فعال قاره‌ای است؛ در جایی‌که خاستگاه گوشته‌ای مذاب با مواد فرورونده، از عنصرهای لیتوفیل با شعاع یونی بزرگ دگرنهاد غنی‌ شده است (Pearce, 1982).

الگوی عنصرهای خاکی کمیاب برای نمونه‌های منطقه که به ترکیب کندریت (Boynton, 1984) بهنجار شده‌ است نشان می‌دهد تغییرات این عنصرها در نمونه‌ها ‌به‌صورت موازی است (شکل 14-B). موازی‌بودن الگوی عنصرهای ‌خاکی کمیاب در سنگ‌ها خاستگاه یکسان و سازوکار تبلور بخشی پیدایش آنها را نشان می‌دهد (Wilson, 1989). در این نمودار عنصرهای ‌خاکی کمیاب نمونه‌ها به ترکیب کندریت پیشنهادیِ بوینتون (Boynton, 1984) غنی‌شدگی نشان می‌دهند. در این میان، عنصرهای ‌خاکی کمیاب سبک (LREE) غنی‌شدگی بیشتری نسبت به عنصرهای ‌خاکی کمیاب سنگین (HREE) نشان می‌دهند (شکل 14-B). در نمودار عنصرهای خاکی کمیاب، نمونه‌ها الگو‌های REE پرشیب ((La/Yb)N: 10/19 تا 05/23) و آنومالی مشخص در Eu نشان می‌دهند (شکل 4) که پیامد جدا‌شدن پلاژیوکلاز از مذاب هنگام جدایش بلوری است. نسبت Eu/Eu* به‌دست‌آمدة این سنگ‌ها از 63/0 تا 72/0 در نوسان است.

بحث

باور بر این است که به‌دنبال برخورد ورقه ایران و عربستان در زمان ائوسن در شمال‌باختری ایران (Ghasemi and Talbot, 2006; Hafkenschied et al., 2006; Horton et al., 2008; Omrani et al., 2008). پاره‌شدن و سقوط ورقة اقیانوسی[2]، بالاآمدن مواد سست‌کره‌ای، آشفتگی دمایی و افزایش شیب دمایی در پی آن، به ذوب پوستة زیرین ستبرشده انجامیده است. ذوب پوستة زیرین مافیک گارنت‌دار و بی پلاژیوکلاز (شبیه اکلوژیت)، توانسته است سنگ‌هایی با ویژگی‌های آداکیتی و نیز سنگ‌هایی اسیدی و سرشار از باریم و استرانسیم را پدید آورد (Lechmann et al., 2018).

شکل 13. نمودارهای رده‌بندی بر پایة ضریب آهن، FSSI، AI، MALI و ASI؛ A) نمودار سیلیس در برابر شاخص آهن؛ B) نمودار ضریب آهک-آلکالی در برابر سیلیس که در آن روند‌های سری‌های آلکالیک، آلکالی-کلسیک، کلسیک-آلکالی و کلسیک نشان داده شده است؛ C) نمودار SiO2 در برابر ASI برای تعیین میزان اشباع‌شدگی از آلومینا؛ D) نمودار FSSI در برابر AI ‌برای تعیین اشباع‌شدگی از سیلیس نمونه‌ها (Frost and Frost, 2008) (نماد نمونه‌ها همانند شکل 12 است).

Figure 13. The classification diagrams of Fe-index, FSSI, AI, MALI and ASI; A) Silica versus Fe-index diagram; B) alkali-lime index versus silica diagram, in which the trends of the alkaline, alkali-calcic, calcic-alkali and calcic series are shown; C) SiO2 versus ASI diagram for determination of alumina saturation; D) FSSI versus AI for determination of silica saturation (Frost and Frost, 2008) (The symbols are the same as in Figure 12).

شکل14. ترکیب سنگ‌های ‌بررسی‌شده در: A) نمودار عنکبوتی بهنجارشده به ترکیب گوشتة اولیه (Wood et al., 1979B) نمودار REE بهنجارشده به ترکیب کندریت (Boynton, 1984).

Figure 14. Composition of studied rocks in: A) The primitive mantle-normalized spider diagram (normalized values from Wood et al. (1979)); B) REE chondrite-normalized diagram (normalized values from Boynton (1984)).

همچنین، ذوب گوشته سنگ‌کره‌ای دگرنهادشده در پی فرایندهای فرورانش، نیز می‌تواند به پیدایش ماگمای مافیک شوشونیتی بیانجامد و تحول آن ماگمای حد واسط و اسیدی را پدید آورده باشد (Ghalamghash et al., 2019). ویژگی زمین‌شیمیایی نمونه‌ها به‌همراه سنگ‌های مشابه از مناطق مختلف آتشفشان سهند نشان‌دهندة پیدایش ماگما‌های آداکیتی و غیرآداکیتی (شامل کالک‌آلکالن عادی و شوشونیتی) است. همان‌گونه‌که نمودارهای Y در برابر Sr/Y و نیز YbN در برابر (La/Yb)N (شکل‌های 15-A و 15-B) نشان می‌دهند، طیف‌های ترکیبی گوناگونی از سنگ‌ها را می‌توان در این منطقه شناسایی کرد که می‌تواند با رخداد جدایش بلورین گسترده ماگما‌های تولیدشده و یا خاستگاه‌های مختلف آنها در ارتباط باشد. با وجود این، نمونه‌ها در محدودة غیرآداکیتی جای گرفته‌اند و ازاین‌رو، چه‌بسا پیدایش آنها با سنگ‌های آداکیتی گوناگون در مجموعة آتشفشانی سهند متفاوت باشد.

شکل 15. ترکیب نمونه‌های داسیتی بررسی‌شده در نمودار‌های متمایزکنندة آداکیت‌ها و سنگ‌های معمولی کمان ماگمایی. A) نمودار Y در برابر Sr/Y (Defant and Drummond, 1990B) نمودار YbN در برابر (La/Yb)N (Martin, 1986).

Figure 15. Composition of the studied dacitic samples in adakites and typical arc rocks discrimination diagrams. A) Y versus Sr/Y plot (Defant and Drummond, 1990; B) YbN versus (La/Yb)N plot (Martin, 1986).

دربارة پیدایش سنگ‌های با گرایش آداکیت در مجموعة سهند، باور بر این است که ذوب‌بخشی پوستة زیرین (Lechmann et al., 2018; Fedele et al., 2023) یا ذوب صفحة فرورونده، پیدایش چنین مذاب‌هایی را به‌دنبال داشته است (Chaharlang et al., 2016; Ghalamghash and Chaharlang, 2014)؛ اما دربارة ماگما‌های غیرآداکیتی باور بر خاستگاه‌گرفتن آنها از گوشتة سنگ‌کره‌ای دگرنهادشده است (Ghalamghash et al., 2019; Fedele et al., 2023). ضخیم‌شدگی پوسته‌ای، پارگی سنگ‌کرة اقیانوسی فرورونده، عقب‌نشینی سنگ‌کرة اقیانوسی فرورونده و ورقة ورقه‌شدن و آشفتگی دمایی از سازوکار‌های تولید مذاب با ترکیب آداکیتی و غیرآداکیتی در این ناحیه دانسته شده‌اند و از نقش فرورانش نئوتتیس در این ارتباط یاد شده است. ازاین‌رو، سنگ‌های بررسی‌شده با ماگماتیسم مرتبط با نئوتتیس ارتباط تنگاتنگی دارند؛ همان‌گونه‌که ترکیب شیمیایی این سنگ‌ها در نمودار‌های متمایز‌کننده محیط‌های زمین‌ساختی (شکل 16) گویای گرایش آنها به محیط‌های زمین‌ساختی کمان قاره‌ای است. با وجود این، همان‌گونه‌که پیشتر گفته شد، شواهد همزمانی پیدایش این سنگ‌ها با فرورانش فعال (در وضعیت پیش از برخورد) وجود ندارد. ازاین‌رو، پیدایش این سنگ‌ها در یک پهنة کمانس پسابرخوردی دور از انتظار نخواهد بود. در نمودارسه‌تایی TiO2/100-La-10Hf که ‌برای تمایز کمان‌های قاره‌ای و پسابرخوردی (CAP+PAP) از کمان‌های اقیانوسی (IOP+LOP) به‌کار برده می‌شود، همة نمونه‌ها در محدوده کمان‌های قاره‌ای و پسابرخوردی جای می‌گیرند (شکل 16-A). در نمودار سه‌تایی 3Zr-50Nb-Ce/P2O5 که ‌برای تمایز کمان پسابرخوردی (PAP) از کمان قاره‌ای (CAP) به‌کار برده می‌شود، نمونه‌ها در محدودة کمان قاره‌ای جای گرفته‌اند (شکل 16-B). در نمودار Zr در برابر (Nb/Zr)n، نمونه‌ها در محدودة سنگ‌های پسابرخوردی جای گرفته‌اند (شکل 16-C). جانمایی نمونه‌ها در نمودار La/Yb در برابر Th/Nb که ‌برای تمایز کمان‌های مختلف از همدیگر به‌کار برده می‌شود، نشان می‌دهد نمونه‌ها به کمان‌های آلکالن گرایش دارند (شکل 16-D).

سنگ‌های ‌بررسی‌شده با توجه به ترکیب شیمیایی کانی‌ها (مانند: سدیک‌بودن پلاژیوکلاز‌ها و داشتن کانی‌های فرومنیزین آبدار (بر پایة بیوتیت) و با در نظرگرفتن ترکیب سنگ کل با سیلیس بیشتر از 60 درصدوزنی ‌به‌عنوان سنگ‌های جدایش‌یافته و تحول‌یافته اسیدی به‌شمار می‌روند. توف‌های داسیتی بررسی‌شده متاآلومینوس، کالک‌آلکالن پتاسیم بالا تا شوشونیتی هستند و الگوی عنصرهای کمیاب آنها، ویژگی‌های شاخص ماگماتیسم نوع کمانی را نشان می‌دهد. نسبت Nb/Ta در نمونه‌ها برابر با 7/14 تا 8/15 است که از مقدار پیش‌بینی‌شده برای پوستة قاره‌ای بیشتر است (Taylor and McLennan, 1985)؛ اما همانند سنگ‌های آتشفشانی کمانی است (Stolz et al., 1996). ویژگی‌های یادشده به‌همراه آنومالی منفی Nb، Ta و Ti و نسبت‌های بالای Ba/La (17 تا 6/20) و Ba/Zr (8/5 تا 7/6) و Ba/Nb >30 (Gill, 1981) همگی شباهت آنها به ماگما‌های مرتبط با فرورانش را نشان می‌دهند. با وجود این، شواهد زمین‌شناسی در منطقه مبنی بر وجود فرورانش فعال در هنگام پیدایش سنگ‌های بررسی‌شده وجود ندارد. ازاین‌رو، گمان می‌رود این ویژگی‌های ترکیبی در ارتباط با خاستگاه باشد تا نشان‌دهندة خاستگاه زمین‌ساختی. فراوانی بالای Sr (540-430 پی پی ام) همراه با آنومالی منفی Eu (73/0 – 62/0 Eu/Eu*=) و تهی‌شدگی نسبی HREE و Y نشان‌دهندة تفالة ذوب با ترکیب گارنت، آمفیبول و در مقدارهای کمتر یا بی پلاژیوکلاز بوده است. چنین شرایطی می‌تواند در پی تبدیل پوستة اقیانوسی فرورونده به گارنت آمفیبولیت یا اکلوژیت و ذوب‌بخشی آنها فراهم شود. البته ذوب‌بخشی اخیر به‌همراه آزاد‌شدن سیال‌ها از ورقة اقیانوسی به نوبه خود نیز می‌توانند گوشتة سنگ‌کره‌ای را دگرنهاد کنند تا مذاب‌هایی با ویژگی‌های زمین‌شیمیایی یادشده پدید آید. گفتنی است سنگ‌های بررسی‌شده با نسبت‌های بالای La/Nb (1/2 تا 4/2) شناخته می‌شوند که اساساً به نواحی با خاستگاه گوشتة سنگ‌کره‌ای زیرقاره‌ای نسبت داده می‌شود (Hawkesworth et al., 1995; Rogers et al., 1995).

غنی‌شدگی سنگ‌های بررسی‌شده از برخی عنصرها مانند Ba، Sr و Rb نیازمند تبلور گسترده، آلایش پوسته‌ای و یا ذوب‌بخشی بسیار کم است. ترکیب ایزوتوپی Sr (704792/0 تا 708505/0) و Nd (512388/0 تا 512802/0) سنگ‌های با سنگ‌شناسی و سن مشابه (Lechmann et al., 2018) در منطقة ‌بررسی‌شده نشان می‌دهد تنها آلایش پوسته‌ای نمی‌تواند این تغییرات گسترده را به‌دنبال داشته باشد. افزون‌بر این، نسبت‌های بالای LILE/HFSE و LREE/HREE می‌تواند پیامد نرخ اندک ذوب‌بخشی دانسته شود.

شکل 16. نمودار‌های متمایز کنندة محیط‌های زمین‌ساختی. A) نمودار سه‌تایی TiO2/100-La-10Hf (Muller and Groves, 1997) ‌برای تمایز کمان‌های قاره‌ای و پسابرخوردی (CAP+PAP) از کمان اقیانوسی (IOP+LOP)؛ B) نمودار سه‌تایی 3Zr-50Nb-Ce/P2O5 (Muller and Groves, 1997) برای تمایز کمان پسابرخوردی (PAP) از کمان قاره‌ای (CAP)؛ C) نمودار Zr در برابر (Nb/Zr)n (Thiéblemont and Tegyey, 1994) که در آن N نماد بهنجارسازی به ترکیب NMORB (Sun and McDonough, 1989) است؛ D) نمودار La/Yb در برابر Th/Nb (Hollocher et al., 2012) برای شناخت کمان‌های گوناگون از یکدیگر.

Figure 16. Tectonic setting discrimination diagrams. A) TiO2/100-La-10Hf ternary diagram (Muller and Groves, 1997) for discrimination the continental and post-collisional arcs (CAP+PAP) from oceanic arcs (IOP+LOP); B) 3Zr-50Nb-Ce/P2O5 ternary diagram (Muller and Groves, 1997) for discrimination between post-collisional arc (PAP) and continental arc (CAP); C) plot of Zr versus (Nb/Zr)n (Thiéblemont and Tegyey, 1994), where n indicates values normalized to NMORB (Sun and McDonough, 1989); D) La/Yb versus Th/Nb diagram (Hollocher et al., 2012) for discrimination various arc magmas.

با توجه به اینکه آتشفشان سهند روی پی‌سنگ متشکل از نهشته‌های رسوبی پالئوزوییک-مزوزوییک، سنگ‌های آتشفشانی ائوسن، مارن‌های میوسن زیرین (سازند قم) و مارن‌ها و لایه‌های سرخ‌رنگ ماسه‌سنگی و ژیپس (سازند قرمز ‌بالایی) به سن میوسن میانی تا پسین جای گرفته است (Abbassi et al., 2021) و از سوی دیگر، کهن‌ترین فعالیت آتشفشان سهند به سن میوسن پسین دانسته می‌شود (سهند قدیم در رده‌بندی قلمقاش و همکاران (Ghalamghash et al., 2019) با سنی نزدیک به 8 میلیون سال پیش و توالی آذرآواری ضخیم در دامنة باختری آتشفشان با نام سازند قرمزی‌گل در رده‌بندی معین وزیری و امین سبحانی (Moine Vaziri and Amine Sobhani, 1977) با سن 9-12 میلیون سال پیش). ازاین‌رو، با در نظرگرفتن موقعیت چینه‌ای لایة توفی یادشده و شباهت ترکیب شیمیایی آن به فوران‌های غیرآداکیتی سهند، گمان می‌رود این لایة توفی پیامد نخستین فعالیت‌های انفجاری سهند در میوسن پایانی، همزمان با پیدایش رسوبات قرمز بالایی (میوسن پسین) باشد. با توجه به رخداد سنگ‌های بررسی‌شده در نزدیکیِ گسل تبریز و با در نظر گرفتن اینکه بیشترین فعالیت این گسل در میوسن پایانی روی داده است (Reichenbacher et al., 2011)، نقش کلیدی آن در ارتباط با فوارن‌های آتشفشانی سهند دور از انتظار نیست.

برداشت

رخداد لایة توفی با ترکیب داسیتی درون سازند قرمز ‌بالایی (میوسن پسین) در نزدیکیِ گسل تبریز در نوع خود جالب توجه است. این ماگماتیسم می‌تواند یکی از نخستین فوران‌های آتشفشانی سهند در میوسن پسین دانسته شود. از دیدگاه سنگ‌نگاری، این سنگ‌ها بافت پورفیری با فنوکریست‌های پلاژیوکلاز، بیوتیت و کوارتز در زمینة شیشه‌ای (ویتروفیری) دارند. ترکیب اسیدی (وجود کوارتز و پلاژیوکلاز‌های سدیک)، تغییرات گستردة فوگاسیتة اکسیژن ماگما با توجه به ترکیب بیوتیت‌ها، دمای کمابیش کم با توجه به دماسنجی پلاژیوکلاز‌ها و داشتن کانی‌های آبدار نشان‌دهندة فوران‌های انفجاری ماگمای سازندة سنگ‌های یادشده است. داده‌های زمین‌شیمیایی، گویای ترکیب کالک‌آلکالن با پتاسیم بالا و شوشونیتی ماگماست و با توجه به متاآلومینوس‌بودن از نظر میزان آلومینیم و شیب و ناهنجاری‌های الگوی چندعنصری، عنصرهای کمیاب و خاکی کمیاب گمان نمی‌رود این سنگ‌ها تنها خاستگاه پوستة قار‌ه‌ای داشته باشند. مقایسة ترکیب شیمیایی سنگ‌های بررسی‌شده با سنگ‌های مجموعه آتشفشانی سهند نشان‌دهندة تشابه آنها به سنگ‌های غیرآداکیتی این مجموعه است و ‌ویژگی‌هایی مانند آنومالی منفی Ta، Nb و Ti، داشتن نسبت‌های بالای Ba/La، Ba/Zr، Ba/Nb و La/Nb و داشتن فراوانی بالای Sr و Ba می‌تواند به خاستگاه گوشتة سنگ‌کره‌ای آنها ارتباط داشته باشد که در یک محیط زمین‌ساختی مرتبط با فرورانش یا ‌حاشیة فعال قاره‌ای پدید آمده‌اند. با توجه به بازة زمانی فعالیت آتشفشان سهند که از میوسن پسین تا پلیستوسن پایانی دانسته شده است و در پی این فعالیت آتشفشانی، حجم بزرگی از مواد آذرآواری در جنوب گسل تبریز پدید آمده است، گمان می‌رود لایة توفی بررسی‌شده، با توجه به موقعیت چینه‌ای و مکانی آن، از فراورده‌های آتشفشانی و انفجاری میوسن پسین سهند باشد که همزمان با رسوبات قرمز بالایی (میوسن پسین) که در بخش شمالی گسل تبریز در حال پیدایش بوده است، نهشته شده است.

سپاس‌گزاری

نگارندگان از معاونت محترم پژوهشی و فنآوری دانشگاه صنعتی سهند تبریز و نیز از سردبیر وداوران گرامی مجلة پترولوژی برای همکاری صمیمانه‌شان و کمک به بهبود سطح علمی مقاله سپاس‌گزاری می‌کنند.

 

[1] Upper Red Formation

[2] slab break-off

Abbassi, N., Khoshyar, M., Lucas, S.G., and Esmaeili, F. (2021) Extensive vertebrate track site from the Upper Red Formation (middle-late Miocene), west Zanjan, northwestern Iran. Fossil Record, 24, 101-116. https://doi.org/10.5194/fr-24-101-2021
Abdel-Rahman, A.M. (1994) Nature of biotites from alkaline, calc-alkaline, and peraluminous magmas. Journal of Petrology, 35, 525–541. https://doi.org/10.1093/petrology/35.2.525
Alavi, M. (1994) Tectonics of the Zagros orogenic belt of Iran: new data and interpretations. Tectonophysics, 229, 211–238. https://doi.org/10.1016/0040-1951(94)90030-2
Alavi, M., and Shahrabi, M. (1975) Geological Quadrangle map of Maragheh 1/100,000.
Allen, M., Jackson, J.A., and Walker, R. (2004) Late Cenozoic reorganization of the Arabia-Eurasia collision and the comparison of short-term and long-term deformation rates. Tectonics, 23, 1-16. https://doi.org/10.1029/2003TC001530
Asadian, O. (1993) Geological Quadrangle map of Tabriz 1/100,000.
Ashrafi, N., Jahangiri, A., Ameri, A., Hasebe, N., and Eby, G.N. (2009) Biotite mineral chemistry of the Bozqush and Kaleybar intrusions, NW Iran. Iranian Journal of Crystallography and Mineralogy, 17, 381–394 (in Persian). http://ijcm.ir/files/site1/user_files_c470dc/niyaiifar-A-10-134-242-d5d9dce.pdf
Axen, G.J., Lam, P.S., Grove, M., Stockli, D.F., and Hassanzadeh, J. (2001) Exhumation of the west-central Alborz Mountains, Iran, Caspian subsidence, and collision-related tectonics. Geology, 29, 559–562. https://doi.org/10.1130/0091-7613(2001)029%3C0559:EOTWCA%3E2.0.CO;2
Behrouzi, A., Amini Fazl, A., and Amini Azar, B. (1997) Geological Quadrangle map of Bostan Abad 1/100,000.
Berberian, M., and King, G.C.P. (1981) Towards a paleogeography and tectonic evolution of Iran. Canadian Journal of Earth Sciences, 18, 210–265. http://doi.org/10.1139/e81-019
Boynton, W.V. (1984) Cosmochemistry of the earth elements: meteorite studies. In: Henderson, R. (Ed.), Rare Earth Element Geochemistry: Developments in Geochemistry, p. 89-92. Elsevier, Amsterdam. https://doi.org/10.1016/b978-0-444-42148-7.50008-3
Chaharlang, R., Ghalamghash, J., Hashem Emami, M., and Omrani, J. (2016) Source of Kandovan Volcanic Rocks (NW of Iran): Using New Data of Geology, Petrography and Geochemistry. Advanced Applied Geology, 5, 34-43 (in Persian). https://aag.scu.ac.ir/article_11897_876b2a0c7a374363de44959e76bcd216.pdf
Deer, W.A., Howie, R.A., and Zussman, J. (1992) An introduction to the rock forming minerals. Second Edition, 696 p. Longman, London.
Defant, M.J., and Drummond, M.S. (1990) Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature, 347, 662–665. https://doi.org/10.1038/347662a0
Fedele, L., Ghazi, J.M., Agostini, S., Ronca, S., Innocenzi, F., and Lustrino, M. (2023) Concurrent adakitic and non-adakitic Late Miocene-Quaternary magmatism at the Sahand volcano, Urumieh-Dokhtar Magmatic Arc (NW Iran), Lithos, 458-458, 107344. https://doi.org/10.1016/j.lithos.2023.107344
Ferguson, A.K., and Cundari, A. (1975) Petrological aspects and evolution of the leucite-bearing lavas from Bufumbira, southwest Uganda, Contribution to Mineralogy and Petrology. 50, 25-46. https://doi.org/10.1007/BF00385219
Fisher, W.B. (1968) Physical geography. In: Fisher, W.B. (Ed.), The Cambridge History of Iran: The Land of Iran, 1. Cambridge University Press, Cambridge, pp. 3–110. https://doi.org/10.1017/CHOL9780521069359
Foster, M.D. (1960) Interpretation of the composition of trioctahedral micas. United States Geological Survey Professional Paper, 354-B, 11–49. https://doi.org/10.3133/pp354B
Fuhrman, M.L., and Lindsley, D.H. (1988) Ternary-feldspar modeling and thermometry. Am. Mineral. 73, 201-215.
Frost, B.R., and Frost, C.D. (2008) A Geochemical Classification for Feldspathic Igneous Rocks. Journal of Petrology, 49, 1955-1969. https://doi.org/10.1093/petrology/egn054
Gansser, A. (1955) New aspects of the geology in central Iran. 4th World Petroleum Congress. Rome, Italy.
Ghadirzadeh, A. (2002) Geological Quadrangle map of Azar Shahr 1/100,000.
Ghalamghash, J., and Chaharlang, R. (2014) Petrogenesis and tectonic setting of Sahand subvolcanic and volcanic domes. Petrology, 18,123-139 (in Persian). https://ijp.ui.ac.ir/article_16168.html?lang=fa
Ghalamghash, J., Schmitt, A.K., and Chaharlang, R. (2019) Age and compositional evolution of Sahand volcano in the context of post-collisional magmatism in northwestern Iran: Evidence for time-transgressive magmatism away from the collisional suture. Lithos, 344-345, 265-279. https://doi.org/10.1016/j.lithos.2019.06.031
Ghasemi, A., and Talbot, C.J. (2006) A new tectonic scenario for the Sanandaj–Sirjan Zone (Iran). Journal of Asian Earth Sciences, 26, 683–693. https://doi.org/10.1016/j.jseaes.2005.01.003
Ghauori, N. (2002) Study of stratigraphy, petrography and geochemistry of ignimbrite sequences of Sahand. M.Sc. thesis, Tarbiat Moalem University, Tehran, Iran (in Persian)
Gill, J.B. (1981) Orogenic Andesites and Plate Tectonics, 385 pp. Springer, New York. https://doi.org/10.1007/978-3-642-68012-0
Hafkenscheid, E., Wortel, M., and Spakman, W. (2006) Subduction history of the Tethyan region derived seismic tomography and tectonic reconstruction. Journal of Geophysical Research Atmospheres, 111, B08401. https://doi:10.1029/2005JB003791
Hawkesworth, C., Turner, S., Gallagher, K., Hunter, A., Gallagher, K., Hunter, A., Bradshaw, T., and Rogers, N. (1995) Calcalkaline magmatism, lithospheric thinning and extension in the Basin and Range. J. Geophys. Res. 100, 10271–10286. https://doi.org/10.1029/94JB02508
Henry, D.J., and Guidotti, C.V. (2002) Titanium in biotite from metapelitic rocks: Temperature effects, crystal-chemical controls, and petrologic applications. American Mineralogist, 87, 375-382. https://doi.org/10.2138/am-2002-0401
Henry, D.J., Guidotti, C.V., and Thomson, J.A. (2005) The Ti-saturation surface for low to medium pressure metapelitic biotites: Implications for geothermometry and Ti-substitution mechanisms. American Mineralogist, 90, 316–328. https://doi.org/10.2138/am.2005.1498
Hollocher, K., Robinson, P., Walsh, E., and Roberts, D. (2012) Geochemistry of amphibolite-facies volcanics and gabbros of the Storen Nappe in extensions west and southwest of Trondheim, western gneiss region, Norway: A key to correlations and paleotectonic settings. American Journal of Science 312, 357–416. https://doi.org/10.2475/04.2012.01
Horton, B.K., Hassanzadeh, J., Stockli, D.F., Axen, G.J., Gillis, R.J., Guest, B., Amini, A., Fakhari, M.D., Zamanzadeh, S.M., and Grove, M. (2008) Detrital zircon provenance of Neoproterozoic to Cenozoic deposits in Iran: implications for chronostratigraphy and collisional tectonics. Tectonophysics, 451, 97–122. https://doi.org/10.1016/j.tecto.2007.11.063
Hyndman, D.W. (1985) Petrology of Igneous and Metamorphic Rocks, Second Edition, 720 p. McGraw-Hill.
Jahangiri, A. (2007) Post-collisional Miocene adakitic volcanism in NW Iran: Geochemical and geodynamic implication. Journal of Asian Earth Sciences, 30, 433-447. https://doi.org/10.1016/j.jseaes.2006.11.008
Jahangiri, A., and Ashrafi, N. (2006) Transition from phreatomagmatic to strombolian eruptions in the Sahand volcano; constrain from pyroclastic studies. 25th Geosciences Congress, Geological Survey of Iran, Tehran, Iran (in Persian).
Janousek, V., Farrow, C.M., and Erban, V. (2006) Interpretation of whole-rock geochemical data in igneous geochemistry: introducing Geochemical Data Toolkit (GCDkit). Journal of Petrology 47, 1255-1259. http://dx.doi.org/10.1093/petrology/egl013
Keskin, M., Pearce J.A., and Mitchell, J.G. (1998) Volcano-stratigraphy and geochemistry of collision-erlated volcanism on the Eezurum-Kars Plateau, northeastern Turkey. Journal of Volcanology and Geothermal Research, 85, 355-404. https://doi.org/10.1016/S0377-0273(98)00063-8
Keskin, M. (2003) Magma generation by slab steepening and breakoff beneath a subduction accretion complex: an alternative model for collision-related volcanism in eastern Anatolia, Turkey. Geophysical Research Letters 30, 1-4. https://doi.org/10.1029/2003GL018019
Khodabandeh, A.A., and Amini Fazl, A. (1995) Geological Quadrangle map of Osku 1/100,000.
Le Bas, M.J., Le Maitre, R.W., Streckeisen, A., and Zanetti, B. (1986). A chemical classification of volcanic rock based on the total alkali–silica diagram. Journal of Petrology, 27, 745–750. https://doi.org/10.1093/petrology/27.3.745
Lechmann, A., Burg, J.P., Ulmer, P., Guillong, M., and Faridi, M. (2018) Metasomatized mantle as the source of Mid-Miocene-Quaternary volcanism in NW-Iranian Azerbaijan: Geochronological and geochemical evidence. Lithos, 304-307, 311-328. https://doi.org/10.1016/j.lithos.2018.01.030
Liégeois, J.P., Navez, J., Hertogen, J., and Black, R. (1998) Contrasting origin of post-collisional high-K calc-alkaline and shoshonitic versus alkaline and peralkaline granitoids. The use of sliding normalization. Lithos, 45, 1-28. http://doi.org/10.1016/S0024-4937(98)00023-1
Martin, H. (1986) Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone magmas. Geology, 14, 753–756. https://doi.org/10.1130/0091-7613(1986)14%3C753:EOSAGG%3E2.0.CO;2
Mohajjel, M., Fergusson, C.L., and Sahandi, M.R. (2003) Cretaceous–Tertiary convergence and continental collision, Sanandaj–Sirjan Zone, western Iran. Journal of Asian Earth Sciences, 4, 397–412. https://doi.org/10.1016/S1367-9120(02)00035-4
Moine Vaziri, H. (1999) An introduction to the magmatism in Iran. 440 p. Kharazmi University Publication, Tehran (in Persian).
Moine Vaziri, H., and Amine Sobhani, E. (1977) Carte Geologique du Sahand 1:100,000 (Ecole normale superieure de Tehran).
Morrison, G.W. (1980) Characteristics and tectonic setting of the shoshonite rock association, Lithos, 13, 97- 108. https://doi.org/10.1016/0024-4937(80)90067-5
Muller, D., and Groves, D.I. (1997) Potassic igneous rocks and associated gold-copper mineralization, 238 p. Springer. https://link.springer.com/book/10.1007/978-3-319-23051-1
Nekvasil, H. (1992) Ternary feldspar crystallization in high-temperature felsic magmas. American Mineralogist, 77, 592-604.
Nachit, H., Razafimahefa, N., Stussi, J.M., and Carron, J.P. (1985) Composition chimique des biotites et typologie magmatique des granitoides. Comptes Rendus Hebdomadaires de l’ Académie des Sciences, 301, 813–818.
Nachit, H., Ibhi, A., Abia, E.H., and Ohoud, M.B. (2005) Discrimination between primary magmatic biotites, reequilibrated biotites and neoformed biotites. Geomaterials (Mineralogy), Copmtes Rendus, Geoscience, 337, 1415–1420. http://dx.doi.org/10.1016/j.crte.2005.09.002
Omrani, J., Agard, P., Whitechurch, H., Benoit, M., Prouteau, G., and Jolivet, L. (2008) Arc magmatism and subduction history beneath the Zagros Mountains, Iran: a new report of adakites and geodynamic consequences. Lithos, 106, 380–398. http://dx.doi.org/10.1016/j.lithos.2008.09.008
Pearce, J.A. (1982) Trace element characteristics of lavas from destructive plate boundary, In: Thorpe, R.S. (Ed.), Andesites, p. 525-548. Wiley, New York. https://orca.cardiff.ac.uk/id/eprint/8625/
Pirmohammadi, F., Ameri, A., Jahangiri, A., Mojtahedi, M., Chen, C., and Kesckin, M. (2011) Mineralogy, petrography and petrogenesis of Sahand Volcanic Suite, Northwest. Iranian Journal of Crystallography and Mineralogy, 19, 83–102 (in Persian). http://ijcm.ir/article-1-469-en.html
Reichenbacher, B., Alimohammadian, H., Sabouri, J., Haghfarshi, E., Faridi, M., Abbasi, S., Matzke-Karasz, R., Fellin, M.G., Carnevale, G., Schiller, W., Vasilyan, D., and Scharrer, S. (2011) Late Miocene stratigraphy, palaeoecology and palaeogeography of the Tabriz Basin (NW Iran, Eastern Paratethys). Palaeogeography, Palaeoclimatology, Palaeoecology. 311, 1–18. https://doi.org/10.1016/j. palaeo.2011.07.009
Rezaei, H., and Fonoudi, H. (2009) Geological Quadrangle map of Ghareh Aghaj 1/100,000.
Rieder, M., Cavazzini, G., Yakonov, Y.D., Frank-Kanetskii, V.A., Gottardi, G., Guggenheim, S., Koval, P.V., Müller, G., Neiva, A.M.R., Radoslovich, E.W., Robert, J.L., Sassi, F.P., Takeda, H., Weiss, Z., and Wones, D.R. (1998) Nomenclature of the micas. Canadian Mineralogist 36, 905–912. http://doi.org/10.1346/CCMN.1998.0460513
Rogers, N.W., Hawkesworth, C.J., and Ormerod, D.S. (1995) Late Cenozoic basaltic magmatism in the Western Great Basin California and Nevada. Journal of Geophysical Research, 100, 10287–10301. https://doi.org/10.1029/94JB02738
Shahabpour, J. (2007) Island-arc affinity of the Central Iranian Volcanic Belt, Journal of Asian Earth Sciences 30, 652–665. https://doi.org/10.1016/j.jseaes.2007.02.004
Soltani Sisi, G.A. (2005). Geological Quadrangle map of Ajab Shir 1/100,000.
Stolz, A.J., Jochum, K.P., Spettel, B., and Hofmann, A.W. (1996) Fluid and melt-related enrichment in the subarc mantle: Evidence from Nb/Ta variations in island-arc basalts. Geology, 24, 587–590. https://doi.org/10.1130/0091-7613(1996)024%3C0587:FAMREI%3E2.3.CO;2
Sun, S.S., and McDonough, W.F. (1989) Chemical and isotopic systematic of oceanic basalts: implication for mantle composition and processes. Geol. Soc. Am. Spec. Publ. 42, 313–345. http://dx.doi.org/10.1144/GSL.SP.1989.042.01.19
Thompson, R.N. (1985) Asthenospheric source of Ugandan ultrapotassic magma? Journal of Geology, 93, 603–608. https://doi.org/10.1086/628984
Taylor, S.R., and McLennan, S.M. (1985) The Continental Crust: Its Composition and Evolution, 312 pp. Blackwell, Oxford. https://doi.org/10.1017/S0016756800032167
Thiéblemont, D., and Tegyey, M. (1994) Une discrimination géochimique des Roches différenciées témoin de la diversité d’origine et de situation Tectonique des magmas calco-alcalins. Comptes Rendus Geoscience, 319, 87-94.
Venturelli, G., Thorpe, R.S., Dal Piaz, G.V., Del Moro, A., and Potts, P.J. (1984) Petrogenesis of calc-alkaline, shoshonitic and associated ultrapotassic Oligocene volcanic rocks from the Northwestern Alps, Italy. Contribution to Mineralogy and Petrology, 86, 209–220. https://doi.org/10.1007/BF00373666
Wang, Q., Xu, J.F., Jian, P., Bao, Z.W., Zhao, Z.H., Li, C.F., Xiong, X.L., and Ma, J.L. (2006) Petrogenesis of Adakitic Porphyries in an Extensional Tectonic Setting, Dexing, South China: Implications for the Genesis of Porphyry Copper Mineralization. Journal of Petrology, 47(1), 119–144. https://doi.org/10.1093/petrology/egi070
Wilson, M. (1989) Igneous Petrogenesis; A Global Tectonic Approach, 466 p. Unwin Hyman, London. https://doi.org/10.1007/978-1-4020-6788-4
Winchester, J.A., and Floyd, P.A. (1977) Geochemical discrimination of different magma series and their differentiation products using immobile elements. Chemical Geology, 20, 325–343. http://doi.org/10.1016/0009-2541(77)90057-2
Wones, D.R., and Eugster, H.P. (1965) Stability of biotite: experiment, theory, and application. American Mineralogist, 50, 1228-1272. http://www.minsocam.org/ammin/AM50/AM50_1228.pdf
Wood, I.A., Joron, J.L., Treuil, M., Norry, L., and Tarney, J. (1979) Elemental and Sr isotope variations in basic lavas from Iceland and surrounding ocean floor. Contribution to Mineralogy and Petrology, 70, 319-339. https://doi.org/10.1007/BF00375360
Yavuz, F. (2003) Evaluating micas in petrologic and metallogenic aspect: I–definitions and structure of the computer program MICA+. Computers & Geosciences, 29, 1203–1213. http://dx.doi.org/10.1016/S0098-3004(03)00142-0
Zhou, Z.X. (1986) The origin of intrusive mass in Fengshangdong, Hubei Province. Acta Petrologica Sinica, 2, 59-70.
 
Volume 15, Issue 4 - Serial Number 60
Petrological Journal, 15th Year, No. 60, Winter 2025
December 2025
Pages 69-96
  • Receive Date: 23 May 2024
  • Revise Date: 03 August 2024
  • Accept Date: 10 August 2024