## IDENTIFIKASI BAHAYA PENGUATAN GEMPA WILAYAH JAWA TAHUN 1974-2020

## Iis Ratna Sari, Haqqie Ramadhan, Supardi, Syamsurijal Rasimeng, Sunanda Yogi
Jurusan Teknik Geofisika, Fakultas Teknik, Universitas Lampung e-mail : iisratna790@gmail.com
Abstrak. Penelitian  mengenai  identifikasi  bahaya  penguatan  gempa  sangatlah  penting,  mengingat dampak  gempabumi  yang  cukup  merugikan.  Data  yang  digunakan  dalam  penelitian  ini,  dianalisis menggunakan program matlab. Adapun tujuannya yaitu untuk mengetahui persebaran gempa, periode munculnya  gempa,  dan  seberapa  besar  gempa  yang  terjadi  serta  frekuensinya.  Lokasi  penelitian berada pada koordinat  5,6 o - 8,5 o LS hingga 105,9 o - 114,6 o BT dengan rentang waktu antara 1 Januari 1974-30 April 2020. Berdasarkan data yang diperoleh menunjukkan bahwa magnitude gempa tertinggi mempunyai  nilai  7,5  Mw  yang  terjadi  tanggal  08  Agustus  2007  dengan  kedalaman  280  km,  dengan magnitude  gampa  terendah  adalah  5  Mw.  Sedangkan,  nilai  b-value  sebesar  0,651.  Hasil  analisa  data gempa dengan program matlab digunakan untuk mengungkapkan  adanya periodogram, dan diperoleh nilai  frekuensi  dominan  sebesar  0,17266.  Selain  itu,  pada  penelitian  ini  juga  dihasilkan  periode munculnya gempa dalam skala bulan, dimana dari hasil analisa yang dilakukan diperoleh nilai periode pelepasan energi sebesar 0,20863 bulan dengan energi sebesar 2,08 x 1034  joule.
Kata Kunci: Frekuensi Dominan; Gempa bumi; Nilai B-value; Matlab; Surfer
Abstract. Research  on  the  identification  of  the  hazards  of  earthquake  amplification  is  very  important, given the adverse impact of earthquakes. The data used in this study were analyzed using the surfer and matlab  programs.  The  aim  is  to  determine  the  distribution  of  earthquakes,  the  period  when  the earthquake occurred, and how big the earthquake occurred and its frequency. The research location is at coordinates 5,6 o - 8,5 o LS to 105,9 o - 114,6 o BT with a time span between January 1, 1974 - April 30, 2020. Based on the data obtained, it shows that the highest earthquake magnitude has a value of 7.5 Mw which occurred  on  August  8,  2007  with  a  depth  of  280  km,  with  the  lowest  earthquake  magnitude  is  5  Mw. Meanwhile,  the  b-value  is  0.651.  The  results  of  earthquake  data  analysis  using  Matlab  to  reveal  the existence  of  a  periodogram,  and  obtained  a  dominant  frequency  value  of  0.17266.  In  addition,  this research also generates the period of the occurrence of earthquakes on the scale of the month, where the results of the analysis obtained the value of the energy release period of 0.20863 months with an energy of 2,08 x 1034   joules.
Keywords: Dominant Frequency; Earthquake; B-value; Matlab; Surfer

## PENDAHULUAN
Indonesia  berada  pada  kawasan  seismik  aktif dengan frekuensi kegempaan yang sangat tinggi, hal ini  disebabkan  oleh  kondisi  tatanan  tektonik  yang begitu  kompleks,  diantaranya  pertemuan  lempeng pasifik  yang  berada  disepanjang  Maluku,  dan pertemuan  lempeng  Indo-Australia  serta  Eurasia yang  membentang  sepanjang  barat  sumatera, selatan Jawa, Bali, dan Nusa Tenggara (Rahmadi Tri Wardana, 2008).
Hal itu lah, yang menyebabkan beberapa wilayah di  Indonesia  menjadi  rentan  akan  terjadinya bencana  gempabumi,  salah  satunya  yaitu  Pulau Jawa.  Aktivitas  tektonik  dari  pertemuan  antar lempeng ini menyebabkan terbentuknya berbagai
struktur  geologi  di  Pulau  Jawa  yaitu  berupa patahan/sesar,  seperti:  Sesar  Cimandiri,  Sesar Lembang,  Sesar  Baribis,  Sesar  Opak,  dan  lain sebagainya (PUSGEN, 2017).
Berdasarkan  data  BMKG,  gempa  &  tsunami berulang  kali  terjadi  di  bagian  Selatan  Jawa  yang mengakibatkan  kerusakan  dan  juga  korban  jiwa. Contohnya, yaitu  gempa  dengan  kekuatan  6.8  Mw pada  tanggal  2  September  2009  yang  melanda bagian  selatan  Jawa  Barat  (Hanifa  dkk.,  2014). Kemudian,  peristiwa  slow  slip  sesar  dengan kekuatan gempabumi sebesar 7.8 Mw pada tanggal 17 Juli 2006 yang menimbulkan tsunami dahsyat di Pangandaran  (Gunawan  dkk.,  2016)  dan  gempa besar lainnya.
Berdasarkan  kondisi  tersebut,  maka  identifikasi bahaya  penguatan  gempa  sangatlah  penting khususnya  untuk  wilayah  yang  memiliki  resiko tinggi terjadinya gempa bumi. Identifikasi penguatan  gempa  dapat  dipelajari  melalui  analisis intensitas  gempa.  Data  yang  digunakan  dalam penelitian  ini,  akan  diolah  menggunakan  program matlab.  Adapun  tujuan  dari  penelitian  ini  yaitu untuk  mengetahui  persebaran  gempa,  periode munculnya  gempa,  dan  seberapa  besar  gempa yang terjadi serta frekuensinya.
Gambar 1. Kondisi tektonik Indonesia (BMKG, 2018)
Gambar 2. Patahan/sesar di Pula Jawa (Ryanto dkk., 2019)
Gambar 3. Peta lokasi penelitian

## METODOLOGI
Penelitian  ini menggunakan  data sekunder  yang diunduh  dari  website  (USGS,  2020)  pada  rentang waktu  antara  1  Januari  1974  hingga  30  April  2020 dengan magnitude 5-10 Mw, dan kedalaman gempa 0-900  km.  Batasan  wilayah  penelitian  yaitu  pada 5,6 o -8.5 o LS  hingga  105,9 o -114,6 o BT.  Jumlah  data yang  digunakan  sebanyak  371  event  gempa.  Alat yang  digunakan  terdiri  dari  laptop,  Microsoft  Excel 2007,  Software  ArcGis,  Notepad++,  dan  Software Matlab R2007b.
Pengolahan data dalam penelitian ini yaitu: a. Pengumpulan  data  gempa berupa data latitude , longitude , origin  time ,    depih,  moment  seismik (Mo),  magnitude  gelombang  badan  (Mb),  dan magnitude gelombang permukaan (Ms).
b. Konversi  jenis  magnitude  Mb  dan  Ms  menjadi
Mw  untuk jenis magnitude.
c. Membuat klasifikasi gempa berdasarkan magnitude beserta log frekuensinya.
d. Menghitung  nilai  b-value  melalui  grafik  dari trendline data.
e. Membuat  peta  sebaran  gempa  pada  daerah penelitian  dengan  menggunakan  software ArcGIS.
Analisa  menggunakan  program  matlab  akan diperoleh  grafik time  series . Time  series adalah serangkaian  nilai-nilai  variabel  yang  disusun berdasarkan  waktu.  Analisis time  series akan mempelajari  pola  gerakan  nilai-nilai  variabel  pada suatu  interval  waktu  yang  diatur.  Kemudian hasilnya  diperoleh  sebuah  grafik  hubungan  antara Julian month dan magnitude (Suwandana, N. 2013). Adapun  langkah  pengolahan  data  menggunakan program matlab yaitu :
a. Membuat Magnitude Time Series Data  yang  diperlukan  berupa  data month dan data magnitude (Mw).
b. Membuat Energi Time Series
Data yang diperlukan yaitu data month dan data energi  gempa.  Energi  gempa  dapat  dihitung dengan menggunakan rumus:
LogEs = 11.8 + 1.5M , (1) c. Membuat Baseline Correction
Data  yang  diperlukan  dalam  pembuatan baseline  corection terdiri  dari  data month ,  data energi rata-rata dan data perubahan energi.
d. Membuat Periodogram Grafik periodogram dibuat dengan
pengembangan  script  matlab  untuk  FFT.  Pada
Gambar 6. Peta sebaran hiposenter gempa wilayah Jawa
pembuatan perodogram ini, data yang digunakan  yaitu  waktu  kejadian  gempa  yang dikonversi ke julian month , dan nilai magnitudo yang telah dikoreksi baseline .
e. Membuat Grafik Power Vs Periode Grafik ini juga dibuat menggunakan pengembangan  script  matlab  untuk  FFT.    Data yang  digunakan  yaitu  data  waktu  kejadian gempa  yang  dikonversi  ke julian  month ,  dan data  perubahan  energi  atau  energi-energi  rata- rata  yang  telah  dikoreksi baseline (Suwandana, 2013).

## HASIL DAN PEMBAHASAN
Data  gempa  yang  diperoleh  dari  katalog  USGS akan diklasifikasikan berdasarkan range magnitude dan  frekuensinya.  Tabel  1  menunjukkan  bahwa gempa dengan kekuatan 5 Mw sering mengguncang  wilayah  Pulau  Jawa  dengan kedalaman  gempa  yang  bervariasi.  Magnitude paling  tinggi  tercatat  sebesar  7,5  Mw  dengan kedalaman  280  km  yang  terjadi  tepatnya  pada tanggal  08 Agustus 2007.
Tabel 1. Klasifikasi gempa No Mw Frekuensi Log Frekuensi 1 5-5.9 343 2.53529412 2 6-6.9 28 1.447158031 3 7-7.9 2 0.301029996 4 8-8.9 0 - 5 9-9.9 0 -
Pada grafik (Gambar 5) dibawah ini, diperoleh: B Value - Data y = -0,651x + 2,812 (2) dengan
m = b-value
sehingga  nilai  b-value  wilayah  Pulau  Jawa  adalah 0,651.
Beberapa  ahli  mengatakan  bahwa  nilai  konstan b-value yaitu sekitar 1. Apabila nilai b-value rendah maka  menunjukkan  bahwa  daerah  tersebut memiliki  tingkat  kerapuhan  batuan  yang  rendah, daya tahan terhadap batuan besar, dan nilai b-value yang  rendah  berhubungan  dengan  tingkat  stress yang tinggi (Hilmi dkk., 2019).
Gambar 5. Grafik nilai b-value
Sedangkan,  jika  nilai  b-value  tinggi  maka menunjukkan  bahwa  daerah  tersebut  memiliki tingkat  kerapuhan  batuan  yang  tinggi,  daya  tahan
terhadap  batuan  kecil,  dan  nilai  b-value  tinggi berhubungan tingkat stress yang rendah (Hilmi dkk., 2019).
Pada peta gambar 6, menunjukan bahwa tingkat seismisitas atau intensitas kegempaan di Pulau Jawa sangatlah  tinggi.  Terlihat  bahwa  gempabumi  yang terjadi sangat dominan pada bagian wilayah selatan Pulau  Jawa  dan  kebanyakan  gempa  berasal  dari lautan.  Hal  ini  karena  bagian  selatan  pulau  jawa terdapat  lempeng  besar,  dan  pergerakan  lempeng inilah  yang  sering  menjadi  pemicu  gempabumi  di wilayah jawa.
Gambar 7. (a) Grafik magnitudo time series, (b) Grafik periodogram, (c) Grafik power vs periode
Dari  grafik time  series gambar  7a,  menunjukan bahwa  gempa  yang  sering  terjadi  di  daerah penelitian didominasi oleh gempa dengan kekuatan menengah yaitu dari 5 sampai 6 Mw.
Pada grafik periodogram (gambar 7b) diperoleh nilai  frekuensi  dominan  di  wilayah  Pulau  Jawa sebesar 0,17266. Selain itu, pada penelitian ini juga dihasilkan  periode  munculnya  gempabumi  dalam skala  bulan  ( month )  seperti  pada  gambar  7c.  Dari hasil  analisa  yang  dilakukan  maka  diperoleh  nilai periode  pelepasan  energi  yaitu  sebesar  0,20863 bulan dengan energi sebesar 2,08 x 10 34 joule.
Frekuensi  dan  Periode  muncul  nya  gempa  di Pulau  Jawa  sangat  sering  terjadi,    sehingga  perlu diwaspadai terutama pada daerah selatan wilayah Jawa.  Kesiapsiagaan  menghadapi  bencana  harus dibangun  dalam  setiap  kelompok  masyarakat untuk  meminimalisir  kerusakan  dan  juga  korban jiwa.  Kesiapsiagaan dapat dimulai dengan memetakan  daerah-daerah  yang  rawan  akan bencana, memperkuat infrastruktur dengan membuat  bangunan  dan  rumah  tahan  gempa, serta memberikan pemahaman mengenai kebencanaan kepada seluruh lapisan masyarakat.

## PENUTUP Simpulan dan Saran
Berdasarkan hasil penelitian, disimpulkan bahwa:  Banyaknya  aktivitas  seismik  terutama  di  bagian selatan  disebabkan  oleh  adanya  zona  subduksi yang  membentang  di  sepanjang  daerah  selatan pulau jawa.
 Gempabumi yang terjadi didominasi oleh  gempa dengan  kedalaman  yang  relatif  dangkal  yaitu pada  0-70  km  dan  berkekuatan  sedang  yaitu  5 hingga  6  Mw.  Gempabumi  dengan  kedalaman dangkal ini akan menyebabkan tingginya potensi kerusakan.
 Energi Gempa dari tahun 1974-2020 mempunyai nilai rata-rata energi  Sebesar  1.48E+21 joule.
 Frekuensi  Dominan  yang  dihasilkan  dari penelitian ini sebesar  0,17266.
 Periode  yang  dihasilkan  yaitu  sebesar  0,20863 sehingga  rangkaian  gempa  akan  terjadi  lagi setiap 0,20863 bulan.
(a)
(b) (c)
Ketika  proses  pengumpulan  dan  pemilihan  data, ketelitian  sangat  diperlukan  agar  pada  saat pengolahan  menggunakan  matlab,    kurva  yang dihasilkan tidak error.

## Ucapan Terima Kasih
Penulis  mengucapkan  terimakasih  kepada  seluruh pihak yang telah membantu sehingga penelitian ini dapat diselesaikan.

## DAFTAR PUSTAKA
BMKG  (2018), Peta  Tektonik  Indonesia .  Diambil  23 Desember 2020, dari http://balai3.denpasar.bmkg.go.id/. Gunawan, E., Meilano, I., Abidin, H.Z., Hanifa, N.R., dan Susilo (2016), "Investigation of the Best Coseismic  Fault  Model  of  the  2006  Java Tsunami Earthquake Based on Mechanisms of  Postseismic  Deformation", Journal  of
Asian  Earth  Sciences ,  Vol.117,  hal.  64–72. http://doi.org/10.1016/j.jseaes.2015.12.00
3.
Hanifa, N.R., Sagiya, T., Kimata, F., Efendi, J., Abidin, H.Z.  dan  Meilano,  I.  (2014),  "Interplate
Coupling Model off the Southwestern Coast of  Java,  Indonesia,  Based  on  Continuous GPS  Data  in  2008–2010", Earth  and Planetary  Science  Letters ,  Vol.401,  hal. 159–171. http://doi.org/10.1016/j.epsl.2014.06.010.
Hilmi,  I.L.,  Sutrisno,  S.  dan  Sunarya,  D.  (2019),
"Analisis  Seismisitas  Berdasarkan  Data Gempa
Bumi Periode
1958-2018 Menggunakan b-Value Pada Daerah Selatan Jawa  Barat  Dan  Banten", Al-Fiziya:  Journal of Materials Science, Geophysics, Instrumentation  and  Theoretical  Physics , Vol.2, No.1, hal. 10–16. http://doi.org/10.15408/fiziya.v2i1.10482. PUSGEN  (2017), Buku  Peta  Sumber  Dan  Bahaya
Gempa  Indonesia  Tahun  2017 .  Diambil  23 Desember 2020, dari
https://www.pu.go.id/berita/view/15309/b uku-peta-sumber-dan-bahaya-gempa- indonesia-tahun-2017. Rahmadi  Tri  Wardana  (2008), Analisa  Statistik Tingkat  Keaktifan  Gempa  Bumi  di  Daerah Jawa  Timur  dan  sekitarnya  Tahun  1973-
2007 STMKG, Diambil dari
http://perpus.stmkg.ac.id/book- detail.php?id=5783.
Ryanto, T.A., Suntoko, H. dan Setiaji, A.B.W. (2019), "Pendugaan  Awal  Patahan  di  Pulau  Jawa Menggunakan  Anomali  Gravitasi  dan
Riwayat Kegempaan", EKSPLORIUM , Vol.40, No.1, hal. 43–52. http://doi.org/10.17146/eksplorium.2019.4 0.1.5470. Suwandana,  N.  (2013), Identifikasi  Gempa  dengan Matlab Diambil dari https://id.scribd.com/document/17133607 7/Identifikasi-Gempa-dengan-Matlab.
USGS  (2020), Data  Gempa .  Diambil  23  Desember 2020,  dari  https://www.usgs.gov/natural- hazards/earthquake-hazards/earthquakes.
-------------------------

## Lampiran 1
Gambar 10. Peta sebaran hiposenter gempa pada kedalama kurang dari 30 km.
Gambar 11. Peta sebaran hiposenter gempa pada kedalaman 30-60 km.
Gambar 12. Peta sebaran hiposenter gempa pada kedalaman 60-100 km.
Gambar 13. Peta sebaran hiposenter gempa pada kedalaman lebih dari 100 km.
Gambar 14. Grafik time series energi gempa
Gambar 15. Grafik perubahan energi terhadap waktu ( month )

## Lampiran 2
Tabel 2 . Data gempa yang digunakan dalam penelitian diunduh dari website USGS
time latitude longitude depth mag magType nst gap dmin rms net id 2020-03-10T10:18:04.505Z -6.7852 106.6813 10 5 mb 86 1.173 1.15 us us60008d3j 2020-02-05T18:12:37.734Z -6.0817 113.0778 592.42 6.2 mww 31 5.479 1.04 us us70007j6z 2019-09-19T07:32:58.793Z -6.0534 111.8968 636.41 5.6 mb 25 1.752 0.69 us us60005kuc 2019-09-19T07:32:01.582Z -6.0586 111.8877 617.56 5.9 mww 19 1.742 0.87 us us60005ku0 2019-09-19T07:06:33.291Z -6.0708 111.8422 610 6.1 mww 19 1.698 0.97 us us60005kta 2019-07-28T14:25:02.625Z -6.9777 106.1423 58.47 5.1 mb 102 1.469 0.82 us us60004wrw 2019-07-24T01:29:12.731Z -8.5946 114.3134 102.51 5 mb 54 0.201 1.05 us us60004ugn 2019-06-29T06:47:04.303Z -8.1929 108.9633 78.89 5.1 mb 122 1.569 0.84 us us7000477c 2019-06-09T09:32:25.075Z -8.1979 108.9471 70.44 5.4 mww 45 1.586 0.89 us us70003xf0 2019-02-13T23:41:54.620Z -6.7647 105.803 66.89 5 mb 40 1.799 1.08 us usd0009e9l 2019-01-08T09:54:47.990Z -7.8778 106.5655 35 5.2 mb 46 1.476 1.11 us us2000j1c6 2018-11-29T20:42:05.420Z -8.7025 109.7932 89.24 5 mb 114 1.066 1.2 us us1000hy1r 2018-10-10T18:44:55.280Z -7.4525 114.4553 9 6 mww 27 1.054 1.15 us us1000ha6q 2018-03-19T11:36:27.050Z -7.41 105.9404 48.55 5 mww 58 3.064 1.2 us us1000d6de 2018-01-23T06:34:54.980Z -7.0924 105.9632 48.19 5.9 mww 27 1.66 1.45 us us2000cmwz 2017-12-16T00:22:30.840Z -7.8842 106.8195 44 5.2 mww 27 1.316 0.99 us us2000c551 2017-12-15T16:47:58.230Z -7.4921 108.1743 90 6.5 mww 28 2.365 1.07 us us2000c4v8 2017-12-12T02:23:10.090Z -7.1269 107.0917 100.75 5 mb 89 3.324 0.92 us us2000c35s 2017-09-25T18:10:43.200Z -5.9298 106.3876 150.05 5.1 mb 35 3.446 1.15 us us2000aufz 2017-09-20T23:59:24.900Z -6.1113 113.024 588 5.7 mww 13 2.597 1.17 us us2000as7r 2017-08-03T16:14:23.470Z -8.1608 108.9879 63.18 5 mww 27 1.54 1.18 us us2000a467 2017-07-08T22:48:31.650Z -7.4094 105.9824 8.54 5.2 mb 25 3.069 1.15 us us1000989m 2017-06-11T23:15:06.150Z -8.3206 106.259 7 5.7 mww 46 2.004 1.46 us us20009lwc 2017-04-23T18:01:13.310Z -7.7404 108.0192 72.22 5.2 mb 28 2.486 1.24 us us10008k20 2017-02-15T18:53:38.570Z -7.5161 107.3878 77.33 5 mb 40 3.063 1.21 us us20008k1a 2017-02-08T08:34:18.360Z -7.5557 106.9498 74.88 5 mb 64 0.982 1.32 us us20008iff 2017-02-07T17:21:47.160Z -8.4287 113.1306 140.72 5.2 mb 49 1.011 1.2 us us20008i9k 2017-01-02T21:02:25.190Z -8.7459 108.7974 38.14 5.2 mb 53 1.897 0.86 us us10007pbq 2016-04-06T14:45:29.620Z -8.2036 107.3857 29 6.1 mww 31 1.39 1.44 us us20005fv5 2016-03-02T06:09:12.940Z -8.6191 112.5355 78.32 5.1 mb 31 1.605 1.22 us us10004u04 2016-02-25T23:35:19.300Z -8.7463 111.5091 57.44 5.1 mb 40 1.28 1.36 us us10004srp 2015-12-19T15:03:36.010Z -8.3684 107.0225 27 5.7 mww 87 2.496 1.4 us us100048d4 2015-11-28T14:47:18.970Z -7.0548 105.9999 51.34 5.4 mww 47 3.421 1.42 us us100041wc 2015-10-10T09:04:01.020Z -8.0486 107.4009 35 5.2 mb 35 2.969 1.13 us us10003mnt 2015-09-04T20:08:31.580Z -8.043 107.2248 37 5.4 mww 50 0.759 1.63 us us20003ghh 2015-07-24T21:44:38.920Z -8.2477 108.9163 48 5.5 mwb 17 1.625 1.33 us us20002zzg 2014-11-03T08:56:30.790Z -7.4091 105.9668 25.05 5 mb 52 1.735 0.74 us usc000stmf 2014-10-17T09:26:00.240Z -6.4963 107.1731 151.61 5.2 mb 19 0.545 1.04 us usb000snbi 2014-06-22T06:38:51.440Z -7.7615 108.0537 72.02 5.1 mb 28 0.312 0.85 us usc000rjxl 2014-01-27T16:14:00.800Z -8.0452 109.2927 87.85 5 mb 23 1.225 0.68 us usb000m8w7 2014-01-25T05:25:16.000Z -8.0009 109.3003 84.51 5.2 mb 44 1.214 0.5 us usc000mkb3 2014-01-25T05:14:18.510Z -7.9855 109.2653 66 6.1 mww 20 1.248 1.45 us usb000m7wd 2013-08-08T10:45:58.200Z -8.6243 110.9786 9.45 5.3 mb 25 0.838 1.43 us usb000iyfk 2013-04-08T18:53:44.370Z -7.091 106.072 56.4 5.3 mwb 410 23.3 0.82 us usb000g2b9 2012-11-10T06:05:27.190Z -7.373 106.666 49.5 5.1 mb 411 17.2 1.13 us usp000jva5 2012-11-01T14:12:01.820Z -6.762 107.529 142.1 5.7 mww 271 23.8 1.08 us usp000juw5 2012-06-12T05:59:38.690Z -5.677 105.494 163.9 5.5 mb 444 15.7 0.93 us usp000jmnp 2012-06-04T11:18:13.450Z -7.692 106.371 50 5.9 mww 542 26.6 0.88 us usp000jmap 2011-11-24T03:55:36.540Z -8.644 109.888 90 5.1 mb 98 101.8 0.86 us usp000jbaz 2011-11-06T04:41:57.140Z -6.132 106.301 163.8 5.2 mb 213 24 0.98 us usp000jadv 2011-08-08T02:35:04.180Z -5.824 111.294 510.3 5 mb 71 31.2 0.96 us usp000j65j 2011-07-24T08:23:36.380Z -7.197 106.458 59.1 5 mwc 121 28.7 1.02 us usp000j5dj 2011-07-17T10:59:09.690Z -6.64 106.489 94.8 5 mb 171 25.9 0.95 us usp000j52k 2011-05-28T20:28:55.840Z -8.756 108.754 45.3 5.1 mb 161 25.8 1.17 us usp000j2ce 2011-04-26T06:39:36.520Z -8.081 108.477 95.4 5.6 mww 215 27 1.03 us usp000j0m1 2011-04-05T04:34:46.260Z -6.93 106.22 82.7 5.2 mwc 155 23 1 us usp000hza4 2011-02-26T06:29:35.440Z -6.556 107.187 187.4 5 mb 68 34.7 1.12 us usp000hv3a 2010-12-21T03:59:38.350Z -8.7 111.197 54.6 5.6 mwb 310 21.4 1.01 us usp000hr88 2010-11-09T12:38:59.620Z -7.865 107.206 56.8 5.1 mwc 174 17.8 1.23 us usp000hpas 2010-11-09T07:03:30.360Z -8.634 110.026 63 5.4 mwb 302 22.9 1.14 us usp000hpad 2010-08-16T22:19:17.250Z -8.546 112.53 73.6 5 mb 138 23.3 1.13 us usp000hjan 2010-08-11T19:10:21.430Z -7.495 107.058 37.4 5.3 mwc 317 18.4 1.01 us usp000hhyz 2010-08-09T03:34:15.980Z -6.082 112.995 598.5 5.3 mwc 320 23 0.95 us usp000hhsy 2010-07-09T10:13:38.760Z -6.04 105.495 67.6 5.1 mb 68 28.2 1.38 us usp000hf2a 2010-06-26T09:50:43.090Z -8.07 108.089 68 5.9 mwc 381 19.2 1.05 us usp000hejb 2010-06-18T14:42:53.310Z -8.58 111.34 73.9 5 mb 103 44.1 1.31 us usp000he85 2010-06-06T16:57:26.980Z -8.572 112.451 87.6 5.1 mb 134 41 0.96 us usp000hdnp 2010-05-18T11:59:57.840Z -7.845 107.22 34.2 5.3 mwc 243 18.7 1.17 us usp000hcws 2010-01-10T00:25:04.230Z -7.907 107.879 65.2 5.1 mb 155 23.3 1.23 us usp000h5wc 2009-09-02T08:07:44.610Z -7.765 107.449 54 5.3 mb 68 37.7 0.99 us usp000h155 2009-09-02T08:02:21.200Z -7.817 107.393 51.8 5.5 mb 53 30.1 0.97 us usp000h153 2009-09-02T07:55:01.050Z -7.782 107.297 46 7 mwc 312 16.3 0.89 us usp000h152 2009-03-27T20:53:26.670Z -6.963 106.154 47.4 5 mwc 90 46.4 1.08 us usp000gvg4 2008-07-20T06:11:08.200Z -8.636 111.225 77.8 5.3 mwc 201 44.1 1.04 us usp000gcpf 2008-04-27T05:02:39.910Z -7.782 107.876 40.7 5 mwc 162 27.6 0.98 us usp000g4ym 2008-02-03T10:55:33.480Z -8.617 111.521 35 5.1 mb 21 168 1.2 us usp000fy25 2007-12-01T23:01:32.900Z -7.068 105.993 87.7 5 mwc 136 44.9 0.92 us usp000ftps 2007-11-10T02:43:46.200Z -8.565 107.854 34.2 5.3 mb 119 31.4 1.04 us usp000fs9w 2007-09-24T17:35:29.110Z -5.875 105.572 35 5 mb 137 59.5 1.11 us usp000fp1a 2007-09-11T16:36:09.930Z -5.73 105.53 10 5 mb 114 47.5 0.82 us usp000fmx7 2007-08-08T17:05:04.920Z -5.859 107.419 280 7.5 mwc 228 30.7 1.13 us usp000fj6k 2007-08-08T17:04:57.840Z -5.926 107.681 291.2 6.1 mb 531 32.1 1.05 us usp000fj6j 2007-01-31T20:31:34.840Z -7.817 107.193 51.7 5.5 mwb 275 36.2 0.93 us usp000f41s 2006-12-31T13:12:49.520Z -8.257 111.672 139.9 5 mb 52 75.2 0.85 us usp000f1e1 2006-12-23T22:59:40.530Z -6.785 105.635 30 5.5 mwb 206 48.5 0.81 us usp000f0uy 2006-11-28T13:54:07.740Z -7.735 107.798 38.4 5.1 mb 47 69.3 0.83 us usp000ez1h 2006-05-26T22:53:58.920Z -7.961 110.446 12.5 6.3 mwb 213 42.8 1.32 us usp000ej1c 2006-05-19T22:19:23.010Z -8.687 105.977 28.5 5 mb 55 49.2 0.89 us usp000ehh3 2006-03-16T05:43:09.700Z -7.366 106.72 20.3 5 mwc 187 47.1 0.84 us usp000ec8v 2006-02-26T14:42:59.610Z -8.66 106.04 3.6 5.4 mwc 117 46.9 0.93 us usp000eax5 2006-01-20T20:26:51.270Z -6.082 113.064 579.6 5.2 mwc 116 39.1 0.68 us usp000e8j5 2005-09-24T04:40:24.720Z -5.845 105.498 57.9 5 mwc 103 51.3 1.05 us usp000e02j 2005-06-22T20:08:23.870Z -7.847 107.463 38.3 5.3 mwc 130 45.5 1 us usp000dtpj 2005-03-14T15:37:01.030Z -8.618 112.716 120.7 5 mwc 145 20 0.98 us usp000dj73 2005-03-07T13:51:03.570Z -8.524 113.297 116.8 5.1 mb 69 95.6 0.77 us usp000dhrc 2004-09-25T16:04:16.970Z -7.753 107.963 86.6 5.1 mwc 102 28.4 1.1 us usp000d51g 2004-08-19T06:33:31.360Z -8.291 109.794 79.9 5.2 mwc 133 29.7 1.04 us usp000d2en 2004-07-17T01:20:33.390Z -6.519 113.367 53.8 5.1 mb 66 65.7 1.2 us usp000d0ab 2004-06-06T09:38:07.280Z -6.04 113.109 579.4 5.9 mwb 417 26.5 0.87 us usp000cx7y 2004-05-09T22:25:31.130Z -6.944 106.178 78.3 5.2 mwc 208 48.4 0.86 us usp000cukk 2004-01-25T09:45:55.350Z -7.994 112.077 115.5 5.2 mwc 77 47.5 0.83 us usp000cjmu 2003-12-03T12:04:04.210Z -6.959 106.274 33 5.4 mwc 153 56 0.94 us usp000cedr 2003-11-24T23:53:57.320Z -7.56 106.203 33 5.4 mwc 144 45.1 1.09 us usp000cdud 2003-09-16T16:41:15.720Z -7.962 111.687 91.3 5 mb 68 52.1 1 us usp000c7qz 2003-09-08T06:26:32.230Z -8.55 110.189 50.6 5.9 mwb 225 45.4 1 us usp000c726 2003-07-19T21:20:37.010Z -8.682 111.227 56.2 5.9 mwb 397 19.4 0.92 us usp000c2x1 2003-05-28T01:34:02.080Z -8.577 112.499 65.3 5.2 mwc 165 57.3 0.97 us usp000bybv 2003-05-14T07:40:36.100Z -8.06 107.315 79.1 6 mwc 172 38.1 1.17 us usp000bx7m 2003-01-06T21:16:01.190Z -8.609 110.645 116.4 5.1 mb 61 0.88 us usp000bn5y 2003-01-05T16:16:35.540Z -6.607 105.645 33 5.2 mwc 32 0.78 us usp000bn45 2002-10-21T12:15:22.700Z -7.374 107.356 65.5 5.1 mwc 103 1.12 us usp000beyz 2002-08-26T17:36:44.700Z -6.398 113.406 33 5.6 mwc 91 1.07 us usp000bap7 2002-08-26T14:55:13.940Z -6.75 105.705 62.9 5.8 mwc 265 0.97 us usp000bap0 2002-06-21T05:07:11.930Z -6.416 113.404 42.2 5.1 mwc 61 0.86 us usp000b6mg 2002-04-07T23:36:59.620Z -7.642 107.868 43.9 5.2 mwc 77 0.94 us usp000b1yj 2002-02-24T13:14:48.390Z -5.7 110.897 545 5.2 mwc 93 0.86 us usp000az3e 2001-10-14T01:10:45.610Z -8.598 110.633 67.2 5.8 mwc 279 0.99 us usp000aqw5 2001-08-07T03:22:35.890Z -8.61 113.832 90.5 5.6 mwc 181 1.01 us usp000akyv 2001-06-28T03:46:28.190Z -6.99 108.275 36.9 5 mb 67 1.3 us usp000ahm8 2001-06-01T23:09:21.650Z -7.742 107.213 90.3 5.3 mwc 134 1.1 us usp000afwf 2001-05-25T05:06:10.680Z -7.869 110.179 143.1 6.3 mwc 283 1.2 us usp000afgk 2001-03-12T23:35:08.370Z -7.206 106.116 33 5.8 mwc 256 0.91 us usp000aba1 2001-03-03T06:08:44.650Z -6.783 106.341 50.8 5 mb 59 0.86 us usp000aatq 2001-01-19T20:26:41.700Z -8.776 112.584 112.3 5.2 mb 54 1.04 us usp000a808 2001-01-14T19:51:26.360Z -8.469 112.635 74.1 5.2 mb 60 0.97 us usp000a7pd 2001-01-07T12:55:46.890Z -8.703 108.893 33 5.4 mwc 99 0.99 us usp000a794 2001-01-04T03:50:11.440Z -8.403 108.374 33 5.3 mb 46 1.01 us usp000a749 2000-10-25T17:42:22.640Z -7.464 107.823 33 5.3 mb 67 1.01 us usp000a2db 2000-10-25T09:32:23.970Z -6.549 105.63 38 6.8 mwc 415 0.89 us usp000a2ck 2000-09-21T14:33:15.520Z -5.711 110.622 560.9 5.4 mwc 187 0.82 us usp000a0mf 2000-09-05T02:58:00.420Z -8.474 112.137 136.9 5.1 mwc 70 0.99 us usp0009zr8 2000-07-12T01:10:42.630Z -6.675 106.845 33 5.4 mwc 103 1.17 us usp0009w9p 2000-06-05T17:43:24.020Z -7.359 106.664 33 5.1 mwc 106 1.13 us usp0009u2j 2000-04-03T11:13:17.970Z -8.652 110.701 82.8 5 mb 1.1 us usp0009qrs 2000-03-10T16:22:47.180Z -8.701 106.379 33 5.5 mwb 1.17 us usp0009phb 2000-02-12T03:11:02.690Z -8.701 112.445 107.9 5.3 mwc 1.33 us usp0009nc4 2000-01-29T08:13:10.730Z -8.633 111.137 60.7 5.4 mwc 1.09 us usp0009mtg 2000-01-08T03:44:54.320Z -6.746 105.646 61.7 5.6 mwc 1 us usp0009kw1 1999-12-21T14:14:57.610Z -6.845 105.555 56 6.5 mwb 1.03 us usp0009k3r 1999-11-12T12:48:19.900Z -7.709 108.043 80.4 5 mb 1.22 us usp0009heh 1999-07-13T15:52:46.640Z -7.094 107.088 33 5.2 mwc 1.17 us usp0009bj9 1999-06-10T01:06:41.860Z -6.582 106.143 109.3 5.2 mwc 0.92 us usp00099ks 1999-05-26T00:37:53.580Z -6.319 106.461 135.7 5.2 mwc 1.05 us usp00098s7 1999-03-07T01:32:27.000Z -5.812 107.558 325.3 5.2 mwc 0.97 us usp0009402 1999-01-26T15:53:26.950Z -8.239 108.724 96 5.3 mwc 0.91 us usp00091rs 1998-12-25T18:15:41.830Z -7.535 107.126 58.5 5.3 mwc 1 us usp000905p 1998-12-12T09:50:48.840Z -7.609 107.863 85.7 5 mwc 0.98 us usp0008zgz 1998-10-17T11:04:04.790Z -6.382 106.872 147.9 5.2 mwc 1.1 us usp0008wnc 1998-09-28T13:34:30.490Z -8.194 112.413 151.6 6.6 mwb 1.39 us usp0008vnn 1998-09-20T21:21:55.120Z -7.77 106.951 67.4 5.5 mwc 1.2 us usp0008v91 1998-08-17T12:43:39.840Z -7.58 107.199 76.6 5.4 mwc 1.23 us usp0008tdu 1998-04-27T23:51:35.720Z -6.082 113.097 590.7 5.7 mwc 1 us usp0008m2a 1998-02-06T21:55:28.470Z -6.007 105.519 33 5 mb 1.11 us usp0008fk2 1997-07-11T09:55:12.580Z -5.697 110.796 574.4 6 mwc 0.99 us usp00084q6 1997-05-10T21:14:16.130Z -5.909 110.676 600 5 mb 0.97 us usp000822j 1997-03-17T08:05:48.410Z -6.614 105.514 33 6.4 mwc 1.1 us usp0007yym 1997-01-10T18:32:41.520Z -5.885 105.584 33 5 mb 1.07 us usp0007vqe 1996-12-13T20:33:50.650Z -6.2 105.574 55.3 5.1 mwc 1.47 us usp0007udb 1996-12-12T18:23:53.020Z -6.03 112.959 606.2 5.1 mwc 1.17 us usp0007uba 1996-12-09T03:54:16.310Z -7.936 107.489 50.9 6.1 mwc 1.25 us usp0007u5k 1996-11-06T17:04:33.840Z -7.761 106.981 33 5.4 mwc 1.15 us usp0007sj5 1996-05-19T18:18:21.000Z -6.037 112.328 612.5 5.4 mwc 1.1 us usp0007hk9 1995-10-08T05:15:53.350Z -6.89 106.318 67.5 5.2 mw 1 us usp00074st 1995-07-27T12:21:32.430Z -8.683 111.206 71.5 5.4 mw 1.2 us usp000711e 1995-05-28T02:42:27.150Z -6.783 107.277 115.8 5.1 mb 1 us usp0006y6b 1995-05-16T08:40:45.450Z -8.614 106.013 31 5.2 mb 1 us usp0006xf3 1995-05-10T19:01:12.470Z -8.327 112.913 142.6 5 mb 0.9 us usp0006x5s 1995-05-05T17:19:19.390Z -8.725 111.034 76.5 5.9 mw 1.2 us usp0006wxp 1995-02-24T06:24:06.750Z -8.109 105.72 38.8 5.1 mb 0.9 us usp0006tex 1995-01-12T14:37:43.550Z -7.199 106.587 78.2 5.2 mw 1.4 us usp0006r82 1994-10-04T09:54:20.770Z -5.958 110.22 652.2 5 mb 0.6 us usp0006kdd 1994-09-28T17:33:58.210Z -5.731 110.364 628.2 6 mw 0.9 us usp0006k58 1994-09-28T16:39:51.670Z -5.786 110.352 637.5 6.6 mw 1 us usp0006k55 1994-07-17T04:42:41.940Z -5.941 107.515 314.1 5.2 mw 1.2 us usp0006ft2 1994-07-01T01:57:57.160Z -7.96 109.137 101.4 5.1 mw 1.2 us usp0006f13 1994-03-16T00:02:36.190Z -7.149 108.542 33 5 mb 1.2 us usp00069rv 1994-03-13T10:50:11.150Z -7.666 111.429 97.2 5 mb 1.3 us usp00069n6 1994-01-18T06:46:10.120Z -8.364 113.187 140.7 5.2 mb 1.1 us usp00066tf 1993-11-05T07:02:06.110Z -7.033 106.101 74.8 5.2 mw 1.1 us usp0006375 1993-10-13T07:08:36.000Z -8.438 106.817 33 5.1 mb 1.5 us usp00061wn 1993-08-26T01:43:10.470Z -8.672 110.358 66.4 5.3 mb 1.1 us usp0005z77 1993-07-17T23:56:36.570Z -6.091 113.09 567.4 5.3 mw 1 us usp0005wv6 1993-04-10T17:58:07.250Z -8.387 111.254 88.2 5 mb 1.1 us usp0005qzf 1993-03-22T05:24:56.030Z -8.664 110.882 52.4 5.1 mb 1.1 us usp0005pxy 1993-03-15T04:15:35.550Z -7.292 106.796 73.5 5.1 mb 1.4 us usp0005phe 1992-12-27T21:49:04.470Z -6.087 113.05 600.9 6 mw 0.9 us usp0005jyn 1992-12-14T19:57:49.480Z -8.205 107.169 33 5 mb 1.4 us usp0005j9z 1992-11-21T03:04:08.050Z -8.643 110.423 47.9 5.5 mw 1.1 us usp0005h4g 1992-11-18T21:26:13.030Z -7.268 106.193 48.2 5.3 mb 1.2 us usp0005h1e 1992-10-30T16:53:52.590Z -7.878 107.05 64.6 5.2 mw 1 us usp0005g34 1992-10-25T14:37:31.770Z -5.83 106.78 104.1 5.3 mb 0.6 us usp0005fvu 1992-10-18T17:42:51.450Z -6.74 105.602 37.4 5.8 ms 1.1 us usp0005fhu

## Lampiran 3
Script pembuatan grafik energi time series : a=dlmread( 'plotjawaenergi.txt' ); x=a(:,1) y=a(:,4) plot(x,y) title( 'Grafik Time Series Energi Gempa') xlabel( 'Month' ) ylabel( 'Energi' )
Script pembuatan grafik  magnitude time series : a=dlmread( 'plotjawamonth.txt' ); x=a(:,1) y=a(:,2) plot(x,y) title( 'Grafik Magnitude vs Month' ) xlabel( 'Month' ) ylabel( 'Magnitude' )
Script pembuatan grafik  Periodogram : infile=dlmread( 'periodogramiislagi.txt' ); n=length(infile);
month=infile(:,1); % waktu kejadian gempa yang dikonversi eq=infile(:,2); %magnitudo yang telah dikoreksi baseline
Y=fft(eq);
Y(1)=[]
n=length(Y); power=abs(Y(1:floor(n/2))).^2; nyq=1/2; freq=(1:n/2)/(n/2)*nyq; plot(freq,power); xlabel( 'Siklus/bulan' ); title( 'Periodogram' )
hold on ; index=find(power==max(power)); mainFreqStr=num2str(freq(index)); plot(freq(index),power(index),'r','Markersize',25); text(freq(index)-0.1,power(index),['Frekuensi Dominan:',mainFreqStr]);
hold off ;

## Lampiran 4
Script pembuatan grafik power vs periode: infile=dlmread( 'periodevspoweriis.txt' ); n=length(infile);
month=infile(:,1); % waktu kejadian gempa yang dikonversi ke julian month eq=infile(:,2); % perubahan energi (energi – energi rata-rata yang telah dikoreksi baseline)
Y=fft(eq);
Y(1)=[]
n=length(Y); power=abs(Y(1:floor(n/2))).^2; nyquist=1/2; freq=(1:n/2)/(n/2)*nyquist;
period=1./freq; plot(period,power); title( 'Power vs Periode' ) ylabel( 'Power' ) xlabel( 'Periode (siklus/bulan)' );
hold on ; index=find(power==max(power)); mainFreqStr=num2str(freq(index)); plot(freq(index),power(index),'r','Markersize',25); text(freq(index)+5,power(index),['Periode:',mainFreqStr]);
hold off ;