Combination of geochemical- statistical- petrological and crystal size distribution (CSD) as a new confidence in determination of magmatic processes in Kodegan Area (South Khorasan)

Document Type : Original Article

Authors

1 Faculty of Geology, University of Tehran, Tehran, Iran

2 Faculty of Basic Sciences, Tarbiat Modares University, Tehran, Iran

3 Islamic Azad University, Science and Research Branch, Department of Earth Sciences, Tehran, Iran

Abstract

Eocene volcanic and plutonic rocks exposed in the 90 km south of Khusf, south Khorasan Province. The texture of these rocks is porphyritic, glomeroporphyritic, seriate and poikilitic. Plagioclase, Alkali feldspar, quartz, pyroxene and hornblende are abundant mineral in these rocks. Plagioclase crystals show sieve texture and zoning. These rocks are calc-alkaline to tholeiiti. They characterized by enrichment of LREE and LILE, depletion of HREE and HFSE. Negative anomalies of Ti, Nb, P, Zr, and positive anomalies of Pb, K, U reflect emplacement of volcanic rocks in an active continental margin. In order to investigation of magma crystallization, crystal size distribution (CSD) technique also used for 10 andesite and basalt rock samples. Increasing of plagioclase crystal size in a semi-logarithmic plot which demonstrate population density versus size, represents fractional crystallization. Concavity of the CSD diagrams resulting from combination of two populations of crystals, that one of them formed during slow cooling and the other formed during a period of pressure reduction and rapid cooling. Additionally, chemical features such as variation trends of Nb/Y vs. Nb and MgO vs. Fe2O3 and two different trends of volcanic and plutonic samples in Harker and dendrogram diagram, regarded as an indicator of magma mixing processes.

Keywords

Main Subjects


 

یکی از روش‏‌‏‌های کاربردی که اطلاعات مهمی دربارة تاریخچه تبلور و سردشدن سنگ‏‌‏‌های ماگمایی ارائه می‌دهد، بررسی کمّی اندازه کانی‏‌‏‌های سنگ‏‌ساز سنگ‏‌‏‌های آذرین و به‏‌ویژه در سنگ‌های آذرین بیرونی است. این روش به‌نام توزیع اندازة بلور‏‌ها یا CSD (Crystal Size Distribution) خوانده می‏‌شود (Marsh, 1988). روش CSD بیشتر متأثر از هسته‏‌بندی، متوسط سرعت و زمان رشد بلورها و به‌دنبال آن، زمان توقف ماگما در آشیانة ماگمایی است.

ازآنجایی‌که تبلور سنگ‏‌‏‌های آذرین در شرایط پیچیده هسته‏‌بندی و رشد انجام می‏‌شود، سردشدگی سریع در سنگ‏‌‏‌های آتشفشانی یا نیمه‏‌ژرف یا سردشدن آرام در سنگ‏‌‏‌های آذرین درونی تغییرات آشکار و قابل دیدنی در بافت آنها پدید می‏آورد. عامل کلیدی در درک چگونگی پیدایش سنگ‏‌‏‌های آذرین و فرایند‏‌های مؤثر در پیش، همزمان و پس از پیدایش آنها، نیازمند تفسیر ترکیب شیمیایی و بافت آنهاست (Renzulli and Santi, 1997; Jerram and Cheadle, 2000; Kuscu and Floyd, 2001). تجزیة کیفی بافت‏‌‏‌ها پایة فهم ویژگی‏‌‏‌های فیزیکی تبلور ماگما به‌شمار می‌رود (Marsh, 1988; Higgins, 2006a; Jerram and Higgins, 2007).

روش CSD برپایة بررسی تعادل میان شمار بلور‏‌ها در یک طیف اندازة مشخص در حجم مشخصی از ماگماست که به نام چگالی تجمعی از آن یاد می‏‌شود. پلاژیوکلاز از فراوان‏‌ترین کانی‏‌‏‌های سازندة سنگ‏‌‏‌های آذرین است و در بیشتر شرایط گوناگون ماگمایی حالت پایدار دارد و نوسان‌های فیزیکو- شیمیاییِ آشیانة ماگمایی را در خود ثبت کند (Gagnevin et al., 2007; Ruprecht and Wörner, 2014). ازاین‌رو، در بیشتر بررسی‌های CSD انجام‌شده در سنگ‏‌‏‌های آذرین، بیشترین تأکید روی اندازه‏‌گیری پارامتر‏‌ها دربارة بلور‏‌های پلاژیوکلاز بوده است.

 Randolph و Larson (1971) از نخستین پژوهشگرانی هستند که به معرفی کاربرد روش توزیع اندازة بلور‏‌ها در گدازه‏‌‏‌ها پرداختند. سپس Marsh (1988) این روش را به‌صورت کاملاً علمی توسعه‌ داد. برپایة بررسی‌های وی، Higgins (2000) نرم‏‌افزار CSD Correction را طراحی کرد.

با تکمیل تئوری CSD روشن شد سرعت رشد بلور‏‌ها به اندازه آنها وابسته نیست؛ اما میزان هسته‏‌بندی بلور‏‌ها دربرابر اندازه آنها روندی وارونه نشان می‌دهد؛ یعنی هرچه میزان هسته‏‌بندی سریع‏‌تر باشد، بلور‏‌های پدیدآمده اندازة کوچک‏‌تری خواهند داشت. اندازه و آرایش فضایی بلور‏‌ها در سیستم‏‌‏‌های ماگمایی نشان‏‌دهندة شرایط بالاآمدن ماگماست، بدین‌گونه‌که میزان میکرولیت‏‌‏‌ها و میکروفنوکریست‏‌‏‌ها در سنگ‏‌‏‌های آذرین بیرونی نشان‌دهندة رشد بلور‏‌ها در مجاری و جمعیت فنوکریست‏‌‏‌ها نشان‏‌دهندة توقف ماگما در آشیانة ماگمایی است.

در این پژوهش، برپایة بررسی توزیع اندازة بلور‏‌های (CSD) پلاژیوکلاز در 10 نمونه از سنگ‏‌‏‌های آتشفشانی محدودة کودگان و تلفیق این اطلاعات با داده‌های به‏‌دست‏‌آمده‏‌ از تجزیة زمین‏‌شیمیایی نمونه‏‌‏‌ها، شرایط تبلور ماگمای سازنده سنگ‏‌‏‌های منطقه را بازسازی و فرایند‏‌های دخیل در تحولات ماگمایی آنها را بررسی می‌شود.

 

زمین‏‌شناسی عمومی

مجموعة کودگان در 90 کیلومتری جنوب خوسف، در استان خراسان جنوبی و خاور پهنه ساختاری لوت در پهنة ایران مرکزی جای دارد (شکل 1). بلوک لوت روند شمالی- جنوبی دارد و با گسل‌های نهبندان در خاور، درونه در شمال، نایبند در باختر و جازموریان در جنوب فراگرفته شده است (Berberian and King, 1981).


 

 

 

شکل 1- جایگاه محدودة کودگان (دایرة سبز رنگ) روی نقشة پهنه‏‌‏‌های ساختاری ایران، برگرفته ‏از Nowroozi (1976)، Aghanabati (2004) و Hessami و همکاران (2006)

 

 

برخی پژوهش‏‌‏‌های انجام‌شده روی توالی افیولیتی در مرز خاوری پهنة لوت با بلوک افغان و فعالیت‏‌‏‌های ماگماتیسم دیده‌شده در بلوک افغان، فرورانش پوستة اقیانوسی به زیر بلوک افغان را پیشنهاد کرده‏‌اند (Saccani et al., 2010)؛ اما برخی پژوهش‏‌‏‌های دیگر برپایة حجم و پراکندگی زمانی و مکانی ولکانیسم و پلوتونیسم در بلوک لوت، فرورانش پوستة اقیانوسی به‏‌سوی باختر و به زیر بلوک لوت را پیشنهاد کرده‏‌‏‌اند (Eftekharnezhad, 1981). بررسی‏‌‏‌های زمین‏‌شناسی و ایزوتوپی جدید روی برخی توده‏‌‏‌های آذرین درونیِ بلوک لوت، فرورانش دوسویه و نامتقارن پوستة اقیانوسی به زیر بلوک لوت و بلوک افغان با سرعت‏‌‏‌های گوناگون را نشان می‌‏دهند (Arjmandzadeh et al., 2011).

فعالیت ماگماتیسم گستردة بلوک لوت از زمان ژوراسیک آغاز و تا ترشیاری ادامه داشته و در این زمان به اوج خود رسیده است (Pang et al., 2013). سن‏‌سنجی بازالت‏‌‏‌های آلکالن میوسن تا کواترنری منطقه لوت- سیستان نشان‏‌دهندة برخورد بلوک لوت با بلوک افغان در زمان کرتاسه پسین است (Pang et al., 2013).

محدودة کودگان در بخش شمالی نقشة زمین‏‌شناسی 1:100.000 کودگان جای دارد. گروه‌‏‌های سـنگی در نقشة کودگان بیشتر از نوع آذرین هستند و سن آنها پالئوسـن تـا ائوسـن (Bolourian, 2010) است. ایــن سنگ‏‌ها دربردارندة تــوده‏‏‌های آذرین درونی ژرف و نیمه‏ژرف گوناگون، سنگ‏‏‌های آتشفشانی و آذرآواری هستند. کهن‏تـرین گروه‌های سنگی که گستردگی کمابیش بسیاری در جنوب منطقه دارند، واحـد‏‌های آگلـومرایی و تـوف هستند. بیشتر سطح محدوده را واحد‏‌های آندزیت (پیروکسن‌آندزیت) و دیوریت فراگرفته‌اند. نمونه‏‌‏‌های سنگی بررسی‌شده در این پژوهش، ترکیب آندزیت، پیروکسن‌آندزیت و بازالت دارند (شکل 2).

 

 

 

شکل 2- نقشة زمین‏‌شناسی محدودة کودگان


 


بیشتر گسل‏‌‏‌های محدودة کودگان روند شمال‏‌باختری- جنوب‏‌خاوری و خاوری- باختری نشان می‌دهند. این گسل‌ها در واحد آندزیتی دیده می‏‌شوند و شواهد رفتار گسل‌های نرمال و چپگرد در آنها دیده می‏‌شود.

 

سنگ‏‌نگاری

الف- آندزیت

آندزیت فراوان‏ترین ترکیب سنگ‏شناسی دیده‌شده در این محدوده است. بیشتر نمونه‏‏‌های برداشت‌شده از این واحد ترکیب آندزیت، پیروکسن‏آندزیت و هورنبلندآندزیت دارند. در این نمونه‏‏‌ها، بافت‏‌‏‌های پورفیری، گلومروپورفیری، پویی‏‌کلیتیک و غربالی از بافت‏‌‏‌های غالب به‌شمار می‌روند. پلاژیوکلاز فراوان‏‌ترین کانیِ این نمونه‏‌‏‌هاست و پس از آن هورنبلند، پیروکسن و الیوین که بیشتر به‏‌صورت شکل‌دار در نمونه‏‌‏‌ها دیده می‌شوند. نشانه‌های دگرسانی سریسیتی، اپیدوتی، کربناتی، بیوتیتی و اکسیداسیون آهن با درجات ضعیف تا کمی متوسط نیز دیده می‏‌شوند. بلور‏‌های پلاژیوکلاز به‏‌صورت شکل‏‌دار تا نیمه‏‌شکل‏‌دار هستند و در اندازه‏‌‏‌های متوسط تا درشت دیده می‌شوند. در این بلور‏‌ها، نشانه‌های ماکل پلی‏‌سینتتیک، زونینگ، بافت غربالی و پلاژیوکلاز دونسلی دیده می‏‌شوند. منطقه‏بندی ترکیبی در این بلور‏‌ها نشان‏‌دهندة تغییر ترکیب ماگما، در پی تزریق مداوم ماگما به آشیانة ماگمایی و تغییر شرایط پیرامون بلور‏‌های پلاژیوکلاز است (Holton et al., 2000) (شکل‌های 3- A تا 3- D و 3- F).

 

ب- بازالت

نمونه‏‏‌های بازالت دیده‌شده در این محدوده بسیار کم هستند؛‌ به‌گونه‌ای‌که برپایة شمار کم آنها واحد جدایی تفکیک نشده است. در این نمونه‏‏‌ها نیز مانند آندزیت‏‌‏‌ها، بافت‌های پورفیری، پویی‏کلیتیک و کمتر تراکیتی دیده می‌شوند. بیشتر بلور‏‌های الیوین دیده‌شده در این نمونه‏‏‌ها به سیلیس دگرسان شده و ایدنگزیتی شده‏اند. در این نمونه‏‏‌ها، نشانه‌هایی از دگرسانی سیلیسی و کربناته با درجه ضعیف تا متوسط دیده شده‌اند. وجود بافت‏‌‏‌های غیرتعادلی مانند منطقه‏‌بندی ترکیبی، بافت غربالی و پلاژیوکلاز دو نسلی معمولاً به متغیر‏‌های مستقل (فشار، دما و ترکیب شیمیایی) وابسته است که حالت تعادل پیشین سیستم را برهم می‏‌زند. این بافت‏‌‏‌ها در پی تغییرات حرارتی و شیمیایی مذابی پدید می‏‌آیند که در تماس با بلور است و چه‌بسا پیامد انتقال بلور به بخش دیگری از ماگما، جریان مذاب میان بلورهایی که زودتر پدید آمده‏‌اند، و جدایش پیشرونده مذاب باشد. دربارة چگونگی پیدایش بافت غربالی باور بر اینست که ناپایداری کم بخشی از بلور در مذاب، باعث ذوب و تبلور درجا (و دوباره) بخش‏‌‏‌های ذوب‌شده می‏شود (Stamatelopoulou et al., 1990; Tsuchiyama, 1985). این بافت‏‌‏‌ها ‏‌‏‌نشان‏‌دهندة نبود تعادل در آشیانة ماگمایی است که چه‌بسا پیامد ورود یا حضور دو یا چند ماگما در آشیانة ماگمایی باشد (شکل 3- C). پیروکسن و هورنبلند فراوان‏‌ترین کانی‏‌‏‌های مافیک دیده‌شده در واحد‏‌های سنگی آذرین بیرونی این محدوده هستند. بیشتر این بلور‏‌ها شکل‌دار هستند و در آنها دگرسانی‏‌‏‌های سریسیتی، کربناتی و اپیدوتی با درجه ضعیف تا متوسط دیده می‏‌شوند. بیشتر بلور‏‌های هورنبلند در نمونه، دگرسان ‌شده‌اند و حاشیة اپاسیتی دارند (شکل‌های 3- A، 3- C، 3- F و 3- E).


 

 

 

شکل 3- A) فنوکریست پلاژیوکلاز با نشانه‌های ضعیف زونینگ و کربناتی‏‌شدن فنوکریست‏‌‏‌های پیروکسن و پلاژیوکلاز (پیروکسن‌آندزیت)؛ B) بلور‏‌های شکل‌دار پلاژیوکلاز با ماکل پلی‏‌سینتتیک و الیوین سرپانتینیتی‌شده در زمینة ریزبلور و شیشه‏‌ای (آندزیت)؛ C) پلاژیوکلاز با بافت غربالی در کنار فنوکریست دگرسان‌شدة پیروکسن (پیروکسن‌آندزیت)؛ D) فنوکریست‏‌‏‌های پلاژیوکلاز با بافت غربالی، الیوین سیلیسی‌شده و ایدینگزیتی‌شده در زمینه‌ای از بافت میکرولیتی (بازالت)؛ E) بافت گلومروپورفیری با تجمع بلور‏‌های شکل‌دار پیروکسن (پیروکسن‌آندزیت)؛ F) زونینگ در فنوکریست پلاژیوکلاز در زمینة میکرولیتی (آندزیت)

 


زمین‏‌شیمی

در این بخش، برای بررسی و تکمیل نتایج به‌دست‌آمده از روش توزیع اندازة بلور‏‌ها و بررسی تحولات زمین‏‌شیمیایی مؤثر در پیدایش سنگ‏‌‏‌های این محدوده، از داده‌های به‌دست‌آمده از تجزیه به روش XRF و ICP-MS (برای 8 نمونة آذرین بیرونی و 7 نمونة آذرین درونی برداشت‌شده از محدودة کودگان) بهره گرفته شده است (جدول‌های 1 و 2).

برپایة جانمایی نمونه‏‌‏‌ها در نمودار‏‌های پیشنهادیِ Cox و همکاران (1979) و Middlemost (1994)، نمونه‏‌‏‌های آذرین بیرونیِ کودگان ترکیب آندزیت و بازالت دارند؛ اما نمونه‏‌‏‌های آذرین درونی در گسترة ترکیبی گرانیت، گرانودیوریت و کوارتزولیت جای می‏‌گیرند (زیرا دچار دگرسانی بالای سیلیسی شده‌اند) (شکل 4). همچنین، در نمودار‏‌های شناسایی سری‏‌‏‌های ماگمایی (شکل 5) نمونه‏‌‏‌های آذرین درونی سرشت توله‏‌ایتی دارند؛ اما نمونه‏‌‏‌های آذرین بیرونی روند کالک‏‌آلکالن نشان می‏‌دهند. ماگما‏‌های کالک‏آلکالن ویژة پهنه‌های فرورانشی هستند؛ اما سری‏‏‌های توله‏‌ایتی (که دربردارندة سنگ‏‏‌های حد واسط و اسیدی نیز هستند) گاه در پهنه‌های فرورانشی دیده شده‏اند. هر دو گروه، نمونه‏‌‏‌های آذرین بیرونی و درونی در جایگاه کمان آتشفشانی مرز قاره پدید آمده‏‌اند (شکل‌های 6 و 7).

 

 

جدول 1- داده‌های به‌دست‌آمده از تجزیة نمونه‏‌‏‌های کودگان به روش XRF برپایة درصدوزنی

Sample No.

97-KOD-85

97-KOD-54

97-KOD-340

97-KOD-56

97-KOD-5

97-KOD-11

97-KOD-309

 

Intrusive Rocks

SiO2

73.62

74.04

75.59

76.87

78.57

80.60

86.11

 

TiO2

0.47

0.47

0.64

0.41

0.45

0.68

0.43

 

Al2O3

11.94

14.76

4.75

11.07

10.64

13.06

5.05

 

Fe2O3

3.56

2.55

2.27

2.91

4.40

2.36

2.20

 

MnO

0.08

0.09

0.05

 

MgO

2.05

0.32

0.98

2.24

1.61

0.19

0.43

 

CaO

4.91

5.27

10.69

5.10

3.32

1.63

5.02

 

Na2O

0.22

0.37

0.23

0.37

0.46

0.21

0.20

 

K2O

0.11

0.43

0.10

0.11

0.16

0.19

0.12

 

P2O5

0.20

0.16

0.07

0.13

0.19

0.26

0.07

 

SO3

2.93

1.56

4.52

0.68

0.14

0.75

0.31

 

Sample No.

97-KOD-86

97-KOD-233B

97-KOD-147

97-KOD-361

97-KOD-401

97-KOD-13

97-KOD-394

97-KOD-196

Volcanic Rocks

SiO2

54.79

54.79

55.57

55.92

56.95

56.36

56.13

56.28

TiO2

0.75

0.76

0.67

0.72

0.76

0.71

0.76

0.66

Al2O3

17.33

17.62

16.58

16.88

18.27

17.69

16.51

15.94

Fe2O3

9.60

8.22

8.28

8.37

7.50

8.10

8.25

7.79

MnO

0.15

0.08

0.11

0.09

0.06

0.13

0.16

0.15

MgO

3.79

3.86

3.82

3.36

2.89

3.60

2.99

4.44

CaO

8.43

9.01

8.99

8.49

8.17

7.96

7.93

8.71

Na2O

2.27

2.67

2.49

2.50

2.73

2.77

3.11

2.33

K2O

2.50

2.55

3.01

3.17

2.37

2.31

3.50

3.20

P2O5

0.24

0.33

0.36

0.36

0.24

0.24

0.42

0.31

SO3

0.08

0.06

0.06

0.06

0.05

0.14

0.12

 

جدول 2- داده‌های به‌دست‌آمده از تجزیة نمونه‏‌‏‌های کودگان به روش ICP-MS (برپایة ppm)

Sample No.

97-KOD-5

97-KOD-11

97-KOD-13

97-KOD-54

97-KOD-56

97-KOD-85

97-KOD-86

97-KOD-147

Li

46

41

13

33

36

31

10

8

Be

0.9

0.6

1

1

0.9

0.9

1.1

1.3

Na

1776

611

19130

1642

1918

988

15882

16618

Mg

8273

1110

17510

1807

10786

10042

17787

18294

Al

46455

55857

79697

60952

46451

48266

76527

73320

P

499

802

783

453

384

602

805

1134

S

661

2698

57

5964

2332

12027

220

138

K

946

1171

16493

2818

659

630

17719

20352

Ca

18223

8836

44131

28762

27562

25886

46847

49440

Sc

9.6

14.7

20.9

7.9

7

9.4

27

23.4

V

72

184

171

64

66

108

202

187

Cr

27

102

24

10

20

32

16

48

Mn

155

271

762

409

523

118

901

708

Fe

25574

14964

47361

15486

16902

20719

54061

48384

Co

11

20

21

6

8

15

22

26

Ni

7

29

10

4

7

24

5

19

Cu

1

52

30

<1

4

16

<1

21

Zn

46

27

87

43

39

42

84

80

As

6.8

>100

0.8

9.5

5.3

12.1

5.5

0.4

Se

2.1

1.8

2.3

2.1

3.9

3.6

1.6

1.3

Rb

<1

<1

41

<1

<1

<1

49

64

Sr

619

667

407

330

441

988

428

595

Y

8.5

8.8

13.1

12.5

9.1

8

12.5

11.7

Zr

80

88

76

55

105

90

69

75

Nb

4.1

3.9

5.1

6.9

4

5.2

4.8

2.6

Mo

0.1

0.4

0.1

0.6

0.4

0.4

<0.1

0.1

Ag

0.1

0.4

0.3

<0.1

0.7

0.2

0.2

<0.1

Cd

<0.1

1.5

<0.1

0.1

<0.1

0.1

0.1

<0.1

In

<0.5

<0.5

<0.5

<0.5

<0.5

<0.5

<0.5

<0.5

Sn

1.3

1

1.1

0.6

1.3

1.2

1.1

1.2

Sb

0.5

32.5

0.5

2.4

0.8

0.9

<0.5

<0.5

Te

0.2

0.1

0.2

0.2

0.1

0.1

0.1

0.2

Cs

0.9

0.9

1.7

0.9

0.7

0.7

3.8

1.5

Ba

384

98

435

233

220

354

423

550

Tl

<0.1

0.2

0.1

<0.1

<0.1

0.4

0.3

0.2

Pb

14

17

18

8

9

15

16

17

Bi

0.3

0.2

0.2

<0.1

0.2

0.3

0.1

0.1

Ti

2276

3592

3767

2206

1954

2322

3928

3558

La

14

17

14

21

15

16

15

18

Ce

26

30

25

36

31

29

26

31

Pr

3.2

3.9

3.3

4.6

3.1

3.6

3.5

4.3

Nd

13

16.3

14.3

16.6

12.7

12.9

14.8

18.9

Sm

2.5

3.1

3.2

3.2

2

2.4

3.3

3.8

Eu

0.7

1

1.1

0.8

0.6

0.7

1.3

1.4

Gd

2.1

2.5

3

2.6

1.8

2.1

2.9

3.4

Tb

0.4

0.5

0.6

0.5

0.4

0.4

0.6

0.6

Dy

2.6

2.9

4

3.5

2.8

2.6

4

3.8

Er

1.8

1.6

2.5

2.3

1.8

1.6

2.4

2.3

Tm

0.3

0.2

0.4

0.4

0.3

0.3

0.4

0.4

Yb

1.2

1.4

2.3

1.6

1.3

1.2

2.3

2.1

Lu

0.2

0.2

0.3

0.3

0.2

0.1

0.3

0.3

Hf

2.4

2.5

2.3

1.8

2.6

2.4

2.2

2.4

Ta

0.6

0.5

0.6

0.9

0.6

0.6

0.6

0.5

W

<1

5.8

37.7

61.3

1.2

6.8

<1

8.1

Th

4.5

6.8

3.3

11

5

5.3

3.3

4.4

U

1.1

3.6

1.1

2.2

9

14.4

1.1

1.3

جدول 2- ادامه

Sample No.

97-KOD-196

97-KOD-233B

97-KOD-309

97-KOD-340

97-KOD-361

97-KOD-394

97-KOD-401

Li

11

7

35

31

8

7

10

Be

1.2

1.1

0.9

0.7

1.4

1.5

1.1

Na

15731

18079

846

1098

17231

20628

19104

Mg

19938

18016

2446

5007

15954

13600

13845

Al

70003

75655

22706

21053

76630

69050

81855

P

1028

1072

327

217

1218

1352

799

S

606

104

2022

17037

216

646

88

K

22427

17373

732

618

22379

22641

16800

Ca

47490

48205

27402

55866

48418

40486

45529

Sc

23.1

18.5

5.2

5.4

24.4

17.3

21.3

V

164

165

64

57

190

79

176

Cr

75

94

23

15

50

26

19

Mn

812

507

245

217

550

923

369

Fe

43564

46089

13430

13395

48015

43381

41979

Co

21

19

9

7

22

18

20

Ni

19

17

4

6

16

11

9

Cu

19

12

<1

<1

28

7

<1

Zn

82

67

36

35

97

83

79

As

0.5

<0.1

35.8

3.3

2.9

3.2

<0.1

Se

1.9

2.4

2.3

3.4

2.9

2.1

2

Rb

64

40

<1

<1

70

64

40

Sr

566

493

284

475

660

603

448

Y

10.7

12.7

7.2

6.8

11.7

11.1

12.4

Zr

73

69

32

66

79

65

75

Nb

5

4

<1

1.5

4.7

2.3

4.3

Mo

<0.1

<0.1

<0.1

0.3

<0.1

<0.1

<0.1

Ag

0.2

0.1

<0.1

0.1

0.2

0.1

0.2

Cd

<0.1

0.1

0.3

<0.1

0.1

<0.1

0.1

In

<0.5

<0.5

<0.5

<0.5

<0.5

<0.5

<0.5

Sn

1.2

1

0.7

1.1

1.1

1.3

1.1

Sb

<0.5

<0.5

4.6

2.4

1

<0.5

<0.5

Te

0.2

0.1

0.3

0.2

0.2

0.3

0.2

Cs

1.5

1.2

0.8

0.8

2

1.7

2

Ba

479

455

449

1084

540

560

438

Tl

0.3

0.2

<0.1

<0.1

0.2

<0.1

0.2

Pb

16

13

6

4

18

17

15

Bi

0.2

0.2

0.1

0.2

0.1

0.1

0.2

Ti

3373

3924

2149

3133

3717

3654

3872

La

17

16

7

11

19

16

14

Ce

29

28

10

19

34

30

25

Pr

4.1

3.9

1.2

2.5

5

4.6

3.3

Nd

17.3

17.6

5.3

8.3

21

19.6

14.4

Sm

3.6

3.7

0.9

1.6

4.7

4.1

3.1

Eu

1.2

1.4

0.4

0.5

1.5

1.5

1.2

Gd

3.2

3.5

0.8

0.5

3.9

3.7

3.4

Tb

0.5

0.7

0.4

0.3

0.6

0.6

0.6

Dy

3.2

4

2.3

2.3

4

3.6

3.7

Er

2.1

2.5

1.5

1.5

2.2

2.1

2.5

Tm

0.4

0.4

0.3

0.2

0.4

0.3

0.4

Yb

1.9

2.2

1

0.9

2.1

1.6

2.2

Lu

0.2

0.3

0.1

0.1

0.3

0.2

0.3

Hf

2.3

2.1

0.8

1.9

2.6

2.2

2.5

Ta

0.7

0.6

0.3

0.5

0.6

0.5

0.6

W

<1

<1

<1

<1

65.5

3.3

3.1

Th

3.9

3.1

1

3.3

4.9

3.1

3.1

U

1.3

1.2

0.6

0.5

1.4

1

1.1

 

 

 

 

 
 

 

 

 

شکل 4- ترکیب نمونه‏‌‏‌های آذرین محدودة کودگان (●: سنگ‌های آذرین بیرونی؛ ▲ : سنگ‌های آذرین درونی) در: A) نمودار Cox و همکاران (1979)؛ B) نمودار Middlemost (1994)

 

 

شکل 6- ترکیب نمونه‏‌‏‌های آذرین درونی محدودة کودگان در: A) نمودار Yb دربرابر Ta Pearce et al., 1984))؛ B) نمودار Y دربرابر Nb Pearce et al., 1984))

 

 

 

 

 

 

 

 

شکل 5- ترکیب نمونه‏‌‏‌های آذرین محدودة کودگان در نمودار سه‏‌تایی (Irvine and Baragar, 1971)

 

 

شکل 7- ترکیب نمونه‏‌‏‌های آذرین محدودة کودگان در نمودار Nb/Yb دربرابر Th/Yb Hollocher et al., 2012))

 


عنصرهای اندازه‏‌گیری‌شده در نمونه‏‌‏‌های بررسی‌شده دربرابر ترکیب گوشتة طبیعی (شکل 8) و گوشتة اولیه (شکل 9) و برپایة ترکیب‌های پیشنهادیِ Sun و McDonough (1989)، بهنجارسازی شده‎اند. غنی‏‌شدگی از LIEE در کنار آنومالی منفی P، Ti، Ba و Nb چه‌بسا پیامد ذوب‏‌بخشی سنگ‏‌های پوستة قاره‏‌ای در کمان ماگمایی مرتبط با فرورانش (Kolb et al., 2013) و یا ذوب گوة ‏‌گوشته‏‌ای غنی‏شده با سیال‌های متاسوماتیسم‌کننده است (Girardi et al., 2012).

حضور نا‏‌هنجاری منفی Nb از ویژگی‌های پهنه‌های فرورانش و مرز‏‌‏‌های فعال قاره‏‌ای و همچنین، از نشانه‏‌‏‌های شرکت پوسته در فرایند‏‌های ماگمایی است (Wilson, 1989). آنومالی منفی Nb و Ti نشان‏دهندة: ویژگی ماگماتیسم مرتبط با فرورانش (Rudnick and Gao, 2003)، مشارکت پوسته در فرایند‏‌های ماگمایی (Nagudi et al., 2003) و آلودگی پوسته یا ذوب دوبارة آن (Swain and Kirby, 2008) است.

آنومالی مثبت Pb گویای متاسوماتیسم گوة ‏‌گوشته‏‌ای با سیال‌های پدیدآمده از سنگ‏‌کرة اقیانوسی فرورو و یا آلایش ماگما با سنگ‏‌کرة قاره‏‌ای است (Atherton and Ghani, 2002). در پهنه‌های فرورانش، به‌دنبال تأثیر سیال‌های پدیدآمده از تختة فرورو، نسبت LILE به HFSE افزایش می‏‌یابد. در حقیقت، این ناهنجاری در سنگ‏‌‏‌های آتشفشانی کالک‏‌آلکالن در پی افزوده‌شدن عنصرهای LILE از تختة فرورو به درون گوة گوشته‌ای بالای آن روی می‏‌دهد.

آنومالی منفی P نشان‏‌دهندة رفتار سازگار فسفر هنگام روند جدایش بلورین کانی‌هاست؛ زیرا این عنصر وارد کانی آپاتیت می‌شود.

 

 

 

شکل 8- نمودار بهنجارسازی عنصرهای کمیاب و خاکی کمیاب دربرابر ترکیب گوشتة طبیعی Sun and McDonough, 1989)) برای نمونه‏‌‏‌های محدودة کودگان

 

شکل 9- بهنجارسازی عنصرهای کمیاب و خاکی کمیاب دربرابر ترکیب گوشتة اولیه Sun and Mc Donough, 1989)) نمونه‏‌‏‌های سنگی محدودة کودگان

 


به‌دست‌آوردن پارامتر‏‌های آماری و رسم نمودار‏‌های توزیع داده‏‌ها

دست‏‌یابی به پارامتر‏‌های آماری، نخستین گزینه‌ای است که به‏‌عنوان مهم‏‌ترین و جامع‏‌ترین اطلاعات آماری در پیش روی داده‏‌پردازان است. این اطلاعات، چگونگی پراکندگی داده‏‌‏‌ها حول میانگین و تمایل به میانگین، میزان چولگی و کشیدگی جامعة عیار نمونه‏‌‏‌ها و تشابه یا تمایز آن با یک توزیع نرمال را به نمایش می‏‌گذارند. همة داده‏‌‏‌های خام به‏‌دست‏‌آمده‏‌ از تجزیة اکسید عنصر‏‌های اصلی در واحد‏‌های سنگی آذرین درونی و بیرونی محدودة کودگان، پس از گروه‏‌بندی مناسب با نرم‏‌افزار SPSS25 بررسی‏‌ آماری شدند. برای این کار، داده‌های به‏‌دست‏‌آمده‏‌ از تجزیة 8 نمونة برداشت‌شده از سنگ‌های آذرین بیرونی و 7 نمونة برداشت‌شده از سنگ‌‏‌های آذرین درونی به‌کار برده شدند. برپایة سرشت متفاوت آنها، پارامتر‏‌های آماری برای هر گروه از واحد‏‌های سنگی جداگانه به‌دست آورده شدند.

در نمونه‏‌‏‌های آذرین بیرونی برداشت‌شده از این محدوده، 12 اکسید عنصر اصلی زیر تجزیه و پردازش شدند:

SiO2, TiO2, Al2O3, Fe2O3, MnO, MgO, CaO, K2O, Na2O, P2O5, BaO, SO3

میانگین، انحراف معیار، چولگی، کشیدگی، و کمترین و بیشترین مقدار داده‏‌‏‌ها از پارامتر‏‌های آماری داده‏‌‏‌های خام اکسید‏‌های عنصرهای گوناگونی است که در جدول 3 نشان داده شده‌اند. گفتنی است این مقدارها و بازه‏‌‏‌های عددی مرتبط با هر عنصر، با نرم‏‌افزار برگزیده و به‌دست آورده می‏‌شوند و سپس در نمودار‏‌های هیستوگرام رسم و نشان داده می‏‌شوند. در شکل‌های 10 و 11، نمودارهای هیستوگرام داده‏‌‏‌های خام عنصر‏‌های اصلی آورده شده‌اند.

 


جدول 3- پارامتر‏‌های آماری داده‏‌‏‌های خام اکسید‏‌های عنصرهای گوناگون در نمونه‏‌‏‌های آذرین بیرونی محدودة کودگان

Statistics

SiO2

TiO2

Al2O3

Fe2O3

MnO

MgO

N

Valid

8

8

8

8

8

 

Missing

0

0

0

0

0

 

Mean

77.9147

0.5066

10.1828

2.8941

0.0527

1.1169

Std. Deviation

4.37373

0.10900

3.85648

0.81311

0.02225

0.85197

Skewness

1.168

1.154

- 0.689

1.271

1.073

0.271

Kurtosis

1.178

- 0.579

- 1.001

0.756

- 0.644

- 2.063

Minimum

73.62

0.41

4.75

2.20

0.04

0.19

Maximum

86.11

0.68

14.76

4.40

0.09

2.24

Statistics

CaO

Na2O

K2O

P2O5

SO3

BaO

N

Valid

8

8

8

 

8

8

Missing

0

0

0

 

0

0

Mean

5.1337

0.2950

0.1755

0.1533

1.5552

0.0617

Std. Deviation

2.78372

0.10168

0.11801

0.06847

1.61258

0.04706

Skewness

1.321

0.719

2.231

0.064

1.256

2.372

Kurtosis

3.249

- 1.266

5.201

- 0.796

0.653

5.816

Minimum

1.63

0.20

0.10

0.07

0.14

0.04

Maximum

10.69

0.46

0.43

0.26

4.52

0.17

 

 

شکل 10- نمودارهای هیستوگرام برای داده‏‌‏‌های خامِ درصدوزنی اکسید‏‌های Al2O3، BaO، CaO، Fe2O3 و K2O در نمونه‏‌‏‌های آذرین بیرونی محدودة کودگان

 

شکل 11- نمودارهای هیستوگرام داده‏‌‏‌های خامِ درصدوزنی اکسید‏‌های MnO، Na2O، P2O5، SiO2 و TiO2 در نمونه‏‌‏‌های آذرین بیرونی محدودة کودگان

 

 

با نگاهی به نمودار هیستوگرام داده‏‌‏‌های خام عنصرها، دیده می‏‌شود اکسید‏‌های عنصر‏‌های اصلی در واحد‏‌های سنگی آذرین بیرونی پراکندگی نرمال دارند. تنها در 2 نمونه، اکسید‏‌های BaO و K2O به اندازة بسیار کم چولگی مثبت نشان می‏‌دهند؛ ازاین‌رو، از نرمال‏‌سازی این اکسید‏‌ها چشم‌پوشی می‌شود.

از 8 نمونة برداشت‌شده از سنگ‏‌‏‌های آذرین درونی نیز مقدار اکسید عنصر‏‌های اصلی (همانند سنگ‌های آذرین بیرونی) در 7 نمونه تجزیه شده است و پارامتر‏‌های آماری داده‏‌‏‌های خام این متغیر‏‌ها در جدول 4 نشان داده شد‌ه‌اند. نمودارهای هیستوگرام داده‏‌‏‌های خامِ اکسید‏‌های عنصر‏‌های اصلی در شکل‌های 12 و 13 نشان داده شده‌اند.

نمودارهای هیستوگرامِ داده‏‌‏‌های خامِ اکسید‏‌های عنصرها برای نمونه‏‌‏‌های آذرین درونی نشان می‌دهند این متغیر‏‌ها در سنگ‏‏‌های آذرین درونی پراکندگی نرمال دارند و تنها در مقدار Fe2O3 به اندازة بسیار کم چولگی مثبت دیده می‏‌شود و ازاین‌رو، از نرمال‏‌سازی مقدار این اکسید در نمونه‏‌‏‌ها چشم‌پوشی شد.

 


جدول 4- پارامتر‏‌های آماری داده‏‌‏‌های خام اکسید‏‌های عنصرهای گوناگون در نمونه‏‌‏‌های آذرین درونی محدودة کودگان

Statistics

SiO2

TiO2

Al2O3

Fe2O3

MnO

MgO

N

Valid

7

7

7

7

7

 

Missing

0

0

0

0

0

 

Mean

55.8508

0.7239

17.1018

8.2623

0.1173

3.5947

Skewness

- 0.359

- 0.606

0.010

1.499

- 0.301

0.087

Kurtosis

- 0.731

- 1.385

- 0.728

3.737

- 1.404

- 0.164

Minimum

54.79

0.66

15.94

7.50

0.06

2.89

Maximum

56.95

0.76

18.27

9.60

0.16

4.44

Statistics

CaO

Na2O

K2O

P2O5

SO3

BaO

N

Valid

7

7

7

 

7

7

Missing

0

0

0

 

0

0

Mean

8.4618

2.6093

2.8267

0.3120

0.0772

0.0670

Skewness

0.075

0.667

0.255

0.285

1.039

1.700

Kurtosis

- 1.496

0.456

- 1.725

- 0.939

- 0.188

2.711

Minimum

7.93

2.27

2.31

0.24

0.04

0.06

Maximum

9.01

3.11

3.50

0.42

0.14

0.09

 

 

شکل 12- نمودارهای هیستوگرام داده‏‌‏‌های خامِ درصدوزنی اکسید‏‌های Al2O3، BaO، CaO، Fe2O3 و K2O در نمونه‏‌‏‌های درونی محدودة کودگان

 

شکل 13- نمودارهای هیستوگرام داده‏‌‏‌های خامِ درصدوزنی اکسید‏‌های MnO، Na2O، P2O5، SiO2 و TiO2 در نمونه‏‌‏‌های آذرین درونی محدودة کودگان

 

 

آنالیز خوشه‏‌ای

آنالیز خوشه‌ای ابزاری قدرتمند برای شناخت و برگزیدن گروه‏‌‏‌های همگن از داده‏‌‏‌های خام است. این روش، برپایة شماری از الگوریتم‏‌‏‌های گوناگون، داده‏‌‏‌های خام را به‌ترتیب میزان ارتباط آنها با یکدیگر رده‏‌بندی می‏‌کند. نتیجه این آنالیز، نموداری است که نمودار درختی خوانده می‏‌شود.

نمودار درختی به‏‌دست‏‌آمده‏‌ از آنالیز نمونه‏‌‏‌های آذرین درونی و بیرونی محدودة کودگان در شکل 14 آورده شده است. همان‏‌گونه‏‌که دیده می‏‌شود، نمونه‏‌‏‌های بررسی‌شده در دو دستة جداگانه A و B جای می‏‌گیرند. خوشة A دربردارندة نمونه‏‌‏‌های آذرین درونی و خوشة B دربرگیرندة نمونه‏‌‏‌های آذرین بیرونی است. ازاین‏‌رو، از یکدیگر جدا هستند. این پدیده شاید پیامد تفاوت رفتار اکسید عنصر‏‌های اصلی در این دو گروه از واحد‏‌های سنگی باشد که جایگیری آنها در دو گروه متفاوت را به‌دنبال داشته است.

 


 

شکل 14- نمودار درختی عنصرهای اندازه‏‌گیری‌شده در نمونه‏‌‏‌های لیتوژئو‏‌شیمیایی سیستماتیک محدودة کودگان

 


روش کار (انجام آنالیز توزیع اندازة بلور)

پس از بررسی‏‌‏‌های صحرایی و نمونه‏‌برداری از سنگ‏‌‏‌های آذرین بیرونی منطقة کودگان، بررسی‌های سنگ‏‌نگاری روی 107 مقطع نازک انجام شد. برای شناسایی تاریخچة تبلور، میزان رشد، زمان توقف ماگمای سازنده و تحولات ماگمایی تودة آذرین درونی، ابعاد بلور‏‌های پلاژیوکلاز در 10 نمونه با سنگ میزبان آندزیت و بازالت اندازه‏‌گیری شد. داده‏‌‏‌‏‌های به‌دست‌آمده که دربردارندة درازا، پهنا و مساحت بلور‏‌های پلاژیوکلاز هستند به نرم‏‌افزار CSD correction 1.3 انتقال یافتند. با به‌کارگیری این نرم‏‌افزار نمودار‏‌های سه‏‌بعدی لازم برای تعبیر و تفسیر داده‏‌‏‌های CSD تهیه شدند. برای تبدیل داده‏‌‏‌های دوبعدی به سه‏‌بعدی از نسبت ابعاد بلور، بعد کوتاه (S)، بعد متوسط (I)، بعد بلند (L) و فاکتور گردشدگی (بیضوی= 1، مکعبی=0) بهره گرفته می‏‌شود.

انواع نمودار‏‌های CSD عبارتند از: 1) نیمه‏‌لگاریتمی یا کلاسیک ؛2) لگاریتم طبیعی؛ 3) فرکتال (مستقل از مقیاس) یا ترکیبی از آنها (Higgins, 2006b) (شکل 15).

 

 

شکل 15- انواع نمودار‏‌های CSD (Higgins, 2006b)

داده‌های به‌دست‌آمده از تجزیة زمین‏‌شیمیایی نمونه‏‏‌های آذرین بیرونی و درونی کودگان به روش XRF و ICP-MS برای شناخت فرایند‏‌های مؤثر در پیدایش این سنگ‌ها به کار برده شدند.

 

توزیع اندازة بلور‏‌ها (CSD)

نتیجة آنالیز توزیع اندازة بلور‏‌ها (CSD) به‏‌صورت نمودارهایی نمایش داده می‏‌شود که در آنها شمار بلور‏‌ها در یک واحد حجمی (n)، به‏‌عنوان تابعی از اندازة بلور (L) هستند. برای به‌دست‌آوردن شیب و محل برخورد از روند خطی تابع ریشة دوم ln(n) دربرابر L (طول) بهره گرفته می‏‌شود. برای به‌دست‌آوردن سهم اشغال‌شده توسط کانی‏‏‌ها، فرمول سوم به‌کار برده می‏شود:

16SLOP= -1/Gt">             16Intercept=I/G">

V1=6σno1C14

برای بررسی توزیع اندازة بلور‏‌های پلاژیوکلاز از سنگ‏‏‌های آذرین بیرونی (آندزیت و بازالت) محدودة کودگان، تصویرهای میکروسکوپی دوبعدی با بزرگنمایی 40 از مقطع‌های نازک واحد میزبان تهیه شده است. این عکس‏‌‏‌ها در نرم‏‌افزار تخصصی JMicroVision 1.3.1 برپایة بزرگنمایی میکروسکوپ کالیبره شده‌اند و سپس ابعاد 1739 بلور پلاژیوکلاز اندازه‏‌گیری شده است. داده‏‌‏‌‏‌های دوبعدی به‏‌دست‏‌آمده‏‌ از این اندازه‏‌گیری‏‌‏‌ها که شامل درازا، پهنا و مساحت بلور‏‌ها هستند به نرم‏‌افزار CSDcorrection 1.3 انتقال و با به‌کارگیری رابطة n=dn/dl (dn: شمار بلور‏‌ها؛ dl: اندازة بلور‏‌‏‌ها)، چگالی تجمعی دانه‏‌‏‌ها (n) به‌دست ‌آورده ‌شد. در پایان، نمودار‏‌های کلاسیک تغییرات اندازه دانه‏‌‏‌‏‌ها (mm) دربرابر چگالی تجمعی (mm-4) رسم شدند (شکل 16).

 

 

 

شکل 16- فرایند‏‌های مؤثر در شکل نمودار‏‌های CSD برگرفته از Higgins (b2006) و تغییراتی از Vient و Higgings (2010)

 

 

دربارة نمونه‏‌‏‌های بررسی‌شده، همة نمودار‏‌های کلاسیک تغییرات اندازة بلور‏‌ها دربرابر چگالی تجمعی شکستکی و تقعر به‏‌سوی بالا نشان می‌دهند (شکل 16). پیچیدگی یا خمیدگی در نمودار‏‌های CSD شاید پیامد فرایندهایی مانند انباشتگی یا کاهش شمار بلور‏‌ها (Marsh, 1998)، تراکم (Boorman et al., 2004)، آمیختگی دو سری جمعیت بلوری یا چندین ماگمای مختلف (Higgins, 1996; Jerram et al., 2003)، درشت‏‌شدگی بافتی (Higgins and Roberge, 2003) و یا همزمانی فرایند درشت‏‌شدگی بافتی و تراکم پدیدآمده در پی تبلور دوباره باشد (Higgins, 2002b) (شکل 17).

درشت‏‌شدگی بافتی فرایندی است که بلور‏‌های کوچک حل می‌شوند و همزمان بلور‏‌های درشت‏‌تر رشد می‏‌کنند؛ ازاین‏‌رو، انرژی سطح کاهش می‏‌یابد (Voorhees, 1992).

برپایة ویژگی‌های سنگ‏‌نگاری، زمین‏‌شیمیایی و توزیع اندازة بلور در نمونه‏‏‌های بررسی‌شده، فرایند ورود و حضور چندین ماگمای مختلف در آشیانة ماگمایی و آمیختگی، شکستگی و خمیدگی در نمودار‏‌های CSD را به دنبال دارند.

 

 

شکل 17- نمودار‏‌های لگاریتم طبیعی اندازه‏‌گیری‌شده برای بلور‏‌های پلاژیوکلاز در محدودة کودگان

 

 

به باور Higgins (a2002)، افزوده‏‌شدن دو نوع CSD صاف و مستقیم با شیب و عرض از مبدأ‏‌های متفاوت، یک CSD با شیب تند برای بلور‏‌های با اندازة کوچک و با شیب آرام برای بلور‏‌های بزرگ‏‌تر پدید می‏‌آورد. ازآنجایی‌که ترکیب این دو CSD خطی لگاریتمی‏ است، محور قائم یک CSD خمیده با دو بخش خطی پدید می‌آورد. حضور متعدد ماگما در آشیانه‏‌‏‌های ماگمایی، خمیدگی CSD را به‌دنبال دارد. برپایة اندازه‏‌‏‌های متفاوت در هر نسل بلوری مشخص می‌شود هر نسل دچار فشار‏‌های متفاوت بوده و ازاین‌رو، در ژرفای متفاوتی پدید آمده است.

همان‏‌گونه‏‌که در شکل 17 دیده می‏‌شود، CSD در کوچک‏‌ترین اندازة بلور‏‌‏‌‏‌ها، بیشترین مقدار را دارد و میزان آن با افزایش اندازة بلور‏‌ها کاهش می‏‌یابد. در سیستم‏‌‏‌های باز که با ورود مکرر ماگمای جدید به درون مخزن تایید شده است، شیب نمودار‏‌ها منفی است (Marsh, 1998). رویداد شیب منفی نشانة افزایش رشد همراه با کاهش هسته‏‌بندی است و افزایش اندازه و کاهش چگالی تجمعی بلور‏‌ها منجر را به‌دنبال دارد.

ازآنجایی‌که پلاژیوکلاز فراوان‏ترین کانی دیده‌شده در این نمونه‏‏‌هاست، برای به‌دست‌آوردن سهم اشغال‌شده با بلور‏‌های پلاژیوکلاز در مقطع‌ها، فرمولV1=6σno1C14 به‌کار برده می‏‌شود. در این فرمول، V سهم حجمی، n میزان هسته‏‌بندی و C میانگین درازای بلور‏‌ها هستند (Marsh, 1988). ازاین‌رو، در نمونه‏‌‏‌های گوناگون، سهم حجمی بلور‏‌های پلاژیوکلاز از 40/5 تا 02/65 متغیر است. زمان رشد و سرعت هسته‏‌بندی رابطة مستقیمی با درصدحجمی هر بلور دارد. هرچه سرعت هسته‏‌بندی و زمان رشد بلور افزایش یابد، بلور‏‌ها درصد بیشتری از حجم سنگ را فرا می‏‌گیرند.

به پیشنهاد Marsh (1988) برای به‌دست‌آوردن متوسط زمان رشد بلور‏‌ها، معادله زیر به‌کار برده می‏‌شود:

Tr=(- 1/G*m)/31536000

در این معادله Tr زمان توقف (سال)، G میزان رشد بلور، m شیب منحنی CSD (لگاریتم خطی) و 31536000 ضریب تبدیل ثانیه به سال است. انتخاب میزان رشد مناسب، فاکتور مهمی در به‌دست‌آوردن زمان توقف است. به پیشنهاد Cashman (1993)، این میزان برابربا 10- 9 mm/s و 10- 10 mm/s دانسته می‏‌شود. در بررسی CSD محدودة کودگان از میزان رشد 10- 10 mm/s برای انجام محاسبه‌ها بهره گرفته شد. اگر میزان هسته‏‌بندی و رشد برابر باشند، CSD افقی و بدون شیب است و هرچه اختلاف این دو بیشتر باشد، میزان شیب نیز بیشتر خواهد شد. برپایة محاسبه‌های انجام‌شده روی نمونه‏‌‏‌های محدودة کودگان، زمان توقف ماگما در آشیانه ماگمایی و تبلور بلور‏‌های پلاژیوکلاز، برپایة شیب به‏‌دست‏‌آمده‏‌ از نمودار‏‌ها، برابربا 79 تا 1141 سال به‌دست آمده است (جدول 5).

جدول 5- نتایج محاسبه‌های انجام‌شده برای بلور‏‌های پلاژیوکلاز

Sample number

Slop

Intercept

Time

G

mms- 1

9A

0.47

1.15

1489

10- 10

12

0.45

0.24

1425

10- 10

14

0.45

0.28

1425

10- 10

147

0.25

0.16

792

10- 10

347

0.6

0.17

1901

10- 10

368

0.2

0.23

633

10- 10

 

خاستگاه

هم‏‌روند‏‌بودن نمونه‏‌‏‌های یک مجموعه روی نمودار، احتمال وابستگی خویشاوندی میان ماگما‏‌ها و روند تحولی ماگما را نشان می‏‌دهد. پیوستگی این روند نشان می‌دهد نمونه‏‌‏‌های یادشده از یک ماگمای مشترک جدایش یافته‌اند و یا پیامد ذوب‏‌بخشی خاستگاهی یکسان هستند. از سوی دیگر، نبود پیوستگی نقاط و قطع‏‌شدگی روند‏‌ها روی نمودار، چه‌بسا نشانة نبود وابستگی خویشاوندی میان ماگما‏‌ها در یک مجموعه است.

نمودار‏‌های تغییرات اکسید‏‌های عنصر‏‌های اصلی و کمیاب دربرابر SiO2 که برای نمونه‏‌‏‌های آذرین بیرونی و درونی محدودة کودگان (شکل 18) رسم شده‌اند، نبود پیوستگی در نمونه‏‌‏‌های یاد شده را نشان می‌دهند؛ به‌گونه‌ای‌که آنها را در دو دستة کاملاً جدا جای می‏‌دهند. نبود پیوستگی و روند‏‌های متفاوت در نمودار‏‌های هارکر نشان‌دهندة وجود دو خاستگاه متفاوت برای سنگ‏‏‌های آذرین بیرونی و درونی این منطقه هستند.

برپایة بررسی‏‏‌های سنگ‏‌نگاری، زمین‏‌شیمیایی و توزیع اندازة بلور، این نمونه‏‏‌ها از خاستگاه ‏‌گوشته‏‌ای و پوستة آغشته‌شده به گوشتة بالایی پدید آمده باشند.

برپایة یافته‌های Yousefzadeh و همکاران (2019)، کانسار‏‌های در محدودة کودگان از نوع کانسار‏‌های IOCG ذخایر سولفیدی با خاستگاه گرمابی هسـتند و کانی‏سازی‌های گوناگونِ برشی، رگه و رگچه‏ای، افشـان و تـود‏‌ه‌ای در آنها دیده مـی‏شـوند. به باور Karimpour و همکاران (2011)، این کانسار‏‌ها از دیدگاه زمین‏ساختی در حوضه‏‏‌های کششی پشت کمربند پهنة فرورانش مرز قاره‏‌ای، ریفـت‏‌‏‌های درون قـاره، نقـاط داغ درون‌ قـاره و حوضــه‏‏‌های کششــی پــسابرخــوردی پدید مــی‏آیند. برپایة پــژوهش Omaraeei و Niroumand (2014)، سنگ‏‏‌های محدودة کودگـان از نـوع کلسـیمی قلیـایی و متـاآلومینوس است و از دیدگاه پهنة زمـین‌سـاختی، الگـوی ماگماتیسـم در ارتبـاط بـا کمان مرز قاره‏‌ای را نشان می‏دهند.

 

 

 

شکل 18- نمودار تغییرات درصدوزنی SiO2 دربرابر اکسید‏‌های اصلی

 

 

روند کاهشی MgO و MnO دربرابر SiO2 و کاهش Cr و Co با افزایش SiO2 روند عادی جدایش بلورین ماگمایی، و نیز نشان‏‌دهندة آلایش یا آمیختگی ماگمایی است (Kumar and Rino, 2006)، همبستگی ضعیف Rb و Sr با SiO2 و پراکندگی بالای آنها چه‌بسا پیامد واکنش با ترکیب‌های پوسته‏‌ای است (Kumar and Rino, 2006) (شکل 19).

نسبت‏‌‏‌های کم Ce/Pb (کمتر از 5) و Nb/La (کمتر از 7/0) و همچنین، نسبت بالای Sr/Ce (بیشتر از 5) نشانة درگیر‏‌بودن پوسته و گوشته در پیدایش ماگمای سازنده است (Jung et al., 2005). این نسبت‏‌‏‌ها در نمونه‏‌‏‌های محدودة کودگان به‌ترتیب برابربا 75/4- 3/1، 36/0- 1/0 و 5/34- 1/14 هستند. همة این‏‌ نسبت‌ها چه‌بسا نشان‏‌دهندة درگیر‏‌بودن پوسته و گوشته در پیدایش ماگمای مادر واحد‏‌های آذرین بیرونی و درونی پدیدآمده در محدودة کودگان باشد.

در نمودار‏‌ MgO دربرابر Fe2O3 و نیز نمودار Nb دربرابر Nb/Y، روند نمونه‏‌‏‌های بررسی‌شده در محدودة کودگان (شکل 20)، با رویداد فرایند آمیختگی ماگمایی همخوانی دارد.


 

 

 

شکل 19- روند تغییرات درصدوزنی SiO2 دربرابر Rb و Sr (برپایة ppm) برای نمونه‌های آذرین بیرونی و درونی کودگان

 

 

شکل 20- A) نمودار درصدوزنی MgO دربرابر Fe2O3 (Zhou et al., 1994)؛ B) نمودار Nb (برپایة ppm) دربرابر Nb/Y (Treuil and Joron, 1975)

 

 

برداشت

محدودة کودگان دربردارندة سنگ‏‌‏‌های آذرین بیرونی، نیمه‏‌ژرف با ترکیب آندزیت، دیوریت و بازالت و واحد‏‌های سنگی آذرین درونی با ترکیب گرانیت تا گرانودیوریت به سن ائوسن است.

واحد آندزیتی که فراوان‏ترین ترکیب در این محدوده را دارد، بازة ترکیبی پیروکسن‌آندزیت و هورنبلندآندزیت نشان می‌دهد. همچنین، ویژگی‌های مهمی مانند بافت غربالی، زونینگ و پلاژیوکلاز دونسلی، شکستگی و خمیدگی در نمودار‏‌های CSD و روند تغییرات Nb دربرابر Nb/Y در بیشتر نمونه‏‏‌های بررسی‌شده دیده شده‌اند.

در نمونه‏‏‌های بازالتی برداشت‌شده از این محدوده که فراوانی کمی نیز داشته‏اند، پلاژیوکلاز فراوان‏ترین کانی دیده‌شده است. همانند نمونه‏‏‌های آندزیتی، در بازالت‌ها نیز بافت غربالی و پلاژیوکلاز دونسلی دیده می‌شود. همچنین، در نمودار‏‌های CSD برای نمونه‏‏‌های بازالت نیز شکستکی و خمیدگی در پی تفاوت اندازة بلور‏‌ها دیده می‏شوند.

ویژگی‌های سنگ‏‌نگاری مانند بافت غربالی، زونینگ و پلاژیوکلاز دونسلی، نشانه‌ای از رویداد شرایط نامتعادل در آشیانة ماگمایی و نشان‏‌دهندة رخداد فرایند آمیختگی ماگمایی در نمونه‏‌ها هستند. حضور بلور‏‌های کوچک پلاژیوکلاز در درون پلاژیوکلاز‏‌های بزرگ‌تر چه‌بسا نشانة آشفتگی در آشیانة ماگمایی و حضور دو نسل گوناگون پلاژیوکلاز در آشیانة ماگمایی است که از دو ماگما در فاصلة زمانی متفاوت متبلور شده‌اند.

نبود پیوستگی در روند تغییرات اکسید‏‌های عنصر‏‌های اصلی و عنصر‏‌های کمیاب دربرابر SiO2 و نمودار خوشه‌ای در نمونه‏‌‏‌های آذرین بیرونی و درونیِ این محدوده، نشان‌دهندة پیدایش نمونه‏‌‏‌ها از دو خاستگاه متفاوت است.

روند‏‌های دیده‌شده در نمودار‏‌های Nb دربرابر Nb/Y و MgO دربرابر Fe2O3 نیز درستیِ یافته‌های به‌دست‌آمده از بخش سنگ‏‌نگاری و توزیع اندازة بلور را نشان می‌دهند و نشانة حضور دو یا چند ماگما در آشیانة ماگمایی و رویداد فرایند آمیختگی ماگمایی در هنگام پیدایش ماگما هستند.

بررسی توزیع اندازة بلور‏‌های پلاژیوکلاز در نمونه‏‌‏‌ها و وجود شکستگی و خمیدگی در نمودار‏‌های نیمه‏‌لگاریتمی نیز نشان‏‌دهندة ورود ماگمای جدید به آشیانة ماگمایی، برهم‌خوردن تعادل دمایی آشیانة ماگمایی، تبلور بلور‏‌های جدید با اندا‏زه‏‏‌های متفاوت و یا نشانة دیگری از رویداد فرایند آمیختگی ماگمایی در نمونه‏‌هاست.

زمان توقف ماگما و تبلور بلور‏‌های پلاژیوکلاز بیشتر، نزدیک به 1000 سال بوده است که نشان‏‌دهندة ژرفای کم جایگیری این توده در محدودة اکتشافی کودگان بوده است.

 

سپاس‌گزاری

این پژوهش با پشتیبانی مادی و معنوی شرکت طرح و توسعة لاله‌گستر نیایش به انجام رسیده است. ازاین‌رو، از مدیریت گرامی شرکت طرح و توسعة لاله‏‌گستر نیایش برای فراهم‌آوردن تسهیلات و امکانات لازم در راستای انجام این پژوهش و همچنین، در اختیارگذاشتن اطلاعات سپاس‌گزاری می‏‌شود.

 

 

Aghanabati, A. (2004) Geology of Iran. Geological Survey of Iran, Tehran, Iran (in Persian).
Arjmandzadeh, R., Karimpour, M. H., Mazaheri, S. A., Santos, J. F., Medina, J. M. and Homam, S. M. (2011) Two-sided asymmetric subduction; implications for tectonomagmatic and metallogenic evolution of the Lut block, eastern Iran. Journal of Economic Geology 1(3): 1-14.
Atherton, M. P. and Ghani, A. A. (2002) Slab breakoff: a model for Caledonian, Late Granite syn-collisional magmatism in the morphotectonic (metamorphic) zone of Scotland and Donegal, Ireland. Lithos 62: 65-85.
Berberian, M. and King, G. (1981) Towards a Paleogeography and Tectonic Evolution of Iran. Canadian Journal of Earth Sciences 18: 210-265.
Bolourian, G. (2010) Descriptions of the geological map of 1: 100000 Koodakan. Geological Survey and Mineral Exploration of Iran, Tehran, Iran.
Boorman, S., Boudreaux, A. and Kruger, F. J. (2004) The Lower Zone-Critical Zone Transition of the Bushveld Complex: A Quantitative Textural Study. Journal of Petrology 45(6): 1209-1235.
Cashman, K. V. (1993) Relationship between plagioclase crystallization and cooling rate in basaltic melts. Contributions to Mineralogy and Petrology 113: 126-142.
Cox, K. G., Bell, J. D. and Pankhurst, R. J. (1979) The Interpretation of Igneous Rocks. George Allen and Unwin, London, UK.
Eftekharnezhad, J. (1981) Tectonic division of Iran with respect to sedimentary basins. Journal of Iranian Petroleum Society 82: 19-28 (in Persian).
Gagnevin, D., Waight T. E., Daly, J. S., Poli, G. and Conticelli, S. (2007) Insights into magmatic evolution and recharge history in Capraia Volcano (Italy) from chemical and isotopic zoning in plagioclase phenocrysts. Journal of Volcanology and Geothermal Research 168(1-4): 28-54.
Girardi, V. A. V., Correa da Costa, P. C. and Teixeira, W. (2012) Petrology and SrNd characteristics of the Nova Lacerda dike swarm, SW Amazonian Craton: new insights regarding its subcontinental mantle source and Mesoproterozoic geodynamics. International Geology Review 54: 165–182.
Hessami, K., Nilforoushan, F. and Talbot, C. J. (2006) Active deformation within the Zagros Mountains deduced from GPS measurements. Journal of Geological Society 163: 143–148.
Higgins, M. D. (1996) Crystal size distributions and other quantitative textural measurements in lavas and tuff from Egmont volcano (Mt. Taranaki), New Zealand. Bulletin of Volcanology 58(2-3): 194-204.
Higgins, M. D. (2000) Measurement of crystal size distributions. American Mineralogist 85(9): 1105-1116.
Higgins, M. D. (2002a) Closure in crystal size distributions (CSD), verification of CSD calculations, and the significance of CSD fans. American Mineralogist 87(1): 171-175.
Higgins, M. D. (2002b) A crystal size-distribution study of the Ki-glapait layered mafic intrusion, Labrador, Canada: evidence for textural coarsening. Contributions to Mineralogy and Petrology 144: 314-330.
Higgins, M. D. (2006a) Quantitative textural measurements in igneous and metamorphic petrology. Cambridge University Press, Cambridge, UK.
Higgins, M. D. (2006b). Use of appropriate diagrams to determine if crystal size distributions (CSD) are dominantly semi-logarithmic, lognormal or fractal (scale invariant). Journal of Volcanology and Geothermal Research 154: 8-16.
Higgins, M. D. and Roberge, J. (2003) Crystal size distribution (CSD) of plagioclase and amphibole from Soufriere Hills volcano, Monteserrat: evidence for dynamic crystallization/textural coarsening cycles. Journal of Petrology 44: 1401-1411.
Hollocher, K., Robinson, P., Walsh, E. and Roberts, D. (2012) Geochemistry of amphibolite facies volcanics of the Storen Nappe in extensi-ons southwest and west from the Trondheim region, western Gneiss Region, Norway: A key to correlations and paleotectonic settings. American Journal of Science 312: 357-416.
Holton, T., Jamtveit, B. and Meakin, P. (2000) Noise and oscillatory zoning of minerals. Geochimica Acta 64: 1893-1904.
Irvine, T. N. and Baragar, W. R. A. (1971) A guide to the chemical classification of the common volcanic rocks. Canadian J Earth Science 8: 523-548.
Jerram, D. A. and Higgins, M. D. (2007) 3D analysis of rock textures: quantifying igneous microstructures. Elements 3(4): 239-245.
Jerram, D. A. and Cheadle, M. J. (2000) On the cluster analysis of grains and crystals in rocks. American Mineralogist 85(1): 47-67.
Jerram, D. A., Cheadle, M. C. and Philpotts, A. R. (2003) Quantifying the building blocks of igneous rocks: are clustered crystal frameworks the foundation? Journal of Petrology 44: 2033-2051.
Jung, S., Pfander, J. A., Brugmann, G. and Stracke, A. (2005) Sources of primitive alkaline volcanic rocks from the Central European Volcanic Province (Rhon, Germany) inferred from Hf, Os and Pb isotopes. Contributions to Mineralogy and Petrology 150: 546-559.
Karimpour, M. H., Stern, C. R., Farmer, L., Saadat, S. and Malekzadeh, A. (2011) Review of age Rb-Sr geochemistry and petrogenesis of Jurassic to Quaternary igneous rocks in Lut Block, Eastern Iran. Geopersia 1(1): 19–54.
Kolb, M., Von Quadt, A., Peytcheva, I., Heinrich, C. A., Fowler, S. J. and Cvetković, V. (2013) Adakite-like and normal arc magmas: distinct fractionation paths in the East Serbian Segment of the Balkan–Carpathian Arc. Journal of Petrology 54: 421-451.
Kumar, S. and Rino, V. (2006) Mineralogy and geochemistry of microgranular enclaves in Palaeoproterozoic Malanjkhand granitoids, central India: Evidence of magma mixing, mingling, and chemical equilibration. Contributions to Mineralogy and Petrology 152(5): 591-609.
Kuşcu, G. G. and Floyd, P. A. (2001) Mineral compositional and textural evidence for magma mingling in the Saraykent volcanics. Lithos 56(2-3): 207-230.
Marsh, B. D. (1988) Crystal size distribution (CSD) in rocks and the kinetics and dynamics of crystallization theory. Contributions to Mineralogy and Petrology 99: 277–291.
Marsh, B. D. (1998) On the interpretation of crystal size distributions in magmatic systems. Journal of Petrology 39(4): 553-599.
Middlemost, E. A. K. (1994) Naming materials in magma/igneous rock system. Earth Science- Reveiws 37: 215–224.
Nagudi, N., Koberl, C. and Kurat, G. (2003) Petrography and Geochemistry of the sigo granite, Uganda and implications for origin. Journal of African earth Sciences 36(1): 1-14.
Nowroozi, A. A. (1976) Seismotectonic provinces of Iran. Bulletin of the Seismological Society of America 66: 1249–1276.
Omaraeei, S. and Niroumand, S. (2014) Geochemistry, mineralogy, alteration and fluid inclusion studies in copper gold vein system Koodakan prospecting area, South Khorasan province. Journal of Advanced Applied Geology 19(6): 34–47 (in Persian with English abstract).
Pang, H., Chen, J., Pang, X., Liu, L., Liu, K. and Xiang, C. (2013) Key factors controlling hydrocarbon accumulations in Ordovician carbonate reservoirs in the Tazhong area, Tarim basin, western China. Marine and Petroleum Geology 43: 88-101.
Pearce, J. A., Harris, N. B. W. and Tindle, A. G. (1984) Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. Journal of Petrology 25: 956-983.
Randolph, A. D. and Larson, M. A. (1971) Theory of Particulate Processes. New York, Academic Press, US.
Renzulli, A. and Santi, P. (1997) Sub-volcanic crystallization at Stromboli (Aeolian Islands, southern Italy) preceding the Sciara del Fuoco sector collapse: evidence from monzonite lithic suite. Bulletin of Volcanology 59(1): 10-20.
Rudnick, R. and Gao, S. (2003) Composition of the Continental Crust. Treatise Geochem 3. Elsevier 3: 1-64.
Ruprecht, P. and Wörner, G. (2014) Variable regimes in magma systems documented in plagioclase zoning patterns: El Misti stratovolcano and Andahua monogenetic cones. Journal of Volcanology and Geothermal Research 165(3-4): 142-162.
Saccani, E., Delavari, M., Beccaluva, L. and Amini, S. A. (2010) Petrological and geochemical constraints on the origin of the Nehbandan ophiolitic complex (eastern Iran): implication for the evolution of the Sistan Ocean. Lithos 117: 209-228
Stamatelopoulou-Seymour, K., Vlassopoulos, D., Pearce, T. H. and Rice, C. (1990) The record of magma chamber processes in plagioclase phenocrysts at Thera Volcano, Aegean Volcanic Arc, Greece. Contributions to Mineralogy and Petrology 104: 73–84.
Sun, S. S. and McDonough, W. F. (1989) Chemical and isotopic systematics of oceanic basalts; implications for mantle composition and processes. In: Magmatism in the ocean basins (Eds. Saunders, A. D. and Norry, M. J.) 42: 313-345. Geological Society of London, London, UK.
Swain, C. J. and Kirby, J. F. (2008) An accuracy assessment of the fan wavelet coherence method for elastic thickness estimation. Geochemistry, Geophysics, Geosystems 9(3): Q03022.
Treuil, M. and Joron, J. L. (1975) Utilisation des elements hygromag- matophiles pour la simplification de la modelisation quantita- tive des processus magmatiques. Exemples de l’Afar et de la Dorsale Medioatlantique Rend. Society Italy Mineralogy and Petrology 31: 125-174.
Tsuchiyama, A. 1985. Dissolution kinetics of plagioclase in the melt system diopside-albite-anorthite, and origin of dusty plagioclase in an- desites. Contributions to Mineralogy and Petrology 89: l-16.
Voorhees, P. W. (1992) Ostwald ripening of tow-phase mixture. Annual Review Materials Science 22: 197-215.
Wilson, M. (1989) Igneous Petrogenesis: A Global Tectonic Approach. Unwin Hyman, London, UK.
Yousefzadeh, M., Rahmani, A. and Mohammad, S. S. (2019) Petrology and tectonomagmatic setting of volcanic and subvolcanic rocks in the east of Khousf (Southwest of Birjand). Iranian Journal of Petrology 37: 1-22 (in Persian).
Zhou, M. F., Robinson, P. T. and Bai, W. J. (1994) Formation of podiform chromitites by melt/rock interaction in the upper mantle. Mintralium Deposita 29: 98-101.