Abstract:
آتشفشان تفتان با ارتفاع حدود ۴۰۰۰ متر در ۴۵ کیلومتری شمال شهر خاش در استان سیستان و بلوچستان واقع شده است. این استراتوولکان یکی از آتشفشان های متعدد تشکیل دهنده قوس آتشفشانی زون فرورانشی مکران است که در آن پوسته اقیانوسی دریای عمان به سمت شمال دچار فرورانش می شود. مطالعات قبلی انجام شده در مورد آن حاکی از خروج گدازه، فوران خاکستر و تولید جریان های پیروکلاستیکی (ابر سوزان) در گذشته بوده است. با توجه به مورفولوژی نه چندان فرسوده این آتشفشان، وجود چشمه های آبگرم و سولفاتار و نیز انجام فرآیند فرورانش در حال حاضر، احتمال فوران مجدد این آتشفشان وجود دارد که در این صورت ممکن است فوران از نوع استرومبولی، ولکانو تا ساب پلینی باشد. لذا به منظور تهیه نقشه های پهنه بندی خطر تفرا، گدازه و ابرهای سوزان این آتشفشان از مدل های ارتفاعی رقومی (DEM)، تصاویر ماهواره ای و نرم افزار های ArcGIS، ENVI و VORIS استفاده گردید. همچنین برای تهیه نقشه خاکستر مدل فرارفت انتشار و برای تهیه نقشه خطر پیروکلاستیک های جریانی مدل مالین و شریدان (۱۹۸۲) به کار گرفته شدند. اطلاعات جوی مورد نیاز از داده های پایگاه مرکز NCEP/NCAR اخذ گردیدند. درنهایت، نقشه های پهنه بندی خطر تهیه گردیدند. بر اساس این مطالعه، خاکسترهای ناشی از فوران به سمت شرق آتشفشان حرکت کرده و تعدادی از روستاهای این محدوده را تهدید خواهند کرد ولی گدازه های خروجی به سمت شمال و جنوب جریان می یابند و به روستاهای تمین، سنگان، گوشه و تمندان می توانند برسند. در صورت فوران از نوع ابر سوزان، آن ها در بعضی سطوح حاشیه ای شمال، شرق و جنوب آتشفشان پخش خواهند شد و روستاهای این محدوده ها را تهدید خواهند کرد. درمجموع نقشه-های تهیه شده محدوده های تحت خطرهای آتشفشانی را مشخص می سازند و برای اقدامات مدیریتی نظیر تغییر کاربری زمین، پیش بینی و هشدار، آمادگی برای مقابله و امداد و نجات می توانند مورد استفاده قرار گیرند.
Introduction : There are about 1500 active volcanoes across the world، among which only about 400 cases have erupted during last century، only one third of which have had recorded historical eruptions. Also، over 10% of world population lives around active and semi-active volcanoes (Tilling، 2005). The majority of hazardous volcanoes are of stratovolcano type. Plio-quaternary stratovolcanoes of Iran are considered as potential volcanic threats in Iran. Hence، it is helpful to evaluate the probability of their reactivation، their kind of threat as well as to determine locations at threat.
Study area: In this regard، we decided to evaluate the hazard aspects of Taftan volcano. Taftan volcano which is 4000 m in height lies at 61°،8ʹ eastern longitude and 28°،36ʹ northern latitude at 45 km north of Khash town in Sistan and Baluchistan province of Iran (Bahrami، Yamani & Alavipanah، 2008). This stratovolcano belongs to the volcanic arc of Makran subduction zone in which the oceanic lithosphere of Oman sea is subducted northward beneath continental lithosphere of eastern Iran and Pakistan (Farhoudi & Karig، 1977). Currently، this volcano exhibits post volcanic activity via sulfataric and fumarolic activities near its summit (Moeinvaziri & Aminsobhani، 1978). Although previously several geological studies have been undertaken on Taftan، but their purpose have been to understand it، so they have not dealt with its hazards and hazard zoning.Since the potassium-Argon and rubidium-strontium dating investigations undertaken by Biabangard (2006) revealed that its latest flows have erupted about 6.95±0.02 thousand years ago، the conclusion is that we deal with a potentially active volcano. Taftan has erupted with strombolian to plinian types of activity in the past (Biabangard & moradian، 2002) and it may erupt with 3-5 intensity on VEI scale in the future. Moreover، it is probable that the eruptions have pyroclastic flows (nuée ardente) or ash (Tephra). Thus، hazard zoning must be conducted for this kind of threat.
Material and methods: After gathering the needed data، Landsat images and digital elevation model (DEM) of the study area were employed to create a 3D model and hill-shaded topographic model to enhance visibility. They also were used to prepare the land use layers of the area. In this study، we employed the VORIS software which is based on advection-diffusion model. In this model several parameters such as height of eruption column (H)، the volume of erupting materials (M)، eruption duration (D)، particle size، velocity and wind direction at different elevations of the atmosphere، air temperature and some other are used. Relevant parameters are needed and were determined for the studied area was introduced to the software and it was run. The output was the ash thickness map (isopach map) (Mastin، 2009).
The contributing factors on tephra movement the air include horizontal velocity of the wind، turbulence and the rate of vertical discendance which vary according the Reynolds number of particles. The atmosphere is divided into several horizontal layers each of them having specific velocity and directions. Ash particles follow the conditions of each atmospheric layer and when entering a new layer are affected by the conditions of this new one. This process continues until the ash particle reaches the ground surface (Mortazavi، Sparks، & Amigo، 2009). The necessary wind data include its direction and velocity in five atmospheric levels (4000، 9000، 16000، 20000 & 35000 m) which in this study were taken from NCEP/NCAR national centers.
In order to determine the route of lava flows from a point (the summit of volcano)، the simulation model of lava flow was employed. This model is based on this presumption that topography plays the main role in the route of flows. This model uses two logics، the first one supposes that lava will flow from one cell (pixel) to the nearby pixel whenever their elevation differences are positive، and secondly، the volume of lava entering from one pixel to the other one depends on the value of positive differences between them. In preparing the zoning map with this model، the maximum distance for moving lava was considered to be 5 kms.
The model used for simulating the maximum potential extent affected by pyroclastic density currents (PDC) is a very simple model proposed by Malin and Sheridan (1982). The principle is that the height of the starting point of the flow (Hc) ratios to the length of the run out (L) as a type of friction parameter is termed the Heim coefficient. The inclination of the energy cone is an angle (αc) defined by arctan (Hc/L). The intersection of the energy cone، originating at the eruptive source، with the ground surface defines the distal limits of the flow (Felpeto، 2009). In this study، the collapse equivalent height of nuée ardente was considered to be 200 m and the collapse equivalent angel which is the angel between collapse height and slope angel was presumed to be 6°.
Conclusion: According to this study، Taftan is a potentially semi-active volcano which may reactive in the future. It has erupted lava flows، tephras and nuée ardentes in the past. According to the presented hazard zonation map، the flows and nuée ardentes threats some nearby villages. The ashes will disperse and extend eastward and will rise and accumulate in some villages. It is suggested that a seismological station be established near this volcano for forecasting its probable eruptions. The kind of threats identified and the hazard zoning maps may be used for emergency evacuation and hazard reduction programs.
Machine summary:
"پهنه بندی خطر تفرا، گدازه و ابرهای سوزان ناشی از فوران احتمالی آتشفشان تفتان احمد خیاط زاده ١- کارشناس ارشد زمین شناسی زیست محیطی- دانشگاه شهید باهنر کرمان – ایران احمد عباس نژاد – دانشیار ژئومورفولوژی – دانشگاه شهید باهنر کرمان – ایران حجت الله رنجبر – استاد- گرایش دورسنجی – دانشگاه شهید باهنر کرمان - ایران تاریخ دریافت : ١٣٩٤/٩/١٥ تاریخ تصویب : ١٣٩٥/٣/٣١ چکیده آتشفشان تفتان با ارتفاع حدود ٤٠٠٠ متر در ٤٥ کیلومتری شمال شهر خاش در استان سیستان و بلوچستان واقع شده است .
در مجموع ، اگرچه بشر در طی زمان تاریخی ناظر بر فوران آن نبوده ، ولی با توجه به فعالیت شدید فومرولیها، سولفاتارها و چشمه های آبگرم و سن آخرین فوران و عدم فرسایش شدید این آتشفشان ، و نیز، با توجه به اینکه هم اکنون نیز فرآیند فرورانش صفحه 1 Hall 2 Keller & DeVecchio 3 Tilling اقیانوسی عمان که عامل ایجاد این آتشفشان است ادامه دارد، میتوان نتیجه گرفت که احتمال فوران مجدد آن وجود دارد و بدین لحاظ لازم است نوع تهدیدهای ناشی از آن و مناطق تحت تهدید مورد شناسایی و ارزیابی واقع شوند.
) شکل ٧ پنجره پارامترهای ورودی مدل جریان های پیروکلاستیک برای وضعیتی که فوران از یک دهانه انتخاب شده پهنه بندی خطر الف ) خطر تفرا (خاکستر آتشفشانی) همانگونه که ذکر شد خطر تولید و پخش خاکستر در اثر فوران احتمالی تفتان وجود دارد."