Mineral chemistry and geochemistry of Kabbarik syenite stock, North of Zanjan

Document Type : Original Article

Authors

1 Assistant Professor, Department of Geology, Faculty of Sciences, University of Mohaghegh Ardabili, Ardabil, Iran,

2 Professor, Earth Science Department, Natural Sciences Faculty, University of Tabriz, Tabriz, Iran,

3 M.Sc. in Exploration Mining Engineering, Arak University, Arak, Iran,

Abstract

The study area, Tarom, is situated in northwestern Iran, within the Zanjan province. Geologically, it falls within the western Alborz zone, a volcanic-plutonic belt, oriented northwest-southeast. This region, extending from Qazvin (west of Taleghan) to the north-northwest of Miane and north of the Manjil embayment, is bounded by the Abhar-Zanjan-Miane axis to the south (Ghorbani, 2009). Within this mountain range, numerous intrusive masses, primarily granodiorite, have intruded into the Eocene volcanic rocks. These masses, including the Bakhtar Takhestan, Khorasanlu, Zaker, and Chal masses, follow a northwest-southeast trend, aligning with the overall orientation of the Tarom mountains. They are arranged in two parallel rows, cutting through the older volcanic formations.
Geology of the area
The most important rock units in the Tarom Mountains consist of Eocene pyroclastic rocks, including andesite, rhyolite, basalt, andesite basalt, and rhyodacite, as well as shear tuff and andesite tuff. The region's main elevations are composed of Eocene volcanic and volcano-clastic rocks, along with Oligocene intrusive masses, which form the bulk of the mountain unit. Oligocene trachytic deposits, which are deposited on the Eocene units at an angle of approximately 23 degrees, are exposed in the northern part of the Kabbarik syenitic porphyroid stock. The syenitic porphyroid stock has intruded into the Eocene volcanic and volcano-clastic rocks.
This study is grounded in field surveys, sampling of the slightly intrusive mass, and the examination of thin sections. Additionally, chemical analyses were conducted on ten samples from Zarazma Company to determine major, minor, and trace elements using inductively coupled plasma mass spectrometry (ICP-MS). Two samples were also analyzed using electron microscopy. To analyze the mineral composition of Kabbarik rocks for major elements, a JEOL JXA-8100 Superprobe electron microprobe (EMP) equipped with an Oxford Instruments INCA EDS system was employed at the Institute of Geology and Geophysics, Chinese Academy of Science. A 5-micron spot size, a 20 nA beam current, and a 15 keV accelerating voltage were utilized during these analyses.
Petrography
The intrusive Kabbarik stock exhibits a porphyritic texture, characterized by euhedral to subhedral phenocrysts of alkali feldspar. Its primary mineral composition includes alkali feldspar (Orthoclase) constituting 65-70 Vol.%, followed by 10-15 Vol.% clinopyroxene (diopside), 5-7 Vol.% anorthoclase, and 3-5 Vol.% secondary minerals such as olivine, altered to iddingsite. Minor amounts of apatite, zircon, and dark minerals are also present. Orthoclase phenocrysts often display Carlsbad and polysynthetic twinning. Megacrysts of Orthoclase may contain tiny inclusions of olivine and clinopyroxene, exhibiting a poikilitic texture. Some crystals exhibit a sieve texture, and in certain cases, overgrowth on plagioclase results in a Rapakivi texture. Clinopyroxene (diopside) crystals are generally euhedral to subhedral, with some containing small dark mineral inclusions distributed along their edges and conforming to their shape. Subhedral crystals of iddingsitized olivine are scattered throughout the rock matrix in small quantities. The overall texture of the microlithic rock remains porphyritic.
Mineral Chemistry
The chemical composition of clinopyroxene and alkali feldspar minerals from the syenitic stock was analyzed using a JEOL JXA-8200 electron microscope at the Department of Geology and Geophysics, Chinese Academy of Sciences.
Based on relevant diagrams, the clinopyroxenes were classified as diopside-type (Eby et al., 1998; Morimoto, 1988). The chemical composition of the pyroxenes indicates a sub-alkaline, volcanic arc magmatic setting. EPMA analysis of the alkali feldspars revealed an anorthoclase composition (Deer et al., 1991).
Whole Rock Geochemistry
The chemical analysis of significant, rare, and rare earth elements (ICP-MS) reveals a calc-alkaline composition for the intrusive mass of the syenite zone and its parent magma. Normalized spider diagrams compared to chondrite demonstrate a clear enrichment of light rare earth elements (LREE) and incompatible elements relative to heavy rare earth elements (HREE) in the examined rocks. Based on tectonic environment diagrams, the studied samples are classified as post-collisional.
Results
The predominant rock formations in the Tarom Mountains consist of Eocene volcanic and volcaniclastic rocks, which constitute the bulk of the mountain range. The investigated syenitic stock, intruding into the overlying rocks, dates back to the Oligocene epoch. Based on spider diagrams, the absence of significant heavy rare earth element (HREE) depletion suggests a mantle that has been metasomatized by subduction-related fluids. Studies indicate that the analyzed samples originated in a setting associated with an arc and a post-collisional environment.
 
 

Keywords

Main Subjects


منطقة کبریک در شمال‏‌باختری ایران و شمال استان زنجان جای دارد. طارم از دیدگاه پهنه‏‌بندی واحدهای زمین‏‌شناسی ایران (Nabavi, 1977) در پهنة البرز باختری جای دارد و بر پایة (Aghanabati, 2005) بخشی از پهنة ایران مرکزی شمرده می‌شود (شکل‌های 1-A و 1-B). این نوار ولکانو- پلوتونیک روند شمال‏‌باختری-جنوب‏‌خاوری دارد. پهنة طارم- هشجین از قزوین (باختر طالقان) آغاز می‏‌شود و تا شمال- شما‏‌ل‏‌باختری میانه و شمال نواحی فروافتادگی منجیل ادامه دارد و از جنوب به محور ابهر- زنجان- میانه می‌رسد. در حقیقت از دید ساختاری، این کمربند از شمال، جنوب و باختر به گسل‏‌های بزرگ سفیدرود (در شمال) و ادامه گسل تبریز- سلطانیه و سلطانیه- تاکستان (در جنوب) می‏‌رسد و از باختر با گسل آستارا- مریوان محدود می‏‌شود (Ghorbani, 2009).

 

 

شکل 1. A) واحدهای اصلی رسوبی و ساختاری ایران (Aghanabati, 2005) و محدوده بررسی‌شدة کبریک که با چهار گوش نشان داده شده است؛ B) نقشة زمین‏‌شناسی محدودة کبریک (مقیاس: 1:25000).

Figure 1. A) The major sedimentary and structural units of Iran and the location of the study area of Kabbarik shown by a quadrangle (Aghanabati, 2005); B) Geological map of the Kabbarik area (Scale: 1:25000).

 

 

بررسی‏‌های پیشین نشان می‏‌دهند سرزمین ایران در بازة زمانی پالئوزوییک فازهای مختلفی از رخدادهای زمین‌ساختی کششی مرتبط با زایش پهنه‏‌های اقیانوسی تتیسی را پشت سر گذاشته است (Berberian and King, 1981; Derakhshi and Ghasemi, 2014; Ghasemi and Dayhimi, 2015; Azizi et al., 2017; Delavari et al., 2019). زمین‏‌شناسی و زمین‏‌ساخت البرز باختری-آذربایجان در برابر دیگر مناطق ایران پیچیدگی‏‌های خاصی دارد و رویدادهای مهمی از مزوزوییک تا سنوزوییک را در خود ثبت کرده است. با فعالیت فاز کوهزایی آلپی از ائوسن تا کواترنری، ماگماتیسم گسترده‏‌ای در پهنة البرز باختری-آذربایجان رخ داده است (Verdel et al., 2011; Allen et al., 2013). در البرز باختری رخداد فاز کششی ائوسن تا الیگوسن آغاز فعالیت‏‌های ماگمایی گسترده‏‌ای را به‌دنبال داشته است. در پی رخداد فاز کوهزایی پیرنه (الیگوسن زیرین) سری‏‌های ائوسن در البرز باختری چین‏‌خورده و گسل‏‌های کششی به موازات محور چین‏‌خوردگی‏‌ها پدید آمده‏‌اند. در آذربایجان‏‌شرقی و طارم، به‌ویژه در شمال زنجان چین‏‌هایی با روند خاوری-باختری پدید آمده‏‌اند (Moayyed, 2001). پهنة ماگمایی البرز- آذربایجان که در بخش باختری پهنة ماگمایی البرز جای دارد، تکاپوی ماگمایی گسترده‏‌ای، به‌ویژه در سنوزوییک و با سرشت کالک‏‌آلکالن داشته است (Moayyed, 2001). سنگ‌های آتشفشانی منطقه زنجان و کوه‏‌های طارم بیشتر به سازند کرج نسبت داده شده‏‌اند (Ghorbani et al., 2013; Nabatian et al., 2014). موید (Moayyed, 2001) با بررسی دقیق سنگ‌های ماگمایی پهنة ماگمایی ولکانو-پلوتونیک ترشیری البرز باختری- آذربایجان، سری سنگ‌های آذرین آتشفشانی و درونی ائوسن- الیگوسن منطقه هشتجین را کالک‏‌آلکالن پتاسیم بالا و شوشونیتی دانسته و محیط زمین‏‌ساختی آن را کمان ماگمایی پسابرخوردی دانسته است.

در کوه‌های طارم، واحدهای پرکامبرین به‌صورت سنگ‌های دگرگونی و گرانیت در چند نقطه بیرون‌زدگی دارند. پالئوزوییک زیرین (کامبرین و اردوویسین)‏‌ و دونین بالایی و کربونیفر در کوه‌های طالش چندین رخنمون دارند. ژوراسیک شامل نهشته‏‌های سازند شمشک و آهک لار و کرتاسه شامل نهشته‏‌های آهک اربیتولین‏‌دار است (نقشة زمین‏‌شناسی 1:250000 زنجان (Stöcklin and Eftekharnezhad, 1969) و نقشة زمین‏‌شناسی 1:100000 ماسوله (Device and Ghoreishi, 1969). در این رشته کوه‏‌ها، توده‏‌های آذرین درونی بسیاری با روند شمال‏‌باختری- جنوب‏‌خاوری به پیروی از روند کوه‏‌های طارم یافت می‌شوند. این توده‏‌ها (تودة باختر تاکستان، تودة خراسانلو، تودة ذاکر و تودة چال) بیشتر گرانودیوریتی هستند و در امتداد هم و در دو ردیف موازی، سنگ‌های آتشفشانی ائوسن را قطع کرده‏‌اند. در پهنة طارم بیشتر بررسی‏‌های انجام‌شده بر پایة کانه‏‌زایی‏‌های موجود بوده‌اند و در آنها کمتر به بررسی سنگ‏‌شناختی و زمین‏‌شیمیایی ماگماتیسم منطقه پرداخته شده است (Ebrahimi et al., 2016). بر پایة پژوهش‏‌های (Kheirkhah and Karimi, 2021)، بررسی‌های سنگ‏‌شناختی، گدازه‏‌های حد واسط الیگوسن منطقه اوچ‏‌بلاغ در شمال زنجان و جنوب‌باختری طالش در بلوک البرز باختری-آذربایجان و بخشی از زیرپهنة طارم - هشجین، ترکیب آندزیت- تراکی‏‌آندزیتی دارند و گوشتة بالایی خاستگاه ماگمای سازندة گدازه‏‌های کالک‏‌آلکالن پتاسیم بالا تا شوشونیتی بوده است (Kheirkhah and Aghaali, 2018) که در یک جایگاه زمین‏‌ساختی کوهزایی مرتبط با فرورانش برون ریخته‏‌اند. توده‏‌های آذرین درونی منطقة ذاکر در خاور زنجان در بلوک البرز باختری-آذربایجان و بخشی از زیرپهنة طارم - هشجین، با ترکیب کوارتزمونزونیت و میکرودیوریت تا کوارتزدیوریت سنگ‌های ولکانیک و ولکانو-کلاستیک ائوسن را قطع کرده‌اند و سنی جوان‌تر از ائوسن از خود نشان می‏‌دهند (Rezanajad, 2013). توده‏‌های آذرین درونی میکرودیوریت تا کوارتزدیوریت دگرسانی پروپیلیتیک و تودة کوارتزمونزونیت پورفیری دگرسانی سریسیتیک از خود نشان می‏‌دهند (Nabatian et al., 2016). سنگ‌های سینیتی یک بخش فرعی از سلسله آذرین هستند که بیشتر آنها در شمال‏‌باختری پهنة سنندج-سیرجان در جاهایی مانند ارومیه (Jafari et al., 2018)، پیرانشهر (Mazhari et al., 2009)، سردشت (Fazlnia, 2019)، سقز (Sepahi and Athari, 2006)، آلموقولاق (Amiri et al., 2017; Jamshidi-Badr, et al., 2018)، گلپایگان (Davoudian et al., 2007) و دهگلان (Sarjoughian et al., 2016) جای دارند. گمان می‌رود بیشتر آنها در یک محیط زمین‌ساختی کششی پدید آمده‏‌اند و یا در ارتباط با گرانیتوییدهای نوع A رخ داده‏‌اند (Sepahi and Athari, 2006; Davoudian et al., 2007; Mazhari et al., 2009; Sarjoughian et al., 2016; Amiri et al., 2017). سنگ‌زایی بقایای سینیت همچنان موضوع بحث برانگیزی است و برخی الگو‏‌ها در تضاد مستقیم با نتایج بررسی‌های تجربی هستند. یکی از پیامدهای این الگو‏‌های زایشی این است که ماگماهای سینیتی را می‏‌توان در محیط‏‌های مختلف زمین‌ساختی در یک چرخة کوهزایی جای داد (Adetunji et al., 2018). با وجود این، به‌طور کلی فرض می‏‌شود که بیشتر ماگماتیسم سینیتی ویژگی محیط‏‌های پس از کوهزایی، درون‌صفحه‏‌ای و کششی است (Kaul and Cordiani, 2000; Martin and Devito, 2005; Yang et al., 2005)؛ اگرچه شواهدی نشان می‏‌دهند سینیت‏‌ها می‏‌توانند در محیطی همزمان با کوهزایی در هنگام فرورانش پدید آیند (Chen et al., 2003; Yan and Jiang, 2019). سینیت‏‌ها بر پایة کانی‏‌شناسی به دو دستة هیپرسولووس (بی پلاژیوکلاز و با یک آلکالی‏‌فلدسپار پرتیتی یا آنتی‏‌پرتیتی) و ساب‏‌سولووس (دارای پلاژیوکلاز و آلکالی‏‌فلدسپار) دسته‌بندی می‏‌شوند. سینیت‏‌ها به دو صورت دیده می‏‌شوند: توده‏‌های نیمه‌ژرف مستقل و کوچک و رخساره‌‏‌های آذرین محلی و کوچک مرتبط با توده‏‌های گرانیتوییدی بزرگ. بر پایة سن‏‌سنجی به روش U-Pb برای سینیت بردکیش (مجموعه پلوتونیک ارومیه) سن 8/0 ± 1/91 میلیون سال پیش (تورونین) برآورد شده است. تجزیة شیمیایی سنگ کل با پیدایش ماگمای سینیتی در یک سیستم کمان با تبلوربخشی یک ماگمای مادر مافیک با مقداری خرج پوسته‏‌ای سازگار است (Jafari et al., 2020). بررسی‌ نخستین گزارش علمیِ این استوک سینیتی در منطقة کبریک و اینکه در نقشة 1:100000 هشجین این استوک شناسایی نشده است؛ همچنین، برای نخستین‌بار در یک نوشتار پژوهشی معرفی می‏‌شود و پیشینه‏‌ای برای این بررسی وجود ندارد تا بتوان یافته‏‌های آن را با بررسی‌های پیشین مقایسه کرد؛ پس لازم است ترکیب سنگ‏‌شناسی بر پایه مودال کانیایی، ویژگی‏‌های زمین‏‌شیمیایی و کانی‏‌شیمی این استوک ارزیابی و شناخته شود و تلاش این پژوهش در این راستاست.

زمین‏‌شناسی منطقه

پی‏‌سنگ کهن پهنة ماگمایی البرز-آذربایجان از سنگ‌های دگرگونی ساخته شده است که در شمال‏‌باختری لاهیجان، جنوب‏‌باختری خوی، شمال تبریز و شمال-شمال‏‌باختری زنجان رخنمون دارند. سن این پی‏‌سنگ، نئوپروتروزوییک بالایی-کامبرین زیرین است (Shafaii Moghadam et al., 2015). به‌طور کلی، روند عمومی ساختارهای زمین‏‌شناسی در منطقة طارم، همسان با دیگر مناطق البرز باختری WNW-ESE است که جایگیری و پراکندگی سطحی توده‏‌های آذرین درونی منطقه نیز از روند گفته شده پیروی می‏‌کند. توده‏‌های آذرین درونی طارم با سن ائوسن بالایی شامل توده زنجان (کوارتزمونزودیوریت)، ذاکر (کوارتزمونزونیت)، مروارید (مونزوگرانیت تا سینوگرانیت) و کوه تبر (میکروکوارتزدیوریت پورفیری) هستند که به درون واحدهای آتشفشانی و آذرآواری سازند کرج نفوذ کرده‏‌اند (Nabatian, 2008). میزبان توده‏‌های آذرین درونی طارم، بیشتر فورآن‌های آتشفشانی زیردریایی (سازند کرج) شامل سنگ‌های آذرآواری و گدازه‏‌های تراکیتی، تراکی‏‌آندزیتی، آندزیتی و الیوین‏‌بازالتی است (Nabatian, 2008). به‌طور کلی، روند عمومی ساختارهای زمین‏‌شناسی در منطقة طارم، همسان با دیگر مناطق البرز باختری روندِ WNW-ESE دارد. جایگیری و پراکندگی سطحی توده‏‌های آذرین درونی منطقةکبریک نیز از روند یادشده پیروی می‏‌کند. مهم‏‌ترین واحدهای سنگی در کوه‏‌های طارم متشکل از واحدهای آذرآواری ائوسن (آندزیت، ریولیت، بازالت، آندزیت بازالت و ریوداسیت همراه با توف برشی و توف آندزیتی) هستند. در کوه‏‌های طارم نیز مانند ایران مرکزی و البرز مرکزی حرکت‌های پالئوسن با بیرون‌ریختن سنگ‌های آتشفشانی زیردریایی پالئوسن ادامه یافته است (Moeinvaziri, 1997). رخنمون‌های منطقه که ارتفاعات اصلی محدوده را ساخته‌اند، متشکل از سنگ‌های آتشفشانی و ولکانوکلاستیک ائوسن و توده‌های آذرین درونی مربوط به الیگوسن (به‌علت نبود انجام سن‏‌سنجی مطلق روی استوک یادشده، از آن‌جایی ‏‌که این پهنه ادامه پهنة ولکانو-پلوتونیک طارم به‌سوی باختر و شمال‏‌باختری است و با توجه به شباهت‌های ترکیبی توده‏‌های گرانیتوییدی پهنة طارم و ادامه این پهنه به‌سوی شمال‏‌باختری ایران (پهنة ارسباران)، این استوک هم‏‌سنِ توده‏‌های این دو منطقه دانسته شده است) هستند که بخش بزرگی از واحد کوهستانی را ساخته‌اند (شکل 2-A).

 

 

شکل 2. A) رخنمون واحدهای آتشفشانی ائوسن در منطقة کبریک؛ B) فرسایش‌پذیری و نفوذپذیری شدید واحدهای توفی نسبت به گدازه‏‌های بازالتی و آندزیتی و رخداد فرسایش تفریقی؛ C) دم‌ تراکیتی مربوط به الیگوسن که به‌صورت دگرشیبی زاویه‏‌دار روی واحدهای توف جای گرفته است؛ D، E) نمایی از استوک پورفیرویید سینیتی کبریک (دید رو به شمال‏‌باختری).

Figure 2. A) Outcrop of Eocene volcanic rocks in Kabbarik area; B) Erodibility and extreme permeability of tuff units compared to basaltic and andesitic lavas and creating differential erosion; C) Oligocene trachytic tail, which is placed on the tuff units as an angular unconformity; D, E) A view of Kabbarik syenitic porphyroid stock (NW view)

 

 

فرسایش‏‌پذیری و نفوذپذیری شدید واحدهای توفی متعلق به سازند کرج به میزان چشمگیری در رخسارة زمین‌ریخت‌شناسی این بخش تأثیر گذاشته است و در نزدیکی آنها گدازه‌های آندزیتی و بازالتی مربوط به ائوسن (بر پایة اطلاعات نقشة 1:100000 هشجین و وجود فسیل‏‌های نومولیتس و دیسکوسیکلینا در لابلای توف‏‌ها و گدازه‏‌های یادشده (Moayyed, 2001)) جای می‌گیرند و بیشتر آنها صخره‌ساز هستند و در رخداد فرسایش تفریقی[1] نقش عمده‌ای دارند (شکل 2-B). در منطقة کبریک به فاصلة کمی از استوک سینیتی، رخنمون‏‌های بزرگی از توده‏‌های کوارتزمونزونیتی تا گرانودیوریتی گلجین-رشت‏‌آباد به سن الیگوسن دیده می‌شوند که درون سنگ‌های آتشفشانی و ولکانو-کلاستیک ائوسن تزریق شده‏‌اند و این ویژگی احتمال جدایش ماگمایی این سینیت‏‌ها از ماگمای مادر این توده‏‌ها را افزایش می‏‌دهد. روانه‌های تراکیتی مربوط به الیگوسن که با دگرشیبی زاویه‏‌دار نزدیک به 23 درجه روی واحدهای ائوسن جای گرفته‏‌اند (Moayyed, 2001) (شکل 2-C)، در بخش شمال استوک پورفیرویید سینیتی کبریک رخنمون دارند. در شکل‏‌های 2-D و 2-E جایگاه استوک پورفیرویید سینیتی نشان داده شده است.

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

این پژوهش بر پایه بررسی‏‌های میدانی، نمونه‏‌برداری از استوک کبریک، مطالعة مقاطع نازک و داده‏‌های به‌دست‌آمده از تجزیة شیمیایی (10 نمونه سنگ) در شرکت زرآزما برای تجزیه عنصرهای اصلی (XRF)، فرعی و کمیاب (ICP-MS) و 2 نمونه سنگی ریزکاو الکترونی انجام شده است.

تجزیه عنصرهای اصلی کانی‏‌های منتخب سازندة سنگ‌های در منطقة کبریک در مؤسسه زمین‏‌شناسی و ژئوفیزیک آکادمی علوم چین با استفاده از ریزکاو الکترونی JEOL JXA-8100 Superprobe (EMP) با سیستم Oxford Instruments INCA EDS با باریکة پرتوی به قطر 5 میکرون تجزیه شدند. جریان پرتو 20 نانوآمپر و ولتاژ شتاب‌دهنده keV15 لحاظ شده است.

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

سنجش عنصرهای اصلی به روش فلورسانس پرتوی ایکس (XRF) روی قرص‏‌های تهیه‌شده از پودر سنگ انجام شد. آماده‌سازی نمونه در روش XRF به‌صورت تهیة نمونه‌های پودری انجام می‌شود که شامل خشک‌کردن، خردایش، نمونه‌گیری و نرمایش است و در پایان، قرص فشرده تهیه می‌شود. فرایند فشرده‌سازی پودر سازندة قرص‌ها شامل آسیاب‌کردن نمونه برای به‌دست‌آوردن یک پودر خوب با حالت ایده‌آل، اندازة ذرات کمتر از 75 میکرومتر و سپس مخلوط‌کردن آنها با چسب مناسب و اعمال فشار نزدیک به 20 تا 30 نیوتن بر مخلوط همگن درون قالب به‌مدت زمان 15 ثانیه برای ساخت قرص نمونه‌ها در قالب یا رینگ‌های مخصوص است. قرص هر نمونه با دستگاه XRF در شرکت زرآزما تجزیه شد و مقدار اکسیدهای عنصرهای اصلی آن اندازه‏‌گیری شد.

مقدار عنصرهای فرعی و کمیاب نیز با دستگاه ICP-MS در شرکت زرآزما به‏‌دست آورده شد. در روش ICP-MS، هضم هر نمونه هم به‌صورت ذوب با کمک لیتیم‏‌متابورات و سپس انحلال در اسیدنیتریک رقیق بود وهم از تجزیة چند اسیدی و به‏‌کارگیری هضم میکروویو بهره گرفته شد. در این روش آستانه آشکارسازی برای همة اکسیدها نزدیک به 05/0 درصدوزنی بود.

سنگ‏‌نگاری

استوک سینیتی کبریک ساخت پورفیری دارد و با فنوکریستال‏‌های شکل‏‌دار تا نیمه‌شکل‏‌دار آلکالی‏‌فلدسپار شناخته می‏‌شود. رنگ آن در نمونة دستی خاکستری مایل به سبز است و یک تودة آذرین درونی نیمه‌آتشفشانی مزوکرات شمرده می‏‌شود. بافت کلی سنگ میکرولیتی پورفیریک است (شکل 3-A). کانی‏‌های اصلی سازندة آن شامل بلورهای شکل‏‌دار تا نیمه‌شکل‏‌دار آلکالی‏‌فلدسپار (ارتوکلاز) به‌صورت فنوکریستال تا مگاکریستال (65 تا 70 درصدحجمی، مجموع میکروکریستال‏‌های ارتوکلاز در خمیره و مگاکریستال‏‌ها)، کلینوپیروکسن (دیوپسید) شکل‏‌دار تا نیمه‌شکل‏‌دار (10 تا 15 درصدحجمی)، آنورتوکلاز (5 تا 7 درصدحجمی)، کانی‏‌های فرعی الیوینِ دگرسان‌شده به ایدنگسیت (3 تا 5 درصدحجمی) به‌همراه مقادیر اندکی آپاتیت، زیرکن (شکل‌های 3-E و 3-F) و کانی کدر است. فنوکریستال‏‌های ارتوکلاز ماکل کارلسباد و پلی‏‌سینتتیک نشان می‌دهند و مگاکریستال‏‌های آن میانبار‏‌های ریزی از الیوین و کلینوپیروکسن دارند که بافت پویی‏‌کیلیتیک را به نمایش گذاشته‏‌اند (شکل 3-B). برخی بلورهای ارتوکلاز بافت غربالی نیز دارند. در برخی بخش‌ها، تجمعات گلومروپورفیریک از ارتوکلاز دیده می‏‌شود.

 

 

شکل 3. تصویرهای میکروسکوپی (در XPL) از A) بافت میکرولیتی پورفیری و فنوکریستال‏‌های ارتوکلاز و کلینوپیروکسن در خمیرة میکرولیتیک؛ B) مگاکریستال ارتوکلاز با میانبارهای ریز کلینوپیروکسن؛ C) فنوکریستال شکل‏‌دار کلینوپیروکسن با خلیج‏‌خوردگی و میانبار‏‌های ریز کانی تیره که از شکل بلور پیروی کرده‏‌اند؛ D) فنوکریستال الیوین که اندکی به ایدنگسیت دگرسان شده است؛ E) بلور آپاتیت در کنار فنوکریستال الیوینی که اندکی به ایدنگسیت دگرسان شده است؛ F) بلور زیرکن در حاشیة فنوکریستال ارتوکلاز (نام اختصاری کانی‌ها از ویتنی و اوانس (Whitney and Evans, 2010) برگرفته شده است).

Figure 3. Photomicrographs (in XPL) of A) porphyric microlithic texture and orthoclase and clinopyroxene phenocrysts in the microlithic groundmass; B) Orthoclase megacrysts with small inclusions of clinopyroxene; C) Euhedral clinopyroxene phenocrysts with corrosion gulf and small dark mineral inclusions following the crystal shape; D) Olivine phenocryst slightly altered to iddingsite; E) Apatite crystal next to olivine phenocryst that is slightly altered to iddingsite; F) Zircon crystal at the edge of an orthoclase phenocryst (Mineral abbreviations are from Whitney and Evans, 2010).

 

بلورهای کلینوپیروکسن (دیوپسید) شکل‏‌دار تا نیمه‌شکل‏‌دار هستند. برخی بلورهای کلینوپیروکسن خلیج‏‌خوردگی دارند و میانبار‏‌های ریزی از کانی‏‌های کدر دارند که موازی حاشیه‌های بلور و به‌پیروی از شکل آن پراکنده شده‏‌اند (شکل 3-C). در برخی بخش‌ها تجمعات گلومروپورفیریک از این کانی دیده می‏‌شود. بلورهای نیمه‌شکل‏‌دار الیوین ایدنگسیتی‌شده با فراوانی اندک در زمینة سنگ پراکنده‏‌اند (شکل 3-D). سینیت کبربک با داشتن ارتوکلاز و آنورتوکلاز یک آلکالی‏‌فلدسپار سینیت هیپرسولووس شمرده می‏‌شود.

شیمی کانی‏‌ها

فرمول ساختاری کانی کلینوپیروکسن بر پایة 6 اتم اکسیژن و 4 کاتیون به‌دست آورده شده است. داده‌های به‌دست‌آمده از تجزیه این کانی در جدول 1 آورده شده‌اند.

 

جدول 1. داده‌های تجزیة نقطه‏‌ای کلینوپیروکسن‏‌ (بر پایة درصدوزنی) در سینیت کبریک به‌همراه فرمول ساختاری به‌دست‌آمده (برپایة apfu) برپایة 6 اتم اکسیژن، توزیع کاتیون‌ها در جایگاه‌های بلوری و سازنده‌های پایانی.

Table 1. The EPMA analytical data of clinopyroxene (in wt%) in Kabbarik syenite, the calculated structural formula (in apfu) based on 6 oxygens, cation distribution in crystal sites, and the end members.

KBR9 core

KBR8 mid

KBR7 core

KBR6 mid

KBR5 rim

KBR4 rim

KBR3 rim

KBR2 core

KBR1 core

Point No.

50.54

50.4

50.43

50.19

50.18

50.14

50.17

50.77

50.65

SiO2

0.70

0.68

0.76

0.53

0.49

0.51

0.51

0.60

0.55

TiO2

1.46

1.44

1.58

1.29

1.29

1.26

1.28

1.47

1.37

Al2O3

14.47

14.44

14.31

14.54

14.75

14.82

14.68

13.51

14.28

FeO

0.77

0.79

0.75

0.83

0.85

0.79

0.81

0.71

0.77

MnO

10.03

9.81

9.91

9.82

9.58

9.65

9.72

10.35

9.99

MgO

20.91

20.92

20.98

20.98

21.03

21.04

21.05

21.09

21.09

CaO

0.51

0.49

0.54

0.63

0.61

0.63

0.59

0.54

0.60

Na2O

0.00

0.00

0.00

0.00

0.00

0.01

0.01

0.00

0.00

K2O

0.00

0.02

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Cr2O3

99.39

98.98

99.25

98.84

98.78

98.86

98.82

99.04

99.29

Sum

1.97

1.97

1.97

1.97

1.97

1.97

1.97

1.98

1.98

Si

0.06

0.02

0.02

0.02

0.01

0.01

0.01

0.02

0.02

Ti

0.07

0.07

0.07

0.06

0.06

0.06

0.06

0.07

0.06

Al

0.47

0.47

0.47

0.47

0.48

0.49

0.48

0.44

0.47

Fe

0.03

0.03

0.02

0.03

0.04

0.02

0.04

0.02

0.02

Mn

0.5

0.57

0.58

0.58

0.56

0.57

0.57

0.60

0.58

Mg

0.89

0.88

0.88

0.88

0.89

0.89

0.89

0.88

0.88

Ca

0.06

0.04

0.04

0.05

0.05

0.05

0.04

0.04

0.05

Na

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

K

0.00

0.001

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Cr

4.05

4.05

4.05

4.06

4.060

4.06

4.06

4.05

4.06

Sum

1.95

1.95

1.94

1.94

1.95

1.94

1.94

1.95

1.95

Si

0.05

0.05

0.06

0.06

0.05

0.06

0.06

0.05

0.05

Al(IV)

2.00

2.00

2.00

2.00

2.00

2.00

2.00

2.00

2.000

Sum

0.01

0.02

0.02

0.00

0.01

0.00

0.00

0.02

0.01

Al(VI)

0.41

0.42

0.40

0.38

0.40

0.39

0.39

0.39

0.39

Fe2+

0.06

0.05

0.06

0.09

0.08

0.09

0.08

0.05

0.07

Fe3+

0.50

0.49

0.50

0.51

0.50

0.51

0.52

0.52

0.51

Mg

0.02

0.02

0.02

0.02

0.01

0.01

0.01

0.02

0.02

Ti

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Cr

1.00

1.00

1.00

1.00

1.00

1.00

1.00

1.00

1.00

Sum

0.04

0.04

0.04

0.05

0.05

0.05

0.04

0.04

0.04

Na

0.86

0.87

0.87

0.87

0.87

0.87

0.87

0.87

0.87

Ca

0.07

0.06

0.06

0.05

0.05

0.05

0.06

0.07

0.06

Mg

0.97

0.97

0.97

0.97

0.97

0.97

0.97

0.98

0.97

Sum

1.85

1.85

1.84

1.82

1.83

1.82

1.83

1.85

1.83

Quad

0.08

0.07

0.08

0.09

0.09

0.09

0.09

0.08

0.09

Jd

0.69

0.94

0.92

0.09

0.30

0.01

0.13

1.16

0.66

Jd%

55.33

54.34

55.13

55.78

54.38

55.89

55.21

57.04

55.54

Di%

38.78

39.87

38.77

37.34

38.94

39.18

38.20

37.50

38.40

Hd%

3.27

2.93

3.26

4.88

4.46

4.92

4.54

2.37

3.98

Ae%

100.00

100.00

100.00

100.00

100.00

100.00

100.00

100.00

100.00

Sum

جدول 1. ادامه.

Table 1. Continued.

KBR18 core

KBR17 rim

KBR16 core

KBR15 core

KBR14 rim

KBR13 core

KBR12 core

KBR1 1 core

KBR10 core

Point No.

50.19

49.99

50.55

50.18

50.14

50.58

50.45

50.54

50.54

SiO2

0.86

0.49

0.75

0.83

0.48

0.69

0.74

0.71

0.71

TiO2

1.73

1.29

1.54

1.69

1.24

1.41

1.58

1.48

1.47

Al2O3

14.34

15.14

14.16

14.38

15.21

14.50

13.97

14.56

14.61

FeO

0.78

0.84

0.78

0.70

0.82

0.82

0.77

0.71

0.83

MnO

9.98

9.28

10.22

9.97

9.39

9.89

10.24

9.97

9.89

MgO

20.75

20.75

20.79

20.83

20.85

20.84

20.84

20.87

20.90

CaO

0.57

0.69

0.50

0.49

0.65

0.53

0.52

0.52

0.55

Na2O

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

K2O

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Cr2O3

99.18

98.46

99.28

99.06

98.77

99.24

99.10

99.35

99.50

Sum

1.96

1.98

1.97

1.97

1.98

1.97

1.97

1.97

1.97

Si

0.02

0.01

0.02

0.02

0.01

0.02

0.02

0.02

0.02

Ti

0.08

0.06

0.07

0.08

0.06

0.06

0.07

0.07

0.07

Al

0.47

0.50

0.46

0.47

0.50

0.47

0.46

0.47

0.48

Fe

0.03

0.03

0.03

0.02

0.03

0.03

0.02

0.02

0.03

Mn

0.58

0.55

0.59

0.58

0.55

0.57

0.59

0.58

0.57

Mg

0.88

0.88

0.87

0.87

0.88

0.89

0.88

0.88

0.87

Ca

0.04

0.05

0.04

0.04

0.05

0.04

0.04

0.04

0.04

Na

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

K

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Cr

4.06

4.06

4.05

4.05

4.060

4.05

4.05

4.05

4.05

Sum

1.94

1.95

1.95

1.94

1.95

1.95

1.94

1.95

1.95

Si

0.06

0.05

0.05

0.06

0.05

0.05

0.06

0.05

0.05

Al(IV)

2.00

2.00

2.00

2.00

2.00

2.00

2.00

2.00

2.00

Sum

0.01

0.01

0.02

0.02

0.01

0.02

0.02

0.01

0.01

Al(VI)

0.39

0.41

0.40

0.41

0.41

0.42

0.39

0.41

0.41

Fe2+

0.09

0.09

0.05

0.06

0.08

0.05

0.06

0.07

0.07

Fe3+

0.49

0.48

0.51

0.49

0.49

0.49

0.51

0.49

0.49

Mg

0.02

0.01

0.02

0.02

0.01

0.02

0.02

0.02

0.02

Ti

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

Cr

1.00

1.00

1.00

1.00

1.00

1.00

1.00

1.00

1.00

Sum

0.04

0.05

0.04

0.04

0.05

0.04

0.04

0.04

0.04

Na

0.86

0.87

0.86

0.86

0.86

0.86

0.86

0.86

0.86

Ca

0.07

0.05

0.08

0.08

0.06

0.07

0.07

0.08

0.07

Mg

0.97

0.97

0.98

0.98

0.97

0.97

0.97

0.98

0.97

Sum

1.826

1.817

1.846

1.843

1.823

1.847

1.842

1.847

1.84

Quad

0.08

0.10

0.07

0.07

0.09

0.08

0.08

0.08

0.08

Jd

0.79

0.47

0.89

0.85

0.27

0.96

0.88

0.82

0.69

Jd%

56.59

53.63

56.99

56.25

54.03

55.38

56.53

54.75

54.71

Di%

38.95

40.95

39.12

39.92

40.85

40.51

37.56

39.33

39.12

Hd%

3.67

4.94

2.99

2.98

4.85

3.15

3.13

3.19

3.56

Ae%

100.00

100.00

100.00

100.00

100.00

100.00

100.00

100.00

100.00

Sum

 

 

 

 

برای رده‌‏‌بندی کلینوپیروکسن‏‌های موجود در استوک پورفیرویید سینیتی منطقه از نمودار دو متغیرة Q-J (Morimoto, 1988) بهره گرفته شده است. بر پایة این رده‌‏‌بندی نمونه‏‌ها در محدودة Quad (آهن- منیزیم-کلسیم) (کلسیک) جای گرفته‏‌اند (شکل 4).

 

شکل 4. نمودار Q-J (Morimoto, 1988).

Figure 4. Q-J diagram (Morimoto, 1988).

بر پایة نمودارهای سه‌تایی Mg-Na-(Fe+2+Mn) (شکل 5-A) و En-Fs-Wo (شکل 5-B)، ترکیب شیمیایی کلینوپیروکسن‏‌ها از نوع دیوپسید ارزیابی می‏‌شود.

برای نام‏‌گذاری فلدسپارها در سنگ‌های منطقة کبریک از نمودار آلبیت NaAlSi3O8 (Ab)-آنورتیت CaAl2Si2O8 (An)- ارتوز KAlSi3O8 (Or) بهره گرفته شده است (Deer et al., 1991). بر پایة این رده‌‏‌بندی، ترکیب آلکالی‏‌فلدسپارها در استوک پورفیرویید سینیتی در گسترة آنورتوکلاز جای می‏‌گیرد (شکل 6). داده‌های به‌دست‌آمده از تجزیة این کانی در جدول 2 نشان داده شده‌اند.

زمین‏‌شیمی سنگ کل

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

بر پایة نسبت Na2O+K2O در برابر SiO2، نمونه‏‌های یادشده در محدودة ترکیبی سینیت جای می‏‌گیرند (شکل 7).

در بررسی الگوی توزیع عنصرهای کمیاب بر پایة نمودارهای عنکبوتی، مقدار عنصرهای کمیاب ناسازگار و خاکی کمیاب نمونه‏‌های منطقه نسبت به ترکیب گرانیت پشتة اقیانوسی (ORG) و کندریت بهنجار شده‏‌اند. الگوی عنصرهای کمیاب بهنجارشده به ORG (گرانیتوییدهای پشتة اقیانوسی) دربارة سنگ‌های منطقه و مقایسه آن با الگوهای رسم‌شده برای سنگ‌های گرانیتوییدی از محیط‏‌های گوناگون زمین‌ساختی (Pearce et al., 1984) نشان می‏‌دهد الگوی ترکیبی سنگ‌های بررسی‌شده با الگوی گرانیتوییدهای درون‌صفحه‏‌ای (WPG) همخوانی خوبی نشان می‏‌دهد (شکل 8-A).

 

شکل 5. ترکیب شیمیایی کلینوپیروکسن‏‌های درون استوک پورفیرویید سینیتی کبریک روی A) نمودار سه‌تایی Mg-Na-(Fe+2+Mn)، (Eby et al., 1998B) نمودار سه‌تایی En-Fs-Wo (Morimoto, 1988).

Figure 5. Chemical composition of clinopyroxenes of Kabbarik syenitic porphyroid stock on A) Mg-Na-(Fe+2+Mn) ternary diagram (Eby et al., 1998); B) En-Fs-Wo ternary diagram (Morimoto, 1988).

 

جدول 2. داده‌های تجزیة نقطه‏‌ای آلکالی‏‌فلدسپارها‏‌ در سینیت کبریک به‌همراه فرمول ساختاری (برپایة 8 اتم اکسیژن) و سازنده‌های پایانی.

Table 2. The EPMA analytical data of alkali feldspars in Kabbarik syenite, the calculated structural formula (based on 6 oxygen) and the end members.

Point No.

PKB1

PKB2

PKB3

SiO2

64.69

64.84

64.69

TiO2

0.06

0.05

0.04

Al2O3

22.39

22.33

22.25

FeO

0.17

0.17

0.18

MnO

0.02

0.00

0.02

MgO

0

0.00

0.01

CaO

1.57

1.37

1.47

Na2O

6.49

6.20

6.2

K2O

5.89

6.32

6.33

Cr2O3

0.01

0.01

0.00

Sum

101.29

101.28

101.20

Si4+

2.86

2.87

2.86

Ti4+

0.00

0.00

0.00

Al3+

1.17

1.16

1.16

Fe2+

0.01

0.01

0.01

Mn2+

0.00

0.00

0.00

Mg2+

0.00

0.00

0.00

Ca2+

0.07

0.06

0.07

Na1+

0.55

0.53

0.53

K1+

0.33

0.36

0.36

Cr3+

0.00

0.00

0.00

Sum

4.99

4.99

4.99

a(ab, Na)

0.82

0.83

0.82

a(an, Ca)

0.1

0.1

0.1

a(Or,K)

0.17

0.15

0.16

 

شکل 6. ترکیب شیمیایی آلکالی‏‌فلدسپارهایِ سینیت کبریک در نمودار سه‌تایی Ab-An-Or (Deer et al., 1991).

Figure 6. Chemical composition of alkali feldspars from Kabbarik syenite in Ab-An-Or ternary diagram (Deer et al., 1991).

 

 

 

 

جدول 3. داده‌های به‌دست‌آمده از تجزیة شیمیایی عنصرهای اصلی (برپایة درصدوزنی)، کمیاب و خاکی کمیاب (برپایة ppm) نمونه‏‌های سینیت کبریک.

Table 3. Results of chemical analysis of major (in wt%), rare, and rare earth (in ppm) elements in the Kabbarik syenite samples.

Sample No.

HAB-01

HAB-02

HAB-03

HAB-04

HAB-05

HAB-06

HAB-07

HAB-08

HAB-09

HAB-10

SiO2

59.84

61.36

59.92

58.94

60.34

59.48

61.12

60.91

58.9

60.52

Al2O3

17.04

15.74

16.42

16.77

15.82

16.78

15.94

15.86

16.72

16.32

Fe2O3

5.72

5.58

5.83

5.62

5.71

5.68

5.71

5.51

5.69

5.78

CaO

2.64

2.47

2.46

2.48

2.08

2.51

2.38

2.34

2.58

2.36

MgO

0.76

0.72

0.66

0.72

0.64

0.75

0.78

0.72

0.84

0.66

Na2O

6.72

6.54

6.67

6.44

6.98

6.64

6.38

6.54

6.45

6.98

K2O

5.72

5.44

5.66

5.62

5.84

5.77

5.32

5.24

5.61

5.57

TiO2

0.83

0.82

0.81

0.86

0.77

0.84

0.75

0.81

0.94

0.85

MnO

0.192

0.193

0.212

0.21

0.214

0.194

0.196

0.214

0.204

0.196

P2O5

0.24

0.23

0.23

0.26

0.21

0.27

0.23

0.22

0.24

0.22

L.O.I.

1.04

0.86

0.94

1.02

0.94

0.85

1.32

1.48

1.88

0.54

Li

12

16

11

15

12

17

15

11

14

16

Be

4.2

4.3

4.4

4.6

3.9

4.2

3.6

3.8

4.4

4.9

S

77

121

78

102

114

88

86

109

92

88

Sc

6.3

6.5

7.4

6.9

6.4

6.8

6.9

7.3

8.1

6.9

V

16

15

17

14

13

17

18

19

14

15

Cr

4

6

3

3

2

5

3

2

2

1

Co

4.3

3.8

4.4

3.9

4.2

4.4

3.7

4.2

5.3

4.5

Cu

34

28

19

31

38

42

37

29

18

31

Ni

5

4

4

6

4

6

5

8

3

5

Zn

108

100

106

98

102

94

102

108

110

112

Ga

161.4

159.8

167.2

173.2

171.4

168.9

170.6

172.4

164.2

174.2

As

4

4.2

4.7

4.4

5.4

4.8

5.3

4.7

5.6

5.1

Se

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

Rb

86

87

90

92

96

88

94

85

92

82

Sr

238.4

204.3

214.3

222.6

234.2

244.6

216.2

225.4

239.4

235.1

 

جدول 3. ادامه.

Table 3. Continued.

Sample No.

HAB-01

HAB-02

HAB-03

HAB-04

HAB-05

HAB-06

HAB-07

HAB-08

HAB-09

HAB-10

Y

37.8

38.2

37.6

38.4

38.3

36.9

37.4

38.2

37.5

37.3

Zr

496

477

489

497

506

492

508

516

485

501

Nb

71.6

75.8

77.1

75.6

82.1

73.2

75.4

79.1

74.9

80.3

Mo

3

5

4

2

5

7

5

4

2

5

Ag

2.6

2.8

2.4

2.4

1.6

1.4

3.4

1.8

1.6

2.8

Cd

0.3

0.2

0.1

0.2

0.1

0.3

0.2

0.2

0.1

0.3

In

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

<0.05

Sn

3.2

3.1

2.8

3.6

2.2

3.1

3.7

2.4

3.6

3.1

Sb

1.3

0.8

1.4

0.9

1.5

0.9

1.1

1.2

1.4

1.6

Te

0.15

0.12

0.14

0.15

0.11

0.18

0.15

0.18

0.12

0.11

Cs

1.8

1.7

2.1

2.3

2.1

2.4

1.7

2.6

2.9

1.8

Ba

476

448

482

471

493

489

482

458

496

477

La

56.2

58.3

58.9

59.1

60.4

57.9

56.8

58.9

60.7

58.1

Ce

134

128

131

129

135

136

129

134

131

130

Pr

12.68

13.45

12.88

12.64

13.28

12.14

13.12

12.95

13.12

12.92

Nd

44.9

47.2

46.5

47.1

44.6

44.9

45.7

46.9

47.1

45.9

Sm

9.12

9.68

9.72

9.76

9.12

9.14

9.36

9.62

9.12

9.58

Eu

2.12

2.26

2.19

2.17

2.12

2.18

2.22

2.13

2.12

2.24

Gd

6.58

7.04

6.66

6.72

6.32

6.45

6.82

6.38

7.04

6.58

Tb

1.23

1.34

1.29

1.32

1.24

1.25

1.26

1.32

1.41

1.32

Dy

7.86

8.51

7.82

8.36

8.44

7.78

8.26

8.48

7.98

8.52

Ho

1.15

1.13

1.16

1.19

1.05

1.17

1.09

1.15

1.26

1.17

Er

5.08

5.36

5.16

5.24

5.32

5.12

5.28

5.36

5.24

5.14

Tm

0.72

0.74

0.76

0.78

0.67

0.73

0.77

0.68

0.76

0.72

Yb

3.94

4.08

4.14

3.96

4.02

3.84

4.08

3.94

4.12

4.14

Lu

0.66

0.65

0.63

0.66

0.71

0.62

0.68

0.64

0.71

0.66

Hf

11.4

11.5

12.3

11.9

11.1

11.4

11.7

10.9

12.4

11.8

Ta

6.61

6.67

6.81

6.54

6.97

6.77

6.58

6.91

6.58

6.62

W

1

1

3

4

2

2

3

6

4

2

Tl

0.33

0.28

0.32

0.24

0.34

0.28

0.22

0.26

0.33

0.35

Pb

3

2

2

5

4

3

5

2

6

4

Bi

<0.1

<0.1

<0.1

<0.1

<0.1

<0.1

<0.1

<0.1

<0.1

<0.1

Th

8.44

8.78

8.34

8.45

8.82

8.34

8.71

8.08

8.56

8.62

U

2.9

3.2

2.8

3.4

3.6

2.8

2.9

3.4

3.9

3.7

 

 

 

شکل 7. ترکیب سینیت‌های منطقة کبریک در نمودار مجموع آلکالی در برابر SiO2  (Middlemost, 1994).

Figure 7. Composition of Kabbarik syenite in total alkali versus SiO2 diagram (Middlemost, 1994).

نمودار عنصرهای خاکی کمیاب بهنجارشده به ترکیب کندریت (Boynton, 1984) برای سنگ‌های سینیتی کبریک در شکل 8-B نشان داده شده‏‌ است. بررسی روند تغییرات میانگین عنصرهای کمیاب بر پایة نمودار‏‌های عنکبوتی بهنجارشده به ترکیب کندریت در سنگ‌های سینیتی کبریک نشان‌دهندة غنی‏‌شدگی عنصرهای خاکی کمیاب سبک (LREE) و عنصرهای ناسازگار نسبت به عنصرهای خاکی کمیاب سنگین (HREE) است.

بحث

برای ارزیابی و بررسی سرشت ماگمایی و جایگاه زمین‌ساختی پیدایش مذاب سازندة استوک پورفیرویید سینیتی منطقة کبریک، ترکیب شیمیایی پیروکسن‏‌ها در نمودارهای دوتایی SiO2-Al2O3 (شکل 9-A) و Ti در برابر Ca+Na (شکل 9-B) بررسی شد. بر پایة این نمودارها، نمونه‏‌ها در گسترة ساب‏‌آلکالن جای دارند.

همچنین، بر پایة نمودار دوتایی Ti + Cr در برابر Ca، نمونه‏‌ها در محدودة کمان آتشفشانی جای دارند (شکل 10).

 

 

شکل 8. نمودار عنکبوتی بهنجارشدة سنگ‌های سینیتی منطقة کبریک به A) ترکیب ORG (Pearce et al., 1984B) ترکیب REE Chondrite (Boynton, 1984).

Figure 8. Normalized spider diagrams for Kabbarik syenitic rocks A) ORG composition (Pearce et al., 1984), B) Chondrite composition (Boynton, 1984).

 

 

شکل 9. ترکیب کلینوپیروکسن‌هایِ نمونه‏‌های سینیتی کبریک در A) نمودار SiO2 در برابر Al2O3 (Le Bas, 1962B) نمودار Ti در برابر Ca+Na (Leterrier et al., 1982).

Figure 9. Composition of clinopyroxenes in the Kabbarik syenitic rocks in A) SiO2 versus Al2O3 diagram (Le Bas, 1962); B) Ti versus Ca+Na diagram (Leterrier et al., 1982).

 

 

 

 

شکل 10. ترکیب کلینوپیروکسن‌هایِ نمونه‏‌های سینیتی کبریک در نمودار Ca در برابر Ti+Cr (Leterrier et al., 1982).

Figure 10. Composition of clinopyroxenes in the Kabbarik syenitic rocks in Ca versus Ti+Cr diagram (Leterrier et al., 1982).

تعیین سری ماگمایی سنگ‌های آذرین از مهم‏‌ترین مباحث در بررسی‏‌های سنگ‌شناسی سنگ‌های آذرین هر منطقه و تعیین محیط زمین‌ساختی آنهاست. همچنین، شناخت نوع سری ماگمایی سازندة سنگ‌ها کلید اساسی در حل مسائل سنگ‌زایی و به‌ویژه جایگاه ژئودینامیکی منطقة کبریک شمرده شود. بر پایة نمودار Th/Yb در برابر Ta/Yb، که برای تمایز سری‏‌های شوشونیتی، کالک‏‌آلکالن و توله‌ایتی پیشنهاد شده است، ماگمای مادر سینیت‏‌های منطقه از نوع کالک‏‌آلکالن شناخته شد (شکل 11).

 

شکل 11. ترکیب نمونه‏‌های سینیتی کبریک در نمودار Th/Yb در برابر Ta/Yb (Pearce, 1983).

Figure 11. Composition of Kabbarik syenitic rocks on Th/Yb versus Ta/Yb diagram (Pearce, 1983).

ترکیب عنصرهای خاکی کمیاب مکملی در کنار دیگر اطلاعات به‌دست‌آمده از مجموعة سنگی به‌شمار می‌رود که امکان پیشنهاد الگوی سنگ‏‌شناختی را برای ما فراهم می‏‌آورد. برپایة این داده‌ها، سنگ‌های سینیتی کبریک مایعات مختلفی هستند که از یک توالی خاص ماگمایی (فرایند تبلوربخشی یا ذوب‌بخشی سنگ‌های خاستگاه) جدایش یافته‌اند؛ زیرا تغییرات و الگوی پراکندگی این عنصرها در سنگ‌های آذرین با ترکیب شیمیایی خاستگاه و درجة ذوب و تعادل بلور- مایع کنترل می‏‌شود که هنگام تکامل ماگما رخ می‏‌دهند. بالابودن عنصرهای LILE در این نمونه‏‌ها و تهی‏‌شدگی از عنصرهای HREE می‏‌تواند پیامد نقش پوسته‏‌ای و یا پیدایش ماگمای سازندة سنگ‌های منطقه از مذاب‌های حاصل از پوستة زیرین و یا پیامد نرخ کم ذوب و درجة جدایش بلوری بسیار بالای ماگمای حاصل از گوشتة بالایی باشد. غنی‏‌شدگی از LREE نسبت به HREE با الگوی کمابیش مسطح HREE نشان‌دهندة ماگماهایی با سرشت کالک‏‌آلکالن هستند (Machado et al., 2005). گفتنی است که شیب نمودار عنکبوتی عنصرهای خاکی کمیاب کاهشیِ شدید نیست و تهی‏‌شدگی شدید از HREE دیده نمی‏‌شود. این موضوع می‏‌تواند گویای گوشتة دگرنهاده توسط سیال‌های فرورانشی باشد. آنومالی منفی Eu نشان‌دهندة جدایش پلاژیوکلاز از ماگما و تحول آن به‌سوی سینیت خواهد بود. بررسی‌های تجربی نشان داده‌اند با خروج 17 درصدحجمی پلاژیوکلاز، 16 درصدحجمی پیروکسن و 2 درصدحجمی ایلمنیت، از یک ماگمای مونزونیتی یک ماگمای سینیتی پدید می‏‌آید (Middlemost, 1994).

برای تعیین محیط زمین‌ساختی سنگ‌های منطقه از نمودارهای دوتایی TiO2 در برابر Al2O3 (شکل 12-A) و Zr/Al2O3 در برابر TiO2/Al2O3 بهره گرفته شد (Muller and Groves, 1993) (شکل 12-B). با توجه به این نمودار‏‌ها، نمونه‏‌های بررسی‌شده در محدودة مرتبط با کمان (AR) جای دارند. همچنین، بر پایة نمودار سه‏تایی Nb*50-Zr*3-Ce/P2O5  (Muller and Groves, 1993)، همة نمونه‏‌های یادشده در محیط پسابرخوردی جای دارند (شکل 12-C).

 

 

 

 

شکل 12. ترکیب نمونه‏‌های سینیتی کبریک در A) نمودار دوتایی TiO2 در برابر Al2O3؛ B) نمودار دوتایی Zr/Al2O3 در برابر TiO2/Al2O3؛ C) نمودار سه‏تایی Nb*50-Zr*3-Ce/P2O5 (Muller and Groves, 1993) (CAP: Continental Arcs; PAP: Postcollisional Arcs).

Figure 12. Composition of Kabbarik syenitic rocks on A) TiO2 versus Al2O3 binary diagram; B) Zr/Al2O3 versus TiO2/Al2O3 binary diagram; C) Nb*50-Zr*3-Ce/P2O5 ternary diagram (Muller and Groves, 1993).

 

برداشت

میزبان توده‏‌های آذرین درونی طارم، بیشتر فوران‌های آتشفشانی زیردریایی (سازند کرج) شامل سنگ‌های آذرآواری و گدازه‏‌های تراکیتی، تراکی‏‌آندزیتی، آندزیتی و الیوین‏‌بازالتی است. به‌طور کلی، روند عمومی ساختارهای زمین‏‌شناسی در منطقة طارم، با دیگر مناطق البرز باختری WNW-ESE همسان است و جای‌گیری و پراکندگی سطحی توده‏‌های آذرین درونی منطقة کبریک نیز از روند یادشده پیروی می‏‌کند. مهم‏‌ترین واحدهای سنگی در کوه‏‌های طارم متشکل از سنگ‌های آتشفشانی و ولکانوکلاستیک ائوسن هستند که بخش بزرگی از واحد کوهستانی را تشکیل داده‌اند. استوک پورفیرویید سینیتی کبریک درون سنگ‌های یادشده تزریق شده است و سنی نزدیک به الیگوسن دارد. بر پایة کانی‏‌شیمی، ماگمای مادر این توده ساب‏‌آلکالن (کالک‏‌آلکالن) بوده است. بر پایة ترکیب شیمیایی پیروکسن‏‌ها، نوع پیروکسن‏‌ها در تودة کبریک دیوپسید و ترکیب آلکالی‏‌فلدسپار، آنورتوکلاز به‌دست آمد. بر پایة نمودارهای عنکبوتی، نبود تهی‏‌شدگی شدید از HREE چه‌بسا گویای گوشتة دگرنهاده توسط سیال‌های فرورانشی باشد. بررسی‏‌ها نشان می‏‌دهند نمونه‏‌های بررسی‌شده در محدودة مرتبط با کمان و در یک محیط پس برخوردی پدید آمده‏‌اند. گمان می‌رود استوک سینیتی کبریک از جدایش بلورین ماگمای مادر توده‏‌های کوارتزمونزونیتی تا گرانودیوریتی گلجین-رشت‏‌آباد (در فاصلة کمتر از استوک سینیتی) پدید آمده باشند. آنومالی منفی Eu نشان‌دهندة جدایش پلاژیوکلاز از ماگما و تحول آن به سینیت خواهد بود. بررسی‌های تجربی نشان داده‌اند با خروج 17 درصدحجمی پلاژیوکلاز، 16 درصدحجمی پیروکسن و 2 درصدحجمی ایلمنیت، از یک ماگمای مونزونیتی یک ماگمای سینیتی پدید می‏‌آید.

 

 

 

1 Differential erosion

Adetunji, A., Olarewaju, V.O., Ocan, O.O., Macheva, L., and Ganevc, V.Y. (2018) Geochemistry and U–Pb zircon geochronology of Iwo quartz potassic syenite, southwestern Nigeria: Constraints on petrogenesis, timing of deformation and terrane amalgamation. Precambrian Research, 307, 125–136. https://doi.org/10.1016/j.precamres.2018.01.015
Aghanabati, A. (2005) Geology of Iran. Geological Survey of Iran (in Persian).
Allen, M.B., Kheirkhah, M., Nill, I., Emami, M.H., and Mcleod, C. (2013) Generation of arc and within-
plate chemical signatures in collision zone magmatism: Quaternary lavas from Kurdistan province,
Iran. Journal of Petrology 54(5), 887- 911. https://doi.org/10.1093/petrology/egs090    
Amiri, M., Ahmad Khalaji, A., Tahmasbi, Z., Santos, S.J., Zarei Sahamieh, R., and Zamanian, H. (2017) Geochemistry, petrogenesis, and tectonic setting of the Almogholagh batholith in the Sanandaj-Sirjan zone, western Iran. Journal of African Earth Sciences, 134, 113–133. https://doi.org/10.1016/j.jafrearsci.2017.06.018
Azizi, H., Kazemi, T., and Asahara, Y. (2017) A- type granitoid in Hasansalaran complex, northwestern Iran: Evidence for extensional tectonic regime in northern Gondwana in the Late Paleozoic. Journal of Geodynamics 108: 56- 72. https://doi.org/10.1016/j.jog.2017.05.003
Berberian, M., and King, G.C.P. (1981) Towards a paleogeography and tectonic evolution of Iran. Canadian Journal of Earth Sciences 18(2): 210- 265. https://doi.org/10.1139/e81-163
Boynton, W.V. (1984) Geochemistry of the rare earth elements: meteorite studies, In: Henderson P., (Ed.), Rare Earth Element Geochemistry. Elsevier 63-114.
Chen, J.F., Xie, Z., Li, H.M., Zhang, X.D., Zhou, T.X., Park, Y.S., Ahn, K.S., and Chen, D.G. (2003) U-Pb zircon ages for a collision-related K-rich complex at Shidao in the Sulu ultrahigh pressure terrane, China. Geochemical Journal, 37(1), 35–46. https://doi.org/10.2343/geochemj.37.35
Davoudian, A., Hamedani, A., Shabanian, N., and Mackizadeh, M. A. (2007) Petrological and geochemical constraints on the evolution of the Cheshmeh-Sefid granitoid complex of Golpaygan in the Sanandaj- Sirjan zone, Iran. Neues Jahrbuch für Mineralogie Abhandlungen, 184(2), 117–129. https://doi.org/10.1127/0077-7757/2007/0085
Deer, W.A., Howie, R.A., and Zussman, J. (1991) An introduction to the Rock-forming minerals. Longman Scientific and Technical.
Delavari, M., Arab Asadi, F., and Mohammadi, A. (2019) Paleozoic magmatism in the southwest of Julfa (northwestern Iran): Geochemical characteristics, U- Pb dating and tectonic setting. Iranian Journal of Petrology, 10(2), 99- 120. https://doi.org/10.22108/ijp.2019.117112.1136  (in Persian).
Derakhshi, M., and Ghasemi, H. (2014) Ordovician- Devonian magmatism in the north of Shahrood: implication for long lived rifting of Paleotethys in eastern Alborz. Iranian Journal of Petrology, 5(18), 105- 122. https://ijp.ui.ac.ir/article_16171_en.html?lang=en (in Persian).
Device, H., and Ghoreishi, A. (1969) Geological map of the Masouleh area 1:100000. Geological Survey and Mineral Exploration of Iran.
Ebrahimi, M., Kouhestani, H., Mokhtari, M., and Feizi, M. (2016) Petrology and geochemistry
of the Aqkand acidic volcanic rocks and perlites, North of Zanjan. Journal of Geoscience,
26(101), 99-110. https://doi.org/10.22071/gsj.2016.41039  (in Persian).
Eby, G.N., Woolley, A.R., Din, V.I.C., and Platt, G. (1998) Geochemistry and petrogenesis of nepheline syenites: Kasungu–Chipala, Ilomba, and Ulindi nepheline syenite intrusions, North Nyasa alkaline province, Malawi. Journal of Petrology, 39(8), 1405-1424. https://doi.org/10.1093/petroj/39.8.1405
Fazlnia, A. (2019) Petrogenesis and tectonic significance of Sardasht syenite–monzonite–gabbro–appinite intrusions, NW Iran. International Journal of Earth Sciences, 108, 49–66. https://link.springer.com/article/10.1007/s00531-018-1641-7
Ghasemi, H., and Dayhimi, M. (2015) Devonian Alkaline Basic Magmatism in Eastern Alborz, North of Shahrood: Evidence for Paleotethys Rifting. Iranian Journal of Geology, 8(32), 1-19. https://dorl.net/dor/20.1001.1.17357128.1393.8.32.3.6  (in Persian).
Ghorbani, M. (2009) Economic geology of Iran. Arian Zamin Publications. https://publish.arianzamin.com (in Persian).
Ghorbani, M.R., Ahmadi, A.R., Tiepolo, M., and Langone, A. (2013) Clinopyroxene
chemistry in Tertiary alkaline volcanic rocks from Taleghan, Central Alborz, Iran:
implications for two parental melts. Chemie der Erde, Geochemistry, 73(4), 565-568. https://doi.org/10.1016/j.chemer.2013.09.003
Jafari, A., Fazlnia, A.N., and Jamei, S. (2018) Geochemistry, petrology and geodynamic setting of the Urumieh plutonic complex, Sanandaj-Sirjan zone, NW Iran: New implication for Arabian and Central Iranian plate collision. Journal of African Earth Sciences, 139, 421–439. https://doi.org/10.1016/j.jafrearsci.2017.11.039
Jafari, A., Moayyed, M., and Williams, I.S. (2020) The age and geochemistry of the Bardkish syenite, northwest Iran: Syenite formation during Neo-Tethyan subduction. Island Arc, 29(1), e12375. https://doi.org/10.1111/iar.12375
Jamshidibadr, M., Collins, A.S., Salomao, G.N., and Costa, M. (2018) U–Pb zircon ages, geochemistry and tectonic setting of felsic and mafic intrusive rocks of Almogholagh complex, NW Iran. Periodico di Mineralogia, 87(1), 21–53. https://doi.org/10.2451/2018PM761
Kaul, P.F.T., and Cordiani, U.G. (2000) Geochemistry of Serra do Mar granitic magmatism and tectonic implications, southern Brazil. Revista Brasileira de Geoẻiencias, 30(1), 115–119. https://repositorio.usp.br/bitstreams
Kherikhah, M., and Karimi, A. (2021) The Geothermobarometric of Uchbolaq trachyandesites lavas, research on the Early-Cenozoic volcanism in Tarom-Hashtjin magmatic zone. Journal of Petrology, 12(4), 101-120. https://doi.org/10.22108/ijp.2022.130813.1254 (in Persian).
Kherikhah, M., and Aghaali, E. (2018) Petrographical investigation of the prismatic columns
rocks in the North of Zanjan (NW Iran). 2nd Trigger International Conference, Tehran,
Iran. https://trgr.ut.ac.ir (in Persian).
Le Bas, M.J. (1962) The role of aluminium in igneous clinopyroxenes with relation to their parentage. American Journal of Science, 260(4), 267-288. https://ui.adsabs.harvard.edu/link_gateway/1962AmJS.260.267L/doi:10.2475/ajs.260.4.267
Leterrier, J., Maury, R.C., Thonon, P., Girard, D., and Marchal, M. (1982) Clinopyroxene composition as a method of identification of the magmatic affinities of palaeo-volcanic series. Earth and Planetary Sciences, 59(1), 139-154. https://doi.org/10.1016/0012-821X(82)90122-4
Machado, A.T., Chemale, Jr.F., Conceicao, R.V., Kawaskita, K., Morata, D., Oteıza, O., and Schmus, W.R.V. (2005) Modeling of subduction components in the Genesis of the Meso- Cenozoic igneous rocks from the South Shetland arc, Antarctica. Lithos, 82(3-4), 435- 453. https://doi.org/10.1016/j.lithos.2004.09.026
Martin, R.F., and Devito, C. (2005) The patterns of enrichment in felsic pegmatites ultimately depend on tectonic setting. The Canadian Mineralogist, 43(6), 2027–2048. https://doi.org/10.2113/gscanmin.43.6.2027
Mazhari, S.A., Bea, F., Amini, S., Ghalamghash, J., Molina, J. F., Montero, P., and Williams, I. S. (2009) The Eocene bimodal Piranshahr massif of the Sanandaj–Sirjan Zone, NW Iran: A marker of the end of the collision in the Zagros orogen. Journal of the Geological Society, 166(1), 53–69. https://doi.org/10.1144/0016-76492008-022
Middlemost, E.A.K. (1994) Naming material in the magma - igneous rock system. Earth Sciences Reviews 37(3-4), 215-224. https://doi.org/10.1016/0012-8252(94)90029-9
Moayyed, M. (2001) Petrological investigations of the Tertiary volcanic-plutonic belt of western Alborz-Azerbaijan with a special view on Hashtjin region. Ph.D. thesis, Shahid Beheshti University, 328 pages, (in Persian).
Moein Vaziri, H. (1997) An introduction to magmatism in Iran. Tarbiat Moalem University Publications, Tehran, 440 pp (in Persian).
Morimoto, N. (1988) The Nomenclature of Pyroxenes. Mineralogical Magazine, 52, 425-433.
Muller, D., and Groves, D.I. (1993) Direct and indirect associations between potassic igneous rocks, shoshonites and gold-copper deposits. Ore Geology Reviews, 8(5), 383-406. https://doi:10.1016/0169-1368(93)90035-W
Nabatian, G. (2008) Geology, geochemistry and genesis of Sorkhe-Dizaj iron oxide apaite deposit,
southeast of Zanjan. MSc thesis, Tarbiat Modares University, Tehran, Iran (in Persian).
Nabatian, G., Ghaderi, M., Neubauer, F., Honarmand, M., Liu, X., Dong, Y., Jiang, S.Y., von Quadt, A., and Bernroider, M. (2014) Petrogenesis of Tarom high-potassic granitoids in the Alborz–Azarbaijan belt, Iran: Geochemical, U–Pb zircon and Sr–Nd–Pb isotopic constraints. Lithos, 184, 324-345. https://doi.org/10.1016/j.lithos.2013.11.002
Nabatian, G., Jiang, S., Honarmand, M., and Neubauer, F. (2016) Zircon U-Pb ages, geochemical and SrNd-Pb-Hf isotopic constraints on petrogenesis of the Tarom-Olya pluton, Alborz magmatic belt, NW Iran. Lithos, 244, 43-58. https://doi.org/10.1016/j.lithos.2015.11.020
Nabavi, M.H. (1977) An introduction to the geology of Iran. Publications of the Geological Organization of the country, 109 pages, (in Persian).
Pearce, J. (1983) Role of the sub-continental lithosphere in magma genesis at active continental magams. Continental Basalts and mantle xenoliths, 1983, 230-249. https://cir.nii.ac.jp/crid/1573387448936656128
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(4), 956-983. https://doi.org/10.1093/petrology/25.4.956
Rezanajad, H. (2013) Petrological investigation of volcanic and plutonic rocks of Kaj Kolah-Zhaker area, Northeast of Zanjan. Conference of Economic Geology Association of Iran.
Sarjoughian, F., Kananian, A., Haschke, M., and Ahmadian, J. (2016) Transition from I-type to A-type magmatism in the Sanandaj–Sirjan Zone, NW Iran: An extensional intra-continental arc. Geological Journal, 51(3), 387–404. https://doi.org/10.1002/gj.2637
Sepahi, A.A., and Athari, S.F. (2006) Petrology of major granitic plutons of the northwestern part of the Sanandaj-Sirjan Metamorphic Belt, Zagros Orogen, Iran: With emphasis on A-type granitoids from the SE Saqqez area. Neues Jahrbuch für Mineralogie Abhandlungen, 183(1), 93–106. https://doi.org/10.1127/0077-7757/2006/0063
Shafaii Moghadam, H., Khademi, M., Hu, Z.C., Stern, R.J., Santos, J.F., and Wu, Y.B. (2015) Cadomian
(Ediacaran- Cambrian) arc magmatism in the ChahJam- Biarjmand metamorphic complex (Iran):
magmatism along the northern active margin of Gondwana. Gondwana Research, 27(1), 439- 452. https://doi.org/10.1016/j.gr.2013.10.014  
Stöcklin, J., and Eftekharnezhad, B. (1969) Geological map of the Zanjan area 1:250000. Geological Survey and Mineral Exploration of Iran.
Verdel, C., Wernicke, B. P., Hassanzadeh, J., and Guest, B. (2011) A Paleogene extensional
arc flare-up in Iran. Tectonics, 30(3), TC3008. https://doi:10.1029/2010TC002809
Whitney, D.L., and Evans, B.W. (2010) Abbreviation for names of rock-forming minerals. American Mineralogist, 95(1), 185-187. DOI: 10.2138/am.2010.3371
Yan, X., and Jiang, S.Y. (2019) Petrogenesis and tectonic implications of Early Cretaceous shoshonitic syenites in the northern Wuyi Mt Range, Southeast China. Journal of Asian Earth Sciences, 180, 103877. https://doi.org/10.1016/j.jseaes.2019.103877
Yang, J.H., Chung, S.L., Wilde, S.A., Wu, F. Y., Chu, M.F., Lo, C.H., and Fan, H.R. (2005) Petrogenesis of post-orogenic syenites in the Sulu Orogenic Belt, East China: Geochronological, geochemical and Nd–Sr isotopic evidence. Chemical Geology, 214(1-2), 99–125. https://doi.org/10.1016/j.chemgeo.2004.08.053