Document Type : Research Paper

Authors

1 Master's Degree in Geophysics-Earthquake Engineering

2 Professor of Malek Ashtar University of Technology

3 Master’s Degree in Artificial Intelligence

Abstract

Extended Abstract
Introduction
Unfortunately, seismic data recorded globally during the last fifty years does not include every type of wave propagation conditions in the environment, types of construction, the rupture process on the fault, and the geometrical relationship between the construction and the fault. This is especially seen in near-field regions. Before the 1999 Chi-chi earthquake in Taiwan and the 1999 Izmit earthquake in Turkey, there were only about 20 records of earthquakes with a magnitude greater than 7 at a distance of less than 20 kilometers from the fault.
The Turkish earthquake added 5 records and the Thai earthquake added 65 records to this collection, but only two fault rupture scenarios were added to our knowledge, while thousands of other possible scenarios may occur. Thus, seismologists and earthquake engineers have tried to estimate parameters related to strong near-field motions of the earth with an acceptable confidence using various experimental and theoretical simulation methods.           
Methods
In earthquake engineering and seismology, earthquake phenomenon and the resulting movements are generally investigated and analyzed using dynamic and kinematic methods. Seismological models and problems are thus divided into two categories: Kinematic models which are based on slip distribution and do not take the state of stress on the fault into account. Dynamic models deal with the physics of fault rupture and its causes. Simulation methods are also divided into three main categories: deterministic (low frequencies), stochastic (high frequencies) and hybrid (broad band) methods.
Generally speaking, simulating strong ground motion plays an important role in the estimation of related parameters especially in regions lacking such data. Accelerographs are used to simulate strong ground motions. The present study has introduced, investigated and validated two methods: decisive simulation models (Discrete-Wave Number and Finite Fault) and Finite Fault models. It also explains how the simulated recording are produced for near-field (less than 20 km to a seismogenic fault) and far-field events, presents attenuation relationships for the Zagros seismotectonics region, and predicts parameters of strong ground motions.
Results & Discussion
Due to the special geological conditions and the existence of many active faults in Iran, our country is considered to be located in an earthquake-prone region. Zagros region is considered to be the most earthquake-prone region of Iran. Finite fault modeling combines various aspects of plate source with the ground motion model based on point source. Since previously mentioned limitations are not naturally present in finite fault modelling, the method takes geometry of the fault and the directivity effect into account. Time delay method and the sum of accelerations recorded in maps of a two-dimensional network are used for simulation in finite fault model. The fault plate is divided into various elements and a minor event is simulated for each one. The overall seismic acceleration equals the sum of the effects of these minor events. The strong ground motions in each micro-fault are calculated using the random point source method and then summed up at the desired point with an appropriate time delay to obtain the ground motion of the entire fault.
Previous geological and seismic studies of each seismic region are used to determine the key parameters of the simulation input. To produce a comprehensive database, a significant number of stations are taken into account around the fault based on different hypotheses and artificial accelerograms are produced in accordance with the seismological parameters of the region. A suitable function is then selected and an attenuation relationship is fitted. The simulation results and the resulting attenuation relationship are then compared with valid global attenuation relationships and their consistency (compliance percentage) is investigated.
 Conclusion
The present study has produced a wide range of simulated records (about 20 thousand records) for Zagros seismotectonics region. Thus, the resulting relationships will hopefully have sufficient accuracy and efficiency to be used in structure designing and urban development. It worth noting that the regression correlation coefficient (R-Square) was above 0.95 in all fits.
These attenuation relationships can provide a new perspective on site selection, and help us in understanding the dynamic behavior of structures, and the development of various infrastructure. They also help urban managers to predict and reduce earthquake damages.

Keywords

Main Subjects

1- حسنخانی، زعفرانی؛ علی، حمید(1393). "پیش‌ بینی خصوصیات پالس‌های جهت‌پذیری حوزه نزدیک با شبیه‌سازی پدیده زلزله". هشتمین کنگره ملی مهندسی عمران، بابل 1393.
2- زعفرانی، نورزاد، برگی؛ حمید، اسدا..، خسرو(1386). "شبیه‌سازی حرکات ثبت شده درزلزله دی‌ماه 1382بم، به روش تصادفی گسل با ابعاد محدود و بررسی کمّی نقش چشمه لرزه‌زا در شکل‌گیری توزیع خرابی مشاهده شده". نشریه دانشکده فنی، جلد41، شماره6، دی‌ماه 1386، از صفحه753 تا صفحه 764.
3- زعفرانی، نورزاد؛ حمید، اسدا... (1392). لرزه‌شناسی مهندسی. تهران: انتشارات دانشگاه تهران.
4- Akkar, S, and Bommer, J. J. (2010). ‘Empirical equations for the prediction of PGA, PGVand spectral accelerations in Europe, the Mediterranean region, and the Middle East, Seismol’. Res. Lett. 81, 195–206.
5- Arora, S., Joshi, A., Kumari, P. et al. Strong ground motion simulation techniques—a review in world context. Arab J Geosci 13, 673 (2020). https://doi.org/10.1007/s12517-020-05583-5
6- Beresnev, I., and G. Atkinson (1998a). Stochastic finite-fault modeling of ground motions from the 1994 Northridge, California earthquake, part I: Validation on rock sites, Bull. Seism. Soc. Am. 88, 1392–1401.
7- Boore, D. M., and W. B. Joyner (1997). ‘Site amplification for generic rock sites’. Bull. Seismol. Soc. Am. 87, 327–341.
8- Cheloni. Daniele, Aybige Akinci .” Source modelling and strong ground motion simulations for the 24 January 2020,Mw6.8Elazığ earthquake,Turkey”Geophysical Journal International, Volume 223, Issue 2, November 2020,Pages1054–1068,https://doi.org/10.1093/gji/ggaa350.
9- Frankel. Arthur (2009). ‘A Constant Stress – Drop Model for Producing Broad Band Synthetic Seismograms: Comparison with the Next Generation Attenuation Relations’. Bull. Seismol. Soc. Am., 99:2A, 664–680.
10- Gail M. Atkinson and David M. Boore. ‘Earthquake Ground-Motion Prediction Equations for Eastern North America’. Bulletin of the Seismological Society of America, Vol. 96, No. 6, pp. 2181–2205, December 2006, doi: 10.1785/0120050245.
11- Gallovič, F., and Brokešová, J., (2004). ‘On strong ground motion synthesis with k-2 slip distributions’. J. Seismol. 8, 211–224.
12- G‌h‌a‌e‌m‌i‌a‌n. M,  H. Zafarani.” s‌e‌i‌s‌m‌i‌c a‌n‌a‌l‌y‌s‌i‌s o‌f d‌a‌m-f‌o‌u‌n‌d‌a‌t‌i‌o‌n-r‌e‌s‌e‌r‌v‌o‌i‌r-f‌a‌u‌l‌t s‌y‌s‌t‌e‌m b‌y h‌y‌b‌r‌i‌d d‌i‌s‌c‌r‌e‌t‌e w‌a‌v‌e‌n‌u‌m‌b‌e‌r-f‌i‌n‌i‌t‌e e‌l‌e‌m‌e‌n‌t m‌e‌t‌h‌o‌d”. Volume 38.2, Issue 1.2 - Serial Number 1June 2022 Pages 39-50.
13- Hanks, T., and R. McGuire (1981). The character of high-frequency strong ground motion, Bull. Seism. Soc. Am. 71, 2071–2095.
14- Hartzell, S. (1978). Earthquake aftershocks as Green’s functions, Geophys. Res. Lett. 5, 1–14.
15- Herrero, A., and Bernard, P., (1994). ‘A kinematic self-similar rupture process for earthquakes’. Bull. Seismol. Soc. Am. 84, 1216–1228.
16-  Hutchings, L., Ioannidou, E., Foxall, W., Voulgaris, N., Savy, J., Kalogeras, I.,    Scognamiglio, L.,andStavrakakis, G., (2007).’ A physically based strong ground-motion prediction methodology; application to PSHA and the 1999 Mw~6.0 Athens earthquake’, Geophys. J. Int. 168,659–680.
17- Irikura, K., 1992, The construction of large earthquake by a superposition of small events, Earthquake Engineering. Tenth world Coference. Balkema. Rotterdam.
18- Jackson, J. A., Fitch, T. J. and McKenzi, D. P., 1981, Active thrusting and the evolution of the Zagros fold belt, in: McClay, K. R. and Price, N. J. (eds.), Thrust and Nappe Tectonics, Geol. Soc. Spec. Pub., 9, 371-379.
19- Karimiparidari, S. , Zare, M. , Memarian,H., Kijko, A., 2013. Iranian earthquakes, a uniform catalog with moment magnitudes. Journal of seismology, 17: 897-911. DOI: 10.1007/s10950-0139360-9
20- Mirzaei, N., Gao, M., Chen, Y. T., 1998, Seismic source regionalization for seismic zoning of Iran: major seismotectonic provinces. J. Earthquake prediction Research, 7, 465-495.
21- Molnar, P. and Chen, W. P., 1982, Seismicity and mountain building, in: mountain building Processes, Hsued, K. J. (ed.), Academic Press, 41-57.
22- Motazedian, D., and G. Atkinson (2005). Stochastic finite-fault model based on dynamic corner frequency, Bull. Seism. Soc. Am. 95, 995–1010.
23- Nicholas J. Gregor, Walter J. Silva, Ivan G. Wong, and Robert R. Youngs (2002). ‘ Ground-MotionAttenuation Relationships for Cascadia Subduction Zone Megathrust Earthquakes Based on a Stochastic Finite-Fault Model’. Bulletin of the Seismological Society of America, Vol. 92, No. 5, pp. 1923–1932, June 2002.
24- N S H Harith, P J Ramadhansyah, M I Adiyanto and N I Ramli(2020), “Ground Motion Observation of Sabah Earthquakes on the Use of Next Generation Attenuation (NGA) Ground-Motion Models” IOP Conference Series: Earth and Environmental Science, Volume 682, 4th National Conference on Wind & Earthquake Engineering 16 - 17 October 2020, Putrajaya, Malaysia.
25- Sepehr, M., and Cosgrove, J. W., 2005. Role of the Kazerun Fault Zone in the formation and deformation of the Zagros Fold-Thrušt Belt, Iran.
 26- Shakiba. Somayeh, Habib Rahimi, Mehrdad Pakzad (2018) “ Simulation of Strong Ground Motion of Bam Earthquake Using Stochastic Finite Fault Method” Proceedings of the 18th Iranian Geophysical Conference, May 2018, pages 642-645.
27- Somerville, P., Irikura, K., Graves, R., Sawada, S., Wald, D., Abrahamson, N., Iwasaki, Y., Kagawa, T., Smith, N., and Kowada, A., (1999). ‘Characterizing earthquake slip models for the prediction of strong ground motion’. Seism. Res. Lett. 70, 59–80.
28- Spudich. Paul and Brian S. J. Chiou (2008). ‘Directivity in NGA Earthquake Ground Motions: Analysis Using Isochrone Theory’. Earthquake Spectra: February 2008, Vol. 24, No. 1, pp. 279-298.
29- Stephen H. Hartzell and Thomas H. Heaton (1983).’ inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 imperial valley, california, earthquake ‘ Bulletin of the Seismological Society of America, Vol 73, No 6, pp 1553-1583, December 1983.
30- Wells, D. L., and K. J. Coppersmith (1994). ‘New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement’. Bull. Seism. Soc. Am. 84, 974–1002.
31- Zafarani .H, A. Noorzad and A. Ansari,(2005). ‘ Generation of near-fault response spectrum for a large dam in Iran ‘. Hydropower and Dams,12, Issue 4. 51-55.
32- Zafarani.H & B. Hassani (2010), ‘Site response and source spectra of S waves in the Zagros region, Iran’. Journal of Seismology ISSN 1383-4649 Volume 17 Number 2 J Seismol (2013) 17:645-666 DOI 10.1007/s10950-012-9344-1.
33- Zafarani .H and M. Soghrat, (2012). ‘Simulation of Ground Motion in the Zagros Region of Iran Using the Specific Barrier Model and the Stochastic Method’. Bulletin of the Seismological Society of America, Vol. 102, No. 5, pp. 2031–2045, October 2012, doi: 10.1785/0120110315
34- Zafarani .H, n, H. Vahidifard, A. Ansari (2012). ‘Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area, Iran ‘. Soil Dynamics and Earthquake Engineering 43, 342–354.
35- Zafarani. Hamid, Hesam Vahidifard, and Anooshirvan Ansar I (2013) .’ Prediction of Broadband Ground-Motion Time Histories: The Case of Tehran, Iran’.Earthquake Spectra 29 (2), 633-660.