PEMETAAN BAHAYA TANAH LONGSOR DENGAN METODE FREQUENCY RATIO DI KECAMATAN PIYUNGAN DAN PLERET, KABUPATEN BANTUL, DAERAH ISTIMEWA YOGYAKARTA

PROCEEDING, SEMINAR NASIONAL KEBUMIAN KE-8 Academia-Industry Linkage 15-16 OKTOBER 2015; GRHA SABHA PRAMANA

PEMETAAN BAHAYA TANAH LONGSOR DENGAN METODE FREQUENCY RATIO DI KECAMATAN PIYUNGAN DAN PLERET, KABUPATEN BANTUL, DAERAH ISTIMEWA YOGYAKARTA Nusantara, Y.P.*, Setianto, A. Jurusan Teknik Geologi Fakultas Teknik Universitas Gadjahmada *corresponding author: [email protected]

ABSTRAK Tanah longsor merupakan salah satu bencana yang sering terjadi di Provinsi Yogyakarta. Tercatat telah terjadi 212 kejadian tanah longsor di Yogyakarta dari tahun 2011-2014 (BPBD DIY, 2014). Kejadian tanah longsor di Indonesia banyak menimbulkan kerugian baik dari segi fisik maupun sosial. Kejadian tanah longsor di daerah penelitian dipicu oleh curah hujan yang tinggi sehingga sangat mungkin bencana ini akan terus terjadi di masa depan. Pemetaan bahaya longsor dapat dilakukan dengan metode kualitatif dan metode kuantitatif. Salah satu metode kuantitatif yang digunakan adalah frequency ratio. Frequency ratio didasarkan kepada hubungan antara lokasi kejadian tanah longsor dan faktor-faktor yang mempengaruhi terjadinya longsor. Metode ini diaplikasikan dengan menggunakan GIS dan data penginderaan jauh seperti citra DEM Terrasar-X dan citra Google Earth serta peta geologi 1:100.000 dan peta topografi 1:25.000. Data tersebut merepresentasikan faktor-faktor yang mempengaruhi longsor yaitu kemiringan lereng, jarak dari drainase, jarak dari patahan, satuan litologi, tata guna lahan ,dan presipitasi. Faktor-faktor tersebut kemudian dikonstuksikan menjadi peta dengan format raster 10x10m dan ditumpangtindih sehingga menghasilkan nilai frequency ratio. Nilai frequency ratio dikalkulasikan menjadi LSI (Landslide Susceptibility Index) yang kemudian diklasifikan menjadi 4 zona bahaya longsor yaitu zona bahaya rendah, zona bahaya sedang, zona bahaya tinggi dan zona bahaya sangat tinggi yang kemudian disajikan dalam bentuk peta bahaya tanah longsor. Hasil dari metode ini menunjukkan bahwa kemiringan lereng merupakan faktor yang berperan penting menyebabkan terjadinya longsor di daerah penelitian. Hasil validasi menunjukkan bahwa metode ini menghasilkan peta bahaya longsor yang cukup memuaskan dengan nilai akurasi sekitar 70% sehingga dapat diterapkan dalam upaya mitigasi bencana longsor.

I.

kerugian (BNPB,2014). Daerah yang terjadi tanah longsor kebanyakan terjadi di daerah dengan lereng yang terjal dan dipicu oleh curah hujan yang tinggi terutama di daerah Kecamatan Piyungan dan Pleret yang menjadi daerah penelitian yang ditunjukkan Gambar 1. Tanah longsor di daerah ini banyak menimbulkan kerugian yaitu rusaknya jalan dan rumah warga. Dalam pemetaan bahaya tanah longsor terdapat 2 metode yaitu metode pemetaan langsung dan metode pemetaan tidak langsung. Dalam penelitian ini menggunakan metode pemetaan tidak langsung yaitu metode frequency ratio. Metode frequency ratio merupakan metode yang dibangun berdasarkan hubungan antara lokasi kejadian tanah longsor dan faktor-faktor yang mengontrol terjadinya tanah longsor (Lee dan

PENDAHULUAN Tanah longsor merupakan bencana alam yang paling banyak menimbulkan kerugian baik materi maupun korban jiwa yang terus menjadi pusat perhatian karena bencana ini terjadi secara kontinu dari tahun ke tahun. Beberapa peneliti telah melakukan penelitian untuk mengembangkan peta bahaya tanah longsor yang dapat digunakan untuk upaya mitigasi bencana tersebut. Peta bahaya tanah longsor memiliki fungsi untuk membagi area bahaya longsor menjadi beberapa zona yang berbeda tergantung tingkat bahaya area itu sendiri. Daerah Istimewa Yogyakarta merupakan daerah di Indonesia yang sering terjadi tanah longsor. Tercatat dari tahun 2004 sampai tahun 2014 telah terjadi 212 kasus tanah longsor dan telah menimbulkan banyak 513

PROCEEDING, SEMINAR NASIONAL KEBUMIAN KE-8 Academia-Industry Linkage 15-16 OKTOBER 2015; GRHA SABHA PRAMANA Pradhan, 2006). Faktor-faktor yang mengontrol terjadinya tanah longsor yaitu kemiringan lereng, jarak dari sesar, jarak dari drainase, litologi, tata guna lahan dan presipitasi. Faktor-faktor tersebut dikonstruksikan dan ditumpangtindihkan menggunakan GIS sehingga menghasilkan peta bahaya tanah longsor.

II.

dikonstruksikan menggunakan GIS (Geographic Information System) seperti yang ditunjukkan Tabel 1. Lokasi kejadian tanah longsor ditentukan berdasarkan laporan peneliti terdahulu dan pengamatan lapangan secara langsung. Pada daerah penelitian, total telah terjadi longsor sebanyak 31 kejadian yang tersebar di seluruh area. Faktor-faktor yang mengontrol tanah longsor terdiri dari 6 faktor yaitu kemiringan lereng, jarak dari sesar, jarak dari drainase, litologi, tata guna lahan dan presipitasi. Faktor-faktor tersebut telah ditransformasikan menjadi data vektor dan raster menggunakan GIS. Faktor kemiringan lereng lereng diketahui dari citra DEM TerraSAR-X dengan resolusi 10 m. Kemiringan lereng di daerah penelitian bervariasi dari 00-440 yang kemudian dibagi menjadi 3 zona berdasarkan yaitu zona 00-200, zona 200-400 dan zona >400 (Karnawati,2005) seperti yang ditunjukkan Gambar 2. Faktor Tata guna lahan diketahui dari interpretasi citra satelit Bingmaps dan dikoreksi dengan peta RBI 1:25.000 Lembar Lembar Timoho, Jabung, dan Imogiri. Tata guna lahan di daerah penelitian terdiri dari pemukiman, persawahan, perkebunan, tanah kosong, Tempat Pembuangan Sampah (TPS), hutan dan semak belukar seperti Gambar 3. Jarak dari sesar dan litologi diketahui dari peta geologi 1:100.000 lembar Yogyakarta dan Surakarta dan hasil pemetaan di lapangan. Daerah penelitian dibagi menjadi 3 satuan yaitu satuan perselingan batupasir tuffan dan tuff dengan sisipan breksi autoklastik, satuan perselingan batupasit tuffan dengan tuff dengan sisipan lapili tuff, dan satuan breksi andesit seperti Gambar 4. Jarak dari sesar di daerah penelitian dibagi menjadi 4 zona yaitu zona 0-100 m, zona 100-200 m, zona 200-300 m dan zona >300 m seperti Gambar 5. Faktor jarak dari drainase diketahui dari Peta Pola Penyaluran. Jarak dari drainase dibagi menjadi 4 zona yaitu zona 0-100 m, zona 100-200 m, zona 200-300 m, dan zona >300 m seperti Gambar 6. Faktor presipitasi diketahui dari data curah hujan selama 10 tahun terakhir

KONDISI GEOLOGI REGIONAL

Secara fisiografi regional, daerah penelitian termasuk ke dalam Zona Pegunungan Selatan Bagian Timur yang terletak di selatan Jawa Tengah dan Jawa Timur, memanjang ke arah timur dari bagian tenggara Provinsi Daerah Istimewa Yogyakarta sampai Pantai Selatan Jawa Timur (Van Bemmelen, 1949). Secara stratigrafi, daerah penelitian berada di stratigrafi Pegunungan Selatan. Adapun urutan stratigrafi dari yang tertua sampai termuda Pegunungan Selatan yang tersingkap di daerah penelitian yaitu (Surono, 2009) :  Formasi Semilir Formasi ini terdiri dari batupasir tuffan, tuff lapili, tuff lempungan dan tuff yang berumur Miosen Awal.  Formasi Nglanggran Formasi ini terdiri dari breksi gunung api dan aglomerat, dengan sisipan tuf dan lava andesit yang berumur Miosen Awal. Struktur geologi yang terdapat di daerah penelitian merupakan struktur dengan arah baratlaut-tenggara atau biasa disebut dengan Pola Sumatra dan struktur dengan arah timurlaut-baratdaya atau disebut dengan Pola Meratus (Surono, 2009). Struktur geologi yang ditemukan berupa sesar turun, sesar geser dekstral dan sesar geser sinistral.

III.

DATA DAN PENELITIAN

METODE

Dalam penelitian ini, spatial database yang menjelaskan kejadian longsor dengan faktorfaktor pengontrol tanah longsor 514

PROCEEDING, SEMINAR NASIONAL KEBUMIAN KE-8 Academia-Industry Linkage 15-16 OKTOBER 2015; GRHA SABHA PRAMANA Badan Meteorologi, Klimatologi, dan Geofisika Yogyakarta. Curah hujan di daerah penelitian bervariasi dari 1699-1709 mm/tahun yang kemudian dibagi menjadi 9 zona seperti Gambar 7.

dengan masing-masing faktor pengontrol longsor. Hubungan tersebut menjelaskan kontribusi setiap faktor terhadap terjadinya tanah longsor. Bila nilai frequency ratio lebih dari 1 maka Nilai frequency ratio masingmasing zona faktor pengontrol longsor ditunjukkan oleh Tabel 2. Nilai frequency ratio setiap zona di daerah penelitian bervariasi dari 0 - 6,19. Dalam faktor kemiringan lereng, lereng dengan kemiringan yang lebih terjal memiliki probabilitas lebih tinggi untuk terjadinya tanah longsor. Kemiringan lereng di bawah 200 memiliki nilai frequency ratio 0,50 yang mengindikasikan probabilitas untuk terjadinya longsor lebih rendah. Kemiringan lereng di 200-400 memiliki nilai frequency ratio 2,83 yang mengindikasikan probabilitas untuk terjadinya longsor lebih tinggi. Semakin tinggi kemiringan lereng, meningkatkan shear stress batuan penyusun lereng sehingga kestabilan lereng terganggu dan berpotensi terjadinya tanah longsor (Pradhan, 2010). Untuk menghitung nilai LHI, setiap nilai frequency ratio masing-masing zona yang bertampalan dihitung dengan rumus (2). Nilai LHI yang terdapat di area penelitian bervariasi dari 2,77 – 18,22. Nilai LHI ini kemudian diklasifikasikan menjadi 3 zona yaitu zona bahaya rendah (22%), zona bahaya menengah (65%0dan zona bahaya tinggi (13%) yang ditunjukkan oleh Gambar 8.

Faktor-faktor kemiringan longsor ini kemudian ditransformasikan menjadi data raster dengan resolusi 10 m seperti yang ditunjukkan Tabel 1. Faktor pengontrol longsor tersebut ditumpangtindihkan dengan titik kejadian longsor. Setiap zona dari faktor memiliki 2 atribut yaitu jumlah piksel area dalam zona dan jumlah longsor yang terdapat di dalam zona. Kedua atribut ini kemudian digunakan untuk menghitung nilai frequency ratio masing-masing zona dengan rumus (1) (Bonham, 1994 dalam Tazik dkk, 2014). 𝐷𝑖/𝐴𝑖 ∑𝑁 𝐷𝑖/ 𝑖=1 𝑖=1 𝐴𝑖

𝐹𝑅 = ∑𝑁

FR = Frequency Ratio Di = jumlah titik longsor pada suatu kelas dalam faktor tertentu Ai = jumlah area pada suatu kelas dalam faktor tertentu Setelah nilai FR masing-masing zona faktor pengontrol longsor dihitung, kemudian masing-masing peta ditumpang-tindihkan sehingga menghasilkan nilai LHI (Landslide Hazard Index) seperti rumus (2) (Pradhan,2010). LHI = FR1 + FR2 + FR3 + ... + FRn

Hasil dari analisis bahaya longsor kemudian divalidasikan menggunakan data tanah longsor yang telah ada. Nilai LHI yang didapatkan kemudian dibagi menjadi 10 kelas dari yang tertinggi sampai terendah (100%-90%, 9080%,..., 10%-0%). Setiap kelas LHI dihitung persentase jumlah longsor dan persentase komulatifnya seperti Kurva perubahan laju dibuat dengan Nilai LHI sebagai sumbu x dan persentase komulatif kejadian longsor sebagai sumbu y seperti Gambar 9. AUC ( Area Under Curve) dihitung dari kurva perubahan laju tersebut yang mengindikasikan nilai akurasi prediksi dari model frequency ratio. Nilai AUC berkisar 0,705 yang berarti nilai akurasi prediksi bernilai 70,5%. Nilai AUC dari kurva

(2)

Kemudian nilai LHI ini diklasifikasikan menjadi 3 zona yaitu zona tingkat bahaya rendah, zona tingkat bahaya menengah dan zona tingkat bahaya tinggi yang disajikan dalam bentuk peta bahaya tanah longsor. Peta bahaya tanah longsor divalidasikan dengan titik kejadian tanah longsor. Validasi ini menunjukkan seberapa baik model dalam memprediksi longsor. Hasil dari validasi ini akan menghasilkan nilai akurasi prediksi berdasarkan AUC (Area Under Curve).

IV.

ANALISIS DATA DAN DISKUSI

Metode frequency ratio didasarkan kepada hubungan antara penyebaran titik longsor 515

PROCEEDING, SEMINAR NASIONAL KEBUMIAN KE-8 Academia-Industry Linkage 15-16 OKTOBER 2015; GRHA SABHA PRAMANA tersebut >0,5 yang menyatakan bahwa model frequency ratio memiliki nilai akurasi prediksi yang baik dan dapat diterapkan di daerah penelitian.

V.

tanah yang paling berpengaruh di daerah penelitian adalah kemiringan lereng 200-400 dengan nilai FR 2,83 dengan pelamparan yang luas dibandingkan faktor-faktor pengontrol longsor yang lain.

KESIMPULAN

Pemetaan bahaya tanah longsor di daerah Kecamatan Piyungan dan Pleret, Kabupaten Bantul, Daerah Istimewa Yogyakarta dengan metode frequency ratio menghasilkan 3 zona bahaya yaitu zona bahaya rendah dengan pelamparan sekitar 22%, zona bahaya menengah dengan pelamparan sekitar 65% dan zona bahaya tinggi dengan pelamparan sekitar 13%. Validasi menghasilkan frequency ratio memiliki akurasi prediksi yang baik dengan nilai akurasi 70,5%. Faktor pengontrol

VI.

ACKNOWLEDGEMENT

Ucapan terima kasih disampaikan kepada : 1. BPBD Daerah Istimewa Yogyakarta atas bantuannya dan ijin pengambilan data kejadian tanah longsor daerah penelitian. 2. BMKG Daerah Istimewa Yogyakarta atas batuannya dan ijin pengambilan data curah hujan daerah penelitian.

DAFTAR PUSTAKA Badan Penanggulangan Bencana Daerah (BPBD) Daerah Istimewa Yogyakarta, 2013, Peta Bahaya Tanah Longsor Daerah Istimewa Yogyakarta, skala 1:250.000, Pemerintah Daerah Istimewa Yogyakarta. Bonham-Carter G.F., 1994, Geographic Information System for Geoscientists, Modelling with GIS, Pergamon, Oxford, 398 p. Karnawati, D., 2005, Bencana Alam Gerakan Massa Tanah di Indonesia dan Upaya Penanggulangannya, Jurusan Teknik Geologi Fakultas Teknik Universitas Gadjah Mada, Yogyakarta. Lee, S., Pradhan, B., 2007, Landslide Hazard Mapping at Selangor, Malaysia Using Frequency Ratio and Logistic Regression Models, Landslide Vol.4, p. 33-41 Margono, U., Surono, dan Kusnama, 2009 (dalam persiapan). Peta Geologi Lembar Wonosari dan Semanu, Yogyakarta, Skala 1:50.000. Pusat Survei Geologi, Bandung Pradhan, B., 2010, Landslide Susceptibility Mapping of a Catchment Area Using Frequemcy Ratio, Fuzzy Logic and Multivariate Logistic Regression Approaches, Journal Indian Society Remote Sensing Vol. 38, 301-320 p. Rahardjo, W., Sukandarrumidi, Rosidi, H.M.D., 1995, Peta Geologi Lembar Yogyakarta, Jawa, skala 1:100.000. Pusat Penelitian dan Pengembangan Geologi, Bandung Sabins, F.F., 1987. Remote Sensing : Principles and Interpretation, Freeman and Company, New York, 449 p. Surono, 2008, Sedimentasi Formasi Semilir di Desa Sendang, Wuryanto, Wonogiri, Jawa Tengah, Jurnal Sumber Daya Geologi Vol.18, p. 29-41 Surono, 2008, Stratigrafi dan Sedimentasi Formasi Kebo dan Formasi Butak di Pegunungan Selatan, Jawa Bagian Selatan. Jurnal Geologi Indonesia Vol.3, p. 183-193 Surono, 2009, Litostratigrafi Pegunungan Selatan Bagian Timur Daerah Istimewa Yogyakarta dan Jawa Tengah, Jurnal Sumber Daya Geologi Vol. 19, p. 31-43 516

PROCEEDING, SEMINAR NASIONAL KEBUMIAN KE-8 Academia-Industry Linkage 15-16 OKTOBER 2015; GRHA SABHA PRAMANA Surono, Toha, B., Sudarno, 1992, Peta Geologi Lembar Surakarta-Giritontro, Jawa, Skala 1:100.000. Pusat Penelitian dan Pengembangan Geologi, Bandung. Sutanto, 1986, Penginderaan Jauh Jilid I, Gadjah Mada University Press, Yogyakarta Tazik, E., Jahantab, Z., Bakhtiari, M., Rezaei, A., Alavipanah, K.S., 2014, Landslide Susceptibility Mapping by Combining the Three Methods Fuzzy Logic, Frequency Ratio and Analytical Hierarchy Process in Dozain Basin, The International Arhives of the Photogrammetry, Remote Sensing and Spatial Information Sciences Vol. 40, p. 267-272 Varnes, D.J., 1978, Slope Movement Types and Processes, Special Report 176; Landslides; Analysis and Control, Eds: R.L. Schuster dan R.J. Krizek, Transport Research Board, National Research Council, Washington, D.C, p 11-33 Bemmelen, Van R.W., 1949, The Geology of Indonesia, Govt. Printing Office, Nijhoff, The Hague, Jawa Tengah

TABEL Tabel 1. Spatial database daerah penelitian beserta sumbernya. Sumber Informasi Tipe data Citra DEM TerraSar-X

Kemiringan lereng

Citra IKONOS dan Peta RBI

Tata guna lahan

Laporan Peneliti Terdahulu

Titik kejadian longsor

Data Curah Hujan

Presipitasi

Peta Pola Penyaluran Peta Geologi dan Pemetaan

Skala

Raster

1 : 50.000

Vektor

1 : 50.000

Vektor

1 : 50.000

Raster

1 : 50.000

Jarak dari drainase

Vektor

1 : 50.000

Litologi Jarak dari sesar

Vektor

1 : 50.000

Tabel 2. Faktor pengontrol longsor beserta nilai FR (frequency ratio) masing-masing zona Faktor Kemiringan Lereng

Litologi

Tata Guna Lahan

Zona

Di

∑ Di

Ai

∑ Ai

FR

0°-20° 20°-40° >40° Perselingan batupasir tuffan dan tuff dengan sisipan breksi autoklastik Perselingan batupasir tuffan dan tuff dengan sisipan lapili tuff Breksi Andesit

259169 72086 1133

332318 332318 332318

12 19 0

31 31 31

0,50 2,83 0

69530

332388

12

31

1,85

256913 5875

332388 332388

16 3

31 31

0,67 5,47

Tanah kosong Hutan

16 3132

330745 330745

0 0

31 31

0 0

517


Jarak dari drainase

Jarak dari sesar

Presipitasi

Tempat Pembuangan Sampah Pemukiman Persawahan Ladang/Tegalan Perkebunan Semak belukar

1304 45575 91470 155598 12981 20669

330745 330745 330745 330745 330745 330745

0 0 2 12 5 12

31 31 31 31 31 31

0 0 0,23 0,82 4,11 6,19

0-100 100-200 200-300 >300 0-100 100-200 200-300 >300 1699-1709 1709-1719 1719-1729 1729-1739 1739-1749 1749-1759 1759-1769 1769-1779 1779-1790

152857 94695 46456 38328 56082 53003 47896 175356 88867 85599 25604 22352 24313 29851 28553 27006 192

332336 332336 332336 332336 332337 332337 332337 332337 332337 332337 332337 332337 332337 332337 332337 332337 332337

12 8 8 3 3 5 6 17 7 2 1 2 4 6 7 2 0

31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31

0,84 0,91 1,85 0,84 0,57 1,01 1,34 1,04 0,84 0,25 0,42 0,96 1,76 2,15 2,63 0,79 0

GAMBAR

Gambar 1. Peta lokasi daerah penelitian. (tidak dalam skala sebenarnya)

518


Gambar 2. Peta geologi daerah penelitian (tidak dalam skala sebenarnya).

Gambar 3. Peta tata guna lahan daerah penelitian (tidak dalam skala sebenarnya).

519


Gambar 4. Peta geologi daerah penelitian (tidak dalam skala sebenarnya)

Gambar 5. Peta jarak dari sesar daerah penelitian (tidak dalam skala sebenarnya).

520


Gambar 6. Peta jarak dari drainase daerah penelitian (tidak dalam skala sebenarnya).

Gambar 7. Peta curah hujan daerah penelitian (tidak dalam skala sebenarnya). 521


Gambar 8. Peta bahaya tanah longsor daerah penelitian (tidak dalam skala sebenarnya)

Gambar 9. Kurva perubahan laju dari peta bahaya tanah longsor dengan metode frequency ratio Nilai AUC (Area Under Curve) berkisar 0,705.

522

PEMETAAN BAHAYA TANAH LONGSOR DENGAN METODE FREQUENCY RATIO DI KECAMATAN PIYUNGAN DAN PLERET, KABUPATEN BANTUL, DAERAH ISTIMEWA YOGYAKARTA

Recommend Documents