زمین‌شناسی، کانه‌نگاری و ژئوشیمی کانسار آهن 16 ب بافق یزد

نوع مقاله : مقاله پژوهشی

نویسندگان

1 زمین‌شناسی اقتصادی پژوهشکده علوم زمین، سازمان زمین شناسی کشور، تهران، ایران

2 زمین شناسی اقتصادی، گروه زمین شناسی، دانشگاه آزاد اسلامی میانه، ایران

چکیده

مقدمه
پراکندگی کانسارهای آهن در ایران از لحاظ زمانی مربوط به دوره خاص زمین ساختی نیست. ذخایر آهن بافق با 39 کانسار و اندیس معدنی با حدود دو میلیارد تن کانسنگ آهن، یکی از مهمترین مناطق استخراج سنگ آهن در ایران است. کانسارهای آهن سه چاهون، چادرملو و چغارت، هر یک به ترتیب 117، 400 و 216 میلیون تن ذخیره دارند کانسار آهن آنومالی 16 ب بافق در زون ساختاری ایران مرکزی واقع شده است.
مواد و روش‌­ها
با توجه به پیمایش­های صحرایی و مطالعه مغزه­های حفاری که در منطقه انجام شده است. 54 نمونه جهت مطالعات سنگ نگاری و کانه­نگاری، 44 (19 نمونه کانسنگ و 25 نمونه از توده نفوذی و سنگ میزبان) نمونه از مغزه­های حفاری جهت تجزیه شیمیایی عناصر اصلی و کمیاب به روش طیف سنجی پلاسمای جفت شده القایی - طیف سنجی جرمی (ICP-MS) مدل Varian برای 44 عنصر در آزمایشگاه شرکت فرآوری کرج و 5 نمونه از سنگ­های نفوذی میزبان رگه­های کانه­دار به روش فلورسانس اشعه ایکس (XRF) ساخت کمپانی PHILIPS مدل  PW1480در آزمایشگاه فرآوری کرج برای اکسیدهای اصلی تجزیه شده است. طبق تقسیم‌بندی واحدهای ساختمانی- رسوبی ایران، منطقه مورد مطالعه در زون ایران مرکزی واقع شده است (Nabavi, 1355). در این زون قدیمی‌ترین سنگ­های دگرگونی پرکامبرین وجود دارد (Aghanabati, 1384). منطقه معدنی بافق بخشی از ریز صفحه ایران مرکزی(خاوری) است که تکامل این ریز صفحه در ارتباط با کوهزایی کاتانگایی و حرکات تکتونیکی آن طی 600 میلیون سال گذشته است (Taghavi, 2007). مهمترین ذخایر مگنتیت‌های توده‌ای نئوپروتروزوئیک-کامبرین پیشین در سنگ­های آتش‌فشانی بافق در دیاترم­ها (بالا‌آمدگی­های مواد گوشته‌ای) در حاشیه فروافتادگی­های کالدراهای ایگنمبریتی واقع شده است (Förster et al, 1994). ذخایر مربوط به دیاترم­ها که در منطقه معدنی بافق مشخص شده‌اند شامل ذخایر چاه­گز، چادر ملو، چغارت و سه چاهون است. کانی­سازی­های پراکنده آهن و عناصر نادر خاکی در منطقه معدنی بافق با ایجاد ریفت درون قاره­ای ارتباط مستقیم دارد. فعالیت­های آتشفشانی، ماگماتیسم و تکتونیک منطقه تا حدی متاثر از فرایند کافتی شدن است که سبب کانی­سازی در سنگ­های آذرین شده است (Samani, 1993).
کانسار 16ب بخشی از محدوده معدنی بافق و ریز قاره ایران مرکزی، خاوری است. بر اساس ویژگی­های سنگ و چینه­شناختی، منطقه معدنی بافق- ساغند به سه بخش غربی، مرکزی و شرقی تقسیم می شود (Ramezani and Tukker, 2003). قدیمی­ترین سنگ­ها در بخش شرقی (نئوپروتروزوئیک-کامبرین پیشین) و جوانترین سنگ­ها در بخش غربی (ائوسن) قرار گرفته است. گسل نی­باز-چاتک جداکننده دو بخش غربی و مرکزی از یکدیگر است. بخش­های مرکزی و شرقی نیز بوسیله پشت بادام از یکدیگر جدا می­شوند.  
نتایج و بحث
در محدوده کانسار انواع سنگ­های آذرین، دگرگونی وجود دارد. مهم­ترین سنگ­های آذرین منطقه گابرو، دیوریت، سینیت، کوارتز مونزوسینیت،گرانیت و سنگ­های بازیک به شدت دگرسان شده (متابازیت) است. سنگ­های دگرگونی منطقه مرمر و سنگ­های اسکارنی است. گابرو، سینیت، گرانیت، کوارتز مونزوسینیت از سنگ­های آذرینی هستند که در حاشیه محدوده و در نتایج حاصل از مغزه­گیری دیده شده است. سنگ گرانیتی محدوده از نوع لوکوگرانیت و با بافت هیپ ایدیومورفیک و گرافیکی(هم­رشدی فلدسپار پتاسیک و کوارتز) از هم­رشدی بلورهای پلاژیوکلاز به صورت تیغه­های منشوری و فلدسپات سدیک به صورت بلورهای بی­شکل و قطور همراه با دانه­های کوارتز بی­شکل تا نیمه­شکل­دار با حواشی گرد شده یا خلیجی است. مگنتیت، هماتیت، آلکالی فلدسپار­ها و بیوتیت جز کانی­های فرعی و سرسیت، کانی­های رسی، اپیدوت و کربنات­ها از کانی­های ثانویه این طیف سنگ­ها هستند. مقدار کانی­های کدر در این نمونه­ها 8 تا 10 درصد است. در برخی مقاطع بلورها، در اثر فشارهای تکتونیکی، شکسته و خرد شده و بافت کاتاکلاستی را به وجود آورده­اند. کانی­زایی فلزی در منطقه مورد مطالعه در ارتباط با سنگ­های سینیتی، گابرویی و اسکارنی رخ داده است، کانی­شناسی محدوده مورد مطالعه بسیار ساده است و شامل دو دسته کانی­های فلزی و غیرفلزی (باطله) می­باشد. کانی­های فلزی مگنتیت، هماتیت، پیریت و کالکوپیریت است که با کوارتز، اکتینولیت، کلسیت و اپیدوت همراهی می­شود و به شکل­های مختلف در داخل سنگ­های میزبان، توده­های نفوذی سینیتی، گابرویی و اسکارنی دیده شده است. کانی­سازی دارای ساخت توده­ای و بافت­های پرکننده فضای خالی، جانشینی (مارتیتیزاسیون)، عدم آمیزش و افشان از بافت­های غالب در نمونه­ها هستند. بررسی نتایج حاصل از مطالعات تجزیه نمونه­های کانسنگ نشان می­دهد که مقدار اکسیدآهن کل در نمونه­های برداشت شده از 25 تا75 درصد تغییر می­کند. مقدار سیلیس در نمونه­ها از (5 تا 45) درصد متغیر است. میزان تیتانیم در نمونه­های کانسنگ منطقه مورد مطالعه از 1/0 تا 5/0 درصد است. با توجه به روند تغییرات بین تیتان و آهن منفی است. تیتانیوم معمولا در ساختار مگنتیت وارد می­شود (Bin et al, 2016). میزان تغییرات اکسید پتاسیم از حدود 1/0 تا 8/1 درصد متغیر است. مقدار فسفر در نمونه­های مورد آنالیز از 02/0 تا 35/0 متغیر است که این موضوع نشاندهنده عدم وجود کانه­های فسفات­دار در کانه­زایی این کانسار بوده و تغییرات مقدار فسفر نسبت به آهن تا حدودی ثابت است. اکسید منیزیم از 1 تا 12 درصد است. افزایش مقدار اکسید منیزیم ناشی از وجود کانی­های فرومنیزیم از جمله آمفیبول و دولومیت است. روند تغییرات اکسید منیزیم نسبت به آهن دارای شیب منفی است که این نشاندهنده عدم وجود احتمالی جانشینی منیزیم به جای آهن با دو بار مثبت است. تغییرات مقدار اکسید آلومینیوم از 2 تا 12 درصد و مقدار اکسید کلسیم از 2 تا نزدیک 26 درصد است. مقدار عناصر کبالت و نیکل در نمونه­های کانسنگ به ترتیب از 3 تا 75 و 1 تا 17 گرم در تن است. مقدار کروم و وانادیم در محدوده مورد مطالعه به ترتیب 10 تا 96 و 40 تا 120 گرم در تن است. رفتار عناصر نادر خاکی برای تعیین نحوه تشکیل، برای کانسارهای آهن با منشا متفاوت آذرین یا رسوبی دارای اهمیت است (Naslund et al, 2000). در این کانسار همبستگی مثبتی بین عناصر نادر خاکی با P2O5 در منطقه اسکارنی (خارج از منطقه آنومالی) دیده نشده است (Zamanian et al, 2014). برای بهتر نمایان شدن همبستگی­ها می­توان اعداد در یک ضریب ثابت ضرب یا تقسیم نمود (Zamanian et al, 2014). نسبت Eu/Eu* و Ce/Ce* که بنا بر معادلات روابط 1 و 2:
رابطه 1)
Eu/Eu*= {(2Eu) sp/ (Eu)ch}/ {(Sm)sp/(Sm)ch} + {(Gd)sp/(Gd)ch}
رابطه 2)
Ce/Ce*= {(2Ce) sp/ (Ce)ch}/ {(La)sp/(La)ch} + {(Pr)sp/(Pr)ch}
محاسبه می­شود و مقیاسی از بی­هنجاری بوده و مقادیر بالاتر از 1 بی­هنجاری مثبت و مقادیر کمتر از 1 بی هنجاری منفی است. مقادیر محاسبه شده برای Eu/Eu* بین 4/1 تا 4/2 و Ce/Ce* بین 67/0 تا 28/3 در تغییر است. این ارتباط بین Ce و REE بیان کننده ثابت شدن Ce به صورت Ce4+ در شرایط اکسیدی و اسیدی است، در حالی که در شرایط بازی عنصر Ce با لیگاندهای کربنات، تشکیل هم­بافت داده و از گستره شسته شده و حرکت می­کند که باعث بی­هنجاری منفی Ce می­گردد (Gasper et al, 2007). رفتار Eu در نمونه­ها قابل اندازه­گیری نبوده است بر همین اساس نمی­توان راجع به نتایج آن بررسی مناسبی نمود. تفریق بیشتر در بخش LREE از نمودار قابل مشاهده است، جایی که میانگین La/Sm 12 می­باشد، در حالی که میانگین نسبت Gd/Lu در بخش HREE برابر با 95/2 است. الگوی تغییرات REE در این کانسار نشان می­دهد که تحت­تاثیر دگرسانی هیدروترمال قرار گرفته است. مقدار مجموع عناصر نادر خاکی بین 1/13 تا 1/375 گرم در تن است که مشابه میزان گزارش شده برای کانسارهای مگنتیت­دار اسکارنی می­باشد (Gasper et al, 2007)، با توجه به نمودار نرمالایز شده نمونه ها دارای غنی شدگی در LREE و تهی شدگی در HREE بوده و همچنین آنومالی مثبت Eu است (مقدار اندازه­گیری شده در نمونه (P-34)، که این وضعیت نشان دهنده کانسارهای تیپ اسکارن است (Bea et al, 1996). به احتمال زیاد غنی شدگی و تهی شدگی نشان می­دهد که این عناصر به صورت جانشینی در شبکه کانی­های گارنت و زیرکن و همچنین کانه مگنتیت می­توانند حضور داشته باشند.
نتیجه­‌گیری
از دیدگاه زمین­شناسی کهن­ترین سنگ­های این ناحیه، واحد­های دگرگون شده پرکامبرین از نوع گنایس، میکاشیست، آمفیبولیت و میگماتیت بوده که پی سنگ منطقه را تشکیل داده است. محدوده کانی­سازی به شکل کامل بوسیله سری رسوبات ترشیاری وکواترنری حوضه بافق پوشیده شده است. توده نفوذی زیر زون کانه­دار دیوریت­ - سینیت­های قلیایی است. در این کانسار انواع سنگ­های آذرین، دگرگونی وجود دارد. از نظر پتروگرافی سنگ­های آذرین محدوده شامل گابرو، سینیت، کوارتز مونزوسینیت، گرانیت و سنگ­های بازیک به شدت دگرسان شده (متابازیت) است. سنگ­های دگرگونی منطقه از مرمر و سنگ­های اسکارنی است. کانی­زایی فلزی در ارتباط با سنگ­های سینیتی، گابرویی رخ داده است. مگنتیت فراوانترین کانه اکسیدی آهن در محدوده مورد مطالعه است و با ساخت توده­ای، نواری، پرکننده فضاهای خالی و افشان دیده می­شود. دانه­های آن شکل­دار تا بی­شکل است. مگنتیت در نزدیکی سطح زمین اکسید شده و به کانی­های هماتیت، گوتیت و دیگر اکسیدهای آهن دگرسان می­شود. دیگر کانی­های فلزی همراه مگنتیت، پیریت و کالکوپیریت است که با کوارتز، اکتینولیت، کلسیت و اپیدوت به شکل­های مختلف در داخل سنگ­های میزبان، توده­های نفوذی سینیتی، گابرویی و اسکارنی دیده شده است. تطابق نمودارهای عناصر نادر در کانی­سازی آهن­ و توده نفوذی نشاندهنده احتمالی یکسان بودن منشا کانی­سازی با توده نفوذی است. عیار اکسید آهن در کانسنگ بین 25 تا 75 درصد تغییر می­کند. عنصر آهن با اکسیدهای تیتانیم، منیزیم، منگنز، فسفر، پتاسیم و سدیم همبستگی منفی نشان می­دهد. بر اساس رابطه کبالت با نیکل، کروم با نیکل، کروم با وانادیم، این کانسار آهن در محدوده کانسارهای با منشا هیدروترمال قرار می­گیرد. با توجه به نسبتAl/Co  و Sn/Ga این کانسار در محدوده کانسارهای تیپ اسکارنی قرار می­گیرد. بر اساس الگوهای پراکندگی عناصر کمیاب کانسار آهن آنومالی 16 ب با کانسارهای تیپ اسکارنی بیشتر شباهت دارد. شواهد زمین­شناسی، کانی­شناسی و کانه­نگاری و همچنین ژئوشیمیایی کانی مگنتیت نشاندهنده منشاء اسکارنی برای کانی­سازی آهن آنومالی 16ب است. آهن توسط سیالات داغ حاصل از توده­های نفوذی، جابجا شده و در حدفاصل بین واحدهای دگرگونی اسکارن و مرمر تجمع حاصل کرده است.

کلیدواژه‌ها


عنوان مقاله [English]

Geology, Mineralography and geochemistry in 16 B Fe mineralization Bafq (Yazd)

نویسندگان [English]

  • pouria salami 1
  • afshin akbarpour 1
  • mohammad lotfi 1
  • arash gourabjiri 2
1 Economic geology at the Research Institute for Earth Sciences (Geological Survey of Iran), Tehran, Iran
2 Economic geology at Mianeh Azad University, Mianeh, Iran
چکیده [English]

Introduction
The XVI-B anomaly iron ore deposit is situated within the Central Iranian structural zone, specifically in the Bafq region. This region is notable for its lack of a direct association between specific tectonic periods and iron ore deposits. The Bafq area contains 39 iron ore deposits with an estimated total reserve of 2 billion tons, making it one of the most significant iron ore extraction regions in Iran.
 
Materials and Methods
Prominent iron ore deposits in the area include Sechahoon (117 Mt), Chadormalu (400 Mt), and Choghart (216 Mt). The Bafq anomaly iron ore deposit is located within the Central Iranian structural zone. This study involved the collection of drilling core samples for ICP-MS analysis (conducted at Karaj Laboratory), thin-polish and thin-section preparations (54 samples), and XRF analysis (5 samples, also conducted at Karaj Laboratory). According to the structural-sedimentary unit classification of Iran, the study area lies within the central Iranian zone (Nabawi, 1976). This zone contains some of the oldest metamorphic rocks in Iran, dating back to the Precambrian (Aghanbati, 2004). The Bafq mineral district, a subset of the central-eastern Iranian microplate, has experienced tectonic evolution influenced by the Katanga orogeny and related movements over the past 600 million years (Taghavi, 2015). The region's Neoproterozoic–Early Cambrian mass magnetite deposits are predominantly found in volcanic rocks associated with mantle diapirism along caldera margins. Notable deposits include Chaghez, Chadormalu, Choghart, and Sechahoon (Torab et al, 2007). The dispersed mineralization of iron and rare earth elements in the Bafq mining area is directly linked to intracontinental rifting. Volcanism, magmatism, and regional tectonics, influenced by rift dynamics, played a significant role in the mineralization of igneous rocks (Samani, 1993).
The XVI-B deposit is part of the Bafq mineral district, located within the central-eastern Iranian microplate. Stratigraphically, the area is divided into western, central, and eastern sections (Ramezani and Tucker, 2003). The oldest rocks (Neoproterozoic–Early Cambrian) are situated in the eastern part, while the youngest rocks (Eocene) are found in the west. Major faults, including the Nybaz-Chatak and Poshte-Badam faults, delineate these sections.
Results and Discussion
The region contains diverse igneous and metamorphic rocks. The igneous rocks include gabbro, diorite, syenite, quartz monzonite, granite, and highly altered basic rocks (metabasites). The metamorphic rocks are dominated by marble and skarn formations. The leucogranite, characterized by idiomorphic and graphic textures, contains plagioclase and sodic feldspar crystals alongside quartz, alkali feldspar, biotite, and secondary minerals such as sericite, clay minerals, epidote, and carbonates. Tectonic activity has resulted in cataclastic textures in some samples. Mineralization in the area is primarily associated with syenite, gabbro, and skarn rocks. The metallic minerals include magnetite, hematite, pyrite, and chalcopyrite, while non-metallic minerals such as quartz, actinolite, calcite, and epidote are also present.
Mineralization and Geochemical Characteristics
Iron mineralization predominantly occurs as magnetite, which is observed in massive, void-filling, and disseminated forms. Near the surface, magnetite undergoes oxidation, resulting in its transformation into hematite, goethite, and other iron oxides. Associated metallic minerals include pyrite and chalcopyrite, often found with quartz, actinolite, calcite, and epidote in host rocks, intrusive syenites, gabbros, and skarns.The total iron oxide content in the collected samples ranges from 25% to 75%, while silica content varies between 5% and 45%. Titanium concentrations are relatively low, between 0.1% and 0.5%, and exhibit a negative correlation with iron content. Potassium oxide levels range from 0.1% to 1.8%, and phosphorus content varies between 0.02% and 0.35%, indicating an absence of phosphate mineralization. Magnesium oxide levels range from 1% to 12%, largely attributable to the presence of ferromagnesian minerals such as amphibole and dolomite. Negative correlations between magnesium oxide and iron suggest minimal substitution of magnesium for iron in the mineral lattice. Aluminum and calcium oxides range from 2% to 12% and 2% to 26%, respectively. The ore samples contain cobalt (3–75 ppm), nickel (1–17 ppm), chromium (10–96 ppm), and vanadium (40–120 ppm). The behavior of rare earth elements (REEs) indicates hydrothermal alteration, with total REE content varying between 13.1 and 375.1 ppm. Enrichment in light REEs (LREEs) and depletion in heavy REEs (HREEs), alongside a positive Eu anomaly, are consistent with skarn-type deposits (Bea et al, 1996). This pattern suggests that REEs may substitute for elements in garnet, zircon, and magnetite structures.
 
Conclusion
The geological evidence indicates that the oldest rocks in the area are Precambrian metamorphic units, including gneiss, mica-schist, amphibolite, and migmatite, which form the bedrock of the mineralization anomaly. The deposit is covered by Tertiary and Quaternary sediments of the Bafq Basin. The alkaline diorite-syenite intrusive units play a significant role in hosting mineralization.The mineralization comprises a variety of igneous and metamorphic rocks, including gabbro, syenite, quartz monzonite, granite, marble, and skarn. Magnetite is the dominant iron oxide ore and is accompanied by pyrite and chalcopyrite. Oxidation near the surface leads to the formation of secondary iron oxides like hematite and goethite. The geochemical data suggest that the deposit is of hydrothermal origin, with characteristics aligning it with skarn-type deposits. Correlations among trace elements such as Co/Ni, Cr/Ni, and Cr/V, along with ratios like Al/Co and Sn/Ga, further support this classification. The REE patterns, including a positive Eu anomaly and LREE enrichment, are consistent with skarn-type mineralization and provide insights into the role of hydrothermal fluids in the formation process. This study highlights the skarn origin of the XVI-B anomaly iron ore deposit, emphasizing the significance of intrusive magmatic activity and hydrothermal processes in its formation. These findings contribute valuable insights for the geological modeling and economic evaluation of the deposit.

کلیدواژه‌ها [English]

  • Fe mineraliztion
  • Geochemistry
  • Skarn
  • Anomaly 16 B
  • Bafq
  • Yazd
Afzal, P., 2003.Genesis of Iran's iron deposits and valuable rare elements and minerals with them, master's thesis in mining engineering, Islamic Azad University, Tehran South Branch, 168 p (In Persian).
Amirkhani, R., Ebrahimi, M., Mokhtari, M.A. and Azimzadeh, M.A., 2016. Mineralogy, geochemistry, fluid inclusion study and genesis of Homeijan magnetite- apatite minerlization, SW Behabad, Yazd province. Scientific Quaternary Jounral, Geoscience, V. 26(103), p. 83-100.
Azizi, H., Mehrabi, B. and Akabarpour, A., 2009. Genesis of Tertiary Magnetite– Apatite Deposits, Southeast of Zanjan, Iran. Resource Geology, v. 59, p. 330-341.
Baiat, A., 2014. Geochemistry of gabbro-diorite masses of south western Behbad (Yazd province) and their possible role in iron-apatite mineralization. Master thesis, University of Zanjan, Zanjan, Iran.
Bakhtiari, S., 2016. Genetic investigation of XIV Anomaly iron based on rare elements and rare earth elements in Bafq Area, Bafq, Iran. Ph.D thesis, University of Olomtahghighat, Tehran, Iran.
Bankuitz, P., Bau, U. and Seltmann, R., 1997. Magnetite Ore bearing belts of protozoic and lower Paleozoic age in Europe-comparison between the Besser Arc (Germany) and the Bergslagen province (Sweden). Mineral deposits, papuhan Rotterdam, 36 p.
Bajwah, Z.U., Secombe, P.K. and Offler, R., 1987. Trace element distribution, Co: Ni ratios and Genesis of the Big Cadiairon-copper deposit, New South Wales, Australia, Mineralium deposita, v. 22, p. 292-300.
Bea, F., 1996. Residence of REE, Y, Th and U in granites and crustal protoliths, implications for the chemistry of crustal melts. Journal of Petrology, v. 37, p. 521-552.
Bonyadi, Z., Davidson, G.J., Mehrabi, B. and Meffre, S., 2011. Significance of apatite REE depletion and monazite inclusions in the brecciated Se-Chahun iron oxide- apatite deposit, Bafq district, Iran: Insights from paragenesis and geochemistry. Chemical Geology, v. 281, p. 253-269.
Bin, Z., Hong-Fu, Z., Xin-Miao, Z. and Yong-Sheng, H., 2016. Iron isotope fractionation during skarn-type alteration Implications for metal source in the Han-Xing iron skarn deposit, Ore Geology Reviews, v. 4, p. 139-150.
Daliran, F., 2002. Kiruna-type iron oxide-apatite ores and apatitites of the Bafq district, Iran, with an emphasis on the REE geochemistry of their apatites, In Porter, T.M., (ed.), Publishing, p. 303-320.
Daliran, F., Stosch, H.G. and Williams, P.J., 2007. Multistage metasomatism and mineralization at hydrothermal Fe oxide_REE_apatite deposits and ‘apatitites’ of the Bafq district, central east Iran. In: Stanely C. J. eds. Digging Deeper, p. 1501_1504. Proceedings 9thBiennial SGA Meeting Dublin, Ireland.
Dare, S.A.S., Barnes, S.J., Beaudoin, G., Méric, J., Boutroy, E. and Potvin-Doucet, C., 2014. Trace elements in magnetite as petrogenetic indicators. Mineralium Deposita, v. 49(7), p. 785-796.
Darvishzadeh, A., 1992. Geological survey of Iran, Tehran (In Persian).
Emadi, B., 2009. The structural-metallurgical evolution of Zarigan-Sfordi mining field and its role in the formation of deposits of radioactive materials, lead and zinc, iron and apatite and rare earth elements, Master's Thesis, Research institute of Earth Sciences, GSI, 175 p (In Persian).
Förster, H. and Jafarzadeh, A., 1994. The Bafq mining district in Central Iran: A highly mineralized InfraCambrian volcanic field. Econ. Geol., v. 89, p. 1697-1721.
Gasper, M., Knaack, Ch., Meinert, D.L. and Moretti, R., 2007. REE in skarn systems: A LA-ICP-MS study of garnets from the Crown Jewel gold deposit, Geochimica et Cosmochimica Acta, v. 79, p. 185-205.
Haghipour, A., 1977. Geological Map of Biabanak –BafqArea.1: 500,000. Geol. Survev. Iran.
Heidarian, H., Alirezaie, S. and Lentz, D., 2017. Chadormalu Kiruna-type magnetite-apatite deposit, Bafq district, Iran: Insights into hydrothermal alteration and petrogenesis from geochemical, fluid inclusion, and sulfur isotope data, Ore Geology Reviews, v. 83, p. 43-62.
Hooshmandzadeh, A.R., 1989. An introduction to the geology of Biabank-Bafaq area, Tehran Ministry of Mines and Metals (In Persian).
Iranmanesh, S., 2013. Investigating the distribution, distribution and mineralogy of rare earth elements in iron ore-sfordi apatite, Bafaq area of ​​Yazd, Master's Thesis, Research institute of Earth Sciences, GSI, Tehran (In Persian).
Jafarzadeh, A., Gorbani, M. and Pezeshkpour, M., 1996. Iron deposits, book editing plan, Geological Organization of the country, v. 26, 213 p (In Persian).
Kouhestani, H., Abrahimi, M. and Shahidi, A., 2016. Investigating the type and origin of iron mineralization in the copper mineral occurrence using lithology, mineralogy and geochemical data, Economic Geology, v. 7(1), p. 111-127 (In Persian).
Moore, F. and Modabberi, S., 2003. Origin of Choghart iron oxide deposit, Bafq District, Central Iran: New isotopic and geochemical evidence. Journal of Sciences, Islamic Republic of Iran, v. 14(3), p. 259-269.
Meinert, L.D., 1992. Skarns and skarn deposits. Geosciences Canada, v. 19(4),p. 145-162.
Moghtaderi, A., 2005. Geology and geochemistry of Chador Melo iron mine, southeast of Yazd, north of Shazaq Bafaq, PhD thesis, Shiraz University (In Persian).
Naslund, H.R., Aguirre, R., Dobbs, F.M., Henriquez, F.J. and Nystrom, J.O., 2000. The origin, emplacement, and eruption of ore magmas, IX Congreso Geologico Chileno, Sociated geologica de chile, v. 2, p. 135-139.
Nabavi, M.H., 1977. An introduction to the geology of Iran, Geological survey of Iran (In Persian).
Rahmani, Sh., Amami, H. and Mokhtari, M., 2004. Mineralization of rare earth elements in Bafaq area, Poshteh Badam, 22nd Earth Sciences conferences. Tehran (In Persian).
Samani, B., 1993. Saghand Formation, a riftogenic unit of upper Precambrian in central Iran. Geosciences: Scientific Quarterly Journal of the Geological Survey of Iran 2 (in Persian with English abstract), 32-45 p.
Torab, F.M. and Lehmann, B., 2007. Magnetite-apatite deposits of the Bafq district, Central Iran, apatite geochemistry and monazite geochronology. Min. Mag, v. 71, p. 347-363.
Taghavi, A., 2007. Study and investigation of mineralogy and geochemistry of northern iron anomaly (North Bafq) and determination of the next iron grade" master's thesis, Tarbiat Moalem University, Tehran (In Persian).
Williams, P.J., Barton, M.D., Johnson, D.A., Fontbote, L., De Haller, A., Mark, G., Oliver, N.H.S. and Marschik, R., 2007. Iron oxide copper-gold deposits: Geology, space-time distribution, and possible modes of origin, Economic Geology, v. 100, p. 371-405.
 Zamanian, H. and Radmard, K., 2014. The geochemistry of REEs in Baba Ali skarn deposit, key point for finding mineralization condition, Journal of Crystallography and Mineralogy of Iran, p. 743-758.
Zhiwei, B., Zhenhua, Z., Jayanta, G., Anthony, E. and Jones, W., 2004. HFSE, REE and PGE geochemistry of three sedimentary rock-hosted disseminated gold deposits in southwestern Guizhou Province, China, Geochemistry Journal, v. 38, p. 363-381.