th The 19 Iranian Geophysical Conference–November, 2020, pages 548-551

Is WNW striking faults a seismic source threatening ?

Abbassi, M.R1, Ghods, A.R2,

1 Intenational Institute of Earthquake Engineering and Seismology, Tehran/, [email protected] 2 Department of Earth Sciences, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan 45137- 66731, Iran

ABSTRACT Tehran lies on the southern flank of the Central Alborz, an active mountain belt characterized by many historical earthquakes, some of which have affected Tehran itself. It is characterized by very low seismic activity. Furthermore a rapid urbanization in last two decades made it difficult to conduct geological study. A relocation of seismic events together with kinematic study allowed locating some new active fault trends, which could be considered as seismic source. The newly recognized faults are striking WNW, forming the main structural grain of western branch of the central Alborz. Based on aerial photos and field observation these fault trends arranged in en- echelon manner occupying a considerable fault zone width, especially in alluvial deposits.

Keywords: Tehran, seismic source, active faulting, WNW-striking fault

INTRODUCTION In last two decades many researchers have shown the importance of E-W striking faults as potential seismic source threatening the Grater Tehran and surrounding area, including former major ancient Iranian city called Ray (Rhaga) lying in southern Tehran. The historical city of Ray has experienced three major earthquakes evidenced by historical data (Ambraseys and Melville, 1982; Berberian et al., 1985; Berberian and Yeats, 1999). It is suggested that events of 855– 856 AD (M~7.1) and 1177 AD (M~7.2) may be associated with the Kahrizak, North Rey and South Rey. The first paleoseismological study in Tehran conducted on Kahrizak fault by De Martini et al. (1998) showed that the fault activity and its related destructive events (855 A.D., 1177 A.D., 312-280 B.C. and 958 A.D). In a newer paleoseismological study Nazari et al., (2010) rejected the events assigned to Kahrizak fault, North Ray fault and South Ray and considered the three faults as features related to ancient lake terraces and ancient beaches. On the other hand Azadi et al., (2010) using three different geophysical methods, seismic refraction, electrical resistivity and magnetic techniques on Kahrizak fault have revealed the geometry and activity of this fault. The inherited WNW-directed faults display various morphological features. The length of these faults, observed on aerial photos, within Tehran's piedmont, varies between 2.6 to 6.5 Km, arranged in an en-echelon manner. Morphological features of the WNW-directed faults display controversial faulting mechanisms (Abbassi and Farbod, 2009) on aerial photos and cannot be explained by a single stress direction. Whereas in rock mass cropping out in northern and eastern Tehran they are well known as Emamzadeh Davoud, Poorkan-Vardij, western termination of Mosha Fault and Ghasr-e Firuzeh. A newer paleoseismological study on WNW striking part of North Tehran fault revealed three historical events (Ritz et al., 2012). However, an alternative data to shed a light on the true nature of fault activity in Tehran could be the instrumental data. Hence due to the extensive urbanization the morphological features of WNW trending faults hinders any geological study. According to the present day stress prevailing in the South Central Alborz obtained by P-axis of focal mechanisms and fault slip data in relative young alluvial deposits (Abbassi, M. R., and E. Shabanian ,1999; Abbassi, M.R. et al. 2003) NE-directed stress proposes a compressive faulting mechanisms associated with WNW-striking faults.

549

Methodology The assessment of seismic hazard depends upon understanding of activity, mechanism and trends of faults affecting an area. The rapid urbanization of Tehran with remarkably concentration of people needs to be studied by new studies considering all possible active fault trends.

Figure 1.SRTM image of studied area. a) WNW faulting trend and their names. Yellow dots show the relocated seismic events. b) Kinematic stations are shown with their related numbers. In black are focal mechanisms from this study, numbers refer the events (table 2), in gray are focal mechanisms from Soltani Moghadam et al. 2018 (for the numbers see table 1).

550

Table1. Source mechanism data.

Date Lat. Lang. Ms Plane Dip Rake Ref. 1957.07.02 36.14 52.70 6.8 119 47 80 M72 2007.06.18 34.47 50.79 5.5 144 66 124 CMT 1 35.67 51.59 2.0 122 48 118 S 2 35.57 51.52 3.7 284 47.9 39.3 S 3 35.46 51.45 2.6 118.5 60 90 S M72 first-motion fault-plane solution from McKenzie (1972). HVD Harvard Centroid Moment Tensor Catalogue. Soltani Moghadam et al., (2018)

Table2. Source mechanism data (this study) Date Lat. Lang. M Plane Dip Rake Numb. of event 2003.03.09 35.787 51.516 3.2 231 36 31 4 2008.05.26 35.21 51.497 3.2 52 67 9 17 2009.10.17 35.569 51.505 3.2 316 40 71 25 2013.08.21 35.372 51.798 3.1 63 38 29 50 The E-W-trending active faulting in Tehran is delineated by their morphological features recognized on aerial photos of 1955. This study examines WNW-trending faults as a possible seismic source responsible for historical earthquakes. The assessment of seismic hazard depends upon understanding of activity, mechanism and trends of faults affecting an area. The rapid urbanization of Tehran with remarkable concentration of people needs to be studied by new studies considering all possible active fault trends. The E-W- trending active faulting in Tehran is delineated by their morphological features recognized on aerial photos of 1955. Generally active faults with mainly vertical movements are not easy to recognize in aerial photos. Another obstacle to trace active faulting in urban area is human activity by demolishing the morphology. To overcome these problems we have used kinematic study conducted in Late Pleistocene deposits and relocated seismic events (Fig.1a&b). To evaluate seismicity in the Tehran region, we have collected a comprehensive data set made of all permanent and temporary networks. We use four permanent networks including the Iranian National Seismic Network (IRSC) belonging to Institute of Geophysics University of Tehran (IGUT) operating since 1996, the Iranian National Seismograph Network (INSN) belonging to International Institute of Earthquake Engineering and Seismology (IIEES) operating since 1998, the Tehran Disaster Mitigation and Management Organization (TDMMO) seismic network operating since 2004, the Iranian Building and Housing Research Center (BHRC) belonging to Road, Housing & Urban Development Research Center operating since 1973. The BHRC stations do not have calibrated timing systems so we utilize differential Sg-Pg phase readings. The temporary network include the East Alborz temporary network (EATN) belonging to International Institute of Earthquake Engineering and Seismology (IIEES) The database has two important aspects, namely long duration (1996-2019) and better spatial coverage. The Central Alborz is a slow slip region (Nilfouroushan et al, 2003; Vernant et al., 2004; Djamour et al, 2010) and consequently the recurrence time of earthquakes are long. By expanding the time period of our catalog, we have increased the number of well-located events and focal mechanisms by as much as possible. By combining the data from different networks, we have improved the azimuthal coverage and decreased the possibility of systematic location biases caused by deviations of the assumed travel time model from the true Earth. All the events shown in Figure 1 have Pg azimuthal gap of less than 180 degree. The focal mechanism shown in Fig. 1 obtained using first polarity method and grouped as quality A (black focals). The aim of this study is in first line, to shed a light on faulting mechanism of WNW striking bordering the southern flank of the Central Alborz, by combination of seismic data and kinematic study conducted in the Late Pleistocene deposits. The morphological evidence for proposed WNW faulting are almost vanished in Tehran City, therefore we relay on seismic and kinematic data.

551

CONCLUSION(S) The compressive faulting mechanism of WNW-trending and moderately dipping faults is in accordance with the NE-directed present-day stress direction. Based on kinematic study conducted in Late Pleistocene alluvial deposit on - fault (Fig.1), these faulting trends tends normally to have a minor left-lateral horizontal component, which is in accordance with seismological data obtained by this and previous studies. However, some of the focal mechanisms show a minor right-lateral horizontal component (Fig. 1a &b; Tables 1 &2). The relocated events (Fig.1a) and kinematic studies measured mainly in Late Pleistocene deposits (Fig.1b) suggest several WNW-trending fault affecting Tehran city, namely as Karaj-Rey, Kan- Rey and Chitgar-Rey. The Karaj-Rey fault is the mountain frontal faulting, with a compressive focal mechanism (event 25, Fig 1a and Table 2). A major part of recent shortening occurs along WNW striking faults in the lower elevated south Central Alborz and the central Iran.

REFERENCES Abbassi, M.R., and Shabanian, E., 1999. Evolution of the stress field in the Tehran region during the Quaternary. In: Third International Conference on Seismology and Earthquake Engineering. Tehran-Iran. Abbassi, M. R., E. Shabanian, Y. Farbod, K. Feghhi, and H. Tabassi (2003). The state of contemporary stress in the southern flank of central Alborz [abstract in English], Rep. 81-2003- 7, Int. Inst. of Earthquake Eng. and Seismol., Tehran. Abbassi, M.R. and Farbod, Y., 2009. Faulting and folding in quaternary deposits of Tehran’s piedmont (Iran). J. Asian Earth Sci. 34, 522–531. https://doi.org/10.1016/j.jseaes. 2008.08.001. Amberaseyas, N.N. and Melville, C.P., 1982. A history of Persian earthquakes. Cambridge University Press. Cambridge, 219 p. Azadi, A., Hessami, K., Javan-Doloei, Gh., 2010. Integrated geophysical methods for determining geometry of the Kahrizak Fault, Tehran, Iran. Nat Hazards 54:813–825. DOI 10.1007/s11069-010-9506-9 Berberian, M., Qorashi, M., Arzhang-Ravesh, B., Mohajer-Ashjai, A.,1985. Recent tectonics, seismotectonics, and earthquake-fault hazard study in the greater Tehran region: contibution to the seismotectonics of Iran, part V. Geological Survey of Iran 56, 316 (in Persian). Berberian, M., and Yeats, R.S., 1999. Patterns of historical earthquake rupture in the Iranian Plateau. Bulletin of the Seismological Society of America, 89, 120-139. DeMartini PM, Hessami K, Pantosti D, D’Addezio G, Alinaghi H, Ghafory-Ashtiani M (1998) A geologic contribution to the evaluation of the seismic potential of the Kahrizak fault (Tehran, Iran). Tectonophysics, 287:187–199. Djamour, Y. et al., 2010. GPS and gravity constraints on continental deformation in the Alborz mountain range, Iran, Geophys. J. Int., 183(3), 1287–1301. McKenzie, D. P., 1972. Active tectonics of the Mediterranean region, Geophys. J. R. astr. SOC., 30,109- 185. Nazari, H., Ritz, J.F., Salamati, R., Shahidi, A., Habibi, H., Ghorashi, M. & Bavandpur, A.K., 2010. Distinguishing between fault scarps and shorelines: the question of the nature of the Kahrizak, North Rey and South Rey features in the Tehran plain (Iran), Terra Nova, 22(3), 227– 237. Nilfouroushan, F. et al., 2003. GPS network monitors the Arabia-Eurasia collision deformation in Iran. Journal of Geodesy, 77: 411–422 DOI 10.1007/s00190-003-0326-5 Ritz, J.F. et al., 2012. Paleoearthquakes of the past 30 000 years along the North Tehran Fault (Iran), J. geophys. Res., 117(B6), doi:10.1029/2012JB009147. Soltani Moghadam, S., Sepanlo, K., Kheyri, M. 2018 Velocity model calculation and seismicity study of last decade on Tehran and high Alborz elevations. Iranian Journal of Geophysics, Vol. 12, No 2, 2018, P. 6 Vernant, P. et al. 2004. Present-day crustal deformation and plate kinematics in theMiddle East constrained by GPS measurements in Iran and northern Oman, Geophys. J. Int., 157(1), 381–398.