RELOKASI HIPOSENTER GEMPABUMI MENGGUNAKAN

METODE MODIFIED JOINT HYPOCENTER

DETERMINATION (MJHD) UNTUK ANALISIS ZONA
SUBDUKSI SUMATERA BAGIAN SELATAN
(Skripsi)

Oleh
Deswita Sari
1315051013

KEMENTERIAN RISET, TEKNOLOGI DAN PENDIDIKAN TINGGI

UNIVERSITAS LAMPUNG
FAKULTAS TEKNIK
JURUSAN TEKNIK GEOFISIKA
2017
 ABSTRACT

HYPOCENTER RELOCATION OF EARTHQUAKE USING MODIFIED

JOINT HYPOCENTER DETERMINATION (MJHD) METHOD FOR
ANALYSIS SUBDUCTION ZONE AT THE PART OF SOUTH
SUMATERA

By

Deswita Sari

The part of south Sumatera is very vulnerable region in case of earthquake

disaster caused by convergent boundary of two tectonic plates Indo-Australian
Plates and Eurasian Plates. Precise hypocenter analysis is needed to understand
about the accurate tectonic setting such as subduction zone in the area.
Hypocenter relocation is used to recalculate earthquake hypocenter to become
more accurate. To produce a more accurate hyposenter this hyposenter relocation
is done by using the method of Modified Joint Hypocenter Determination
(MJHD). Relocation using the Modified Joint Hypocenter Determination (MJHD)
method uses IASP91 wave velocity which assumes that the inner structures are
heterogeneous. In this study, used data P-wave and S-wave arrival time in the
period January 2010 s.d December 2016 with coordinates -3.5º s.d -9º LS - 99º s.d
106.5º BT. The results of the relocation using MJHD showed a change of
earthquake hypocenter shown by RMS (Root Mean Square) value ranging from
0.2 s.d 0.5. There are three subduction of the part in south sumatra. The
subduction zone formed in Bengkulu is about 26.78º, the subduction zone of
Lampung is around 30.225º and the subduction of the Sunda Strait is about
52.53º. Subduction zone of Bengkulu at depth of 250 km, Lampung and Sunda
Strait at depth 400 km.

Keywords: Hypocenter relocation, Modified Joint Hypocenter Determination

(MJHD), Subduction Zone

i
 ABSTRAK

RELOKASI HIPOSENTER GEMPABUMI MENGGUNAKAN METODE

MODIFIED JOINT HYPOCENTER DETERMINATION (MJHD) UNTUK
ANALISIS ZONA SUBDUKSI SUMATERA BAGIAN SELATAN

Oleh

Deswita Sari

Daerah Sumatera bagian Selatan merupakan daerah yang rawan terhadap bencana
gempabumi karena adanya aktifitas tumbukan lempeng tektonik yaitu Lempeng
Indo-Australia dan Lempeng Eurasia. Untuk memahami kondisi tektonikyang
tepat seperti pola zona subduksi diperlukan analisis hiposenter di daerah tersebut.
Sehingga dilakukanlah relokasi hiposenter untuk menentukan ulang hiposenter
gempabumi menjadi lebih akurat. Untuk menghasilkan hiposenter yang lebih
akurat ini dilakukanlah relokasi hiposenter dengan menggunakan metode
Modified Joint Hypocenter Determination (MJHD). Relokasi dengan metode
Modified Joint Hypocenter Determination (MJHD) menggunakan kecepatan
gelombang IASP91 yang mengasumsikan bahwa struktur dalam bumi bersifat
heterogen. Pada penelitian ini data yang digunakan berupa data arrival time
gelombang P dan S pada rentang waktu Januari 2010 s.d Desember 2016 dengan
koordinat -3.5º s.d -9º LS – 99º s.d 106.5º BT. Hasil dari relokasi menggunakan
MJHD menunjukkan adanya perubahan hiposenter gempabumi yang ditunjukkan
dengan nilai RMS (Root Mean Square) berkisar 0.2 s.d 0.5. Terdapat 3 sudut
penunjaman didaerah Sumatera bagian Selatan. Sudut penunjaman yang terbentuk
di Bengkulu sekitar 26.78º, sudut penunjaman Lampung sekitar 30.225º dan sudut
penunjaman Selat Sunda sekitar 52.53º. Masing-masing kedalaman penunjaman
daerah Bengkulu yaitu sekitar 250 km, Lampung dan Selat Sunda sekitar 400 km.

Kata kunci: Relokasi hiposenter, Modified Joint Hypocenter Determination

(MJHD), Zona Subduksi

ii
RELOKASI HIPOSENTER GEMPABUMI MENGGUNAKAN
METODE MODIFIED JOINT HYPOCENTER
DETERMINATION (MJHD) UNTUK ANALISIS ZONA
SUBDUKSI BAGIAN SELATAN

Oleh

DESWITA SARI

Skripsi

Sebagai salah satu syarat untuk mencapai gelar

SARJANA TEKNIK

pada

Jurusan Teknik Geofisika

Fakultas Teknik Universitas Lampung

KEMENTRIAN RISET TEKNOLOGI DAN PENDIDIKAN TINGGI

JURUSAN TEKNIK GEOFISIKA
FAKULTAS TEKNIK
UNIVERSITAS LAMPUNG
2017
 RIWAYAT HIDUP

Penulis dilahirkan di Kotabumi, Lampung Utara pada

tanggal 25 Desember 1994, sebagai anak pertama dari

dua bersaudara dari pasangan Bapak Margiono dan Ibu

Yuliati. Penulis memulai pendidikan di TK Nurul

Muttaqin pada tahun 2000, kemudian melanjutkan

Sekolah Dasar di SDN 6 Kotabumi dari tahun 2001

hingga 2007. Penulis melanjutkan sekolah menegah pertama di SMP N 3 Cilacap

pada tahun 2007 hingga 2008 kemudian pindah ke SMP N 1 Kotabumi dan lulus

pada tahun 2010. Selanjutnya penulis meneruskan sekolah di SMAN 3 Kotabumi

pada tahun 2010 hingga akhirnya lulus tahun 2013.

Pada tahun 2013 penulis terdaftar sebagai mahasiswa Universitas Lampung

Fakultas Teknik Jurusan Teknik Geofisika melalui jalur SNMPTN. Selama

menjadi mahasiswa, penulis aktif di beberapa organisasi kemahasiswaan seperti

Himpunan Mahasiswa Teknik Geofisika (HIMA TG) pada periode 2014/2015 dan

2015/2016 sebagai anggota Danus (Dana dan Usaha). Selain organisasi internal

kampus, penulis juga aktif dalam organisasi eksternal yaitu sebagai anggota

vii
AAPG SC Unila, anggota SEG (Society of Exploration Geophysicsts) Unila dan

anggota HAGI (Himpunan Ahli Geofisika Indonesia).

Pada bulan Juli hingga Agustus 2016 penulis melaksanakan Kerja Praktik di Pusat

Sumber Daya Geologi (PSDG), Bandung dengan melakukan survei lapangan dan

dengan judul “Akuisisi Data Metode Magnetotellurik Daerah Panasbumi ‘DS’

Sumatera Barat” selama kurang lebih satu bulan. Pada tahun 2017 penulis

melaksanakan Tugas Akhir (TA) selama kurang lebih 2 bulan di Badan

Meteorologi Klimatologi dan Geofisika (BMKG), Kotabumi dengan judul

“Relokasi Hiposenter Gempabumi Menggunakan Metode Modified Joint

Hypocenter Determination (MJHD) Untuk Analisis Zona Subduksi Sumatera

Bagian Selatan”. Hingga akhirnya penulis berhasil menyelesaikan pendidikan

sarjananya pada 12 Juni 2017.

viii
Aku persembahkan karya kecil ini untuk:

Allah SWT

Ayahanda tercinta, Bapak Margiono

Ibunda terkasih, Ibu Yuliati

Saudara kandungku,
Agas Indra Kusuma
dan Keluarga besarku

Teknik Geofisika UNILA 2013

Keluarga Besar Teknik Geofisika UNILA
Almamater Tercinta UNILA

Sahabat-sahabatku tercinta
Serta Calon Imamku kelak.
 ﴾٧٢﴿ ً‫ضيَّة‬
ِ ْ‫اضيَةً َّمر‬
ِ ‫ك َر‬ ْ ‫يَا أَيَّتُهَا النَّ ْفسُ ْال ُم‬
ِ ِّ‫﴾ ارْ ِج ِعي إِلَ ٰى َرب‬٧٢﴿ ُ‫ط َمئِنَّة‬

﴾٠٣﴿ ‫﴾ َوا ْد ُخلِي َجنَّتِي‬٧٢﴿ ‫فَا ْد ُخلِي فِي ِعبَا ِدي‬

Hai jiwa yang tenang, kembalilah kepada Tuhanmu dengan hati yang
puas lagi diridhai; lalu masuklah ke dalam jemaah hamba-hamba-Ku,
dan masuklah ke dalam surga-Ku (QS al-Fajr [89]: 27-30).

Dipersimpangan jalan aku berdiri, aku hilang arah, Hampa

Hanya Tuhanku, Allah yang mampu menyelamatkan

Orang bodoh seringkali beralasan sabar terhadap segala sesuatu yang

sebenarnya dia mengalah dengan kedaan tanpa pernah berusaha

(Albert Einstein)

Tugas Sains adalah untuk menemukan keindahan alam

Ku tengadahkan tangan, Ku basuhkan doa, Ku niatkan karya kecil ini untuk

orang-orangku tersayang
 KATA PENGANTAR

Alhamdulillahirobbilalamin, segala puji dan syukur bagi ALLAH SWT yang telah
memberikan nikmat, karunia dan perlindungan-Nya sehingga penulis dapat
menyelesaikan Skripsi yang berjudul “RELOKASI HIPOSENTER
GEMPABUMI MENGGUNAKAN METODE MODIFIED JOINT
HYPOCENTER DETERMINATION (MJHD) UNTUK ANALISIS ZONA
SUBDUKSI BAGIAN SELATAN” sebagai salah satu bagian dari kurikulum
dan salah satu syarat bagi penulis untuk menyelesaikan studi sebagai Sarjana
Teknik pada Jurusan Teknik Geofisika, Fakultas Teknik, Universitas Lampung.
Skripsi ini merupakan hasil kegiatan Tugas Akhir di BMKG (Badan Meteorologi
Klimatologi dan Geofisika) Kotabumi. Namun demikian, penulis menyadari
masih banyak ketidaksempurnaan dan banyak kelemahan dalam laporan Tugas
Akhir ini. Oleh karena itu penulis mengharapkan kritik dan saran yang
membangun agar dapat memperbaiki dan menyempurnakan nya. Semoga Skripsi
ini dapat bermanfaat bagi kita semua.

Bandar lampung, Juli 2017

Penulis,

Deswita Sari

xi
 SANWACANA

Skripsi yang berjudul: “Relokasi Hiposenter Gempabumi Menggunakan

Metode Modified Joint Hypocenter Determination (MJHD) Untuk Analisis
Zona Subduksi Sumatera Bagian Selatan” adalah salah satu syarat untuk
memperoleh gelar Sarjana Teknik, Jurusan Teknik Geofisika, Fakutas Teknik,
Universitas Lampung.

Penyusunan skripsi ini tidak akan terwujud tanpa adanya dukungan, bantuan dan
bimbingan dari berbagai pihak. Oleh karena itu penulis mengucapkan terimakasih
kepada :

1. Allah SWT yang selalu mengawasi dan meridhoi setiap proses sampai skripsi

ini selesai, sehingga tiada alasan bagi penulis untuk berhenti bersyukur

“Alhamdulillah”.

2. Nabi Muhammad saw yang memberikan teladan kepada seluruh umatnya,

terutama penulis untuk selalu ingin menjadi orang yang lebih baik lagi.

3. Kedua orang tua ku tercinta, Ibu Yuliati dan Bapak Margiono serta adik

Agas Indra Kusuma yang telah memberikan kasih sayang dan kesabaran

serta doa yang terus di berikan tanpa rasa lelah.

4. Bapak Dr. Ahmad Zaenudin, S.Si., M.T selaku Ketua Jurusan Teknik

Geofisika sekaligus pembimbing utama yang telah membantu, memberikan

saran dan bimbingan serta memberikan motivasi selama penelitian hingga

penulisan skripsi.

xiii
5. Bapak Dr. Nandi Haerudin, S.Si., M.Si selaku pembimbing kedua yang

telah memberikan arahan dan bimbingan penulisan skripsi

6. Bapak Rustadi, S.Si., M.T sebagai penguji yang telah memberi kritik, saran

dan bimbingan dalam perbaikan-perbaikan skripsi

7. Bapak Syamsurijal Rasimeng, M.Si sebagai pembimbing akademik yang

telah membantu dan mendukung terselesaikan skripsi ini.

8. Bapak Prof. Drs. Suharno, M.S., M.Sc., Ph.D selaku Dekan Fakultas

Teknik Universitas Lampungmasukan dan nasehat, baik untuk skripsi

ataupun untuk masa depan penulis.

9. Seluruh Staf Tata Usaha Jurusan Teknik Geofisika Unila, Pak Marsuno dan

Mbak Dewi yang telah memberi banyak bantuan dalam proses administrasi;

10. Badan Meteorologi Klimatologi dan Geofisika (BMKG) Kotabumi sebagai

institusi yang telah memberi kesempatan untuk melaksanakan Tugas Akhir.

11. Bapak Joharman, S.H selaku pimpinan BMKG Kotabumi yang telah

mengizinkan dan membantu penulis dalam melaksanakan Tugas Akhir.

12. Mba Fadiah dan Mas Gatut yang telah banyak membantu dan memotivasi

serta menjadi mentor yang baik bagi penulis selama melaksanakan tugas

akhir di BMKG Kotabumi.

13. Mba Juwita, Mba Ferina, Mba Dewi, Mba Vibri, Mas Devid, Pak Tris,

Pak Agung, Pak Rudi serta seluruh pegawai BMKG Kotabumi yang telah

banyak membantu penulis saat melaksanakan Tugas Akhir.

14. Bapak dosen Jurusan Teknik Geofisika atas didikan, bimbingan, serta ilmu

pengetahuan yang telah diberikan.

xiii
15. BMKG pusat yang telah mengizinkan dan memberikan data untuk

menyelesaikan Tugas Akhir.

16. Teman seperjuangan selama melaksanakan tugas akhir di BMKG Kotabumi

yaitu Pipit dan Ulfa yang telah berbagi ilmu dan memotivasi penulis.

17. Teman-teman Teknik Geofisika 2013 yang telah memberikan dukungan dan

keceriaan setiap harinya.

18. Keluarga Besar Teknik Geofisika Unila angkatan 2007, 2008, 2009, 2010,

2011, 2012, 2014, 2015 dan 2016 yang telah memberikan dukungan, do’a dan

semangat untuk penulis.

19. Sahabat perjuangan pada saat KKN di Sidomakmur, PenawarTama, Tulang

Bawang yaitu Ayu Novita Sari, Desy Desmania, Yunita Elsa Pane,

Rafinko Anggriawan dan Varga Desnar Zendya tempat berbagi suka dan

duka.

20. Yang Terakhir namun tak kalah pentingnya, Danu Wahyu Purnawan yang

selalu setia memberikan motivasi tiada hentinya dalam suka maupun duka.

Semoga segala kebaikan dan bantuan yang telah diberikan kepada penulis dibalas

oleh Allah SWT dan bermanfaat bagi kita semua.

Bandar Lampung, Juli 2017

Penulis,

Deswita Sari

xiv
 DAFTAR ISI

Halaman
ABSTRACT ....................................................................................................... i

ABSTRAK ......................................................................................................... ii

HALAMAN JUDUL ......................................................................................... iii

HALAMAN PERSETUJUAN ......................................................................... iv

HALAMAN PENGESAHAN ........................................................................... v

HALAMAN PERNYATAAN ........................................................................... vi

RIWAYAT HIDUP ........................................................................................... vii

HALAMAN PERSEMBAHAN ....................................................................... ix

MOTTO ............................................................................................................. x

KATA PENGANTAR ....................................................................................... xi

SANWACANA .................................................................................................. xii

DAFTAR ISI ...................................................................................................... xv

DAFTAR GAMBAR .........................................................................................xvii

DAFTAR TABEL .......................................................................................... xix

I. PENDAHULUAN
A. Latar Belakang ...................................................................................... 1
B. Tujuan Penelitian .................................................................................. 4
C. Batasan Masalah ................................................................................... 4
D. Manfaat Penelitian ................................................................................ 5

xv
II. TINJAUAN PUSTAKA
A. Daerah Penelitian .................................................................................. 6
B. Letak Stasiun Seismik BMKG Kotabumi ............................................. 7
C. Fisiografi ............................................................................................... 8
D. Cekungan Sumatera Bagian Selatan ..................................................... 9
E. Sistem Sesar Sumatera dan Sejarah Kegempaan Lampung .................. 10

III. TEORI DASAR

A. Klasifikasi Gempabumi ........................................................................ 12
B. Teori Elastic Rebound........................................................................... 15
C. Pergerakan Lempeng Tektonik ............................................................. 16
D. Zona Subduksi ...................................................................................... 19
E. Tektonik Sumatera ................................................................................ 22
F. Gelombang Seismik ............................................................................. 27
G. Parameter Sumber Gempa ................................................................... 31
H. Modified Joint Hypocenter Determination (MJHD)............................. 34
I. Menentukan Kelengkungan dan Sudut Penunjaman ........................... 43

IV. METODE PENELITIAN

A. Lokasi dan Tempat Penelitian ............................................................... 45
B. Data Penelitian dan Perangkat Lunak ................................................... 45
C. Pengolahan Data ................................................................................... 46
D. Diagram Alir ......................................................................................... 49
E. Time Schedule ....................................................................................... 50

V. HASIL PENGAMATAN DAN PEMBAHASAN

A. HasilPengamatan................................................................................... 51
1. Sebelum Relokasi ............................................................................ 51
2. Sesudah Relokasi ............................................................................. 52
B. Pembahasan........................................................................................... 59
1. Distribusi Event Gempa dan Pola Penunjaman ............................... 59
2. Sudut Penunjaman ........................................................................... 69

VI. KESIMPULAN DAN SARAN

A. Kesimpulan ........................................................................................... 73
B. Saran ..................................................................................................... 74

DAFTAR PUSTAKA

LAMPIRAN

xvi
 DAFTAR GAMBAR

Gambar Halaman
1. Daerah penelitian ........................................................................................... 6
2. Letak stasiun seismik BMKG ........................................................................ 7
3. Peta fisiografi daerah penelitian ...................................................................... 8
4. Model Elastic Rebound .................................................................................. 15
5. Pergerakan lempeng tektonik ......................................................................... 18
6. Batas-batas lempeng tektonik ........................................................................ 18
7. Zona subduksi ................................................................................................ 19
8. Bentuk zona subduksi di Indonesia ................................................................ 20
9. Peta tektonik Indonesia .................................................................................. 23
10. Sumatera Fault Zone ...................................................................................... 25
11. Zona subduksi Megathrust Sumatera ............................................................. 26
12. Simulasi gerakan gelombang P dan S ............................................................ 29
13. Gelombang love dan gelombang rayleigh ..................................................... 31
14. Diagram Wadati ............................................................................................. 32
15. Jarak hiposenter .............................................................................................. 33
16. Ilustrasi relokasi menggunakan MJHD ........................................................... 41
17. Model kecepatan IASP91................................................................................ 42
18. Ilustrasi garis singgung kurva ........................................................................ 44
19. Data dari katalog BMKG ............................................................................... 46
20. Data dalam format MJHD ............................................................................. 46
21. Tampilan input program station ..................................................................... 47
22. Tampilan parameter iterasi yang digunakan .................................................. 47
23. Contoh plotting gempabumi dengan GMT .................................................... 48

xvii
24. Diagram alir ................................................................................................... 49
25. Event gempabumi sebelum relokasi ............................................................... 51
26. Distribusi stasiun ............................................................................................ 53
27. Event gempabumi sesudah relokasi ............................................................... 54
28. Sebaran hiposenter dalam 3D sebelum relokasi ............................................ 55
29. Sebaran hiposenter dalam 3D sesudah relokasi ............................................. 56
30. Irisan penampang vertikal .............................................................................. 59
31. Pola penunjaman daerah Bengkulu sebelum relokasi ..................................... 60
32. Pola penunjaman daerah Bengkulu sesudah relokasi ..................................... 61
33. Pola penunjaman daerah Lampung sebelum relokasi .................................... 63
34. Pola penunjaman daerah Lampung sesudah relokasi ..................................... 64
35. Pola penunjaman daerah Selat Sunda sebelum relokasi ................................ 67
36. Pola penunjaman daerah Selat Sunda sesudah relokasi ................................. 68
37. Pola kelengkungan penunjaman Bengkulu .................................................... 70
38. Pola kelengkungan penunjaman Lampung .................................................... 70
39. Pola kelengkungan penunjaman Selat Sunda ................................................. 71
40. Model slab 3D (USGS)5 ................................................................................ 72

xviii
 DAFTAR TABEL

Tabel Halaman
1. Data dan parameter sumber gempa Sumatera Fault ........................................ 24
2. Data dan parameter sumber gempa subduksi (Megathrust) ............................ 27
3. Time schedule penyusunan skripsi ................................................................... 50
4. Nilai RMS sebelum relokasi ........................................................................... 57
5. Nilai RMS sesudah relokasi ............................................................................ 58

xix
 BAB I. PENDAHULUAN

A. Latar Belakang

Indonesia berada pada zona tektonik yang sangat aktif karena terletak di daerah

dengan tingkat aktifitas gempabumi tinggi, hal tersebut akibat dari bertemunya

tiga lempeng besar dunia yaitu Lempeng Indo-Australia, Eurasia dan Lempeng

Pasifik. Lempeng Indo-Australia bertabrakan dengan Lempeng Eurasia di lepas

pantai Sumatera, Jawa dan Nusa Tenggara. Sedangkan Lempeng Pasifik di utara

Irian dan Maluku Utara.

Keberadaan interaksi antar lempeng-lempeng ini menempatkan wilayah Indonesia

sebagai wilayah yang sangat rawan terhadap gempabumi. Sumatera merupakan

salah satu pulau yang aktif terjadi gempa, karena adanya penunjaman

(subduction) Lempeng Indo-Australia yang masuk ke dalam Lempeng Eurasia,

sehingga membentuk jalur gempa di Laut Sumatera Indonesia. Sebelah Barat

Pulau Sumatera dikenal dengan jalur gempa mediteran dan di darat Pulau

Sumatera menimbulkan Patahan Besar Sumatera (Great Sumatera Fault) yang

membujur sepanjang Bukit Barisan dan membentang dari Aceh, Sumatera Utara,

Sumatera Barat sampai ke Teluk Lampung, sehingga wilayah Sumatera bagian

 2

Selatan memiliki tingkat seismisitas tinggi terancam oleh bencana alam

gempabumi.

Gempabumi yang terjadi di Pulau Sumatera merupakan implikasi geodinamika

dari deformasi aktif di sekitar Sunda dan Java Trench. Kejadian gempa yang

terjadi di dekat batas pertemuan antara Lempeng Samudera yang menunjam

masuk ke bawah Lempeng Benua diklasifikasikan sebagai Zona Subduksi. Pada

penelitian ini penulis menggunakan data gempa yang terjadi di wilayah Bengkulu,

Lampung dan sekitaran Selat Sunda.

Menurut sudut pandang ilmu kebumian, Sumatera bagian Selatan mempunyai

aktifitas tumbukan dua lempeng. Akibat proses tumbukan dua lempeng yang

berlangsung hingga saat ini menyebabkan adanya zona sesar lokal. Selain itu

wilayah Sumatera bagian Selatan memiliki tingkat ancaman gempabumi yang

cukup tinggi. Hal ini dikarenakan adanya Zona Subduksi yang merupakan batas

antar Lempeng India-Australia yang menunjam kedalam Lempeng Eurasia. Zona

ini berpotensi menimbulkan gempabumi dengan magnitudo yang relatif lebih

besar (Madlazim, 2013).

Salah satu parameter gempabumi yang dapat dihitung yaitu hiposenter

gempabumi. Penentuan hiposenter gempabumi sangatlah penting didalam dunia

seismologi. Hal ini sangat diperlukan dalam analisis struktur tektonik secara

detail, misalnya untuk identifikasi zona patahan maupun pola Zona Subduksi.

Namun, parameter hiposenter yang dihasilkan masih kurang optimal karena hanya
 3

untuk memberikan informasi sesegera mungkin kepada masyarakat tentang

bahaya gempabumi. Sehingga, perlu dilakukan studi lebih lanjut untuk merelokasi

parameter hiposenter gempabumi yang telah dihasilkan sebelumnya.

Penentuan parameter hiposenter sebelumnya ditentukan dengan menggunakan

waktu tiba gelombang P dan S dari beberapa stasiun menggunakan pendekatan

Single Event Determination (SED) yang dapat menghasilkan origin time dari

setiap event gempabumi. Dari pendekatan SED ini hanya menghasilkan hiposenter

yang belum akurat karna hanya menggunakan kecepatan bumi 1D yang belum

termodelkan. Sehingga untuk menyelesaikan permasalahan ini dilakukanlah

relokasi hiposenter gempabumi dengan menggunakan metode Modified Joint

Hypocenter Determination (MJHD) yang dikembangkan oleh Hurukawa dan

Imoto, (1992). MJHD ini menggunakan model kecepatan global IASP91. Model

kecepatan IASP91 merupakan suatu model kecepatan global yang dihasilkan dari

International of Seismology and Physics of the Earth’s Interior (IASPEI) dengan

memperhatikan keheterogenan suatu lapisan bumi. Model kecepatan IASP91 ini

merupakan hasil pemodelan dengan menggunakan ribuan data gelombang P dan

gelombang S yang direkam oleh ribuan seismometer diseluruh dunia.

Relokasi hiposenter merupakan koreksi dari lintang, bujur, dan kedalaman dari

gempabumi. Dalam studi ini, penulis melakukan relokasi hiposenter gempabumi

di Sumatera bagian selatan dengan menggunakan metode Modified Joint

Hypocenter Determination (MJHD). Metode MJHD dipilih karena metode ini

mampu merelokasi posisi gempabumi relatif lebih akurat dengan struktur dalam
 4

bumi yang heterogen dan distribusi stasiun yang tidak merata. Metode ini dapat

menghitung banyak data gempabumi secara simultan (Maung, 2009).

Kelebihan dalam relokasi menggunakan metode MJHD ini adalah adanya

penambahan koreksi stasiun pada proses penentuan ulang hiposenter gempabumi.

Selain itu juga dengan menggunakan kecepatan IASP91 akan menghasilkan

hiposenter yang lebih akurat karena asumsi dari kecepatan yang digunakan berupa

kecepatan bumi global dengan asumsi struktur bumi yang heterogen.

B. Tujuan Penelitian

Secara garis besar tujuan dari dilakukannya penelitian ini antara lain:

1. Melakukan relokasi data hiposenter gempabumi wilayah Sumatera bagian

Selatan untuk mengetahui perbedaan hasil sebaran kegempaan sebelum

dilakukan relokasi dan sesudah dilakukan relokasi

2. Membuat pemodelan hiposenter gempabumi secara 3D

3. Menganalisis pola dan sudut Zona Subduksi di wilayah Sumatera bagian

Selatan berdasarkan distribusi hiposenter hasil dari relokasi

C. Batasan Masalah

Batasan masalah pada penelitian ini yaitu:

1. Melakukan relokasi gempabumi Sumatera bagian Selatan periode waktu

Januari 2010 s.d Desember 2016 dengan data yang digunakan adalah data

arrival time gelombang P dan gelombang S

 53

2. Metode yang digunakan dalam penelitian ini adalah metode relokasi

hiposenter Modified Joint Hypocenter Determination (MJHD)

D. Manfaat Penelitian

Manfaat yang dilakukan dalam penelitian ini adalah:

1. Merelokasi hiposenter gempabumi untuk wilayah Sumatera bagian Selatan

pada rentang waktu Januari 2010 s.d Desember 2016

2. Mengetahui pola Zona Subduksi wilayah Sumatera bagian Selatan

berdasarkan penyebaran hiposenter hasil relokasi

 BAB II. TINJAUAN PUSTAKA

A. Daerah Penelitian

Daerah penelitian terletak dengan koordinat -3.5º s.d -9º LS dan 99º - 106.5º BT

yang mencakup wilayah Sumatera bagian Selatan yaitu Bengkulu, Lampung dan

Selat Sunda. Dengan kelompok geologi yang ada di daerah tersebut berupa

kelompok sedimen tersier, sedimen kuarter, plutonik tersier, sedimen paleozoikum

dan malihan paleozoikum. Berikut ini adalah peta dan posisi daerah penelitian yang

dapat dilihat pada Gambar 1.

Gambar 1. Daerah penelitian (Sukamto, dkk., 1996)

 7

B. Letak Stasiun Seismik BMKG Kotabumi

Gambar 2. Letak stasiun seismik BMKG (BMKG Kotabumi, 2016)

Berdasarkan letak stasiun yang terdapat pada Gambar 2, BMKG Kotabumi

bertanggung jawab atas 8 stasiun seismik yang ditempatkan dibeberapa daerah

yaitu:

1. Stasiun Seismik KLI yang berada di Kotabumi Lampung Utara

2. Stasiun Seismik LWLI yang berada di Liwa Lampung Barat

3. Stasiun Seismik BLSI yang berada di Bandar Lampung

4. Stasiun Seismik KASI yang berada di Kotaagung Tanggamus

5. Stasiun Seismik MDSI yang berada di Muara Dua Sumatera Selatan

6. Stasiun Seismik LHSI yang berada di Lahat Sumatera Selatan

7. Stasiun Seismik KLSI yang berada di Sungkai Lampung Utara

8. Stasiun Seismik PMBI yang berada di Palembang Sumatera Selatan

 8

C. Fisiografi

Secara umum daerah penelitian ini dapat dibagi menjadi tiga satuan morfologi,

yaitu dataran rendah di bagian timur dan timurlaut, pegunungan bergelombang di

bagian tengah dan barat daya, dan daerah pantai berbukit sampai dataran. Daerah

dataran rendah menempati lebih dari 60% luas lembar dan terdiri dari endapan

vulkanoklastika tersier-kuarter dan aluvium dengan ketinggian puluhan meter di

atas muka laut. Pegunungan Bukit Barisan menempati 25-30% luas lembar, terdiri

dari batuan beku dan malihan serta batuan gunungapi muda. Lereng-lereng

umumnya curam dengan ketinggian 500 s.d 1.680 m di atas muka laut (Mangga,

dkk., 1993). Daerah pantai bertopografi beraneka ragam dan seringkali terdiri dari

pebukitan kasar, mencapai ketinggian 500 m di atas muka laut dan terdiri dari

batuan gunungapi tersier dan kuarter serta batuan terobosan yang terdapat pada

Gambar 3.

Gambar 3. Peta fisiografi daerah penelitian (Mangga, dkk., 1993)

 9

D. Cekungan Sumatera Bagian Selatan

Cekungan Sumatera bagian Selatan adalah suatu hasil kegiatan tektonik yang

berkaitan dengan penunjaman Lempeng Indo-Australia, yang bergerak ke arah

utara hingga timur laut terhadap Lempeng Eurasia yang relatif diam. Zona

penunjaman lempeng meliputi daerah sebelah barat Pulau Sumatera dan Selatan

Pulau Jawa. Beberapa lempeng kecil yang berada diantara zona interaksi tersebut

turut bergerak dan menghasilkan zona konvergensi dalam bentuk dan arah.

Penunjaman Lempeng Indo-Australia tersebut dapat mempengaruhi keadaan

batuan, morfologi, tektonik dan struktur di Sumatera Selatan.

Tumbukan tektonik lempeng di Pulau Sumatera menghasilkan jalur busur depan,

magmatik, dan busur belakang. Cekungan Sumatera Selatan telah mengalami tiga

kali proses orogenesis, yaitu :

• Mesozoikum Tengah

• Kapur Akhir sampai Tersier Awal

• Pilo-Plistosen

Secara fisiografis, Cekungan Sumatera Selatan merupakan cekungan Tersier

berarah barat laut-tenggara, yang dibatasi Sesar Semangko dan Bukit Barisan di

sebelah barat daya, Paparan Sunda di sebelah timur laut, Tinggian Lampung di

sebelah tenggara yang memisahkan cekungan tersebut dengan Cekungan Sunda,

serta Pegunungan Dua Belas dan Pegunungan Tiga Puluh di sebelah barat laut yang

memisahkan Cekungan Sumatera Selatan dengan Cekungan Sumatera Tengah.

Posisi Cekungan Sumatera Selatan sebagai Cekungan Busur Belakang, daerah

cekungan ini merupakan Cekungan Busur Belakang berumur tersier yang terbentuk
 10

sebagai akibat adanya interaksi antara Paparan Sunda (sebagai bagian dari

Lempeng Kontinen Asia) dan Lempeng Samudera Hindia. Daerah cekungan ini

meliputi daerah seluas 330 x 510 km2. Kenampakan struktur yang dominan adalah

struktur yang berarah barat laut – tenggara sebagai hasil orogenesa plio – plistosen.

Dengan demikian pola struktur yang terjadi dapat dibedakan atas pola tua yang

berarah utara-selatan dan barat laut – tenggara serta pola muda yang berarah barat

laut – tenggara yang sejajar dengan Pulau Sumatera (Sieh, 2000).

E. Sistem Sesar Sumatera dan Sejarah Kegempaan Lampung

Di Pulau Sumatera, pergerakan Lempeng Indo dan Australia yang mengakibatkan

kedua lempeng tersebut bertabrakan dan menghasilkan penunjaman menghasilkan

rangkaian busur pulau depan (forearch islands) yang non-vulkanik (seperti: P.

Simeulue, P. Banyak, P. Nias, P. Batu, P. Siberut hingga P. Enggano), rangkaian

pegunungan Bukit Barisan dengan jalur vulkanik di tengahnya, serta sesar aktif

’The Great Sumatera Fault’ yang membelah Pulau Sumatera mulai dari Teluk

Semangko hingga Banda Aceh. Sesar Besar ini menerus sampai ke Laut Andaman

hingga Burma. Patahan aktif ini (Sesar Besar Sumatra) diperkirakan bergeser

sekitar sebelas sentimeter per tahun dan merupakan daerah rawan gempa bumi dan

tanah longsor.

Penunjaman Lempeng Indo – Australia juga mempengaruhi geomorfologi Pulau

Sumatera menjadikan bagian barat Pulau Sumatera terangkat, sedangkan bagian

timur relatif turun. Hal ini menyebabkan bagian barat mempunyai dataran pantai

yang sempit dan kadang-kadang terjal. Pada umumnya, terumbu karang lebih
 11

berkembang dibandingkan berbagai jenis bakau. Bagian timur yang turun akan

menerima tanah hasil erosi dari bagian barat (yang bergerak naik), sehingga bagian

timur memiliki pantai yang datar lagi luas. Di bagian timur, gambut dan bakau lebih

berkembang dibandingkan terumbu karang.

Keadaan Pulau Sumatera menunjukkan bahwa kemiringan penunjaman,

punggungan busur muka dan cekungan busur muka telah terfragmentasi akibat

proses yang terjadi. Kenyataan menunjukkan bahwa adanya transtensi (trans-

tension) Paleosoikum tektonik Sumatera menjadikan tatanan tektonik Sumatera

menunjukkan adanya tiga bagian pola (Sieh, 2000). Bagian selatan terdiri dari

lempeng mikro Sumatera, yang terbentuk sejak 2 juta tahun lalu dengan bentuk,

geometri dan struktur sederhana, bagian tengah cenderung tidak beraturan dan

bagian utara yang tidak selaras dengan pola penunjaman.

Lampung Sendiri merupakan salah satu zona rawan bencana gempabumi. Salah

satu gempa merusak di wilayah Lampung yaitu gempa di Liwa Lampung Barat

yang berkekuatan 6,6 SR pada 15 Februari 1994. Akibat gempa tersebut terjadi

kerusakan parah, sekitar 196 jiwa dari beberapa desa dan kecamatan di Lampung

meninggal dunia, sementara jumlah korban luka mencapai 2000 orang. Rata-rata

mereka tewas dan terluka akibat tertimpa reruntuhan bangunan dan jumlah

penduduk yang kehilangan tempat tinggal sekitar 75.000 orang.

 12

BAB III. TEORI DASAR

A. Klasifikasi Gempabumi

Gempabumi merupakan goncangan pada permukaan bumi yang dihasilkan dari

gelombang seismik akibat pelepasan energi secara tiba-tiba dari dalam bumi.

Berdasarkan sumber terjadinya gempa, Hoernes (Subardjo dan Ibrahim, 2004)

mengelompokan menjadi:

1. Gempabumi Tektonik, yaitu gempabumi yang disebabkan oleh aktivitas

pergerakan lempeng tektonik misalnya adanya tumbukan antar lempeng

pembentuk kulit bumi.

2. Gempabumi Vulkanik, yaitu gempabumi yang terjadi karena adanya aktivitas

vulkanik misalnya desakan magma dari gunungapi ke permukaan. Gempabumi

ini memiliki kekuatan kurang dari 4 SR dan termasuk gempabumi sedang.

3. Gempabumi Runtuhan, yaitu gempabumi yang terjadi karena adanya runtuhan

atau longsoran dari massa batuan. Gempabumi ini memiliki kekuatan yang

sangat kecil sehingga getarannya tidak bisa terasa dan hanya bisa terdeteksi oleh

seismograf. Gejala ini disebut dengan tremor dan banyak terjadi di pegunungan.

4. Gempabumi Buatan, yaitu gempabumi yang sengaja dibuat oleh manusia, seperti

ledakan dinamit atau ledakan nuklir untuk mencari bahan tambang.

 13

Berdasarkan dalamnya sumber gempa, Bath (Subardjo dan Ibrahim, 2004)

mengelompokan gempabumi menjadi :

1. Gempabumi dangkal (kedalaman 0 - 60 km)

Gempabumi dangkal yang terjadi di daratan biasanya berasosiasi dengan

patahan-patahan besar yang bergeser atau bergerak akibat pergerakkan lempeng.

Gempabumi dangkal memiliki daya rusak kontsruksi yang sangat kuat. Hal ini

disebabkan jarak hiposenter relatif dekat dengan permukaan sehingga getaran

gempa terasa sangat kuat dipermukaan. Walaupun skala gempanya hanya

memiliki amplitudo 4 - 5 Skala Richter namum mampu merusakkan bangunan.

Sedangkan gempabumi dangkal yang terjadi di laut bisa mengakibatkan tsunami

bila kekuatan gempa ≥ 6 SR.

2. Gempabumi menengah (kedalaman 60 – 300 km)

Bumi terdiri dari tiga lapisan utama yaitu crust, mantle dan core. Crust atau

lithosfer adalah lapisan yang paling terluar dari bumi berbentuk padat dengan

ketebalan lapisan mencapai 100 km. Lithosfer terdiri dari kerak bumi dan bagian

atas selubung. Gempa-gempa menengah terjadi pada kedalaman dibawah kerak

bumi. Sehingga digolongkan sebagai gempagempa yang mungkin tidak

berasosiasi dengan penampakan patahan di permukaan. Namun gempa-gempa

ini masih dapat di perkirakan mekanisme terjadinya.

3. Gempabumi dalam ( > 300 km)

Gempa dalam ini sebenarnya relatif sering terjadi. Namun karena berada pada

kedalaman dibawah 300 km maka manusia tidak bisa merasakan getarannya.

Gempa dalam merupakan gempa-gempa yang disebabkan oleh pergerakan kerak

benua.
 14

Berdasarkan kekuatan gempanya, Hagiwara (Subardjo dan Ibrahim, 2004)

mengklasifikasikan gempabumi menjadi:

1. Gempa sangat besar, M> 8,0

2. Gempa besar, 7,0 < M < 8,0

3. Gempa sedang, 4,5 < M < 7,0

4. Gempa mikro, 1,0 < M < 4,5

5. Gempa ultra mikro, M < 1,0

Pola mekanisme terjadinya gempa diatas tergantung pada keadaan struktur kulit

bumi dan distribusi gaya/stress yang bekerja. Stress yang bekerja pada gempa

tektonik yang terjadi umumnya adalah seragam/uniform, sehingga perbedaan

yang terjadi adalah merupakan perbedaan keadaan struktur/medium daerah

bersangkutan.

Menurut Mogi (1967) pola umum terjadinya gempabumi dibedakan dalam 3

(tiga) jenis :

1. Tipe I : yaitu gempabumi utama (main shock) tanpa didahului gempa

permulaan (fore shock), tetapi diikuti dengan banyak gempabumi susulan

(after shock). Gempabumi tipe ini biasanya terjadi di daerah yang mempunyai

medium homogen dengan stress yang bekerja hampir merata (uniform).

Sebagian besar gempabumi tektonik yang terjadi di bumi tergolong jenis ini.

2. Tipe II : yaitu gempabumi utama (main shock) didahului gempa-gempa

pendahuluan (fore shock) dan kemudian diikuti gempa susulan (after shock)

yang cukup banyak jumlahnya. Gempabumi tipe ini terjadi pada daerah
 15

dengan struktur batuan yang tidak seragam dengan distribusi stress yang

bekerja tidak seragam.

3. Tipe III : yaitu gempabumi dimana tidak terdapat gempa utama (main

shock), biasa disebut gempabumi “swarm”. Gempabumi tipe ini terjadi

dalam daerah yang terbatas, biasanya terjadi didaerah gunung api.

Gempabumi ini terjadi pada daerah yang struktur mediumnya tidak seragam

dengan stres yang bekerja terkonsentrasi pada area yang terbatas.

B. Teori Elastic Rebound

Elastic Rebound Theory yang dikemukakan oleh seorang seismologist Amerika

bernama Reid (Bolt, 1976) menjelaskan proses terjadinya gempabumi yang

terdapat pada Gambar 4.

Gambar 4. Model Elastic Rebound (Subardjo dan Ibrahim, 2004)

Pada keadaan I menunjukkan suatu lapisan yang belum terjadi perubahan bentuk

geologi. Karena di dalam bumi terjadi gerakan yang terus-menerus, maka akan

terdapat stress yang lama kelamaan akan terakumulasi dan mampu merubah bentuk

geologi dari lapisan batuan. Keadaan II yakni suatu lapisan batuan telah

mengandung stress dimana telah terjadi perubahan bentuk geologi. Untuk daerah A

mendapat stress ke atas, dan daerah B mendapat stress ke bawah. Proses ini berjalan

terus hingga stress yang dikandung di daerah ini cukup besar untuk merubahnya

menjadi gesekan antara daerah A dan daerah B. Suatu ketika karena lapisan batuan
 16

sudah tidak mampu lagi untuk menahan stress, maka akan terjadi suatu perpindahan

massa batuan secara tiba-tiba berupa patahan yang melepaskan gelombang seismik.

Pada keadaan III menunjukan lapisan batuan yang sudah patah. Gerakan perlahan-

lahan sesar ini akan berjalan terus, sehingga seluruh proses diatas akan diulangi lagi

dan sebuah gempa akan terjadi lagi setelah beberapa waktu lamanya.

C. Pergerakan Lempeng Tektonik

Pergerakan lempeng tektonik terbagi atas 3 zona (Awaludin, 2011) yang terdapat

pada Gambar 5, yaitu :

1. Zona Divergen

Zona divergen adalah pergerakan dua buah lempeng tektonik atau lebih yang

bergerak saling menjauh satu sama lainnya yang mengakibatkan material mantel

naik keatas atau terjadi pergerakan mantel (mantle convection) membentuk

lantai samudra (sea floor spreading). Pada zona ini juga terdapat pegunungan

bawah laut (mid oceanic ridge). Pergerakan mantel ini terjadi karena adanya

pendinginan dari atas dan pemanasan dari bawah sehingga mantel akan bergerak

keatas. Aktivitas semacam ini menimbulkan gempa tektonik dangkal dan gempa

vulkanik.

2. Zona Konvergen

Zona konvergen merupakan pergerakan dua lempeng tektonik yang bergerak

relatif saling mendekati. Zona konvergen juga ditandai dengan adanya

penghancuran meteri-materi lempeng, sehingga zona ini disebut zona destruktif.

Zona konvergen terbagi dua, yaitu :

 17

a. Zona Tumbukan

Zona tumbukan merupakan pertemua dua lempeng dengan berat jenis sama

yang bergerak relatif saling mendekati. Tumbukan ini menghasilkan

pegunungan lipatan seperti Pegunungan Himalaya dan Pegunungan Andes.

Aktifitas lempeng seperti ini menimbulkan gempa tektonik dangkal dan

gempa vulkanik.

b. Zona Subduksi

Zona Subduksi merupakan pertemuan dua lempeng tektonik yang

mempunyai berat jenis berbeda dan bergerak relatif saling mendekati

sehingga lemepeng yang lebih berat menyusup atau menujam ke bawah

lempeng yang lebih ringan. Zona ini ditandai dengan adanya palung laut atau

trench sebagai batas pertemuan kedua lempeng. Selain itu, pada zona

subduksi juga terdapat rangkaian gunung api yang sejajar trench sebagai

akibat dari melelehnya lempeng yang menujam pada kedalaman 100-400

km. aktifitas ini mengakibatkan terjadinya gempa tektonik dangkal,

menengah dan dalam serta gempa vulkanik.

3. Zona Transform

Zona transform merupakan daerah singgungan dua lempeng yang bergerak

relatif sejajar dan berlawanan arah sehinga pada batas kedua lempeng ini terjadi

gesekan. Aktivitas ini sering menimbulkan gempa dangkal dan bersifat merusak.
 18

Gambar 5. Pergerakan lempeng tektonik (Subardjo dan Ibrahim, 2004)

Batas Divergen

Batas Konvergen

Batas Transform

Gambar 6. Batas-batas lempeng tektonik (Subardjo dan Ibrahim, 2004)

 19

D. Zona Subduksi

Zona Subduksi merupakan daerah zona seismik aktif, sebagian besar gempabumi

terjadi di daerah ini, baik dangkal, menengah maupun dalam sebagai hasil dari

penunjaman lithosfer. Penyebaran titik–titik pusat gempa (hypocenter)

menunjukkan struktur dari penunjaman lithosfer dari mulai batas palung sampai ke

dalam mantel bumi. Terjadinya zona subduksi dimulai dengan dua lempeng (benua

dan samudera) yang bergerak mendekat satu dengan yang lainnya. Kemudian

terjadi konvergen, plat samudera melengkung dan terdorong ke bawah plat benua

yang lebih tebal dan lebih stabil. Proses penunjaman ini memberikan bentuk dari

permukaan bumi (Malik, 2009).

Selain dari terbentuknya palung – palung yang dalam, yang merupakan batas

penunjaman lithosfer. Penunjaman lithosfer juga menyebabkan terjadinya

deformasi dan kompressi di daerah lempeng benua yang dekat dengan daerah

tumbukan, mengakibatkan terjadinya intrusi dari magma sehingga terbentuknya

busur vulkanik. Gambar zona subduksi dapat dilihat pada Gambar 7

Gambar 7. Zona subduksi (Subardjo dan Ibrahim, 2004)

 20

Secara umum di Indonedia terdapat 4 bentuk zona subduksi (Subardjo dan Ibrahim,

2004) yang terdapat pada Gambar 8.

Gambar 8. Bentuk zona penunjaman di Indonesia (Subardjo dan Ibrahim, 2004)

1. Zona Penunjaman pendek

Di sepanjang Sumatera sampai Jawa Barat, kedalaman penunjaman sejauh 180

km, yang berarti disini zona subduksi menunjam sejauh 180 km. Arahnya dari

barat daya ke timur dengan sudut penunjaman 25º, jenis subduksi miring

(oblique fault). Kedalaman palung laut Sumatera sekitar 4500 meter dan palung

laut Jawa mencapai 7000 meter. Lempeng menunjam rata-rata dengan kecepatan

6,8 cm pertahun.

2. Zona penunjaman diskontinu.

Bentuk ini ditemui mulai dari Jawa Tengah sampai Flores. Kedalaman

maksimum 650 km Tetapi pada kedalaman antara 260 – 542 km di Jawa Tengah

dan kedalaman antara 280 – 360 km di Flores Barat terdapat diskontinu

lempeng. Kecepatan penunjaman lempeng sekitar 7,5 cm pertahun. Gunung api

 21

terdapat pada lokasi yang berkaitan dengan kedalaman gempabumi antara

100 dan 200 km.

3. Zona penunjaman berbentuk permukaan cekung.

Bentuk zona ini terdapat mulai dari Alor sampai kepulauan Kai (daerah Laut

Banda). Di wilayah Laut Banda bentuk penunjaman lempeng lebih komplek.

Sebelah selatan terdapat Palung Timor dan di utara adalah Palung Seram. Kedua

palung ini melingkar membentuk setengah lingkaran mulai dari selatan pulau

Timor, Tanimbar, berbelok ke atas di sebelah timur Kepulauan Kai dan

kemudian berbalik ke arah barat di sebelah utara Pulau Seram dan Buru.

Kedalaman Palung Timor sekitar 2500 meter, Palung Seram antara 4000 – 5000

meter, dan Basin Weber mencapai kedalaman 7000 meter. Zona subduksi di

daerah Laut Banda adalah berbentuk suatu permukaan cekung, lempeng-

lempeng tektonik menunjam dari arah utara dan dari arah selatan yang bertemu

di Laut Banda. Kedalamannya berkurang dari arah barat ke timur, di sebelah

barat dekat Pulau Alor penunjaman zona subduksi 650 km dan di sebelah timur

dekat Pulau Tanimbar penunjamannya 96 km. Sudut penunjaman juga berkurang

dari arah barat ke timur, di palung sebelah Selatan (Timor) dari 74º sampai

dengan 16º dan di palung sebelah utara dari 57º sampai 14º terdapat diskontinu

pada lempeng di sebelah selatan.

4. Zona penunjaman berbentuk permukaan cembung.

Ditemui didaerah Maluku, Pada daerah Maluku zona subduksi lebih rumit

bentuknya. Di daerah ini terdapat beberapa palung, yaitu Palung Maluku yang

bersambung dengan Palung Filipina diutaranya. Palung Sangihe memanjang dari

Sulawesi Utara sampai selatan Mindanao, dan Palung Cotabato di bagian barat
 22

Mindanao. Zona penunjaman di daerah Maluku membentuk suatu permukaan

cembung disebabkan terdapat lempeng-lempeng yang masing-masing

menunjam ke arah barat dan kearah timur. Kecepatan penunjaman lempeng di

daerah ini 7 cm pertahun. Kedalaman penunjaman di sebelah barat mencapai 625

km dan disebelah timur 275 km. Sudut penunjamanan di sebelah barat adalah

32º– 51º sedangkan di timur antara 34º - 51º. Disamping itu terdapat beberapa

penunjaman lempeng pendek di bagian barat dan timur.

E. Tektonik Sumatera

Pulau Sumatera merupakan sebagian dari Lempeng Eurasia yang bergerak relatif

ke arah barat daya dan berinteraksi dengan Lempeng India-Australia yang terletak

di sebelah barat Pulau Sumatera yang bergerak relatif ke arah utara dengan

kecepatan 6 cm/tahun. Zona pertemuan antara kedua lempeng tersebut membentuk

zona subduksi. Berdasarkan analisis mekanisme sumber (focal mechanism)

kemiringan subduksinya antara 1º sampai 10º dengan dip dominan di bagian bawah

wilayah Sumatera (Ardiansyah, 2012).

Katili dalam Supartoyo, dkk., (2014) dalam penelitiannya menyebutkan bahwa

sistem busur subduksi Sumatera dibentuk oleh penyusupan Lempeng Benua.

Lempeng Benua tebal dan tua ini meliputi busur vulkanik, kapur, dan tersier.

Sedimen elastis sangat tebal menyusup di subduksi Sumatera dan sedimen yang

tebal didorong ke atas membentuk rangkaian kepulauan. Jalur subduksi ini

membujur sepanjang pantai barat Sumatera, tidak terkecuali pantai barat Bengkulu.

Hal ini pula yang menyebabkan di Daerah Bengkulu dan sebagian besar wilayah
 23

Sumatera rentan terhadap bahaya gempabumi. Peta tektonik kepulauan Indonesia

ditunjukan oleh Gambar 9.

Gambar 9. Peta tektonik Indonesia (Irsyam, dkk., 2010)

a. Zona Patahan di Sumatera

Akibat tumbukan lempeng terbentuklah patahan-patahan di Sumatera dan Sesar

Mentawai. Patahan di Sumatera dari Aceh hingga Teluk Semangko. Provinsi

Lampung yang dikenal dengan nama Sesar Besar Sumatera. Sedangkan Sesar

Mentawai terletak di laut, yaitu antara cekungan muka dan zona prismatik akresi

di sebelah barat Pulau Sumatera (Hidayati, dkk., 2010). Sesar Sumatera memiliki

aktivitas yang tinggi sementara Sesar Mentawai hanya sebagiannya saja yang

memiliki aktivitas yang cukup tinggi (Mustafa, 2010).

 24

Beberapa patahan aktif yang terdapat di Sumatera dan parameter gempanya

ditunjukan oleh Tabel 1 dan Gambar 10 sebagai berikut:

Tabel 1. Data dan parameter sumber gempa Sumatera Fault

Fault Slip-Rate
Sense Dip Top Bottom L Mmax
No Nama mm/yr weight Mechanism (km)

1 Aceh 2 1 Strike-slip 90 3 20 230 7,7

2 Seulimeun 2,5 1 Strike-slip 90 3 20 120 7,5

3 Tripa 6 1 Strike-slip 90 3 20 180 7,7

4 Renun 27 1 Strike-slip 90 3 20 220 7,8

5 Toru 24 1 Strike-slip 90 3 20 95 7,4

6 Angkola 19 1 Strike-slip 90 3 20 160 7,6

7 Barumun 4 1 Strike-slip 90 3 20 125 7,5

8 Sumpur 23 1 Strike-slip 90 3 20 35 6,9

9 Sianok 23 1 Strike-slip 90 3 20 90 7,3

10 Sumani 23 1 Strike-slip 90 3 20 60 7,2

11 Suliti 23 1 Strike-slip 90 3 20 95 7,4

13 Dikit 11 1 Strike-slip 90 3 20 60 7,2

14 Ketaun 11 1 Strike-slip 90 3 20 85 7,3

15 Musi 11 1 Strike-slip 90 3 20 70 7,2

16 Manna 11 1 Strike-slip 90 3 20 85 7,3

17 Kumering 11 1 Strike-slip 90 3 20 150 7,6

18 Semangko 5 1 Strike-slip 90 3 20 65 7,2

19 Sunda 5 1 Strike-slip 90 3 20 150 7,6

 25

Gambar 10. Sumatera Fault Zone (SFZ) (Mustafa, 2010)

b. Zona Subduksi Megathrust Sumatera

Zona subduksi Sumatera merupakan wilayah yang paling sering melepaskan

energi gempabumi. Dalam sejarah kegempabumian tercatat banyak gempabumi

yang terjadi dengan magnitudo di atas 8 SR. Di sebelah selatan khatulistiwa,

gempabumi besar pernah terjadi tahun 1833 (M8,9 SR) dan pada tahun 1797

(M8,3-8,7 SR). Kedua gempabumi ini menyebabkan tsunami besar yang

menyapu perairan Sumatera Barat dan Bengkulu. Wilayah zona subduksi di

selatan ini biasa dikenal dengan Segmen Mentawai. Pada bulan September 2007

segmen ini kembali melepaskan energinya sebesar Mw 8,4 (Setyonegoro, dkk.,

 26

2012). Zona subduksi dangkal di Sumatera yang terdiri dari empat zona yaitu :

Zona Subduksi Megathrust Andaman Sumatera, Zona Megathrust Mid-1 (Nias)

Sumatera, Zona Megathrust M2 (Siberut) Sumatera, dan Zona Megathrust

Southern Sumatera (Santoso dan Soehaemi, 2010).

Beberapa zona subduksi yang terdapat di Sumatera dan parameter gempanya

ditunjukan oleh Gambar 11 dan Tabel 2 sebagai berikut:

Gambar 11. Zona subduksi Megathrust Sumatera (Santoso dan Soehaemi,

2010)
 27

Tabel 2. Data dan parameter sumber gempa subduksi (Megathrust)

MMax
No Megathrust MMax History b-val a-val (Desain)
GR Char
1 Andaman-Sumatera 9,2 (26-12-2004) 0,826 4,69 8,0 9,2

2 Nias (Mid-1 Sumatera) 8,7 (28-03-2005) 0,878 4,71 8,7 8,7

3 Siberut (Mid-2 Sumatera) 8,5 (12-09-2007) 0,970 5,35 8,5 8,5

4 Southern Sumatera 7,9 (04-06-2000) 1,050 5,76 8,2 8,2

F. Gelombang Seismik

Gelombang seismik adalah rambatan energi yang disebabkan karena adanya

gangguan di dalam kerak bumi, misalnya adanya patahan atau adanya ledakan.

Energi ini akan merambat ke seluruh bagian bumi yang dapat terekam oleh

seismometer. Efek yang ditimbulkan oleh adanya gelombang seismik ini adalah apa

yang kita kenal sebagai fenomena gempabumi. Perambatan gelombang seismik

tergantung dari sifat elastisitas yang dimiliki oleh suatu batuan. Gelombang seismik

yang melalui bidang interior bumi disebut Gelombang Badan sedangkan

gelombang seismik yang merambat melalui luar bumi disebut Gelombang

Permukaan (Munadi, 2002).

a. Gelombang Badan (Body Wave)

Gelombang badan adalah gelombang seismik yang merambat hingga ke dalam

bumi. Berdasarkan gerak partikel pada media dan arah penjalarannya,

gelombang badan dapat dibedakan atas gelombang P (P-wave) dan gelombang

S (S-wave).
 28

 P-wave atau Gelombang Primer Gelombang P atau disebut juga gelombang

primer adalah gelombang yang pertama kali tercatat di alat seimometer. Hal

ini disebabkan karena gelombang ini memiliki kecepatan yang paling cepat

dibandingkan gelombang seismik yang lain. Bentuk penjalaran gelombang P

adalah longitudinal dimana arah gelombang sejajar dengan arah

penjalarannya. Gelombang P dapat menjalar di semua medium. Arah

getarannya kedepan dan kebelakang sehingga materi yang dilaluinya

mengalami tekanan dan peragangan seperti spiral. Oleh karena itu, sering

disebut dengan Push-Pull Wave atau Compressional Wave. Persamaan dari

kecepatan gelombang P adalah,

𝜆+2µ
𝑉𝑝 = 𝛼 = √ (1)
𝜌

dimana 𝑉𝑝 adalah kecepatan gelombang P, 𝜆 adalah parameter Lame, 𝜇

adalah modulus geser dan adalah densitas batuan.

 S-wave atau Gelombang Sekunder

Gelombang S atau gelombang sekunder adalah gelombang seismik yang

hanya merambat di permukaan bumi. Merupakan gelombang transversal yang

memiliki arah getar tegak lurus dengan arah penjalarannya. Gelombang ini

memiliki waktu perambatan yang lebih lama dari pada gelombang P sehingga

akan tercatat setelah gelombang P pada alat seismometer. Gelombang S tidak

dapat merambat di medium cair. Persamaan dari kecepatan gelombang S

adalah

µ
𝑉𝑝 = 𝛽 = √𝜌 (2)
 29

dimana 𝑉𝑠 adalah kecepatan gelombang S, adalah modulus geser dan

adalah densitas batuan.

Simulasi pergerakan gelombang badan dapat dilihat pada Gambar 12.

Gelombang ini menyebabkan gerakan partikel–partikel media dalam arah

tangensial terhadap arah penjalaran gelombang. Bila arah getar gelombang S

terpolarisir pada bidang vertikal maka gelombang tipe ini disebut gelombang

SV. Sedangkan bila arah getarnya terpolarisir pada bidang horisontal maka

gelombang ini disebut gelombang SH (Munadi, 2002).

Gambar 12. Simulasi gerakan gelombang P dan S (Munadi, 2002)

b. Gelombang Permukaan (Surface wave)

Gelombang permukaan adalah jenis gelombang seismik yang hanya merambat di

permukaan bumi. Amplitudo gelombang ini akan semakin melemah jika semakin

masuk ke dalam bumi. Gelombang ini dapat disamakan dengan gelombang air yang

mengalir diatas permukaan bumi, gerakannya lebih lambat dibandingkan

 30

gelombang badan. Ada dua tipe Gelombang Permukaan, yaitu gelombang Love dan

Gelombang Rayleigh.

 Gelombang Rayleigh

Gelombang Rayleigh adalah gelombang yang merambat pada batas permukaan

saja dan hanya dapat merambat pada media padat serta arah getarannya

berlawanan arah dengan arah perambatannya. Gelombang Rayleigh memiliki

gerakan partikel yang merupakan kombinasi dari gerakan partikel Gelombang P

dan S. Gerakan partikel gelombang ini terpolarisasi elips dengan faktor

amplitude yang mempunyai tanda berkebalikan sehingga gerakan partikelnya

mundur (Santoso, 2004).

 Gelombang Love

Gelombang love adalah gelombang yang hanya merambat pada batas lapisan

saja dan bergerak pada bidang yang horisontal saja (Putri, 2012). Gelombang

Love adalah Gelombang Permukaan yang gerakan partikelnya mirip dengan

Gelombang S, yaitu terjadi secara transversal. Gelombang Love merupakan

gelombang permukaan yang terbesar dan amplitudonya meningkat seiring

dengan bertambahnya kedalaman.

Gelombang Love merambat lebih cepat dibandingkan Gelombang Rayleigh.

Simulasi perambatan gelombang permukaan dapat dilihat pada Gambar 13.

 31

Gambar 13. Gelombang love dan gelombang Rayleigh (Munadi, 2002)

G. Parameter Sumber Gempa

Setiap kejadian gempabumi akan menghasilkan informasi seismik berupa rekaman

sinyal berbentuk gelombang yang setelah melalui proses manual atau non manual

akan menjadi data bacaan fase. Informasi seismik selanjutnya mengalami proses

pengumpulan, pengolahan dan analisis sehingga menjadi parameter gempabumi.

Parameter gempabumi tersebut meliputi :

a. Waktu terjadinya gempa (Origin time)

Origin time atau waktu terjadinya gempabumi merupakan waktu dimana

pelepasan energi pertama kali terjadi pada lempeng tektonik bumi yang

mengalami tekanan akibat tumbukan atau gesekan dan dinyatakan dalam hari,

tanggal, bulan, tahun, jam, menit, detik dalam satuan UTC (Universal Time

Coordinated). Pada umumnya, Origin time ditentukan dari perpotongan garis Tp

ketika Ts – Tp sama dengan nol yang secara sederhana dapat menggunakan

Diagram Wadati yang ditunjukkan pada Gambar 14.

 32

𝑣𝑝
tan 𝜃 = 𝑖 = −1
𝑣𝑠

Gambar 14. Diagram Wadati (Hurukawa, dkk., 2008)

Dimana Tp adalah waktu tiba gelombang P, Ts adalah waktu tiba gelombang S,

To adalah origin time, Vp adalah kecepatan gelombang P dan Vs adalah

kecepatan gelombang S. Pada diagram wadati ini tidak memperhatikan model

struktur dalam bumi, sehingga hanya mengasumsikan struktur bumi yang

bersifat homogen. Hasil origin time ini menunjukkan hiposenter yang kurang

akurat sehingga harus dilakukan relokasi ulang untuk mengetahui hiposenter

yang lebih akurat dengan struktur bumi yang heterogen. Umumnya, di BMKG

menggunakan pendekatan SED (Single Event Determination) yang

mengasumsikan origin time dekat dengan stasiun pencatat yang kemudian

dikembangkan menjadi metode MJHD (Modified Joint Hypocenter

Determination).

b. Hiposenter

Hiposenter merupakan pusat gempabumi yang berada di dalam permukaan

bumi. Untuk memudahkan terkadang hiposenter diasumsikan sebagai sebuah

titik, namun pada kenyataannya hiposenter merupakan sebuah bidang yang

luasnya tergantung pada besarnya energi yang dilepaskan. Penentuan hiposenter

 33

juga dapat menggunakan diagram Wadati, dengan mengasumsikan bahwa

lapisan bumi adalah homogen yang ditunjukan pada Gambar 15.

Gambar 15. Jarak hiposenter (Hurukawa, dkk., 2008)

Dari Gambar 15, D adalah jarak hiposenter dengan stasiun pencatat, dan dapat

dinyatakan dengan persamaan berikut:

𝐷 = 𝑇𝑝𝑜 ∗ 𝑉𝑝

𝐷 = 𝑇𝑠𝑜 ∗ 𝑉𝑠 = (𝑇𝑠 – 𝑇o) ∗ 𝑉𝑠

= {(𝑇𝑠 − 𝑇𝑝) + (𝑇𝑝 – 𝑇o)} ∗ 𝑉𝑠 (3)

= (𝑇𝑠𝑝 + 𝑇𝑝o) ∗ 𝑉𝑠

dimana Tsp adalah Ts – Tp.

Ts merupakan waktu tiba gelombang S, Tp merupakan waktu tiba gelombang

P, Vp merupakan kecepatan gelombang P dan Vs merupakan kecepatan dari

gelombang S. Dari persamaan 3 diatas, kita dapat mengetahui jarak dari stasiun

ke pusat gempabumi.

c. Episenter

Episenter merupakan titik di permukaan bumi yang merupakan reflektifitas

tegak lurus dari hiposenter. Lokasi episenter dibuat dalam koordinat kartesian
 34

bola bumi atau sistem koordinat geografis dan dinyatakan dalam derajat lintang

dan bujur.

d. Magnitudo

Ukuran dari kekuatan gempa disebut magnitudo, yaitu parameter gempa yang

mengukur besarnya energi gempa yang dilepaskan dari sumbernya. Jadi

pengukuran magnitudo yang dilakukan di tempat yang berbeda harus

menghasilkan harga yang sama walaupun gempa yang dirasakan di

tempattempat tersebut tentu berbeda. Satuan yang dipakai adalah Skala Richter.

H. (Modified Joint Hypocenter Determination) MJHD

Metode (Modified Joint Hypocenter Determination) MJHD yang dikembangkan

oleh Hurukawa dan Imoto (1992) merupakan pengembangan dari metode Joint

Hypocenter Determination (JHD) yang sebelumnya dikembangkan oleh Douglas

(1967). Prinsip JHD dan MJHD adalah secara simultan menginversikan waktu

tempuh dari sekelompok data gempabumi untuk mendapatkan lokasi hiposenter

yang lebih baik dibandingkan dengan metode single event determination (SED).

Pada kedua metode ini dimasukkan faktor koreksi stasiun yang bertujuan untuk

untuk mereduksi kesalahan akibat adanya variasi kecepatan secara lateral

(Hurukawa dan Imoto, 1992).

Pada metode SED standar, umumnya digunakan model kecepatan 1-D yang harus

ditentukan sebelumnya. Pemilihan model kecepatan 1-D ini umumnya didasarkan

pada model kecepatan yang diperoleh dari hasil pemodelan sebelumnya atau jika
 35

tidak tersedia biasanya merujuk pada model global. Alternatif pemilihan model cara

kedua yang lebih sering digunakan dalam metode SED standar.

Metode ini dikembangkan oleh Geiger (1910) yang merupakan iterasi numerik

dengan optimasi Gauss-Newton. Metode SED merupakan suatu metode

pengolahan data gempa mikro yang digunakan untuk menentukan hiposenter.

Hiposenter merupakan lokasi fisik berdasarkan koordinat lintang, bujur dan

kedalaman tempat terjadinya gempa. Metode SED melakukan iterasi minimum

dengan teori Geiger Adaptive Damping (GAD).

Berdasarkan propagasi gelombang dibedakan menjadi dua jenis yaitu gelombang

P dan gelombang S. Gelombang P memiliki waktu tempuh yang lebih singkat

dibandingkan dengan gelombang S sehingga memiliki residual waktu kedatangan

kedua gelombang. Langkah pertama adalah menebak hiposenter dan origin time

(x0, y0, z0, t0). Dalam kasus event dekat atau dalam jaringan stasiun, ini dapat

diselesaikan menggunakan lokasi dekat stasiun dengan waktu tiba pertama dan

menggunakan waktu tiba tersebut sebagai t0. Diasumsikan bahwa hiposenter

sebenarnya cukup dekat dengan nilai tebakan sehingga waktu tempuh residual itu

di hiposenter percobaan adalah fungsi linier dari koreksi yang kita punya untuk

membuat jarak hiposenter.

Dimulai dari Tobs adalah waktu tiba pertama gelombang seismik di setiap stasiun

ke-i (xi, yi, zi) dari hiposenter, Tcal adalah waktu tempuh kalkulasi berdasarkan

model kecepatan 1 dimensi bawah permukaan. Waktu residual rij untuk stasiun i
 36

adalah selisih antara waktu tiba observasi dan waktu tiba kalkulasi yang secara

matematis dirumuskan sebagai berikut:

rij = 𝑇𝑖𝑗𝑜𝑏𝑠 − 𝑇𝑖𝑗𝑐𝑎𝑙 (4)

Pendekatan yang paling umum untuk menggunakan solusi least squares untuk

menemukan jumlah residual kuadrat minimum e dari n observasi:

𝑒 = ∑𝑛𝑖=1(𝑟𝑖 )2 (4)

Lokasi hiposenter akan menjadi titik yang paling sesuai antara waktu kalkulasi dan

observasi dengan cara e terkecil. Residual root mean squared (RMS) digambarkan

𝑒
sebagai √ . RMS diberikan dalam hampir semua program lokasi dan umumnya
𝑛

digunakan sebagai panduan untuk akurasi lokasi.

Selanjutnya residual dapat dituliskan:

𝑟𝑖 = (𝑡𝑖 − 𝑡0 ) − 𝑇𝑖 (5)

𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗

𝑟𝑖 = ∆𝑥 + ∆𝑦 + ∆𝑧 + ∆𝑡 (6)
𝜕𝑥𝑖 𝜕𝑦𝑖 𝜕𝑧𝑖

Dengan,

𝑇𝑖𝑗 = waktu tempuh kalkulasi di stasiun i

√(𝑥−𝑥𝑖)2 +(𝑦−𝑦𝑖)2 +(𝑧−𝑧𝑖)2

𝑇𝑖𝑗 = (7)
𝑣

Persamaan yang sama juga digunakan untuk y dan z.

Dalam bentuk matriks kita dapat menuliskan:

J ∆m = ∆d (8)

[JT] [J] ∆m = [J] T ∆d (9)

[JTJ] -1 JT ∆m = [JTJ] -1 JT ∆d (10)

[∆m]= [JTJ] -1 JT ∆d (11)

Dengan elemen penyusun matriks sebagai berikut:

 37

𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗

𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡
… … … …
𝐽= … … … … (12)
𝜕𝑇𝑖𝑗𝑛 𝜕𝑇𝑖𝑗 𝑛 𝜕𝑇𝑖𝑗 𝑛 𝜕𝑇𝑖𝑗 𝑛
[ 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 ]

∆𝑥
∆𝑦
∆m = [ ] (13)
∆𝑧
∆𝑡
𝑟1
…
∆d = [ … ] (14)
𝑟𝑛

Matriks J berupa matriks kernel (jacobian), berisi derivative parsial residual waktu

tempuh tiap stasiun terhadap parameter hiposenter berukuran n x 4, dengan n

adalah jumlah stasiun pengamat dan 4 adalah representasi jumlah parameter

hiposenter yang dicari. Matriks ∆m adalah model yang diinginkan dan ∆d

merupakan matriks residual berisikan parameter residual waktu tiba yang diperoleh

dari data. Sehingga didapatkan perkiraan koordinat setelah dikoreksi, diberikan

oleh persamaan:

X0 + ∆x = x

Y0 + ∆y = y (15)

Z0 + ∆z = z

` t0 + ∆t = t

Dengan,

x0, y0, z0, t0 : lokasi dan origin time dugaan

∆x, ∆y, ∆z, ∆t : selisih hiposenter dugaan hasil kalkulasi

x, y, z, t : lokasi dan origin time hasil kalkulasi

 38

Selanjutnya Metode Geiger ini dikembangkan menjadi persamaan untuk metode

JHD yang dikembangkan oleh Douglas (1967) dan Fredman (1967) diberikan

penambahan besaran koreksi stasiun untuk mereduksi kesalahan akibat variasi

kecepatan lateral yang tidak diperhitungkan dalam model kecepatan 1D (Pujol,

2000).

Penambahan besaran koreksi stasiun menghasilkan formulasi residu waktu tempuh

yang diperoleh pada stasiun ke-i sebagai berikut :

𝑟𝑖𝑗 = 𝑇𝑖𝑗𝑜𝑏𝑠 − 𝑇𝑖𝑗𝑐𝑎𝑙 + 𝑠𝑖 (16)

Tobs adalah waktu tempuh gelombang seismik dari pusat gempa bumi ke stasiun

yang diperoleh dari selisih waktu tiba gelombang gempa bumi dengan waktu

kejadian (origin time), Tcal adalah waktu tempuh kalkulasi dan si adalah koreksi

stasiun. Dengan menggunakan deret Taylor, persamaan (16) dapat diuraikan

menjadi :

𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗

𝑑𝑟𝑖 = 𝑑𝑥𝑗 + 𝑑𝑦𝑗 + 𝑑𝑧𝑗 + 𝑑𝑇𝑜𝑗 + 𝑑𝑆𝑖 (17)
𝜕𝑥𝑗 𝜕𝑦𝑗 𝜕𝑧𝑗

dimana 𝑑𝑥𝑗 , 𝑑𝑦𝑗 , 𝑑𝑧𝑗 adalah koreksi untuk perkiraan awal hiposenter dan 𝑑𝑇𝑜𝑗

adalah origin time untuk gempa ke-j, sementara 𝑑𝑆𝑖 adalah koreksi untuk stasiun

𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗 𝜕𝑇𝑖𝑗

ke–i. Koefisien , , dihitung berdasarkan model kecepatan yang
𝜕𝑥𝑗 𝜕𝑦𝑗 𝜕𝑧𝑗

digunakan.

Jika persamaan (17) disusun menjadi sebuah matriks untuk semua stasiun maka

persamaan matriksnya adalah sebagai berikut:

 39

(18)

𝑟𝑗 berisi residual waktu tempuh tiap gempa, 𝐴𝑗 adalah matriks yang berisi partial

derivative dari residual waktu tempuh terhadap parameter hiposenter, 𝑑𝑥𝑗 berisi

vektor pertubasi atau perubahan parameter hiposenter terhadap parameter dugaan

(awal), dan ds adalah koreksi stasiun.

Nilai dx, dy, dz, dTo, dan dS merupakan perturbasi parameter-parameter model yang

ingin ditentukan, dalam hal ini dihimpun dalam suatu vektor m. Vektor m diperoleh

menggunakan metode optimasi kuadrat terkecil (least square atau disingkat LSQ)

dengan meminimalkan suatu fungsi objektif berupa nilai kuadrat dari residual :

𝑓(𝑚) = ∑(𝑂 − 𝐶)2 → 𝑚𝑖𝑛𝑖𝑚𝑢𝑚 (19)

Dalam metoda LSQ dibutuhkan model inisial sebagai perkiraan awal, dalam hal ini

digunakan hasil dari metoda SED. O merupakan waktu observasi dan C merupakan

waktu hasil dari kalkulasi. Setelah diperoleh nilai dx, dy, dz, dan dTo, nilai model

inisial diperbarui menjadi parameter model yang baru : x1 = x0 + dx, y1 = y0 + dy,

z1 = z0 + dz, dan t1 = t0 + dTo. Penyelesaian ini dilakukan secara iteratif dengan

mengasumsikan model baru yang diperoleh sebagai model inisial. Proses iteratif ini

dilakukan hingga perubahan yang diperoleh sudah tidak signifikan lagi atau dibatasi

oleh jumlah iterasi maksimum.

 40

Metode JHD ini digunakan untuk memperbaiki lokasi gempabumi dengan

memperhitungkan adanya kesalahan dari model bumi yang digunakan sebelumnya.

Sehingga metode ini mampu merelokasi suatu kumpulan gempabumi secara

simultan menggunakan koreksi stasiun. Koreksi ini dihitung untuk setiap stasiun

pengamat yang merekam kumpulan kejadian gempa, sehingga inversi JHD bisa

menghasilkan lokasi hiposenter yang relatif lebih baik.

Namun, apabila gelombang gempabumi melewati medium yang sangat heterogen

dan sebaran stasiunnya kurang baik, maka solusi JHD menjadi tidak stabil dan

kurang dapat diandalkan. Hal ini disebabkan adanya trade – off (loosing quality)

antara nilai koreksi stasiun dengan kedalaman fokus gempabumi (Hurukawa, dkk.,

2008). Dengan memodifikasi metode JHD, Hurukawa dan Imoto (1992)

mengembangkan metode Modified Joint Hypocenter Determination (MJHD)

dengan menambahkan beberapa constrain sebagai berikut:

𝑛

∑ 𝑆𝑖 𝐷𝑖 = 0
𝑖=1

∑ 𝑆𝑖 = 0
𝑖=1
(20)
𝑛

∑ 𝑆𝑖 cos 𝜃𝑖 = 0
𝑖=1

∑ 𝑆𝑖 sin 𝜃𝑖 = 0
𝑖=1

Dimana Si adalah koreksi stasiun pada stasiun ke – i, Di adalah jarak antara stasiun

ke – i dengan pusat cluster, 𝜃𝑖 adalah azimut dari stasiun ke - i terhadap pusat

cluster, dan n adalah nomor stasiun. Pemberian constrain ini berdampak pada nilai
 41

koreksi stasiun yang tidak bergantung pada jarak dan azimuth dari pusat cluster ke

stasiun yang ilustrasinya dapat dilihat pada Gambar 16.

Gambar 16. Ilustrasi relokasi menggunakan MJHD (Hurukawa dan Imoto, 1992)

Meskipun metode ini mengakibatkan perubahan pada lokasi absolut gempabumi,

namun solusi dari MJHD ini menjadi lebih stabil yang dapat kita lihat pada Gambar

16 diatas. Apabila gempabumi yang terjadi dalam, maka solusi menggunakan

metode MJHD akan lebih stabil dibandingkan dengan menggunakan metode JHD.

Keunggulan metode MJHD ini adalah tidak diharuskannya memilih master event,

dimana sangat efektif pada saat terjadi kasus gempabumi yang tidak diamati dengan

jelas di semua stasiun (Hurukawa, dkk., 2008)

MJHD merupakan metode untuk merelokasi posisi pusat gempabumi (lintang,

bujur, kedalaman, dan origin time). Perubahan yang signifikan pada umumnya

terjadi pada parameter kedalaman hiposenter. Relokasi pusat gempa menggunakan

metode MJHD mempunyai kelebihaan yaitu karena adanya penambahan koreksi

stasiun, sehingga meskipun model kecepatan yang digunakan sama seperti yang

digunakan untuk mendapatkan hiposenter awal, namun hasilnya akan berbeda dan

menjadi lebih akurat. Demikian pula dengan adanya penambahan faktor azimuth
 42

dan jarak relatif stasiun membuat metode ini bisa mengadaptasi variasi kecepatan

lateral yang cukup heterogen serta distribusi stasiun yang umumnya tidak merata

(Aswad, dkk., 2012).

Dalam metode MJHD ini menggunakan model kecepatan global IASP91 yang

merupakan model kecepatan bumi satu dimensi yang dihasilkan oleh International

Association of Seismology and Physics of the Earth’s Interior (IASPEI) yang

merupakan hasil pemodelan dengan menggunakan ribuan data gelombang P dan S

yang direkam oleh ribuan seismometer di seluruh dunia. Model IASP91 ini hanya

berisikan parameter kecepatan gelombang P dan S. Pada kecepatan IASP91 ini

memerhitungkan dua lapisan diskontinuitas pada lapisan kerak bumi dikedalaman

20 sampai 365 km, laisan diskontinuitas dimantl bumi pada kedalaman 410 sampai

660 km. Berikut merupakan model kecepatan IASP91 yang digunakan terdapat

pada Gambar 17.

Gambar 17. Model kecepatan IASP91 (Putri, 2012)

 43

Kelebihan dalam relokasi menggunakan metode MJHD ini adalah adanya

penambahan koreksi stasiun pada proses penentuan ulang hiposenter gempabumi.

Selain itu juga dengan menggunakan kecepatan IASP91 akan menghasilkan

hiposenter yang lebih akurat karena asumsi dari kecepatan yang digunakan berupa

kecepatan bumi global yang mengasumsikan struktur bumi yang heterogen.

I. Menentukan Kelengkungan dan Sudut Penunjaman

Kelengkungan penunjaman ditentukan dengan menggunakan metode regresi.

Regresi merupakan hubungan antara dua variabel atau lebih yang dinyatakan dalam

bentuk fungsi. Regresi yang digunakan pada penulisan ini yaitu regresi polinomial

yang merupakan sebuah variabel bebas dengan pangkat terurut. Orde yang

digunakan yaitu orde tiga yang mempunyai korelasi terbaik dari data hiposenter

gempa yang digunakan. Berikut bentuk umum dari persamaan orde tiga:

𝑦 = 𝑎𝑥 3 + 𝑏𝑥 2 + 𝑐𝑥 + 𝑑 (21)

Dimana

Y = fungsi kedalaman (km)

X = jarak pusat gempa dari penunjaman terdangkal (km)

a,b,c,d = konstanta

Besar sudut penunjaman dihitung dengan menggunakan garis singgung kurva pada

salah satu titik di sepanjang kurva. Titik tersebut merupakan titik awal terjadinya

perubahan sudut yang signifikan. Setelah menentukan garis singgung kurva, maka

sudut penunjaman dapat ditentukan. Ilustrasi penentuan sudut dapat dilihat pada

Gambar 18 berikut ini:

 44

Gambar 18. Ilustrasi garis singgung kurva (Listyaningrum, 2015)

Hubungan dari garis singgung kurva tersebut dapat dinyatakan dalam persamaan:

𝑦 = 𝑓(𝑥) (22)

Kemiringan garis singgung atau gradien dapat dinyatakan dengan turunan pertama

dari kurva y=f(x)

𝑚 = 𝑓′(𝑥) (23)

Kemudian untuk menentukan persamaan garis singgung dapat dinyatakan dengan:

𝑦 − 𝑦1 = 𝑚(𝑥 − 𝑥1 ) (24)

Besar sudut (α) didapatkan dengan menentukan terlebih dahulu titik 𝑥1 ′ pada garis

singgung saat y=0.

𝑦
𝛼 = arctan(𝑥 1′) (25)
1
 IV. METODE PENELITIAN

A. Lokasi dan Tempat Penelitian

Tugas Akhir dilakukan di “Badan Meteorologi Klimatologi dan Geofisika” dengan

waktu pelaksanaan selama 2 bulan. Tema dari Tugas Akhir ini adalah “Relokasi
Hiposenter Gempabumi Menggunakan Metode MJHD (Modified Joint
Hyopocenter Determination) Untuk Analisis Zona Subduksi Sumatera Bagian
Selatan”.

B. Data Penelitian dan Perangkat Lunak

Data dan perangkat lunak yang digunakan dalam penelitian ini adalah:

1. Data arrival time gelombang P dan gelombang S periode waktu Januari 2010

s.d Desember 2016 untuk wilayah Sumatera bagian Selatan dengan batasan

wilayah -3.5º s.d -9º LS dan 99º s.d 106.5º BT yang diperoleh dari katalog

BMKG (Badan Meteorologi Klimatologi dan Geofisika)

2. Software Cygwin untuk melakukan relokasi data gempabumi dengan

menggunakan metode MJHD

3. Software GMT untuk melakukan plotting peta sebelum relokasi maupun hasil

dari relokasi

4. Software Excel untuk menentukan pola dan sudut zona subduksi daerah

penelitian
 46

C. Pengolahan Data

Dalam penelitian ini, dilakukan beberapa tahap pengolahan data sebagai berikut:

1. Mengubah data gempabumi dari katalog BMKG dengan format *.txt yang telah

diubah terlebih dahulu menjadi format *.data kedalam format MJHD agar bisa

dilakukan relokasi. Data awal yang didapat dari katalog BMKG berupa data

dalam format arrival time gelombang P dan S yang dapat dilihat pada Gambar

19.

Gambar 19. Data dari katalog BMKG

Data tersebut kemudian diubah kedalam format MJHD yang dapat dilihat pada

Gambar 20.

Gambar 20. Data dalam format MJHD

 47

Pada format data MJHD ini, parameter yang diambil adalah origin time,

lintang, bujur, kedalaman, magnitudo dan jumlah fase dari setiap event

gempabumi yang terjadi, serta keterangan tentang stasiun pencatat.

2. Menentukan nilai Minimum Number of Earthquake (MEQ) dan Minimum

Number of Station (MNST) yang ditunjukkan dengan Gambar 21.

Gambar 21. Tampilan input program station

Nilai MEQ dan MNST digunakan sebagai input pada program station dengan

bahasa fortran yang merupakan input parameter dalam program MJHD. Nilai

MEQ yang digunakan pada penelitian ini yaitu 10 dan MNST yang digunakan

yaitu 5. MEQ merupakan jumlah minimal gempabumi yang dicatat oleh satu

stasiun, sedangkan MNST merupakan jumlah minimal stasiun yang dipakai

untuk satu gempabumi.

3. Melakukan inversi menggunakan MJHD dengan melakukan beberapa iterasi

untuk mendapatkan hasil hiposenter yang lebih akurat. Berikut merupakan

tampilan dari iterasi yang digunakan pada software Cygwin yang ditunjukan

pada Gambar 22.

Gambar 22. Tampilan parameter iterasi yang digunakan

 48

Parameter dari inputan iterasi tersebut berupa iform atau tipe data input, nama

file masukan, nama file keluaran, longitude dan latitude kedalaman fix (ZFIX)

yaitu batas kedalaman gempabumi dalam km, residu maksimum travel time

(RESS), jumlah iterasi maksimum (ITRT), standar deviasi (STD1), jumlah

stasiun yang tidak digunakan (NAST), jumlah gempabumi pada data yang tidak

digunakan (NEXC), akurasi pembacaan (RANKAB), RMAX yaitu batasan

derajat dimana jika sebuah stasiun ada dalam batasan RMAX dari pusat area

maka stasiun tersebut akan dilibatkan dalam perhitungan, pilihan mengenai

hasil yang akan dicetak, magnitude minimum (AMGM), dan nilai SLOPE.

4. Melakukan plotting dengan GMT untuk mengetahui sebaran hiposenter

sebelum dan sesudah direlokasi. Berikut merupakan sebaran event gempabumi

yang digunakan dalam penelitian ini yang ditunjukan pada Gambar 23.

Gambar 23. Contoh plotting gempabumi dengan GMT

5. Menentukan zona subduksi Sumatera bagian Selatan dan mengetahui sudut

penunjamannya dengan menggunakan Ms. Excel.

 49

D. Diagram Alir

Diagram alir yang digunakan dalam penelitian ini terdapat pada Gambar 24.

Mulai

Arrival time
gelombang P dan S

Ubah format data

Data dalam
format MJHD

Inversi MJHD

Tidak

RMS
<1

Ya

Hiposenter Akhir

Plotting hasil dengan GMT

Pola Zona Subduksi Sumatera Bagian Selatan

Selesai

Gambar 24. Diagram alir

 50

E. Time Schedule

Time Schedule dari penyusunan skripsi ini tertera pada Tabel 3.

Tabel 3. Time schedule penyusunan skripsi

Jadwal kegiatan
Kegiatan
Februari Maret April Mei Juni
Studi literatur
Pemilihan data
Pemilihan software
Pengolahan data
Penyusunan laporan
Seminar proposal
Seminar hasil
Uji komprehensif
 BAB VI. KESIMPULAN

A. Kesimpulan

Dari hasil penelitian mengenai relokasi hiposenter gempabumi didapatkan beberapa

kesimpulan antara lain:

1. Relokasi hiposenter gempabumi menggunakan metode MJHD ini menunjukkan

peningkatan kualitas yang lebih baik dilihat dari nilai RMS yang kurang dari 1

yaitu antara 0.2 s.d 0.5.

2. Relokasi menggunakan MJHD menghasilkan distribusi hiposenter yang lebih

dominan mengikuti pola penunjaman dengan distribusi gempa yang semakin

dalam kearah timur.

3. Pola penunjaman yang dihasilkan menunjukan sudut zona subduksi Bengkulu

sebesar 26.78º, Lampung sebesar 30.22º dan Selat Sunda sebesar 52.53º.

4. Zona Seismik Ganda ditemukan sejajar dengan pola penunjaman didaerah

Lampung dengan jarak 50 km dengan kedalaman antara 50 – 150 km.

5. Kedalaman slab daerah Bengkulu berkisar 250 km sedangkan kedalaman slab

daerah Lampung dan Selat Sunda berkisar 400 km.

 74

B. Saran

Saran yang dapat diberikan yaitu perlunya penelitian lanjutan dengan menggunakan

metode lain seperti metode Hypo DD untuk memperkuat hasil dari penelitian ini.
 DAFTAR PUSTAKA

Ardiansyah, S., 2012, Eartquake Potential Energy in the Musi Segment,

Kepahiang-Bengkulu Area, Stasiun Geofisika Kepahiang.

Aswad, S., Ahmad, D., dan Budiati, M.R., 2012, Relokasi Gempa di Sepanjang
Sesar Palu Koro Menggunakan Metode Modified Joint Hyipocenter
Deterministik dan Double Difference. Prodi Geofisika, Jurusan Geofisika,
Fakultas MIPA Universitas Hasanuddin.

Awaludin, A., 2011, Penentuan Waktu Berakhirnya Gempa Susulan Untuk

Gempabumi Biak 16 Juni 2010, Skripsi, Program Studi Fisika, Fakultas Sains
dan Teknologi, Universitas Islam Negeri Syarif Hidayatullah Jakarta.

BMKG Kotabumi, 2016, Peran dan Kapasitas BMKG Kotabumi Dalam

Mendukung Sistem Peringatan Dini, Presentation.

Douglas, A., 1967, Joint Epicenter Determination, Nature, 215 : 47 – 48.

Fauzi, A., 2015, Relokasi Hiposenter Gempabumi Daerah Jawa dan Sekitarnya
Tahun 2009-2014 Menggunakan Metode MJHD, Skripsi, ITB.

Fredman, H.W., 1967, A Statistical Discussion of P Residual From Explosion, Part

II, Bulletin of The Seismological Soecity of America, 57 : 545-546.

Hasegawa, A., Umino, N., dan Takagi, A., 1978. Double-Planed Structure of th Dee
Seismic Zone in the Notheastern Japan Arc. Tectonophysics.

Hidayati, S., Sumaryono, dan Eka, S., 2010, Tsunami Mentawai 25 Oktober 2010,
Bulletin Vulkanologi dan Bencana Geologi, Vol 5 No 3 : 1-11.

Hurukawa, N., dan Imoto, M., 1992, Subducting Oceanic Crust of the Philippine
Sea and Pacific Plates and Weak-Zone-Normal Compression in Kanto
District, Japan, Geophys. J. Int., 109: 639652.

Hurukawa, N., Popa, M., dan Radulian, M. 2008. Relocation of Large

IntermediateDepth Earthquakes in The Vrancea Region, Romania, Since
1934 and a Seismic Gap. Earth, Planets and Space. 60(6):565-572.
Irsyam, M., Sengara, W., Aldiamar, F., Widiantoro, S., Triyoso, W., Hilman, D.,
Kertapati, E., Meilano, I., Suhardjono, Asrusifak, dan Ridwan, M., 2010,
Ringkasan Hasil Studi Tim Revisi Peta Gempabumi Indonesia 2010,
Kementerian Pekerjaan Umum.

Kaluku, A., 2015, Relokasi Gempabumi di Wilayah Selat Sunda Menggunakan

Metode Double Different (hypo DD), Skripsi, Geofisika, STMKG,
Tangerang Selatan.

Listyaningrum, Z., 2015, Studi Hubungan Antara Sudut Penunjaman Lempeng

Indo-Australia - Eurasia Terhadap Tatanan Tektonik Overriding Plate,
Skripsi, Geofisika, STMKG, Tangerang Selatan.

Madlazim, 2013, Kajian Awal Tentang b Value Gempabumi di Sumatera Tahun

1964-2013, Jurnal Penelitian Fisika dan Aplikasinya (JPFA).

Malik, 2009, Analisa Pola Subduksi Daerah Bengkulu dengan Metode Segmen
Irisan Vertikal, Skripsi, Akademi Meteorologi dan Geofisika.

Mangga, S.A., Amirudin, T., Suwarti, S., Gafoer dan Sidarto, 1993, Peta Geologi
Tanjungkarang, Sumatera, Pusat Penelitian dan Pengembangan Geologi,
Bandung.

Maung, P.M., 2009, Relocation of Earthquakes in Myanmar by MJHD Method:

Aftershocks of Large Earthquakes and Seismicity Along the Sagaing
Fault, IISEE, Tsukuba, Ibaraki, Japan.

Munadi, S, 2002, Pengolahan Data Seismik Prinsip Dasar dan Metodologi,

Program Studi Geofisika, Jurusan Fisika, Universitas Indonesia: Depok.

Mustafa, B., 2010, Analisis Gempa Nias dan Gempa Sumatera Barat dan
Kesamaannya yang Tidak Menimbulkan Tsunami, Jurnal Ilmu Fisika (JIF),
Vol 2 No 1.

Pesicek, J.D., Thurber, C.H., Zhang, H., Engdahl, .R., dan Widiyantoro, S.
Teleseismic Doubl-Different Rlocation Arthuakes Along the Sumatera-
Andan Subduction Zone Using 3D Model. Journal Geophysical Research.
Vol 115.

Pujol, J., 2000, Joint Event Location – The JHD Technique and Applications to
Data From Local Seismic Networks, Advances in Seismic Location, 163-204.

Putri, Y.T., 2012, Relokasi Gempabumi Utama Dan Gempabumi Susulan

Menggunakan Metode MJHD (Studi Kasus Gempabumi Mentawai 25
Oktober 2010), Skripsi, Depok: Universitas Indonesia.

Santoso, D., 2004, Volkanologi. Bandung: Institut Teknologi Bandung Press.

Santoso dan Soehaemi, A., 2010, Analisis Bahaya Gempabumumi Lengan Utara
Sulawesi, Pusat Survei Geologi, Vol 20.

Setyonegoro, W., Sunardi, B., Sulastri, Nugraha, J., dan Susilanto, P., 2012,
Analisis Sumber Gempabumi pada Segmen Mentawai (Studi Kasus:
Gempabumi 25 Oktober 2010), Jurnal Meteorologi dan Geofisika (JMG),
Vol 13 No 2.

Sieh, 2000, Nanotectonics of the Sumatran fault, Indonesia, Journal of geophysical

research, Vol 105

Sieh, K., dan Natawidjaja, D.H., 2000, Neotectonics of th Sumateran Fault,

Indonesia. Journal of Gophysical Research. Vol 105.

Subardjo, dan Ibrahim, G., 2004, Pengetahuan Seismologi, Jakarta: Badan

Meteorologi dan Geofisika.

Sukamto, R., Ratman, N., dan Simandjuntak, T.O., 1996, Peta Geologi Indonesia,
Pusat Penelitian dan Pengembangan Geologi.

Supartoyo, Surono, dan Putranto, E.T., 2014, Katalog Gempabumi Merusak di

Indonesia Tahun 1612-2014, Pusat Vulkanologi dan Mitigasi Bencana
Geologi.