Main Article Content
Abstract
In Indonesia, the choice between area source and smoothed seismicity approaches has been driven by data availability rather than scientific evaluation, and no systematic comparison of both methods at the national scale has been conducted. This study compared both approaches across Indonesia to evaluate their differences in PGA estimation and assess the forecasting skill of the smoothed seismicity model. Both models were applied to the same declustered earthquake catalog with a Magnitude of Completeness of 4.7, b-value of 0.89, and a-value of 8.753, using identical Ground Motion Prediction Equations within a probabilistic seismic hazard analysis framework. Both models consistently identified North Sulawesi as the highest hazard region and Kalimantan as the lowest. The smoothed seismicity model produced higher maximum PGA values of 2.13 g compared to 1.1597 g and 1.5901 g from the area source model at 10% and 2% probability of exceedance in 50 years, respectively. Molchan Diagram validation yielded an Area Skill Score of 0.88 and R-score of 0.656, confirming strong forecasting skill. The two methods are complementary, and their integration within a logic tree framework is recommended for future national-scale PSHA in Indonesia.
Keywords
Article Details

This work is licensed under a Creative Commons Attribution 4.0 International License.
References
- [1] S. Dey Rahmawati. (2025). Spatial Modeling of Earthquake Risk in Sulawesi and Maluku Based on Geological Factors, ESTIMASI: Journal of Statistics and Its Application, vol. 6, no. 2, pp. 2721–379.
- [2] N. R. Ananda and D. Pujiastuti. (2024). Studi Bahaya Seismik dengan Metode PSHA (Probabilistic Seismic Hazard Analysis) di Nusa Tenggara Barat Menggunakan Data Gempa Tahun 1900 - 2023, Jurnal Fisika Unand, vol. 13, no. 5, pp. 665–670.
- [3] M. M. Gallahue et al. (2022). A study on the effect of site response on California seismic hazard map assessment, Front. Earth Sci. (Lausanne)., vol. 10.
- [4] R. I. Hielmy and M. L. Rajagukguk. (2025). Analisis Probabilistik Bahaya Seismik Di Denpasar Dan Sekitarnya Berdasarkan Pendekatan Psha, Kel. Tanah Tinggi, vol. 11, no. 1, pp. 21–28.
- [5] R. Wulandari, C.-H. Chan, J.-C. Gao, and D. H. Natawidjaja. (2024). Probabilistic Seismic Hazard Assessment of Sumatra, Indonesia, Aug. 28.
- [6] Daswita, A., Pujiastuti, D., & Anggono, T. (2023). Studi Bahaya Seismik dengan Metode Probabilistic Seismic Hazard Analysis di Sumatera Barat. Jurnal Fisika Unand, 12(3), 445-451.
- [7] Frankel, A. (1995). Mapping seismic hazard in the central and eastern United States. Seismological Research Letters, 66(4), 8-21.
- [8] Taroni, M., & Akinci, A. (2021). A new smoothed seismicity approach to include aftershocks and foreshocks in spatial earthquake forecasting: application to the global Mw≥ 5.5 seismicity. Applied Sciences, 11(22), 10899.
- [9] Rahman, M. M., Bai, L., Li, H., & Liu, C. (2026). Probabilistic Seismic Hazard Map for Bangladesh Including the Smoothed Background Seismicity and Local Site Effects. Pure and Applied Geophysics, 183(2), 335-361.
- [10] Chan, C. H., Ma, K. F., Shyu, J. B. H., Lee, Y. T., Wang, Y. J., Gao, J. C., ... & Rau, R. J. (2020). Probabilistic seismic hazard assessment for Taiwan: TEM PSHA2020. Earthquake Spectra, 36(1_suppl), 137-159.
- [11] Hanks, T. C., & Kanamori, H. (1979). A moment magnitude scale. Journal of Geophysical Research: Solid Earth, 84(B5), 2348-2350.
- [12] Reasenberg, P. (1985). Second‐order moment of central California seismicity, 1969–1982. Journal of Geophysical Research: Solid Earth, 90(B7), 5479-5495.
- [13] Hanafi, R., Dewi, I. K., & Ngatijo, N. (2024). Pemetaan Magnitude Of Completeness (Mc) Untuk Gempa Di Wilayah Bengkulu. JGE (Jurnal Geofisika Eksplorasi), 10(2), 121-138.
- [14] Irsyam, M., Cummins, P. R., Asrurifak, M., Faizal, L., Natawidjaja, D. H., Widiyantoro, S., ... & Syahbana, A. J. (2020). Development of the 2017 national seismic hazard maps of Indonesia. Earthquake Spectra, 36(1_suppl), 112-136
- [15] Lewerissa, R., Rumakey, R., Syakur, Y. A., & Lapono, L. (2021). Completeness magnitude (M c) and b-value characteristics as important parameters for future seismic hazard assessment in the West Papua province Indonesia. Arabian Journal of Geosciences, 14(23), 2588.
- [16] Ritz, V. A., Rinaldi, A. P., & Wiemer, S. (2022). Transient evolution of the relative size distribution of earthquakes as a risk indicator for induced seismicity. Communications Earth & Environment, 3(1), 249.
- [17] Geffers, G. M., Main, I. G., & Naylor, M. (2023). Accuracy and precision of frequency–size distribution scaling parameters as a function of dynamic range of observations: example of the Gutenberg–Richter law b-value for earthquakes. Geophysical Journal International, 232(3), 2080-2086.
- [18] Kijko, A., & Smit, A. (2017). Estimation of the frequency–magnitude Gutenberg–Richter b‐value without making assumptions on levels of completeness. Seismological Research Letters, 88(2A), 311-318.
- [19] Taroni, M., & Akinci, A. (2021). A new smoothed seismicity approach to include aftershocks and foreshocks in spatial earthquake forecasting: application to the global Mw≥ 5.5 seismicity. Applied Sciences, 11(22), 10899.
- [20] Lapajne, J. K., Motnikar, B. S., Zabukovec, B., & Zupancic, P. (1997). Spatially smoothed seismicity modeling of seismic hazard in Slovenia. Journal of seismology, 1(1), 73-85.
- [21] Molchan, G., & Keilis-Borok, V. (2008). Earthquake prediction: probabilistic aspect. Geophysical Journal International, 173(3), 1012-1017.
- [22] Zechar, J. D., & Jordan, T. H. (2008). Testing alarm-based earthquake predictions. Geophysical Journal International, 172(2), 715-724.
- [23] Biondini, E., D'Orazio, F., Lolli, B., & Gasperini, P. (2025). Pseudo-prospective earthquake forecasting experiment in Italy based on temporal variation of the b-value of the Gutenberg–Richter law. Geophysical Journal International, 240(3), 1755-1772.
- [24] Moschetti, M. P., Aagaard, B. T., Ahdi, S. K., Altekruse, J., Boyd, O. S., Frankel, A. D., ... & Withers, K. B. (2024). The 2023 US national seismic hazard model: Ground-motion characterization for the conterminous United States. Earthquake Spectra, 40(2), 1158-1190.
- [25] Kaviris, G., Zymvragakis, A., Bonatis, P., Kapetanidis, V., Spingos, I., Mavroulis, S., ... & Voulgaris, N. (2023). A logic-tree approach for probabilistic seismic hazard assessment in the administrative region of Attica (Greece). Applied Sciences, 13(13), 7553.
- [26] Boore, D. M., Stewart, J. P., Seyhan, E., & Atkinson, G. M. (2014). NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes. Earthquake Spectra, 30(3), 1057-1085.
- [27] Campbell, K. W., & Bozorgnia, Y. (2014). NGA-West2 ground motion model for the average horizontal components of PGA, PGV, and 5% damped linear acceleration response spectra. Earthquake Spectra, 30(3), 1087-1115.
- [28] Chiou, B. S. J., & Youngs, R. R. (2014). Update of the Chiou and Youngs NGA model for the average horizontal component of peak ground motion and response spectra. Earthquake Spectra, 30(3), 1117-1153.
- [29] Pusat Studi Gempa Nasional (Indonesia). (2017). Peta sumber dan bahaya gempa Indonesia tahun 2017. Pusat Penelitian dan Pengembangan Perumahan dan Permukiman, Badan Penelitian dan Pengembangan, Kementerian Pekerjaan Umum.
- [30] Rahman, A. U., Najam, F. A., Zaman, S., Rasheed, A., & Rana, I. A. (2021). An updated probabilistic seismic hazard assessment (PSHA) for Pakistan. Bulletin of Earthquake Engineering, 19(4), 1625-1662.
- [31] Danciu, L., Giardini, D., Weatherill, G., Basili, R., Nandan, S., Rovida, A., ... & Wiemer, S. (2024). The 2020 European seismic hazard model: overview and results. Natural Hazards and Earth System Sciences, 24(9), 3049-3073.
- [32] Woo, G. (1996). Kernel estimation methods for seismic hazard area source modeling. Bulletin of the Seismological Society of America, 86(2), 353-362.
- [33] Petersen, M. D., Shumway, A. M., Powers, P. M., Field, E. H., Moschetti, M. P., Jaiswal, K. S., ... & Witter, R. C. (2024). The 2023 US 50-state national seismic hazard model: Overview and implications. Earthquake Spectra, 40(1), 5-88.
