PEMETAAN TINGKAT RESIKO TSUNAMI DI KABUPATEN SIKKA NUSA TENGGARA TIMUR DENGAN MENGGUNAKAN SISTEM INFORMASI GEOGRAFIS

PEMETAAN TINGKAT RESIKO TSUNAMI DI KABUPATEN SIKKA NUSA TENGGARA TIMUR DENGAN MENGGUNAKAN SISTEM INFORMASI GEOGRAFIS

Oleh : Ernawati Sengaji C64103064

DEPARTEMEN ILMU DAN TEKNOLOGI KELAUTAN FAKULTAS PERIKANAN DAN ILMU KELAUTAN INSTITUT PERTANIAN BOGOR 2009

PERNYATAAN MENGENAI SKRIPSI DAN SUMBER INFORMASI Dengan ini saya menyatakan bahwa Skripsi yang berjudul : ”PEMETAAN TINGKAT RESIKO TSUNAMI DI KABUPATEN SIKKA NUSA TENGGARA TIMUR DENGAN MENGGUNAKAN SISTEM INFORMASI GEOGRAFIS” adalah benar merupakan hasil karya sendiri dan belum diajukan dalam bentuk apapun kepada perguruan tinggi manapun. Semua sumber data dan informasi yang berasal atau dikutip dari karya yang diterbitkan maupun tidak diterbitkan dari penulis lain telah disebutkan dalam teks dan dicantumkan dalam Daftar Pustaka di bagian akhir skripsi ini.

Bogor, Januari 2009

Ernawati Sengaji C64103064

RINGKASAN ERNAWATI SENGAJI. PEMETAAN TINGKAT RESIKO TSUNAMI DI KABUPATEN SIKKA NUSA TENGGARA TIMUR DENGAN MENGGUNAKAN SISTEM INFORMASI GEOGRAFIS. DIBIMBING OLEH BISMAN NABABAN dan SYAMSUL BAHRI AGUS. Penelitian dengan judul “Pemetaan Tingkat Resiko Tsunami di Kabupaten Sikka Nusa Tenggara Timur dengan Menggunakan Sistem Informasi Geografis” dilaksanakan dalam beberapa tahap yaitu penyusunan basis data awal pada bulan Agustus 2007, survei lapang bulan Oktober 2007, pengolahan citra, penyusunan basis data akhir dan pengolahan data spasial sampai bulan November 2008. Fokus daerah penelitian adalah wilayah Kabupaten Sikka yang berada pada daratan Pulau Flores. Secara administratif Kabupaten Sikka terletak pada 08º22'-08º50' Lintang Selatan dan 121º55'40"-122º41'30º Bujur Timur. Penelitian ini bertujuan untuk memetakan tingkat resiko tsunami di wilayah Kabupaten Sikka dengan menggunakan sistem informasi geografis. Metode yang diterapkan adalah metode Cell Based Modelling (CBM). Analisis spasial pada data raster merupakan dasar dari metode CBM. Pada penelitian ini, pemetaan resiko tsunami ditentukan oleh faktor kerawanan (hazard) dan kerentanan (vulnerability). Parameter kerawanan yang digunakan adalah data run up Tsunami Flores Tahun 1992. Sedangkan parameter kerentanan terhadap tsunami meliputi elevasi daratan, slope, morfometri pantai, landuse, jarak dari garis pantai dan jarak dari sungai. Berdasarkan metode CBM keseluruhan parameter tersebut diubah ke dalam format data raster menjadi layer-layer raster. Selanjutnya seluruh data tersebut dioverlay dengan metode Weighted Overlay sehingga didapatkan model peta tingkat resiko tsunami Kabupaten Sikka. Pada penelitian ini, dibuat lima (5) kelas tingkat resiko berdasarkan studi literatur dan konsultasi pakar serta pembimbing. Selang kelas untuk masing-masing tingkat resiko tsunami adalah sebagai berikut : kelas resiko sangat tinggi (R5) = 4,201-5,000; kelas resiko tinggi (R4) = 3,4014,200; kelas resiko sedang = 2,601-3,400; kelas resiko rendah (R2) = 1,8012,600; dan kelas resiko sangat rendah (R1) = 1,000-1,800. Berdasarkan hasil analisis spasial dengan metode CBM dapat diketahui bahwa luasan daerah yang beresiko sangat tinggi dan tinggi terhadap bahaya tsunami adalah relatif lebih kecil dari total luas wilayah Kabupaten Sikka (1,7086% pada run up 10 m dan 4% pada run up 20 m). Akan tetapi, lokasi kemungkinan kejadian tsunami tersebut berada pada daerah pemukiman padat penduduk dan banyak terdapat infrastruktur penting sehingga perlu dikembangkan mitigasi tsunami yang komprehensif di Kabupaten Sikka, terutama pada wilayah-wilayah yang beresiko tinggi dan sangat tinggi terhadap tsunami.

© Hak cipta milik Ernawati Sengaji, tahun 2009 Hak cipta dilindungi Dilarang memperbanyak tanpa izin tertulis dari Institut Pertanan Bogor, sebagian atau seluruhnya dalam bentuk apapun, baik cetak, fotokopi, mikrofilm dan sebagainya.

PEMETAAN TINGKAT RESIKO TSUNAMI DI KABUPATEN SIKKA NUSA TENGGARA TIMUR DENGAN MENGGUNAKAN SISTEM INFORMASI GEOGRAFIS

SKRIPSI Sebagai Salah Satu Syarat untuk Memperoleh Gelar Sarjana Perikanan pada Fakultas Perikanan dan Ilmu Kelautan Institut Pertanian Bogor

Oleh: Ernawati Sengaji C64103064

DEPARTEMEN ILMU DAN TEKNOLOGI KELAUTAN FAKULTAS PERIKANAN DAN ILMU KELAUTAN INSTITUT PERTANIAN BOGOR 2009

SKRIPSI Judul Skripsi : PEMETAAN TINGKAT RESIKO TSUNAMI DI KABUPATEN SIKKA NUSA TENGGARA TIMUR DENGAN MENGGUNAKAN SISTEM INFORMASI GEOGRAFIS Nama : Ernawati Sengaji NRP : C64103064

Disetujui,

Pembimbing I,

Pembimbing II,

Dr. Ir. Bisman Nababan, M.Sc. NIP. 131 953 477

Syamsul Bahri Agus, S.Pi., M.Si. NIP. 132 311 912

Mengetahui, Dekan Fakultas Perikanan dan Ilmu Kelautan

Prof. Dr. Ir. Indra Jaya, M.Sc. NIP. 131 578 849

Tanggal Lulus : 27 Januari 2009

KATA PENGANTAR Puji dan syukur penulis haturkan ke hadirat Allah SWT karena atas rahmat dan hidayah-Nya, penulis bisa menyelesaikan skripsi dengan judul “Pemetaan Tingkat Resiko Tsunami di Kabupaten Sikka Nusa Tenggara Timur dengan Menggunakan Sistem Informasi Geografis”. Selama penyusunan skripsi ini, penulis banyak mendapat bantuan dari berbagai pihak. Oleh karena itu, penulis ingin mengucapkan terima kasih yang sebesar-besarnya kepada : 1. Bapak, Mama, Ade Nia dan Ade Muchlas yang telah memberikan kasih sayang, semangat, motivasi dan dukungan tak terhingga kepada penulis hingga detik ini. 2. Dr. Ir. Bisman Nababan, M.Sc. dan Syamsul Bahri Agus, S.Pi., M.Si. selaku pembimbing yang banyak sekali memberikan bimbingan, arahan, masukan dan motivasi semangat kepada penulis selama penyusunan skripsi ini. 3. Muhammad Helmi, M.Si. dan Wiwin, S.Si. atas bimbingan, arahan dan kritik mengenai Cell Based Modelling. 4. Dr. Ir. I Wayan Nurjaya, M.Sc. dan Prof. Dr. Ir. Setyo Budi Susilo, M.Sc. atas masukan ilmu tentang tsunami bagi penulis. 5. Dr. Ir. Subandono Diposaptono, M.Eng. dan Firdaus Agung yang banyak memberikan masukan dalam metode penentuan resiko tsunami. 6. Drs. M. Taufik Gunawan, Dipl.Seis., Drs. Budi Waluyo, Dipl.Seis., Indra Gunawan, S.Si. serta segenap staf Earthquake and Tsunamy BMG atas masukan ilmu kegempaan dan tsunami serta data run up tsunami. 7. Ir. Djoko Hartoyo, Dr. Ir. Fadhli Syamsudin, M.Sc., Dr. Ir. Nawa Suwedi, M.Sc., Yudi Wahyudi, DEA., Dr. Ir. Agus Wibowo, M.Sc. dan Dr. Ir.

Udrekh, M.Sc. dari BPPT Jakarta serta Dr. Ir. Widjongkoko, M.Sc. dan Dr. Ir. Chaeroni, M.Sc. dari BPDP-BPPT Yogyakarta atas masukan ilmu tsunami bagi penulis. 8. A.A.G. Conterius S. Sos. dan Mauritius Minggo, ST., MT. dari Bappeda Sikka atas arahan dan bantuan data spasial Kab.Sikka kepada penulis. 9. Kapten Dede Yuliadi, Serka Edi Gunadi dari Dishidros TNI AL atas penjelasan ilmu batimetri dan data pasut. 10. Teman-teman Geologi UNPAD atas literatur geologi bagi penulis. 11. Seluruh dosen, segenap staf Tata Usaha, dan teman-teman (khususnya ITK Angkatan 40) atas bantuan, motivasi, dan kebersamaannya selama ini. Akhir kata, semoga skripsi ini dapat bermanfaat bagi kita semua.

Bogor, Januari 2009

Ernawati Sengaji

DAFTAR ISI

Halaman

DAFTAR TABEL ........................................................................................................ xii DARTAR GAMBAR .................................................................................................. xiii DAFTAR LAMPIRAN ................................................................................................ xiv 1. PENDAHULUAN .....................................................................................................1 1.1 Latar belakang ...................................................................................................1 1.2 Tujuan ................................................................................................................2 2. TINJAUAN PUSTAKA ............................................................................................3 2.1 Tsunami .............................................................................................................3 2.1.1 Karakteristik tsunami .................................................................................3 2.1.2. Penyebab tsunami ....................................................................................5 2.1.2.1 Gempa tektonik .............................................................................5 2.1.2.2 Letusan vulkanik ...........................................................................6 2.1.2.3 Longsoran (Landslide) ..................................................................7 2.1.3 Klasifikasi tsunami .....................................................................................8 2.1.3.1 Tsunami lokal (near field/local field tsunami) .................................8 2.1.3.2 Tsunami jarak jauh (far field tsunami) ...........................................8 2.2 Penilaian tingkat resiko tsunami (tsunami tisk assesment) .................................9 2.2.1 Analisis daerah rawan tsunami .................................................................11 2.2.2 Analisis kerentanan terhadap tsunami ......................................................12 2.2.2.1 Batimetri wilayah pantai ...............................................................12 2.2.2.2 Morfometri pantai .........................................................................14 2.2.2.3 Keberadaan pulau-pulau penghalang ..........................................16 2.2.2.4 Ekosistem pesisir .........................................................................16 2.2.2.5 Penggunaan lahan .......................................................................17 2.2.2.6 Kepadatan pemukiman ................................................................18 2.2.3 Pengurangan resiko bencana tsunami .....................................................18 2.3 Pemetaan resiko tsunami dengan metode Cell Based Modelling ....................19 2.3.1 Penginderaan jauh sebagai penunjang riset tsunami ...............................19 2.3.2 Aplikasi Cell Based Modelling untuk kajian resiko tsunami .......................22 2.3.2.1 Sistem informasi geografis (SIG) .................................................22 2.3.2.2 Pemodelan spasial dengan Cell Based Modelling ........................24 2.4. Keadaan umum lokasi penelitian ......................................................................26 3. METODE PENELITIAN ..........................................................................................27 3.1 Waktu dan lokasi penelitian ...............................................................................27 3.2 Alat dan bahan ..................................................................................................28 3.2.1 Alat ...........................................................................................................28 3.2.2 Bahan .......................................................................................................28 3.3 Survei lapang ....................................................................................................29 3.4 Metode penelitian ..............................................................................................29

3.5 Pengolahan citra satelit .....................................................................................30 3.5.1 Pengolahan citra awal (Pre processing) ...................................................30 3.5.2 Penajaman citra (Image enhancement) ....................................................32 3.5.2.1 Penajaman citra untuk terumbu karang .......................................32 3.5.2.2 Penajaman citra untuk ekosistem mangrove ................................34 3.6 Klasifikasi citra .................................................................................................34 3.7 Metode penentuan tingkat resiko tsunami .........................................................35 3.7.1 Penentuan kerawanan tsunami ................................................................35 3.7.1.1 Pemetaan seismisitas ..................................................................35 3.7.1.2 Pemetaan daerah rawan tsunami ................................................35 3.7.2 Penentuan tingkat kerentanan terhadap tsunami .....................................37 3.8 Pengolahan data dengan Cell Based Modelling ..............................................40 3.9 Matriks tingkat resiko tsunami Kabupaten Sikka .............................................43 4. HASIL DAN PEMBAHASAN ..................................................................................47 4.1 Pengolahan data spasial dari citra satelit ..........................................................47 4.1.1 Pengolahan citra awal ..............................................................................47 4.1.2 Ekstraksi data citra ...................................................................................48 4.2 Analisis kerawanan tsunami di Kabupaten Sikka ..............................................52 4.2.1 Seismisitas di Kabupaten Sikka ................................................................52 4.2.2 Tingkat kerawanan tsunami di Kabupaten Sikka ......................................54 4.3 Analisis tingkat kerentanan tsunami di Kabupaten Sikka ...................................58 4.3.1 Elevasi daratan ........................................................................................59 4.3.2 Kemiringan daratan (slope) ......................................................................63 4.3.3 Morfometri pantai .....................................................................................66 4.3.4 Penggunaan lahan ...................................................................................69 4.3.5 Jarak dari garis pantai ..............................................................................74 4.3.6 Jarak dari sungai ......................................................................................76 4.4 Faktor-faktor pendukung tingkat kerentanan tsunami ........................................78 4.4.1 Batimetri ...................................................................................................78 4.4.2 Ekosistem pesisir alami Kabupaten Sikka ................................................80 4.4.3 Pulau-pulau kecil penghalang ..................................................................80 4.5 Analisis spasial tingkat resiko tsunami di Kabupaten Sikka ...............................81 4.6 Pemodelan resiko tsunami pada run up tsunami 20 Meter ................................91 5. KESIMPULAN DAN SARAN ..................................................................................95 5.1 Kesimpulan .......................................................................................................95 5.2 Saran ................................................................................................................95 DAFTAR PUSTAKA ...................................................................................................96 LAMPIRAN ................................................................................................................101 RIWAYAT HIDUP ......................................................................................................113

DAFTAR TABEL Halaman 1. Kajian resiko tsunami .................................................................................... 10 2. Efektivitas peredaman tsunami hutan pantai ................................................. 17 3. Karakteristik Citra Landsat 7 ETM+ .............................................................. 21 4. Hubungan tingkat kerusakan dengan tinggi run up tsunami .......................... 37 5. Matriks resiko tsunami .................................................................................. 43 6. Luasan ekosistem pesisir Kabupaten Sikka dari citra satelit ......................... 51 7. Luasan daerah kerawanan tsunami .............................................................. 58 8. Luasan wilayah tingkat kerentanan elevasi ................................................... 62 9. Luasan wilayah kerentanan slope ................................................................. 65 10. Luasan jenis landuse di Kabupaten Sikka ................................................... 70 11. Luasan tingkat kerentanan landuse ............................................................ 73 12. Jumlah sel per parameter hasil klasifikasi ................................................... 83 13. Jumlah sel hasil weighted overlay ............................................................... 83 14. Tingkat resiko tsunami per kecamatan ........................................................ 86 15. Luasan wilayah tingkat resiko tsunami pada run up 20 meter ..................... 93

DAFTAR GAMBAR Halaman 1. Hubungan panjang, kecepatan, dan tinggi tsunami di laut ............................. 4 2. Proses terjadinya tsunami akibat gempa tektonik ........................................... 5 3. Lingkar api di daerah Pasifik .......................................................................... 7 4. Tinggi gelombang tsunami saat memasuki teluk ........................................... 13 5. Contoh kawasan pesisir yang menyerupai gigi gergaji .................................. 15 6. Pengurangan resiko bencana ....................................................................... 18 7. Konfigurasi Satelit Landsat 7/ETM+ .............................................................. 20 8. Perbandingan antara data raster dan data vektor ......................................... 23 9. Ilustrasi operasi piksel ................................................................................... 25 10. Peta daerah penelitian ................................................................................ 27 11. Diagam alir penelitian ................................................................................. 30 12. Diagam alir pengolahan citra untuk pemetaan ekosistem pesisir ................ 33 13. Desain penentuan tingkat resiko tsunami .................................................... 42 14. Peta ekosistem pesisir Kabupaten Sikka .................................................... 51 15. Peta seismisitas Kabupaten Sikka .............................................................. 53 16. Peta tinggi run up Tsunami Flores tahun 1992 di berbagai lokasi di Pantai Utara Kabupaten Sikka ............................................................................... 55 17. Peta kerawanan tsunami berdasarkan data run up Tsunami Flores tahun 1992............................................................................................................. 56 18. Peta kerentanan elevasi terhadap tsunami ................................................. 61 19. Peta kerentanan slope terhadap tsunami .................................................... 64 20. Peta kerentanan morfologi terhadap tsunami .............................................. 68 21. Peta jenis landuse di Kabupaten Sikka ........................................................ 70

22. Peta kerentanan landuse terhadap tsunami ................................................ 72 23. Peta jarak dari garis pantai ......................................................................... 75 24. Peta jarak dari sungai ................................................................................. 77 25. Peta 3D batimetri Kabupaten Sikka ............................................................ 79 26. Peta pulau penghalang yang ada di Kabupaten Sikka ................................ 81 27. Peta tingkat resiko tsunami Kabupaten Sikka ............................................. 85 28. Grafik resiko tsunami per kecamatan .......................................................... 87 29. Peta resiko infrastruktur terhadap tsunami .................................................. 88 30. Peta resiko pemukiman terhadap tsunami .................................................. 90 31. Peta tingkat resiko tsunami pada run up 20 meter ...................................... 92

DAFTAR LAMPIRAN

Halaman 1. Posisi pengamatan di lapanga ...............................................................................101 2. Foto-foto hasil survei lapang ..................................................................................102 3. Data run up Tsunami Flores tahun 1992 ................................................................103 4. Root mean square error koreksi Citra Landsat 7/ETM+ 2006 ................................108 5. Hasil analisis statistik Transformasi Lyzenga .........................................................109 6. Peta tingkat resiko tsunami per kecamatan ............................................................110

1.PENDAHULUAN 1.1. Latar belakang Indonesia merupakan salah satu negara dengan kondisi geologis yang secara tektonik sangat labil karena merupakan daerah pertemuan Lempeng Eurasia, Lempeng Indo-Australia, Lempeng Pasifik, dan Lempeng Laut Filipina (Diposaptono dan Budiman, 2005). Kondisi ini menyebabkan wilayah Indonesia memiliki tingkat kejadian gempa yang tinggi di dunia dan sangat rawan mengalami tsunami. Hampir 90% kejadian tsunami di Indonesia disebabkan oleh gempa tektonik dan sekitar 85% dari kejadian tsunami tersebut terjadi di wilayah Indonesia Timur, termasuk daerah Kabupaten Sikka (Diposaptono dan Budiman, 2006). Kabupaten Sikka yang termasuk dalam wilayah Provinsi Nusa Tenggara Timur (NTT) merupakan kawasan dengan tingkat resiko tsunami yang cukup tinggi karena daerah ini berada dekat dengan zona subduksi lempeng tektonik Australia dan Eurasia serta dipengaruhi oleh sesar-sesar aktif di sepanjang Pulau Flores (Fauzi, 2006). Bencana tsunami merupakan bencana yang bersifat destruktif dan menimbulkan banyak kerugian. Misalnya Tsunami Flores 12 Desember 1992 telah menyebabkan kerugian materi hingga milyaran rupiah dan korban jiwa sekitar 2100 orang. Oleh karena itu, penting sekali dilakukan suatu upaya mitigasi bencana tsunami, yaitu proses mengupayakan berbagai tindakan preventif untuk meminimalkan dampak negatif bencana tsunami yang diperkirakan akan terjadi. Salah satu langkah mitigasi adalah dengan membuat peta resiko tsunami (Soegiharto, 2006). Pemetaan tingkat resiko tsunami harus dilakukan dengan pendekatan multikriteria sesuai dengan daerah kajian. Untuk itu, diperlukan suatu perangkat

analisis yang tepat untuk membuat peta tersebut. Sistem informasi geografis (SIG) merupakan perangkat yang memiliki kemampuan untuk memvisualisasikan tingkat resiko tsunami. Salah satu metode aplikasi SIG yang sering digunakan saat ini adalah metode Cell Based Modelling (CBM). Analisis spasial pada data raster adalah dasar dari CBM, sehingga luasan area hasil analisis dengan metode ini cukup akurat karena tidak mengalami generalisasi (ESRI, 2002). Kajian mengenai pemetaan tingkat resiko tsunami di Kabupaten Sikka dengan metode CBM belum pernah dilakukan. Oleh karena itu, hasil penelitian ini diharapkan dapat digunakan sebagai bahan masukan dan pertimbangan dalam pengambilan keputusan (Decision support system) oleh stake holders dalam kerangka mitigasi bencana tsunami di Kabupaten Sikka.

1.2. Tujuan Tujuan dari penelitian ini adalah memetakan dan menentukan tingkat resiko tsunami dengan menggunakan sistem informasi geografis di Kabupaten Sikka, Provinsi Nusa Tenggara Timur.

2. TINJAUAN PUSTAKA 2.1 Tsunami 2.1.1 Karakteristik tsunami Tsunami merupakan gelombang laut dengan periode panjang yang ditimbulkan oleh gangguan impulsif pada medium laut yang bersifat transien (gelombangnya bersifat sesar) dan nondispersive (fasa gelombang tidak bergantung pada panjang gelombang). Gangguan impulsif tersebut disebabkan oleh aktivitas gempa tektonik, erupsi vukanik atau longsoran bawah laut (Diposaptono dan Budiman, 2006). Tsunami bergerak keluar dari daerah pembangkitannya dalam bentuk serangkaian gelombang. Kecepatannya bergantung pada kedalaman perairan, akibatnya gelombang tersebut mengalami percepatan atau perlambatan sesuai dengan bertambah atau berkurangnya kedalaman dasar laut. Proses ini menyebabkan arah pergerakan gelombang berubah dengan energi gelombang bisa menjadi terfokus atau menyebar (UNESCO-IOC, 2006). Tsunami bersifat unik karena bentuk gelombangnya memanjang sampai ke seluruh kolom air yaitu dari permukaan sampai ke dasar lautan. Karakteristik inilah yang menjadi penyebab gelombang tsunami berbeda dengan gelombang laut yang terjadi karena terpaan angin dan pasang surut yang hanya mengganggu permukaan laut. Hal ini juga yang merupakan penyebab dasar besarnya jumlah tenaga yang dibentuk oleh suatu gelombang tsunami (UNESCO-IOC, 2006). Pada laut dalam, gelombang tsunami dapat bergerak pada kecepatan 500 sampai 1.000 kilometer per jam. Periode tsunami (waktu untuk siklus satu gelombang) bisa berkisar dari beberapa menit hingga satu jam. Saat mendekati pantai, kecepatan gelombang melambat menjadi beberapa puluh kilometer per

jam dan tinggi tsunami bisa mencapai hingga puluhan meter pada garis pantai seperti yang terlihat pada Gambar 1 (UNESCO-IOC, 2006).

Gambar 1. Hubungan panjang, tinggi, dan kecepatan tsunami di laut (UNESCO-IOC, 2006) Tinggi tsunami tersebut disebabkan karena terjadi konversi energi kinetik gelombang menjadi energi potensial. Artinya, kehilangan energi akibat berkurangnya kecepatan ditransfer ke dalam bentuk pembesaran tinggi gelombang (run up) (Diposaptono dan Budiman, 2006). Kecepatan run up ke daratan bisa mencapai 25-100 km/jam. Kecepatan gelombang tsunami ini yang menyebabkan jatuhnya korban jiwa, rusaknya lahan pertanian, sarana prasarana wilayah, dan dataran rendah menjadi tergenang membentuk lautan baru. Kembalinya air ke laut setelah mencapai puncak gelombang (run-down) juga bersifat merusak karena menyeret segala sesuatu kembali ke laut (Diposaptono dan Budiman, 2006).

2.1.2 Penyebab tsunami Secara umum faktor penyebab terjadinya tsunami adalah gempa tektonik letusan vulkanik, dan longsoran (landslide) (BMG, 2007; Fauzi dan Wandono, 2005; Diposaptono dan Budiman, 2006). 2.1.2.1 Gempa tektonik Gempa tektonik merupakan gempa yang diakibatkan oleh aktivitas tektonik, yaitu pergerakan, pergeseran dan tumbukan lempeng tektonik (Diposaptono dan Budiman, 2006). Gempa ini merupakan penyumbang terbesar terjadinya tsunami.yang biasanya terjadi di zona subduksi atau zona tumbukan antarlempeng tektonik. Pada zona ini, lempeng yang mempunyai berat jenis yang lebih tinggi (lempeng samudera) akan menyusup di bawah lempeng yang berat jenisnya lebih ringan (lempeng benua). Apabila akumulasi tegangan lempeng benua di sekitar tumbukan mencapai batas maksimum maka ujung lempeng samudera akan melenting ke atas yang mengakibatkan terjadinya patahan (sesar) dan selanjutnya menimbulkan tsunami (UNESCO-IOC, 2006; Gambar 2).

Gambar 2. Proses terjadinya tsunami akibat gempa tektonik (UNESCO-IOC, 2006)

Secara umum karakteristik gempa yang dapat menimbulkan terjadinya tsunami adalah sebagai berikut : 1. Lokasi pusat gempa (episenter) berada di laut dengan magnitudo lebih besar dari 6.0 SR. 2. Kedalaman pusat gempa (hiposenter) kurang dari 60 km dari permukaan bumi (gempa dangkal). 3. Mekanisme patahan gempa tektonik bertipe sesar naik (reverse fault) atau sesar turun (normal fault) (BMG, 2007; International Tsunami Information Center, 2006; Himpunan Ahli Geologi Indonesia, 2006 dan Diposaptono dan Budiman, 2006). Wilayah NTT merupakan kawasan aktif gempa tektonik karena berada dekat lempeng Australia yang menyusup ke lempeng Eurasia. Selain itu, terdapat sesar busur muka dan sesar busur belakang di bagian utara Pulau Flores dan sesar lokal yang memotong Pulau Flores. Sesar busur muka (fore arc) merupakan daerah yang berbatasan langsung dengan zona tumbukan atau sering di sebut sebagai zona aktif akibat patahan yang biasa terdapat di darat maupun di laut. Sesar busur belakang (back arc) merupakan bagian paling belakang dari rangkaian busur tektonik. Sesar atau patahan merupakan perpanjangan gaya yang ditimbulkan oleh gerakan-gerakan lempeng utama yang akan menghasilkan gempa bumi di daratan dan tanah longsor (Wah, 2006). 2.1.2.2 Letusan vulkanik Aktivitas erupsi vulkanik dapat juga menyebabkan tsunami, namun frekuensinya sangat jarang terjadi karena kekuatan aktivitas vulkanik biasanya tidak terlalu besar dan hanya bersifat lokal (Diposaptono dan Budiman, 2006). Untuk Indonesia, daerah yang berpotensi mengalami tsunami akibat letusan

vulkanik merupakan bagian dari lingkar api pasifik (The Pasific Ring of Fire) yang bermula dari Kamchatka Alaska, Jepang, Sumatera, Jawa, Bali, Lombok, Flores, Sulawesi dan berakhir di Filipina (Gambar 3; USGS, 2006).

Gambar 3. Lingkar api di daerah Pasifik (USGS, 2006) Menurut data historis, sejak tahun 1814 sampai tahun 2006, dari 25 kali kejadian tsunami di wilayah NTT, hanya dua kejadian yang dipicu oleh letusan gunung berapi yakni meletusnya Gunung Rokatenda tahun 1928 dan longsornya Gunung Werung tahun 1979 (Wah, 2006). 2.1.2.3 Longsoran (Landslide) Longsoran (landslide) dapat mengakibatkan terjadinya tsunami dan biasanya longsoran ini terjadi akibat gempa bumi yang kuat, letusan vulkanik, longsor di dasar laut dan longsor di atas permukaan laut (sea level) (Diposaptono dan Budiman, 2006).

2.1.3 Klasifikasi tsunami Berdasarkan jaraknya, tsunami diklasifikasikan menjadi 2, yaitu: tsunami lokal (near field/local field tsunami) dan tsunami jarak jauh (far field tsunami) (Puspito, 2007 dan ITIC, 2006). 2.1.3.1 Tsunami lokal (near field/local field tsunami) Tsunami lokal adalah tsunami yang terjadi bilamana jarak antara pusat gempa dan daerah bencana tsunami kurang dari 100 km (ITIC, 2006). Dari segi waktu terjadinya tsunami, tsunami lokal terjadi antara 5 sampai 40 menit setelah gempa utamanya. Hal ini menyebabkan bahwa secara teoritis kejadian tsunami lebih mudah diprediksi dibandingkan dengan kejadian gempa (Puspito, 2007). Adanya tenggang waktu antara terjadinya gempa dan tibanya tsunami di pantai memungkinkan tindakan untuk dapat menganalisis karakteristik apakah suatu gempa dapat menimbulkan tsunami atau tidak. Secara umum tsunami yang terjadi di Indonesia adalah tsunami lokal dan mengingat sistem informasi di Indonesia belum memadai maka biasanya sebelum informasi kejadian tsunami sampai ke masyarakat, gelombang tsunami telah menyapu pantai. Hal ini menyebabkan Indonesia belum bisa memaksimalkan sistem peringatan dini tsunami (Tsunami Early Warning System) (Puspito, 2007). 2.1.3.2 Tsunami jarak jauh (far field tsunami) Tsunami jarak jauh (far field tsunami) adalah tsunami yang diakibatkan oleh gempa laut yang jaraknya ribuan kilometer dari pantai (ITIC, 2006). Waktu datang tsunami berkisar antara beberapa jam sampai 24 jam setelah gempa utamanya. Contoh tsunami jarak jauh ini adalah tsunami Aceh yang terjadi tahun 2004, dimana gelombang tsunami tersebut merambat menyebrangi Samudera Hindia sampai ke Pantai Afrika Selatan.

2.2

Penilaian tingkat resiko tsunami (Tsunami Risk Assesment) Hakekat dari mitigasi bencana tsunami adalah menekan hingga seminimal

mungkin resiko bencana tsunami. Pada dasarnya, resiko sebuah bencana memiliki tiga variabel, yaitu : (1) aspek jenis ancaman, (2) aspek kerentanan, dan (3) aspek kemampuan menanggulangi (Diposaptono dan Budiman, 2006). Dewasa ini banyak terminologi yang digunakan untuk menjelaskan pengertian rawan, rentan dan resiko bencana. Menurut Undang-Undang Republik Indonesia Nomor 24 Tahun 2007 tentang Penanggulangan Bencana, ancaman adalah suatu kejadian atau peristiwa yang bisa menimbulkan bencana. Istilah ini sering disebut juga sebagai kerawanan. Kerentanan (vulnerability) adalah sekumpulan kondisi dan atau suatu akibat keadaan (faktor fisik, sosial, ekonomi dan lingkungan) yang berpengaruh buruk terhadap upaya-upaya pencegahan dan penanggulangan bencana (Badan Koordinasi Nasional Penanganan Bencana, 2007). Resiko bencana adalah potensi kerugian yang ditimbulkan akibat bencana pada suatu wilayah dan kurun waktu tertentu yang dapat berupa kematian, luka, sakit, jiwa terancam, hilangnya rasa aman, mengungsi, kerusakan atau kehilangan harta dan gangguan kegiatan masyarakat (UU RI No.24 Tahun 2007) Hubungan antara kerawanan (hazard), kerentanan (vulnerability), kapasitas penanggulangan (capacity) dan resiko (risk) dirumuskan pada Persamaan 1 di bawah ini (Diposaptono dan Budiman, 2006) :

Re siko ( R ) 

Kerawanan ( H ) x Keren tan an (V ) ........................ (1) Kapasitas Penanggulangan (C )

dimana : R = Resiko; H = kerawanan; V = Kerentanan; dan C = Kapasitas penanggulangan.

Parameter kapasitas penanggulangan tidak dikaji dalam penelitian ini karena terkait dengan kemampuan pemerintah dan sumber daya masyarakat dalam penanggulangan bencana tsunami. Resiko berbanding lurus dengan ancaman atau bahaya (kerawanan) dan tingkat kerentanan terhadap tsunami, serta berbanding terbalik dengan kemampuan (kapasitas) dalam menghadapi tsunami. Semakin besar kerawanan dan kerentanan terhadap tsunami, serta semakin rendah kemampuan penanggulangan dalam menghadapi tsunami, maka akan semakin besar resiko tsunami yang timbul. Bilamana jenis kerawanan tsunaminya sama antara satu daerah dengan daerah lainnya, namun jika tingkat kerentanan dan kapasitas penanggulangan yang berbeda-beda, akan mengakibatkan dampak tsunami yang berbeda, antara satu daerah dengan daerah lainnya (Siahaan, 2008). Salah satu contoh analisis resiko tsunami selengkapnya disajikan pada Tabel 1. Tabel 1. Kajian Resiko Tsunami (Latief, 2005) Risk =

Hazard

Catatan sejarah tsunami

Vulnerability

x

Lingkungan bangunan

Surat kabar (berita)

Infrastruktur

Tide gauge

Peta topografi

Tinggi gelombang

Peta hidrografi

Level run up Stastistik

Potensi kerusakan Kerusakan ikutan (collateral damage) Sosial ekonomi

Tsunami katalog

Perencanaan bencana (Disaster management)

Sumber pembangkit tsunami

Emergency management

Peta waktu tiba tsunami

Local response

TSUNAMI RISK ASSESMENT - TSUNAMI ZONATION MAP MITIGATION POLICY FOR TSUNAMI WARNING

Tsunami termasuk salah satu jenis ancaman yang bersifat alami (natural hazard) dan mempunyai kemungkinan untuk terjadi lagi pada suatu periode tertentu dan pada suatu area geografis tertentu (Damen, 2005 in Melianawati, 2006). Oleh karena itu, untuk melakukan penilaian tingkat resiko tsunami di suatu wilayah tertentu harus dilakukan analisis potensi bahaya tsunami (kerawanan) dan kerentanan terhadap tsunami. 2.2.1 Analisis daerah rawan tsunami Menurut Diposaptono dan Budiman (2006), identifikasi daerah yang berpotensi mengalami bencana tsunami dilakukan dengan beberapa cara yaitu identifikasi jalur pertemuan lempeng (tectonic setting plate) yang berpotensi menyebabkan gempa dan tsunami baik near field maupun far field tsunami, analisis aspek historis kejadian gempa yang berpotensi tsunami, analisis historis kejadian tsunami dan pemodelan tsunami terutama pemodelan run up tsunami. Run up tsunami adalah elevasi air laut vertikal yang dapat dicapai oleh tsunami ke arah darat diukur dari muka air laut rata-rata (mean sea level) atau dari garis pantai pada saat tsunami (Diposaptono dan Budiman, 2006). Data run up merupakan data yang penting sebagai parameter kerawanan dalam kajian resiko tsunami. Pembuatan peta resiko tsunami dapat dibuat melalui tiga pendekatan (Diposaptono dan Budiman, 2006), yaitu : 1. Menggunakan data historis genangan dan run up tsunami yang pernah terjadi sebelumnya. Metode ini digunakan jika ada data historis yang lengkap. Sebagai alternatif, digunakan data citra detail sebelum dan sesudah tsunami untuk membandingkan topografi dan landuse sebelum dan sesudah tsunami. 2. Menggunakan simulasi pemodelan matematik untuk pembangkitan, penjalaran, run up dan inundasi tsunami di wilayah pesisir. Pendekatan

matematik ini masih mengalami kendala utama yaitu menyangkut ketelitian data batimetri dan topografi yang kurang detail. 3. Menggunakan asumsi gelombang tsunami yang mencapai pantai mempunyai ketinggian yang sama diukur dari permukaan laut (metode kesamaan ketinggian). Cara terakhir inilah yang paling mudah ditempuh. Berdasarkan analisis kontur wilayah, maka distribusi luas dan tinggi genangan secara spasial dapat diperoleh dengan mudah. Hasilnya kita bisa menentukan daerah potensi rawan tsunami yang kemudian dioverlay dengan peta infrastruktur, ekosistem pesisir, tataguna lahan dan kondisi sosial ekonomi, sehingga bisa diperoleh peta resiko tsunami seperti yang pernah diterapkan untuk daerah Padang. 2.2.2 Analisis kerentanan terhadap tsunami Kerentanan tsunami dapat dikaji dari berbagai faktor antara lain kerentanan fisik (misalnya jenis bahan dan kekuatan struktur bangunan), kerentanan lingkungan (ketinggian dan morfologi), kerentanan infrastruktur (sarana dan prasarana penting), kerentanan sosial kependudukan (jumlah penduduk dan kepadatan penduduk, struktur penduduk lanjut usia dan balita), kerentanan sosial ekonomi (jumlah/proporsi penduduk miskin dan pengangguran), dan kerentanan kelembagaan (Diposaptono dan Budiman, 2006). Pada penelitian ini, analisis kerentanan yang dikaji adalah lingkungan dan sosial kependudukan. Parameter yang dikaji tingkat kerentanannya terhadap tsunami secara rinci dijelaskan sebagai berikut. 2.2.2.1 Batimetri wilayah pantai Kondisi batimetri sangat mempengaruhi pembelokan atau penerusan gelombang tsunami yang menjalar ke pantai. Pada saat memasuki perairan dangkal, gelombang mengalami gesekan dengan dasar laut.

Hal ini menyebabkan kecepatan dan energi gelombang menurun drastis dengan berkurangnya kedalaman tetapi tinggi gelombang semakin meningkat (Gambar 4; Diposaptono dan Budiman, 2006). Tsunami yang menjalar ke pantai yang sempit dan dangkal akan mengalami proses yang kompleks, yaitu shoaling, refraksi, difraksi dan refleksi. Shoaling adalah proses pembesaran tinggi gelombang karena pendangkalan dasar laut. Tsunami yang terjadi di laut dalam selama penjalarannya ke pantai atau laut yang lebih dangkal akan teramplifikasi karena adanya efek shoaling (Diposaptono dan Budiman, 2006). Refraksi merupakan transformasi (perubahan) gelombang akibat adanya perubahan geometri dasar laut yaitu perubahan kedalaman laut. Sedangkan difraksi adalah transformasi gelombang yang diakibatkan oleh adanya bangunan atau struktur yang menghalangi rambatan tsunami. Rintangan tersebut dapat berupa bangunan pantai atau adanya pulau-pulau kecil (Diposaptono dan Budiman, 2006).

Gambar 4. Tinggi gelombang tsunami saat memasuki teluk (Diposaptono dan Budiman, 2006)

Tinggi gelombang tsunami berubah sesuai dengan perubahan kedalaman. Jika tinggi gelombang lebih kecil daripada panjang gelombangnya maka perubahan tinggi gelombang akan sesuai dengan Hukum Green yaitu perbandingan tinggi tsunami pada posisi di pantai yang lebih dangkal (H1) dengan tinggi tsunami di pantai yang lebih dalam (H0) sama dengan pangkat ¼ dari perbandingan perubahan kedalaman air laut (h0 / h1) dikalikan dengan akar dari perbandingan lebar teluk di titik yang lebih dalam (b0) dan lebar teluk di titik yang lebih dangkal (b1). Secara matematis ditulis sebagai berikut (Diposaptono dan Budiman, 2006) : 1/ 4

H1  h0    H 0  h1 

1/ 2

 b0     b1 

............................(2)

Keterangan : Tinggi gelombang tsunami pada laut dangkal H1 : Tinggi gelombang tsunami pada laut dalam H0 ho/h1 : Perbandingan perubahan kedalaman air laut b0 : Lebar teluk di titik yang lebih dalam : Lebar teluk di titik yang lebih dangkal b1 2.2.2.2 Morfometri pantai Morfometri pantai menunjukkan bentuk dan topografi garis pantai. Morfometri merupakan parameter yang sangat penting untuk dikaji karena menentukkan daerah pemusatan atau penyebaran energi gelombang tsunami (Diposaptono dan Budiman, 2006). Gelombang tsunami yang yang bergerak ke pantai akan mengalami perubahan tinggi gelombang yang semakin besar. Perubahan gelombang tersebut sangat dipengaruhi oleh morfometri pantai sebagai berikut (Hendratno, 2005). 

Geometri pantai (arah horizontal) Secara umum wilayah pesisir Indonesia memiliki teluk yang berbentuk V (V-shape bays) yang berasosiasi dengan tanjung dan muara sungai yang

sangat banyak dan berderet satu dengan yang lainnya, sehingga menyerupai gigi gergaji (sawtooth) (Gambar 5). Morfologi pantai yang berbentuk teluk akan mempengaruhi refraksi gelombang dan menyebabkan peningkatan kecepatan dan energi gelombang tsunami (Latief, 2005).

Gambar 5. Contoh kawasan pesisir yang menyerupai gigi gergaji (Diposaptono dan Budiman, 2006) Gelombang tsunami yang menjalar ke teluk, tidak dapat keluar lagi karena sebagian atau seluruh gelombang tersebut akan dipantulkan oleh dinding teluk. Akibatnya gelombang tsunami akan meninggi dan interaksi gelombang tersebut berlangsung dalam waktu yang lama. 

Kelandaian pantai (arah vertikal) Kelandaian pantai sangat mempengaruhi run-up maksimum ke daratan. Pada pantai yang terjal, tsunami tidak terlalu jauh mencapai daratan karena tertahan dan dipantulkan kembali oleh tebing pantai. Sementara pada pantai relatif landai, gelombang tsunami yang tiba akan semakin membesar dan bertambah cepat (Hendratno, 2005).



Kekasaran pantai Kekasaran pantai dapat berfungsi sebagai peredam energi gelombang tsunami. Pantai yang kasar (berbukit-bukit, berbatu cadas, berkarang, dan atau ditutupi oleh vegetasi) dapat meredam energi gelombang tsunami, sedangkan pantai yang halus (berupa endapan aluvial dan pasir berukuran sedang sampai halus) kurang efektif dalam meredam energi gelombang tsunami (Hendratno, 2005).

2.2.2.3 Keberadaan pulau-pulau penghalang Pulau-pulau penghalang berperan penting dalam meredam energi gelombang tsunami. Energi gelombang tsunami akan tereduksi dengan adanya pendangkalan kedalaman perairan pada pulau penghalang. Sehingga, energi gelombang yang akan dirasakan oleh pulau-pulau di belakangnya akan menjadi lebih kecil dari energi awal (Hajar, 2006). 2.2.2.4 Ekosistem pesisir Ekosistem pesisir terutama mangrove dan hutan pantai memegang peranan sebagai greenbelt pelindung pantai dalam mereduksi energi gelombang tsunami. Hutan pantai sangat efektif dalam meredam energi gelombang tsunami seperti dilaporkan Harada dan Imamura (2003) (Tabel 2). Berdasarkan Tabel 2 dapat diketahui bahwa semakin lebat hutan maka tingkat peredaman energi gelombang tsunami semakin besar. Contohnya gelombang tsunami setinggi 3 m yang menerjang hutan pantai dengan lebar 50 m, maka jangkauan run-up yang masuk ke daratan tinggal 81%, tinggi genangan tinggal 82%, arus setelah melewati pantai tinggal 54 % dan gaya hidrolis tsunami tinggal 39%, sehingga daerah pantai yang mempunyai vegetasi pesisir yang rapat dan tebal, mempunyai resiko kerusakan lebih kecil dibandingkan daerah

pantai yang mempunyai kerapatan dan ketebalan yang lebih rendah (Harada dan Imamura, 2003). Keberadaan terumbu karang di wilayah pesisir juga sangat penting. Selain memiliki fungsi ekologis, terumbu karang juga berfungsi sebagai pelindung pantai dari hempasan gelombang dan arus kuat yang berasal dari laut (Diposaptono dan Budiman, 2006). Tabel 2. Efektivitas peredaman tsunami oleh hutan pantai (Harada dan Imamura, 2003) Tinggi tsunami (m)

1

Hutan pantai (Shuto, 1985)

Jarak run up

Lebar hutan

Tinggi genangan

Lebar hutan

Arus

Lebar hutan

Gaya hidrolis

Lebar hutan

50 m 100 m 200 m 400 m 50 m 100 m 200 m 400 m 50 m 100 m 200 m 400 m 50 m 100 m 200 m 400 m

2

3

Mitigasi kerusakan, menghentikan benda yang hanyut dan meredam tsunami 0,98 0,86 0,81 0,83 0,80 0,71 0,79 0,71 0,64 0,78 0,65 0,57 0,86 0,86 0,82 0,76 0,74 0,66 0,46 0,55 0,50 0,11 0,18 0,71 0,58 0,54 0,57 0,47 0,44 0,56 0,39 0,34 0,31 0,24 0,53 0,48 0,39 0,33 0,32 0,17 0,01 0,13 0,08 0,02 0,01

2.2.2.5 Penggunaan lahan Kawasan pesisir termasuk dalam kerentanan yang tinggi terhadap bencana tsunami. Konsep penggunaan lahan harus melihat jarak dari garis pantai agar dapat melindungi daratan dari hantaman gelombang tsunami. Penggunaan lahan pada kawasan pesisir harus memperhitungkan antara tata guna lahan dan tingkat hunian penduduk dengan besar kecilnya resiko bahaya tsunami.

Penggunaan lahan yang berkaitan dengan berbagai aktivitas manusia seharusnya tidak dibangun pada daerah yang sangat rawan tsunami (Irfani, 2005). 2.2.2.6 Kepadatan pemukiman Sebaran dan kepadatan pemukiman menjadi salah satu faktor yang mempengaruhi resiko bencana tsunami yang akan terjadi. Pemukiman penduduk menggambarkan tingkat kepadatan penduduk dan sebaran tempat hunian yang akan mempengaruhi tingkat kerugian akibat tsunami baik dari segi kerugian jiwa maupun harta benda. Penempatan area pemukiman pada zona paling aman dari bahaya tsunami adalah prioritas utama sehingga diletakkan jauh dari laut (Irfani, 2005). Pratikno (2007) menyatakan bahwa pembangunan area pemukiman kurang padat yang efektif adalah minimal dalam radius 500 m dari garis pantai. 2.2.3 Pengurangan resiko bencana tsunami Bencana tsunami tidak dapat dicegah dan mempunyai kemungkinan kejadian berulang, sedangkan kerentanan merupakan aspek yang relatif dapat dilakukan perubahan. Oleh karena itu, untuk mengurangi resiko bencana, hanya dapat dilakukan dengan cara memperkecil kerentanan (Gambar 6).

Resiko 1 > Resiko 2

Gambar 6. Pengurangan Resiko Bencana (Bakornas Penanggulangan Bencana, 2007)

Pembuatan dan analisis tingkat resiko tsunami di suatu daerah merupakan masukan penting dalam rancangan tata ruang wilayah pesisir. Di Indonesia, pedoman resmi untuk pembuatan peta resiko tsunami belum ada. Akan tetapi, pada tanggal 14 April 2008 yang lalu, telah diadakan “Rapat Pedoman Pembuatan Peta Resiko Tsunami” yang dihadiri oleh wakil-wakil dari LIPI, BPPT, PKSPL-IPB, LAPAN, KLH, BMG, ITB, DEPKOMINFO, DEPDAGRI, Artwork dan RISTEK. Salah satu hasil penting dalam rapat tersebut adalah segera menyelesaikan “Pedoman Pembuatan Peta Resiko Tsunami” agar dapat segera didistribusikan kepada pemerintah daerah melalui Departemen Dalam Negeri (Pusat Riset Informasi Bencana Alam, 2008). 2.3 Pemetaan resiko tsunami dengan metode Cell Based Modelling 2.3.1 Penginderaan jauh sebagai penunjang riset tsunami Penginderaan jauh adalah ilmu, seni dan teknologi untuk memperoleh informasi tentang objek, daerah atau fenomena melalui analisis data yang diperoleh dengan alat tanpa kontak langsung dengan objek, daerah atau fenomena yang dikaji (Lillesand dan Kiefer, 1990). Untuk riset tsunami, citra satelit baik secara global, visual, digital dan multi temporal dapat memberikan informasi mengenai dinamika yang terjadi di daerah pesisir, baik sebelum, sewaktu maupun sesudah tsunami (Diposaptono dan Budiman, 2006). Oleh karena itu, penginderaan jauh (remote sensing) merupakan salah satu alat mutakhir yang sangat menunjang kegiatan riset tsunami, terutama jika diintegrasikan dengan SIG. Data penginderaan jauh seperti citra satelit merupakan input yang terpenting bagi SIG karena ketersediaannya secara berkala dan mencakup area yang relatif luas. Kita bisa memperoleh berbagai jenis citra satelit untuk beragam tujuan pemakaian karena

saat ini terdapat bermacam-macam satelit di ruang angkasa dengan spesifikasinya masing-masing (GIS Konsorsium, 2007). Satelit Landsat 7 ETM+ merupakan satelit milik National Atmospheric and Space Administration (NASA) Amerika yang sistem informasinya dilakukan oleh pihak swasta yaitu EOSAT (Earth Observation Satellite). Landsat 7 ETM+ (Gambar 7) diluncurkan pada tanggal 15 April 1999. Konfigurasi dasar satelit berbentuk kupu-kupu dengan tinggi ± tiga m, bergaris tengah 1.5 m dan panel matahari yang melintang kurang dari empat m (NASA, 2005).

Gambar 7. Konfigurasi Satelit Landsat 7 ETM+ (NASA, 2005) Satelit Landsat 7 ETM+ merupakan satelit dengan orbit polar sun synchronous yang memotong garis khatulistiwa ke arah selatan pada waktu lokal pukul 10.00 pagi dengan lebar sapuan 185 km dan resolusi spasial 30 x 30 m. Khusus untuk band 6 resolusinya 60 x 60 m dan band 8 adalah 15x15 m. Satelit ini berada pada ketinggian 705 km dengan periode edar 99 menit dan resolusi temporalnya 16 hari (NASA, 2005). Sistem pada Landsat 7 ETM+ yang dirancang khusus untuk mengumpulkan energi pantulan elektromagnetik

dilakukan oleh saluran 1-5, saluran 7 dan saluran 8, sedangkan untuk memancarkan energi dilakukan oleh saluran 6 (Purwadhi, 2001). Sensor Landsat 7 ETM+ mengkonversi energi pantulan matahari yang diterimanya menjadi suatu radians. Radians adalah fluks energi per satuan sudut ruang yang meninggalkan satu satuan area permukaan pada arah tertentu. Nilai radians akan dikuantifikasi menjadi nilai kecerahan (brightness values) citra yang tersimpan dalam format digital. Citra satelit Landsat ETM 7+ merupakan citra multispektral sehingga untuk mendapat ketepatan interpretasi citra perlu dipilih saluran yang paling sesuai dengan bidang kajian keadaan ekosistem pesisir (BIOTROP Training and Information Center, 2005). Karakteristik dari Landsat 7 ETM+ selengkapnya disajikan pada Tabel 3. Tabel 3. Karakteristik citra Landsat 7 ETM+ (NASA, 2005) Resolusi spektral (µm)

Resolusi spasial pada nadir (m)

Kanal

Spektrum

1

Sinar tampak violet-biru Sinar tampak hijau

0,45 – 0,52 0,52 – 0,60

30x30

2 3

Sinar tampak merah

0,63 – 0,69

30x30

4

Infra merah dekat

0,79 - 0,90

30x30

5

Infra merah menengah

1.55-1.75

30x30

6

Infra merah termal Infra merah jauh

10.40-12.50

60x60

2.08-2.35

30x30

Pankromatik

0.52-0.90

15x15

7 8 Teknologi sensor Altitude Swath widht (lebar sapuan) Resolusi temporal

30x30

Scanning mirror spectrometer 705 km 185 km 16 hari

Resolusi radiometrik Perekaman data

Best 8 of 9 bits 250 image per hari pada 31.450 km

Orbit dan inklinasi

Sun synchronus; 98,2

Peluncuran/operasi

15-April-1999/ 6 tahun

2

0

Pada penelitian ini digunakan juga data Digital Elevation Model Shuttle Radar Topographic Mission ( DEM SRTM). Shuttle Radar Topographic Mission

yang diluncurkan tanggal 11 Februari 2000 merupakan proyek internasional yang dipelopori oleh National Geospatial-Intelegence Agency (NGA) dan National Atmospheric and Space Administration (NASA). SRTM digunakan untuk menghasilkan data topografi digital dan elevasi dari 80% permukaan daratan bumi melalui sensor radar (Hajar, 2007). Permukaan bumi yang direkam oleh sensor radar tersebut berada pada 600 LU sampai 560 LS. Keunggulan sistem radar SRTM adalah dapat beroperasi siang dan malam karena menggunakan sensor aktif yang tidak tergantung oleh cahaya matahari dan tidak terpengaruh oleh cuaca dan awan. Data hasil rekaman sensor ini adalah data Digital Elevation Model atau DEM (USGS, 2007). Data DEM merupakan data digital berformat raster yang memiliki informasi koordinat posisi (x,y) dan elevasi (z) pada setiap pixel atau sel. Data DEM dapat menggambarkan kondisi topografi wilayah penelitian. Hasil pengolahan data DEM harus didukung dengan survei lapang sehingga kenampakan topografi wilayah yang direkam akurat. Akurasi data DEM sangat tergantung pada sumber titik tinggi dan resolusi spasial. Resolusi spasial DEM pada penelitian ini adalah 90 x 90 m. Semakin tinggi resolusi spasial maka semakin tinggi pula akurasi data yang diperoleh. Data DEM merupakan salah satu data penting dalam analisis daerah rawan bencana tsunami (Hajar, 2007). 2.3.2 Aplikasi Cell Based Modelling untuk kajian resiko tsunami 2.3.2.1 Sistem informasi geografis (SIG) Sistem informasi geografis (SIG) adalah kumpulan yang terorganisir dari perangkat keras, perangkat lunak, data geografi dan personal yang dirancang secara efisien untuk memperoleh, menyimpan, memperbaharui, memanipulasi, menganalisis dan menyajikan semua bentuk informasi yang berefensi geografis (Prahasta, 2002).

Sebagian besar data yang ditangani dalam SIG merupakan data spasial yaitu sebuah data yang berorientasi geografis, memiliki sistem koordinat tertentu sebagai dasar referensinya dan mempunyai dua bagian penting yang membuatnya berbeda dari data lain, yaitu informasi lokasi (spasial) dan informasi deskriptif (atribute) (GIS Konsorsium, 2007). Data spasial ini dapat dibagi menjadi dua format yaitu data raster dan data vektor. Data vektor merupakan objek geografis yang direpresentasikan ke dalam kumpulan garis, area (daerah yang dibatasi oleh garis yang berawal dan berakhir pada titik yang sama), titik dan nodes (merupakan titik perpotongan antara dua buah garis). Pada data raster, obyek geografis direpresentasikan sebagai struktur sel grid yang disebut dengan pixel (picture element). Perbandingan visualisasi antara data raster dan vektordisajikan pada Gambar 8 (GIS Konsorsium, 2007).

Data Vektor

Data Raster

Gambar 8. Perbandingan antara data raster dan data vektor Keuntungan data vektor adalah ketepatan dalam merepresentasikan fitur titik, batasan dan garis lurus, akan tetapi, tidak mampu dalam mengakomodasi perubahan gradual. Data raster sangat baik untuk merepresentasikan batasbatas yang berubah secara gradual, seperti jenis tanah dan vegetasi. Keterbatasan utama dari data raster adalah besarnya ukuran file; semakin tinggi

resolusi gridnya semakin besar pula ukuran filenya dan sangat tergantung pada kapasistas perangkat keras yang tersedia (GIS Konsorsium, 2007). 2.3.2.2 Pemodelan spasial dengan Cell Based Modelling Salah satu analisis spasial dalam SIG yang saat ini banyak digunakan untuk memodelkan keadaan di alam adalah Cell Based Modelling (ESRI, 2001). Secara umum suatu model mempresentasikan kekompleksitasan dan interaksi di alam dengan suatu penyederhanaan. Permodelan tersebut akan menolong kita untuk mengerti, menggambarkan, dan memprediksi banyak hal di alam. Ada dua model yang dikenal dalam analisis spasial, yaitu model yang merepresentasikan objek atau kenampakan di alam (representation models) dan model yang mensimulasikan proses di alam (process models). Keseluruhan model tersebut akan lebih efisien bila dilakukan pada data raster (ESRI, 2001). Analisis spasial pada data raster ini disebut metode Cell Based Modelling karena bekerja berdasarkan sel atau piksel. Operasi piksel pada Cell Based Modeling dibagi menjadi lima kelompok (ESRI, 2001) : 1. Local Function adalah operasi piksel yang hanya melibatkan satu sel. Nilai piksel output ditentukan oleh satu piksel input. 2. Focal Function adalah operasi piksel yang hanya melibatkan beberapa sel terdekat. 3. Zonal Function adalah operasi piksel yang melibatkan suatu kelompok sel yang memiliki nilai atau keterangan yang sama. 4. Global Function yang melibatkan keseluruhan sel dalam data raster dan gabungan antara keempat kelompok tersebut. 5. Application Function adalah gabungan dari keempat operasi di atas yang meliputi Local Function, Focal Function, Zonal Function dan Global Function.

Ilustrasi dari keempat operasi Cell Based Modeling dapat dilihat pada Gambar 9.

Gambar 9. Ilustrasi operasi piksel (Sumber : ESRI, 2002) Metode Cell Based Modelling memiliki beberapa keunggulan. Menurut Meaden dan Chi (1996), analisis overlay, pembuatan jarak dan pengkelasan parameter lebih mudah dilakukan secara cepat dan teratur pada data raster. Keunggulan lain yaitu struktur data raster lebih sederhana sehingga memudahkan dalam pemodelan dan analisis, kompatibel dengan data satelit serta memiliki variabilitas spasial yang tinggi dalam merepresentasikan suatu kondisi di alam. Metode ini juga memiliki beberapa kelemahan, yaitu diantaranya membutuhkan space komputer yang besar dan secara spasial memiliki tampilan yang kurang estetik karena berupa data raster yang berbentuk sel (ESRI, 2001).

2.4. Keadaan umum lokasi penelitian Kabupaten Sikka adalah salah satu kabupaten di Provinsi Nusa Tenggara Timur dengan ibu kotanya adalah Maumere. Secara geografis Kabupaten Sikka terletak antara 08022'- 08050' LS dan 121055'40"- 122041'’30" BT. Secara administrasi Kabupaten Sikka berbatasan dengan Laut Flores di sebelah utara, Kabupaten Flores Timur di timur, Laut Sawu di sebelah selatan dan Kabupaten Ende di sebelah barat. Luas wilayah Kabupaten Sikka adalah 2.253,24 km2 dengan rincian luas wilayah daratan adalah 1.731, 91 km2 yang terdiri dari daratan di Pulau Flores seluas 1.618,89 km2 dan daratan pulau-pulau kecil seluas 113,02 km2; luas laut 521.33 km2 dan panjang garis pantai adalah 379 km (Poling, 2005). Penduduk Kabupaten Sikka berjumlah sekitar 300.000 jiwa dan tersebar di 12 kecamatan. Kawasan berpenduduk padat adalah di kawasan utara yang berbatasan dengan Laut Flores, sedang kawasan selatan yang berbatasan dengan Laut Sawu atau Lautan Hindia berpenduduk jarang. Konsentrasi penduduk perkotaan ada di Kota Maumere (Kecamatan Alok) dan kawasan Geliting di Kewapante (Poling, 2005). Kondisi topografis Kabupaten Sikka secara umum merupakan daerah perbukitan. Topografis daerah didominasi oleh wilayah dengan ketinggian >500 m, yakni 42,91% dari luas wilayah daratan. Kondisi kemiringan tanah (lereng) sangat bervariasi, berkisar dari 0 s/d >40% dan didominasi oleh kemiringan tanah yang lebih besar 40% dengan luas 81,641 Ha (Poling, 2005). Pada 12 Desember 1992 Maumere dilanda gempa tsunami 7,8 SR (USGS, 2007) dan mengakibatkan kerusakan harta benda dan lebih dari 2000 penduduk meninggal dunia. Gempa tsunami tersebut disebabkan oleh aktivitas sesar busur belakang yang terletak di sisi utara Maumere, yakni Laut Flores (Poling, 2005).

3. METODE PENELITIAN 3.1 Waktu dan lokasi penelitian Lokasi penelitian terletak di Kabupaten Sikka, Propinsi Nusa Tenggara Timur yang terletak antara 08º22'-08º50' Lintang Selatan dan 121º55'40"-122º41'30º Bujur Timur (Gambar 10). Pada penelitian ini, fokus daerah penelitian adalah daerah Kabupaten Sikka yang berada pada gugus Pulau Flores (hillshade).

Gambar 10. Peta daerah penelitian Penelitian ini dimulai dengan penyusunan basis data awal pada bulan Agustus 2007, survei lapang bulan Oktober 2007 dan pengolahan citra, penyusunan basis data akhir serta pengolahan data spasial sampai bulan November 2008 yang dilakukan di Laboratorium Komputer, Departemen Ilmu dan

Teknologi Kelautan, Fakultas Perikanan dan Ilmu Kelautan, Institut Pertanian Bogor. 3.2 Alat dan bahan 3.2.1 Alat Peralatan yang digunakan pada penelitian ini adalah sebagai berikut. 1. Perangkat lunak ER Mapper 6.4 2. Perangkat lunak ArcView 3.2 3. Perangkat lunak ArcGIS 9.1 4. Perangkat lunak Global Mapperv10.00 5. Perangkat lunak Golden Software Surfer 8.0 6. Perangkat lunak Mapsource 6.3 7. Perangkat lunak Slava ITBD (Integrated Tsunami Data Base) 2006 8. GPS (Global Positionig System) Garmin eTrex v3.0 9. Kamera digital untuk dokumentasi di lapangan 3.2.2 Bahan Bahan yang digunakan dalam penelitian ini yaitu : 1. Citra Landsat 7/ETM+ Kabupaten Sikka akuisisi tanggal 27 September 2001 dan 27 Oktober 2006 path/row 112/66 BTIC BIOTROP 2. Peta batimetri Kab.Sikka 1 : 200.000 tahun 2001 dari DISHIDROS TNI AL 3. Peta rupabumi Kab.Sikka 1 : 25000 tahun 2006 dari Bakosurtanal 4. Data tinggi run up Tsunami Flores 1992 dari BPPT 5. Data seismik Kab.Sikka tahun 1904-2007 dari USGS dan BMG 6. Data kependudukan dari BPS Kab.Sikka tahun 2006 7. Data spasial Kab.Sikka dari Bappeda Sikka tahun 2006 8. Data Digital Elevation Model (DEM) SRTM dari NASA tahun 2005 yang didapat dari : http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp

3.3 Survei lapang Survei lapang yang dilakukan pada tanggal 20-26 Oktober 2007 bertujuan untuk mengetahui kondisi daerah penelitian sekaligus verifikasi pemetaan. Survei dilakukan pada 10 daerah pesisir utara Kabupaten Sikka. Posisi dan koordinat titik pengamatan selengkapnya disajikan pada Lampiran 1. Titik pengamatan tersebut merupakan daerah bekas tsunami dengan tingkat kerugian materi dan korban jiwa yang cukup besar akibat Tsunami 1992. Pengamatan yang dilakukan waktu survei lapang adalah melihat dan melakukan pengamatan terhadap kondisi pemukiman pesisir, kondisi sarana prasarana penting, ekosistem mangrove dan kondisi perairan yang dapat dilihat pada Lampiran 2. 3.4 Metode penelitian Pada penelitian ini dilakukan pengintegrasian data penginderaan jauh dan SIG. Alur penelitian ini meliputi input data (data citra, survei lapang dan data sekunder yang terkumpul), pemrosesan dan analisis.

Selanjutnya dilakukan

penyusunan basis data (baik spasial maupun non spasial) berdasarkan data yang terkumpul, lalu dibuat peta-peta tematik baik untuk parameter kerawanan dan parameter kerentanan. Keseluruhan parameter tersebut berformat raster dan kemudian dioverlay dengan menggunakan metode Weighted Overlay. Selanjutnya dilakukan proses reklasifikasi untuk menentukan tingkat kelas resiko tsunami dengan menggunakan metode Cell Based Modelling seperti terlihat pada diagram alir penelitian (Gambar 11).

Mulai

Pengumpulan data

Data Batimetri

Citra Landsat 7/ETM 2006

Data spasial Kab. Sikka

Data kegempaan

Cropping citra

Koreksi Radiometrik Koreksi Geometrik

Komposit citra

Vegetasi mangrove

Terumbu karang dan derivat

Verifikasi & editing

Konsultasi pakar

Diterima

Basis data spasial

Parameter resiko tsunami

Tidak Diterima

Pemodelan spasial Cell Based Modelling

Peta tingkat resiko tsunami

Selesai

Gambar 11. Diagam alir penelitian 3.5 Pengolahan citra satelit 3.5.1 Pengolahan citra awal (Pre processing) Langkah awal dalam pengolahan citra adalah melakukan penggabungan band dilanjutkan dengan pemotongan citra (cropping/subset) dengan tujuan untuk membatasi area penelitian. Setelah pemotongan citra, dilakukan proses pemulihan citra (image restoration) yang bertujuan untuk memulihkan data citra yang mengalami distorsi akibat proses-proses di atmosfer. Distorsi yang terjadi pada citra dapat dipulihkan dengan melakukan koreksi radiometrik dan koreksi geometrik.

Koreksi radiometrik dilakukan untuk menghilangkan faktor-faktor yang menurunkan kualitas citra. Metode koreksi radiometrik yang digunakan adalah penyesuaian histogram (histogram adjustment). Teknik ini didasarkan pada pengurangan nilai digital number sebesar bias dari masing-masing band. Nilai bias adalah nilai digital number minimum pada setiap band. Nilai bias diasumsikan sama dengan besarnya pengaruh atmosfer terhadap gelombang cahaya. Pada metode ini ditetapkan bahwa respon spektral terendah pada setiap band nilainya adalah nol, sesuai dengan bit coding sensor. Apabila nilai terendah pixel bukan nol, maka nilai penambah tersebut dianggap sebagai hamburan atmosfer. Oleh karena itu, dilakukan pengurangan nilai digital setiap pixel pada semua band sesuai dengan besarnya nilai penambah, sehingga nilai minimum pada semua band sama yaitu nol. Secara matematis, koreksi pengaruh atmosfer dengan pengaturan histogram dapat dilihat pada Persamaan 3 berikut.

DN ijk ( output

terkoreksi )

 DN ijk

( input asli )

 bias ................................. (3)

dimana: DN = Nilai digital number i = Baris j = Kolom k = Input nilai Bias = Nilai digital terendah Selanjutnya dilakukan koreksi geometrik yang bertujuan untuk memperbaiki distorsi posisi atau letak objek agar mempunyai koordinat yang benar. Distorsi ini dapat disebabkan oleh beberapa faktor diantaranya adalah variasi tinggi satelit, ketegakan dan kecepatan satelit (Lillesand dan Kiefer, 1990). Koreksi geometrik dilakukan dengan metode rektifikasi, yaitu menggunakan Ground Control Points (GCPs) dengan distribusi penyebaran titik lokasi harus merata.

Pada koreksi ini, perubahan posisi pixel juga mencakup perubahan informasi spektralnya. Untuk mengatasi hal ini maka diperlukan interpolasi nilai spektral selama transformasi geometri. Proses interpolasi ini disebut resampling. Algoritma pada koreksi ini meliputi algoritma relokasi pixel sekaligus algoritma interpolasi nilai spektral. Untuk relokasi pixel, algoritma berupa fungsi polynomial dan untuk interpolasi nilai spektral dipakai nearest neighbour karena tidak mengubah nilai spektral. Setelah rektifikasi akan diperoleh posisi citra yang sesuai dengan koordinat sebenarnya di bumi (real world coordinate system). Citra Landsat yang digunakan pada penelitian ini sudah dikoreksi baik secara geometri maupun radiometrik. Akan tetapi, untuk keakuratan data dilakukan koreksi ulang terhadap citra tersebut. 3.5.2 Penajaman citra (Image Enhancement) Proses penajaman citra dapat didefinisikan sebagai pemilihan manipulasi kontras kenampakan suatu citra sehingga informasi tersebut dapat lebih mudah diinterpretasikan untuk suatu tujuan tertentu. Penajaman citra bertujuan untuk meningkatkan informasi yang ada dalam citra dengan meningkatkan kekontrasan citra. Penajaman citra dilakukan dengan menggunakan komposit warna (colour composite) dan algoritma yang sesuai dengan objek yang ingin diperjelas informasinya. Pada penelitian ini penajaman citra dilakukan untuk mengetahui kondisi ekosistem pesisir. 3.5.2.1 Penajaman citra untuk terumbu karang Terumbu karang merupakan ekosistem pesisir yang penting sebagai peredam energi gelombang tsunami. Pemetaan terumbu karang (Gambar 12) dilakukan dengan pendekatan algoritma Lyzenga (1978) yang dikenal dengan ”standart exponential attenuation model” yang sering digunakan dalam pemetaan-pemetaan terumbu karang (Engel, 1988 in Siregar et al., 1995).

Algoritma Lyzenga menggunakan band 1 dan band 2 karena kedua band ini diasumsikan memiliki penetrasi yang baik ke kolom air. Persamaan algoritma Lyzenga tersebut yaitu :

Y  ln TM 1  ki / kj * ln TM 2 .................................................. (4) dimana : Y = citra hasil ekstraksi TM1 = band 1 Landsat 7/ETM+ TM2 = band 2 Landsat 7/ETM+ ki/kj = koofisien atenuasi, yang diperoleh dari : = a  a2  1

dengan, a 

VarTM 1  VarTM 2 ......…………………. (5) 2 * Co var TM 1 TM 2

Citra Landsat 7/ETM+ 2006 (corrected and cropped)

   

Komposit 421

Komposit 543

Training area

Training area

Hitung nilai varian band 1 dan band 2 Hitung covar band 1/2 Hitung nilai a Hitung nilai ki/kj

Penajaman dengan algoritma Lyzenga : Ln TM1+ki/kj ln TM2

Citra klasifikasi supervised

Vegetasi mangrove

Peta tematik vegetasi mangrove

Citra hasil penerapan algoritma

Klasifikasi

Peta tematik terumbu karang dan derivatnya

Gambar 12. Diagram alir pengolahan citra untuk pemetaan ekosistem pesisir

3.5.2.2 Penajaman citra untuk ekosistem mangrove Pada penelitian ini dilakukan juga pemetaan ekosistem mangrove, karena mangove merupakan vegetasi pantai yang berperan penting dalam meredam energi gelombang tsunami. Pada komposit citra 453, vegetasi mangrove akan terlihat berwarna merah. Hal ini karena klorofil dalam daun mangrove menyerap dengan kuat sinar merah dan memantulkan kuat sinar infra merah (Earth Observatory NASA, 2007). Setelah dilakukan training area, terbentuk kelas mangrove. Berdasarkan pengolahan data dapat diketahui sebaran dan juga luasan mangrove di Kabupaten Sikka. Selain data terumbu karang dan mangrove, secara visual pada citra dengan komposit citra 321 bisa didapatkan juga data keberadaan pulau-pulau kecil penghalang. Ekosistem pesisir dan pulau penghalang merupakan data pendukung penting dalam mengurangi tingkat resiko tsunami di suatu daerah. 3.6 Klasifikasi citra Citra yang telah ditransformasikan dengan algoritma yang sesuai kemudian diklasifikasi. Klasifikasi merupakan suatu proses pengelompokkan nilai reflektansi dari setiap objek ke dalam kelas-kelas tertentu sehingga mudah dikenali. Pada penelitian ini klasifikasi yang digunakan untuk terumbu karang dan mangrove adalah klasifikasi supervised dengan metode maximum likehood karena metode ini mampu menganalisis fungsi peluang multidimensional untuk menentukan suatu pixel tertentu lebih berpeluang masuk ke dalam suatu kelas tertentu (Pasek, 2007). Tahap ekstraksi keseluruhan data citra dapat dilihat secara lengkap pada Gambar 12.

3.7 Metode penentuan tingkat resiko tsunami 3.7.1 Penentuan kerawanan tsunami Penentuan tingkat kerawanan tsunami dilakukan dengan tujuan mengetahui daerah-daerah yang berpotensi terhadap bencana tsunami. Penentuan tingkat kerawanan dilakukan dengan dua tahap yaitu : (1) membuat peta seismisitas dan (2) membuat peta rawan tsunami. 3.7.1.1 Pemetaan Seismisitas Peta seismisitas adalah peta yang menggambarkan tectonic setting plate dan sebaran titik-titik gempa di suatu wilayah. Peta ini merupakan peta pendukung penting untuk melihat sebaran titik gempa dan tsunami yang pernah terjadi di suatu wilayah, dalam hal ini adalah Kabupaten Sikka. Pada penelitian ini, data yang dipakai untuk membuat peta seismisitas adalah data kegempaan dari BMG selama 107 tahun terakhir. Berdasarkan kumpulan data tersebut, kemudian dipetakan titik-titik gempa dan plate tectonic setting sehingga dapat diketahui dan dianalisis sebaran pusat-pusat gempa bumi dan tsunami yang pernah terjadi di Kabupaten Sikka. Seismisitas merupakan parameter penting dalam pemetaan resiko tsunami karena merupakan preliminary analysis daerah-daerah yang berpotensi terkena tsunami serta untuk mengetahui hubungan antara tektonik setting, kegempaan dan kejadian tsunami. 3.7.1.2 Pemetaan daerah rawan tsunami Untuk analisis tingkat kerawanan tsunami harus dilakukan pemetaan daerah rawan tsunami yaitu dengan memetakan data tinggi run up tsunami. Peta run up merupakan parameter penting dalam menentukan tingkat resiko tsunami. Peta run up yang ideal adalah peta run up yang ditunjang dengan data lapang yang lengkap dan data hasil pemodelan. Akan tetapi, pada penelitian ini data run up yang digunakan adalah hanya dari data survei lapang tsunami Flores tahun 1992

dari BPPT (Lampiran 3). Data tersebut juga sudah dibandingkan dengan data run up dari BMG dan DKP. Data run up yang diperoleh adalah data titik ketinggian tsunami. Untuk mendapatkan data spasial run up, maka harus dibuat poligon daerah run up, yang nantinya dapat dikelaskan dan dioverlay untuk mendapatkan peta resiko tsunami. Adapun pengolahan data run up dilakukan secara sebagai berikut. -

Buat data spasial run up yaitu berupa point ketinggian run up tsunami.

-

Buat poligon area dengan ketinggian topografi 10 m dari permukaan laut (sesuai dengan tinggi run up maksimum tsunami di Sikka, yaitu 10 m)

-

Display dan autolabel data titik run up pada poligon area tersebut, sehingga titik-titik run up akan menempel di sepanjang garis pantai daerah penelitian

-

Antara dua titik run up tarik garis sejajar yang menghubungkan kedua titik tersebut dan tentukan batas tengah dari kedua titik run up tersebut.

-

Selanjutnya, tarik garis tegak lurus ke arah darat sehingga terbentuk poligon daerah run up tsunami

-

Untuk mendapatkan data run up yang lebih baik, pada waktu membuat poligon run up, dioverlay dengan data DEM sehingga bisa dilihat daerah yang flat dan daerah yang tinggi. Ketinggian data di DEM dapat dilihat dengan menggunakan tools identify.

-

Poligon run up tersebut di atas kemudian diberi skor dan dikonvert ke raster dan di reklasifikasi sesuai matriks. Sehingga akhirnya dapat diperoleh data spasial run up tsunami.

Menurut Iida (1963) hubungan tingkat kerusakan dengan tinggi run up disajikan pada Tabel 4.

Tabel 4. Hubungan tingkat kerusakan dengan tinggi run up tsunami No. 1 2 3 4 5

Tinggi run up (m) >16 6-16 2-6 0.75-2 <0.75

Daya rusak

Skala

Sangat besar Besar Menengah Kecil Sangat kecil

5 4 3 2 1

Klasifikasi di atas adalah acuan yang digunakan untuk mengkelaskan tingkat kerawanan tsunami. Pada pembuatan peta run up di atas memang tidak menggambarkan secara detail, kondisi run up yang sebenarnya. Akan tetapi, setelah dioverlay dengan algoritma, baru menunjukkan suatu peta resiko tsunami yang komprehensif. 3.7.2 Penentuan tingkat kerentanan terhadap tsunami Analisis kerentanan yang dikaji pada penelitian ini adalah kerentanan lingkungan, infrastruktur dan kerentanan sosial kependudukan. Akan tetapi, parameter yang dikelaskan dalam matriks, hanya parameter kerentanan lingkungan karena dapat dispasialkan. Sedangkan parameter kerentanan infrastruktur dan sosial kependudukan, tidak dikelaskan dalam matriks tetapi akan dioverlay dengan peta resiko tsunami untuk mendapatkan peta resiko kependudukan terhadap tsunami dan peta resiko infrastruktur terhadap tsunami. Parameter yang digunakan untuk menentukan tingkat kerentanan lingkungan terhadap tsunami adalah elevasi daratan, slope (kemiringan), morfometri pantai, penggunaan lahan, jarak dari garis pantai dan jarak dari sungai. Dasar pengambilan parameter tersebut ditentukan berdasarkan penelitian sebelumnya dengan melihat parameter penentu tingkat kerentanan di suatu wilayah yang kemudian dimodifikasi sesuai dengan konsultasi pakar dan pembimbing berdasarkan kondisi daerah penelitian. Beberapa hasil penelitian terdahulu yang dijadikan acuan yaitu penelitian yang pernah dilakukan oleh

Bakosurtanal dan Pusat Studi Bencana Alam UGM (PSBA UGM) (2002); Papathoma et al., (2003), Diposaptono dan Budiman (2005), dan Agung (2006). Setiap parameter memilki kontribusi yang berbeda terhadap tingkat kerentanan dan resiko bencana tsunami. Oleh karena itu, dalam penelitian ini dilakukan sistem penentuan bobot dan skor untuk masing-masing parameter sesuai dengan besarnya pengaruh parameter tersebut. Parameter yang dianggap memiliki pengaruh paling besar terhadap resiko tsunami, diberikan bobot paling besar. Jadi, semakin besar bobot parameter pemicu resiko bencana tsunami, semakin besar kontribusinya terhadap resiko tsunami dan sebaliknya. Penjelasan masing-masing parameter adalah sebagai berikut. 1. Elevasi daratan Elevasi merupakan salah satu faktor penting yang mempengaruhi resiko tsunami yang akan terjadi. Daerah dengan elevasi rendah, akan mudah dihantam gelombang tsunami, dan sebaliknya. Kelas ketinggian daratan menurut Bappeda Kab.Sikka (2003) adalah 0-25 m, 25-100 m, 100-500 m, 500-1000 m, dan > 1000 m. Akan tetapi, dengan melihat data run up maksimum yang terjadi waktu tsunami Flores adalah 26.2 m dan untuk Kabupatan Sikka adalah 10 m serta untuk mempermudah dalam pengolahan dan analisis data raster (karena cakupan daerah yang luas), maka dalam penelitian ini dilakukan pengkelasan ulang elevasi sebagai berikut : ≤10 m; 10-25 m; 25-50 m; 50-100 m dan >100 m. 2. Kemiringan Kemiringan daratan akan mempengaruhi tinggi run up tsunami di wilayah pesisir. Semakin besar kemiringan daratan (curam), maka tinggi run up semakin rendah, dan sebaliknya. Pengkelasan serta pembobotan kemiringan daratan dalam penelitian ini mengacu pada pembagian

kemiringan wilayah oleh Bappeda Kabupaten Sikka yang kemudian dimodifikasi yaitu 10%; >10-20%; >20-30%; >30-40%; dan >40%. 3. Morfometri pantai Morfologi pantai sangat berpengaruh besar terhadap tingkat energi tsunami yang akan terhempas ke daratan. Tipe pantai teluk akan mengalami amplifikasi energi gelombang dan tipe tanjung akan mereduksi energi gelombang. Pada penelitian ini digunakan lima kelas tipe morfometri pantai hasil modifikasi berdasarkan penelitian dari Bakosurtanal-PSBA UGM (2002), konsultasi pakar dan pembimbing yaitu tipe teluk V, tipe teluk U, tipe tanjung, tipe pantai lurus dan non teluk atau tanjung. 4. Penggunaan lahan (landuse) Penggunaan lahan merupakan salah satu parameter yang mempengaruhi tingkat resiko tsunami. Tsunami dapat menyebabkan perubahan tata guna lahan, contohnya Tsunami Aceh 2004 telah menyebabkan beberapa wilayah pemukiman maupun lahan potensial, sekarang berubah menjadi lahan kosong. Oleh karena itu, perlu penataan ruang dengan baik dalam rangka mengurangi resiko tsunami. Acuan penggunaan lahan dalam penelitian ini dibagi berdasarkan klasifikasi Bakosurtanal-PBSA UGM (2002) dan Diposaptono dan Budiman (2006) serta berdasarkan konsultasi pakar. Pengkelasan untuk landuse dapat dilihat pada matriks resiko tsunami. 5. Jarak dari garis pantai Tsunami merupakan fenomena alam yang bersifat merusak, sehingga perlu memperhatikan adanya kawasan penyangga (buffer zone). Oleh karena itu, pembangunan kawasan untuk pemukiman dan pusat-pusat kegiatan penting lainnya harus memperhatikan jarak dari garis pantai guna mengurangi resiko tsunami. Acuan dasar untuk pembuatan jarak (buffer) merujuk pada UU RI No.27 tahun 2007 tentang pengelolaan wilayah pesisir dan pulau-pulau kecil

yaitu sempadan pantai adalah daratan sepanjang tepian yang lebarnya proporsional dengan bentuk dan kondisi fisik pantai, minimal 100 (seratus) meter dari titik pasang tertinggi ke arah darat. Pada penelitian ini jarak dari garis pantai diklasifikasi menjadi lima kelas yaitu <200 m; 200-500 m; 5001000; 1000-1500 m dan >1500 m. 6. Jarak dari sungai Jarak dari sungai merupakan parameter yang mempengaruhi tingkat resiko tsunami. Tsunami yang memasuki kanal banjir/sungai akan mengakibatkan kerusakan yang lebih besar karena adanya pemusatan energi tsunami sehingga semakin mendorong tsunami masuk lebih jauh ke daratan. Contohnya Tsunami California yang melewati kanal-kanal pengendali banjir dapat masuk ke daratan sampai 1 mil (1,802 km). Oleh karena itu, perlu dilakukan buffer dari sungai. Pada penelitian ini buffer dari sungai dilakukan pada jarak <100 m; 100-200 m; 200-300 m; 300-500 m; dan >500 m. 3.8 Pengolahan data dengan cell based modelling Cell based modelling (selanjutnya disingkat CBM) merupakan metode yang digunakan dalam pengolahan data pada penelitian ini. Oleh karena itu, berbagai parameter kerentanan tersebut di atas harus dikonversi ke bentuk raster. Untuk parameter topografi, input layer data berupa data line contour yang selanjutnya diolah menjadi data set TIN (Triangulated Irregular Network). Data set TIN adalah suatu struktur data yang digunakan untuk model permukaan seperti ketinggian yang berbentuk jejaring triangular (ESRI, 2002). Besar kecilnya ukuran segitiga pada data TIN tergantung dari banyaknya data. Semakin banyak data, semakin rapat/kecil ukuran segitiga yang menunjukkan interpolasi yang detail. Hasil TIN berupa topografi yang smooth dengan kenampakan informasi yang cukup detail. Selanjutnya data TIN tersebut

dikonversi ke raster dengan fungsi convert feature to raster pada menu 3D Analyst. Selanjutnya dilakukan pengkelasan sesuai matriks untuk mendapatkan kelas elevasi daratan. Untuk data kemiringan langsung diturunkan dari data elevasi dengan fungsi Slope Surface Analyst pada menu Spatial Analyst sehingga sudah berformat raster. Selanjutnya tinggal diklasifikasi berdasarkan matriks. Untuk morfologi pantai, penulis tidak mendapatkan data spasial digital morfologi pantai. Oleh karena itu, input data morfologi diambil dengan cara melakukan on screen digitize pada peta batimetri Dishidros TNI AL skala 1:200.000. Sebelum dilakukan proses digitasi, terlebih dahulu dilakukan konversi peta analog ke format digital melaui proses scanning, registrasi peta (menggunakan software Global Mapper 10), transformasi koordinat dan kemudian baru dilakukan proses digitasi. Pada pembuatan peta morfologi, dilakukan buffering area satu (1) km dari garis pantai. Hal ini dilakukan dengan pertimbangan bahwa batas limpasan tsunami ke daratan tidak sampai sejauh 1 km. Digitasi dilakukan dengan teknik visual pada kenampakan daerah pantai bertipe pantai teluk V, teluk U, tanjung, pantai lurus dan non teluk/tanjung. Selanjutnya data tersebut dikonversi ke format raster dan dikelaskan berdasarkan kriteria yang ditetapkan. Input data penggunaan lahan sudah dalam bentuk digital spasial. Untuk mendapatkan data yang benar, dilakukan proses re-editing data landuse (mengacu juga pada kenampakan pada citra satelit). Tahap selanjutnya adalah mengkonversi data ke dalam bentuk raster dan dilakukan pengkelasan. Pada pembuatan buffer zone untuk jarak dari garis pantai dan sungai, dilakukan buffering dan reclassify. Selanjutnya dikonversi ke format data raster dan dikelaskan berdasarkan matriks.

Selain pengolahan di atas juga dilakukan pemetaan batimetri dari peta batimetri dengan cara digitize on screen dan outputnya berupa point kedalaman. Selanjutnya diinterpolasi dengan metode Natural Neighbours. Metode ini merupakan metode interpolasi yang paling efektif jika sel inputnya cukup banyak, sehingga akan dihasilkan peta batimetri yang mirip dengan kondisi real di alam dan hasilnya juga lebih smooth. Pada pengolahan data raster di ArcGIS 9.1, seluruh output cell size adalah 30x30 m. Pemilihan resolusi spasial tersebut mengikuti resolusi spasial citra Landsat7/ETM+ supaya mudah dalam overlay raster. Gambaran desain fisik penentuan tingkat resiko tsunami selengkapnya disajikan pada Gambar 13. Elevasi

Slope Kerentanan lingkungan

Morfometri

Analisis seismisita s

Landuse Parameter Kerentanan

Jarak dari garis pantai

Parameter Kerawanan

Jarak dari sungai Overlay

Pemukiman & infrastruktur

Kerentanan sosial kependudukan & infrastruktur

Peta Resiko Tsunami Kab. Sikka

Analisis resiko tsunami

Gambar 13. Desain penentuan tingkat resiko tsunami

Daerah rawan tsunami

Analisis spasial yang digunakan untuk penentuan resiko tsunami menggunakan metode CBM, baik untuk pengkelasan maupun untuk overlay setiap parameter yang diperoleh. Seluruh parameter yang digunakan akan berformat grid yang terdiri atas sekumpulan sel. Setiap sel memiliki nilai tertentu yang besarnya tergantung dari besarnya nilai masing-masing parameter. Sel-sel tersebut kemudian dikelompokkan berdasarkan nilainya ke dalam lima kelas (zona), yaitu kelas resiko sangat tinggi, resiko tinggi, resiko sedang, resiko rendah dan resiko sangat rendah. Pengelompokkan sel dalam data raster ini mengikuti operasi zonal function karena setiap sel akan dikodekan berdasarkan kriteria yang membentuk suatu zona. Setiap zona akan memiliki kisaran nilai parameter sebagaimana yang terdapat dalam Tabel 5. Selanjutnya dilakukan proses overlay seluruh parameter penentu resiko tsunami dengan sistem pembobotan dan skoring. Tabel 5. Matriks Resiko Tsunami No

Parameter

Bobot* (%)

Resiko Sangat Tinggi

Skor

Resiko Tinggi

1

Run up tsunami (m)

20

>16

5

>6-16

Skor

Resiko Sedang

Skor

Resiko Rendah

Skor

Resiko Sangat Rendah

Skor

>2-6

3

>0.75-2

2

<0.75

1

Kerawanan 4 Kerentanan 2

Elevasi (m)

15

<10

5

>10-25

4

>25-50

3

>50-100

2

>100

1

3

Slope (%)

15

<10

5

10-20

4

>20-30

3

>30-40

2

>40

1

4

Morfometri pantai Landuse

15

Teluk V

5

Teluk U

4

Tanjung

3

Lurus

2

1

10

Pemukiman, sawah, mangrove, hutan rawa

5

Kebun

4

Ladang/ tegalan

3

2

15

0-200

5

>200-500

4

>5001000

3

Semak belukar, rumput/ tanah kosong >1000-1500

2

Non teluk/ tanjung Hutan, batuan cadas dan gamping >1500

10

0-100

5

>100-200

4

>200-300

3

>300-500

2

>500

1

5

6

Jarak dari garis pantai (m) 7 Jarak dari sungai (m) Jumlah bobot x skor

100

5

4

3

2

(Sumber : UU RI No.27 Tahun 2007; Bappeda Kab. Sikka 2006; Bakosurtanal-Pusat Studi Bencana Alam UGM 2002; Diposaptono dan Budiman 2006; *modifikasi berdasar konsultasi pakar dan pembimbing)

3.9 Matriks tingkat resiko tsunami Kabupaten Sikka Pembuatan matriks resiko ini dimulai dengan menentukan parameter apa saja yang berpengaruh terhadap tingkat kerawanan dan kerentanan wilayah dan

1

1

1

penduduk Sikka terhadap tsunami. Penentuan kriteria pun harus disesuaikan dengan kondisi lingkungan di Kabupaten Sikka. Kriteria yang digunakan dalam penelitian ini berasal dari beberapa sumber pustaka hasil penelitian terdahulu, hasil konsultasi dengan pembimbing dan beberapa pakar tsunami-kegempaan. Parameter yang digunakan dalam penentuan resiko tsunami adalah tinggi run up tsunami, topografi, morfologi pantai, penggunaan lahan, jarak dari garis pantai dan jarak dari sungai. Berdasarkan parameter tersebut, maka dapat disusun matriks untuk menentukan tingkat resiko tsunami seperti yang disajikan pada Tabel 5. Pembuatan matriks ditentukan dengan pembobotan dan skor. Pemberian skor dimaksudkan untuk menilai faktor pembatas pada setiap parameter, sedangkan pembobotan setiap parameter berdasarkan pada dominannya pengaruh parameter tersebut dalam penentuan tingkat resiko tsunami. Komponen pemicu yang dianggap memiliki pengaruh yang paling besar terhadap tsunami, diberikan bobot dan skor yang paling besar, dan sebaliknya. Pemberian bobot pada penelitian ini adalah berkisar dari 10-20% dan pemberian skor dalam kisaran 1-5 yang menunjukkan kelas tingkat resiko tsunami. Nilai tiap kelas didasarkan pada perhitungan rumus berikut (Pasek, 2007) :

N   Bi x Si …………………………………………. (6) dimana : N = Total bobot nilai Bi = Bobot pada tiap kriteria Si = Skor pada tiap kriteria i = Parameter Selang tiap-tiap kelas diperoleh dari jumlah perkalian nilai maksimum dari tiap bobot dan skor dikurangi jumlah perkalian nilai minimumnya yang kemudian

dibagi dengan jumlah parameter yang digunakan. Secara matematis selang kelas tingkat resiko dirumuskan sebagai berikut (Pasek, 2007) :

Lebar selang kelas ( L) 

 (B S ) i

i

max

  ( Bi S i ) min n

....................... (7)

dimana : L = Lebar selang kelas n = Jumlah kelas Berdasarkan perhitungan di atas akan diperoleh nilai lebar selang kelas tingkat resiko tsunami sebesar 0.800 dengan nilai Nminimum sebesar satu (1) dan nilai Nmaksimum sebesar lima (5). Kelas resiko sangat rendah (R1) didapatkan dari 1 ditambah dengan 0,800. Nilai kelas resiko rendah (R2) didapat dari selang kelas maksimum R1 yaitu 1,800 ditambah 0,800. Nilai selang kelas resiko sedang (R3) didapatkan dari selang maksimum R2 yaitu 2,600 ditambah dengan 0,800. Nilai selang kelas resiko tinggi (R4) didapatkan dari selang maksimum R3 yaitu 3,400 ditambah dengan 0,800. Nilai selang kelas sangat tinggi (R5) didapatkan dari selang maksimum R4 yaitu 4,200 ditambah dengan 0,800. Secara singkat selang kelas masing-masing kelas resiko dapat ditetapkan sebagai berikut. Kelas R5 (Resiko Sangat Tinggi)

= jika 4,201  N  5,000

Kelas R4 (Resiko Tinggi )

= jika 3,401 < N  4,200

Kelas R3 (Resiko Sedang)

= jika 2,601  N  3,400

Kelas R2 (Resiko Rendah)

= jika 1,801  N  2,600

Kelas R1 (Resiko Sangat Rendah)

= jika 1,000  N  1,800

Pada waktu pengolahan data dengan ArcGIS 9.1, secara otomatis nilai tersebut akan dideskripsikan berupa model tingkat resiko tsunami wilayah

Kabupaten Sikka dengan lima kelas resiko tsunami. Berdasarkan peta resiko tsunami tersebut kemudian dibuat juga peta resiko tsunami per kecamatan. Untuk melihat tingkat resiko sosial kependudukan dan infrastruktur terhadap tsunami, maka peta resiko dioverlay dengan data sosial kependudukan dalam hal ini adalah data pemukiman dan untuk melihat tingkat resiko infrastruktur terhadap tsunami, dioverlay data infrastruktur dengan resiko tsunami. Pada penelitian ini dilakukan juga pemodelan peta resiko tsunami Kabupaten dengan run up tsunami 20 meter sehingga dapat diketahui perbedaaan luasan dengan pemodelan sebelumnya dan sebagai salah satu upaya antisipasi dalam rangka mitigasi bencana tsunami di wilayah Kabupaten Sikka.

4. HASIL DAN PEMBAHASAN 4.1 Pengolahan data spasial dari citra satelit 4.1.1 Pengolahan citra awal Citra satelit yang digunakan adalah Citra Landsat 7/ETM+ path/row 122/066 akuisisi 27 Oktober 2006. Penulis menggunakan data citra tersebut karena selain harganya terjangkau, resolusi spasialnya cukup tinggi yaitu 30x30 m, sehingga bisa mewakili pemetaan resiko tsunami untuk skala satu kabupaten. Perbandingan citra satelit sebelum dan sesudah tsunami 1992 tidak dapat dilakukan karena tidak tersedianya data citra pra tsunami 1992. Data citra yang diperoleh dalam penelitian ini sudah terkoreksi baik secara radiometrik maupun geometrik. Akan tetapi, untuk mendapatkan data dan hasil yang lebih akurat, dilakukan pengkoreksian ulang secara radiometrik dan geometrik. Koreksi radiometrik citra dilakukan untuk mengurangi pengaruh hamburan atmosfer yang mempengaruhi nilai spektral citra. Metode yang digunakan dalam koreksi ini adalah metode penyesuaian histogram (histogram adjustment). Metode ini dipilih karena dapat dilakukan dengan mudah, tidak menggunakan algoritma yang rumit serta memberikan hasil yang baik (Susilo dan Lumban Gaol, 2008). Pada metode ini nilai respon terendah dari setiap band adalah nol, sehingga dilakukan pengurangan nilai digital setiap piksel pada semua band, sehingga nilai minimumnya sama yaitu nol. Pada waktu pengolahan data citra 2006, keseluruhan band sudah memiliki nilai spektral minimum nol dan nilai spektral maksimum 255, sehingga pengolahan langsung dilanjutkan dengan koreksi geometrik.

Koreksi geometrik citra dilakukan dengan metode rektifikasi (image to image rectification) dengan acuan citra akuisisi tanggal 27 September 2001 path/row yang sama. Metode rektifikasinya adalah polynomial dan tipe resampling citra adalah nearest neighbour dengan mengambil beberapa titik kontrol bumi (Ground Control Point). GCP yang diambil adalah sebanyak 30 titik yang menyebar merata pada permukaan objek daratan pada citra. Kualitas pada GCP yang digunakan cukup baik dengan nilai root mean square error (RMSE) yang kecil yaitu 0.3 (Lampiran 4). Ini berarti simpangan kesalahan yang dapat ditoleransi adalah 9 m (didapatkan dari perkalian RMSE dengan piksel Citra Landsat). RMSE merupakan sebuah standar pengukuran statistika yang menjelaskan perbedaan antara lokasi titik yang sebenarnya di alam (actual point location) dan lokasi titik yang diestimasi secara matematik (mathematically estimated point location) (Earth Resource Mapping Ltd., 2008). Nilai RMSE sebaiknya kurang dari satu. Hal ini berarti bahwa rata-rata kesalahan di titik lokasi (X, Y) kurang dari satu sel citra. 4.1.2 Ekstraksi data citra Data yang diekstrak dari Citra Landsat 7/ETM+ adalah data ekosistem pesisir yang meliputi ekosistem terumbu karang, lamun, dan mangrove. Ekosistem pesisir merupakan ekosistem penting yang berperan sebagai peredam energi gelombang tsunami, sehingga bisa meminimalkan dampak resiko tsunami. Pada penelitian ini, ekosistem pesisir termasuk ke dalam parameter pendukung sehingga tidak dikelaskan dalam matriks resiko tsunami. Pendugaan awal ekosistem pesisir dilakukan dengan penajaman citra. Penajaman citra untuk terumbu karang adalah komposit 421 dan untuk mangrove adalah komposit 453. Penajaman tersebut akan memberikan

gambaran yang interpretatif tentang kenampakan atau keberadaan ekosistem di wilayah pesisir. Untuk mendapatkan data ekosistem mangrove, setelah dilakukan training area, akan didapatkan data kelas mangrove. Berdasarkan pengolahan data citra, diketahui bahwa ekosistem mangrove umumnya terdapat di bagian pesisir utara Kabupaten Sikka, sedangkan di bagian selatan tidak ditemukan adanya ekosistem mangrove. Hal ini diduga karena wilayah pesisir selatan Sikka merupakan daerah bertebing dan berarus cukup besar, sehingga mangrove sulit tumbuh dan berkembang biak. Mangrove merupakan vegetasi pesisir pantai yang hidup pada daerah yang terlindung seperti muara sungai. Sebaran mangrove tidak ditemukan pada semua wilayah kecamatan di Kabupaten Sikka. Mangrove dominan terdapat di sekitar pulau-pulau kecil, terutama di Pulau Besar. Selain itu, terdapat juga di Kecamatan Alok dan Talibura dan serta sedikit terdapat pada Kecamatan Magepanda dan Waigete. Luasan vegetasi mangrove yang diekstrak dari citra adalah 249,30 Ha. Berdasarkan laporan dari Dinas Kehutanan dan Perkebunan Kabupaten Sikka diketahui bahwa luas vegetasi mangrove tahun 2007 adalah 219,74 Ha (dengan perincian mangrove yang baik adalah 73,88 Ha dan mangrove yang rusak adalah 145,76 Ha). Hasil wawancara dengan penduduk pesisir, diketahui bahwa semakin banyak hutan mangrove yang rusak di Kabupaten Sikka. Hal ini disebabkan karena adanya kegiatan destruktif masyarakat yang mengganggu keberadaan ekosistem mangrove, seperti pembabatan liar. Ekosistem mangrove yang dominan di Kabupaten Sikka adalah jenis Rhizopora (Dinas Kehutanan dan Perkebunan Kabupaten Sikka, 2007). Rhizopora merupakan jenis mangrove yang paling mampu meredam abrasi dan tinggi gelombang tsunami dibandingkan jenis lainnya. Hal ini disebabkan karena

Rhizopora memiliki formasi akar, batang, ranting, dan daun yang lebih rapat (Whitten et al., 1987 in Yuwono, 2005). Selain ekosistem mangrove, hampir di keseluruhan daerah pesisir Kabupaten Sikka terdapat vegetasi pantai lain yang berperan penting dalam meredam energi gelombang tsunami yaitu pohon kelapa, enau, dan lontar. Selanjutnya untuk mendapatkan data ekosistem terumbu karang dilakukan dengan menggunakan Transformasi Lyzenga. Setelah penulis melakukan training area dan mengekstrak nilai digital Band 1 dan Band 2 didapat nilai koofisien atenuasi perairan (ki/kj) sebesar 1,1527 (Lampiran 5) sehingga algoritma Lyzenga yang digunakan adalah Y= ln (TM1) + 1,1527 + ln (TM2). Selanjutnya nilai spektral citra hasil transformasi tersebut diklasifikasikan dengan metode supervised. Hasil pengolahan yang didapat adalah data ekosistem terumbu karang, karang mati, lamun dan pasir. Terumbu karang jarang hidup pada kedalaman sampai 40-60 m. Oleh karena itu, pada pemetaan terumbu karang, data terumbu karang yang diambil hanya sampai pada kedalaman 25 m saja. Berdasarkan pemetaan ekosistem terumbu karang dan derivatnya (Gambar 14) terlihat bahwa terumbu karang hidup hampir menyebar di sepanjang utara Kabupaten Sikka, kecuali di daerah Kecamatan Kewapante. Sedangkan karang mati banyak tersebar di sepanjang pantai selatan Kabupaten Sikka dan terlihat juga pada pulau-pulau kecil. Hal ini diduga karena pengaruh gelombang yang cukup besar di daerah pantai selatan. Ekosistem lamun tersebar pada semua wilayah tapi dalam jumlah yang sangat sedikit.

Gambar 14. Peta ekosistem pesisir Kabupaten Sikka Ekosistem mangrove, lamun dan terumbu karang yang sudah diekstrak tersebut kemudian digabung menjadi satu data spasial ekosistem pesisir Kabupaten Sikka dengan luasan seperti disajikan pada Tabel 6. Tabel 6. Luasan ekosistem pesisir Kabupaten Sikka dari citra satelit No Ekosistem

Luasan (Ha)

1

Mangrove

249,3

2

Terumbu karang 1.125,18

3

Karang mati

2.453,4

4

Lamun

489,24

Berdasarkan Tabel 6 dapat dilihat bahwa ekosistem terumbu karang merupakan ekosistem yang paling luas, akan tetapi didominasi oleh karang mati. Hal ini diduga karena pengaruh faktor lingkungan dan aktivitas masyarakat yang

menangkap ikan dengan cara-cara yang tidak ramah lingkungan, sehingga banyak karang yang mati. Resiko tsunami umumnya paling parah menimpa daerah pesisir. Keberadaan ekosistem pesisir merupakan ekosistem penting sebagai penghalang pertama untuk meredam energi gelombang tsunami. 4.2 Analisis kerawanan tsunami di Kabupaten Sikka 4.2.1 Seismisitas di Kabupaten Sikka Tsunami merupakan salah satu fenomena bencana alam yang tidak dapat dihindari ataupun dicegah dan memiliki kemungkinan untuk terjadi lagi, sehingga yang dapat dilakukan adalah mengurangi kemungkinan resikonya. Langkah awal penentuan resiko tsunami adalah menentukan tingkat kerawanan suatu wilayah terhadap tsunami. Kabupaten Sikka merupakan wilayah yang berpotensi tinggi dilanda tsunami. Berdasarkan Gambar 15 dapat dilihat bahwa Kabupaten Sikka yang berada di tengah Pulau Flores, diapit oleh Florest Thrust di utara, Wetar Thrust di bagian Timur, Banda Arc di bagian tenggara, dan Sunda Arc di bagian barat daya yang merupakan daerah patahan lempeng. Daerah tersebut adalah zona seismik aktif yang merupakan sumber gempa dan tsunami. Hal ini menyebabkan Kabupaten Sikka sering mengalami gempa bumi, baik gempa dalam, menengah maupun gempa dangkal (yang dapat menyebabkan tsunami). Adanya sesarsesar di daratan juga menyebabkan gempa-gempa yang bersifat merusak. Besarnya tingkat kegempaan dapat diketahui dengan melihat penyebaran pusatpusat gempa bumi seperti pada Gambar 15 berikut.

Gambar 15. Peta seismisitas Kabupaten Sikka Peta seismisitas Kabupaten Sikka dibuat berdasarkan kumpulan data kegempaan dari tahun 1900 sampai Oktober 2007. Gempa dalam (>300 m) ditunjukkan dengan node bulat hijau, gempa menengah (60-300 m) ditunjukkan dengan node bulat kuning, dan gempa dangkal (<60 m) ditunjukkan dengan node bulat merah. Besarnya ukuran node menunjukkan besarnya magnitudo gempa. Berdasarkan Gambar 15 di atas dapat dilihat bahwa sebagian besar gempa bumi terjadi di laut dan berada di daerah dangkal dengan kedalaman kurang dari 60 km dari dasar laut. Gempa bumi yang berpusat di laut dangkal merupakan pemicu timbulnya tsunami. Pusat gempa bumi dangkal di laut tersebut menyebar merata dan dominan di bagian utara dan selatan Pulau Flores. Pada bagian utara berpusat pada patahan naik Flores (Flores thrust) dan bagian selatan berpusat pada daerah subduksi Sunda Arc dan Banda Arc.

Frekuensi gempa dangkal relatif sering terjadi di zona-zona aktif gempa tersebut. Hampir sebagian besar gempa yang terjadi di Kabupaten Sikka, dipicu oleh gempa tektonik yang bermagnitudo lebih dari lima SR berpusat di laut. Besarnya frekuensi dan intensitas gempa tersebut menunjukkan bahwa Kabupaten Sikka rawan terhadap gempa dan tsunami. Kelompok pantai yang rawan terhadap tsunami adalah kelompok pantai yang langsung menghadap ke laut karena langsung berhadapan dengan daerah patahan subduksi yang merupakan sumber gempa besar dan sumber pembangkit tsunami. 4.2.2 Tingkat kerawanan tsunami di Kabupaten Sikka Berdasarkan sejarah kegempaan, diketahui Tsunami Flores 1992 merupakan tsunami terbesar yang pernah melanda wilayah NTT, khususnya Kabupaten Sikka. Oleh karena itu, untuk kajian resiko tsunami, data run up Tsunami Flores 1992 dijadikan sebagai acuan parameter kerawanan tsunami. Pada saat Tsunami Flores, ketinggian gelombang maksimum saat mencapai pantai adalah 10 m dengan inundasi lebih dari 600 m. Tinggi tsunami di setiap lokasi tidaklah sama (Gambar 16). Hal ini sangat ditentukan oleh kondisi fisik wilayah itu sendiri. Peta run up tsunami merupakan parameter penting dalam analisis resiko tsunami karena dijadikan sebagai input utama dalam parameter kerawanan tsunami. Peta run up ini dibuat berdasarkan data survei lapang Tsunami Flores 1992 (Lampiran 3). Pengolahan data run up tersebut dilakukan dengan metode manual untuk menspasialkan tingkat kerawanan daerah berdasarkan kelas tinggi run up tsunami.

Gambar 16. Peta tinggi run up Tsunami Flores tahun 1992 di berbagai lokasi Pantai Utara Kabupaten Sikka Berdasarkan data survei yang diperoleh diketahui bahwa tinggi run up tsunami di setiap tempat berbeda-beda. Tinggi run up maksimum yang terjadi di Kabupaten Sikka adalah 10 m, sedangkan yang terendah adalah 1.7 m. Run-up maksimum terjadi pada pantai yang landai dan berteluk. Run-up tsunami tidak hanya ditentukan oleh magnitudo gempa, tetapi juga sangat ditentukan oleh morfologi dasar laut, bentuk pantai serta karakteristik penjalaran gelombangnya (Fauzi, 2006). Selanjutnya, berdasarkan peta run up hasil survei lapang tersebut kemudian dibuat peta kerawanan tsunami berdasarkan matriks resiko tsunami (Gambar 17). Pada dasarnya peta tersebut menunjukkan tingkat kerawanan suatu daerah terhadap bencana tsunami.

Gambar 17. Peta kerawanan tsunami berdasarkan data run up Tsunami Flores tahun 1992

Setelah dilakukan pengolahan spasial tinggi run up tsunami, didapatkan tiga kelas kerawanan tsunami di Kabupaten Sikka berdasarkan skala Iida yaitu kelas kerawanan tinggi, kelas kerawanan sedang dan kelas kerawanan rendah. Berdasarkan peta run up tsunami (Gambar 17) di atas dapat diketahui bahwa hanya wilayah utara Kabupaten Sikka saja yang tergenang oleh limpasan tsunami, sedangkan daerah selatan tidak terkena limpasan tsunami. Hal ini disebabkan karena pusat tsunami terletak di wilayah utara Flores. Akan tetapi, dampak dari gempa dan tsunami tersebut dirasakan ke seluruh wilayah Kabupaten Sikka. Kelas kerawanan tinggi memiliki tinggi run up >6-16 m yang meliputi daerah Kecamatan Alok dan sebagian kecil wilayah Kecamatan Talibura. Selanjutnya, kelas kerawanan sedang memiliki tinggi run up >2-6 m, meliputi wilayah Kecamatan Waigete, Magepanda dan sebagian Kecamatan Talibura, Kewapante, dan Alok. Kerawanan rendah memiliki ketinggian >0.75-2 m, meliputi sebagian besar wilayah Kecamatan Kewapante. Kelas kerawanan tinggi merupakan daerah yang beresiko tinggi terhadap tsunami. Berdasarkan pengukuran dengan menggunakan tools measure diketahui jarak jangkauan tsunami (inundasi) pada wilayah kerawanan tinggi adalah ± 1500 m ke darat. Besarnya kemungkinan jangkauan limpasan tsunami yang masuk ke daratan tersebut menunjukkan pula batas kawasan rawan tsunami. Setiap kecamatan memiliki tinggi run up dan inundasi tsunami yang berbeda-beda. Hal ini sangat dipengaruhi oleh kondisi fisik lingkungan di masing-masing wilayah tersebut. Jika kerentanan lingkungannya tinggi, maka akan mudah untuk terpapar tsunami, sehingga resikonya pun akan lebih besar. Wilayah pesisir Sikka bagian utara adalah daerah yang memiliki potensi tinggi terkena tsunami, sedangkan daerah selatan relatif lebih aman.

Hal ini terutama disebabkan karena topografi wilayah utara Sikka lebih rendah daripada daerah selatan. Pengolahan data pada ArcGIS 9.1 dapat diketahui luas daerah kerawanan tsunami yaitu pada Tabel 7 berikut. Tabel 7. Luasan daerah kerawanan tsunami No. Tingkat kerawanan Jumlah sel Luas (m2)

Luas (Ha)

1

Sangat Tinggi

0

0

0

2

Tinggi

5.666

5.099.400

509,94

3

Sedang

34.102

30.691.800 306,18

4

Rendah

4.035

3.631.500

363,15

5

Sangat Rendah

0

0

0

Berdasarkan Tabel 7 di atas, diketahui bahwa luasan daerah kerawanan tinggi adalah 509,94 Ha; luasan daerah kerawanan sedang adalah 306,91 Ha; dan luasan daerah kerawanan sedang adalah 363,15 Ha. Berdasarkan Skala Iida (Tabel 4), diketahui bahwa sebagian besar tinggi Tsunami Flores termasuk dalam skala menengah. Skala Iida ini merupakan skala yang digunakan sebagai acuan parameter kerawanan dalam kajian resiko tsunami. Untuk mendapatkan tingkat resiko tsunami, maka parameter kerawanan akan diintegrasikan dengan data-data parameter kerentanan. 4.3 Analisis tingkat kerentanan tsunami di Kabupaten Sikka Daerah yang berpotensi tinggi rawan tsunami, belum tentu memiliki tingkat resiko yang tinggi terhadap tsunami. Besar kecilnya resiko tsunami sangat tergantung dari kondisi sosial, ekonomi, budaya, dan fisik lingkungan sekitarnya (Diposaptono dan Budiman, 2006). Pada penelitian ini, parameter kerentanan yang dikelaskan adalah parameter lingkungan, sedangkan parameter yang tidak dikelaskan adalah parameter sosial kependudukan dan infrastruktur wilayah.

Kajian parameter kerentanan lingkungan yang mempengaruhi tingkat resiko tsunami di Kabupaten Sikka adalah sebagai berikut. 4.3.1 Elevasi daratan Elevasi daratan merupakan parameter penting yang mempengaruhi tingkat kerentanan tsunami. Semakin rendah suatu daratan, akan mudah terkena limpasan tsunami dan sebaliknya, semakin tinggi suatu daratan, maka jangkauan tsunami pun akan semakin rendah. Jadi, tinggi dan rendah suatu daratan, sangat mempengaruhi jarak dan tinggi run up tsunami. Oleh karena itu, pada penelitian ini parameter topografi diberikan bobot sebesar 15%. Topografi Kabupaten Sikka terdiri atas daratan yang sempit di pesisir utara dan makin ke selatan wilayahnya berbukit-bukit dan bergunung-gunung, sehingga banyak terdapat pantai curam dan terjal. Pemetaan topografi dalam penelitian ini dilakukan secara digital menggunakan data topografi TIN yang disudah diubah ke dalam format raster. Pemetaan topografi pun dapat dilakukan dengan menggunakan data Digital Elevation Model (berformat raster). Untuk mempermudah dalam analisis topografi, penulis menggunakan tampilan data raster dengan pemodelan cahaya (sun angle shadding) atau hillshade. Hillshade adalah suatu tampilan data raster dengan pencahayaan matahari dari berbagai arah dan ketinggian sudut matahari. Efek dari tampilan hillshade adalah sisi permukaan yang terkena cahaya dan sisi yang tidak terkena pencahayaan, sehingga mempermudah melihat cepat profil ketinggian (ESRI, 2002). Pada penelitian ini, pemodelan cahaya dilakukan dengan perangkat lunak ArcGIS 9.1 dengan azimuth (sudut arah pencahayaan matahari) 315º dan elevation (ketinggian pencahayaan matahari) 45º dengan efek transparansi 67%.

Peta topografi yang sudah dikelaskan berdasarkan matriks, kemudian dioverlay dengan peta hillshade tersebut. Pemodelan hillshade ini diterapkan pada semua peta dalam penelitian ini, seperti yang sudah diaplikasikan pada peta daerah penelitian. Pemetaan elevasi dapat dilihat pada Gambar 18. Berdasarkan Gambar 18 dapat diketahui Kabupaten Sikka merupakan daerah yang memiliki topografi yang kompleks mulai dari dataran rendah sampai dengan topografi terjal yang mencapai ketinggian 1700 m di atas permukaan laut. Akan tetapi, secara umum dapat dilihat dari peta bahwa sebagian besar wilayah Kabupaten Sikka berbukit-bukit (topografi kasar). Hal ini ditandai dengan banyak terdapat dataran tinggi yang berbukit dan bergunung-gunung. Daerah utara Sikka merupakan daerah yang cukup landai dibandingkan dengan daerah selatan. Hal ini ditandai yaitu warna-warna ketinggiaan rendah lebih banyak terdapat di daerah utara daripada selatan. Kabupaten Sikka bagian utara terdiri atas daratan yang sempit dan makin ke selatan wilayahnya berbukitbukit dan bergunung. Untuk kajian resiko tsunami, topografi dibagi dalam lima kelas yaitu kelas kerentanan sangat tinggi (<10 m), kelas kerentanan tinggi (>10-25 m), kerentanan sedang (>25-50 m), kerentanan rendah (>50-100 m) dan kerentanan sangat rendah (>100 m).

Gambar 18. Peta kerentanan elevasi terhadap tsunami

Berdasarkan hasil perhitungan jumlah sel, dapat diketahui luasan wilayah per kelas kerentanan elevasi seperti pada Tabel 8. Tabel 8. Luasan wilayah tingkat kerentanan elevasi 2

No

Tingkat Kerentanan

Jumlah sel

Luas (m )

Luas (Ha)

1

Sangat Tinggi

51.324

46.191.600

4.619,16

2

Tinggi

74.295

66.865.500

6.686,55

3

Sedang

65.731

59.157.900

5.915,79

4

Rendah

125.148

112.633.200

11.263,32

5

Sangat Rendah

1.438.329

1.294.496.100

129.449,60

1.754.827

1.579.344.300

157.934,43

Total

Berdasarkan hasil pemetaan klasifikasi topografi berdasarkan matriks resiko tsunami, dapat diketahui bahwa daerah yang ketinggiannya kurang dari 10 m memiliki luas 4.649,16 Ha yang menyebar sepanjang pantai utara meliputi Kecamatan Talibura, Waigete, Kewapante, Alok, Maumere, dan Nita serta sebagian kecil terlihat di sepanjang pesisir selatan Kecamatan Paga dan Bola. Wilayah dengan ketinggian >10-25 m memiliki luas 6.686,55 Ha juga masih dominan berada pada pesisir utara dan sebagian di pesisir selatan. Untuk wilayah dengan ketinggian di atas 25-50 m dan 50-100 m masih dominan terdapat di wilayah utara Kabupaten Sikka, sedangkan ketinggian >100 m menyebar rata pada bagian tengah Kabupaten Sikka. Secara umum, dapat diketahui bahwa wilayah pesisir Kabupaten Sikka yang berelevasi rendah membuat tingkat kerentanan tsunami di daerah ini lebih tinggi daripada wilayah selatan. Semakin rendah elevasi suatu daerah, maka tingkat kerentanan terhadap bahaya tsunami semakin besar. Semakin besar tingkat kerentanan, maka semakin besar resikonya, dan sebaliknya. Oleh karena itu, peta topografi merupakan peta penting dalam kajian resiko tsunami.

4.3.2 Kemiringan daratan (slope) Kemiringan merupakan parameter penting dalam menentukan tingkat kerentanan tsunami di suatu daerah. Kemiringan daratan akan mempengaruhi tinggi run up tsunami yang akan terjadi. Semakin curam suatu daratan, maka tinggi run up akan semakin rendah. Satuan kemiringan daratan yang digunakan adalah dalam persentase (%). Range slope dalam persen berkisar dari 0-200%. Nilai kemiringan 0% mengindikasikan flat area/no slope (area datar). Nilai kemiringan 100% mengindikasikan kemiringan area 45º dan nilai 200% menunjukkan vertical slope (Earth Resource Mapping Ltd, 2008). Peta slope (Gambar 19) merupakan peta yang diturunkan dari peta topografi. Pada waktu processing data, data topografi dijadikan input dalam algoritma matematis, yang dapat mengubah setiap nilai elevasi menjadi sebuah nilai baru yang menggambarkan kemiringan lahan daratan dengan menggunakan fungsi surface analyst pada menu spatial analyst. Pada kajian resiko tsunami ini, parameter kemiringan daratan juga dibagi ke dalam lima kelas kerentanan yaitu kelas kerentanan sangat tinggi (2%), kerentanan tinggi (>2-10%), kerentanan sedang (>10-15%), kerentanan rendah (>15-40%) dan kerentanan sangat rendah (>40%). Hasil pemetaan slope (Gambar 19) menunjukkan bahwa kondisi kemiringan tanah (kelerengan) di wilayah Kabupaten Sikka cukup bervariasi. Berdasarkan laporan dari Kantor Bappeda Kabupaten Sikka (2005), kisaran kemiringan tanah di Kabupaten Sikka adalah dari 0 sampai dengan 70% dan didominasi oleh kemiringan tanah yang lebih besar dari 40% dengan luas 81,167 Ha atau sekitar 46,87% dari total luas wilayah Kabupaten Sikka.

Gambar 19. Peta kerentanan slope terhadap tsunami

Berdasarkan pengolahan data dengan ArcGIS 9.1 didapat luasan wilayah sebagaimana Tabel 9 di bawah ini. Tabel 9. Luasan wilayah kerentanan slope No Tingkat kerentanan Jumlah sel Luas (m2)

Luas (Ha)

1

Sangat Tinggi

74.995

67.495.500

6.749,55

2

Tinggi

187.314

168.582.600

16.858,26

3

Sedang

143.023

128.720.700

12.872,07

4

Rendah

710.134

639.120.600

63.912,06

5

Sangat Rendah

639.361

575.424.900

57.542,49

1.754.827

1.579.344.300 157.934,43

Total

Berdasarkan Tabel 9 di atas dapat dilihat bahwa daerah yang kerentanan sangat tinggi memiliki luas 6.749,55 Ha dan daerah yang kerentanannya tinggi memiliki luas 16.858,26 Ha. Daerah yang berada pada kelas kerentanan sangat tinggi dan tinggi, sebagian besar berada pada wilayah pesisir utara. Wilayahnya meliputi sebagian besar wilayah Kecamatan Alok, Maumere, Kewapante, Waigete, dan Talibura. Daerah ini akan merupakan daerah yang resiko tsunaminya tinggi jika tidak ada upaya penataan ruang yang baik terutama yang menyangkut area padat penduduk dan basis ekonomi penting. Berdasarkan hasil konsultasi dan studi literatur, diketahui bahwa penempatan kawasan pemukiman dapat dibagi dalam empat kelas kemiringan yaitu : 0-8º, 825º, 25-35º dan lebih dari 35º. Daerah yang memiliki tingkat kerentanan sangat tinggi dan tinggi umumnya terletak pada wilayah pesisir kabupaten dan dominan terlihat pada wilayah pesisir utara. Luas wilayah kerentanan sangat tinggi dan tinggi masing-masing adalah 6.749,55 Ha dan 16.858,26 Ha.

Daerah kerentanan sedang meliputi daratan tengah sebagian Kecamatan Lela, Bola, Talibura, Waigete, Kewapante dan Alok. Daerah kerentanan sangat rendah dan rendah, tersebar merata pada bagian tengah wilayah Kabupaten Sikka. Sehingga dapat diketahui bahwa daerah kerentanan rendah dan sangat rendah adalah daerah yang hampir mendominasi wilayah Kabupaten Sikka yang dicirikan dengan wilayah perbukitan dan pegunungan. Slope merupakan faktor penting yang menentukan sejauh mana limpasan tsunami ke daratan. Daerah yang memiliki slope yang landai merupakan daerah yang beresiko tinggi terhadap tsunami. Oleh karena itu, slope dijadikan sebagai salah satu parameter penting dan diberikan bobot sama seperti elevasi yaitu 15%. 4.3.3 Morfometri pantai Bentuk morfometri pantai sangat berpengaruh besar terhadap tingkat energi tsunami yang akan terhempas ke daratan. Meskipun besaran gempa tektonik yang mengakibatkan gelombang tsunami relatif kecil, tetapi jika morfometrinya mendukung untuk memberikan penguatan terhadap limpasan tsunami, maka resiko korban jiwa dan kerusakan akan semakin besar. Hal ini juga didukung bilamana sedimen pantainya berupa pasir halus. Hal ini juga akan menambah penguatan run up dan mengakibatkan jarak run up yang semakin jauh ke daratan (Istiyanto et al., 2005). Oleh karena itu peta morfometri ini diberikan bobot sebesar 15%. Pada penelitian ini, bentuk morfometri pantai diklasifikasikan dalam lima kelas yaitu teluk V, teluk U, tanjung, pantai lurus dan non teluk atau tanjung. Pemetaan morfometri pantai dilakukan dengan metode on screen digitize berdasarkan hasil visualisasi manual kenampakan topografi daerah pesisir.

Untuk keakuratan data, hasil pemetaan tersebut juga dibandingkan dengan data citra satelit. Pada pemetaan ini dilakukan teknik pembufferan sejauh 1.000 m dari garis pantai. Hal ini dilakukan dengan asumsi bahwa daerah limpasan tsunami tidak akan sampai sejauh satu km dari garis pantai. Peta morfometri pantai disajikan pada Gambar 20. Morfometri pantai berpengaruh terhadap besar kecilnya energi limpasan tsunami ke arah darat. Kabupaten Sikka memiliki morfometri yang unik. Berdasarkan Gambar 20 dapat dilihat bahwa daerah utara Sikka umumnya didominasi oleh teluk, sedangkan daerah selatan merupakan daerah yang didominasi oleh tanjung. Beberapa teluk yang terdapat di Kabupaten Sikka diantaranya yaitu Teluk Maumere, Teluk Rung, Teluk Paga, Teluk Wodong, dan Teluk Pedang. Masingmasing teluk tersebut, memiliki bentuk yang berbeda, dimana ada yang berbentuk U dan ada yang berbentuk V. Sedangkan beberapa tanjung yang terdapat di Kabupaten Sikka yaitu diantaranya Tanjung Wokar, Tanjung Watuntou, Tanjung, dan Tanjung Watunkelahi. Selain teluk dan tanjung, terdapat juga daerah pesisir dengan morfometri pantai lurus. Setiap morfometri pantai memberikan dampak yang berbeda terhadap limpasan gelombang tsunami di pantai. Pantai teluk V memiliki tingkat kerentanan yang lebih tinggi dibandingkan pantai teluk U. Hal ini karena, pantai teluk V lebih memusatkan energi tsunami. Akan tetapi, secara umum pantai yang berteluk memiliki tingkat kerentanan yang tinggi daripada tipe pantai yang lain. Hal ini karena morfometri pantai yang berbentuk teluk akan mempengaruhi refraksi gelombang tsunami sehingga kecepatan dan energi gelombang tsunami bertambah. Akibatnya, akumulasi massa air yang terjadi akan meningkatkan ketinggian serta kecepatan gelombang tsunami di pantai.

Gambar 20. Peta kerentanan morfologi terhadap tsunami

Jika tidak ada penghalang alami dan topografi pantainya landai, maka gelombang tsunami akan jauh masuk ke daratan. Berbeda dengan pantai berteluk, pada pantai bertanjung akan menyebarkan energi gelombang tsunami, sehingga penjalaran dan tinggi run up pun lebih rendah. Pada pantai yang lurus, energi merambat secara lurus tanpa ada pembelokan. Pantai tanjung termasuk ke dalam kerentanan sedang dan pantai lurus termasuk ke dalam kerentanan rendah. Daerah di luar buffer satu km termasuk ke dalam daerah non teluk atau tanjung. Daerah ini dianggap sudah tidak dipengaruhi oleh gelombang tsunami. Selain bentuk pantai, tipe pantai juga sangat mempengaruhi tingkat kerentanan terhadap tsunami karena tipe pantai (berpasir, berlumpur, berbatu, berkarang, atau berawa) turut mempengaruhi run up tsunami. Pada penelitian ini, penulis tidak mendapatkan data geomorfologi daerah Sikka. Akan tetapi, berdasarkan pengamatan di lapangan, sebagian besar daerah pesisir utara dan selatan Kabupaten Sikka pantainya berupa batuan, pecahan karang, tetapi terdapat juga pasir halus sampai kasar dan ada yang berasal juga dari pecahanpecahan karang (rubble). 4.3.4 Penggunaan lahan Pemetaan penggunaan tanah pada penelitian ini dilakukan dengan menggunakan data spasial dari Bappeda Kabupaten Sikka tahun 2005. Berdasarkan data yang didapat diketahui bahwa jenis penggunaan tanah yang terdapat di wilayah Kabupaten Sikka terdiri dari beberapa jenis penggunaan tanah yaitu batuan cadas, batuan gamping, hutan, hutan rawa, kebun, mangrove, pemukiman, rumput/tanah kosong, sawah irigasi, semak belukar, dan tegalan/ladang. Jenis penggunaan lahan di Kabupaten Sikka dapat dilihat pada Gambar 21.

Gambar 21. Peta jenis landuse di Kabupaten Sikka

Berdasarkan Gambar 21 di atas dapat diketahui luasan dari masing-masing landuse tersebut seperti Tabel 10 berikut. Tabel 10. Luasan jenis landuse di Kabupaten Sikka 2

Landuse

Jumlah sel

Luas m

Luas (Ha)

Batuan cadas

4.465

4.018.500

401,85

Batuan gamping

331

297.900

29,79

Hutan

232.081

208.872.900

20.887,29

Hutan rawa

2.974

2.676.600

267,66

Kebun

165.019

148.517.100

14.851,71

Mangrove

3.884

3.495.600

349,56

Pemukiman

29.593

26.633.700

2.663,37

Rumput/tanah kosong

64.949

58.454.100

5.845,41

Sawah irigasi

17.012

15.310.800

1.531,08

Semak/belukar

722.044

649.839.600

64.983,96

Tegalan/ladang

511.711

460.539.900

46.053,99

Total

1.754.063

1.578.656.700

157865,67

Dari Tabel 10 dapat diketahui bahwa tipe penggunaan tanah di wilayah Kabupaten Sikka didominasi oleh jenis semak belukar, tegalan/ladang, hutan belukar dan kebun. Dampak yang ditimbulkan oleh bencana tsunami terhadap masing-masing landuse tidak sama. Hal ini karena masing-masing landuse memiliki tingkat reduksi tertentu saat terkena gelombang tsunami. Misal untuk sawah irigasi. Sawah irigasi penting diketahui tingkat kerentanannya karena sawah merupakan sumber ekonomi pokok masyarakat. Apabila sawah irigasi terkena limpasan tsunami, maka areal sawah tersebut akan tergenang air laut dan tanah sawah yang terkena air asin akan menjadi tanah mati. Dampaknya adalah areal sawah tersebut tidak dapat digunakan lagi untuk bercocok tanam. Artinya, bencana tsunami dapat menyebabkan terjadi perubahan lahan. Oleh karena itu, perlu dilihat tingkat kerentanan landuse terhadap tsunami. Pemetaan kerentanan landuse di Kabupaten Sikka dapat dilihat pada Gambar 22.

Gambar 22. Peta kerentanan landuse terhadap tsunami

Berdasarkan Gambar 22 dapat diketahui bahwa daerah yang termasuk ke dalam kerentanan sangat tinggi dan tinggi banyak terdapat di wilayah pesisir, terutama pesisir utara Kecamatan Alok dan Kewapante. Hal ini karena pada daerah tersebut, banyak dimanfaatkan sebagai daerah pemukiman. Wilayah dengan tingkat kerentanan sedang, rendah dan sangat rendah sebagian besar terdapat pada bagian tengah Kabupaten Sikka. Luas dari masing-masing kelas kerentanan dapat terlihat pada Tabel 11 berikut. Tabel 11. Luasan tingkat kerentanan landuse 2

No

Tingkat Kerentanan

Jumlah sel

Luas (m )

Luas (Ha)

1

Sangat Tinggi

53.463

48.116.700

4.811,67

2

Tinggi

165.019

148.517.100

14.851,71

3

Sedang

511.711

460.539.900

46.053,99

4

Rendah

786.953

708.257.700

70.825,77

5

Sangat Rendah

236.917

213.225.300

21.322,53

1.754.063

1.578.656.700

157.865,67

Total

Tabel 11 di atas menunjukkan bahwa luasan tingkat kerentanan tinggi dan sangat tinggi berturut-turut adalah 14.851,71 Ha dan 4.811,67 Ha. Pada daerah ini, tingkat resikonya akan semakin besar. Berdasarkan Tabel 11 dapat dilihat juga tingkat penggunaan lahan di Kabupaten Sikka sebagian besar berada pada tingkat kerentanan sedang. Untuk kajian resiko tsunami, area pemukiman merupakan area yang paling rentan. Sebagian besar daerah pemukiman terletak di daerah pesisir yang berpotensi besar terhadap bahaya tsunami. Penggunaan lahan yang tidak banyak melibatkan manusia seperti hutan lebat, berada pada daerah yang aman. Oleh karena itu, penggunaan lahan pada kawasan pesisir harus memperhatikan konsep penataan ruang yang berbasis bencana alam, dalam hal ini adalah bencana tsunami.

4.3.5 Jarak dari garis pantai Mengingat tsunami bersifat merusak, maka dalam penataan ruang harus memperhatikan kawasan penyangga (buffer zone). Jarak dari garis pantai merupakan parameter penting dalam kajian resiko tsunami, sehingga diberi bobot 15%. Berdasarkan kejadian Tsunami Flores, banyak warga pesisir yang meninggal dunia dan berbagai sarana penting mengalami kerusakan karena berada di wilayah yang mudah terpapar tsunami. Oleh karena itu, pembuatan jarak dari garis pantai merupakan salah satu faktor penting dalam analisis kerentanan tsunami. Sebagian besar pemukiman dan sarana penting biasanya berada pada daerah pesisir yang sangat rentan terhadap tsunami. Hampir 140 juta penduduk (60%) penduduk Indonesia tinggal di wilayah pesisir dan sekitar 80% dari industri Indonesia memanfaatkan sumberdaya pesisir (Diposaptono dan Budiman, 2006). Oleh karena itu, penataan ruang wilayah pesisir harus melihat konsep jarak dari garis pantai. Pada penelitian ini pembuatan jarak dari garis pantai dilakukan dengan menggunakan data spasial garis pantai yang kemudian dilakukan buffering dan dikelaskan sesuai dengan matriks resiko tsunami. Peta jarak dari garis pantai dapat dilihat pada Gambar 23.

Gambar 23. Peta jarak dari garis pantai

Berdasarkan Gambar 23 di atas dapat diketahui bahwa daerah yang berwarna merah tua menunjukkan daerah yang paling rentan terhadap tsunami yang berada pada radius 200 m dari garis pantai. Daerah yang semakin dekat dengan pantai merupakan daerah yang paling rentan, dan sebaliknya. Berdasarkan pengamatan di lapangan, pemukiman yang rentan sekali terhadap tsunami yaitu pemukiman di daerah Wuring (Kecamatan Alok). Desa Wuring terletak dalam kawasan Teluk Maumere. Masyarakat pada daerah tersebut membangun rumahnya ke arah laut dan bahkan ada yang di laut. Areal pemukiman tersebut semakin lama, semakin bertambah banyak, sehingga semakin menjorok ke laut. Selain Wuring, areal pemukiman di daerah Desa Nangahale (Kecamatan Talibura) juga daerah-daerah pemukiman di pulau-pulau kecil Kabupaten Sikka, umumnya dekat dengan pantai. Dasar pondasi untuk pemukiman di daerah pesisir tersebut banyak yang terbuat dari rubble sehingga sangat rentan bila terhempas gelombang tsunami (Lampiran 2). Oleh karena itu, penting sekali menerapkan penataan ruang yang baik untuk mengurangi resiko tsunami khususnya di daerah pesisir. 4.3.6 Jarak dari sungai Jarak dari sungai (Gambar 24) juga merupakan parameter penting dalam kajian resiko tsunami. Pada umumnya tsunami yang melewati sungai akan menimbulkan kerusakan yang besar. Pada daerah yang menyempit seperti sungai, akan terjadi peningkatan kecepatan dan ketinggian muka air karena dengan debit massa air yang sama harus menjalar melalui celah yang sempit. Oleh karena itu, dalam penataan ruang harus memperhatikan jarak dari sungai.

Gambar 24. Peta jarak dari sungai

Pemetaan jarak dari sungai ini dibatasi sampai ketinggian 10 m di daratan sesuai dengan run up Tsunami Flores (10 m) dan pemetaan dilakukan pada sungai-sungai besar saja. Pada pemetaan di atas dapat dilihat bahwa Kabupaten Sikka memiliki sedikit sungai-sungai besar. Pada umumnya sungai-sungai besar tersebut terdapat pada daerah selatan dan jarak antara sungai pun sangat berjauhan. Berdasarkan literatur ilmiah, pada saat limpasan tsunami ke daratan, jika jarak antara dua sungai saling berdekatan akan menimbulkan kerusakan yang besar karena terjadi akumulasi energi gelombang tsunami dan massa air. Penempatan area pemukiman padat pada zona paling aman dari bahaya tsunami merupakan prioritas utama, sehingga harus diletakkan pada daerah yang jauh dari perkiraan jangkauan tsunami. Berdasarkan Gambar 23 di atas dapat dilihat bahwa daerah yang berwarna merah merupakan daerah yang memiliki tingkat kerentanan tinggi terhadap tsunami. Area padat penduduk dan area akonomi penting, sebaiknya berada pada jarak >100 m dari sungai. 4.4 Faktor-faktor pendukung tingkat kerentanan tsunami Parameter-parameter yang sudah dibahas sebelumnya merupakan parameter utama penentu tingkat resiko tsunami yang dikelaskan dalam matriks. Selain parameter di atas, terdapat parameter lain yang menentukan tingkat resiko tsunami yang dikaji dalam penelitian ini, yaitu batimetri, ekosistem pesisir, dan pulau penghalang. Parameter tersebut tidak dibobotkan dan dikelaskan dalam matriks sehingga tidak dioverlay pada pemodelan spasial tingkat resiko tsunami. 4.4.1 Batimetri Batimetri merupakan parameter yang penting dikaji karena mempengaruhi refraksi gelombang dan tinggi run up. Kecepatan dan energi gelombang tsunami

akan menurun sejalan dengan berkurangnya kedalaman, akan tetapi tinggi gelombang semakin meningkat. Batimetri Kabupaten Sikka dapat dilihat pada Gambar 25. Secara umum dapat dilihat bahwa batimetri Kabupaten Sikka termasuk dalam perairan dalam. Kisaran kedalaman perairan laut adalah 0-5200 m. Kedalaman dekat pantai umumnya dangkal dan semakin ke tengah laut, perairan semakin bertambah dalam. Perairan bagian utara yang merupakan wilayah perairan Flores, batimetrinya kasar karena banyak terdapat pulau-pulau, sedangkan pada bagian selatan batimetrinya halus.

0 -200 -400 -600 -800 -1000 -1200 -1400 -1600 -1800 -2000 -2200 -2400 -2600 -2800 -3000 -3200 -3400 -3600 -3800 -4000 -4200 -4400 -4600 -4800 -5000 -5200

t Lau

Saw

u

Gambar 25. Peta 3D batimetri Kabupaten Sikka Berdasarkan Gambar 25 di atas dapat diketahui bahwa daerah pesisir utara Sikka memiliki tingkat resiko yang lebih besar daripada daerah selatan. Hal ini karena pada daerah utara, batimetrinya perairan dekat pantai dalam dan curam. Jika batimetri curam, maka jarak daerah pecah gelombang dengan pesisir akan semakin kecil. Akan tetapi, energi yang dibawa oleh tsunami sangat besar.

Akibatnya pada daerah utara, run up tsunami akan semakin besar. Pada daerah selatan batimetrinya cenderung dalam dan landai. Pada daerah yang batimetrinya landai, jarak gelombang pecah dengan pesisir jauh, sehingga dampak kerusakan yang ditimbulkan akan lebih rendah dibandingkan dengan perairan yang curam. Berdasarkan pemetaan di atas dapat diketahui bahwa daerah pesisir utara Sikka lebih rentan terhadap tsunami daripada daerah selatan. 4.4.2 Ekosistem pesisir alami Kabupaten Sikka Keberadaan ekosistem pesisir merupakan pelindung alami terhadap ancaman gelombang tsunami. Yang termasuk dalam ekosistem pesisir alami adalah mangrove, lamun, dan terumbu karang. Data ekosistem pesisir ini, keseluruhannya diekstrak dari citra seperti yang sudah dijelaskan pada ekstraksi data citra. Ekosistem pesisir di Kabupaten Sikka dapat dilihat pada Gambar 14. Terumbu karang merupakan penghalang pertama di laut terhadap ancaman tsunami. Ekosistem terumbu karang banyak dijumpai di pesisir Kabupaten Sikka. Akan tetapi, ekosistem yang paling dominan adalah karang mati. Ekosistem lamun dan mangrove, meskipun jumlahnya sedikit tapi masih dijumpai di daerah pesisir. Ekosistem ini berperan penting dalam meredam energi gelombang tsunami Dalam kerangka mitigasi tsunami, ekosistem pesisir sering disebut “soft structure mitigation”. 4.4.3 Pulau-pulau kecil penghalang Keberadaan pulau penghalang (Gambar 26) terbukti sangat efektif dalam meredam energi gelombang tsunami. Berdasarkan kejadian Tsunami Flores 1992, gelombang datang dari arah barat laut dan dalam waktu tiga menit sudah menghantam pesisir utara Sikka. Daerah yang terhalangi oleh pulau penghalang, mengalami kerusakan yang tidak terlalu parah. Jarak, bentuk, dan ukuran pulau

sangat menentukan tingkat resiko tsunami. Sebagai contoh, bencana Tsunami Flores paling parah dialami oleh Pulau Bater yang berjarak lima km dari pesisir utara Flores. Pulau tersebut berbentuk konikal. Bentuk pulau ini menyebabkan gelombang yang datang mengalami refraksi dan difraksi serta terperangkapnya energi gelombang di pulau tersebut (wave trapping) (Diposaptono dan Budiman, 2006).

Gambar 26. Peta pulau penghalang 4.5 Analisis spasial tingkat resiko tsunami di Kabupaten Sikka Analisis resiko tsunami merupakan analisis yang penting dalam kerangka mitigasi bencana alam. Karena mitigasi baru akan diambil setelah diketahui tingkat resikonya. Analisis resiko tsunami yang ideal mencakup fungsi sistematis mulai dari mengidentifikasi perilaku alami, lokasi, intensitas, dan kemungkinan berulangnya suatu kejadian bencana, menentukan keberadaan dan tingkat

kerentanan, mengidentifikasi tingkat kapasitas dan sumberdaya yang tersedia, serta menentukan tingkat resikonya yang biasanya disajikan dalam bentuk petapeta dan analisis kawasan spesifik (site specific). Kebijakan yang diambil berdasarkan GIS juga banyak diyakini banyak pihak lebih akurat dan tepat sasaran (Diposaptono dan Budiman, 2006). Pada penelitian ini, penulis memetakan secara spasial tingkat resiko tsunami di Kabupaten Sikka dengan metode Cell Based Modeling (CBM). Analisis spasial pada data raster merupakan dasar dari pemodelan spasial berbasis sel. Hal ini karena setiap sel memiliki nilai tertentu, sehingga akan mempermudah dalam analisis spasial. Pemodelan tingkat resiko tsunami dilakukan dengan menspasialkan terlebih dahulu keseluruhan parameter utama yang mempengaruhi resiko tsunami. Setiap parameter utama yang sudah berformat raster, direklasifikasi menjadi kelas kerawanan dan kerentanan. Pengelompokkan setiap parameter tersebut mengikuti operasi zonal function karena setiap parameter akan mengelompok berdasarkan kesamaan nilai sel tersebut. Setelah dikelompokkan, pengkodean sel (calculation) dilakukan secara otomatis oleh perangkat lunak ArcGIS 9.1 menurut selang nilai parameter yang ditentukan. Apabila keseluruhan parameter sudah disiapkan, maka dilakukan proses overlay berdasarkan pembobotan sesuai Tabel 5. Oleh karena itu, overlay yang digunakan dalam penelitian ini adalah overlay dengan sistem pembobotan (weighted overlay). Setelah direklasifikasi jumlah sel untuk setiap parameter disajikan pada Tabel 12.

Tabel 12. Jumlah sel per parameter hasil klasifikasi Parameter

Jumlah sel Sangat tinggi

Tinggi

Sedang

Total Rendah

Sangat rendah

Elevasi

51.324

74.295

65.731

125.148

1.438.329

1.754.827

Slope

74.995

187.314

143.023

710.134

639.361

1.754.827

Morfologi

18.691

66.022

64.631

1.531.539

73.987

1.754.870

Landuse

53.463

165.019

511.711

786.953

236.917

1.754.063

Jarak dari gp

47.354

67.432

108.554

10.6175

1.425.355

1.754.870

Jarak dari sungai

1.911

2.207

2.464

5.793

1.742.495

1.754.870

Run up tsunami

0

5.666

34.102

4.035

0

43.803

Overlay data raster dilakukan dengan menggunakan fungsi raster calculator pada menu spatial analyst pada perangkat lunak ArcGIS 9.1. Setiap layer parameter akan dikalikan dengan bobotnya masing-masing. Secara matematis, overlay setiap layer parameter dapat ditulis sebagai berikut.

[( run up * 0,20)  ( elevasi * 0,15)  ( slope * 0,15)  (morfometri * 0,15)  (landuse * 0,10)  ( jarak dari garis pantai * 0,15)  ( jarak dari sungai * 0,10)] Setelah proses overlay, akan didapatkan lima kelas tingkat resiko tsunami , yaitu resiko sangat tinggi, resiko tinggi, resiko sedang, resiko rendah dan resiko sangat rendah. Luasan masing-masing kelas resiko ini disajikan pada Tabel 13. Berdasarkan Tabel 13 dapat diketahui bahwa luasan daerah resiko sangat tinggi, tinggi, sedang, rendah dan sangat rendah berturut-turut adalah 34,04 Ha; 2.733,93 Ha; 3.615,3 Ha; 12.825 Ha; dan 138.447,6 Ha. Tabel 13. Jumlah sel hasil weighted overlay No

Kelas resiko

Jumlah sel

Luas (Ha)

1

Sangat Tinggi

356

32,04

2

Tinggi

30.377

2.733,93

3

Sedang

40.170

3.615,3

4

Rendah

142.500

12.825

5

Sangat Rendah

1.538.307

138.447,6

1.751.710

157.653,9

Total

Daerah resiko sangat tinggi terdapat di dua wilayah pesisir utara kecamatan yaitu Kecamatan Alok dan Magepanda. Daerah resiko tinggi terdapat sebagian besar di pesisir utara. Daerah resiko sedang, rendah dan sangat rendah terdapat di sepanjang bagian tengah wilayah utara dan selatan Kabupaten Sikka. Berdasarkan Tabel 13 dapat diketahui luasan daerah resiko sangat tinggi adalah 32,04 Ha dari total luas wilayah Kabupaten Sikka dan luasan resiko tinggi adalah 2.733,93 Ha dari total luas wilayah Kabupaten Sikka. Pada pemodelan spasial daerah resiko tsunami dapat dilihat bahwa pada area di tengah Kabupaten Sikka masih terdapat warna biru yang menunjukkan daerah resiko rendah. Hal tersebut disebabkan oleh pengaruh dari parameter slope, dimana pada daerah Kecamatan Alok sampai Maumere slopenya masih landai. Sehingga area tersebut masih termodelkan sebagai daerah yang beresiko rendah. Hal ini menunjukkan bahwa faktor pembobotan sangat penting kontribusinya dalam memetakan suatu model yang ideal yang sesuai pada tingkat lokal. Peta tingkat resiko tsunami Kabupaten Sikka dapat dilihat pada Gambar 27.

Gambar 27. Peta tingkat resiko tsunami Kabupaten Sikka

Berdasarkan Gambar 27 juga dapat dilihat bahwa tingkat resiko tsunami di setiap wilayah kecamatannya berbeda-beda. Hal ini sangat tergantung dari tingkat kerentanan yang ada di wilayah tersebut. Tingkat resiko tsunami setiap kecamatan dapat dilihat pada Tabel 14. Tabel 14. Tingkat resiko tsunami per kecamatan Kecamatan

Persentase tingkat resiko (%) Sangat tinggi

Tinggi

Sedang

Rendah

Total (%) Sangat rendah

Alok

0.44

9.24

6.06

17.92

66.33

100.00

Bola

0.00

0.00

0.72

8.71

90.57

100.00

Kewapante

0.00

3.86

5.06

13.39

77.69

100.00

Lela

0.00

0.23

5.45

24.32

70.00

100.00

Magepanda

0.01

2.93

3.60

11.66

81.80

100.00

Maumere

0.00

0.00

0.00

6.72

93.28

100.00

Mego

0.00

0.01

1.23

3.13

95.62

100.00

Nita

0.00

0.00

0.00

1.79

98.21

100.00

Paga

0.00

0.02

0.58

4.96

94.44

100.00

Talibura

0.00

2.38

2.27

6.28

89.06

100.00

Waigete

0.00

1.62

3.23

8.69

86.45

100.00

Berdasarkan Tabel 14 dapat diketahui bahwa kecamatan yang beresiko paling tinggi terhadap tsunami terdapat pada Kecamatan Alok yaitu sebesar 0,44% dari total luas Kecamatan Alok. Daerah dengan tingkat resiko tinggi yang terbesar terdapat pada Kecamatan Alok juga yaitu sebesar 9,24% dari total luas Kecamatan Alok. Hal ini disebabkan wilayah Kecamatan Alok berada pada elevasi dan slope yang rendah. Dalam rangka mitigasi tsunami, daerah resiko sedang, rendah dan sangat rendah dijadikan sebagai area evakuasi jika terjadi tsunami. Secara visual tingkat persentase resiko tsunami dapat dilihat pada Gambar 28 berikut. Pemetaan resiko per kecamatan dapat dilihat pada Lampiran 6.

Gambar 28. Grafik resiko tsunami per kecamatan Untuk melihat tingkat resiko tsunami terhadap berbagai sarana dan prasarana sosial ekonomi masyarakat, dilakukan pemetaan tingkat resiko kerusakan infrastruktur fisik dan ekonomi dengan mengoverlay antara peta resiko dengan data infrastruktur yang ada di Kabupaten Sikka. Hal ini penting dilakukan terkait dengan pentingnya keberadaan sarana prasarana penting guna mendukung aktivitas penduduk. Sarana dan prasarana penting yang dipetakan yaitu pusat listrik, sarana kesehatan, sarana pendidikan, sarana peribadatan, sarana ekonomi, pelabuhan dan lapangan udara. Berdasarkan overlay dengan peta infrastruktur dapat diketahui sarana dan prasarana yang berpotensi terkena dampak tsunami (Gambar 29).

Gambar 29. Peta resiko infrastruktur terhadap tsunami

Berdasarkan Gambar 29 dapat diketahui bahwa sebagian besar infrastruktur penting berada pada daerah beresiko sangat tinggi dan tinggi. Dengan demikian, berdasarkan hasil pemetaan di atas, diharapkan dapat dijadikan sebagai masukan bagi pengambil kebijakan untuk mengambil upaya dan strategi mitigasi yang sesuai pada daerah yang bersangkutan. Untuk melihat tingkat resiko tsunami terhadap data sosial kependudukan dengan mengoverlay peta resiko tsunami dengan sebaran pemukiman. Sebaran dan kepadatan pemukiman menjadi salah satu faktor yang mempengaruhi resiko bencana tsunami yang akan terjadi. Pemukiman penduduk menggambarkan tingkat kepadatan penduduk dan sebaran tempat hunian yang akan mempengaruhi tingkat kerugian akibat tsunami baik dari segi kerugian jiwa maupun harta benda. Tingkat resiko penduduk tersebut dapat dilihat pada Gambar 30. Berdasarkan Gambar 30 dapat dilihat bahwa sebagian besar masyarakat tinggal di daerah dekat pantai, daerah yang bertopografi rendah serta kemiringannya landai. Pemukiman penduduk tersebar di keseluruhan wilayah, tetapi kepadatan pemukiman di tiap wilayah berbeda-beda. Kecamatan yang terpadat adalah Kecamatan Alok dengan kepadatan penduduk 761,76 jiwa/km2. Lalu diikuti dengan kecamatan Kewapante dengan kepadatan penduduk 460,52 jiwa/km2. Sedangkan kepadatan penduduk terendah terdapat di Kecamatan Talibura yaitu 62,62 jiwa/km2 (Bappeda Sikka,2006). Oleh karena itu, penduduk yang tinggal di wilayah Kecamatan Alok dan Kewapante memiliki resiko terhadap tsunami yang lebih tinggi dibandingkan dengan penduduk di kecamatan lain.

Gambar 30. Peta resiko pemukiman terhadap tsunami

Secara fisik lingkungan daerah-daerah padat penduduk yang berada dekat dengan pantai, aliran sungai, serta mempunyai elevasi yang rendah, akan sangat rentan terhadap bahaya tsunami. Sehingga tingkat resikonya pun akan semakin besar. Oleh karena itu, dengan mengetahui tingkat kerentanan di setiap wilayah dapat dilakukan prioritas upaya pengurangan resiko bencana secara lebih spesifik untuk masing-masing variabel. Mitigasi dapat dilakukan jika tingkat kerentanan dan resiko tsunami sudah ditentukan. Berdasarkan hasil wawancara dengan pakar, acuan yang jelas tentang peta resiko tsunami masih sedang dibuat. Selain itu, penulis juga memperoleh informasi bahwa, sekarang ini peta resiko bencana dijadikan dasar untuk semua perencanaan keruangan yang berada di wilayah pesisir. Rencana Tata Ruang Wilayah (RTRW) ataupun Rencana Tata Ruang Detail Kota (RTDK) baik di tingkat propinsi , kabupaten, kecamatan harus berbasis bencana. Oleh karena itu, kajian mendalam tentang analisis resiko tsunami penting untuk dilakukan di suatu wilayah tertentu. 4.6. Pemodelan resiko tsunami pada run up tsunami 20 meter Pada penelitian ini dilakukan juga pemodelan tingkat resiko tsunami di Kabupaten Sikka dengan menggunakan run up tsunami 20 meter. Hal ini dilakukan sebagai salah satu upaya mitigasi terhadap tsunami dengan membuat perkiraan daerah rawan adalah 2 kali run up maksimum yang dicapai saat Tsunami Flores. Perkiraan daerah resiko berdasarkan run up 20 m dapat dilihat pada Gambar 31.

Gambar 31. Peta tingkat resiko tsunami pada run up 20 meter

Berdasarkan Gambar 31 dapat dilihat bahwa daerah yang beresiko tinggi terhadap tsunami terutama adalah di pesisir utara Kecamatan Alok, Talibura dan Magepanda. Luasan daerah tingkat resiko pun bertambah karena daerah yang rawann semakin luas. Besarnya luasan area per tingkat resiko tsunami dapat dilihat pada Tabel 15. Tabel 15. Luasan wilayah tingkat resiko tsunami pada run up 20 meter No

Kelas resiko

Jumlah sel

Luas (Ha)

1

Sangat Tinggi

6.154

553,86

2

Tinggi

65.800

5.922

3

Sedang

44.375

3.993,75

4

Rendah

86.974

7.827,66

5

Sangat Rendah

1.548.680

139.381,2

1.751.983

157.678,5

Total

Berdasarkan Tabel 15 di atas dapat diketahui bahwa luasan wilayah untuk tiap kelas resiko berbeda-beda. Luasan daerah resiko sangat tinggi, tinggi, sedang, rendah dan sangat rendah berturut-turut adalah 553,86 Ha; 5.922 Ha; 3.993,75 Ha; 7.827,66; dan 139.381,2 Ha. Apabila dibandingkan dengan pemetaan resiko pada run up 10 meter dapat dilihat bahwa luasan daerah resiko tinggi dan sangat tinggi pada run up 20 meter lebih luas daripada run up 10 meter. Secara matematis diketahui bahwa luasan daerah beresiko sangat tinggi dan tinggi pada run up 10 m adalah 1,7086% dari total luas wilayah Kabupaten Sikka. Sedangkan luasan daerah beresiko sangat tinggi dan tinggi pada run up 20 m adalah sebesar 4 % dari total luas wilayah Kabupaten Sikka. Hal ini menunjukkan bahwa luasan daerah yang beresiko tinggi dan sangat tinggi terhadap bahaya tsunami adalah relatif lebih kecil dari total luas wilayah Kabupaten Sikka. Akan tetapi, lokasi kemungkinan kejadian tsunami tersebut

berada pada daerah pemukiman padat penduduk dan terdapat banyak infrastruktur penting. Oleh karena itu, pada daerah-daerah yang beresiko tinggi dan sangat tinggi inilah, upaya mitigasi penting untuk diterapkan. Pengembangan mitigasi tsunami yang komprehensif di Kabupaten Sikka dapat dilakukan secara struktural diantaranya yaitu dengan membangun breakwater. Selain itu, dapat juga dilakukan upaya mitigasi non struktural, yaitu salah satunya dengan membuat peta evakuasi terutama pada wilayah-wilayah yang beresiko tinggi dan sangat tinggi terhadap tsunami.

5. KESIMPULAN DAN SARAN 5.1 Kesimpulan Berdasarkan analisis SIG dengan metode CBM diperoleh bahwa tingkat resiko tsunami di setiap wilayah kecamatan di Kabupaten Sikka berbeda-beda. Daerah resiko sangat tinggi terdapat di dua wilayah pesisir utara kecamatan yaitu Kecamatan Alok dan Magepanda dengan luasan daerah resiko sangat tinggi adalah 32,04 ha dan luasan resiko tinggi adalah 2.733,93 ha. Secara umum dapat diketahui bahwa luasan yang beresiko sangat tinggi dan tinggi terhadap bahaya tsunami adalah relatif lebih kecil (1,7086% pada run up 10 m dan 4% pada run up 20 m) dari total luas wilayah Kabupaten Sikka. Akan tetapi, lokasi kemungkinan kejadian tsunami tersebut berada pada daerah pemukiman padat penduduk sehingga perlu dikembangkan mitigasi tsunami yang komprehensif di Kabupaten Sikka, terutama pada wilayah-wilayah yang beresiko tinggi dan sangat tinggi terhadap tsunami. 5.2 Saran Saran yang dapat diberikan dalam penelitian ini yaitu : 1. Perlu dilakukan penelitian lebih lanjut dengan memasukan data pemodelan tsunami dengan model pusat gempa tsunami dari arah yang berbeda untuk memetakan daerah resiko tsunami yang lebih detail. 2. Penerapan ”soft mitigation dan hard mitigation structure” penting untuk dilakukan pada daerah yang beresiko tinggi terhadap bahaya tsunami. 3. Perlu dibuat “evacuation zone route map” dan database bencana alam tsunami sebagai input penting untuk perkembangan daerah berbasis bencana.

DAFTAR PUSTAKA Badan Meteorologi dan Geofisika. 2007. Tsunami. http://www.bmg.go.id [6 Juni 2007] Badan Koordinasi Nasional Penanggulangan Bencana. 2007. Fasilitasi Rencana Penanggulangan Bencana : Materi 2. Bakornas PB. Jakarta www.bakornas_pb.go.id [19 Juni 2007] Bakornas Penanggulangan Bencana. 2008. Operasionalisasi Program Penanganan Bencana Alam Bidang Penataan Ruang. www.bakornas_pb.go.id [19 Juni 2008] Badan Koordinasi Survey dan Pemetaan Nasional. 2002. Inventarisasi Kerawanan Bencana Alam Wilayah Pesisir dan Laut Indonesia. Bakosurtanal. Cibinong, Bogor. BIOTROP Training and Information Center. 2005. Training Course on Remote Sensing Satellite Data Processing For Land Feature Mapping and Monitoring. BTIC BIOTROP. Bogor. Diposaptono, S dan Budiman. 2006. Tsunami. Buku Ilmiah Populer. Bogor. Diposaptono, S. 2006. Mitigasi Bencana Wilayah Pesisir Berbasis Ekosistem Mangrove. Ditjen Kelautan, Pesisir, dan Pulau-pulau Kecil DKP. Jakarta. Earth Observatory. 2007. Measuring Vegetation : Normalized Difference Vegetation Index (NDVI). http://[email protected] [26 April 2008] Environmental System Research Institute. 2001. Using ArcGis Spatial Analyst. ESRI Institute, Inc. New York. USA. Environmental System Research Institute. 2002. Using ArcGis Spatial Analyst. ESRI Institute, Inc. New York. USA. Earth Resource Mapping Inc. 2008. User Guide : Slope. www.ermapper.com [8 Agustus 2008]

Fauzi. 2006. Daerah Rawan Gempa Tektonik di Indonesia. http://www.bmg.go.id [5 Juni 2007] Fauzi dan Wandono. 2005. Penentuan Episenter Untuk Analisa Resiko Limpasan Tsunami in Penerapan Hasil Riset Untuk Penanggulangan Bencana Tsunami di Indonesia (Sadikin et al. eds). Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on Indonesian Coastal Region (2002-2003-2004). BPPT Press. Jakarta. Pp 163-169. GIS Konsorsium Aceh Nias. 2007. Modul Pelatihan ArcGIS Tingkat Dasar. Pemerintah Kota Banda Aceh. Nangroe Aceh Darussalam Hajar, M. 2006. Penentuan Kawasan Rawan Bencana Tsunami di Kota Padang, Sumatera Barat. Skripsi (Tidak Dipublikasikan). Program Studi Ilmu dan Teknologi Kelautan. Fakultas Perikanan dan Ilmu Kelautan. Institut Pertanian Bogor. Bogor. Harada, K dan Imamura, F. 2003. Study on The Evaluation of Tsunami Reducing by Coastal Control Forest for Actual Conditions. Asia and Pasific Coast. Japan Hendratno, A. 2005. Peta Resiko Limpasan Tsunami Dalam Perbaikan Tata Ruang Pantai in Penerapan Hasil Riset Untuk Penanggulangan Bencana Tsunami di Indonesia (Sadikin et al. eds). Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on Indonesian Coastal Region (2002-2003-2004). BPPT Press. Jakarta. Pp. 235-242. Himpunan Ahli Geofisika Indonesia (HAGI). 2007. Pengenalan Tsunami. www.hagi.or.id [5 Januari 2007] Iida, K. 1963. Magnitude, Energy And Generation Mechanism Of Tsunami And A Catalogue Of Earthquacke Associated With Tsunami. International Union Geodesy and Geophysics Monograph. 24 : 7-17 ITIC. 2006. Tsunami Glossary. http://www.shoa.cl/oceano/itic/frontpage/html [12 Februari 2006] Istiyanto, D. C., K. S. Utomo, Suranto dan M. Z. Jauzi. 2005. Studi Model Fisik Perambatan dan Limpasan Tsunami pada Berbagai Variasi Bentuk Pantai. dalam Penerapan Hasil Riset Untuk Penanggulangan Bencana Tsunami di Indonesia (Sadikin et al. eds). Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on

Indonesian Coastal Region (2002-2003-2004). BPPT Press. Jakarta. Pp. 81-92. Irfani, M. 2005. Pola Lansekap di Kawasan Rawan Tsunami in Penerapan Hasil Riset Untuk Penanggulangan Bencana Tsunami di Indonesia (Sadikin et al. eds). Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on Indonesian Coastal Region (20022003-2004). BPPT Press. Jakarta. Pp. 265-274. Latief. 2005. Prakiraan Pembangkitan Tsunami Menyusul Kejadian Gempa (Studi Kasus Tsunami Flores 1992) in Penerapan Hasil Riset Untuk Penanggulangan Bencana Tsunami di Indonesia (Sadikin et al. eds). Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on Indonesian Coastal Region (2002-2003-2004). BPPT Press. Jakarta. Pp. 27-52. Lillesand, T. M. dan F. W. Kiefer. 1990. Penginderaan Jauh dan Interpretasi Citra. Alih Bahasa : Dulbahri, Prapto Suharsono, Hartono, Suharyadi; Sutanto (penyunting). Gadjah Mada University Press, Yogyakarta. Meaden, G. J. dan T. D. Chi. 1996. Geographical information Systems : Application to Marine Fisheries. FAO Fisheries Technical Paper No. 356. Rome. Melianawati, B. D. 2006. Tsunamy Risk Assesment Using Probabilistic Approach and GIS Based SAW Method : Case Study Coastal of Padang, West Sumatera, Indonesia. Graduated School Bogor Agriculture University. Bogor. National Atmospheric and Space Administration (NASA). 2005. Landsat 7 Science Data Users Handbook. http://landsathandbook.gsfc.nasa.gov/handbook.html [6 Agustus 2007] Pasek, I. M. R. S. 2007. Penentuan Zona Potensial Budidaya Mutiara (Pinctada spp.) dengan Cell Based Modelling : Studi Kasus Perairan Sekotong, Lombok Barat NTB. Skripsi (Tidak Dipublikasikan). Program Studi Ilmu dan Teknologi Kelautan. Fakultas Perikanan dan Ilmu Kelautan. Institut Pertanian Bogor. Bogor. Poling, P. 2005. Sikka Dalam Kenangan 12 Desember 1992 dalam Penerapan Hasil Riset Untuk Penanggulangan Bencana Tsunami di Indonesia (Sadikin et al. eds). Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on Indonesian Coastal Region (2002-2003-2004). BPPT Press. Jakarta. Pp 3-14.

Prahasta, E. 2002. Konsep-Konsep Dasar Sistem Informasi Geografis. Informatika. Bandung. Pratikno, W. A. 2007. Disusun, Tata Ruang Kota Pesisir Rawan Bencana. www.bppt.go.id [7 Januari 2008] Puntodewo, A. 2003. Sistem Informasi Geografis untuk Pengelolaan Sumber Daya Alam. Center for International Forestry Research. Bogor. Purwadhi, F.S.H. 2001. Interpretasi Citra Digital. Gramedia. Jakarta. Puspito, N. 2007. Indonesia memang Rawan Tsunami. www.bppt.go.id [12 April 2007] Pusat Informasi Riset Bencana Alam. 2008. Pembuatan Peta Resiko. http://www.pirba.ristek.go.id/main.php [14 Oktober 2008] Pusat Riset Teknologi Kelautan. 2004. Penanganan Abrasi, Erosi dan Tsunami dengan Optimasi Vegetasi. Badan riset Kelautan dan Perikanan Departemen Kelautan dan Perikanan. Jakarta. Siahaan, F. R. 2008. Perlu Manajemen Resiko Bencana Tsunami Di Mentawai. http://www.puailiggoubat.com [15-31 Mei 2008] Siregar, V., D. Soedharma dan J. L. Gaol. 1995. Pemetaan Tematik Terumbu Karang dengan Menggunakan Data Landsat TM. Fakultas Perikanan IPB. Bogor. Soegiharto, R. 2006. Mitigasi Bencana Di Kampung Nelayan : Upaya Sistematis Mengurangi Kerugian Jiwa, Harta Benda dan Kerusakan Lingkungan. www.dkp-banten.go.id [12 April 2007] Susilo, S. B. 2000. Penginderaan Jauh Kelautan Terapan. Fakultas Perikanan dan Ilmu Kelautan IPB. Bogor. Susilo, S. B. dan J.L. Gaol. 2008. Dasar-dasar Penginderaan Jauh Kelautan. Fakultas Perikanan dan Ilmu Kelautan IPB. Bogor. UNESCO-IOC. 2006. Tsunami Glossary: IOC Information document No. 1221. Paris. www.shoa.cl [26 April 2008]

United State Geologyc Survey (USGS). 2007. Seamless Shuttle Radar Topography Mission (SRTM). http://rockyweb.cr.usgs.gov/nmpstds/demstds.html [1 Mei 2007] United State Geologyc Survey (USGS). 2007. Data Tsunami. http://neic.usgs.gov/neic/epic/epic_rect.html [28 Agustus 2007] Wah. 2006. NTT Rawan Gempa dan Tsunami. http://www.kompas.co.id [12 Januari 2007] Yuwono. 2005. Karakteristik Hidraulik Hutan Bakau dalam Meredam Energi Tsunami. Prosiding Seminar Tsunami dalam Kerangka Research on Tsunami Hazard and Its Effect on Indonesian Coastal Region (20022003-2004). BPPT Press. Jakarta. Pp 323-334.

LAMPIRAN

Lampiran 1. Posisi pengamatan di lapangan Waypoint 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

Longitude

Latitude

Elevasi (m)

122.21978 122.21942 122.21942 122.21942 122.21942 122.21942 122.7133 122.63601 122.30133 122.20351 122.20351 122.20227 122.21902 122.21902 106.73486

-8.61468 -8.61507 -8.6151 -8.6151 -8.6151 -8.6151 -8.52242 -8.50106 -8.63926 -8.59848 -8.59848 -8.601 -8.61777 -8.61777 -6.55469

5.49 2.74 4.69 4.69 4.69 4.69 39.77 307 148 6.13 6.13 3.97 1.57 1.57 200

Lampiran 2. Foto-foto hasil survei lapang

Pemukiman di laut (Wuring)

Jarak rumah <5 m dari garis pantai (Nangahale)

Pondasi rumah dari rubble (Nangahale)

Ekosistem mangrove (Wuring)

Vegetasi pantai di P. Babi (Barier Island)

Stasiun Pasut Bakosurtanal di Maumere

Lampiran 3. Data run up Tsunami Flores tahun 1992

Lampiran 4. Root mean square error koreksi Citra Landsat 7/ETM+ 2006 # GCPs for dataset : D:\Processing data\#DATA\iMAGE_2006\Citra_2006.ers # Total number of GCPs: 30 # Number turned on : 30 # Warp order :1 # GCP CORRECTED map projection details: # Map Projection : SUTM51 # Datum : WGS84 # Rotation : 0.000 # RMS error report: # Warp Type - Polynomial # -----ACTUAL----- ---PREDICTED--# Point Cell-X Cell-Y Cell-X Cell-Y RMS # "1" 1887.083 1393.253 1886.674 1393.103 0.4361 # "2" 1971.893 1548.266 1971.943 1548.409 0.1509 # "3" 2000.504 1536.869 2000.600 1537.252 0.3949 # "4" 1591.964 1392.319 1591.979 1392.283 0.0392 # "5" 2226.264 1350.979 2226.316 1351.364 0.3886 # "6" 2578.870 1234.771 2578.692 1234.813 0.1824 # "7" 2463.114 1392.603 2463.096 1392.605 0.0182 # "8" 2167.035 1927.095 2166.903 1927.228 0.1875 # "9" 2148.665 1940.020 2149.278 1940.417 0.7302 # "10" 1513.977 2118.522 1513.958 2118.437 0.0879 # "11" 1343.114 2122.766 1343.120 2122.818 0.0521 # "12" 689.840 1731.502 690.249 1731.398 0.4216 # "13" 604.233 1710.468 604.676 1710.355 0.4574 # "14" 3011.829 2119.217 3012.197 2118.813 0.5462 # "15" 2196.901 2444.678 2196.269 2444.246 0.7654 # "16" 1974.220 2491.784 1974.549 2491.088 0.7697 # "17" 1910.677 2526.427 1910.654 2526.407 0.0308 # "18" 1718.464 2557.584 1718.286 2557.757 0.2485 # "19" 1145.117 2542.399 1145.141 2542.648 0.2508 # "20" 611.243 2591.961 611.028 2592.237 0.3497 # "21" 445.050 2728.221 444.854 2728.362 0.2416 # "22" 1623.259 2144.108 1623.220 2144.102 0.0391 # "23" 598.763 2541.554 598.838 2541.703 0.1666 # "24" 2765.451 2132.513 2765.344 2132.467 0.1163 # "25" 2922.702 2092.173 2922.523 2092.223 0.1858 # "26" 1407.323 2564.805 1407.378 2565.066 0.2663 # "27" 1084.917 2549.076 1085.037 2549.016 0.1340 # "28" 208.831 1428.068 207.984 1427.798 0.8890 # "29" 415.492 1591.059 415.578 1590.946 0.1420 # "30" 518.425 1708.782 518.856 1708.478 0.5270 # Average RMS error : 0.307 # Total RMS error : 9.216 # End of GCP details

Lampiran 5. Hasil analisis statistik Transformasi Lyzenga Means Summary Report Class/Region

Band1

Band2

Band3

Band4

Band5

Band6

Band7

------

-----

-----

-----

-----

-----

-----

-----

1

124.5

109.25

90.75

20.25

13.25

11

41.75

2

115.667

98.667

73.333

21

16

11.333

42.333

3

117.667

97

69

17.333

18.667

16.333

38

4

116.9

100.9

82.6

19.6

13.6

12.2

41.5

5

137.167

126

92.667

19.833

13.167

11.833

46.333

6

144

137.222

114.852

20.667

17.37

15.37

53.296

7

139.071

122.25

91.071

18.179

18.036

15.929

44.357

8

138.86

127.035

96.368

17.018

15.649

14.211

45.825

9

139.4

133.9

113

19.6

14.2

12.6

52.3

10

117.917

103.667

82.417

19.083

18

15.917

39.833

11

121.214

107.857

78.643

15.5

15.5

13.429

37

12

124.875

105.875

76.75

20.5

23.625

20.5

39.625

13

121.222

101.111

78.778

26.667

35.667

30.333

44.667

14

119.667

98.444

73.889

24.556

32

27.778

41.556

15

115

99.6

72

18.2

20

16.6

36.6

16

128.4

114.6

83.6

16

14

12.8

40.8

17

119

109.5

96

21

14

13.5

41.5

18

128

120.333

94.667

16.667

13

12.333

42.667

19

125.667

115.5

95.833

18.167

16

13.5

44.5

20

128

113.5

93.25

21.75

19.75

17.25

40.5

21

124.5

105

77

19.5

22.5

20.5

40.5

22

127.333

104.667

66.333

18.333

21.667

19

35

23

127

101

65.5

19

22.5

20

36.5

24

135

111.5

71

18.5

22

18.5

39

25

123.333

102

69

15.667

17

15.333

37

26

157.167

142.833

97.533

15.867

15.5

14.2

47.3

27

122

106.75

87

16.5

15.75

14

42

28

127.8

114.6

98.4

21.6

14.2

12.4

47.4

29

134.778

110.556

67.556

17.667

20

17.556

39

30

149

125

72.5

17

20.5

17.5

38

31

130.286

105.571

65.571

16.714

18.143

15.571

36.143

86.416

56.363

47.596

32.969

39.689

27.899

34.285

Darat

79.76

64.339

60.87

81.076

97.941

57.781

62.769

Laut

87.307

52.466

41.621

16.797

19.527

17.336

24.268

All

  

Varian band 1 = 86.416 Varian band 2 = 56.363 Covar band 1 2= 105.409 else null a = (var b1-var b2)/(2*cov b1 b2) = 0.1426 ki/kj = a + √(a2+1) = 1.1527

Nilai batas darat laut adalah 23 Formula Lyzenga = if band 4 <=23 then log band1 + 1.1527 + log band 2

Lampiran 6. Peta tingkat resiko tsunami per kecamatan

113

RIWAYAT HIDUP Penulis dilahirkan di Sikka 24 Juli 1985. Penulis merupakan anak pertama dari tiga bersaudara pasangan Bapak Manduli Swa Sengaji, Bc.Ku. dan Ibu Nurhayati Bonda S.Pd.I. Pada tahun 2003 penulis diterima di Departemen Ilmu dan Teknologi Kelautan Institut Pertanian Bogor melalui jalur Seleksi Penerimaan Mahasiswa Baru (SPMB). Selama menjalani perkuliahan, penulis aktif dalam berbagai kegiatan kemahasiswaan yaitu staf Hubungan Luar dan Komunikasi Himpunan Mahasiswa Ilmu dan Teknologi Kelautan (Hublukom HIMITEKA) tahun 2004/2005 dan sebagai staf Pengembangan Sumber Daya Manusia HIMITEKA pada tahun 2005/2006. Pada tahun yang sama, penulis juga menjadi anggota Badan Eksekutif Mahasiswa Fakultas Perikanan dan Ilmu Kelautan IPB (BEM FPIK) dan anggota Himpunan Mahasiswa Ilmu dan Teknologi Kelautan seIndonesia (HIMITEKINDO). Pada bidang akademis, penulis merupakan asisten luar biasa mata kuliah Ekologi Perairan (2004/2005, 2005/2006 dan 2009) , asisten mata kuliah Ekologi Laut Tropis tahun 2007, asisten mata kuliah Oseanografi Umum tahun 2007/2008, asisten mata kuliah Sistem Informasi Geografis tahun 2008/2009. Untuk memperoleh gelar sarjana, penulis menyelesaikan tugas akhirdengan judul “Pemetaan Tingkat Resiko Tsunami di Kabupaten Sikka Nusa Tenggara Timur dengan Menggunakan Sistem Informasi Geografis”.