• صفحه اصلی
  • ارزیابی و پهنه‌بندی خطرهای مرتبط با فوران احتمالی آتشفشان سبلان

اشتراک گذاری

آدرس مقاله


کد مقاله : 1402070944252 بازدید : 1545 صفحه: 68 - 80

20.1001.1.17357128.1401.17.66.5.5

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

ارزیابی و پهنه‌بندی خطرهای مرتبط با فوران احتمالی آتشفشان سبلان

محورهای موضوعی :

احمد عباس نژاد 1 , احمد خیاط زاده 2 , حجت الله رنجبر 3 , حمید احمدی پور 4 , بهنام عباس نژاد 5

1 - دانشگاه شهید باهنر کرمان
2 - گروه زمین شناسی، دانشگاه شهید باهنر کرمان، کرمان، ایران
3 - گروه مهندسی معدن، دانشگاه شهید باهنر کرمان، کرمان، ایران
4 - گروه زمین شناسی، دانشگاه شهید باهنر کرمان، کرمان، ایران
5 - گروه زمین شناسی، دانشگاه شهید باهنر کرمان، کرمان، ایران

تاریخ دریافت : 1402/07/09 تاریخ پذیرش : 1402/07/09 تاریخ انتشار : 1402/06/06

کلید واژه: آتشفشان سبلان, خطر لاهار, خطر گدازه, خطر ابر سوزان, خطر خاکستر آتشفشانی,

چکیده مقاله :

استراتوولکان سبلان یکی از آتشفشان‌های نیمه فعال ایران است و احتمال فوران آن در آینده منتفی نیست، زیرا که دارای چشمه‌های آبگرم فراوان و مورفولوژی سرپا (با فرسایش کم) است و از آخرین فوران آن، در مقیاس زمین‌شناسی، زمان اندکی طی شده است. براساس رژیم فوران این نوع آتشفشان‌ها، و نیز، فوران‌های قبلی حاصل از مطالعات زمین‌شناسی ، فوران احتمالی این آتشفشان می‌تواند با تهدیدهای خاکستر (تفرا)، گدازه، پیروکلاستیک‌های جریانی و لاهار همراه باشد. شدت فوران احتمالی آن در مقیاس VEI در حد سه تا پنج پیش‌بینی می‌شود. درصورتی¬که این آتشفشان نشانه‌های شروع فوران را بروز دهد باید براساس نقشه‌های پهنه¬بندی مربوط به خطرهایی که می‌تواند ایجاد کند اقدامات اضطراری صورت گیرد. بنابراین، تهیه این نقشه‌ها برای اقدامات اضطراری ضروری است. در این مطالعه، از مدل ارتفاعی رقومی آتشفشان، تصاویر ماهواره‌ای و نرم‌افزارهای Arc GIS، ENVI و VORIS و داده‌های جوی پایگاه NCEP/NCAR استفاده شد. برای تهیه نقشه خطر پیروکلاستیک‌های جریانی از مدل مالین و شریدان (Malin and Sheridan, 1982) استفاده شد. در تهیه نقشه گسترش و موقعیت جریان‌های گدازه مدل شبیه‌سازی جریان گدازه مورد استفاده قرار گرفت، و بالاخره، برای تعیین مسیرهای احتمالی و موقعیت جریان‌های لاهاری از تلفیق مدل ارتفاعی رقومی آتشفشان و تصویرهای ماهواره‌ای موزائیک شده، آبراهه‌های اصلی شناسائی و بافرگذاری شدند. براساس این مطالعه، خاکستر آتشفشان به سمت شرق حرکت کرده و شهر اردبیل و آبادی‌های دیگری را مورد تهدید قرار خواهد داد. لاهارهای آن هم کاربری‌های متعددی در اطراف مخروط را تهدید می‌کنند و امکان ورود به شهرهای مشکین‌شهر و اردبیل را هم دارند. جریان‌های گدازه اکثر جان¬پناه‌ها، تأسیسات استفاده از آبگرم و پیست اسکی آلوارس را تهدید خواهند کرد. همچنین، مشکین‌شهر، روستای موئیل، پیست‌های اسکی و جان¬پناه‌ها در معرض پیروکلاستیک‌های جریانی (ابر سوزان) هم قرار خواهند داشت.

چکیده انگلیسی:

Sabalan is one of dormant stratovolcanoes of Iran with likely eruption, because there are many thermal springs around it. It has a slightly- eroded cone and, geologically, a short time has passed from its last eruption. According to eruption regime of such volcanoes, as well as its former activities, Sabalan volcanic eruptions may happen along with such hazards as tephra, lava, pyroclastic flows and lahars. Its eruption intensity on VEI scale is anticipated to be 3 to 5. Whenever this volcano shows awakening evidence, it would be necessary to take appropriate actions according to hazard zoning maps which are the subject of this study. That is, these maps are necessary for taking emergency actions. In this study, digital elevation data (DEM); satellite pictures; Arc GIS, ENVI and VORIS softwares; as well as atmospheric data of NCEP/NCAR center were used. For construction of pyroclastic density current spreads, the Malin and Sheridan (1982) model was employed. For preparation of lava flow hazard map, simulation model was used, and for preparation of lahar hazard zones, major valleys on the cone and its surroundings were identified and buffered using DEM and satellite data. According to this study, volcanic ashes will move eastward and threat Ardabil city and several nearby villages. Lahars would threaten many surrounding land uses and probably Meshkinshahr and Ardabil cities. Lava flows would damage many mountaineering shelters and ski facilities of Alvares. Additionally, Meshkinshahr, Moil, ski resorts as well as many shelters are at the threat of nuee ardentes.

منابع و مأخذ:

اسفندیاری درآباد، ف. و خیام، م.، 1386. تحلیلی بر اثرات ژئومورفیک برفاب در دامنه شرقی سبلان، پژوهش‌های جغرافیائی، 6، 83-97.
امامی، م.، 1379. ماگماتیسم در ایران، انتشارات سازمان زمین‌شناسی و اکتشافات معدنی کشور، 608.
باباخانی، ع.، 1387. کوه آتشفشان سبلان، گزارش سازمان زمین‌شناسی کشور، مرکز تبریز، 122.
بیاتی خطیبی، م.، 1386. تحلیل و بررسی نقش عوامل توپوگرافی و دینامیک رودخانه‌ای بر اندازه مخروط واریزه‌ای؛ مطالعه موردی: دامنه شمال غربی سبلان (شمال غرب ایران)، پژوهش‌های جغرافیائی، 6، 157- 175.
خیاط زاده، ا. و عباس نژاد، ا.، 1395. کاربرد روش‌های تصمیم‌گیری EN- SAW و ANP در مطالعات زمین‌شناسی: ارزیابی و اولویت‌بندی خطر فعال شدن استراتوولکان‌های ایران به‌عنوان مطالعه موردی، فصلنامه علوم زمین، (102)26، 137- 146.
خیاط زاده، ا.، عباس نژاد، ا. و رنجبر، ح.، 1395. پهنه‌بندی خطر تفرا، گدازه و ابرهای سوزان احتمالی آتشفشان تفتان، فصلنامه جغرافیا و مطالعات محیطی، (18)5، 17- 30.
خیاط زاده، ا.، عباس نژاد، ا.، رنجبر، ح.، احمدی پور، ح. و عباس نژاد، ب.، 1402. ارزیابی و پهنه‌بندی خطر تفرا، گدازه و ابرهای سوزان ناشی از فوران احتمالی آتشفشان بزمان. زمین‌شناسی کاربردی پیشرفته, 13(2)، 455-468.‎
خیاط زاده، ا.، عباس نژاد، ا.، رنجبر، ح.، احمدی پور، ح. و عباس نژاد، ب.، 1402. پهنه‌بندی خطر تفرا، گدازه، ابر سوزان و لاهار ناشی از فوران احتمالی آتشفشان دماوند. فصلنامه علمی علوم زمین, 32(4)، 355-366.‎
دلال اوغلی، ع. و رجبی، م.، 1394، بررسی اختلافات ساختاری و فرسایشی در مخروط‌های استراتوولکانوسهند و سبلان، اولین همایش ملی علوم زمین و توسعه شهری.
دلال اوغلی، ع.، 1382. بررسی مورفولوژی و نحوه فعالیت یخچال‌های سنگی دامنه شمالی کوه سبلان، پژوهش‌های جغرافیائی، 45، 1- 12.
رحیمی، ع.، 1375. زمین‌شناسی میدان زمین گرمائی شمال غرب سبلان (جنوب مشکین‌شهر) و ژئوشیمی دگرسانی و ته نشست‌های گرمابی منطقه، پایان‌نامه کارشناسی ارشد دانشگاه تهران، 284.
زماني، ر.، امامي، م.، وثوقي عابديني،م. و كريم زاده ثمرين، ع.، 1394. مطالعه شيمي بلور و منطقه بندي کلينوپيروکسن هاي موجود در سنگ‌های آتشفشاني آلکالن شمال غرب مشکین‌شهر، ايران، فصلنامه زمین‌شناسی ایران، (33)9، 31-44.
شاه بیک، ا.، 1372. آب‌های معدنی و گرم ایران، سازمان زمین‌شناسی کشور، 212.
فتح الهی، م. و خیرخواه، م.، 1394. منشاء و جایگاه تکتونو ماگمائی سنگ‌های آتشفشانی کواترنری سبلان، فصلنامه کواترنری ایران، جلد 1، 2، 125- 136.
فتح الهی، م.، امامی، م. و خیرخواه، م.،1392. پتروگرافی و پترولوژی سنگ‌های آتشفشانی ماگماتیسم نهایی سبلان، فصلنامه زمین‌شناسی محیط‌زیست (7)23، 51-68.
فتح الهی، م.، امامی، م. و خیرخواه، م.، 1389. تاثیر فرایند اختلاط ماگمائی و آلایش پوسته‌ای در پیدایش سنگ‌های آتشفشانی سبلان، بیست و نهمین گردهمایی علوم زمین.
قربانی، م.، 1382. آتشفشان شناسی با نگرشی بر آتشفشان‌های ایران، نشر آرین زمین، 362.
قلمقاش، ج.، کتابی، ز.، اصفهانی، ا. و موسوی، ز.، 1398. مقایسه زمین‌شناسی، پتروگرافی، زمین‌شیمی و منشاء ایگنیمبریت در آتشفشان‌های سبلان و سهند، زمین‌شناسی کاربردی پیشرفته، (131)9، 17- 24.
کمالی نژاد، م. و عباس نژاد، ا.، 1398. مطالعه چشمه‌های آبگرم و معدنی ایران (ویژگی‌های شیمیائی، درمانی و گردشگری)، انتشارات سنجش و دانش، 435.
مبشر گرمی، م.، احمدزاده، غ. و آقازاده، ا.، 1397، برآورد ویژگی‌های فیزیکوشیمیایي و محیط تکتونیکی تشکیل توده‌های بازیک جنوب شهرستان گرمی (استان اردبيل) بر اساس شيمي بلورهای آمفیبول و پلاژیوکلاز، )، فصلنامه زمین‌شناسی ایران، (12) 48، 81-98.
یادگاری، ن.، علوی، غ. و مؤید، م.، 1402، بررسی زمین‌ساخت، داده‌های گسلی و ارتباط آن‌ها با کانه‌زایی و دگرسانی در گستره اسبخان هریس (استان آذربایجان‌شرقی - شمال‌غرب ایران)، فصلنامه زمین‌شناسی ایران، (17)65، 33-47.
Abdollahzadeh Bina, F., 2009. Geothermal resource assessment of the NW- Sabalan Field, Iran, Through well testing, Geothermal training programe, United Nations University (UNU), 6, 15- 44.
Alatorre-Ibargüengoitia, M.A., Morales-Iglesias, H., Ramos-Hernández, S.G., Jon-Selvas, J. and Jiménez-Aguilar, J.M., 2016. Hazard zoning for volcanic ballistic impacts at El Chichón Volcano (Mexico). Natural Hazards, 81(3), 1733-1744.
Amini, B., 1994. Geological map of Meshkinshahr, Geological Survey of Iran.
Blong, R. J., 1984. Volcanic hazards: a sourcebook on the effects of eruptions, Elsevier.
Calder, E., Wagner, K. and Ogburn, S., 2015. Volcanic hazard maps. Global volcanic hazards and risk, 335-342.
Capra, L., Norini, G., Groppelli, G., Macías, J.L. and Arce, J.L., 2008. Volcanic hazard zonation of the Nevado de Toluca volcano, México. Journal of Volcanology and Geothermal Research, 176(4), 469-484.
Charlton, D., 2018. New approaches to volcanic hazard mapping at Campi Flegrei, Southern Italy (Doctoral dissertation, UCL (University College London)).
Chevrel, M.O., Favalli, M., Villeneuve, N., Harris, A.J., Fornaciai, A., Richter, N., Derrien, A., Boissier, P., Di Muro, A. and Peltier, A., 2021. Lava flow hazard map of Piton de la Fournaise volcano. Natural Hazards and Earth System Sciences, 21(8), 2355-2377.
Connor, C. B., Hill, B. E., Winfred, B., Franklin, N. W. and La Femina, P. C., 2001. Estimation of volcanic hazards from tephra fallout, Natural Hazards Review, 2, 33- 42.
Didon, J. and Gemain, Y. A., 1976. La Sabalan, volcan plio- Quaternaire de l, Azerbaijan oriental Iran, Etude Geologique et Scientifique et Medical de Grenoble, France.
Emami, M., 1994. Geological Map of Iran, Meshkinshahr quadrangle (1: 100000), Geological Survey of Iran.
Felpeto, A., 2009. VORIS, a GIS based tool for volcanic hazard assessment, User Guide, Observatorio Geofisico Central, IGN.
Felpeto, A., Marti, J. and Ortiz, R., 2007. Automatic GIS- based system for volcanic hazard assessment, Journal of volcanology and Geothermal Research, 166, 106- 116.
Fotouhi, M., 1995. Geothermal development in Sabalan, Iran, in: proceedings of world geothermal congress, Italy, 191- 196.
Ghaedrahmati, R., Moradzadeh, A., Fathianpour, N., Lee, S. K. and Porkhial, S., 2013. 3D inversion of MT data from Sabalan geothermal field, Ardabil, Iran, Journal of Applied Geophysics, 93, 12- 24.
Ghalamghash, J., Mousavi, S. Z., Hassanzadeh, J. and Schmitt, A. K., 2016. Geology, zircon geology, and petrogenesis of Sabalan volcano (northwestern Iran), Journal of Volcanology and Geothermal Research, 327, 192- 207.
Haynes, K., Barclay, J. and Pidgeon, N., 2008. The issue of trust and its influence on risk communication during a volcanic crisis. Bulletin of Volcanology, 70, 605-621.
Latter, J. H. (ED.), 2013. Volcanic hazards: assessment and monitoring, Springer Science and Business Media.
Loughlin, S. C., Sparks, R. S. J., Sparks, S., Brown, S. O K., Jenkins, S. F. and Vye- Brown, C., (Eds). 2015. Global volcanic Hazards and risk, Cambridge University Press.
Malin, M. C. and Sheridan, M. F., 1982. Computer- assisted mapping of pyroclastic surges, Science, 217 (4560), 637- 640.
Mastin, L. G., Guffanti, M., Servrankx, R., Webley, P., Barsotti, S. and Dean, K. et al., 2009. A multidisuplinary efforf to assign realistic source parameters to models of volcanic ash- cloud transport and dispersion during eruptions, Journal of volcanology and Geothermal Research, 186, 10- 21.
Mc Guire, W. J., 1998. volcanic hazards and their mitigation, Geological Society, London, Engineering Geology Special Publications, 15 (1), 79- 95.
Michaud-Dubuy, A., Carazzo, G. and Kaminski, E., 2021. Volcanic hazard assessment for tephra fallout in Martinique. Journal of Applied Volcanology, 10(1), 1-20.
Mousavi, S. Z., Darvishzadeh, A., Ghalamghash, G. and Vosoughi Abedini, M., 2014. Volcanology and geochronology of Sabalan volcano, the hieghest stratovolcano in Azerbaijan region, NW Iran, Nautilus, 128, 85- 98.
Nakasuji, A. and Satake, j., 2004. Volcanic Hazard Map An Introduction and Overseas Cases. Journal of the Japan Society of Engineering Geology, 44(6), 341-348.
Negro, C.D., Cappello, A., Neri, M., Bilotta, G., Hérault, A. and Ganci, G., 2013. Lava flow hazards at Mount Etna: constraints imposed by eruptive history and numerical simulations. Scientific Reports, 3(1), 1-8.
Neri, M., Le Cozannet, G., Thierry, P., Bignami, C. and Ruch, J., 2013. A method for multi-hazard mapping in poorly known volcanic areas: an example from Kanlaon (Philippines). Natural Hazards and Earth System Sciences, 13(8), 1929-1943.
Salvatici, T., Di Roberto, A., Di Traglia, F., Bisson, M., Morelli, S., Fidolini, F., Bertagnini, A., Pompilio, M., Hungr, O. and Casagli, N., 2016. From hot rocks to glowing avalanches: Numerical modelling of gravity-induced pyroclastic density currents and hazard maps at the Stromboli volcano (Italy). Geomorphology, 273, 93-106.
Tarigan, A.P.M., Suwardhi, D., Fajri, M.N. and Fahmi, F., 2017. March. Mapping a volcano hazard area of Mount Sinabung using drone: preliminary results. In IOP Conference Series: Materials Science and Engineering 180(1), 012277.
Thompson, M. A., Lindsay, J. M. and Leonard, G. S., 2018. More than meets the eye: volcanic hazard map design and visual communication. Observing the volcano world: volcano crisis communication, 621-640.
Tilling, R. J., 1984. Volcanic hazards: a sourcebook on the effects of eruptions, Elsevier.
Wagner, K., Ogburn, S.E. and Calder, E.S., 2015. Volcanic Hazard Map Survey.
Walker, G. P., 1982. Volcanic hazards, Interdisciplinary Science Reviews, 7 (2), 148- 157.
Xu, J., Wa, Y., Wa X., Pan, B., Ya, H., Zh. B. and Ya W., 2022. Review on the development of volcanic hazard zonation in china. Geology and Resources, 1;31(3).
متن کامل:


 

ارزیابی و پهنه‌بندی خطرهای مرتبط با فوران احتمالی آتشفشان سبلان

 

احمد عباس نژاد(*و1)، احمد خیاط زاده2، حجت اله رنجبر3، حمید احمدی پور4، بهنام عباس نژاد5

1. استاد گروه زمین‌شناسی، دانشکده علوم، دانشگاه شهید باهنر کرمان، کرمان، ایران

2. دانش­آموخته کارشناسی ارشد، دانشکده علوم، گروه زمین‌شناسی، دانشگاه شهید باهنر کرمان، کرمان، ایران

3. استاد گروه مهندسی معدن، دانشکده فنی و مهندسی، دانشگاه شهید باهنر کرمان، کرمان، ایران

4. استاد گروه زمین‌شناسی، دانشکده علوم،  دانشگاه شهید باهنر کرمان، کرمان، ایران

5- پژوهشگر پسادکترا، گروه زمین‌شناسی، دانشکده علوم، دانشگاه شهید باهنر کرمان، کرمان، ایران

چکیده

استراتوولکان سبلان یکی از آتشفشان‌های نیمه فعال ایران است و احتمال فوران آن در آینده منتفی نیست، زیرا که دارای چشمه‌های آبگرم فراوان و مورفولوژی سرپا (با فرسایش کم) است و از آخرین فوران آن، در مقیاس زمین‌شناسی، زمان اندکی طی شده است. براساس رژیم فوران این نوع آتشفشان‌ها، و نیز، فوران‌های قبلی حاصل از مطالعات زمین‌شناسی ، فوران احتمالی این آتشفشان می‌تواند با تهدیدهای خاکستر (تفرا)، گدازه، پیروکلاستیک‌های جریانی و لاهار همراه باشد. شدت فوران احتمالی آن در مقیاس VEI در حد سه تا پنج پیش‌بینی می‌شود. درصورتی­که این آتشفشان نشانه‌های شروع فوران را بروز دهد باید براساس نقشه‌های پهنه­بندی مربوط به خطرهایی که می‌تواند ایجاد کند اقدامات اضطراری صورت گیرد. بنابراین، تهیه این نقشه‌ها برای اقدامات اضطراری ضروری است. در این مطالعه، از مدل ارتفاعی رقومی آتشفشان، تصاویر ماهواره‌ای و نرم‌افزارهای Arc GIS، ENVI  و VORIS و داده‌های جوی پایگاه NCEP/NCAR استفاده شد. برای تهیه نقشه خطر پیروکلاستیک‌های جریانی از مدل مالین و شریدان (Malin and Sheridan, 1982) استفاده شد. در تهیه نقشه گسترش و موقعیت جریان‌های گدازه مدل شبیه‌سازی جریان گدازه مورد استفاده قرار گرفت، و بالاخره، برای تعیین مسیرهای احتمالی و موقعیت جریان‌های لاهاری از تلفیق مدل ارتفاعی رقومی آتشفشان و تصویرهای ماهواره‌ای موزائیک شده، آبراهه‌های اصلی شناسائی و بافرگذاری شدند. براساس این مطالعه، خاکستر آتشفشان به سمت شرق حرکت کرده و شهر اردبیل و آبادی‌های دیگری را مورد تهدید قرار خواهد داد. لاهارهای آن هم کاربری‌های متعددی در اطراف مخروط را تهدید می‌کنند و امکان ورود به شهرهای مشکین‌شهر و اردبیل را هم دارند. جریان‌های گدازه اکثر جان­پناه‌ها، تأسیسات استفاده از آبگرم و پیست اسکی آلوارس را تهدید خواهند کرد. همچنین، مشکین‌شهر، روستای موئیل، پیست‌های اسکی و جان­پناه‌ها در معرض پیروکلاستیک‌های جریانی (ابر سوزان) هم قرار خواهند داشت.

 

واژه­های کلیدی:  آتشفشان سبلان، خطر لاهار، خطر گدازه، خطر ابر سوزان، خطر خاکستر آتشفشانی

 


مقدمه

پهنه­­­بندی مهم­ترین و اساسی‌ترین اقدام برای کاهش ریسک خطرهای زمین‌شناختی است، زیرا براساس آن وجود یا فقدان خطر و محدوده‌های با درجه­های مختلف خطر مشخص می‌شوند. مقایسه یا همپوشانی این نقشه‌ها با کاربری زمین موقعیت و نوع عناصر در معرض ریسک را مشخص می‌کند. درواقع، این نقشه‌ها پیش نیاز اصلی برای هرگونه اقدامی در مورد خطرها می‌باشند. این نکته در مورد آتشفشان‌ها هم صدق می‌کند. آن‌ها به شکل‌های مختلفی می‌توانند باعث خسارات جانی و مالی شوند و در منابع مختلف از جمله Walker (1982), Blong (1984), Tilling (1989), McGuire (1998), Latter (2013), Loughlin et al., (2015) معرفی شده‌اند. در هر حال، ازجمله مهم­ترین تهدیدهای ناشی از آن‌ها می‌توان به خطر خاکستر و گازهای آتشفشانی، خطر ابر سوزان، ورود گدازه، خطر لاهار و خطر تغییرات آب و هوایی اشاره کرد. براساس بررسی‌های خیاط زاده و عباس نژاد (1395) می‌توان آتشفشان‌های دماوند، سبلان، بزمان و تفتان را خاموش یا نیمه فعال2 فرض کرد، امکان فعال شدن آن‌ها در هر لحظه وجود دارد. در صورت فعال شدن ناگهانی، نیاز به داشتن اطلاعات در مورد محدوده‌های گسترش انواع تهدیدهای آتشفشانی جهت انجام اقدامات مدیریتی است و باید با توجه به نقشه‌های پهنه­بندی انجام شوند.

ناکاسوجی و ساتاکه (Nakasuji and Satake, 2004) نقشه پهنه­بندی خطر آتشفشان را نقشه‌ای تعریف کرده‌اند که موقعیت محدوده‌هایی که در اثر فعالیت آتشفشان آسیب می‌بینند را مشخص سازد. تهیه نقشه‌های پهنه‌بندی خطر آتشفشان‌های فعال یا بالقوه فعال گام مهمی به سمت کاهش ریسک ناشی از آن‌ها است (Tilling, 2005). از آنجا که در دوران تاریخی در ایران هیچ فوران آتشفشانی رخ نداده است، در مقایسه با دیگر خطرها (زلزله و سیل)، توجه کافی به این نوع تهدید نشده است (خیاط زاده و همکاران، 1401 و 1402). در 150 سال گذشته صدها نقشه پهنه‌بندی خطر متعلق به آتشفشان‌ها و پهنه­های آتشفشانی مختلف دنیا تهیه و ارائه شده است که نقش مهمی در مراوده3 این خطر داشته‌اند (Charlton, 2018). جزئیات این نقشه‌ها متأثر از نوع و شرایط آتشفشان، نوع داده‌های مورد استفاده و هدف از تهیه آن‌ها بوده است (Thompson et al., 2017). در هر حال، اسامی بعضی از نقشه‌های پهنه­بندی این نوع خطر در کشورهای مختلف با ذکر منبع و مشخصات نقشه توسط ناکاسوجی و ساتاکه (Nakasuji and Satake, 2004) اعلام شده است.

این نقشه‌ها ابزار مهمی برای مراوده ریسک آتشفشانی بین دانشمندان، مسئولان و توده مردم می‌باشند (Haynes et al., 2008) و نقش مهمی در پاسخ اضطراری و کاهش خطر آتشفشان‌ها دارند (Xu et al., 2022). در شرایط بحرانی، می‌توان آن‌ها را از طریق اینترنت به‌سرعت و به­طور وسیع انتشار داد تا در اختیار نیازمندان و افراد ساکن در اطراف آتشفشان قرار گیرند (Thompson et al., 2017). انواع مختلفی از نقشه پهنه­بندی خطر آتشفشان وجود داشته و قابل تهیه است که توسط کالدر و همکاران (Calder et al., 2015) معرفی شده‌اند. هینس و همکاران (Haynes et al., 2008) نقشه‌های پهنه­بندی خطر آتشفشان بوکورون (Boqueron) در کشور سان­سالوادور را براساس سه وضعیت فرضی مختلف (فوران خفیف، متوسط و شدید) تهیه کرده‌اند. به‌عنوان نوعی از نقشه تک خطره می‌توان به نقشه خطر گدازه آتشفشان اتنا (Negro et al., 2013) اشاره کرد. همچنین، نری و همکاران (Neri et al., 2013) با استفاده از تکنیک درخت حوادث، نقشه پهنه­بندی خطر آتشفشان کانلائون4 فیلیپین را تهیه کرده‌اند. همچنین، حتی نقشه پهنه‌بندی خطر قطعات بزرگ پرتابی از دهانه آتشفشان هم تهیه شده‌ است (Alatorre- Ibargueng et al., 2016).

کاپرا و همکاران (Capra et al., 2008) اقدام به پهنه­بندی خطر آتشفشان نوادو د تولوکا 5واقع در 70 کیلومتری غرب مکزیکوسیتی براساس معیارهای زمین­شناسی و شبیه­سازی کامپیوتری با استفاده از نرم‌افزارهای FLOW3D، TTTAN2D، LAHARL وHAZMAP  کرده‌اند. اگرچه این آتشفشان از 33000 سال قبل فورانی نداشته، ولی مشابه با سبلان به لحاظ قرارگیری در یک کمربند فعال و با توجه به فوران‌های گذشته، خطرناک تشخیص داده شده است. براساس این مطالعه، جریان‌های آواری آن می‌توانند تا 15 کیلومتر از دهانه دور شده و خسارت وارد سازند. خاکستر ناشی از فوران نوع پلینی آن تا 70 کیلومتر از محل فوران توانائی وارد ساختن خسارت را دارد. بدین ترتیب، مکزیکوسیتی هم در معرض تهدید آن قرار دارد. لاهارهای آن هم در سال‌های متمادی پس از فوران می‌توانند تولید شوند و امکان خسارت به روستاهای موجود در دره‌های اطراف آتشفشان را خواهند داشت.

سالواتیچی و همکاران (Salvatici et al., 2016) با استفاده از مدل­سازی رقومی اقدام به تهیه نقشه پهنه­بندی خطر جریان‌های آذرآواری در آتشفشان استرومبولی ایتالیا کرده‌اند. همچنین، تاریگان و همکاران (Tarigan et al., 2017) نقشه پهنه‌بندی خطر آتشفشان سینابونگ6  را با استفاده از پهباد تهیه کرده‌اند. چورل و همکاران (Chevrel et al., 2021) نقشه پهنه­بندی خطر آتشفشان سپری جزیره رونیون فرانسه را تهیه کرده‌اند. آن‌ها به مواردی نظیر موقعیت نقطه خروج، زمان تکرار فوران، طول جریان و مسیر جریان توجه داشته­اند.

واگنر و همکاران (Wagner et al., 2015) تعداد 19 نقشه پهنه­بندی متعلق به آتشفشان‌های مختلف دنیا را جمع‌آوری کرده و از نظر نحوه تهیه، نوع خطر و نحوه ارائه مورد مقایسه قرار داده‌اند. ناکاسوجی و ساتاکه (Nakasuji and Satake, 2004) مراحل تهیه نقشه‌های پهنه­بندی را توضیح داده‌اند. تهیه نقشه با استفاده از GIS امکان اصلاح و تغییر آن در هر زمان را فراهم می‌کند. در عین حال که با استفاده از GIS می‌توان نقشه هم‌زمان با فوران را نیز به‌سرعت تهیه کرد و برای کمک به مدیران و افراد درگیر از طریق اینترنت منتشر ساخت.

همان­گونه که ذکر شد، براساس مطالعات خیاط زاده و عباس نژاد ( 1395)، سبلان همراه با تفتان، دماوند و بزمان ازجمله آتشفشان‌های ایرانی مستعد به فوران در آینده است. پیش­تر نقشه‌های پهنه­بندی خطر آتشفشان تفتان (خیاط زاده و همکاران، 1395)، بزمان ( خیاط زاده و همکاران، 1401) و دماوند ( خیاط زاده و همکاران، 1401) تهیه شده‌اند. بنابراین، هنوز چنین نقشه‌هایی برای آتشفشان سبلان تهیه نشده‌اند. اگرچه مطالعات متعددی روی آتشفشان سبلان صورت گرفته است، اما هدف آن‌ها شناخت آن از جنبه آتشفشان‌شناسی است و جنبه‌های خطری آن مورد مطالعه قرار نگرفته‌اند. نبود نقشه پهنه­بندی در هنگام بروز نشانه‌های فعال شدن و در حین فوران می‌تواند خسارات و تلفات ناشی از فوران را به‌شدت افزایش دهد.

 

منطقه مورد مطالعه

استراتوولکان سبلان در 40 کیلومتری باختر اردبیل و در جنوب مشکین‌شهر قرار دارد(شکل 1). ارتفاع آن از سطح دریا 4820 متر است و محدوده‌ای به وسعت حدود 1200 کیلومترمربع را اشغال کرده است. این آتشفشان دارای سه قله است که بلندترین آن‌ها سبلان سلطان نام دارد. یکی دیگر از قله‌ها که هرم داغ نام دارد دهانه‌ای به قطر 200 متر دارد و از بمب، لاپیلی و بلوک پر شده است. این آتشفشان به­طور معمول از آذر ماه به طور کامل توسط برف پوشیده می‌شود.

آتشفشان سبلان روی یک واحد کوهستانی با امتداد خاور- باختر و متشکل از گدازه‌های ائوسن و توده‌های نفوذی قرار دارد (باباخانی، 1387) و دارای یک کالدرای بزرگ با قطر حدود 12 کیلومتر است. این کالدرا توسط فوران‌های بعدی تا حدی پر شده است (شکل 2). پیدایش این کالدرا را ناشی از فرونشست زمین دانسته‌اند. آتشفشان سبلان توسط افراد مختلفی ازجمله دیدون و ژمین (Didon and Gemain, 1976)، امامی (Emami, 1994)، قلمقاش و همکاران (Ghalamghash et al., 2016)، موسوی و همکاران (Mousavi et al., 2014)، امامی (1379)، فتح الهی و همکاران (1392)، باباخانی (1387)، فتح الهی و همکاران (1389)، فتح الهی و خیرخواه (1394) و قلمقاش و همکاران (1398) مورد مطالعه قرار گرفته است. سنگ‌شناسی بعضی از واحدهای ائوسنی که این آتشفشان روی آنها قرار گرفته توسط زمانی و همکاران (1394)، مبشر گرمی و همکاران (1397) و یادگاری و همکاران (1402) مورد مطالعه قرار گرفته است.

 

img0 

شکل 1. موقعیت آتشفشان سبلان و مهم­ترین شهرهای اطراف آن

اولین فعالیت‌های فورانی این آتشفشان از ائوسن شروع شده‌اند، ولی شکل واقعی آن محصول فعالیت در پلیوسن است. البته، در کواترنر هم فعالیت داشته است. بنابراین، در طی یک دوره زمانی به نسبت طولانی شکل گرفته است. فتح الهی و خیرخواه (1394) آن را به‌عنوان عضو جوانی از مجموعه آتشفشانی سنوزوئیک متعلق به کمان ماگمائی البرز ذکر کرده‌اند و معتقدند جوان‌ترین کالدرای آتشفشانی در شرق فلات ایرانی- ترکی است. براساس مطالعات فتح الهی و خیرخواه (1394)، ماگماهای اولیه آن به حاشیه فعال قاره‌ای متاثر از فرورانش اقیانوس تتیس زیر صفحه ایران مرکزی تعلق دارند. باباخانی (1387) فعالیت آتشفشانی آن را هم ارز کواترنر پیشین در نظر گرفته‌ است. قلمقاش و همکاران (Ghalamghash et al., 2016) با استفاده از سن سنجی زیرکن به روش U- Pb زمان فوران سنگ‌های آتشفشانی قدیمی آن را 5/4 تا 3/1 میلیون سال قبل و سن جوان‌ترین سنگ‌های خروجی از آن را 48000 تا 188000 سال قبل به دست آورده است، ولی براساس داده‌های منتشر نشده حسن‌زاده و همکاران (به نقل از فتح الهی و خیرخواه (1394))، جدیدترین فوران سبلان حدود 75000 سال قبل رخ داده است.

در این آتشفشان سه سری فعالیت شناسائی شده است که عبارتند از: (1) سری زیرین، پیش از پیدایش مخروط تشکیل شده و شامل گدازه‌های لاتیتی- آندزیتی میوسن است، (2) سری میانی، مخروط اصلی را ساخته و در طی آن ابتدا گدازه‌های لاتیتی- آندزیتی خارج شده و در پایان جای خود را به داسیت داده‌اند و (3) سری فوقانی، پس از پیدایش کالدرا تشکیل شده ، ترکیب عمده آن آندزیتی- داسیتی است و سن آن پلیوکواترنر است.

در جدول 1 واحدهای زمین‌شناختی تشکیل‌دهنده سبلان به ترتیب سن معرفی شده‌اند. این واحدها حاکی از این می‌باشند که سبلان در فعالیت‌های خود لاهار، گدازه و ابر سوزان (ایگنیمبریت) هم تولید کرده است. بدیهی است در فوران‌های این استراتوولکان خاکستر هم تولید شده و به‌طور عمده در محل­های دور از مخروط نهشته شده‌اند. بنابراین، براساس نوع مواد خروجی، در صورت فوران، می‌تواند خطرهای از نوع گدازه، تفرا (خاکستر)، لاهار و ابر سوزان را ایجاد کند.

جدول 1. واحدهای زمین‌شناسی تشکیل‌دهنده آتشفشان سبلان (با اقتباس از نقشه زمین‌شناسی مشکین‌شهر، Amini, 1994)

دوره فعالیتی 

جنس سنگ‌های تشکیل­دهنده واحدهای مربوطه  

فرسایشی 

شامل:لاهار و پادگانه‌های آبرفتی قدیمی‌تر که خود شامل واریزه‌ها و مورن‌های یخچالی‌اند

مرحله چهارم 

شامل: گدازه‌های آندزیتی- داسیتی و درنهایت گدازه‌های تراکیتی- تراکی آندزیتی

مرحله سوم

شامل: لاهار و آگلومرا، ایگنیمبریت‌های داسیتی، روانه‌های خاکستر داسیتی- ریوداسیتی و گدازه‌ها و گنبدهای ریولیتی- تراکیتی

مرحله دوم

شامل: گدازه و گنبدهای داسیتی- ریوداسیتی در بلندی‌های سبلان

مرحله اول

شامل: جریان‌های گدازه آندزیتی، جریان‌های گدازه تراکی آندزیتی و گدازه‌های آندزیتی پرفیری و داسیت پرفیری

 

دیدون و ژمین(Didon and Gemain, 1976) پنج فاز فورانی در گدازه‌های آن شناسایی کرده‌اند که سه مورد از آن‌ها غالب است و واحدهای عمده را تشکیل داده­اند. این سه مورد عبارتند از: آندزیت‌های زیرین ، تراکی آندزیت‌ها و گدازه‌های داسیتی انتهائی. 

براساس مطالعات باباخانی (1387) و امامی (1379) فرو­­ریزی این آتشفشان و تشکیل کالدرای آن همراه با فعالیت اصلی انفجاری است. در اثر آن، جریان­های آذرآواری، خاکستر، ابر سوزان، برش‌های تیپ پله، روان‌های ایگنیمبریتی و بمب‌های آتشفشانی تولید شده‌اند.

پس از تشکیل کالدرا، جریان‌های کوتاه گدازه خارج شده و گنبدهای آتشفشانی را شکل داده‌اند (شکل 2). آن‌ها بخش‌های مرکزی مرتفع این آتشفشان را اشغال کرده و ترکیبی آندزیتی تا داسیتی دارند. آخرین فعالیت سبلان به‌صورت تشکیل دو مخروط مرکب کوچک و سه جریان گدازه بوده است (امامی، 1379). ناگفته نماند، درحالی­که گدازه‌های این آتشفشان به­طور عمده ترکیب آندزیتی تا داسیتی دارند، ایگنیمبریت‌های آن ترکیب اسیدی‌تر (ریولیتی تا داسیتی) نشان می‌دهند (باباخانی، 1387؛ امامی، 1379 و قلمقاش و همکاران، 1398).

فتح­الهی و خیرخواه (1394) مواد خروجی آن را محصول تبلور تفریقی ماگمای والد بازالتی با خصوصیات کالک آلکن پتاسیم بالا ذکر کرده‌اند و پوسته لیتوسفری غنی شده را منشاء آن دانسته‌اند. این مواد به کمان آتشفشانی حاشیه قاره‌ای تعلق دارند، ولی در پوسته بالایی دچار آلودگی شده‌اند. قلمقاش و همکاران (1398) منشاء آن را ذوب بخشی پوسته قاره‌ای پایینی با ترکیب آمفیبولیت گارنت دار دانسته‌اند.

img0 

شکل 2. برش زمین‌شناختی از آتشفشان سبلان (باباخانی و همکاران، 1367-به نقل از قربانی، 1382)

 

در بعضی از حواشی سبلان (به­طور مثال، نزدیک اردبیل، مشکین‌شهر، سراب و نیر) جریان‌های لاهاری به سن احتمالی پلیستوسن تولید شده‌اند و در سطح مخروط هم به چشم می‌خورند. ضخامت آن‌ها 10 تا 300 متر است. ژئومورفولوژی این آتشفشان توسط اسفندیاری و خیام (1386)، بیاتی خطیبی (1386)، دلال اوغلی (1382 و 1383) و دلال اوغلی و رجبی (1394) مورد مطالعه قرار گرفته است.

بالا بودن دمای ژئوترمال در محدوده این آتشفشان خود مؤید پتانسیل بالای فوران احتمالی آن در آینده می‌باشد. مطالعات زمین گرمائی این آتشفشان توسط افراد مختلف ازجمله رحیمی (1375)، Fotouhi (1995)،  Abdolahzadeh Bina (2009) Ghaedrahmati et al., (2013) انجام شده است.

 

روش مطالعه 

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

 از تصاویر DEM 30 متر به‌عنوان اطلاعات ارتفاعی رقومی استفاده شد و تصاویر رقومی موزائیک تصویری آن تهیه شدند. تصاویر سه بعدی آتشفشان با استفاده از نرم‌افزار ENVI و با داشتن اطلاعات ارتفاعی رقومی به دست آمدند. این تصاویر سه بعدی مبنای تهیه نقشه‌های پهنه‌بندی قرار گرفتند. همچنین، تصاویر لندست ETM+ برای شناسایی کاربری زمین و عناصر تحت خطر مورد استفاده قرار گرفتند.

در تهیه نقشه‌های پهنه­بندی موقعیت نقطه خروج مواد بسیار اهمیت دارد. در این رابطه، همواره نقطه خروج در استراتوولکان‌ها دهانه فعلی در نظر گرته می­شود (Nakasuji and Satake, 2004). بنابراین، در این مطالعه هم نقطه خروج در دهانه فعلی آتشفشان در نظر گرفته شد. یک جنبه دیگر این مطالعه اندازه فوران و مشخصات فوران است و مشابه با سایر استراتوولکان‌های ایران (خیاط زاده و همکاران، 1401 و 1395) عمل شد.

 

نحوه تهیه نقشه پهنه­بندی خطر خاکستر

برای تهیه نقشه گسترش خاکسترهای ناشی از فوران احتمالی این آتشفشان از مدل فرارفت- انتشار (Mastin et al., 2009) استفاده شد. این مدل اغلب برای هشدار خطر خاکستر و مقابله با آن استفاده می­شود (Connor et al., 2001). برای استفاده از این مدل باید پارامترهای تأثیرگذارنده مختلفی؛ نظیر ارتفاع ستون فوران، مقدار مواد فورانی، مدت زمان فوران، اندازه ذرات و سرعت و جهت وزش باد در ارتفاعات مختلف معرفی شوند. از نرم‌افزار VORIS که براساس مدل فوق طراحی شده، استفاده شد. این نرم‌افزار می‌تواند نقشه ضخامت و گسترش خاکستر ناشی از فوران فرضی را نشان دهد .( Felpeto, 2007) اطلاعات جوی مورد نیاز از پایگاه NCEP/NCAR که سرتاسر جهان را پوشش می‌دهد کسب شدند. در این پایگاه، آمار مربوط به‌سرعت و جهت باد در محدوده‌های 5/2 در 5/2 درجه طول و عرض جغرافیائی از سال 1948 ارائه می‌شوند. براساس آن، مقادیر سرعت و جهت وزش باد در ارتفاعات مختلف منطقه سبلان در جدول 2 ارائه شده‌اند. در ضمن، به تبعیت از مرتضوی و همکاران (Mortazavi et al., 2009) و خیاط زاده و همکاران (1395) حجم فوران 5/0 کیلومتر مکعب و اندازه ذرات Md=4/5 و Md= 3/0 میکرون در نظر گرفته شدند.

 

جدول 2. سرعت و جهت وزش باد در سطوح مختلف جو در پهنه سبلان (برگرفته از داده‌های مرکز NCER/NCAR)

 

ارتفاع  (Km)

فشار  (mb)

جهت وزش باد (deg)

سرعت وزش باد (m/s)

1

4820 

500

255

11

2

9000 

300

260

21

3

16000

100

260

18

4

20000

50

251

6.75

5

25000

20

250

6.25

 

نحوه تهیه نقشه پهنه­بندی خطر گدازه 

برای تهیه نقشه پهنه­بندی خطر گدازه از مدل شبیه­سازی جریان گدازه استفاده شد. اساس آن این فرض است که توپوگرافی نقش اصلی را در تعیین مسیر گدازه ایفا می‌کند. در این مدل دو فرض منطقی به کار گرفته می‌شود. اول اینکه گدازه در صورتی از یک نقطه (پیکسل) به نقطه دیگر وارد می‌شود که تفاضل مقدار ارتفاع آن‌ها مثبت باشد. براساس فرض دوم مقدار ورود آن از یک نقطه (پیکسل) به دیگری تابع مقدار عددی تفاضلات مثبت است. در این مدل، با انتخاب یک نقطه به‌عنوان مبداء، مسیرهای حرکت با استفاده از میانگین‌های الگوریتم مونت کارلو7 مشخص می‌شوند. مقادیر ارتفاع براساس مدل ارتفاعی رقومی8 محل تعیین می­شوند. حداکثر طول جریان گدازه در هر مسیر توسط یک پارامتر ثابت برای همه محاسبات کنترل می‌شود   (Felpeto et al., 2007) . البته، با وارد کردن مقادیر انتخابی برای هر پارامتر، پارامترهای کلی (مانند حداکثر طول جریان، تصحیح ارتفاعی و تعداد از سرگیری محاسبات) باید برای هر نقطه خروجی فرضی گدازه به‌صورت عدد وارد شوند. مختصات نقطه شروع فوران به‌صورت دستی با کلیک روی نقطه مورد نظر مشخص می‌شوند. در این مطالعه، حداکثر طول جریان 5 کیلومتر در نظر گرفته شد. همچنین، ارتفاع روانه برای بی‌اثر کردن موانع موجود در مسیر و یا خطای مدل ارتفاعی رقومی سه متر فرض شد. در ضمن، به‌عنوان شروع حرکت گدازه 5000 نقطه در اطراف دهانه انتخاب شدند.

 

نحوه تهیه نقشه پهنه­بندی خطر پیروکلاستیک‌های جریانی

در این مطالعه، برای تعیین حداکثر مسافت گسترش پیروکلاستیک‌های جریانی از مدل ارائه شده توسط مالین و شریدان (Malin and Sheridan, 1982) استفاده شد. اساس این مدل در نظر گرفتن ارتفاع نقطه شروع جریان (Hc) و نسبت آن با طول مسافت جریان (L) به‌عنوان یک پارامتر اصطکاکی است و در آن به زاویه میل جریان (cα) نیاز است. این زاویه از طریق نسبت  تعیین می‌شود. پس، در آن، محدوده حرکت با محل تلاقی جریان با سطح زمین (Felpeto et al., 2007) برابر است. این مدل به‌خوبی توسط فلپتو (Felpeto, 2009) توضیح داده شده است. در عمل، برای اجرای آن مدل ارتفاعی رقومی این آتشفشان مورد استفاده قرار گرفت و با توجه به موقعیت دهانه، مختصات نقطه شروع مشخص شد. پس از آن، مقادیر پارامترهای مورد نیاز در پنجره شبیه‌سازی آن به‌صورت دستی وارد شدند. با فرمان اجرا9، محاسبات انجام شده و لایه اطلاعاتی جدیدی به دست آمد و با استفاده از نرم‌افزار Arc GIS به نقشه پهنه­بندی پیروکلاستیک‌های جریانی تبدیل شد. در ضمن، در این مطالعه ارتفاع سقوط ابر جریانی 200 متر و زاویه معادل سقوط (زاویه بین ارتفاع سقوط و شیب توپوگرافی)  شش درجه در نظر گرفته شدند.

 

نحوه تهیه نقشه پهنه­بندی خطر لاهار

برای تهیه این نقشه، با استفاده از تصویر سه بعدی تهیه شده از تلفیق مدل ارتفاعی رقومی آتشفشان و تصویرهای ماهواره‌ای موزائیک شده آن، آبراهه‌های اصلی مشخص و بافرگذاری شدند. آن‌ها به‌عنوان معابر اصلی حرکت لاهار عمل می‌کنند. به‌عنوان محدوده خطر لاهار، بافر 500 متر برای آبراهه‌های اصلی در نظر گرفته شد و نقشه مسیرهای احتمالی جریان‌های لاهاری به دست آمد. ناگفته نماند، در این روش مسیر لاهار مشخص می­شود، ولی مسافت طی شده تابع شرایط زمان و مکان فوران است.

بحث

احتمال فعالیت 

سبلان در دوره تاریخی فاقد هرگونه فوران بوده است. چشمه‌های آب­گرم متعدد اطراف این مخروط آتشفشانی که توسط شاه بیک (1372) و کمالی نژاد و عباس نژاد (1398) مطالعه و معرفی شده‌اند؛ گزارش خروج گازهای گوگردی در برهه‌های زمانی کوتاه توسط افراد محلی؛ فعالیت شدید تکتونیکی پهنه (رحیمی و همکاران، 1375) و نیز، به نسبت سرپا بودن آن و فعالیتش در دوران چهارم زمین‌شناسی، همگی مؤید نیمه فعال بودن آن و احتمال فوران مجدد آن می‌باشند. براساس مطالعات خیاط زاده و عباس نژاد (1395)، این آتشفشان در میان استراتوولکان‌های ایران دارای پتانسیل فوران در آینده است.

 

اندازه فوران احتمالی

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

 

نوع خطرهای احتمالی

همان­گونه که ذکر شد، سبلان در گذشته با ایجاد فعالیت‌های از نوع خروج گدازه، تولید ابر سوزان، فوران خاکستر و ایجاد لاهار همراه بوده است. توجه به استراتوولکان بودن آن و نیز، نوع تهدیدهایی که در گذشته داشته حاکی از این است که خطرهای فوق (گدازه، لاهار، ابر سوزان و خاکستر) در فوران‌های بعدی آن باید جدی گرفته شوند.

 

عناصر تحت خطر

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

 

محدوده‌های تحت خطر تفرا (خاکستر)

نقشه گسترش خاکسترهای آتشفشانی ناشی از فوران احتمالی سبلان در شکل 3 نشان داده شده است. با توجه به این نقشه، تحت تأثیر بادهای غالب مرتفع در این عرض جغرافیائی که غربی- شرقی می‌باشند خاکسترها به سمت شرق تا اندکی متمایل به سمت شمال‌شرق حرکت کرده و درنهایت در سطح زمین نهشته می‌شوند. شهر اردبیل و روستاهای شرق مخروط در محدوده گسترش پلوم آن قرار داشته و می‌توانند دچار مشکل شوند. براساس فرضیات انجام شده برای تعیین پلوم خاکستر، بستگی به فاصله و موقعیت محل، ضخامت خاکستر می‌تواند در مواردی به حدود سه متر هم برسد. در اثر خاکستر، مشکلاتی از نظر تنفس، دید در هنگام رانندگی، عبور هواپیماها، آلودگی منابع آب و آسیب به محصولات کشاورزی عارض می‌شوند. در جدول 3 اثرات ناشی از ضخامت‌های مختلف خاکستر آتشفشانی درج شده است.

جدول 3. اثرات ضخامت‌های مختلف خاکستر (Michaud-Dubuy et al., 2021- با تغییرات)

ردیف

ضخامت خاکستر

اثرات

1

در حد میلی‌متر

توقف فعالیت فرودگاه، ایجاد مشکلات برای شبکه‌های تأمین

2

در حد سانتیمتر

آسیب به فرودگاه، آسیب به پوشش گیاهی، ایجاد مشکلات برای شبکه‌های تأمین (جاده و غیره)

3

حدود 20 سانتیمتر

غیرقابل عبور شدن جاده‌ها، ریزش سقف خانه‌های چوبی

4

بیش از 1 متر

تخریب کامل زیرساخت‌ها

 

C:\Users\user\Desktop\سبلان\Sabalan\tephra.jpg 

شکل 3. نقشه گسترش خاکسترهای آتشفشانی در اثر فوران احتمالی سبلان

 

 

محدوده‌های تحت خطر جریان گدازه

براساس فرضیات و روش ارائه شده در مبحث «مواد و روش‌ها»، نقشه پهنه­بندی خطر گدازه ناشی از فوران احتمالی سبلان در شکل 4 ارائه شده است. با توجه به این نقشه، گدازه‌ها در مجاری حفر شده در شمال و جنوب مخروط حرکت کرده و تنها جریان شمالی می‌تواند به مشکین­شهر نزدیک شود. براساس این نقشه، اکثر جان­پناه‌ها، تأسیسات استفاده از آب­گرم و پیست اسکی آلوارس در مسیر گدازه‌های احتمالی قرار خواهند داشت، ولی احتمال رسیدن گدازه‌ها به شهرها و روستاهای اطراف مخروط کم است.

C:\Users\user\Desktop\سبلان\Sabalan\lava flow.jpg شکل 4. نقشه گسترش گدازه‌های خروجی از دهانه اصلی آتشفشان سبلان

 

محدوده‌های تحت خطر پیروکلاستیک‌های جریانی

نقشه گسترش ابر سوزان و سرج در شکل 5 نشان داده شده است. پیش­تر (در مبحث مواد و روش‌ها) فرضیه­های مرتبط با تهیه این نقشه توضیح داده شدند. براساس این نقشه، مشکین‌شهر، روستای موئیل، محل‌های استقرار عشایر شاهسون و قطورسوئی (در شمال مخروط) و هر دو پیست اسکی اوج و آلوارس، روستای آلوارس در جنوب، اکثر جان­پناه‌ها و آب­گرم قطورسوئی در معرض تهدید قرار دارند. به عقیده قلمقاش و همکاران (1398) آتشفشان سبلان استعداد فوران‌های دمای بالا، سریع و پر گسترش ابر سوزان را دارد.

C:\Users\user\Desktop\سبلان\Sabalan\PDCs.jpg 

شکل 5. نقشه گسترش پیروکلاستیک های جریانی (نویی آردنت) در اثر فوران احتمالی سبلان

 

محدوده‌های تحت خطر لاهار

نقشه مسیرهای احتمالی و گسترش لاهار در شکل 6 نشان داده شده است. این نقشه مسیرهای احتمالی لاهار را نشان می‌دهد، ولی مسافت طی شده توسط آن‌ها متأثر از حجم آب موجود در مخروط (ناشی از ذوب یخ و برف و بارش باران) است و در ایام مختلف سال می‌تواند تفاوت زیادی داشته باشد. بنابراین، ابعاد لاهارها تابع اوقات فوران در سال هم هست، به­طوری­که در زمستان به علت تراکم بیشتر برف لاهارهای بیشتری می‌توانند تولید شوند و تا مسافت طولانی‌تری پیش می‌روند.

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

 

 

C:\Users\user\Downloads\2_144176761431115765.jpg 

 شکل 6. نقشه مسیر لاهارهای ایجاد شده از فوران فرضی سبلان

 

نتیجه‌‌گیری

براساس آنچه در مورد شرایط آتشفشان سبلان توضیح داده شد، احتمال فوران آن در آینده مطرح است. در این صورت، ممکن است فوران با شدت سه تا پنج در مقیاس VEI صورت گیرد. در اثر آن، تهدیدهای ناشی از خاکستر، گدازه، جریان‌های آذرآواری و لاهار امکان بروز دارند. این تهدیدها می‌توانند به کاربری‌های سطح مخروط و حاشیه آن آسیب برسانند (به‌ویژه لاهارها و پیروکلاستیک‌های جریانی). خاکسترهای این آتشفشان تحت تأثیر جریان بادهای غالب مرتفع به سمت شرق منتشر شده و رسوب می‌کنند. آن­ها می‌توانند مشکلاتی برای شهر اردبیل و روستاهای اطراف ایجاد کنند.

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

در این مطالعه، نقشه‌های پهنه­بندی براساس محتمل‌ترین وضعیت برای هر یک از خطرها تهیه شده‌اند، ولی بهتر است در مطالعات بعدی نقشه‌های پهنه­بندی هر یک از خطرها برای وضعیت­های احتمالی مختلف هم تهیه شوند. البته، از آنجا که راهی برای صحت­سنجی این نوع نقشه‌ها وجود ندارد، شاید بهتر باشد انواعی از نقشه با روش‌ها و تکنیک‌های مختلف تهیه شوند. در این حالت، وضعیت‌های مورد تائید نقشه‌های متعدد، قابل قبول‌تر و صحیح‌تر فرض می‌شوند.

 

منابع

اسفندیاری درآباد، ف. و خیام، م.، 1386. تحلیلی بر اثرات ژئومورفیک برفاب در دامنه شرقی سبلان، پژوهش‌های جغرافیائی، 6، 83-97. ##امامی، م.، 1379. ماگماتیسم در ایران، انتشارات سازمان زمین‌شناسی و اکتشافات معدنی کشور، 608. ##باباخانی، ع.، 1387. کوه آتشفشان سبلان، گزارش سازمان زمین‌شناسی کشور، مرکز تبریز، 122. ##بیاتی خطیبی، م.، 1386. تحلیل و بررسی نقش عوامل توپوگرافی و دینامیک رودخانه‌ای بر اندازه مخروط واریزه‌ای؛ مطالعه موردی: دامنه شمال غربی سبلان (شمال غرب ایران)، پژوهش‌های جغرافیائی، 6، 157- 175. ##خیاط زاده، ا. و عباس نژاد، ا.، 1395. کاربرد روش‌های تصمیم‌گیری EN- SAW و ANP در مطالعات زمین‌شناسی: ارزیابی و اولویت‌بندی خطر فعال شدن استراتوولکان‌های ایران به‌عنوان مطالعه موردی، فصلنامه علوم زمین، (102)26، 137- 146. ##خیاط زاده، ا.، عباس نژاد، ا. و رنجبر، ح.، 1395. پهنه‌بندی خطر تفرا، گدازه و ابرهای سوزان احتمالی آتشفشان تفتان، فصلنامه جغرافیا و مطالعات محیطی، (18)5، 17- 30. ##خیاط زاده، ا.، عباس نژاد، ا.، رنجبر، ح.، احمدی پور، ح. و عباس نژاد، ب.، 1402. ارزیابی و پهنه‌بندی خطر تفرا، گدازه و ابرهای سوزان ناشی از فوران احتمالی آتشفشان بزمان. زمین‌شناسی کاربردی پیشرفته, 13(2)، 455-468.‎ ##خیاط زاده، ا.، عباس نژاد، ا.، رنجبر، ح.، احمدی پور، ح. و عباس نژاد، ب.، 1402. پهنه‌بندی خطر تفرا، گدازه، ابر سوزان و لاهار ناشی از فوران احتمالی آتشفشان دماوند. فصلنامه علمی علوم زمین, 32(4)، 355-366.‎ ##دلال اوغلی، ع. و رجبی، م.، 1394، بررسی اختلافات ساختاری و فرسایشی در مخروط‌های استراتوولکانوسهند و سبلان، اولین همایش ملی علوم زمین و توسعه شهری. ##دلال اوغلی، ع.، 1382. بررسی مورفولوژی و نحوه فعالیت یخچال‌های سنگی دامنه شمالی کوه سبلان، پژوهش‌های جغرافیائی، 45، 1- 12. ##رحیمی، ع.، 1375. زمین‌شناسی میدان زمین گرمائی شمال غرب سبلان (جنوب مشکین‌شهر) و ژئوشیمی دگرسانی و ته نشست‌های گرمابی منطقه، پایان‌نامه کارشناسی ارشد دانشگاه تهران، 284. ##زماني، ر.، امامي، م.، وثوقي عابديني،م. و كريم زاده ثمرين، ع.، 1394. مطالعه شيمي بلور و منطقه بندي کلينوپيروکسن هاي موجود در سنگ‌های آتشفشاني آلکالن شمال غرب مشکین‌شهر، ايران، فصلنامه زمین‌شناسی ایران، (33)9، 31-44. ##شاه بیک، ا.، 1372. آب‌های معدنی و گرم ایران، سازمان زمین‌شناسی کشور، 212. ##فتح الهی، م. و خیرخواه، م.، 1394. منشاء و جایگاه تکتونو ماگمائی سنگ‌های آتشفشانی کواترنری سبلان، فصلنامه کواترنری ایران، جلد 1، 2، 125- 136. ##فتح الهی، م.، امامی، م. و خیرخواه، م.،1392. پتروگرافی و پترولوژی سنگ‌های آتشفشانی ماگماتیسم نهایی سبلان، فصلنامه زمین‌شناسی محیط‌زیست (7)23، 51-68. ##فتح الهی، م.، امامی، م. و خیرخواه، م.، 1389. تاثیر فرایند اختلاط ماگمائی و آلایش پوسته‌ای در پیدایش سنگ‌های آتشفشانی سبلان، بیست و نهمین گردهمایی علوم زمین.##قربانی، م.، 1382. آتشفشان شناسی با نگرشی بر آتشفشان‌های ایران، نشر آرین زمین، 362. ##قلمقاش، ج.، کتابی، ز.، اصفهانی، ا. و موسوی، ز.، 1398. مقایسه زمین‌شناسی، پتروگرافی، زمین‌شیمی و منشاء ایگنیمبریت در آتشفشان‌های سبلان و سهند، زمین‌شناسی کاربردی پیشرفته، (131)9، 17- 24. ##کمالی نژاد، م. و عباس نژاد، ا.، 1398. مطالعه چشمه‌های آبگرم و معدنی ایران (ویژگی‌های شیمیائی، درمانی و گردشگری)، انتشارات سنجش و دانش، 435. ##مبشر گرمی، م.، احمدزاده، غ. و آقازاده، ا.، 1397، برآورد ویژگی‌های فیزیکوشیمیایي و محیط تکتونیکی تشکیل توده‌های بازیک جنوب شهرستان گرمی (استان اردبيل) بر اساس شيمي بلورهای آمفیبول و پلاژیوکلاز، )، فصلنامه زمین‌شناسی ایران، (12) 48، 81-98. ##یادگاری، ن.، علوی، غ. و مؤید، م.، 1402، بررسی زمین‌ساخت، داده‌های گسلی و ارتباط آن‌ها با کانه‌زایی و دگرسانی در گستره اسبخان هریس (استان آذربایجان‌شرقی - شمال‌غرب ایران)، فصلنامه زمین‌شناسی ایران، (17)65، 33-47. ##Abdollahzadeh Bina, F., 2009. Geothermal resource assessment of the NW- Sabalan Field, Iran, Through well testing, Geothermal training programe, United Nations University (UNU), 6, 15- 44. ##Alatorre-Ibargüengoitia, M.A., Morales-Iglesias, H., Ramos-Hernández, S.G., Jon-Selvas, J. and Jiménez-Aguilar, J.M., 2016. Hazard zoning for volcanic ballistic impacts at El Chichón Volcano (Mexico). Natural Hazards, 81(3), 1733-1744. ##Amini, B., 1994. Geological map of Meshkinshahr, Geological Survey of Iran. ##Blong, R. J., 1984. Volcanic hazards: a sourcebook on the effects of eruptions, Elsevier. ##Calder, E., Wagner, K. and Ogburn, S., 2015. Volcanic hazard maps. Global volcanic hazards and risk, 335-342. ##Capra, L., Norini, G., Groppelli, G., Macías, J.L. and Arce, J.L., 2008. Volcanic hazard zonation of the Nevado de Toluca volcano, México. Journal of Volcanology and Geothermal Research, 176(4), 469-484. ##Charlton, D., 2018. New approaches to volcanic hazard mapping at Campi Flegrei, Southern Italy (Doctoral dissertation, UCL (University College London)). ##Chevrel, M.O., Favalli, M., Villeneuve, N., Harris, A.J., Fornaciai, A., Richter, N., Derrien, A., Boissier, P., Di Muro, A. and Peltier, A., 2021. Lava flow hazard map of Piton de la Fournaise volcano. Natural Hazards and Earth System Sciences, 21(8), 2355-2377. ##Connor, C. B., Hill, B. E., Winfred, B., Franklin, N. W. and La Femina, P. C., 2001. Estimation of volcanic hazards from tephra fallout, Natural Hazards Review, 2, 33- 42. ##Didon, J. and Gemain, Y. A., 1976. La Sabalan, volcan plio- Quaternaire de l, Azerbaijan oriental Iran, Etude Geologique et Scientifique et Medical de Grenoble, France. ##Emami, M., 1994. Geological Map of Iran, Meshkinshahr quadrangle (1: 100000), Geological Survey of Iran. ##Felpeto, A., 2009. VORIS, a GIS based tool for volcanic hazard assessment, User Guide, Observatorio Geofisico Central, IGN. ##Felpeto, A., Marti, J. and Ortiz, R., 2007. Automatic GIS- based system for volcanic hazard assessment, Journal of volcanology and Geothermal Research, 166, 106- 116. ##Fotouhi, M., 1995. Geothermal development in Sabalan, Iran, in: proceedings of world geothermal congress, Italy, 191- 196. ##Ghaedrahmati, R., Moradzadeh, A., Fathianpour, N., Lee, S. K. and Porkhial, S., 2013. 3D inversion of MT data from Sabalan geothermal field, Ardabil, Iran, Journal of Applied Geophysics, 93, 12- 24. ##Ghalamghash, J., Mousavi, S. Z., Hassanzadeh, J. and Schmitt, A. K., 2016. Geology, zircon geology, and petrogenesis of Sabalan volcano (northwestern Iran), Journal of Volcanology and Geothermal Research, 327, 192- 207. ##Haynes, K., Barclay, J. and Pidgeon, N., 2008. The issue of trust and its influence on risk communication during a volcanic crisis. Bulletin of Volcanology, 70, 605-621. ##Latter, J. H. (ED.), 2013. Volcanic hazards: assessment and monitoring, Springer Science and Business Media. ##Loughlin, S. C., Sparks, R. S. J., Sparks, S., Brown, S. O K., Jenkins, S. F. and Vye- Brown, C., (Eds). 2015. Global volcanic Hazards and risk, Cambridge University Press. ##Malin, M. C. and Sheridan, M. F., 1982. Computer- assisted mapping of pyroclastic surges, Science, 217 (4560), 637- 640. ##Mastin, L. G., Guffanti, M., Servrankx, R., Webley, P., Barsotti, S. and Dean, K. et al., 2009. A multidisuplinary efforf to assign realistic source parameters to models of volcanic ash- cloud transport and dispersion during eruptions, Journal of volcanology and Geothermal Research, 186, 10- 21. ##Mc Guire, W. J., 1998. volcanic hazards and their mitigation, Geological Society, London, Engineering Geology Special Publications, 15 (1), 79- 95. ##Michaud-Dubuy, A., Carazzo, G. and Kaminski, E., 2021. Volcanic hazard assessment for tephra fallout in Martinique. Journal of Applied Volcanology, 10(1), 1-20. ##Mousavi, S. Z., Darvishzadeh, A., Ghalamghash, G. and Vosoughi Abedini, M., 2014. Volcanology and geochronology of Sabalan volcano, the hieghest stratovolcano in Azerbaijan region, NW Iran, Nautilus, 128, 85- 98. ##Nakasuji, A. and Satake, j., 2004. Volcanic Hazard Map An Introduction and Overseas Cases. Journal of the Japan Society of Engineering Geology, 44(6), 341-348. ##Negro, C.D., Cappello, A., Neri, M., Bilotta, G., Hérault, A. and Ganci, G., 2013. Lava flow hazards at Mount Etna: constraints imposed by eruptive history and numerical simulations. Scientific Reports, 3(1), 1-8. ##Neri, M., Le Cozannet, G., Thierry, P., Bignami, C. and Ruch, J., 2013. A method for multi-hazard mapping in poorly known volcanic areas: an example from Kanlaon (Philippines). Natural Hazards and Earth System Sciences, 13(8), 1929-1943. ##Salvatici, T., Di Roberto, A., Di Traglia, F., Bisson, M., Morelli, S., Fidolini, F., Bertagnini, A., Pompilio, M., Hungr, O. and Casagli, N., 2016. From hot rocks to glowing avalanches: Numerical modelling of gravity-induced pyroclastic density currents and hazard maps at the Stromboli volcano (Italy). Geomorphology, 273, 93-106. ##Tarigan, A.P.M., Suwardhi, D., Fajri, M.N. and Fahmi, F., 2017. March. Mapping a volcano hazard area of Mount Sinabung using drone: preliminary results. In IOP Conference Series: Materials Science and Engineering 180(1), 012277. ##Thompson, M. A., Lindsay, J. M. and Leonard, G. S., 2018. More than meets the eye: volcanic hazard map design and visual communication. Observing the volcano world: volcano crisis communication, 621-640. ##Tilling, R. J., 1984. Volcanic hazards: a sourcebook on the effects of eruptions, Elsevier. ##Wagner, K., Ogburn, S.E. and Calder, E.S., 2015. Volcanic Hazard Map Survey. ##Walker, G. P., 1982. Volcanic hazards, Interdisciplinary Science Reviews, 7 (2), 148- 157. ##Xu, J., Wa, Y., Wa X., Pan, B., Ya, H., Zh. B. and Ya W., 2022. Review on the development of volcanic hazard zonation in china. Geology and Resources, 1;31(3).##

 


 

 

Assessment and zoning of hazards related to likely eruption of Sabalan volcano

 

Ahmad Abbasnejad1*, Ahmad Khaiatzadeh2, Hojjatolah Ranjbar3, Hamid Ahmadipour4, Behnam Abbasnejad5

1. Professor, Department of Geology, Faculty of Sciences, Shahid Bahonar University of Kerman

2. M.Sc. graduate, Department of Geology, Faculty of Sciences, Shahid Bahonar University of Kerman

3. Professor, Department of Mining Engineering, Faculty of Engineering, Shahid Bahonar University of Kerman

4. Professor, Department of Geology, Faculty of Sciences, Shahid Bahonar University of Kerman

1. Post Doc Researcher, Department of Geology, Faculty of Sciences, Shahid Bahonar University of Kerman

 

 

 Abstract

Sabalan is one of dormant stratovolcanoes of Iran with likely eruption, because there are many thermal springs around it. It has a slightly- eroded cone and, geologically, a short time has passed from its latest eruption. According to eruption regime of such volcanoes, as well as its former activities, Sabalan volcanic eruptions may happen along with such hazards as tephra, lava, pyroclastic flows and lahars. Its eruption intensity on VEI scale is anticipated to be 3 to 5. Whenever this volcano shows awakening evidence, it would be necessary to take appropriate actions according to hazard zoning maps which are the subject of this study. That is, these maps are necessary for taking emergency actions. In this study, digital elevation data (DEM); satellite pictures; Arc GIS, ENVI and VORIS softwares; as well as atmospheric data of NCEP/NCAR center were used. For construction of pyroclastic density current spreads, the Malin and Sheridan (1982) model was employed. For preparation of lava flow hazard map, simulation model was used, and for preparation of lahar hazard zones, major valleys on the cone and its surroundings were identified and buffered using DEM and satellite data. According to this study, volcanic ashes will move eastward and threat Ardabil city and several nearby villages. Lahars would threaten many surrounding land uses and probably Meshkinshahr and Ardabil cities. Lava flows would damage many mountaineering shelters and ski facilities of Alvares. Additionally, Meshkinshahr, Moil, ski resorts as well as many shelters are at the threat of nuee ardentes.

 

 Keywords: Sabalan, Lahar hazard, Lava hazard, Nuee ardente hazard, Tephra hazard.

 

 

17

 

[1] * نویسنده مرتبط:  Aabbas@uk.ac.ir

[2]  Dormant

[3]  Communication

[4]  Kanlaon

[5]  Nevado de Toluca

[6]  Sinabung

[7]  Mont Carlo Method

[8]  DEM

[9]  RUN


رایمگ

سامانه رایمگ تمامی فرآیندهای دریافت، ارزیابی و داوری، ویراستاری، صفحه‌آرایی و انتشار الکترونیکی نشریات علمی را به انجام می‌رساند

پیوندهای سایت

مراکز مرتبط

پشتیبانی

صفحات رسمی

حقوق این وب‌سایت متعلق به سامانه مدیریت نشریات رایمگ است.
حق نشر © 1403-1396

صفحه اصلی| عضویت/ ورود| درباره رایمگ| تماس با ما|
[English] [العربية] [en] [ar]