Zoning of Iran region based on changes in seismicity coefficient ratio (a/b)

Abstract

 
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 Gutenberg-Richter (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 b-Value and a-Value 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 Gutenberg-Richter (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 a-Value, b-Value and a/b-Value have been provided. The changes in the value of seismic coefficients indicate different tectonic situations in the region. According to the b-Value and a-Value zoning maps, the regions with various seismicity coefficients values can be separated. In the study and comparison of the maps, the zoning map a/b-Value 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.
 
2-Methodology
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 a-value, b-value, and a/b-value 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 b-values 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 b-value. 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 Gutenberg-Richter line was drawn and a- and b-values 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.
 
3-Result and discussion
Based on zoning maps, the a- and b-values 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 a-values, 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 b-value 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 high-coefficient regions surround regions including parts of the Central and Eastern Iran with low seismic coefficients. The increase and decrease in b-value 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 a-value zoning map across Iran; Left, the b-value zoning map across Iran
 
 
Figure 2. The a/b-value 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 (low-density colors) or those with a low seismic potential (Fig. 2).
 
4-Conclusion
Based on zoning maps of a-value, b-value, 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 sedimentary-structural 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: 1-trend northwest to southeast of Zagros, 2- northern Strait of Hormoz with a northern southern trend corresponding to the Nayband-Sistan faults belt, 3- northeastern part of the Makran coast, 4- widespread sections of Azerbaijan, 5- Western Alborz, 6-Eastern Alborz, 7 Kopet Dagh, and 8-parts 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 Paleo-Tethys remnants in northeastern Iran. Geological Society of America Bulletin, 103(8), 983-992.
Ansari, A., Noorzad, A., Zafarani, H., 2009. Clustering analysis of the seismic catalog of Iran. Computers & Geosciences, 35(3), 475-486.
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), 653-666.
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: Zagros-Hindu Kush-Himalaya Geodynamic Evolution 3, 33-69.
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), 1363-1367.
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), 157-178.
Mirzaei, N., Mengtan, G., Yuntai, C., 1998. Seismic source regionalization for seismic zoning of Iran: major seismotectonic provinces. Journal of earthquake prediction research 7, 465-495.
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), 15081-15090.
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, 1249-1276.
Scholz, C.H., 1968. The frequency-magnitude relation of microfracturing in rock and its relation to earthquakes: Bulletin of the seismological society of America, 58(1), 399-415.
Stocklin, J., 1968. Structural history and tectonics of Iran: a review. AAPG Bulletin, 52(7), 1229-1258.
Tavakoli, B., Ghafory-Ashtiany, 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).
 

Keywords

Main Subjects


آقانباتی، ع.،  1383، زمین شناسی ایران: انتشارات سازمان زمین شناسی و اکتشافات معدنی کشور، تهران، 606 ص.
افتخار نژاد، ج.، 1359،پهنه ساختاری ایران در ارتباط با حوضه های رسوبی: انجمن نفت ایران، شماره 82، صفحه 19-28.
نبوی، م. ح.، 1355، دیباچه ای بر زمین شناسی ایران: سازمان زمین شناسی کشور، گزارش شماره 38، صفحه 110.
نیری، ع.  و همکاران، 1375، استان های لرزه زمین ساخت ایران زمین: کمیته ملی سدهای بزرگ ایران، (IRCOLD)،نشریه 12.
Ahmadi,G., Mostaghel,N. and Nowroozi,A.A.,1989, Earthquake Risk analysis of Iran-V: probabilistic seismic risk for various peak  ground accelerations: Iranian Journal of Science and Technology,13:115-156.
AlaviNainiM.,1972, Etudegeologiquedelaergiondedjam: Geological Survey of Iran Reports23, pp.1–45 (in French).
 Alavi, M., 1991,Sedimentary and structural characteristics of the Paleo-Tethys remnants in northeastern Iran: Geological Society of America Bulletin, 103(8), pp.983-992
Allen, C.R., St. Amand, P., Richter, C.F. and Nordquist, J., 1965, Relationship between seismicity and geologic structure in the southern California region: Bulletin of the Seismological Society of America, 55(4), pp.753-797.
Ansari, A., Noorzad, A. and Zafarani, H., 2009, Clustering analysis of the seismic catalog of Iran: Computers & Geosciences, 35(3), pp.475-486.
Ambraseys, N. N. and Melville, C. P., 1982, A history of persian earthquakes: Cambridge university press.
Barton, D.J., Foulger, G.R., Henderson, J.R. and Julian, B.R., 1999, Frequency–magnitude statistics and spatial correlation dimensions of earthquakes at Long Valley caldera, California: Geophysical Journal International, 138(2), pp.563-570.
Bayrak, Y., Yılmaztürk, A. and Öztürk, S., 2002, Lateral variations of the modal (a/b) values for the different regions of the world: Journal of Geodynamics, 34(5), pp.653-666.
Berberian, M., 1981, Active faulting and tectonics of Iran: Zagros-Hindu Kush-Himalaya Geodynamic Evolution, 3, pp.33-69.
Berberian, M., 1976a, Contribution to Seismotectonics of Iran (part II): Geological Survey of Iran Report No. 39, 518 pp.
Berberian, M., 1976b, Seismotectonic map of Iran (1:2 500 000): Geological Survey of Iran Report No. 39.
Berberian, M., 1995, Master “blind” thrust faults hidden under the Zagros folds: active basement tectonics and surface morphotectonics: Tectonophysics, 241(3-4), pp.193-224.
Bhattacharya, P.M., Majumdar, R.K. and Kayal, J.R., 2002, Fractal dimension and b-value mapping in northeast India: Current Science, pp.1486-1491.
Boccaletti, M., Dainelli, P., 1982, Il sistema regmatico neogenicoquaternario nell’area mediterranea: esempio di deformazione plastico-rigida post-collisionale: Mem. Soc. Geol. Ital. 24, 465– 482.
Bonini, M., Corti, G., Sokoutis, D., Vannucci, G., Gasperini, P. and Cloetingh, S., 2003, Insights from scaled analogue modelling into the seismotectonics of the Iranian region: Tectonophysics, 376(3-4), pp.137-149.
Cao, A. and Gao, S.S., 2002, Temporal variation of seismic b‐values beneath northeastern Japan island arc: Geophysical research letters, 29(9), pp.48-1.
DeMets, C., Gordon, R.G., Argus, D.F. and Stein, S., 1990, Current plate motions: Geophysical journal international, 101(2), pp.425-478.
Gardner, J.K. and Knopoff, L., 1974, Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian?: Bulletin of the Seismological Society of America, 64(5), pp.1363-1367.
Gutenberg, B. and Richter, C.F., 1954, Seismicity of the earth and Associated Phenimena: Princeton Univ. Press. Princeton, New Jersey, U.S.A.
Hatzidimitriou, P.M., Papadimitriou, E.E., Mountrakis, D.M. and Papazachos, B.C., 1985, The seismic parameter b of the frequency-magnitude relation and its association with the geological zones in the area of Greece: Tectonophysics, 120(1-2), pp.141-151.
Hashemi, N., 2009, An assessment of the Spatial Variation of the Seismic b_value across Iran: Scientific Quarterly Journal, Geoscience, Vol.18,No 72, Summer 2009.
Jackson, M.P.A. ed., 1990, Salt diapirs of the Great Kavir, central Iran (Vol. 177): Geological Society of America.
Jackson, J., Haines, J. and Holt, W., 1995, The accommodation of Arabia‐Eurasia plate convergence in Iran: Journal of Geophysical Research: Solid Earth, 100(B8), pp.15205-15219.
Jackson, J. and McKenzie, D., 1984, Active tectonics of the Alpine—Himalayan Belt between western Turkey and Pakistan: Geophysical Journal International, 77(1), pp.185-264.
Jackson, J. and McKenzie, D., 1988, The relationship between plate motions and seismic moment tensors, and the rates of active deformation in the Mediterranean and Middle East: Geophysical Journal International, 93(1), pp.45-73.
Jackson, J.A., 2002, Strength of the continental lithosphere: time to abandon the jelly sandwich?: GSA today, 12, pp.4-10.
Jestin, F., Huchon, P. and Gaulier, J.M., 1994, The Somalia plate and the East African Rift System: present-day kinematics: Geophysical Journal International, 116(3), pp.637-654.
Kalyoncuoglu, U.Y., 2007, Evaluation of seismicity and seismic hazard parameters in Turkey and surrounding area using a new approach to the Gutenberg–Richter relation: Journal of Seismology, 11(2), pp.131-148.
Khan, P.K., 2005, Mapping of b-value beneath the Shillong Plateau: Gondwana Research, 8(2), pp.271-276.
Mandal, P. and Rastogi, B.K., 2005, Self-organized fractal seismicity and b value of aftershocks of the 2001 Bhuj earthquake in Kutch (India): pure and applied geophysics, 162(1), pp.53-72.
McKenzie, D., 1972, Active tectonics of the Mediterranean region: Geophysical Journal of the Royal Astronomical Society, 30(2), pp.109-185.
Mirzaei, N., Mengtan, G. and Yuntai, C., 1998, Seismic source regionalization for seismic zoning of Iran: major seismotectonic provinces: Journal of earthquake prediction research, 7, pp.465-495.
Mogi, K., 1962, Magnitude-frequency relation for elastic shocks accompanying fractures of various materials and some related problems in earthquakes: Bull. Earthq. Res. Inst., 40, pp.831-853.
Mori, J. and 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), pp.15081-15090.
NogolSadat,M.A.A.,1993, SeismotectonicmapofIran(scale1:1000000): Geological Survey of Iran.
Nowroozi, A., 1976, Seismotectonic Provinces of Iran: Bull. Seism. Soc. Am.66, pp1249-1276.
Nowroozi, A.A. and Ahmadi, G., 1986, Analysis of earthquake risk in Iran based on seismotectonic provinces: Tectonophysics, 122(1-2), pp.89-114.
Priestley, K., Baker, C. and Jackson, J., 1994, Implications of earthquake focal mechanism data for the active tectonics of the South Caspian Basin and surrounding regions: Geophysical Journal International,  118(1), pp.111-141.
Scholz, C.H., 1968, The frequency-magnitude relation of microfracturing in rock and its relation to earthquakes: Bulletin of the seismological society of America, 58(1), pp.399-415.
Stocklin, J., 1968, Structural history and tectonics of Iran: a review. AAPG Bulletin, 52(7), pp.1229-1258.
Tavakoli, B. and Ghafory-Ashtiany, M., 1999, Seismic hazard assessment of Iran: Annals of Geophysics, 42(6), 1013–1021.
Tchalenko, J.S. and Braud, J., 1974, Seismicity and structure of the Zagros (Iran): The Main Recent Fault between 33 and 35 N. Phil. Trans. R. Soc. Lond. A, 277(1262), pp.1-25.
Zamani, A. and Asadi,A.,1995, Distribution of a and b values and quantitative seismicity of Iran: Proceeding of the Second International Conference on Seismology and Earthquake Engineering, Tehran,Iran,187-196.
Zamani, A.,Nedaei,M.,Boostani,R.,2009, Tectonic zoning of Iran based on self- organizing map: J. Appl. Sci. 9, 4099 –4114.
Zare, M. and Memarian, H., 2000, Simulation of earthquakes intensity in Iran: Research Report of Iranian Red Crescent, Tehran, Iran, 150 pp. (in Persian).
Wang, J.H., 1988, b values of shallow earthquakes in Taiwan: Bulletin of the Seismological Society of America, 78(3), pp.1243-1254.
Wiemer, S. and Wyss, M., 1997, Mapping the frequency‐magnitude distribution in asperities: An improved technique to calculate recurrence times? : Journal of Geophysical Research: Solid Earth, 102(B7), pp.15115-15128.
Wiemer, S. and Wyss, M., 2002, Mapping spatial variability of the frequency-magnitude distribution of earthquakes: In Advances in geophysics (Vol. 45, pp. 259-V). Elsevier.