RINGKASAN SKRIPSI S-1
ATENUASI R ESPON HUB UN GAN YA DEN GAN GEM PA B UM I SUBD UKSI SUMATER A
R ESPONSE ATTENUATION RELATIONSHIPS F OR SU MATRAN -SU BDUCTION EARTHQUAKE
Diajukan unt uk me menuhi s alah satu syarat mempe role h de rajat Sarjana Sains Ilmu Fis ika
WISNU SLAM ET PR IYAN TO 07/253936/PA/11750
PROGRAM STUD I GEOFISIKA JURUSAN FISIKA D EPAR TEM EN PENDIDIKAN N ASION AL UNIVERSITAS GADJAH M AD A FAKULTAS M ATEM ATIKA D AN ILM U PENGETAHUAN ALAM YOG YAKARTA 2012
i
RINGKASAN SKRIPSI S-1
ATENUASI R ESPON HUB UN GAN YA DEN GAN GEM PA B UM I SUBDUKSI SUM ATER A
RESPONSE SPECTR AL ATTENUATION RELATIONSHIPS F OR SU MA TRANSUBDUCTION EARTHQU AKE
Diajukan unt uk me menuhi s alah satu syarat mempe role h de rajat Sarjana Sains Ilmu Fis ika
WISNU SLAM ET PR IYAN TO 07/253936/PA/11750
PROGRAM STUD I GEOFISIKA JURUSAN FISIKA D EPAR TEM EN PENDIDIKAN N ASION AL UNIVERSITAS GADJAH M AD A FAKULTAS M ATEM ATIKA D AN ILM U PENGETAHUAN ALAM YOG YAKARTA 2012
ii
KATA PENGANTAR
Assalamualaikum Warahmatullahi Wabarakatuh, Puji syukur penulis panjatkan kepada Tuhan Yang Maha Esa yang telah memberikan rahmat dan hidayah Nya, sehingga penulis dapat menyelesaikan Laporan Tugas Akhir “Atenuasi Respon Hubunganya Dengan Gempa Bumi Subduksi” Tugas Akhir ini merupakan syarat yang wajib ditempuh di Prodi Geofisika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Gadjah Mada. Dengan selesainya laporan Tugas Akhir ini tidak terlepas dari bantuan banyak pihak yang telah memberikan masukan- masukan kepada penyusun. Untuk itu penyusun mengucapkan banyak terima kasih kepada : 1. Tuhan Yang Maha Esa yang telah memberikan rahmat dan karunianya. 2. Kepada kedua Orang tua saya, keempat adik saya dan keluarga besar yang selalu memberikan support
yang lain sehingga penyusun dapat
menyelesaikan Tugas Akhir ini dengan baik. 3. Bapak Drs. rer. nat. Wiwit Suryanto, Msi selaku dosen pembimbing yang dengan sabar membantu dan memberi diskusi, serta ilmu sehingga Tugas Akhir ini dapat terselesaikan dengan baik. 4. Bapak Drs.h. Suparwoto, M.sc. selaku dosen penguji yang telah menguji dan membantu dan memberikan diskusi mengenai hasil dari penelitian ini 5. Bapak Sudarmadji M.si selaku dosen penguji yang telah memberikan arahan dan masukan dalam pengujian penelitian ini. 6. Bapak Drs. Ari setiawan Msc. Selaku ketua jurusan Fisika. 7. Bapak Drs. Sismanto, selaku ketua Program Studi Geofisika UGM. 8. Kepada teman-teman geofisika’07 yang sudah melewati masa-masa bersama selama hampir 4 tahun. 9. Kepada teman-teman satu kontrakan OPC yang sudah lima tahun mengontrak bersama, banyak hal yang dilakukan bersama dan semua pihak yang terlibat. Penulis menyadari bahwa masih banyak kekurangan dari laporan ini, baik dari penelitianmateri maupun teknik penyajianya, mengingat kurangnya pengetahuan
iii
dan pemahaman penulis. Oleh karena itu, kritik dan saran yang membangun sangat penulis harapkan. Wassalamualaikum Warahmatullahi Wabarakatuh Yogyakarta, Juni 2012
Penulis
iv
Daftar Isi Halaman Halaman Sampul ........................................................................................................i Halaman Judul.......................................................................................................... ii Kata Pengantar ........................................................................................................ iii Daftar Isi ...................................................................................................................v Daftar Gambar....................................................................................................... viii Intisari ......................................................................................................................xi Abstract .................................................................................................................. xii BAB I PENDAHULUAN 1.1 LatarBelakang ............................................................................................... 1 1.2 Batasan Masalah ...........................................................................................2 1.3 TujuanPenelitian ...........................................................................................3 1.4 ManfaatPenelitian .........................................................................................3 1.5 Lokasi dan Waktu Penelitian .......................................................................3 BAB II TINJAUAN PUSTAKA 2.1 KeteranganUmum Daerah Penelitian............................................................4 2.2 Fisiografiprovinsi Sumatera Barat ................................................................4 2.3 TektonikLempeng .........................................................................................5 2.3.1 PetaTektonikGeologi Regional Sumatera ...........................................6 2.4 StudiKasus ....................................................................................................7 2.4.1 Suheri (2009).......................................................................................7 2.4.2 Saputro (2007).....................................................................................8 2.4.3 Megawati (2004) ..............................................................................10 2.5 Peta Daerah Penelitian ................................................................................12 BAB III DASAR TEORI 3.1 Seismik Hazard ...........................................................................................13 3.1.1 Analisis Hazard Kegempaan .............................................................13 3.2 FungsiAtenuasidan Yang Mempengaruhinya .............................................14 3.2.1 MekanismeGempa ............................................................................15 3.2.2 JarakEpisenter ...................................................................................15 3.2.3 Kondisi Tanah Lokal.........................................................................16 3.3 Peak Ground Velocity .................................................................................16 3.4 Perhitungan Jarak Episenter........................................................................18 3.5 FungsiAtenuasi............................................................................................19 BAB IV METODOLOGI PENELITIAN 4.1 Data Penelitian..............................................................................................21 4.2FasilitasPendukung .........................................................................................21 4.3Diagram AlirPengolahan ................................................................................22 4.4Editing Data Gelombang ................................................................................23 4.5Filtering Data Gelombang ..............................................................................24
4.6Pemilihan Event Data Gelombang .................................................................25 4.7Konversi Data Gelombang .............................................................................26 4.8Picking Maksimum Amplitude danMenentukan PGV ...................................27 BAB V HASIL DAN PEMBAHASAN 5.1 Picking Event.................................................................................................30 5.2 Grafik Peak Ground Velocity VsJarakEpisentrum........................................31 BAB VI PENUTUP 6.1 Kesimpulan ....................................................................................................42 6.2 Saran ..............................................................................................................42 DAFTAR PUSTAKA LAMPIRAN
DAFTAR GAMBAR hal Gambar 2.1 Gambar 2.2 Gambar 2.3 Gambar 2.4 Gambar 2.5 Gambar 2.6 Gambar 2.7
Gambar 2.8 Gambar 3.6 Gambar 4.1 Gambar 4.2 Gambar 4.3
Gambar 4.4a
Gambar 4.4b
Gambar 4.5 Gambar 4.6
Gambar 4.7 Gambar 4.8 Gambar 4.9 Gambar 4.10 Gambar 4.11a
(Bakosurtanal, 2005)…. Peta Fisiografi Kota Padang, Sumatera Bara t(I Made Sandy, 1985) …………..…………………………………………….. Distribusi Lempeng Tektonik Indonesia…………………… Peta Tektonik Geologi Regional Sumatera………………....
5 6 7
Pet a ting kat ba ha ya gemp ab u mi (d ala m g al) un tu k period e 100 t ah un (Suheri, 2009) ………………………………………
8
Pet a Kontu r Peak Gro un d A cc ele rat ion lin esou rce & Po int So u rce (Sap ut ro, 2007) ………………………………………….
9
Pet a kota Pad an g , Sumat era barat
4
Peak ground velocity simulasi horizontal untuk Gempa Bumi subduksi Sumatera dengan range moment magnitude 5 s/d 6 untuk kedalaman antara 15-33 km (Megawati et. al,
2004).………....................................................................... Peta Daerah Penelitian …………..……….............................. Notasi teori Koordinat bola…………………………......... Diagram Alir Penelitian………..………………………....... Tampilan Open File pada Geopsy Data Pada Tanggal 6 April 2011………………………………………………...... Contoh tampilan ed iting ge lombang untuk komponen vertikal dari mas ing- mas ing stas iun kedalam gra fik pada data tangga l 6 April 2011………………….….. Contoh data sebelum F ilter ing untuk data pada data tangga l 6 April 2011 untuk lima stasiun penga matan dengan kompone n vertika l.………………………….... Contoh data setelah Filtering untuk data pada data tanggal 6 April 2011 untuk lima stasiun pengamatan dengan komponen vertikal.…………………………………….......
11 12 18 22
Menu untuk pemotongan sinyal.…………………………… Kotak dengan garis merah digunakan untuk memperbesar dan memperjelas sinyal untuk contoh untuk data pada data tanggal 6 April 2011………………………………….......... Contoh Sinyal yang telah diperjelas untuk data pada data tanggal 6 April 2011………………………………………. Type File Eksport…………..……………………………... Type File Eksport …………….…………………………... Menu Open File Pada Seisgram2k………………………... Contoh Picking gelombang pada lima stasiun pengamatan mencari nilai amplitude maksimum pada data tanggal 6
25
23
24
24
25
26 26 27 27 28
Gambar 4.11b
Gambar 5.1 Gambar 5.2a
Gambar 5.2b
Gambar 5.2c
Gambar 5.2d
Gambar 5.2e
Gambar 5.3a
Gambar 5.3b
Gambar 5.3c
Gambar 5.3d
Gambar 5.3e
April 2011.…………………………………......................
28
Contoh Gambar diperbesar untuk P-P max merupakan amplitude maksimum Stasiun Pengamatan 2 pada data tanggal 6 April 2011.……………..…………………….....
29
Gempa pada tanggal 25 Februari 2011 pada pukul 22.09 tercatat di lima stasiun.……………………………………. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 2 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 4 s/d 5……………………................................................. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 4 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 4 s/d 5…………………………………………………….….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 5 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 4 s/d 5………………………………………………………... Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 6 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 4 s/d 5…………………………………………………….….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 7 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 4 s/d 5………………………………………………………... Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 2 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 5 s/d 6………………………………………………………... Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 4 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 5 s/d 6………………………………………………………... Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 5 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 5 s/d 6……………………………………………………….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 6 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 5 s/d 6……………………………………………………….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 7 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 5 s/d
31
33
33
34
34
35
35
36
36
37
Gambar 5.4a
Gambar 5.4b
Gambar 5.4c
Gambar 5.4d
Gambar 5.4e
6……………………………………………………….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 2 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 6 s/d 7……………………………………………………….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 4 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 6 s/d 7…………………………………………………….….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 5 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 6 s/d 7……………………………………………………….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 6 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 6 s/d 7……………………………………………………….. Grafik Peak ground veloc ity s imulas i hor izonta l Stasiun penga matan 7 untuk ge mpa subduksi Sumatera barat dengan range moment magnitude 6 s/d 7………………………………………………………...
37
38
38
39
39
40
INTISARI
Atenuasi Respon Hubunganya Dengan Gempa Bumi Subduksi Sumatera. Oleh Wisnu Slamet Priyanto 07/253936/PA/11750
Kota Padang merupakan daerah gempa yang memiliki nilai seismisitas tinggi di Asia Tenggara. Pulau Sumatera yang terletak dikawasan yang terdapat lempeng tektonik aktif, subduksi lempeng di Sumatera berada pada zona Sunda arc. Penelitian dilakukan untuk mengetahui respon atenuasi dari pergerakan tanah di Kota Padang. Dalam penelitian ini menggunakan metode Regresi dengan acuan persamaan Molas dan Yamazaki. Pada penelitian kali ini dengan mendapatkan nilai baru dari hubungan atenuasi dari Peak Ground Velocity (PGV) pada gempa bumi subduksi di Sumatera. Hasil dari penelitian ini dapat memperkirakan seismik hazard Kota Padang atau Indonesia secara umum dengan daerah jangkauan dengan jarak dari 20-350 km dan moment magnitude dengan jangkauan 4-7. Berdasarkan perhitungan hasil regresi respon atenuasi pergerakan kecepatan tanah maksimum, dari perhitungan didapatkan nilai koefisien regresi. Nilai koefisien tersebut menggambarkan kemungkinan tingkat kerusakan yang terjadi akibat kejadian subduksi besar dengan nilai magnitudo tinggi yang mungkin dapat terjadi di kota Padang, Sumatera Barat untuk kejadian dengan jangkauan jarak 350 km.
xi
ABSTRACT
Response Attenuation Relationships for Sumatran-Subduction Earthquakes
Wisnu Slamet Priyanto 07/253936/PA/11750
Padang City representing an earthquake region that located in highseismicity of Southeast Asia. The island of sumatra contained active plate, the Sumatran subduction zone lies within the sunda arc. Investigation have been done by knowing response spectral attenuation from ground motions in Padang City. In this research using regression methods that relations with Molas and Yamazaki Equations. The present research particularly addressed issues in deriving new set of attenuation relationships of Peak Ground Velocity (PGV) for the Sumatran subduction earthquakes. To be relevant for seismic hazard assessment of the remote metropolises, the attenuation relationships cover a long distance range 20350 km and magnitude moment range 4.0 – 7.0. Based on calculation result of regretion response attenuation Peak Ground Velocity, the value obtained from the calculation of the regression coefficients. That value was described a larger event with great Mw in this particular subduction zone may capable of generating destructive ground motions in Padang City, even at distance 350 km.
xii
BAB I PENDAHULUAN 1.1. Latar Belakang Indonesia merupakan negara kepualauan yang secara tektonik berada pada tiga lempeng tektonik besar, yaitu lempeng Eurasia, lempeng Indo-Australia dan lempeng Pasifik. Hal tersebut menyebabkan Indonesia menjadi negara yang berpotensi terjadi bencana kebumian. Selain itu kondisi topografi dan geologi permukaan sangatlah mempunyai pengaruh yang tinggi terhadap tingkat kerusakan yang diakibatkan oleh gempa bumi (Daryono, 2011). Beberapa bencana yang sangat berpotensi terjadi dan telah terjadi di Sumatera Barat diantaranya adalah gempa bumi. Peraturan Gempa Indonesia (SNI‐1726, 2002) menempatkan Provinsi Sumatera Barat sebagai salah satu provinsi yang memiliki percepatan gempa maksimum (PGA) tertinggi di Indonesia. Hal ini menunjukkan bahwa Provinsi Sumatera Barat bisa dipastikan adalah daerah yang rawan terhadap bencana gempa bumi. Wilayah Sumatera Barat yang paling rentan terhadap bencana ini adalah Kepulauan Mentawai yang berada pada zona subduksi lempeng Indo-Australia. Setelah
kejadian-kejadian
gempabumi
bermunculan
maka
akan
memunculkan pemikiran mengenai bagaimana mengetahui hubungan distribusi sumber gempa suatu wilayah serta seberapa besar kemungkinan kerusakan yang terjadi pada suatu daerah dengan radius pada jarak tertentu. Perhitungan-perhitungan mengenai tingkat bahaya gempa bumi telah banyak dilakukan untuk mengetahui seberapa besar tingkat bahaya yang diakibatkan apabila terjadi gempabumi. Perhitungan tersebut melibatkan parameter-parameter yang menyangkut getaran tanah maupun intensitas gempa bumi. Perhitungan dapat dilakukan dengan metode deterministik maupun probabilistik. Untuk menghubungkan antara data gempa dengan model yang akan digunakan dalam analisis bahaya gempa diperlukan identifikasi dan karakterisasi semua potensi sumber gempa yang mungkin memberikan pengaruh signifikan pada ground motion di lokasi gempa. 1
2
Sumber gempa dapat diidentifikasi dari geologi, tektonik, sejarah dan dari hasil pencatatan gempa. Analisis penilaian bencana gempa melibatkan estimasi kuantitatif dari karakteristik ground motion pada site lokasi. Pendekatan analisis yang digunakan dalam penentuan analisis bencana gempa tersebut terdiri atas dua pendekatan yaitu secara deterministik (Deterministic Seismic Hazard Analysis (DSHA)) dan secara probabilistik (Probabilistic Seismic Hazard Analysis (PSHA)). Metoda deterministik merupakan metoda dengan menggunakan input data skenario gempabumi dari satu sumber energi patahan gempa bumi yang paling berpotensi untuk menimbulkan bencana di wilayah yang bersangkutan. Perhitungan dengan metode probabilistik mulai dikenal pertama kali oleh Cornell (1968). Berkembang kemudian untuk menghitung bahaya gempabumi pada daerah yang memiliki bangunan yang vital, Mc Guire (1976) membuat program EQRISK untuk menghitung bahaya berdasarakan paparan metode Cornell (1968). Asumsi dasar yang digunakan dalam perhitungan tingkat bahaya gempabumi adalah bahwa aktifitas gempabumi akan terjadi kembali disuatu area pada masa yang akan datang. Kejadian ini akan terjadi secara tektonik berpotensi menghasilkan gempabumi dan dianggap sebagai suatu yang bersifat probabilistik. Untuk itu perhitungan tingkat bahaya gempabumi ini secara probabilistik seismik hazard. Dalam PSHA tingkat dari bahaya gempabumi H(V) didefinisikan sebagai gempabumi tahunan yang menghasilkan amplitudo getaran tanah a melampaui nilai harapan V pada suatu area tertentu. Variabel bebas V umumnya menggambarkan kuantitas bahaya yang terkait dengan getaran tanah,kali ini PGV
1.2. Batasan Masalah Penelitian ini memiliki beberapa batasan masalah antara lain : 1. Data yang digunakan adalah perekaman gelombang di 5 stasiun pengamatan di daerah Padang Sumatera barat selama 3 bulan dan data katalog USGS dari bulan Febuari-April 2011. 2. Data yang dimasukan dalam pengolahan ini adalah amplitudo, magnitude dan jarak episentrum.
3
3. Metode yang digunakan oleh analisa dalam menentukan peak ground velocity adalah metode regresi.
1.3. Tujuan Penelitian Tujuan dari penelitian ini adalah : 1. Menentukan Peak Ground Velocity (PGV) di Kota Padang. 2. Mendapatkan fungsi atenuasi baru untuk gempabumi subduksi di wilayah Padang, Sumatera Barat.
1.4. Manfaat Penelitian Penelitian ini diharapkan dapat memberikan manfaat, antara lain : 1. Penelitian ini dijadikan sebagai salah satu syarat untuk memperoleh derajat Sarjana S-1. 2. Penelitian ini memberikan kesimpulan yang dapat digunakan untuk study Seismik Hazard Analisis di wilayah Padang atau Indonesia secara umum .
1.5. Lokasi dan Waktu Penelitian Penelitian dilakukan di Laboratorium Geofisika Universitas Gadjah Mada. Waktu penelitian dimulai dari bulan 1 September sampai 30 November 2011.
BAB II TINJAUAN PUSTAKA 2.1. Keterangan Umum Daerah Penelitian Provinsi Sumatera Barat (Ga mbar 2.1 ) terletak pada koordinat 3º 50' LS - 1º 20' LU dan 98º 10' - 102º 10' BT dengan memiliki total luas wilayah sekitar 42.297,30 km2 .
Gambar 2.1. Peta kota Padang, Sumatera barat ( Sumber: Hasil Modifikasi dari “Peta RBI Sumatera Barat” (Bakosurtanal, 2005))
Wilayah Sumatera Barat sebelah bawah terdapat lempeng samudera yang bergerak
relatif tegak lurus terhadap bentuk Pulau Sumatera dengan
kecepatan 6 – 7 cm/tahun khususnya menjadi faktor sering terjadinya gempa bumi setelah periode tertentu yang berada pada bidang batas kontak lempeng samudera bertemu dengan lempeng kerak kepulauan. 2.2. Fisiografi Provinsi Sumatera Barat Sumatera Barat merupakan provinsi yang memiliki kondisi fisiografi sangat 4
5
kompleks. Menurut Sandy (1985) di Sumatera Barat dapat ditemui tiga wilayah fisiografi utama, yaitu: wilayah pegunungan vulkanik, perbukitan lipatan tersier, dan wilayah dataran rendah (Gambar 2.2).
Gambar 2.2. Peta Fisiografi kota Padang, Sumatera barat (Sumber peta: hasil modifikasi dari “Peta Wilayah Fisiografi” (I Made Sandy, 1985))
Wilayah pegunungan vulkanik membujur pada bagian tengah provinsi ini dari utara sampai selatan, dengan patahan semangko pada tengahnya. Sedangkan perbukitan lipatan tersier membentang dibagian timur pegunungan vulkanik tersebut. Perbukitan tersier ini di beberapa tempat mengandung deposit batubara dengan medan yang curam. Sementara itu pada bagian barat provinsi ini terdapat dataran rendah. 2.3. Tektonik Lempeng Teori tektonika Lempeng (Plate Tectonics) adalah teori dalam bidang geologi yang dikembangkan untuk memberi penjelasan terhadap adanya buktibukti pergerakan skala besar yang dilakukan oleh litosfer bumi. Teori ini telah mencakup dan juga menggantikan Teori Pergeseran Benua yang lebih dahulu dikemukakan pada paruh pertama abad ke-20 dan konsep seafloor spreading yang dikembangkan pada tahun 1960-an (Gambar 2.3).
6
Gambar 2.3 Distribusi Lempeng Te ktonik Indonesia
2.3.1 Peta Tektonik Geologi Regional Sumatera Pulau Sumatra terletak dibaratdaya dari Kontinen Sundaland dan merupakan jalur konvergensi antara Lempeng Hindia-Australia yang menyusup di sebelah barat Lempeng Eurasia/Sundaland. Konvergensi lempeng menghasilkan subduksi sepanjang Palung Sunda dan pergerakan lateral ke arah kanan dari Sistem Sesar Sumatra.
7
Gambar 2.4 Peta tektonik geologi regional Su maeta (Sumber: Wikimapia “Map Sumatera”)
2.4 Studi Kasus 2.4.1 Suheri (2009) Penelitian juga telah dilakukan mengenai pergerakan tanah oleh Suheri (2009) dan Saputro (2007). Penelitian yang dilakukan oleh Suheri (2009) yaitu menghitung tingkat bahaya gempabumi di Jawa Tengah secara probabilistik. Suheri
(2009)
meneliti
mengenai
tingkat
bahaya
gempabumi
dengan
menggunakan percepatan pergerakan tanah maksimum (PGA) yang diolah dengan menggunakan metode Horizontal to Vertical Spectral Ratio (HVSR). Metode ini pertama kali diperkenalkan oleh Nogoshi dan Iragashi (1971) yang menyatakan adanya hubungan antara perbandingan komponen horisontal dan vertikal terhadap kurva eliptisitas pada gelombang Rayleigh. Metode ini kemudian disempurnakan
8
oleh Nakamura (1989) yang menyatakan perbandingan spektrum H/V sebagai fungsi frekuensi berhubungan erat dengan fungsi waktu perjalanan untuk gelombang S ke stasiun. Model seismisitas zona sumber gempabumi beserta parameter seismisitasnya pada penelitian milik Suheri (2009) dibuat dan dihitung berdasarkan data seismisitas dan geologinya. Kemudian nilai percepatan maksimumnya dihitung menggunakan metode Tong dan Katayama (1988), Boore et. al. (1994) dan Mcguire (1976). Hasil dari penelitian yang dilakukan oleh Suheri (2009) berhasil memperoleh nilai periode dominan tanah di Jawa Tengah, serta peta tingkat bahaya gempabumi (Gambar 1.1) yang dibuat.
Gambar 2.5 Peta tingkat bahaya gempabumi (dalam gal) untuk periode 100 tahun (Suheri, 2009)
2.4.2 Saputro (2007) Penelitian dilakukan oleh Saputro (2007) dengan memetakan percepatan maksimum pergerakan tanah atau dikenal juga peak ground acceleration (Gambar 1.2) akibat dari gempa utama Yogyakarta 27 Mei 2006 dan susulanya dengan menggunakan metode Kanai (1966) yaitu pendekatan garis untuk sumber garis untuk gempa utama dan sumber titik untuk gempa susulan. Metode penurunan garis (Line Source Derivation Method) yaitu dengan menganggap sumber gempa
9
sebagai garis karena merupakan suatu sesar penyebab gempa yang diteliti dan perhitunganya untuk site/stasiun adalah tegak lurus dengan arah strike atau fault tersebut. Metode penurunan titik (Point Source Derivation Method) merupakan metode yang menggambarkan sumber gempa/source pada perhitungan percepatan getaran tanah dianggap sebagai point yang tersebar disuatu region penelitian sehingga nantinya diambil sebagai harga maksimum dari masing- masing site. Perioda dominan pada batuan sedimen dihitung menggunakan frekuensi resonansi dari perambatan gelombang sekunder mulai dari paling atas dari lapisan dasar hingga permukaan tanah. Hasil dari penelitian yang dilakukan oleh Saputro (2007) menunjukan bahwa terdapat 4 range persebaran nilai percepatan getaran tanah maksimum. Penelitian yang dilakukan oleh Kanai (1966), yaitu menentukan percepatan getaran tanah maksimum dengan menyertakan tingkat bahaya gempa bumi dari geologi maupun lithologi terutama adanya faktor amplifikasi dan diamplifikasi suatu wilayah (Kanai, 1966).
Gambar 2.6 Peta Kontur Peak Ground Acceleration linesource & Point Source (Saputro, 2007)
Penelitian kali ini dilakukan untuk mengetahui respon atenuasi dari pergerakan kecepatan tanah maksimum (Peak Ground Velocity) di Kota Padang, serta memperkirakan seismik hazard yang relevan secara umum daerah jangkauan dengan jarak dan moment magnitude tertentu.
10
Peak Ground Velocity merupakan amplitudo puncak dari kecepatan pergerakan tanah. PGV memberikan korelasi terbaik mengenai akibat kerusakan dan PGV hanya menyatakan puncak pertama integrasi dari catatan percepatan, sehingga kita dapat mengetahui dampak paling besar akibat dari sebuah event gempa. Dalam studi gempa bumi Peak Ground Velocity masih jarang dipakai, kebanyakan peneliti meneliti mengenai Peak Ground Acceleration (PGA) yang menitik beratkan pada respon percepatan yang terjadi apabila terjadi suatu event. Dalam penelitian ini penulis ingin meneliti pertama kali mengenai respon kecepatan maksimum pergerakan tanah (Peak Ground Velocity) untuk daerah Padang, Sumatera Barat. 2.4.3 Megawati (2004) Megawati et. al. (2004) melakukan penelitian mengenai atenuasi respon spektral untuk subduksi gempa pulau Sumatra dan pengaruhnya terhadap Singapura dan kota Kuala Lumpur. Penelitian dilakukan pada daerah yang memiliki nilai seismisitas yang rendah dan jarak yang jauh yaitu dari stasiun pengamatan. Dalam penelitianya Megawati (2004) menggunakan acuan dari persamaan Molas dan Yamazaki (1995) serta persamaan Atkinsoon dan Boore (1997) yang kemudian diplotkan dan dicari hasil regeresinya kedalam grafik Peak Ground Velocity. Hasil dari regeresi grafik tersebut didapatkan persamaan atenuasi baru untuk Singapura dan Kuala lumpur. Gambar 2.5 merupakan Peak ground velocity simulasi horizontal untuk subduksi sumatera yang diambil dari penelitian yang dilakukan oleh Megawati et. al. (2004) data diatas di plot berdasarkan range moment magnitude 5 s/d 6 untuk kedalaman antara 15-33 km untuk data yang diambil dari Singapura dan Kuala Lumpur dari sumber subduksi gempa di Sumatera pada periode Bulan Juni tahun 2000. Gambar 2.5 menjelaskan mengenai metode regeresi yang diambil dengan cara mendekatkan antara grafik data nilai peak ground velocity oleh Megawati et. al. (2004) dengan persamaan Atkinson dan Boore (1997) serta
Molas dan
Yamazaki (1995). Dari gambar tersebut maka didapatkan hasil nilai dari koefisien atenuasi untuk peak ground velocity untuk Singapura dan kota Kuala Lumpur.
11
Hasil persamaan dari rumus Megawati et. al. (2004) didapatkan dari hasil modifikasi atenuasi Atkinson dan Boore (1997) serta Molas dan Yamazaki (1995) yang kemudian digabungkan dan diberi koefisien baru.
Gambar 2.7 Peak ground velocity simu lasi horizontal untuk Gempa Bu mi subduksi Sumatera dengan range mo ment magnitude 5 s/d 6 untuk kedalaman antara 15-33 km. (Megawati et. al., 2004)
Dari koefisien persamaan tersebut kemudian Megawati et. al. (2004) dapat menjelaskan mengenai seismik hazard yang menggambarkan kemungkinan tingkat kerusakan yang terjadi akibat kejadian subduksi dengan nilai magnitudo tinggi yang mungkin dapat terjadi di Singapura dan Kuala Lumpur, Malaysia untuk observasi dengan jangkauan jarak 200-1000 km.
12
2.5 Peta Daerah Penelitian Data gelombang direkam pada 5 stasiun pengamatan di kota Padang, Sumatera barat , yaitu stasiun pengamatan 2 ( Sikabu ), stasiun pengamatan 4 (Air Dingin ), stasiun pengamatan 5 ( Sungai Lareh ), stasiun pengamatan 6 ( Bukit Gado-Gado ), stasiun pengamatan 7 ( Tarantang ). Hal tersebut sebagaimana digambarkan pada gambar peta daerah penelitian (Gambar 2.6).
Gambar 2.8 Peta Vs 3 daerah penelit ian dengan koordinat stasiun
BAB III DASAR TEORI 3.1 Seismik Hazard Tujuan seismic hazard analysis (SHA) adalah untuk mengkuantifikasi peluang terjadinya dan terlampauinya berbagai tingkat percepatan maupun kecepatan maksimum tanah akibat gempa bumi yang dapat terjadi di suatu daerah. Hasil dari studi seismik hazard ini berupa percepatan maksimum pergerakan tanah, kecepatan maksimum pergerakan tanah dan amplitudo dengan respon spektral. Nilai peak ground velocity (PGV) dan besaran karakteristik lainnya digunakan sebagai parameter perhitungan pergerakan tanah dalam analisis penelitian kali ini. Sumatera Barat merupakan daerah yang memiliki tingkat kerentanan sangat tinggi terhadap kejadian bencana alam. Tingkat kerentanan yang tinggi tersebut terjadi karena banyak hal, seperti lokasi yang berada di kawasan zona subduksi lempeng indo-australia, keberadaan patahan besar semangko yang sampai saat ini masih aktif, keberadaan beberapa gunung berapi aktif, kondisi pantai yang landai, banyaknya sungai besar yang mengalir di wilayah ini, dan masih banyak lagi alasan lain yang menyebabkan propinsi ini sangat rentan terhadap kejadian bencana alam. 3.1.1. Analisis hazard kegempaan Analisis hazard kegempaan dimulai dengan mengembangkan model matematik yang akan digunakan untuk memperkirakan kemungkinan kejadian gempa dalam level skala magnitude atau intensitas tertentu pada interval periode ulang untuk suatu daerah tertentu. Analisis ini menghasilkan parameter desain seismik seperti percepatan maksimum dan kecepatan maksimum yang dapat terlampaui untuk probabilitas serta periode ulang tertentu. Pada studi ini, kecepatan maksimum gempa diperoleh dari hasil analisis yang dilakukan dengan metode regeresi menggunakan software Geopsy dan Seisgram2k. Nilai amplitude maksimum didapat dari picking data gelombang yang merupakan respon dari
13
14
kecepatan maksimum dari sebuah event yang terjadi dan terekam pada stasiun pengamatan. Analisis untuk kemungkinan magnitude dan jarak dari site ke sumber gempa yang akan memberikan hazard terbesar pada site tidak terlihat dengan jelas dalam seismic
hazard
analisis.
Dengan
menggunakan
metode regeresi
menggunakan komponen magnitude (M), jarak episenter (D) dan kecepatan maksimum pergerakan tanah (PGV) dapat digambarkan da lam suatu fungsi persamaan. Konsep ini dapat memberikan gambaran mengenai ukuran pergerakan tanah maksimum pada gempa dengan magnitude dan jarak untuk sumber gempa tertentu, yang kemungkinan besar besar destruktif terhadap suatu daerah. Metode regresi dapat mengkarakterisasi dengan suatu rata-rata magnitude dan jarak dengan nilai amplitude maksimum. Dengan melakukan plotting seperti ini, maka perbedaan hazard dalam suatu bentang magnitude untuk suatu periode spektral akan dapat terlihat. Gempa bumi yang dominan dalam memberikan kerusakan yang berarti terhadap suatu wilayah. 3.2 Fungsi atenuasi dan Faktor yang Mempengaruhinya Beberapa fungsi atenuasi telah dipublikasikan oleh sejumlah peneliti berdasarkan rekaman percepatan gempa yang pernah terjadi dan kondisi site kajiannya. Prediksi hubungan empiris untuk parameter gempa yang melemah (berkurang) sejalan dengan bertambahnya jarak, seperti percepatan puncak (Peak Ground acceleration) dan kecepatan puncak (Peak Ground Velocity), dikenal sebagai fungsi atenuasi (attenuation relationship atau attenuation function). Fungsi atenuasi adalah suatu fungsi yang menggambarkan korelasi antara intensitas (i) gerakan tanah setempat, magnitude (M) dan jarak (R) dari sumber titik dalam daerah sumber gempa. Memperkirakan fungsi atenuasi untuk gerakan tanah akibat gempa, telah menjadi bahasan yang menarik dalam penelitian bidang kegempaan. Fungsi atenuasi merupakan alat yang penting
dalam
mengaplikasikan
resiko
kegempaan
bangunan tahan gempa. Faktor-faktor yang mempengaruhi fungsi atenuasi adalah :
dalam perencanaan
15
3.2.1 Mekanisme gempa Gempa- gempa besar biasanya
terjadi karena
pergeseran tiba-tiba
lempeng tektonik yang mengakibatkan terlepasnya energi yang sangat besar. Pergeseran tektonik ini biasa terjadi pada daerah subduction, ataupun pada patahan yang tampak di permukaan bumi, seperti patahan semangko di sumatera. Gempa yang terjadi pada daerah subduction biasanya merupakan gempa dalam yang mempunyai kandungan frekuensi yang berbeda dengan gempa dangkal. Gempa dalam biasanya mempunyai gelombang permukaan yang lebih sedikit, sehingga memberikan spektrum respon yang lebih rendah pada periode tinggi. 3.2.2 Jarak episenter Respon spektrum dari gempa yang tercatat pada batuan mempunyai bentuk yang berbeda tergantung jarak episenternya (near field, mid field, dan far field). Gempa near field memberikan respon yang tinggi pada perioda yang rendah tapi mengecil secara drastis dengan bertambah perioda. Di lain hal, gempa far field pada perioda rendah tetapi responnya terlihat konstan sampai perioda sekitar satu detik.
Hal
ini
menunjukkan
adanya
perubahan
kandungan frekuensi gempa dengan semakin jauhnya suatu daerah yang ditinjau ke episenter. Ada beberapa sistem untuk menyatakan besarnya sudut, diantaranya yaitu : a. Sistem Seksagesimal. Dalam sistem seksagesimal keliling lingkaran dibagi dalam 360 bagian yang disebut derajat. 1o (1 derajad) = 60’ (60 menit) dan 1’ = 60” (60 detik). b. Sistem Sentisimal Dalam sistem sentisimal keliling lingkaran dibagi dalam 400 bagian yang disebut grade. 1g (1 grade) = 100c (100 centigrade) dan 1c = 100cc (100 centicentigrade). Sistem Radial Dalam sistem radial keliling lingkaran dibagi dalam bagian yang disebut dengan satu radial. c . Sistem Waktu,
16
Sistem waktu digunakan dalam pengukuran astronomi. Dimana, 360o = 24 jam; 1 jam =15o 3.2.3 Kondisi tanah lokal Kondisi tanah lokal mempunyai peran yang sangat penting dalam menentukan respon
suatu
daerah
terhadap
gelombang
gempa.
Respon
gempa yang tiba dibatuan dasar bisa diperkuat, diperlemah atau berubah kandungan frekuensinya karena tersaringnya getaran berfrekuensi tinggi. 3.3 Peak Ground Velocity Kecepatan getaran tanah maksimum atau peak ground velocity (PGV) pada saat ini masih mendapat perhatian yang jauh lebih sedikit dalam literatur teknis daripada parameter digunakan yang lebih luas, seperti percepatan getaran tanah maksimum atau peak ground acceleration (PGA) dan koordinat respon spektral. Namun, ada banyak contoh dari penggunaan PGV dalam rekayasa seismologi dan rekayasa
gempa,
termasuk
sebagai
parameter
yang
digunakan
untuk
memperkirakan intensitas makroseismik dan kerusakan struktural. Dalam hal ini peak ground velocity sangat baik dalam mengkorelasi akibat kerusakan dan PGV hanya menyatakan puncak pertama integrasi dari catatan percepatan. Dalam pengolahan data kali ini menggunakan persamaan yang beberapa persamaan hubungan atenuasi dengan gempa bumi subduksi pada beberapa wilayah. Atkinsoon dan Boore (1997) menerangkan hubungan preliminary dengan persamaan atenuasi untuk daerah Kaskadia, dengan hanya mencakup kejadian dengan magnitude kurang dari 7.
Molas dan Yamazaki
(1995) menerangkan mengenai kombinasi dari pergerakan tanah dari Shallow Crustal dan zona subduksi gempabumi menjadi satu paket yang mengacu pada hubungan atenuasi di Jepang. Megawati et. al. (2004) melakukan penelitian di daerah yang memiliki nilai seismisitas yang rendah dan dilakukan di daerah gempa bumi subduksi Sumatera untuk jarak yang jauh yaitu dari stasiun pengamatan di Singapura dan Kuala Lumpur, Malaysia serta menggunakan acuan dari persamaan Molas dan Yamazaki (1995) serta persamaan Atkinsoon dan Boore (1997).
17
Tabel 2.1 Referensi persamaan ground motion Referensi
Notasi
Data yang digunakan
Mw
R (km)
model Atkinson and
AB97
Boore (24)
Nilai prediksi
Simulasi stokhastik
4.0-
100-
pergerakan tanah akibat
8.25 400
gempabumi subduksi di
PGA, RSA, PGV
Kaskadia Molas
and
Rekaman shallow crustal
4.0-
Yamazaki
dan zona gempabumi
7.0
(23)
subduksi di jepang
Kusnowidjaja Megawati
MY95
MG04
20-600
PGA, PGV
Respon spektral atenuasi
5.0-
200-
PGA,
dan hubunganya dengan
6.5
1000
RSA,
subduksi di Sumatra dan
PGV
pengaruhnya terhadap Singapura dan Kuala Lumpur
Pada dasarnya persamaan Molas dan Yamazaki (1995), Atkinsoon dan Boore (1997) dan Megawati et. al (2004). dapat menjelaskan hubungan dari mengenai Peak Ground Velocity, yang memungkinkan untuk memprediksi hubungan atenuasi antara pergerakan tanah (Ground motion) dan jarak gempabumi subduksi di Sumatera Barat. Hal tersebut didasarkan bahwa kesimpulanya persamaan Molas dan Yamazaki (1995), Atkinsoon dan Boore (1997) dan Megawati et. al (2004) dapat memprediksi data rekaman gempabumi subduksi dengan cukup akurat dari respon amplitudo yang terekam pada stasiun pengamatan. Tetapi semua tergantung terhadap jarak episentrum yang digunakan.
18
3.4 Perhitungan Jarak episenter Dalam perhitungan jarak diantara dua koordinat bumi biasanya dihitung menggunakan satuan derajat tetapi dapat juga dalam kilometer ( 1 0 =111.1 km). bumi sebenarnya tidak benar-benar bulat. maka dari itu biasanya dalam menghitung jarak episentrum dan azimuth mengunakan spherical trigonometry, kita harus merubah latitude
menjadi geocentric latitude
, rumusnya : (3.7)
dimana
merupakan konstanta perataan tanah.
Gambar 3.6 Notasi teori Koordinat bola
(Handoko, 2004) Untuk menghitung jarak dan azimuth akan tepat bila menggunakan cosines garis geosentris dari stasiun dan episenter. Kita anggap
( positif untuk belahan
bumi bagian utara, dan negative nilainya untuk belahan bumi bagian selatan) menjadi latitude dan
menjadi longitude ( dihitung positif dari arah Greenwich ke
timur dan bervariasi nilainya antara 0 o hingga 360o ), berikut adalah cosines garis geosentris
(Gambar 3.6): (3.8)
Menghitung Jarak Episentrum dalam derajat(∆):
19
(3.9) Dimana
merupakan episentrum dan
adalah stasiun, sedangkan
menghitung jarak dalam kilometer yaitu : (3.9) 3.5 Fungsi Atenuasi Berkembangnya penelitian mengenai kegempaan juga memunculkan banyak literature sebagai pedoman dan biasa menjadi alat konseptual dalam melakukan penelitian, pekembangan mengenai penelitian pergerakan tanah telah berlangsung sejak abad 19 hingga kini, hal tersebut di jelaskan dalam jurnal milik Jhon Douglas. Dibawah ini beberapa persamaan atenuasi untuk perge rakan tanah yang akan digunakan dalam penelitian ini sebagai referensi dasar. a. Persamaan Rumus Atkinsoon and Boore 1997 yang digunakan: (3.10) dengan: = Kecepatan maksimum pergerakan tanah
= Moment magnitude = Jarak episenter = Koefisien regeresi Tabel 3.1 Koefisien Regeresi yang digunakan: Detik(S) PGV
4.903
1.223
0.000
0.00253
b. Persamaan Rumus Molas and Yamazaki 1995 yang digunakan: (3.11) dengan: = Kecepatan maksimum pergerakan tanah
= Moment magnitude
20
= Jarak episenter = Koefisien regeresi = Kedalaman Tabel 3.2 Koefisien Regeresi yang digunakan: Detik(S) PGV
0.206
0.477
-0.0144
-1
0.00311
c. Persamaan Rumus Megawati et. al 2004 yang digunakan: (3.12) dengan: = Kecepatan maksimum pergerakan tanah
= Moment magnitude = Jarak episenter = Koefisien regeresi = Kedalaman =
Koefisien batuan
Tabel 3.3 Koefisien Regeresi yang digunakan: Detik (S) PGV
-12.155
3.3842
-0.085810
-0.9953
-0.000597
-0.08551
0.4776
BAB IV METODOLOGI PENELITIAN 4.1. Data Penelitian Penelitian menggunakan data gelombang yang merupakan data perekaman gempa di daerah Padang, Sumatera Barat selama tiga bulan, yaitu selama bulan Febuari, Maret dan April tahun 2011. Daerah pengamatan ada 5 stasiun pengamatan dengan koordinat posisi seperti tabel 4.1 dibawah ini: Tabel 4.1. Data Koordinat Posisi Stasiun No.
Stasiun
LATITUDE LONGITUDE
1
Sikabu ( SP2 )
00 41’19”S
1000 19’00”E
2
Air Dingin ( SP4 )
00 49’41”S
1000 23’28”E
3
Sungai Lareh ( SP5 )
00 51’46”S
1000 24’04”E
4
Bukit Gadi- Gado ( SP6 ) 00 58’47”S
1000 22’02”E
5
Tarantang ( SP7 )
00 57’94”S
1000 25’60”E
Data perakaman gelombang gempa ini terdiri atas tiga macam sensor, yaitu utara, timur dan vertikal. Dalam penelitian ini digunakan sensor vertikal untuk menentukan respon kejadian gempa pada tiap-tiap stasiun pengamatan. untuk data gelombang berbentuk data perekaman selama 1 hari yang dibagi sesuai dengan waktu 1 hari normal (24 jam). Data perekaman gempa bertipe .seed, data melalui proses
pengeditanan agar memudahkan dalam proses picking.
Keuntungan data melalui proses pengeditan bentuk gelombang yaitu dapat mempermudah dalam mengenali sebuah event. Untuk pengeditan, filtering dan pemilihan event gempa menggunakan software geopsy setelah itu mencari nilai dari amplitude maksimumnya dengan menggunakan software Seisgram2k. 4.2. Fasilitas Pendukung Fasilitas pendukung untuk pengolahan data penelitian ini menggunakan software Geopsy yang digunakan untuk pengeditan, filtering, pemilihan event dan mengkonversi data gelombang, kemudian untuk pengepick-an Amplitudo
21
22
maksimum menggunakan software Seisgram2k, software lainya juga Microsoft office 2007 dan Surfer. 4.3. Diagram Alir Pengolahan Alur kerja dari penelitian ini dapat dijelaskan secara garis besar dengan diagram alir seperti gambar 4.1 di bawah ini: MULAI
Data Perekaman Gelombang Gempa
Pengeditan Data
Pemilihan Kejadian
Informasi Data Katalog
Picking Maksimum
Penentuan Peak Ground
Penentuan Distance
Grafik Peak Ground Velocity Vs Distance Episenter
SELESAI
Gambar 4.1 Diagram alir penelit ian
23
4.4. Editing data gelombang Pengolahan data dimulai dengan editing data perekaman gelombang, yaitu memilih data untuk komponen vertikalnya saja pada masing- masing stasiun dan menyatukanya ke dalam sebuah grafik, hal tersebut dilakukan dengan tujuan mempermudah pengolahan dalam proses pemilihan gelombang. Untuk proses editing ini digunakan software Geopsy. Data awal sudah berbentuk .seed maka dapat langsung menjadi input pada software geopsy. Untuk membuka file cukup dengan meng-klik file pada menu bar geopsy kemudian pilih import signals kemudian file.
Gambar 4.2 Tampilan Open File pada Geops y data pada tanggal 6 April 2011.
Setelah file kita buka misalnya untuk stasiun pengamatan 2 data pada tanggal 1 April masih menjadi 3 komponen, yaitu komponen North, East dan vertical (Gambar 4.2). Untuk pengolahan data kali ini kita memisahkan komponen tersebut dengan mengeplot pada grafik untuk
komponen vertikal pada setiap
stasiun. Untuk pengeplotan dapat dengan cara meng-klik view pada menu bar geopsy kemudian klik grafik dan lalu kita tarik komponen vertikal pada table menuju grafik (Gambar 4.3).
24
Gambar 4.3 Contoh tampilan ed iting gelo mbang untuk ko mponen vertikal dari masingmasing stasiun kedalam grafik pada data tanggal 6 April 2011.
4.5. Filtering data gelombang Tahap selanjutnya yaitu filtering, hal ini dilakukan apabila gelombang banyak noise dan sulit untuk melakukan picking. Filtering dilakukan untuk memisahkan gelombang yang yang tidak di inginkan (Noise). Filtering dilakukan dengan menggunakan bandpass filter antara 1-5 hz. Hasilnya dapat terlihat pada gambar 4.4a untuk sebelum filtering dan gambar 4.4b setelah filtering.
Gambar 4.4a Contoh data sebelum Filtering untuk data pada data tanggal 6 April 2011 untuk lima stasiun pengamatan dengan komponen vertikal.
25
Gambar 4.4b Contoh data setelah Filtering untuk data pada data tanggal 6 April 2011 untuk lima stasiun pengamatan dengan komponen vertikal.
4.6.
Pemilihan event data gelombang Pada pemilihan data event ini dilakukan dengan memperbesar time dan
juga amplitudenya sehingga event terlihat lebih jelas. Untuk memperbesarnya tinggal menggeser time dan amplitude pada windows yang terletak dibawah gelombang (Gambar 4.6) apabila event yang dilihat sudah tepat, seperti yang dijelaskan pada teori, yaitu setiap stasiun pengamatan akan memberikan respon yang sama. Setelah event yang dipilih dirasa tepat maka data tersebut dipotong dengan meng-klik waveform pada menubar kemudian cut (Gambar 4.5), lalu pilih this time. Hasil dari pemilihan event (Gambar 4.7) ini akan mempermudah nantinya dalam proses picking.
Gambar 4.5 Menu untuk pemotongan sinyal
26
Gambar 4.6 Kotak dengan garis merah digunakan untuk memperbes ar dan memperjelas sinyal untuk contoh untuk data pada data tanggal 6 April 2011.
Gambar 4.7 Contoh Sinyal yang telah diperjelas untuk data pada data tanggal 6 April 2011.
4.7.
Konve rsi data gelombang Setelah melakukan pemotongan tahap berikutnya, yaitu mengeksport data
dari type file .seed menjadi file dengan type .mseed. untuk mengeksport file tersebut digunakan software geopsy, open software geopsy, kemudian file-import
27
signal-file lalu open sinyal tersebut. Setelah membuka file eksport data dapat dilakukan dengan cara file-eksport atau langsung dengan mengklik ctrl E secara bersamaan. Untuk type file yang dieksport dapat diubah menjadi beberapa macam seperti .arr, .saf, .su, .mseed dll (Gambar 4.8). Untuk penelitian kali ini menggunakan type file .mseed.
Gambar 4.8 Type file eksport
Selain itu juga terdapat beberapa pilihan juga dapat menggunakan nama dari file sebelumnya tanpa mengetik ulang, yaitu dengan mencentang pilihan use original base name (Gambar 4.9).
Gambar 4.9 Type file eksport.
4.8.
Picking amplitude mkasimum Untuk pengepick-an data gempa
menggunakan software seisgram2k,
caranya meng-klik file pada menu bar software seisgram2k kemudian open file (Gambar 4.10) lalu buka file hasil pemotongan pada pemilihan data yang sudah diubah ekstensinya menjadi .mseed kemudian klik open.
28
Gambar 4.10 Menu open file pada software Seisgram2k
Setelah file tebuka picking dapat dilakukan dengan meng-klik pick pada menu bar kemudian dapat pula mengepicknya secara manual atau otomatis tinggal meng-klik P-P max pada menu bar. P-P max (Gambar 4.11a)disini untuk mencari nilai amplitude maksimumnya dalam rangka mengetahui nilai Peak Ground Velocity.
Gambar 4.11a Contoh Picking gelo mbang pada lima stasiun pengamatan mencari nilai amp litude maksimu m pada data tanggal 6 April 2011.
29
Amplitude maksimum dan label waktu
Gambar 4.11b Contoh Gambar diperbesar untuk P-P max merupakan a mplitude maksimu m Stasiun Pengamatan 2 pada data tanggal 6 April 2011.
Setelah mendapatkan amplitude maksimumnya beserta label waktu (Gambar 4.11b) kemudian menyamakan dengan data informasi katalog, kemudian kita dapat
menghitung
nilai Peak
Ground
Velocitynya, yaitu dengan
mengkonversinya menjadi cm/s, yaitu dengan mengkalikanya dengan 10-9 .
BAB V PEMBAHASAN 5.1 Picking Event Proses picking data dilakukan dengan menggunakan software Seisgram2k, dengan menggunakan data rekaman gelombang dari 5 stasiun pengamatan. datadata yang terekam merupakan data-data yang sebagian besar akibat dari pergerakan tanah pada zona subduksi yang terekam oleh stasiun pengamatan di kota Padang, Sumatera Barat. Pergerakan tanah tersebut berupa gempa- gempa baik akibat gejala tektonik maupun vulkanik. Proses picking disini berbeda dengan biasanya karena bukan p icking pada first time arrival pada fase gelombang P maupun fase gelombang S, melainkan pengepick-an dilakukan pada maksimum amplitudonya. Picking dilakukan dengan memperhatikan waktu event pada setiap stasiun untuk jarak stasiun yang relatif dekat (Gambar 5.1) maka retang selisih waktunya sekitar orde kurang dari 1 menit. Misalnya saja pada stasiun pengamatan 2 (Sikabu) 20110225 2209 53.880, Stasiun pengamatan 4 (Air Dingin) 20110225 2209 46.990, Stasiun pengamatan 5 (Sungai Lareh) 20110225 2209 44.290, Stasiun Pengamatan 6 (Bukit Gado- gado) 20110225 2209 39.400, Stasiun pengamatan 7 (Tarantang) 20110225 2209 41.530. Proses picking tersebut menghasilkan nilai amplitude maksimum yang masih dalam satuan nanometer kemudian di konversi menjadi Peak Ground Velocity (PGV). Jarak episenter didapatkan dari hasil kalkulasi antara koordinat stasiun pengamatan dengan koordinat sumber gempa dari sumber informasi Katalog global dari bulan Februari sampai dengan April 2011. Dalam perangkat lunak seisgram2k terdapat beberapa tools yang digunakan untuk melakukan picking amplitude maksimum secara otomatis, diantaranya P-P_max dan 0-P_max. P-P_max digunakan untuk picking amplitudo maksimal dari lembah sampai puncak gelombang atau disebut dengan peak to peak amplitude. Sedangkan 0-P_max digunakan untuk picking amplitudo maksimal yang berasal dari posisi 0 sampai puncak gelombang, ini berarti pada
30
31
tombol 0-P_max, amplitude yang diambil adalah amplitudo dengan puncak tertinggi. Sedangkan pada mode P-P_max, amplitudo yang diambil adalah amplitudo puncak ditambah amplitudo lembah maksimal.
Gambar 5.1 Gempa pada tanggal 25 Februari 2011 pada pukul 22.09 tercatat di lima stasiun.
Selain dilakukan secara otomatis, picking juga dapat dilakukan dengan cara manual. Cara manual ini dilakukan jika hasil picking dengan cara otomatis tidak bisa dipercaya atau kurang meyakinkan. 5.2 Grafik Peak Ground Velocity Vs Jarak Episentrum Setelah melalui tahap picking kemudian kita menginterpretasikanya ke dalam sebuah grafik, yaitu peak ground velocity vs jarak episentrum. Untuk jarak episentrum tiap stasiun pengamatan terhadap event-event tidaklah sama karena koordinat setiap stasiun berbeda walau jika di lapangan jaraknya tidak jauh. Perbedaan ini pula menjadikan nilai peak ground velocity di tiap stasiun memiliki nilai yang berbeda-beda, juga waktunya juga berbeda, tetapi dalam selisih waktu pada setiap stasiun tidak terpaut jauh yaitu kurang 1 menit.
32
Pada penelitian ini menggunakan metode regresi untuk mendapatkan persamaan pergerakan kecepatan tanah maksimum. Dengan menggunakan rumus Molas dan Yamazaki (1995) sebagai acuan, hal tersebut karena dalam penelitian ini jarak episentrum gempa dengan stasiun pengamatan berada di kisaran 20 s/d 350 Km. Sedangkan untuk rumus Megawati et. al.(2004) sebenarnya rumusnya dapat dikategorikan mendekati karena untuk penelitian zona subduksi di Pulau Sumatera, tetapi dalam pengukuran dilakukan dari singapura dan untuk jarak ketelitian berada pada jangkauan 200-600 km, dalam data penelitian bisa saja menggunakan rumus Megawati et. al.(2004) tetapi menjadi kurang valid, karena sebagian besar data berada pada jarak kurang dari 250 km. Pada runus Megewati et. al. (2004) untuk data yang berada kurang dari 250 km tidak bisa di kalkulasi. Pada persamaan Megawati et. al.(2004) juga digunakan untuk gempa subduksi dengan moment magnitude jangkauan 5.0 s/d 6.5 jadi untuk data dengan nilai magnitude diatas maupun dibawah nilai tersebut masih belum dapat di jangkau. Faktor lain bahwa penelitian Megawati et. al.(2004) berada di Singapura dan Kuala Lumpur yang merupakan daerah dengan nilai seismisitas yang rendah, sedangkan untuk penelitian kali ini berada pada daerah dengan nilai seismisitas yang tinggi. Berdasarkan faktor tersebut karena jarak sangat mempengar uhi maka tidak digunakanya rumus tersebut. Demikian halnya untuk persamaan Atkinson dan Boore (1997) alasan tidak digunakanya karena ketelitian jarak yang dapat dijangkau berada pada jangkauan 100-400km. Pada penelitian Atkinson dan Boore (1997) ini berada di daerah Kaskadia, dengan jangkauan moment magnitude 4.0 s/d 8.5 . Molas dan Yamazaki (1995) meneliti untuk gempa bumi subduksi dengan jarak episenter dengan jangkauan 20 s/d 600 km dan jangkauan moment magnitude 4.0-7.0. Penelitian tersebut dilakukan di Jepang yang merupakan daerah yang mempunyai nilai seismisitas yang tinggi. Berikut grafik (Gambar 5.2a, Gambar 5.2b, Gambar 5.2c, Gambar 5.2d, Gambar 5.2e, Gambar 5.3a, Gambar 5.3b, Gambar 5.3c, Gambar 5.3d, Gambar 5.3e, Gambar 5.4a, Gambar 5.4b, Gambar 5.4c, Gambar 5.4d, Gambar 5.4e) hasil
33
dari pengolahan data gempa yang terekam oleh stasiun pengamatan di kota Padang, Sumatera Barat pada periode bulan Febuari sampai dengan April 2011 :
10 4<M<=5 SP2
1 1.0
10.0
100.0
1000.0
LOG PGV ( cm/s )
0.1
4<M<=5 Hasil Regeresi 4<M<=5 MY 95
0.01
SP 2
0.001 0.0001
SP 2 Hasil Regresi
0.00001
MY 95
0.000001
Distance ( Km )
Gambar 5.2a Grafik Peak Ground Velocity simu lasi horizontal Stasiun pengamatan 2 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 4 s/d 5 .
10
4<M<=5 SP4 1
LOG PGV ( cm/s )
1.0
10.0
100.0
0.1
1000.0
4<M<=5 Hasil Regeresi 4<M<=5 MY 95
0.01
SP 4 0.001
SP 4 Hasil Regresi
0.0001
MY 95
0.00001
Distance ( Km )
Gambar 5.2b Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 4 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 4 s/d 5 .
34
10 4<M<=5 SP5 1 LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
0.1
4<M<=5 Hasil Regeresi 4<M<=5 MY 95
0.01
SP 5 0.001
SP 5 Hasil Regresi 0.0001
MY 95
0.00001
Distance ( Km )
Gambar 5.2c Gra fik Peak Ground Velocity simu lasi horizontal Stasiun pengamatan 5 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 4 s/d 5 .
100 4<M<=5 SP6 10
LOG PGV ( cm/s )
1
1.0
10.0
100.0
0.1
SP 6
0.01 0.001
SP 6 Hasil Regresi
0.0001 0.00001
1000.0
4<M<=5 Hasil Regeresi 4<M<=5 MY 95
MY 95 Distance ( Km )
Gambar 5.2d Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 6 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 4 s/d 5 .
35
10 4<M<=5 SP7
1 LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
0.1
4<M<=5 Hasil Regeresi
4<M<=5 MY 95
0.01
SP 7 0.001 SP 7 Hasil Regresi
0.0001
MY 95
0.00001
Distance ( Km ) Gambar 5.2e Grafik Peak Ground Velocity simu lasi horizontal Stasiun pengamatan 7 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 4 s/d 5 .
10 5<M<6 SP2 1
LOG PGV ( cm/s )
1.0
10.0
100.0
0.1
1000.0
5<M<6 Hasil Regeresi
5<M<6 MY 95 0.01 SP 2
0.001 SP 2 Hasil Regresi
0.0001 0.00001
MY 95 Distance ( Km )
Grafik 5.3.a Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 2 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 5 s/d 6.
36
10
5<M<6 SP4
1 LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
0.1
5<M<6 Hasil Regeresi 5<M<6 MY 95 SP 4
0.01
SP 4 Hasil Regresi
0.001
MY 95 0.0001
Distance ( Km )
Grafik 5.3.b Grafik Peak ground velocity simulasi horizontal Stasiun pengamatan 4 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 5 s/d 6.
10 5<M<6 SP5 1
LOG PGV ( cm/s )
1.0
10.0
100.0
0.1
1000.0
5<M<6 Hasil Regeresi 5<M<6 MY 95
0.01
SP 5 0.001
SP 5 Hasil Regresi
0.0001 0.00001
MY 95 Distance ( Km )
Grafik 5.3.c Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 5 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 5 s/d 6.
37
10 5<M<6 SP6
1 LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
0.1
5<M<6 Hasil Regeresi 5<M<6 MY 95 SP 6
0.01
SP 6 Hasil Regresi
0.001
MY 95 0.0001
Distance ( Km )
Grafik 5.3.d Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 6 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 5 s/d 6.
10 5<M<6 SP7 1 LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
5<M<6 Hasil Regeresi
5<M<6 MY 95
0.1
SP 7 0.01 SP 7 Hasil Regresi
0.001
MY 95 0.0001
Distance ( Km )
Grafik 5.3.e Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 7 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 5 s/d 6.
38
10 6<M<7 SP4 1
LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
0.1
6<M<7 Hasil Regeresi 6<M<7 MY 95
0.01
SP 4 0.001
SP 4 Hasil Regresi
0.0001
MY 95 0.00001
Distance ( Km )
Grafik 5.4.a Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 2 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 6 s/d 7.
10
6<M<7 SP4
1
LOG PGV ( cm/s )
1.0
10.0
100.0
0.1 0.01
6<M<7 Hasil Regeresi 6<M<7 MY 95 SP 4
0.001
SP 4 Hasil Regresi
0.0001 0.00001
1000.0
MY 95
Distance ( Km )
Grafik 5.4.b Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 4 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 6 s/d 7.
39
10
6<M<7 SP5
1
LOG PGV ( cm/s )
1.0
10.0
100.0
1000.0
0.1
0.01
6<M<7 Hasil Regeresi 6<M<7 MY 95
SP 5
0.001
SP 5 Hasil Regresi
0.0001
MY 95
0.00001
Distance ( Km )
Grafik 5.4.c Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 5 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 6 s/d 7.
10 6<M<7 SP6 1
LOG PGV ( cm/s )
1.0
10.0
100.0
0.1
1000.0
6<M<7 Hasil Regeresi 6<M<7 MY 95
0.01
SP 6 0.001 SP 6 Hasil Regresi
0.0001
MY 95 0.00001
Distance ( Km )
Grafik 5.4.d Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 6 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 6 s/d 7.
40
10 6<M<7 SP7
1 1.0
10.0
100.0
LOG PGV ( cm/s )
0.1
6<M<7 Hasil Regeresi 6<M<7 MY 95
0.01
SP 7
0.001 0.0001
SP 7 Hasil Regresi
0.00001 0.000001
1000.0
MY 95 Distance ( Km )
Grafik 5.4.e Grafik Peak Ground Velocity simulasi horizontal Stasiun pengamatan 7 untuk gempa subduksi Sumatera barat dengan range mo ment magnitude 6 s/d 7.
Grafik diatas merupakan grafik peak ground velocity simulasi horizontal untuk gempa subduksi Sumatera barat stasiun pengamatan 2 ( Sikabu ), stasiun pengamatan 4 (Air Dingin ), stasiun pengamatan 5 ( Sungai Lareh ), stasiun pengamatan 6 ( Bukit Gado-Gado ), stasiun pengamatan 7 ( Tarantang ) serta grafik hasil regresi menggunakan persamaan Molas dan Yamazaki (1995). Dari grafik diatas berada pada jarak episenter dengan jangkauan 20 s/d 350 km dan moment magnitude dengan jangkauan 4.0 s/d 7. Masing-masing dari stasiun di plotkan dalam warna garis berbeda untuk respon peak ground velocity dari tiap masing- masing stasiun pengamatan. Nilai respon peak ground velocity tiap stasiun berbeda semua, tetapi semuanya menggambarkan respon yang kurang lebih sama, hal tersebut terlihat dalam grafik. Respon peak ground velocity dari tiap stasiun menunjam dari arah kanan gambar menuju arah kiri gambar, kemudian garis hitam lurus menggambarkan nilai respon peak ground velocity dari hasil menggunakan rumus Molas dan Yamazaki, sedangkan garis hitam putus-putus menggambarkan nilai respon peak ground velocity yang telah di regresi. Berdasarkan grafik diatas hasil dari metode regresi maka didapat persamaan atenuasi pergerakan tanah sebagai berikut :
41
(5.1) Dari hasil persamaan diatas didapatkan mendapatkan nilai baru dari hubungan atenuasi dari Peak Ground Velocity (PGV) untuk gempa bumi subduksi di Kota Padang atau Indonesia secara umum. Untuk memperkirakan seismic hazard yang relevan secara umum mempunyai daerah jangkauan dengan jarak dari 20 s/d 350 km dan moment magnitude dengan jangkauan 4 s/d 7. Berdasarkan perhitungan
hasil regeresi respon atenuasi pergerakan kecepatan tanah
maksimum, dari perhitungan didapatkan nilai koefisien regeresi yang di aplikasikan pada persamaan dan menghasilkan persamaan (14) diatas. Nilai koefisien tersebut menggambarkan kemungkinan tingkat kerusakan yang terjadi akibat kejadian subduksi besar dengan nilai magnitudo yang tinggi yang mungkin dapat terjadi di kota Padang, Sumatera Barat untuk kejadian pada jarak 350 km.
BAB VI KESIMPULAN DAN SARAN 6.1 Kesimpulan Dalam penelitian ini berdasarkan pada grafik peak ground velocity simulasi horizontal untuk gempa subduksi Sumatera Barat SP 2 ( Sikabu ), SP 4 (Air Dingin ), SP 5 ( Sungai Lareh ), SP 6 ( Bukit Gado-Gado ), SP 7 ( Tarantang ) didapatakan persamaan atenuasi pergerakan tanah maksimum sebagai berikut :
(6.1) Data perekaman gelombang gempa dari rumus tersebut didapat pada periode Februari s/d April 2011. Persamaan 6.1 dapat digunakan untuk mendapatkan informasi dari respon atenuasi pergerakan tanah dengan jarak episenter dengan jangkauan 20 s/d 350 km dan moment magnitude dengan jangkauan 4.0 s/d 7. Dari rumus diatas dapat digunakan untuk menentukan besar bahaya yang terjadi apabila terjadi gempa subduksi. 6.2 Saran 1. Dalam penelitian alangkah baiknya apabila menghitung juga jarak dan magnitude gempa dari hasi pengepick-an dan membandingkanya dengan data yang diperoleh dari informasi kalalog gempa, sehingga didapat perbedaanya dan dapat diperhitungkan kembali nilai ralat erornya. 2. Hal lain juga alangkah baiknya jika data yang diolah banyak, sehingga jangkauan magnitude dan jangkauan jarak episenternya mungkin dapat bertambah lebar.
42
DAFTAR PUSTAKA Abdullah, A. 2007. Ensiklopedi Seismik (Online). Tersedia : http://ensiklopediseismik.blogspot.com/(28januari2010). Anonymous. (2007). Arahan Pemanfaatan Ruang Wilayah Rawan Bencana Sumatera Barat Tahun 2003-2019. Sumatera Barat: Bappeda Smatera Barat. 17 Maret 2011. Adrian Rodrigues,Marek. Near-fault Seismic Site Response. Engineering-Civil Engineering, University of California, Berkeley. Atkinson GM, Boore DM. Stochastic point-source modeling of ground motions in the Cascadia Region. Seismological Res Lett 1997; 68(1):74–85. Bakosurtanal. (2005). Peta Rupa Bumi Provinsi Sumatera Barat. Boore, D. M., Joyner, W. B., & Fumal, T. E. 1994a. Estimation of response spectra and peak accelerations from western North American earthquakes: An interim report. Part 2. Open-File Report 94-127. U.S. Geological Survey. Cornell, C.A. (1968). Engineering seismic risk analysis, Bull. Seism. Soc. Am., 58, 1583-1606. Darman, H., and Sidi, F.H., eds., 2000, An outline of the geology of Indonesia; Indonesian Association of Geologists, 192 p. Douglas, J., 2011. BRGM. Ground Motion Prediction Equation 1964-2010. Pacific Earthquake Engineering Research Center College of Engineering University of California, Berkeley. Handoko, E.Y.,2004, Pendidikan dan Pelatihan (Diklat)Teknis Pengukuran dan Pemetaan Kota. Fakultas Teknik Sipil dan Perencanaan, Institut Teknologi Sepuluh November Irsyam, M., Sengara, W.I, Aldlamar, F., Widiyantoro, S., Triyoso, W., Hilma, D., Kertapati, E., Mellano, I., Asrurifka, M., Ridwan, M., Suhardjono. 2010. Ringkasan Hasil Studi Tim Revisi Peta Gempa Indonesia 2010 ( Edisi 2). Bandung, Juni 2010. Kanai, K. 1966. Improved empir ical formula for character istics of stray [sic] ear thquake mo-tions. Pages 1–4 of: Proceedings of the Japanese Ear thq uake Sym posium . N ot seen. Repor ted in Tr ifunac & Brady (1975). McGuire, R., 1976, Fortran Computer Program for Seismic Risk Analysis, Open-File Report 76-67, U.S. Geological Survey. Megawati, K., 2004, Response Spectral Attenuation Relationships for Sumatran-Subduction Earthquakes and The Seismic Hazard Implications to Singapore and Kuala Lumpur. Centre for Earthquake Engineering Research, Department of Civil Engineering, The University of Hongkong. Milsom, J., 2002, Field Geophysics:The Geological Field Guide Series. Citabook; The atrium Gate, Chicester, West Sussex, England.
Molas GL, Yamazaki F. Attenuation of earthquake ground motion in Japan including deep focus events. Bull Seismological Soc Am 1995; 85(5):1343–58. Nakamura, Y., (1989). A Method for Dynamic Characteristics Estimation of Subsurface usingMicrotremor on the Ground Surface. Quarterly Report of RTRI, Railway Technical Research Institute(RTRI), 30(1). Nogoshi, I., Igarashi, T.,(1971). On the amplitude characteristics of microtremor (Part 2), Jour. Seism. Soc. Japan, 24, 24-60 (in Japanese with English abstract). Pusat Vulkanologi Dan Mitigasi Bencana Geologi.2008. Wilayah Potensi Gerakan Tanah Di Provinsi Sumatra Barat Bulan Oktober 2008.Bandung. http://www.bgl.esdm.go.id/dmdocuments. Reid, H.F., The Mechanics of the Earthquake: The California Earthquake of April 18, 1906., Report of the State investigation Committee., Caregie Institution of Washington, D.C, 1910; Vol. 2 Salwan, S.,(2009). Perhitungan tingkat bahaya gempabumi di Jawa Tengah secara probabilistic. Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Gadjah Mada. Sandy, I Made. (1996). Republik Indonesia: Geografi Regional. Depok: Departemen Geografi, Universitas Indonesia. Saputro, N.E.,(2007). Pemetaan peak ground acceleration akibat dari gempa utama Yogyakarta 27 Mei 2006 dan susulanya dengan menggunakan metode Kanai (1966) pendekatan garis untuk sumber garis untuk gempa utama dan sumber titik untuk gempa susulan. Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Gadjah Mada. Tong, H., & Katayama, T. 1988. Peak acceleration attenuation by eliminating the ill-effect of the correlation between magnitude and epicentral distance. Pages 349–354 of: Proceedings of Ninth World Conference on Earthquake Engineering, vol. II. Wallace, R. L. & R. A. Colburn, 1989. Phylogenetic relationships within phylum Rotifera: orders and genus Notholca. Hydrobio-logia 186/187: 311–318. Wassermann, J., 2002, IASPEI, New Manual Of Seismological Observatory, Daleh.Id.Au, Diakses 10 Agustus 2011.
