Journal of the Earth and Space Physics

Journal of the Earth and Space Physics

Probabilistic Seismic Hazard Analysis of Hamadan Province Using Monte Carlo Simulation and Evaluation of Standard 2800

Document Type : Research Article

Authors
Ahmad Kiani
Ahmad Sadidkhouy
Department of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran.
10.22059/jesphys.2026.410880.1007759
Abstract
The aim of this study is to carry out a detailed Probabilistic Seismic Hazard Assessment (PSHA) for Hamadan Province in western Iran, with particular emphasis on the seismogenic potential of the Nahavand Fault, and to compare the results with the design provisions of Iranian Seismic Code Standard No. 2800. The region lies within the active Zagros fold-and-thrust belt, where complex fault systems and heterogeneous seismicity patterns necessitate site-specific hazard analyses, especially for critical facilities. In this context, we assembled a homogenized earthquake catalog for the period 1990–2025, converted to moment magnitude using empirical relations developed for the Iranian plateau, and delineated seismogenic sources based on neotectonic, geological, and seismological evidence. Special attention was given to the geometry and segmentation of the Nahavand Fault, including its southern termination, where negative flower structures and distributed deformation indicate ongoing transpressional activity.
Ground motion parameters were estimated through Monte Carlo simulation of 50,000 stochastic earthquake scenarios combined with a logic-tree framework to account for epistemic uncertainties in source characterization and ground motion prediction equations (GMPEs). The logic tree incorporates alternative source models for the Nahavand Fault (linear fault source versus areal seismogenic zone), different magnitude–frequency distributions, and multiple regionally compatible GMPEs. For each realization, synthetic ground motion fields were generated and uniform hazard spectra (UHS) were derived at representative sites across the province for return periods of 475 and 2475 years, corresponding to standard and maximum-considered earthquake levels, respectively.
The results for a 475-year return period indicate a distinctly non-uniform spatial distribution of seismic hazard in Hamadan Province. The highest levels are obtained in Nahavand County, where the maximum peak ground acceleration (PGA) reaches 0.466g. This value is approximately 33.1% higher than the upper bound of the “very high relative hazard” zone defined in Standard 2800 for the region. Elevated hazard levels are also observed in areas located in close proximity to the mapped trace of the Nahavand Fault and its adjacent segments, reflecting the dominant contribution of this structure to seismic risk in the southern part of the province. For the 2475-year return period, the estimated spectral accelerations in the short- and intermediate-period ranges at sites situated near the fault trace exceed the corresponding code-based design values prescribed in Standard 2800, particularly for fundamental periods relevant to low- to mid-rise buildings.
These results should, however, be interpreted in light of the assumptions and limitations inherent in site-specific PSHA. The hazard estimates are conditioned on the adopted source geometry, maximum magnitude, recurrence parameters, selected GMPEs, and the quality and completeness of the input data. Epistemic uncertainties are explicitly treated through the logic-tree approach, while aleatory variability is represented by the dispersion in ground motion predictions; nevertheless, residual model uncertainty remains. Consequently, the findings are not intended to imply a general inadequacy or invalidation of the national seismic zonation framework or the philosophy of Standard 2800. Rather, they underscore the importance of conducting detailed local hazard assessments for critical structures, lifelines, and strategic facilities located in the vicinity of major active faults such as the Nahavand Fault. In such contexts, reliance solely on broad-scale national zoning maps may underestimate the true site-specific hazard, whereas complementary PSHA can provide more realistic input motions for performance-based seismic design and risk mitigation.
Keywords
Probabilistic Seismic Hazard Assessment (PSHA)
Monte Carlo simulation
logic tree
Nahavand Fault
Hamadan Province
Standard No. 2800
Uniform Hazard Spectrum
Subjects
ایزدی کیان، ل.؛ علیپور، ر. و شهریاری، ی. (1403). بررسی ساختاری و فرکتالی پایانه جنوبی گسل نهاوند. زمین شناسی کاربردی پیشرفته، 14(3)، 578-597.
بیت‌اللهی، ع. (۱۳۹۹). گزارش زمین‌لرزه‌های مهرماه ۱۳۹۹ نهاوند (استان همدان). گزارش فنی، تهران: بخش زلزله‌شناسی مهندسی و خطرپذیری، مرکز تحقیقات راه، مسکن و شهرسازی.
حیدری ر. و میرزایی ن. (1388). الگوی لرزه‌زمین‌ساختی گسل اصلی عهد حاضر زاگرس بین 33 تا 35 درجه عرض‌ شمالی، مجله فیزیک زمین و فضا، 35(3)، 83-96.
مرکز تحقیقات راه، مسکن و شهرسازی (۱۳۹۳). آیین‌نامه طراحی ساختمان‌ها در برابر زلزله (استاندارد ۲۸۰۰) (ویرایش چهارم). نشریه شماره 485. تهران: نشر مرکز تحقیقات راه، مسکن و شهرسازی.
Aki, K. (1965). Maximum likelihood estimate of b in the formula log N= a-bM and its confidence limits. Bull. Earthquake Res. Inst., Tokyo Univ., 43, 237-239.‏‏
Akkar, S., Sandıkkaya, M. A., & Bommer, J. J. (2014). Empirical ground-motion models for point-and extended-source crustal earthquake scenarios in Europe and the Middle East. Bulletin of earthquake engineering, 12(1), 359-387.‏
Alipoor, R., & Aghayani, E. (2025). The NW extension of the Nahavand Fault along the Main Recent Fault, west Iran. PLoS One, 20(6), e0326603.‏
Bender, B. (1983). Maximum likelihood estimation of b values for magnitude grouped data. Bulletin of the Seismological Society of America, 73(3), 831-851.‏
Bommer, J. J., & Stafford, P. J. (2016). Seismic hazard and earthquake actions. In Seismic design of buildings to eurocode 8 (pp. 21-54). CRC Press.‏
Cornell, C. A. (1968). Engineering seismic risk analysis. Bulletin of the seismological society of America, 58(5), 1583-1606.‏
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.‏
Ghasemi, H., Zare, M., Fukushima, Y., & Koketsu, K. (2009). An empirical spectral ground-motion model for Iran. Journal of seismology, 13(4), 499-515.‏
Kroese, D. P., Brereton, T., Taimre, T., & Botev, Z. I. (2014). Why the Monte Carlo method is so important today. Wiley Interdisciplinary Reviews: Computational Statistics, 6(6), 386-392.‏
Mousavi-Bafrouei, S. H., Mirzaei, N., & Shabani, E. (2014). A declustered earthquake catalog for the Iranian Plateau. Annals of geophysics, 57(6).‏
Reasenberg, P. (1985). Second‐order moment of central California seismicity, 1969–1982. Journal of Geophysical Research: Solid Earth, 90(B7), 5479-5495.‏
Talebi, M., Zare, M., Peresan, A., & Ansari, A. (2017). Long-term probabilistic forecast for M≥ 5.0 earthquakes in Iran. Pure and Applied Geophysics, 174(4), 1561-1580.‏
Zafarani, H., Luzi, L., Lanzano, G., & Soghrat, M. R. (2018). Empirical equations for the prediction of PGA and pseudo spectral accelerations using Iranian strong-motion data. Journal of Seismology, 22(1), 263-285
Wells, D. L. & Coppersmith, K. J. (1994). New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement, Bulletin of the Seismological Society of America (1994) 84 (4), 974–1002.
Journal of the Earth and Space Physics
Volume 52, Issue 2
online publication date: 5 September 2026
Summer 2026
Pages 259-270