Penerapan Model Epidemic Type Aftershock Sequence (ETAS) pada Data Gempa Bumi Jawa Barat

Authors

  • Sofi Sopiah Statistika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Islam Bandung
  • Sutawanir Darwis Statistika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Islam Bandung

DOI:

https://doi.org/10.29313/bcss.v4i2.14956

Keywords:

Model ETAS, Fungsi Intensitas, Likelihood maksimum

Abstract

Abstract. Earthquake intensity mapping is an important step in mitigating the risk of natural disasters, especially earthquakes. Earthquake activity is still being studied from both seismological and stochastic aspects. One of the stochastic processes that can explain earthquake activity is the point process. The point process is a stochastic process that can explain natural phenomena that are random in both space and time. One of the models in the point process is the epidemic type aftershock sequence (ETAS) model. The ETAS model is expressed by a conditional intensity function that is useful for knowing the chances of earthquake occurrence. Therefore, the purpose of this study is to apply the ETAS model to earthquake data in West Java. Data obtained from the United State Geological Survey (USGS) was declustered with the Gardner Knopoff method and there were 49 clusters and took 1 cluster with the largest earthquake data, namely the earthquake on September 02, 2009 in Tasikmalaya. The maximum likelihood estimation method is used to obtain parameter estimates of the ETAS model. The null model and full model are used to show the magnitude of the exponential distribution, while the full model shows the magnitude of the gamma distribution. With the maximum likelihood method, the parameter estimates of the basic intensity function of the null model and full model were obtained. The earthquake in Tasikmalaya in 2009 obtained the best model, the full model. The results of the parameter estimation are the base seismicity rate of 0.04402, the aftershock productivity of 4.673611e-9, the efficiency of an earthquake with a certain magnitude to produce aftershocks of 59.4750, the time scale of aftershock decay rate of 0.0083, and the aftershock decay rate of 1.3825.

Abstrak. Pemetaan intensitas gempa merupakan langkah penting dalam upaya mitigasi risiko bencana alam, khususnya gempa bumi. Aktivitas gempa bumi masih terus dikaji baik dari aspek seismologi maupun aspek stokastik. Salah satu proses stokastik yang dapat menjelaskan aktivitas gempa bumi adalah proses titik. Proses titik merupakan suatu proses stokastik yang dapat menjelaskan fenomena alam dimana sifatnya acak baik dalam ruang maupun waktu. Salah satu model dalam proses titik adalah model epidemic type aftershock sequence (ETAS). Model ETAS dinyatakan dengan fungsi intensitas bersyarat yang berguna untuk mengetahui peluang kemunculan terjadinya gempa bumi. Oleh karena itu tujuan penelitian ini adalah menerapkan model ETAS pada data gempa bumi di Jawa Barat. Data diperoleh dari United State Geological Survey (USGS) di decslustering dengan metode Gardner Knopoff  terdapat 49 cluster dan mengambil 1 cluster data gempa bumi yang terbesar yaitu gempa pada tanggal 02 September 2009 yaitu di Tasikmalaya. Metode estimasi likelihood maksimum digunakan untuk memperoleh estimasi parameter model ETAS. Digunakan null model dan full model untuk menunjukkan magnitudo berdistribusi eksponensial, sedangkan full model menunjukkan magnitudo berdistribusi gamma. Dengan metode likelihood maksimum diperoleh estimasi parameter fungsi intensitas dasar dari null model dan full model. Gempa di Tasikmalaya pada tahun 2009 memperoleh model terbaik yaitu full model. Hasil estimasi parameter tersebut yaitu laju kegempaan dasar sebesar 0.04402, produktivitas gempa susulan sebesar 4.673611e-9, efisiensi gempa bumi dengan magnitudo tertentu menghasilkan gempa susulan sebesar 59.4750, skala waktu laju peluruhan gempa susulan sebesar 0.0083, dan laju peluruhan gempa susulan sebesar 1.3825.

References

Console, R., Murru, M., & Falcone, G. (2010). Retrospective forecasting of M ≥ 4.0 earthquakes in New Zealand. Pure and Applied Geophysics, 167(6), 693–707.

Daley, D. J., & Vere-Jones, D. (2007). An Introduction to the Theory of Point Processes: Volume I: Elementary Theory and Methods. Springer.

Fisher, R. A. (1922). On the mathematical foundations of theoretical statistics. Philosophical transactions of the Royal Society of London. Series A, containing papers of a mathematical or physical character, 222(594-604), 309-368.

Harte, D. (2010). PtProcess: An R package for modelling marked point processes indexed by time. Journal of Statistical Software, 35, 1-32.

Mutiah, S., Pratiwi, H., & Handajani, S. S. (2019). Penerapan Model Epidemic Type Aftershock Sequence (ETAS) pada Data Gempa Bumi Sulawesi dan Jawa.

Ogata, Y. (1988). Statistical models for earthquake occurrences and residual analysis for point processes. Journal of the American Statistical association, 83(401), 9-27.

Ogata, Y. (1998). Space-time point-process models for earthquake occurrences. Annals of the Institute of Statistical Mathematics, 50,379- 402.

Ogata, Y. (1999). Seismicity analysis through point-process modeling: A review. Seismicity patterns, their statistical significance and physical meaning, 471-507.

Pusat Studi Gempa Nasional. (2017). Peta Sumber dan Bahaya Gempa Indonesia Jakarta: Pusat Penelitian dan Pengembangan Perumahan dan Pemukiman.

Supartoyo, dan Surono. (2009). Kegempaan di Wilayah Jawa Barat dan Kejadian Gempa bumi Jawa Barat Selatan Tanggal 2 September 2009, Jurnal Gunungapi dan Mitigasi Bencana Geologi, 1(2).

Schoenberg, F. P., Brillinger, D. R., & Guttorp, P. (2012). Point Processes, Spatial‐Temporal. Encyclopedia of Environmetrics.

Utsu, T. (1969). Instructions for use Aftershocks and Earthquake Statistics (I). Journal of the Faculty of Science, Hokkaido University, 3(3), 129–195.

Wijaya, A. R. (2023). Model ETAS Spatio-Temporal pada Analisis Pemetaan Intensitas Kegempaan di Wilayah Sumatera. Jambura Journal of Mathematics, 5(1), 179–188.

Zhuang, J. (2000). Statistical modelling of seismicity patterns before and after the 1990 Oct 5 Cape Palliser earthquake, New Zealand. New Zealand Journal of Geology and Geophysics, 43(3), 447–460.

Ikhsanudin, Istiyono, E., & Syamsuddin, S. (2023). Efek Besaran Simpangan Baku Sampel terhadap Nilai Residu dalam Analisis Regresi Berganda Tiga Variabel Bebas. Statistika, 23(2), 88–95. https://doi.org/10.29313/statistika.v23i2.1661

Karimuse, W. Y., Nohe, D. A., & Siringoringo, M. (n.d.). Pendekatan Regresi Nonparametrik Kernel pada Data IHSG Periode Januari 2020-Desember 2021 (Vol. 23, Issue 1). https://www.idx.co.id/

Meidianingsih, Q., Wardani, D. E., Salsabila, E., Nafisah, L., & Mutia, A. N. (n.d.). Perbandingan Performa Metode Berbasis Support Vector Machine untuk Penanganan Klasifikasi Multi Kelas Tidak Seimbang (Vol. 23, Issue 1).

Downloads

Published

2024-08-14