پهنه بندی گستره ایران بر اساس تغییرات نسبت ضرایب لرزه خیری a/b

ضرایب لرزه خیزیb, a در رابطه گوتنبرگ ریشتر از مهمترین پارامتر ها در مطالعه الگوی لرزه خیزی یک منطقه می باشد. تهیه نقشه پهنه بندی ضرایب لرزه خیزی bvalue و a value و همچنین نسبت آن ها a/b value برای گستره ایران هدف اصلی در این پژوهش است. ابتدا یک کاتالوگ از زمین لرزه های رخ داده از سال های 1900 تا 2007 در گستره ایران تهیه شده است، سپس کل ایران به پهنه هایی بسیار کوچک تقسیم شده و ضرایب لرزه خیزی b, a برای هر سلول از این شبکه به طور جداگانه محاسبه شده است. بر اساس مقدار عددی محاسبه شده برای هر سلول در گستره ایران نقشه های پهنه بندی a value ، b value و a /bvalue تهیه شده است. تغییرات در مقدار ضرایب لرزه خیزی بیانگر شرایط متفاوت تکتونیکی در منطقه است. در مقایسه و بررسی نقشه ها با یکدیگر نقشه پهنه بندی a /bvalue اطلاعات جامع تری را در مورد وضعیت لرزه خیزی و تکتونیکی منطقه نشان میدهد. بر این اساس پهنه ایران به سه منطقه کلی تقسیم می شود؛ مناطق با پتانسیل لرزه ای پایین ، مناطق با پتانسیل لرزه ای متوسط ، و مناطق با پتانسیل لرزه ای بالا.

zoning of iran region based on changes in seismicity coefficient ratio (a/b)

1 introduction the geological, structural, and seismic characteristics are not the same in various regions in iran. the difference in the magnitude and frequency of seismic events in this area indicates such variations. different researchers have provided various studies and maps of the geologic, tectonic, and seismotectonic situation in this area over the past years (stocklin, 1968; berberian, 1981; alavi, 1991; alavi naini,1972; nowroozi, 1976; berberian,1976; nogol sadat, 1993; mirzaei et al., 1998; tavakoli et al., 1999; zare and memarian, 2000; ansari et al., 2009). the similarity of such maps indicates the close correlation of the geologic, tectonic, seismotectonic situation in this area. among the seismic parameters, the seismicity coefficients such as a and b in the gutenbergrichter (1954) equation are the most important parameters in determining the seismic pattern as well as tectonic and geologic characteristics of different regions. determining these coefficients is necessary in the studies on the risk analysis and earthquake hazard zoning, and the value of these coefficients are assumed to be constant in the seismotectonic states. therefore, providing a zoning map based on these coefficients will play an important role in better identification of the seismic characteristics of the region and applied studies. providing a zoning map of the seismic coefficients of bvalue and avalue and the ratio a/b value for iran region is the main objective of this research. in order to calculate these coefficients and provide zoning maps, initially a complete catalogue of earthquakes occurred from 1900 to 2007 in iran region has been prepared. then, the entire iran region is divided into very small and regular zones, and the seismicity coefficients (b, a) from the gutenbergrichter (1954) equation are calculated individually from this network for each cell. based on the numerical value calculated for each cell in the range of iran, the zoning maps of avalue, bvalue and a/bvalue have been provided. the changes in the value of seismic coefficients indicate different tectonic situations in the region. according to the bvalue and avalue zoning maps, the regions with various seismicity coefficients values can be separated. in the study and comparison of the maps, the zoning map a/bvalue shows more comprehensive information on the seismicity and tectonic situations of the region. as a matter of fact, using this map, the effect of coefficients such as a, b can be seen together in a map. accordingly, iran zone is divided into three general regions; low seismic potential, moderate seismic potential, and high seismic potential regions.   2methodology the first step in examining the seismicity of any region is collecting earthquakes that have previously occurred in that region. in this study, valid domestic and foreign sources were used for preparing the earthquake catalog. before processing the seismic data, it is essential to refine the data of dependent events in order to obtain a poisson distribution of the data. in the present study, the window elimination method which is a standard method based on logarithmic time drawing of aftershocks based on the magnitude of earthquakes (gardner and knopoff, 1974). finally, a catalog with 8090 earthquakes recorded from 1900 to 2007 (107 years) was prepared, serving as the basis and preliminary data of this study.  zoning of iran based on seismic coefficients avalue, bvalue, and a/bvalue is the main objective of this study. these coefficients were prepared for the entire zone of the country, and then regions with equal values were zoned together. thus, to perform numerical calculations and determine seismic coefficients, the entire zone of the country and adjacent regions were divided into 2°*2° cells with 1.5° overlapping. this overlapping among cells plays a vital role in the continuity of data. a total of 1354 cell were resulted, for each the a and bvalues were calculated and attributed to the central point of that cell. in fact, across the entire zone of iran, there are 1354 points with 0.5° distance from one another, and each point has its specific a and bvalue. due to the large number of earthquakes (8090 records) and the large volume of computations (1354 cell keys), a computer program was written using visual basic using which: 1) the recorded seismic data were read; 2) the data belonging to each cell were separated; 3) seismicity computations were performed for that cell, and 4) results were saved in a separate file. therefore, for all the 1354 cells specified across iran and for the data constrained in each cell, the gutenbergrichter line was drawn and a and bvalues were computed for them. finally, after entering the preliminary data (seismic coefficients) in arc gis, zoning maps were prepared using the inverse distance weighting (idw) method which is a common and frequently used method.   3result and discussion based on zoning maps, the a and bvalues of regions with different seismic activities and seismic coefficients were separated. in each map, nine zones with varying seismic coefficients are evident (fig.1). in the zoning map of avalues, the largest numerical value of the frequency of earthquakes belongs to the zagros region, east alborz and kopet dag, and parts of western and northwestern regions, with the numerical value of 5.6 to 8.1. the bvalue zoning map also showed the largest numerical value for the zagros structural trend, east alborz and kopet dag, and western and northwestern regions with the numerical value of 0.49 to 1.7.  these highcoefficient regions surround regions including parts of the central and eastern iran with low seismic coefficients. the increase and decrease in bvalue of various regions of iran indicate the different tectonic conditions and behaviors of this zone (scholz, 1968; mori and abercombie, 1997; manakou and tsapanos, 2000).   figure 1. right, the avalue zoning map across iran; left, the bvalue zoning map across iran     figure 2. the a/bvalue zoning map across iran next, a zoning map of the a/b ratio was also prepared (fig.2). according to bayrak et al. (2002), this map demonstrates the seismic condition and regions with seismic potential better than a and b maps. based on this map, the iran zone includes regions with high (6.4 to 7.4), moderate (5.9 to 6.3), and low (0.0 to 5.8) seismic activity. the high density of colors (high numerical value of a/b) in zagros, strait of hormuz, some eastern parts, kopet dag, eastern and western alborz, and azarbaijan indicate connected seismic belts of tectonically active regions surrounding regions with a low numerical value of a/b (lowdensity colors) or those with a low seismic potential (fig. 2).   4conclusion based on zoning maps of avalue, bvalue, and a/b value prepared for iran, regions with varying seismic coefficients can easily be distinguished from one another. these seismic zones can be compared with sedimentarystructural zones of iran separated based on geological characteristics. this can also indicate the close relationship between seismic and tectonic features in iran. accordingly, iran zone is divided into three general regions; low seismic potential, moderate seismic potential, and high seismic potential regions. the areas having high seismic potential have surrounded low seismic areas; of course, in form of a large belt of high seismic areas such as: 1trend northwest to southeast of zagros, 2 northern strait of hormoz with a northern southern trend corresponding to the naybandsistan faults belt, 3 northeastern part of the makran coast, 4 widespread sections of azerbaijan, 5 western alborz, 6eastern alborz, 7 kopet dagh, and 8parts of the east of the country in the upper part of the lut desert (khorasan).   reference alavi naini, m., 1972. etude geologique delaergiondedjam. geological survey of iran, report no. 23, 1–45 (in french). alavi, m., 1991. sedimentary and structural characteristics of the paleotethys remnants in northeastern iran. geological society of america bulletin, 103(8), 983992. ansari, a., noorzad, a., zafarani, h., 2009. clustering analysis of the seismic catalog of iran. computers geosciences, 35(3), 475486. bayrak, y., yılmaztürk, a., öztürk, s., 2002. lateral variations of the modal (a/b) values for the different regions of the world. journal of geodynamics, 34(5), 653666. berberian, m., 1976. seismotectonic map of iran (1:2 500 000). geological survey of iran, report no. 39. berberian, m., 1981. active faulting and tectonics of iran: zagroshindu kushhimalaya geodynamic evolution 3, 3369. gardner, j.k., knopoff, l., 1974. is the sequence of earthquakes in southern california, with aftershocks removed, poissonian?. bulletin of the seismological society of america 64(5), 13631367. gutenberg, b., richter, c.f., 1954. seismicity of the earth and associated phenimena. princeton university press, princeton, new jersey, u.s.a. manakou, m.v., tsapanos, t.m.,2000. seismicity and seismic hazard parameters evaluation in the island of crete and the surrounding area inferred from mixed data files. tectonophys, 321(1), 157178. mirzaei, n., mengtan, g., yuntai, c., 1998. seismic source regionalization for seismic zoning of iran: major seismotectonic provinces. journal of earthquake prediction research 7, 465495. mori, j., abercrombie, r.e., 1997. depth dependence of earthquake frequency‐magnitude distributions in california: implications for rupture initiation. journal of geophysical research: solid earth, 102 (b7), 1508115090. nogol sadat, m.a.a., 1993. seismotectonic map of iran (scale1:1000000). geological survey of iran. nowroozi, a., 1976. seismotectonic provinces of iran. bulletin of the seismological society of america 66, 12491276. scholz, c.h., 1968. the frequencymagnitude relation of microfracturing in rock and its relation to earthquakes: bulletin of the seismological society of america, 58(1), 399415. stocklin, j., 1968. structural history and tectonics of iran: a review. aapg bulletin, 52(7), 12291258. tavakoli, b., ghaforyashtiany, m., 1999. seismic hazard assessment of iran. annals of geophysics, 42 (6), 1013–1021. zare, m., memarian, h., 2000. simulation of earthquakes intensity in iran. research report of iranian red crescent, tehran, iran, 150 pp. (in persian).