TESIS-RG092999
ANALISA PENURUNANAN TANAH (LANDSUBSIDENCE) PADA DAERAH SEMBURAN LUMPUR SIDOARJO DENGAN DATA SATELIT GLOBAL POSITIONING SYSTEM (GPS)
WISNU PRIBADI 3512201006 Dosen Pembimbing: Dr. Ing. Ir. Teguh Hariyanto, MSc Ira Mutiara Anjasmara, ST. M.Phil. Ph.D
PROGRAM MAGISTER BIDANG KEAHLIAN MANAJEMEN SURVEY DAN PEMETAAN JURUSAN TEKNIK GEOMATIKA FAKULTAS TEKNIK SIPIL DAN PERENCANAAN INSTITUT TEKNOLOGI SEPULUH NOPEMBER SURABAYA 2014
THESIS-RG092999
ANALYSIS OF LAND SUBSIDENCE ON SIDOARJO MUD FLOW USING GLOBAL POSITIONING SYSTEM SATELITE DATA
WISNU PRIBADI 3512201006 Supervisor: Dr. Ing. Ir. Teguh Hariyanto, MSc Ira Mutiara Anjasmara, ST. M.Phil. Ph.D
MASTER PROGRAM SURVEYING AND MAPPING MANAGEMENT EXPERTISE GEOMATIC ENGINEERING DEPARTMENT FACULTY OF CIVIL ENGINEERING AND PLANNING SEPULUH NOPEMBER INSTITUTE OF TECHNOLOGY SURABAYA 2014
ANALISA PENURUNANAN TANAH (LANDSUBSIDENCE) PADA DAERAH SEMBURAN LUMPUR SIDOARJO DENGAN MENGGUNAKAN DATA GLOBAL POSITIONING SYSTEM (GPS) Nama NRP Pembimbing Co-pembimbing
: : : :
Wisnu Pribadi 3512 201 006 Dr.Ing.Ir.Teguh Hariyanto, M.Sc Ira Mutiara Anjasmara, ST. M.Phil. Ph.D
ABSTRAK Bencana lumpur di Sidoarjo yang terjadi di Kecamatan Porong sejak tahun 2006 lalu mengakibatkan dampak yang luar biasa bagi masyarakat sekitar maupun bagi aktivitas perekonomian di Jawa Timur. Hal ini dikarenakan akibat yang ditimbulkan oleh semburan lumpur menyebabkan beberapa perubahan pada formasi tanah. Perubahan yang terjadi pada permukaan dapat dikategorikan menjadi dua, yaitu perubahan secara horizontal dan perubahan secara vertikal. Salah satu fenomena deformasi permukaan bumi secara vertikal adalah penurunan tanah (landsubsidence) dan penaikkan tanah (uplift). Pada penelitian ini akan dilakukan analisa mengenai penurunan tanah yang diakibatkan oleh keluarnya air tanah bersama dengan lumpur dan gas di daerah kecamatan Porong, Kabupaten Sidoarjo, dengan menggunakan metode GPS untuk mendapatkan nilai penurunan tanah di Sidoarjo. Titik – titik kontrol yang diamati ada 20 titik yang tersebar di kecamatan Porong, Kabupaten Sidoarjo. Pengamatan dilakukan 3 kala yaitu pada bulan Maret, April dan Mei 2014. Hasil dari pengamatan titik yang mengalami penurunan terendah adalah sebesar -0.091 m, yaitu titik BT03 pada perubahan kala 2 dan kala 3, serta terjadi penaikan tanah tertinggi sebesar 0.171 m yaitu pada titik VK23 pada perubahan kala 1 dan kala 2. Dapat disimpulkan bahwa tanah disekitar semburan lumpur sangat tidak stabil, karena terjadi perubahan signifikan pada titik-titik yang berada dekat di sekitar lumpur. Kata kunci
: Penurunan tanah, GPS, kabupaten Sidoarjo
iii
“Halaman ini sengaja dikosongkan”
iv
ANALYSIS OF LAND SUBSIDENCE ON SIDOARJO MUD FLOW USING GLOBAL POSITIONING SYSTEM SATELITE DATA Name NRP Supervisor Co-supervisor
: : : :
Wisnu Pribadi 3512 201 006 Dr.Ing.Ir.Teguh Hariyanto, M.Sc Ira Mutiara Anjasmara, ST. M.Phil. Ph.D
ABSTRACT The Sidoarjo mud disaster that occurred in Porong subdistrict since 2006 ago gave a tremendous impact for human life and economic activity in East Java. It caused by some changes in the land formations. The changes on the surface can be categorized into two models, horizontal deformation and vertical deformation. One of them is called landsubsidence and uplift surface in vertical deformation. This research will be conducted analysis of landsubsidence caused by groundwater and gas turn out to be mud vulcano on the surface around Porong subdistrict, Sidoarjo regency, using GPS method to get subsidence value. The control points that observed was spread in 20 point area in Porong, Sidoarjo Observation will be doing in 3 times, March, April and May 2014. The result of observation declined the lowest point is -0,091 m at BT03 point when taken between on stage 2 and stage 3. As well as the highest uplift occurs at 0.171 m when taken between on stage 1 and stage 2. We could concluded that the soil around mudflow is very unstable and critical damaged, so the significant of changes at the points that are close around the mud.
Keyword: Landsubsidence, GPS, and Sidoarjo
v
“Halaman ini sengaja dikosongkan”
vi
KATA PENGANTAR
Puji syukur kehadirat Tuhan Yang Maha Esa atas limpahan rahmat, taufik dan hidayah-Nya sehingga penulis dapat menyelesaikan penyusunan tesis yang berjudul Analisa Penurunan Tanah (Landsubsidence) Pada Daerah semburan Lumpur Sidoarjo Dengan Menggunakan Data Satelit Global Positioning System (GPS) sebagai syarat kelulusan program Magister. Selama penyelesaian tesis ini penulis banyak mendapatkan bantuan dari berbagai pihak yang turut memberikan kontribusi dalam bentuk saran, kritik dan doa, untuk itu penulis mengucapkan banyak terima kasih kepada : 1. Tuhan Yang Maha Esa atas limpahan Rahmat dan hidayah-Nya. 2. Orang tua yang selalu memotivasi, mendoakan, serta membiayai kami selama penelitian ini. 3. Dr-Ing. Ir. Teguh Hariyanto, M.Sc sebagai dosen pembimbing I atas segala bantuan dan nasihatnya. 4. Ira Mutiara Anjasmara, ST, M.Phil, Ph.D sebagai dosen pembimbing II atas segala bantuan dan nasihatnya. 5. Seluruh dosen pengajar dan staff tata usaha Jurusan Teknik Geomatika ITS atas segala bantuannya. 6. Para sahabat di kampus beserta semua teman-teman S1 Teknik Geomatika ITS. 7. Semua pihak yang telah membantu dan mendukung selama pelaksanaan penelitian dan penyusunan laporan tesis Penulis menyadari dalam laporan tesis
ini masih terdapat banyak
kekurangan, baik dari segi penulisan maupun isi laporan. Oleh karena itu, penulis sangat mengharapkan kritik dan saran yang membangun untuk perbaikan di masa depan. Semoga laporan tesis ini dapat memberikan manfaat pada kalangan mahasiswa dan pembaca.
Surabaya, Juli 2014
vii
(Halaman ini sengaja dikosongkan)
viii
DAFTAR ISI HALAMAN JUDUL........................................................................................ i LEMBAR PENGESAHAN ............................................................................. ii ABSTRAK ....................................................................................................... iii KATA PENGANTAR ..................................................................................... v DAFTAR ISI .................................................................................................... vii DAFTAR GAMBAR ....................................................................................... ix DAFTAR TABEL ............................................................................................ xi DAFTAR LAMPIRAN .................................................................................... xiii
BAB 1 PENDAHULUAN .............................................................................. 1 1.1 Latar Belakang ...................................................................................... 1 1.2 Rumusan dan Batasan Masalah ............................................................. 2 1.3 Tujuan Penelitian ................................................................................... 3 1.4 Manfaat Penelitian ................................................................................. 3
BAB II KAJIAN PUSTAKA DAN DASAR TEORI................................... 5 2.1 Penurunan Tanah ................................................................................... 5 2.2 Teknik Pemantauan Land Subsidence Dengan GPS ............................. 7 2.3 Global Positioning System (GPS) ......................................................... 8 2.3.1 Pengamatan GPS ......................................................................... 9 2.3.2 Penentuan Posisi dengan GPS ..................................................... 10 2.3.3 Penentuan Tinggi dengan GPS .................................................... 11 2.3.4 Kesalahan dan Bias ..................................................................... 12 2.4 Metode Radial ...................................................................................... 15 2.5 Analisa Regresi...................................................................................... 15 2.5.1 Regresi Linier .............................................................................. 16 2.5.2 Regresi Polinomial ...................................................................... 18 2.6 Koefisien Korelasi ................................................................................ 19 2.7 Koefisien Determinasi ........................................................................... 19
vii
2.8 Uji Statistik Student-t ............................................................................ 20 2.9 Lumpur Sidoarjo .................................................................................. 21 2.9.1 Latar Belakang Lumpur .............................................................. 21 2.9.2 Geomorfologi Sidoarjo ................................................................ 21 2.10 Studi Hasil Penelitian Sebelumnya .................................................... 22
BAB III METODA PENELITIAN ............................................................... 25 3.1 Lokasi Penelitian .................................................................................. 25 3.2 Data dan Peralatan ................................................................................. 27 3.2.1 Data ............................................................................................. 26 3.2.2 Peralatan ...................................................................................... 27 3.3 Metodologi Penelitian ........................................................................... 28 3.3.1 Tahap Persiapan .......................................................................... 29 3.3.2 Tahap Pengolahan Data .............................................................. 29 3.3.3 Tahap Analisa Data ..................................................................... 30
BAB IV HASIL PENELITIAN DAN PEMBAHASAN.............................. 31 4.1 Hasil....................................................................................................... 31 4.2 Uji T-Student ......................................................................................... 36 4.3 Analisa Regresi pada Setiap Pengukuran .............................................. 37
BAB V KESIMPULAN DAN SARAN ......................................................... 61 5.1 Kesimpulan ............................................................................................ 61 5.2 Saran ...................................................................................................... 62
DAFTAR PUSTAKA LAMPIRAN BIOGRAFI
viii
DAFTAR GAMBAR Gambar 2.1 Sistem tinggi................................................................................. 11 Gambar 2.2 Kesalahan ephemeris .................................................................... 13 Gambar 2.3 Prinsip Dasar Penentuan Posisi Satelit GPS ................................ 14 Gambar 2.4 Metode pengamatan radial ........................................................... 15 Gambar 2.5 Grafik untuk b positif ................................................................... 17 Gambar 2.6 Grafik untuk b negatif .................................................................. 17 Gambar 2.7 Citra satelit busur vulkanik selatan Pulau Jawa ........................... 22 Gambar 3.1 Daerah pengamatan dalam penelitian .......................................... 25 Gambar 3.2 Diagram Tahap Pengolahan Data ................................................. 28 Gambar 4.1. Sebaran titik pengamatan ............................................................ 31 Gambar 4.2 Grafik regresi di titik BPNRI-53 .................................................. 37 Gambar 4.3 Grafik regresi di titik BT01 .......................................................... 38 Gambar 4.4 Grafik regresi di titik BT03 .......................................................... 39 Gambar 4.5 Grafik regresi di titik BW13 ........................................................ 40 Gambar 4.6 Grafik regresi di titik DG09 ......................................................... 41 Gambar 4.7 Grafik regresi di titik DG11 ......................................................... 42 Gambar 4.8 Grafik regresi di titik DG18 ......................................................... 43 Gambar 4.9 Grafik regresi di titik DG24 ......................................................... 44 Gambar 4.10 Grafik regresi di titik DG34 ....................................................... 45 Gambar 4.11 Grafik regresi di titik DGIJK ..................................................... 46 Gambar 4.12 Grafik regresi di titik PLNG....................................................... 47 Gambar 4.13 Grafik regresi di titik PT06 ........................................................ 48 Gambar 4.14 Grafik regresi di titik PT11 ........................................................ 49 Gambar 4.15 Grafik regresi di titik TTG1304 ................................................. 50 Gambar 4.16 Grafik regresi di titik TTG1305 ................................................. 51 Gambar 4.17 Grafik regresi di titik TTG1307 ................................................. 52 Gambar 4.18 Grafik regresi di titik VK09 ....................................................... 53 Gambar 4.19 Grafik regresi di titik VK13 ....................................................... 54 Gambar 4.20 Grafik regresi di titik VK14 ....................................................... 55
ix
Gambar 4.21 Grafik regresi di titik VK23 ....................................................... 56 Gambar 4.22 Hasil sebaran titik beserta keadaan tanaah selama pengamatan 57 Gambar 4.23 Peta daerah yang mengalami uplift dan subsidence berdasarkan penelitian terdahulu .......................................................................................... 58
x
“Halaman ini sengaja dikosongkan”
xi
DAFTAR TABEL Tabel 3.1. Koordinat sebaran titik-titik kontrol pengamatan ........................... 26 Tabel 4.1 hasil pengolahan data kala 1 bulan Maret 2014 ............................... 32 Tabel 4.2 hasil pengolahan data kala 2 bulan April 2014 ................................ 32 Tabel 4.3 hasil pengolahan data kala 3 bulan Mei 2014 .................................. 33 Tabel 4.4 Nilai Subsidence Maupun Uplift Setiap Titik .................................. 35 Tabel 4.5 Tabel Standar deviasi hasil perhitungan uji T-student ..................... 36 Tabel 4.6 Perbandingan dengan penelitian terdahulu ...................................... 59
xi
“Halaman ini sengaja dikosongkan”
xii
DAFTAR LAMPIRAN Raw Data Pengukuran GPS di lapangan.......................................................... 67 Grafik pengamatan GPS................................................................................... 71 Foto Dokumentasi ............................................................................................ 78
xiii
“Halaman ini sengaja dikosongkan”
xiv
BAB 1 PENDAHULUAN 1.1 Latar Belakang Bencana lumpur yang terjadi di Kecamatan Porong sejak tahun 2006 lalu mengakibatkan dampak yang luar biasa bagi masyarakat sekitar maupun bagi aktivitas perekonomian di Jawa Timur. Peristiwa ini berawal ketika banjir lumpur panas mulai menggenangi areal persawahan, pemukiman penduduk dan kawasan industri. Volume lumpur yang keluar dari titik sumber bencana itu diperkirakan sekitar 5.000 hingga 50 ribu meter kubik perhari (setara dengan muatan penuh 690 truk peti kemas berukuran besar). Semburan lumpur tersebut pertama kali terjadi pada tanggal 27 Mei 2006 di dekat sumur bor Banjar Panji milik PT Lapindo Brantas Inc., di Kecamatan Porong, Kabupaten Sidoarjo. Kemudian pada tanggal 1 Juni 2006, muncul lagi sumber baru di rumah penduduk Balong Nongo. Beberapa sumber telah berhenti, namun ada satu yang masih menyemburkan lumpur sampai sekarang yaitu di dekat sumur bor. (Erif, 2011) Perubahan yang terjadi pada permukaan tanah (deformasi) dapat dikategorikan menjadi dua, yaitu perubahan secara horizontal dan perubahan secara vertikal. Fenomena deformasi permukaan bumi secara vertikal adalah penurunan tanah (landsubsidence) dan penaikan tanah (uplift). Penurunan tanah dapat terjadi baik secara lokal maupun regional. Kondisi tersebut dikarenakan oleh beberapa faktor penyebab, antara lain sebagai berikut, penurunan muka tanah alami (natural subsidence) yang disebabkan oleh proses-proses geologi seperti aktifitas vulkanik dan tektonik, siklus geologi, adanya rongga di bawah permukaan tanah dan sebagainya. Penurunan muka tanah yang disebabkan oleh pengambilan bahan cair dari dalam tanah seperti air tanah. Penurunan muka tanah yang disebabkan oleh adanya beban-beban berat diatasnya seperti struktur bangunan. Penurunan muka tanah akibat pengambilan bahan padat dari tanah (aktifitas penambangan). Pada penelitian ini tinjauan terhadap faktor penyebab penurunan tanah dilakukan terhadap perubahan muka air tanah, kondisi tutupan
1
lahan terkait dengan beban bangunan dan tingkat konsolidasi yang terjadi. (Whittaker and Reddish, 1989) Salah satu dampak yang ditimbulkan luapan lumpur Sidoarjo adalah terjadinya land subsidence (penurunan tanah). Land
subsidence sendiri
didefinisikan sebagai penurunan muka tanah sebagai fungsi dari waktu, atau dapat juga disebut deformasi vertikal, yang artinya perubahan kedudukan vertikal dari permukaan tanah terhadap bidang referensi yang dianggap tetap (Hamdani, 2004). Adapun penyebab land subsidence adalah pengambilan air tanah yang berlebihan (groundwater over exploitation), keluarnya gas dan lumpur secara massive dari dalam tanah ke permukaan, penurunan karena beban bangunan, penurunan karena adanya konsolidasi alamiah dari lapisanlapisan tanah, serta penurunan karena gaya-gaya tektonik (Cahyono, 2010). Informasi penurunan tanah tersebut memberikan berbagai manfaat untuk diantaranya: perencanaan tata ruang, perencanaan pembangunan sarana dan prasarana, pelestarian lingkungan, pengendalian intruksi air laut, pengendalian pengambilan air tanah, pengendalian banjir. (Abidin, 2007) Pada penelitian ini akan dilakukan analisa mengenai penurunan tanah yang diakibatkan oleh penggunanaan air tanah yang keluar bersama lumpur dan gas pada daerah sekitar lumpur lapindo Sidoarjo, dengan menggunakan metode GPS untuk mendapatkan nilai penurunan tanah di Sidoarjo. 1.2 Rumusan dan Batasan Masalah Perumusan masalah dalam penelitian ini adalah bagaimana menganalisa pergerakan tanah yang diakibatkan oleh keluarnya lumpur dan gas di kecamatan Porong, Kabupaten Sidoarjo. Batasan masalah dari penelitian ini adalah: a. Penelitian dilakukan di daerah Kecamatan Porong, Kabupaten Sidoarjo. b. Data Primer yang digunakan adalah data pengamatan GPS pada titik-titik kontrol/pengamatan penurunan tanah di Sidoarjo.
2
c. Metode yang digunakan dalam penelitian ini adalah pengukuran GPS Statik dengan pengukuran radial untuk mendapatkan nilai penurunan tanah di titiktitik kontrol. d. Hasil penelitian adalah nilai penurunan tanah akibat dari keluarnya air dan lumpur ke permukaan tanah secara terus menerus. e. Peta kerentanan gerakan tanah akibat keluarnya lumpur dan gas di daerah kecamatan Porong, kabupaten Sidoarjo.
1.3 Tujuan Penelitian Penelitian penurunan tanah akibat dari
keluarnya air dan lumpur ke
permukaan tanah secara terus menerus yang terjadi di Sidoarjo ini bertujuan : 1. Memberikan informasi mengenai penurunan tanah yang terjadi di titik pengamatan dari data pengamatan pada titik-titik pengamatan dengan menggunakan metode GPS. 2. Mengetahui penurunan tanah dan penaikkan tanah pada titik-titik pengamatan akibat dari keluarnya lumpur dan gas ke permukaan tanah secara terus menerus.
1.4 Manfaat Penelitian Manfaat dari penelitian penurunan tanah ini adalah: 1.
Dapat digunakan untuk rujukan bahwa data pengamatan GPS dapat digunakan untuk pemantauan pergerakan tanah secara vertical dengan hasil yang akurat.
2.
Informasi ada tidaknya penurunan tanah di daerah sekitar Lumpur Lapindo dalam kurun waktu yang singkat dapat digunakan untuk berbagai keperluan penelitian, mitigasi bencana dan keperluan rencana pembangunan.
3.
Memberikan informasi perubahan penurunan tanah yang terbaru dengan data pengamatan GPS.
3
4
“Halaman ini sengaja dikosongkan”
5
BAB II KAJIAN PUSTAKA DAN DASAR TEORI 2.1 Penurunan Tanah (Land Subsidence) Penurunan tanah (landsubsidence) merupakan salah satu fenomena deformasi permukaan bumi secara vertikal disamping terjadi fenomena uplift. Penurunan tanah merupakan suatu fenomena alam yang banyak terjadi di kotakota besar. Penurunan tanah alami terjadi secara regional yaitu meliputi daerah yang luas atau terjadi secara lokal yaitu hanya sebagian kecil permukaan tanah. Hal ini biasanya disebabkan oleh adanya rongga di bawah permukaan tanah, biasanya terjadi didaerah yang berkapur (Whittaker and Reddish, 1989). Secara garis besar penurunan tanah bisa disebabkan oleh beberapa hal (Whittaker and Reddish, 1989), salah satunya yaitu penurunan tanah alami. Penurunan tanah alami (natural subsidence) yang disebabkan oleh proses – proses geologi seperti siklus geologi, sedimentasi daerah cekungan dan sebagainya. Beberapa penyebab terjadinya penurunan tanah alami bisa digolongkan menjadi: a. Siklus Geologi Penurunan muka tanah terkait dengan siklus geologi. Proses – proses yang terlihat dalam siklus geologi adalah: pelapukan (denuation), pengendapan (deposition), dan pergerakan kerak bumi (crustal movement). Adapun keterkaitannya yaitu pelapukan bisa disebabkan oleh air seperti pelapukan batuan karena erosi baik secara mekanis maupun kimia, oleh perubahan temperature yang mengakibatkan terurainya permukaan batuan, oleh angin terutama di daerah yang kering dan gersang karena pengaruh glacial dan oleh gelombang yang biasanya terjadi di daerah pantai (abrasi). b. Sedimentasi Daerah Cekungan Biasanya daerah cekungan terdapat di daerah – daerah tektonik lempeng terutama di dekat perbatasan lempeng. Sedimen yang terkumpul di cekungan semakin lama semakin banyak dan menimbulkan beban yang
5
bekerja semakin meningkat, kemudian proses kompaksi sedimen tersebut menyebabkan terjadinya penurunan pada permukaan tanah. Sebagian besar penurunan muka tanah akibat faktor ini adalah : 1. Adanya gaya berat dari beban yang ditimbulkan oleh endapan dan juga ditambah dengan air menyebabkan kelenturan pada lapisan kerak bumi. 2. Aktivitas internal yang menyebabkan naiknya temperature kerak bumi dan kemudian mengembang menyebabkan kenaikan pada permukaan pada permukaan tanah. Setelah itu proses erosi dan pendinginan kembali menyebabkan penurunan muka tanah. 3. Karakteristik deformasi dari lapisan tanah yang berkaitan dengan tekanan – tekanan yang ada. Beberapa faktor penyebab terjadinya penurunan tanah (Bakti, 2010), antara lain: 1. Pengambilan air tanah yang berlebihan 2. Penurunan karena beban bangunan 3. Konsolidasi alamiah lapisan tanah 4. Gaya-gaya tektonik 5. Ekstraksi gas dan minyak bumi 6. Penambangan bawah tanah 7. Ekstraksi lumpur 8. Patahan kerak bumi 9. Kontraksi panas bumi di lapisan litosfer Menurut Bakti (2010), penurunan tanah di suatu wilayah bisa dipelajari dengan beberapa metode, misalnya metode hidrogeologis, pengamatan level muka air tanah dan pengamatan dengan ekstensometer dan piezometer, maupun metode geodetik seperti survei sipat datar (leveling), survei Global Positioning System (GPS) dan Interferometric Syntetic Aperture Radar (InSAR). Karena laju penurunan tanah umumnya relatif lambat, pemantauan subsidensi lebih efisien dilakukan episodik atau periodik. Karena data dan informasi tentang penurunan muka tanah akan sangat bermanfaat bagi aspek-aspek pembangunan seperti untuk perencanaan tata ruang 6
(di atas maupun di bawah permukaan tanah), perencanaan pembangunan sarana /prasarana, pelestarian lingkungan, pengendalian dan pengambilan air tanah, pengendalian instrusi air laut serta perlndungan masyarakat (linmas) dari dampak penurunan tanah (seperti terjadinya banjir), maka sudah sewajarnya bahwa informasi tentang karakteristik penurunan tanah ini perlu diketahui dengan sebaikbaiknya dengan pemantauan secara berkesinambungan (Bakti, 2010). Dari penelitian yang dilakukan oleh Abidin, dkk (2010), disebutkan bahwa penurunan tanah di sekitar lumpur Sidoarjo diakibatkan oleh: 1. Beban lumpur yang ada di atasnya, 2. Naiknya lumpur yang ada dibawah permukaan, 3. Lahan pemukiman yang ada diatasnya karena adanya pekerjaan permukaan, seperti pembuatan tanggul. Data dan
informasi tentang penurunan
muka tanah akan sangat
bermanfaat bagi aspek- aspek pembangunan seperti untuk perencanaan tata ruang (di atas maupun di bawah permukaan tanah), perencanaan pembangunan sarana/prasarana, pelestarian lingkungan, pengendalian dan pengambilan airtanah, pengendalian intrusi air laut, serta perlindungan masyarakat (linmas) dari dampak penurunan tanah (seperti terjadinya banjir), sehingga informasi tentang karakteristik penurunan
tanah ini perlu diketahui dengan sebaik-baiknya,
dipelajari dan dipantau secara berkesinambungan. Teknik pemantauan land subsidence dari suatu wilayah dapat dipantau dengan menggunakan beberapa metode, baik itu metode-metode hidrogeologis (e.g. pengamatan level muka air tanah serta pengamatan dengan ekstensometer dan piezometer yang diinversikan kedalam besaran penurunan muka tanah) dan metode geoteknik, maupun metode-metode geodetik seperti survei sipat datar (leveling), survei gaya berat mikro, survei GPS dan InSAR. (Abidin dkk, 2008)
2.2 Teknik Pemantauan Land Subsidence Dengan GPS GPS adalah sistem satelit navigasi dan penentuan posisi yang berbasiskan pada pengamatan satelit-satelit Global Positioning System (Abidin, 2000; Hofmann-Wellenhof et al., 1997). Prinsip studi penurunah tanah dengan metode
7
survei GPS yaitu dengan menempatkan beberapa titik pantau di beberapa lokasi yang dipilih, secara periodik untuk ditentukan koordinatnya secara teliti dengan menggunakan metode survei GPS. Dengan mempelajari pola dan kecepatan perubahan koordinat dari titik-titik tersebut dari survei yang satu ke survei berikutnya, maka karakteristik penurunan tanah akan dapat dihitung dan dipelajari lebih lanjut. Studi permukaan dengan metode survey GPS mempunyai beberapa keunggulan dan keuntungan, diantara lain yaitu: a. GPS memberikan nilai vektor pergerakan tanah dalam tiga dimensi (dua komponen horisontal dan satu komponen vertikal). Jadi disamping memberikan informasi tentang besarnya penurunan muka tanah, GPS juga sekaligus memberikan informasi tentang pergerakan tanah dalam arah horisontal. b. GPS memberikan nilai vektor pergerakan dan penurunan tanah dalam suatu sistem koordinat referensi yang tunggal. Dengan itu maka GPS dapat digunakan untuk memantau pergerakan suatu wilayah secara regional secara efektif dan efisien. c. GPS dapat memberikan nilai vektor pergerakan dengan tingkat presisi sampai beberapa mm, dengan konsistensi yang tinggi baik secara spasial maupun temporal. Dengan tingkat presisi yang tinggi dan konsisten ini maka diharapkan besarnya pergerakan dan penurunan tanah yang kecil sekalipun akan dapat terdeteksi dengan baik. d. GPS dapat dimanfaatkan secara kontinyu tanpa tergantung waktu (siang maupun malam), dalam segala kondisi cuaca. Dengan karakteristik semacam ini maka pelaksanaan survei GPS untuk pemantauan pergerakan dan penurunan muka tanah dapat dilaksanakan secara efektif dan fleksibel.
2.3 Global Positioning System (GPS)
8
GPS merupakan suatu metode penetuan posisi dengan menggunakan satelit GPS yang dikelola oleh Amerika Serikat (Abidin, 1995). Sistem ini sudah mulai banyak digunakan banyak kalangan baik untuk keperluan navigasi maupun untuk penentuan posisi dan bisa digunakan dalam kondisi segala cuaca maupun waktu. GPS terdiri dari 3 segmen, yaitu konstelasi satelit, master kontrol, dan segmen pengguna. Sistem konstelasi satelit GPS terdiri dari 24 satelit yang terletak pada 6 bidang orbit yang berpusat ke bumi, dengan jumlah satelit 4 buah pada setiap bidang orbitnya. Satelit GPS memiliki periode orbit 11 jam 58 menit, dengan inklinasi 55 derajat terhadap bidang ekuator. 2.3.1
Pengamatan GPS Pengamatan yang digunakan dalam penentuan posisi antara lain
menggunakan pseudorange dan carrier phase. Pseudorange merupakan jarak antara satelit dengan receiver pada epok tertentu yang didapatkan dari penerimaan sinyal dari satelit GPS. Persamaan pada data pengamatan GPS melalui jarak semu (pseudorange) pada persamaan 1 dan jarak fase (phase range) pada persamaan 2. (Abidin, 2007) : ...(1) .
.....(2) Dimana: Pi
: pseudorange pada frekuensi fi (m), (i=1,2)
Li
: jarak fase pada frekuensi fi (m), (i=1,2)
ρ
: jarak geometris antara pengamat dengan satelit
c
: cepat rambat gelombang : panjang gelombang
dρ
: kesalahan jarak yang diakibatkan kesalahan orbit
dtrop
: bias yang disebabkan oleh refraksi troposfer (m)
dion
: bias yang disebabkan oleh refraksi ionosfer (m)
dt,dT
: kesalahan pada jam receiver dan jam satelit (m) 9
M Pi, Mci : efek multipath pada hasil pengamatan Pi dan Li (m) N1,N2
: ambiguitas fase dari sinyal L1 dan L2 (dalam jumlah gelombang)
,
: gangguan (noise) dalam hasil pengamatan Pi dan Li (m)
2.3.2 Penentuan Posisi dengan GPS Pada dasarnya penentuan posisi dengan GPS adalah pengukuran jarak secara bersama-sama ke beberapa satelit (yang koordinatnya telah dikeahui) sekaligus. Untuk menentukan suatu koordinat suatu titik di bumi, receiver setidaknya membutuhkan 4 satelit yang dapat ditangkap sinyalnya dengan baik. Secara default posisi atau koordinat yang diperoleh bereferensi ke global datum yaitu World Geodetic System 1984 atau disingkat WGS ’84. Secara garis besar penentuan posisi dengan GPS ini dibagi menjadi dua metode yaitu metode absolut dan relatif. 1. Metode absolut atau juga dikenal dengan point positioning, merupakan metode untuk menentukan posisi hanya berdasarkan pada satu pesawat penerima (receiver) saja dan tipe receiver yang digunakan untuk keperluan ini adalah tipe navigasi. Ketelitian posisi yang diperoleh sangat tergantung pada tingkat ketelitian data serta geometri satelit. Metode ini tidak digunakan untuk penentuan posisi yang teliti. Aplikasi utama metode ini adalah untuk keperluan navigasi atau aplikasi-aplikasi lain yang memerlukan informasi posisi yang tidak perlu terlalu teliti tetapi tersedia secara instan (real time), seperti untuk keperluan reconnaissance (orientasi lapangan) dan ground trothing (pengecekan lapangan). 2. Metode relatif atau sering disebut differential positioning, merupakan metode untuk menentukan posisi dengan menggunakan lebih dari sebuah receiver. Satu GPS dipasang pada lokasi tertentu di muka bumi dan secara terus menerus menerima sinyal satelit dalam jangka waktu tertentu dijadikan sebagai referensi bagi yang lainnya. Metode ini menghasilkan posisi berketelitian tinggi dan diaplikasikan untuk keperluan survei Geodesi ataupun pemetaan yang memerlukan ketelitian tinggi.
10
2.3.3 Penentuan Tinggi dengan GPS Ketinggian titik yang diberikan GPS adalah ketinggian titik di atas permukaan ellipsoid, yaitu ellipsoid WGS ’84. Tinggi ellipsoid (h) tidak sama dengan dengan tinggi orthometrik yang merupakan tinggi di atas geoid yang diukukr sepanjang garis gaya berat melalui titik tersebut.
H h
H = tinggi ortometrik (bereferensi ke geoid) N = tinggi (undulasi) geoid di atas ellipsoid h = tinggi ellipsoid (bereferensi ke ellipsoid)
N
H=h-N
Gambar 2.1 Sistem tinggi (Hofmann-Wellenhof, et al., 2006) Hasil penelitian menunjukan bahwa ketelitian komponen tinggi GPS umumnya 2-3 kali lebih rendah dibandingkan ketelitian komponen horisontalnya. Misalkan saat ini ketelitian horisontal sudah dapat mencapai 2 milimeter, maka komponen tingginya menjadi 3 kali lipatnya, yaitu mencapai 6 milimeter. Ada dua penyebab utama dalam hal ini, yaitu (Abidin, 2007): a. Satelit-satelit GPS yang bisa diamati hanya yang berada diatas horison (one-sided geometry). Pada dasarnya hal ini secara geometrik tidak optimal, juga karena tidak ada satelit diatas dan dibawah pengamat, maka tidak akan ada efek pengeliminasian kesalahan seperti halnya dalam kasus komponen horisontal. Pada komponen horisontal, adanya satelit disebelah barat dan timur 11
ataupun di utara dan selatan dari pengamat akan memungkinkan adanya pengeliminasian tersebut. b. Efek kesalahan dan bias (ionosfer, troposfer, dan orbit) umumnya adalah pada jarak, yaitu menyebabkan jarak menjadi semakin panjang atau dapat juga menjadi semakin pendek. Dalam hal ini maka yang paling terpengaruh adalah komponen tinggi
2.3.4 Kesalahan dan Bias Pada saat pengambilan data menggunakan GPS, sinyal satelit akan dipengaruhi oleh beberapa kesalahan dan bias. Kesalahan dan bias pada dasarnya dapat dikelompokkan menjadi: 1. Satelit , seperti kesalahan ephemeris, jam satelit, dan selective availability (SA) 2. Medium propagansi, seperti bias ionosfer dan bias troposfer 3. Receiver GPS, seperti kesalahan jam receiver, kesalahan yang terkait dengan antena, dan noise 4. Data pengamatan, ambiguitas fase dan cycle slips 5. Lingkungan sekitar GPS receiver seperti multipath dan imaging Kesalahan dan bias GPS harus diperhitungkan secara benar dan baik, karena besar dan karakteristik dari kesalahan dan bias tersebut akan mempengaruhi ketelitian informasi (posisi, kecepatan, percepatan, waktu) yang diperoleh serta proses penentuan ambiguitas fase dari sinyal GPS. Berikut ini karakteristik dari sebagian kesalahan dan bias akan dijelaskan secara umum. (Abidin, 2007). Kesalahan Ephemeris (orbit) adalah kesalahan dimana orbit satelit yang dilaporkan oleh ephemeris satelit tidak sama dengan orbit satelit yang sebenarnya. Dengan kata lain, posisi satelit yang dilaporkan tidak sama dengan posisi yang sebenarnya. Penyebab kesalahan ephemeris diantaranya: Proses perhitungan orbit satelit oleh stasiun-stasiun pengontrol yang tidak teliti
12
Kesalahan pada saat memprediksi orbit untuk periode waktu setelah record (perekaman) ke satelit
Ada beberapa cara untuk mengatasi efek kesalahan ini antara lain: Menerapkan metode differensial positioning Memperpendek panjang baseline Memperpanjang waktu pengamatan Tentukan parameter kesalahan orbit dalam proses estimasi Menggunakan informasi orbit yang lebih teliti seperti precise ephemeris
Posisi satelit yang dilaporkan
Posisi satelit sebenarnya
Gambar 2.2 Kesalahan ephemeris (Abidin, 2007) Pada dasarnya, konsep penentuan posisi dengan GPS adalah pengikatan ke belakang, yaitu dengan pengukuran jarak secara simultan ke beberapa satelit GPS yang diketahui koordinatnya. Dalam hal ini receiver GPS harus mempunyai 2 data utama yaitu: 1. Jarak antara receiver GPS dengan setiap satelit GPS yang diamati. 2. Posisi/koordinat dari setiap satelit GPS yang diamati. Posisi yang diberikan GPS adalah posisi tiga dimensi yang dinyatakan dalam Datum WGS 1984, selain itu GPS juga memberikan informasi waktu (t). Secara vektor, prinsip dasar penentuan posisi dengan GPS diperlihatkan oleh gambar 2.5. Dalam hal ini, parameter yang akan ditentukan adalah vektor posisi geosentrik pengamat (R). Untuk itu, karena vektor posisi satelit GPS (r) telah diketahui, maka yang perlu ditentukan adalah vektor posisi toposentris satelit terhadap pengamat (ρ).
13
Gambar 2.3 Prinsip Dasar Penentuan Posisi Satelit GPS (Abidin, 2007) Jarak dari receiver GPS ke satelit GPS dapat dihitung dari sinyal yang diterima dari satelit GPS. Sedangkan koordinat satelit GPS diperoleh dari informasi yang dikirimkan oleh satelit GPS tersebut bersama dengan sinyal yang diterima oleh receiver GPS. Untuk dapat menghitung koordinat receiver GPS, paling sedikit harus ada 4 satelit yang teramati. Posisi yang diberikan oleh GPS adalah posisi 3 dimensi. Pada penentuan posisi menggunakan GPS diperlukan adanya minimal 6 GPS / GLONASS dengan mask angle 15 derajat di ketinggian terhadap prosedur yang direkomendasikan oleh TPS dalam manual alat. Di daerah-daerah multipath tinggi, selama periode PDOP tinggi dan selama periode kinerja aktivitas ionosfer tinggi dapat terdegradasi. Prosedur pemeriksaan Kuat sangat dianjurkan di daerah multipath ekstrim atau di bawah dedaunan lebat. DOP (Dilution of Precision) adalah istilah yang menggambarkan kekuatan akurasi dari konfigurasi geometri satelit. Ketika satelit yang muncul berdekatan di atas langit, geometrinya menjadi lemah dan DOP-nya tinggi. Namun ketika berjauhan maka geometrinya menjadi kuat dan nilai DOP menjadi rendah. Sehingga dengan rendahnya nilai DOP maka kekakuratan posisi GPS menjadi lebih baik karena adanya separasi sudut yang luas dari satelit-satelit yang digunakan dalam memperhitungakan posisi unit GPS. Nilai PDOP harus dibawah 5 bagi para pengguna GPS Geodetik.
14
2.3.5
Bias Ionosfer Ionosfer merupakan bagian dari lapisan atas atmosfer dimana
terdapat sejumlah elektron dan ion bebas yang mempengaruhi perambatan gelombang radio. Sinyal GPS memancarkan sinyalnya ke reciever melewati ionosfer, akibatnya ion-ion tersebut akan mempengaruhi propagasi sinyal satelit terutama pada kecepatan sinyal. Untuk menghilangkan efek ionosfer, bisa menggunakan kombinasi linier bebas. Untuk mereduksi kesalahan bias ionosfer, hal-hal yang dilakukan adalah:
Menggunakan data GPS dari dua frekuensi, L1 dan L2
Melakukan pengurangan (differencing) data pengamatan
Memperpendek panjang baseline
Menggunakan model ionosfer (model Bent atau Klobuchar)
Menggunakan parameter koreksi yang dikirimkan oleh sistem WADGPS (Wide Area Differential GPS)
2.3.6
Bias Troposfer Ketika melalui lapisan troposfer, sinyal GPS akan mengalami
refraksi, yang menyebabkan perubahan kecepatan dan arah sinyal GPS. Bias troposfer ini akan mempengaruhi kecepatan sehingga akan menghasilkan ukuran jarak yang kurang teliti. Lapisan troposfer ini memperlambat data waktu dan data fase. Untuk memodelkan bias troposfer bisa menggunakan metode Goad dan Goodman (1974) yang telah dimodifikasi model Hopfield dengan mengasumsikan bahwa suhu menurun secara linear dengan meningkatnya ketinggian di troposfer, tetapi itu tetap konstan di stratosfer (dua lapisan atmosfer), asumsi yang sama 15
digunakan oleh Saastamoinen (1973) untuk menurunkan fungsi pemetaan nya. Fungsi pemetaan Goad dan Goodman kadang-kadang disebut dalam sastra sebagai modifikasi Hopfield.
2.4 Metode Radial Metode radial merupakan metode yang digunakan dalam pengambilan data di lapangan dengan menggunakan satu base station yaitu CORS yang telah memiliki waktu pengamatan yang sangat panjang dan tempat yang cenderung lebih stabil dari pada titik receiver pengamatan di lokasi penelitian. Penerapan metode radial harus diperhitungkan dalam pergerakan receiver di lapangan, karena hal ini sangat berpengaruh pada operasional serta kekuatan jaring.
Gambar 2.4. Metode pengamatan radial (Abidin, 2007)
2.5 Analisa Regresi Analisis regresi adalah suatu proses memperkirakan secara sistematis tentang apa yang paling mungkin terjadi dimasa yang akan datang berdasarkan informasi yang sekarang dimiliki agar memperkecil kesalahan. Analisis regresi dapat juga diartikan sebagai usaha memprediksi perubahan. Perubahan nilai suatu variabel dapat disebabkan karena adanya perubahan pada variabel-variabel lain yang mempengaruhinya. Dalam fenomena alam banyak sekali kejadian yang saling berkaitan sehingga perubahan pada variabel lain berakibat pada perubahan
16
variabel lainnya. Teknik yang digunakan untuk menganalisis ini adalah analisis regresi. Menurut Kurniawan (2008) dalam bukunya yang berjudul “Regresi Linear”, analisis regresi setidak-tidaknya memiliki 3 kegunaan, yaitu untuk tujuan deskripsi dari fenomena data atau kasus yang sedang diteliti, untuk tujuan kontrol, serta untuk tujuan prediksi. Regresi mampu mendeskripsikan fenomena data melalui terbentuknya suatu model hubungan yang bersifatnya numerik. Regresi juga dapat digunakan untuk melakukan pengendalian (kontrol) terhadap model regresi yang diperoleh. Selain itu, model regresi juga dapat dimanfaatkan untuk melakukan prediksi untuk variabel terikat. Namun yang perlu diingat, prediksi di dalam konsep regresi hanya boleh dilakukan di dalam rentang data dari variabelvariabel bebas yang digunakan untuk membentuk model regresi tersebut. Dalam hal ini analisa regresi digunakan untuk membuat permodelan pada deformasi. Analisa tersebut digunakan untuk memprediksi besar dari deformasi tanah apakah yang terjadi pada tanah tersebut. Ada beberapa macam analisa regresi, diantaranya: 2.5.1 Regresi Linier Regresi linier adalah metode statistika yang digunakan untuk membentuk model hubungan antara variabel terikat (dependen; respon; Y) dengan satu atau lebih variabel bebas (independen, prediktor, X). Apabila banyaknya variabel bebas hanya ada satu, disebut sebagai regresi linier sederhana, sedangkan apabila terdapat lebih dari 1 variabel bebas, disebut regresi linier berganda. Untuk menunjukkan hubungan keterkaitan antara satu atau beberapa variabel yang nilainya sudah diketahui dengan satu variabel lain yang nilainya belum diketahui adalah dengan menggunakan persamaan regresi. Bentuk umum regresi linier sederhana adalah: y = a + bx
(2.1)
dengan : y : variabel tergantung x : variabel bebas a : titik potong y. Merupakan nilai perkiraan bagi y ketika x = 0 17
b : kemiringan garis atau perubahan rata-rata pada y untuk setiap satu unit perubahan (naik atau turun) pada variabel bebas x. Nilai b dapat positif (+) dapat pula negatif (-). Jika b bernilai positif, maka grafiknya:
y y = a + bx
b
x
a
Gambar 2.5 Grafik untuk b positif (Abidin, 2010)
Jika b bernilai negatif, maka grafiknya: y a
b y = a - bx x
Gambar 2.6 Grafik untuk b negatif (Abidin, 2010)
Penentuan nilai a dan b dengan menggunakan persamaan berikut:
∑
∑
1
∑
1
2
∑
1
∑
1 2
1
(2.2)
dengan
18
n : banyak pasangan data yi : nilai peubah tak bebas y ke i xi : nilai peubah bebas x ke i a = y - bx sehingga a=
∑
1
∑
1
(2.3) Algoritma yang digunakan dalam regresi linier adalah 1. Tentukan titik data yang diketahui dalam(xi,yi) untuk i = 1,2,3,..,N 2. Hitung nilai a dan b dengan menggunakan formulasi dari regresi linier 3. Tampilkan fungsi linier 4. Hitung fungsi linier tersebut dalam range x dan step dx tertentu 5. Tampilkan hasil tabel (xn,yn) dari hasil fungsi linier tersebut
2.5.2 Regresi Polinomial Regresi
Polinomial
merupakan
metode
statistika
yang
digunakan untuk menentukan fungsi polinomial yang paling sesuai dengan kumpulan titik data (xn , yn) yang diketahui. Bentuk umum regresi polinomial adalah: 0
1
2
2
⋯
(2,4)
Secara umum, a1 dapat dihitung dengan m+1 dengan persamaan sebagai berikut : 0
⋯
1 1
0
1
2
1 1
1
1
1
⋯ 1
1
⋮
19
0
1
1 1
2
⋯
1
1
1
(2,5) Dari persamaan diatas, dibentuklah matriks persamaan: 1
1
2
1
1 1
1
1
1
1 2
…
1
1
1
1
2
2
…
1
1
1
1
2
1 2
…
1
2
1
2
…
1
1
0
… … … … … 2
1
…
… 1
1
1
(2,6)
2.6 Uji Statistik t-Student Dalam melakukan uji hipotesis, ada banyak faktor yang menentukan, seperti apakah sampel yang diambil berjumlah banyak atau hanya sedikit; apakah standar deviasi populasi diketahui; apakah varians populasi diketahui; metode parametrik apakah yang dipakai, dan seterusnya. Langkah-langkah untuk menguji suatu hipotesis adalah: 1. Menentukan hipotesa nol (H0) dan hipotesa alternatif. H0 merupakan hipotesis nilai para meter dengan dibandingkan dengan hasil perhitungan dari sampel. H0 ditolak hanya jika hasil perhitungan dari sampel tidak mungkin memiliki kebenaran terhadap hipotesis yang ditentukan terjadi. Ha diterima jika H0 ditolak. 2. Menentukan
tingkat
signifikansi
yang
digunakan.
Tingkat
signifikansi merupakan standar statistik yang digunakan untuk menolak H0. Jika ditentukan tingkat signifikansi 5% (α = 0.05). H0
20
ditolak haya jika hasil perhitungan dari sampel sedemikian berbeda dengan nilai dugaan (yang dihipotesakan). 3.
Memilih uji statistik.
4. Menentukan statistik tabel. Nilai staistik tabel dipengaruhi oleh: a. Tingkat kepercayaan b. Derajat kebebasan c. Jumlah sampel yang didapat. 5. Menentukan statistik hitung 6. Mengambil keputusan Uji-t termasuk dalam golongan statistika parametrik. Statistik uji ini digunakan dalam pengujian hipotesis. Uji-t digunakan ketika informasi mengenai nilai variance (ragam) populasi tidak diketahui. Asumsi pengambilan uji t pada penelitian ini antara lain: 1. Datanya berdistribusi normal 2. Data kurang dari 30 sampel 3. Kedua kelompok data dependen (berpasangan) 4. variabel yang dihubungkan berbentuk numerik dan kategorik (dengan hanya 2 kelompok). Untuk menguji kebenaran nilai penurunan yang didapatkan, dilakukan uji dengan menggunakan general linear hypothesis. Uji hipotesa tersebut antara lain: Hipotesa nol :
H0 : ΔHi-j = 0
Hipotesa tandingan :
Ha : ΔHi-j ≠ 0 ∆ ∆
Uji statistik : t
= nilai pengamatan hasil hitung dalam uji t
dhij = varian nilai pengamatan H0 = titik tidak mengalami penurunan Ha = titik mengalami penurunan Uji statistik dilakukan dengan student-t test, hipotesa nol akan ditolak bila | |
, /2
21
Dimana df merupakan derajat kebebasan, dan α adalah derajat kesalahan. H0 diterima, jika thitung > ttabel H0 ditolak, jika thitung ≤ ttabel Dari uji tersebut, H0 diterima jika nilai t hitung lebih kecil dari pada nilai distibusi t pada tabel sehingga dari pernyataan ini, titik tersebut tidak mengalami penurunan maupun kenaikan. H0 ditolak jika nilai uji t lebih besar dari nilai distribusi t pada tabel sehingga dari pernyataan ini, titik tersebut mengalami penurunan ataupun kenaikan. (Caspary, 1987)
2.7 Lumpur Sidoarjo 2.7.1 Latar Belakang Lumpur Banjir
Lumpur
Panas
Sidoarjo (Lusi) adalah
peristiwa
menyemburnya lumpur panas di lokasi pengeboran PT Lapindo Brantas di Kecamatan Porong, Kabupaten Sidoarjo, Jawa Timur, pada tanggal 27 Mei 2006. Sampai saat ini, semburan tersebut belum ada tanda-tanda akan berhenti dan volumenya pun semakin bertambah. Menurut Widodo (2007) dalam bukunya yang berjudul “Memahami Bencana Gunung Lumpur Kasus Lumpur Panas Sidoarjo”, pemunculan dan peningkatan volume lumpur panas ini akan mengakibatkan konsekuensi kejadian beruntun seperti “Tabrakan Karambol” yaitu tenggelamnya sawah, pemukiman, jalan tol, pabrik, rel kereta api, jaringan pipa gas industri, jaringan
telekomunikasi, jaringan PLN, jaringan pipa air dan
infrastruktur lainnya. Konsekuensi selanjutnya akan mengakibatkan terjadinya
pengungsian,
pengangguran,
rusaknya
jaringan
telekomunikasi, jaringan PLN, dan jaringan jalan sampai terganggunya perekonomian Jawa Timur.
2.7.2 Geomorfologi Sidoarjo Geomorfologi daerah Sidoarjo dapat dibedakan menjadi dua macam, yaitu daerah pegunungan yang menempati bagian selatan dan daerah dataran yang menempati bagian utara. Daerah pegunungan dibagi 22
lagi menjadi dua yaitu daerah pegunungan vulkanik dan daerah pegunungan lipatan. Daerah pegunungan vulkanik merupakan bagian dari gunung Penanggungan
yang
tersusun
dari
endapan
dari
pegunungan
Penanggungan dan batuan endapan gunung Arjuno Purba. Daerah pegunungan lipatan berarah timur barat mencakup daerah gunung Bang dan sekitar Raci. Di daerah ini dijumpai batuan-batuan sedimen yang telah terlipat, membentuk antiklinal-antiklinal. Daerah dataran terdapat dibagian utara denga ketinggian mencapai 15-18 meter dpl merupakan delta yang dikenal sebagai delta Brantas yang yang terbentuk oleh sungai Surabaya yang mengalir dibagian utara ke Kota Surabaya dan Sungai Porong. Kedua anak ini merupakan anak-anak sungai Brantas. Di daerah dataran tersebut, terdapat semburan aktif di desa Siring pada ketinggian sekitar 3 meter dpl, dijumpai dua gunung lumpur yaitu di desa Pulungan dengan ketinggian 7 meter dpl dan di desa Kalanganyar dengan ketinggian sekitar 13 meter dpl.
Gambar 2.7 Citra satelit busur vulkanik selatan Pulau Jawa. (Mazzini,dkk.2007) 2.8 Studi Hasil Penelitian Sebelumnya Menurut Kurniawan (2011), dalam penelitiannya perkembangan kota memberikan dampak yang besar terhadap terjadinya penurunan tanah, selain 23
faktor fisik seperti tektonisme, dan konsolidasi tanah. Kota-kota seperti Jakarta, Semarang dan Surabaya mengalami pembangunan yang pesat, pusatpusat bisnis, perdagangan, hiburan dan tempat tinggal berkembang sangat pesat. Perkembangan inilah yang juga dapat mempengaruhi penurunan tanah. (Abidin,dkk. 2010). Pengamatan penurunan tanah di Kota Surabaya dibutuhkan untuk melihat kecenderungan pada pola geometric dan pola fisik yang terjadi. Teknologi Global Positioning System (GPS) dapat memberikan informasi terkait posisi 3 dimensi(x,y dan z) secara akurat, apabila sejak awal pengamatan menggunakan metode yang tepat. Dalam penelitian ini didapatkan kesimpulan bahwa penelitian ini berhasil memberikan informasi adanya penurunan tanah yang terjadi di Kota Surabaya. Penurunan tanah yang terjadi dibeberapa wilayah Kota Surabaya memiliki kecepatan penurunan tanah yang berbeda. Kecepatan penurunan tanah diketahui terendah sebesar -0,322 cm/tahun, kecepatan tertinggi sebesar -2,792 cm/tahun, dan rata-rata sebesar -1,21 cm/tahun. Dalam Pengamatan GPS sebanyak 2 kala (tahun 2007-2010) menunjukkan adanya potensi penurunan tanah pada beberapa lokasi dikawasan Kota Surabaya. Penurunan terbesar terjadi pada titik pengamatan RUKT yang terletak di kawasan industri rungkut. Penurunan kemungkinan terjadi akibat penggunaan air bawah tanah oleh disekitar titik pengamatan secara besar. Penurunan yang terjadi pada titik pengamatan KOTA dan WARU kemungkinan terjadi akibat besarnya beban oleh bangunan dan/ atau penggunaan air tanah untuk kebutuhan pusat-pusat bisnis, perkantoran, pusat perbelanjaan, hotel dan hiburan yang tersebar disekitar titik pengamatan. Penurunan tanah pada titik pengamatan KLNK kemungkinan akibat konsolidasi tanah alluvial yang belum stabil dan/atau besar jalan dari kendaraan-kendaraan berat setiap harinya. Penurunan tanah di titik KJRN terjadi meskipun relatif kecil dibanding titik yang lainnya, hal ini kemungkinan pada wilayah tersebut belum banyak berdiri bangunanbangunan tinggi penggunaan air bawah tanah yang tidak terlampau banyak. Pada titik SBRT perubahan elevasi cenderung tidak ada / stabil, hal ini dikarenakan secara Geologi lokasi titik tersebut terletak pada batuan keras di 24
Formasi Lidah (Sukardi, 1992) relatif lebih stabil dibandingkan titik yang lainnya.
Menurut Bukhori (2011) penelitian penurunan tanah di sekitar lumpur dengan menggunakan GPS yang tersebar di 18 titik yang berbeda disekitar lumpur. Tujuan dari penelitian ini untuk mendapatkan besar penurunan dan kenaikan tanah yang terjadi di sekitar Lumpur Lapindo dengan metode polinomial. Dari pengukuran GPS tersebut, didapatkan bentuk grafik polinomial yang dapat diketahui apakah terjadi kenaikan atau penurunan (uplift atau subsidence) di sekitar Lumpur Lapindo. Prawoko, Eko Yuli, Taufik (2008) melakukan penelitian tentang penurunan tanah disekitar Lumpur Lapindo dengan menggunakan GPS dengan metode analisa arah gerakan tanah dan memprediksi perubahan tinggi titik pada tiap kala. Hasil dari pengukuran GPS menunjukkan pergeseran vertikal yang terjadi rata-rata antara -0,1275 m sampai dengan 0,2574 m. Hutomo, Eko Yuli, (2010) melakukan studi tentang pemanfaatan CORS GPS milik ITS dalam hal ini di gedung Teknik Geomatika. Penelitian ini memanfaatkan CORS GPS sebagai alternatif stasiun referensi dalam pengukuran GPS. Nilai penurunan tanah terbesar adalah sebesar -0.018 m di titik BT03 dan nilai kenaikan tanah terbesar adalah sebesar 0.0115 m di titik 1304. Nilai korelasi dan determinasi yang terbesar adalah pada regresi polinomial ordo 3 yaitu dengan nilai korelasi rata-rata sebesar 0.893 dan nilai determinasi sebesar 0.810.
25
26
“Halaman ini sengaja dikosongkan”
27
BAB III METODA PENELITIAN 3.1 Lokasi Penelitian Lokasi penelitian ini mengambil daerah studi di Kabupaten Sidoarjo terletak antara 112 5’ dan 112 9’ Bujur Timur dan antara 7 3’ dan 7 5’ Lintang Selatan. Batas sebelah utara adalah Kotamadya Surabaya dan Kabupaten Gresik, sebelah selatan adalah Kabupaten Pasuruan, sebelah timur adalah Selat Madura dan sebelah barat adalah Kabupaten Mojokerto
Gambar 3.1 Daerah pengamatan dalam penelitian.
27
Tabel 3.1. Koordinat sebaran titik-titik kontrol pengamatan Nama Titik
Timur (m)
Utara (m)
BPNRI-53
9168497
691251
BT01
9170380
687275
BT03
9168700
687447
BW13
9166787
687841
DG09
9165308
689811
DG11
9168770
690071
DG18
9165896
684993
DG24
9169010
692360
DG34
9166171
683191
DGIJK
9168341
684425
PT06
9170158
692979
PLNG
9168181
688091
PT11
9171031
685757
TTG-1304
9162606
686566
TTG-1305
9165632
687353
TTG-1307
9171609
688853
VK09
9167321
690298
VK13
9165694
688738
VK14
9165749
692869
VK23
9169335
688400
3.2 Data Data yang digunakan pada penelitian ini terdiri dari data spasial dan data tabular yang merupakan data sekunder dan data lapangan yang merupakan data primer. Data primer yang digunakan dalam penelitian ini adalah data pengamatan GPS pada titik-titik pengamatan penurunan tanah di Kota Surabaya tahun 2010 dan 2011 Data yang digunakan dalam penelitian ini adalah data pengamatan 28
dengan menggunakan GPS pada 20 lokasi Benchmark yang tersebar diwilayah Kecamatan Porong, Kecamatan Jabon, dan Kecamatan Tanggulangin di Kabupaten Sidoarjo serta Kecamatan Gempol di Kabupaten Pasuruan. Pengamatan ini dilakukan selama 3 kala yang diambil setiap 1 bulan sekali mulai dari bulan Maret – Mei 2014 dengan lama pengamatan setiap titiknya adalah 6 jam. Data pengamatan Stasiun GPS yang digunakan adalah yang terletak di Teknik Geomatika ITS yaitu CORS ITS.
3.3 Metodologi Penelitian Tahapan dalam pengolahan data ini adalah:
Gambar 3.2 Diagram Tahap Pengolahan Data
29
Berikut adalah penjelasan diagram alir pengolahan data: 3.3.1 Tahap Persiapan Pada tahap ini, kegiatan-kegiatan yang dilakukan adalah : Studi Literatur Pada tahap ini pengumpulan literatur mengenai penurunan tanah, info-info terkini mengenai Bencana Lumpur Sidoarjo serta segala hal mengenai GPS yang akan digunakan dalam penelitian kali ini, CORS sebagai teknologi alternatif stasiun referensi. Literatur yang digunakan dalam bentuk apapun, semisal buku, jurnal ilmiah, majalah, surat kabar, internet dan lain-lain. Orientasi Lapangan Orientasi lapangan wajib dilakukan untuk mengetahui kondisi di lapangan sebenarnya. Dari orientasi lapangan ini dapat membuat keputusan bagaimana rencana pengukuran yang akan dilakukan. 3.3.2 Tahap Pengolahan Data Tahapan dalam Pengolahan Data dan analisa adalah sebagai berikut:
Pengambilan data langsung dari lapangan dengan GPS Topcon HyperPro selama kurang lebih 6 jam tergantung panjang baseline dan data di download dengan software PCCDU (software bawaan TOPCON TOOL)
Data GPS harian diambil dari CORS ITS dan data ephemeris satelit di International GNSS Service. Data precise ephemeris digunakan untuk mereduksi semua kesalahan karena kesalahan satelit GPS sendiri. Sehingga diharapkan mendapat hasil yang lebih baik dalam pengolahan data nantinya.
Pengolahan data dengan software Topcon Tools yang diunduh langsung dari receiver GPS Topcon HyperPro. Data yang dimasukkan antara lain adalah data hasil pengukuran di lapangan, data baseline dari CORS ITS dengan titik-titik BM yang ada dalam penelitian. Dimana data tersebut diolah untuk didapatkan data
30
elevasi titik tersebut. Data-data tersebut diolah sehingga didapatkan koordinat masing-masing titik dengan ketelitian pada software TOPCON TOOLS adalah 3mm untuk x & y dan 5mm untuk z Pengolahan data di TOPCON TOOLS menggunakan Metode perataan kuadrat terkecil (Least Square)
Pengolahan Data Tinggi titik per kala pengamatan Data yang diambil setelah pengolahan data adalah data tinggi per titik yang hasilnya dibandingkan dengan tinggi sebelumnya
3.3.3 Tahap Analisa Data Analisa elevasi tiap titik menggunakan analisa t-student Dari hasil pengolahan, dilakukan uji t-student untuk menguji nilai kebenaran penurunan tanah. Dari hasil uji masing-masing beda tinggi dibandingkan dengan tabel distribusi t yang sesuai. Pada pengujian statistik ini digunakan pengujian dengan tingkat kepercayaan 95% (α =5%), dan 99% (α = 1% ). Pembuatan Model Regresi Polinomial Elevasi setiap titik digambarkan dengan model regresi polinomial ordo dua atau quadratic. Model ini digunakan untuk memprediksi apakah titik pengukuran mengalami penurunan atau kenaikan tanah. Setelah dibuat modelnya, dilakukan analisa menggunakan koefisien korelasi dan determinasi. Analisa ini digunakan untuk menilai apakah model regresi yang digunakan untuk memodelkan uplift dan subsidence.
31
”Halaman ini sengaja dikosongkan”
32
BAB IV HASIL PENELITIAN DAN PEMBAHASAN
4.1 Hasil Data penelitian didapatkan langsung dari pengukuran di lapangan mulai bulan Maret – Mei 2014 dengan menggunakan alat reciever GPS yang meliputi 20 lokasi/titik.
Gambar 4.1. Sebaran titik pengamatan
Data yang diambil dari lapangan kemudian diolah dengan menggunakan software TOPCON TOOLS. Hasil pengolahan dapat dilihat dalam tabel berikut:
33
Tabel 4.1 hasil pengolahan data kala 1 bulan Maret 2014 NO
Koordinat
POSISI Y (m)
X(m)
Z(m)
tinggi antena(m)
Standart Deviasi N(m)
E(m)
U(m)
1 BPNRI-53
9168497.123
691251.778
30.602
1.334
0.01
0.038
0.019
2 BT01
9170380.973
687275.599
32.285
1.161
0.007
0.029
0.017
3 BT03
9168700.731
687447.326
32.727
1.323
0.046
0.115
0.106
4 BW13
9166787.789
687841.592
33.482
1.24
0.029
0.071
0.086
5 DG09
9165308.179
689811.42
34.803
1.307
0.005
0.014
0.016
6 DG11
9168770.289
690071.219
31.347
1.22
0.014
0.038
0.032
7 DG18
9165896.403
684993.558
35.934
1.12
0.009
0.011
0.026
8 DG24
9169010.760
692360.482
31.123
1.128
0.063
0.089
0.122
9 DG34
9166171.182
683191.065
36.447
1.041
0.01
0.03
0.019
10 DGIJK
9168341.289
684425.374
33.903
1.032
0.06
0.136
0.126
11 PT06
9170158.495
692979.909
30.497
1.217
0.026
0.051
0.05
12 PLNG
9168181.325
688091.594
31.408
1.214
0.014
0.044
0.042
13 PT11
9171031.650
685757.177
34.142
0.675
0.011
0.037
0.025
14 TTG-1304
9162606.096
686566.847
43.258
1.168
0.006
0.025
0.012
15 TTG-1305
9165632.400
687353.520
38.992
1.195
0.014
0.044
0.028
16 TTG-1307
9171609.592
688853.560
32.399
1.185
0.012
0.038
0.039
17 VK09
9167321.956
690298.337
31.733
1.095
0.027
0.052
0.052
18 VK13
9165694.08
688738.846
37.762
1.138
0.007
0.009
0.018
19 VK14
9165749.683
692869.647
36.890
1.175
0
0
0
20 VK23
9169335.059
688400.679
32.771
1.14
0.014
0.041
0.028
34
Tabel 4.2 hasil pengolahan data kala 2 bulan April 2014 NO
Koordinat
POSISI Y (m)
X(m)
Z(m)
tinggi antena(m)
Standart Deviasi N(m)
E(m)
U(m)
1 TTG-1307
9171609.551
688853.604
32.273
1.221
0.005
0.016
0.010
2 PT11
9171031.623
685757.094
34.202
0.675
0.019
0.051
0.038
3 BT01
9170380.942
687275.739
32.373
1.161
0.028
0.065
0.055
4 VK23
9169334.852
688400.489
33.099
1.045
0.021
0.051
0.045
5 PLNG
9168181.305
688091.626
31.158
1.22
0.094
0.168
0.198
6 BPNRI-53
9168497.102
691251.551
30.685
1.247
0.026
0.079
0.043
7 DGIJK
9168341.398
684425.251
33.856
1.235
0.048
0.085
0.097
8 DG18
9165896.395
684993.617
36.058
1.21
0.013
0.024
0.032
9 PT06
9170158.548
692979.912
30.577
1.324
0.029
0.092
0.067
10 BT03
9168700.255
687447.343
32.696
1.439
0.063
0.120
0.106
11 VK13
9165694.066
688738.864
37.830
1.082
0.006
0.021
0.013
12 VK09
9167322.043
690298.449
31.625
1.051
0.014
0.033
0.022
13 VK14
9165749.049
692869.567
36.411
1.1
0.015
0.027
0.038
14 TTG-1305
9165632.518
687353.247
39.024
1.13
0.011
0.037
0.025
15 DG24
9169011.012
692360.517
30.974
1.273
0.010
0.022
0.016
16 DG11
9168770.548
690071.205
31.553
1.15
0.029
0.062
0.068
17 TTG-1304
9162606.097
686566.782
43.291
1.093
0.007
0.020
0.015
18 DG09
9165308.248
689811.304
34.972
1.152
0.012
0.032
0.031
19 BW13
9166787.838
687841.623
33.534
1.14
0.006
0.020
0.015
20 DG34
9166171.198
683191.13
36.425
1.305
0.016
0.040
0.030
35
Tabel 4.3 hasil pengolahan data kala 3 bulan Mei 2014 NO
Koordinat
POSISI Y (m)
1 BPNRI-53
X(m)
Z(m)
tinggi antena(m)
Standart Deviasi N(m)
E(m)
U(m)
9168497.13
691251.716
30.581
1.342
0.037
0.067
0.089
2 BT01
9170381.108
687275.761
32.355
1.181
0.019
0.010
0.019
3 BT03
9168700.867
687447.333
32.705
1.416
0.016
0.021
0.05
4 BW13
9166787.841
687841.622
33.535
1.139
0.007
0.018
0.009
5 DG09
9165308.324
689811.332
34.823
1.126
0.01
0.012
0.03
6 DG11
9168770.556
690071.213
31.364
1.230
0.019
0.010
0.019
7 DG18
9165896.414
684993.592
36.010
1.201
0.022
0.011
0.022
8 DG24
9169011.047
692360.501
31.069
1.243
0.007
0.016
0.018
9 DG34
9166171.210
683191.125
36.441
1.236
0.025
0.013
0.026
10 DGIJK
9168341.429
684425.414
33.891
1.156
0.053
0.044
0.055
11 PLNG
9168181.304
688091.633
30.521
1.297
0.006
0.008
0.019
12 PT06
9170158.407
692979.947
31.351
1.347
0.005
0.005
0.019
13 PT11
9171031.671
685757.001
34.189
0.703
0.005
0.02
0.011
14 TTG-1304
9162606.102
686566.79
43.262
1.096
0.014
0.037
0.032
15 TTG-1305
9165632.517
687353.509
39.049
1.173
0.013
0.018
0.038
16 TTG-1307
9171609.584
688853.603
32.469
1.190
0.004
0.015
0.008
17 VK09
9167321.993
690298.528
31.784
1.143
0.006
0.007
0.021
18 VK13
9165694.084
688738.843
37.802
1.206
0.008
0.018
0.010
19 VK14
9165749.571
692869.543
36.854
1.181
0.006
0.018
0.008
20 VK23
9169335.019
688400.655
32.828
1.230
0.007
0.017
0.009
36
Berdasarkan data di atas, maka didapatkan perbedaan tinggi dari masing-masing titik dari tiap kala, yaitu: Tabel 4.4 Nilai Subsidence Maupun Uplift Setiap Titik Nama Titik
Beda Tinggi ∆H
Perambatan kesalahan
(m)
(m)
kala 1-2
kala 1-3
kala 1-2
kala 1-3
BPNRI-53
0.083
-0.021
0.021
0.091
BT01
0.088
0.070
0.042
0.025
BT03
-0.031
-0.022
0.119
0.117
BW13
0.052
0.053
0.097
0.086
DG09
0.007
0.020
0.199
0.034
DG11
0.088
0.017
0.054
0.037
DG18
0.063
0.076
0.100
0.034
DG24
-0.024
-0.054
0.126
0.123
DG34
-0.022
-0.006
0.070
0.032
DGIJK
-0.054
-0.012
0.165
0.137
PLNG
-0.060
0.024
0.052
0.053
PT06
0.121
-0.057
0.047
0.046
PT11
0.060
0.047
0.045
0.027
TTG-1304
0.041
0.004
0.028
0.034
TTG-1305
0.032
0.057
0.032
0.047
TTG-1307
-0.055
0.070
0.078
0.040
VK09
-0.004
0.051
0.054
0.056
VK13
0.062
0.040
0.036
0.021
VK14
-0.479
-0.036
0.015
0.008
VK23
0.171
0.057
0.041
0.029
Pada tabel di atas, nilai negatif (-) menunjukkan bahwa lokasi tersebut mengalami penurunan tanah (land subsidence), sedangkan nilai positif (+) menunjukkan bahwa daerah tersebut mengalami kenaikan tanah (land uplift). Hasil dari pengamatan titik yang mengalami penurunan terkecil adalah -0.060 m, yaitu pada titik PLNG serta terjadi penaikan tanah terbesar sebesar 0.171 m yaitu pada titik VK23. Untuk kala 2 mengalami penurunan tanah terbesar pada titik 37
BT03 yaitu sebesar -0.091 m serta penaikan terbesar pada titik VK14 yaitu sebesar 0.093 m.
4.2 Uji T-Student Hasil perhitungan uji statistik dengan menggunakan uji t-student dapat dilihat dalam table dibawah ini. Tabel 4.5 Tabel standar deviasi Titik
Kala 1 Kala 2 Kala 3
BPNRI-53
0.019
0.010
0.089
BT01
0.017
0.038
0.019
BT03
0.106
0.055
0.05
BW13
0.086
0.045
0.009
DG09
0.016
0.198
0.03
DG11
0.032
0.043
0.019
DG18
0.026
0.097
0.022
DG24
0.122
0.032
0.018
DG34
0.019
0.067
0.026
DGIJK
0.126
0.106
0.055
PLNG
0.05
0.013
0.019
PT06
0.042
0.022
0.019
PT11
0.025
0.038
0.011
TTG-1304
0.012
0.025
0.032
TTG-1305
0.028
0.016
0.038
TTG-1307
0.039
0.068
0.008
VK09
0.052
0.015
0.021
VK13
0.018
0.031
0.010
VK14
0
0.015
0.008
VK23
0.028
0.030
0.009
38
Pada pengujian statistik ini digunakan pengujian dengan tingkat kepercayaan 95% (α =5%), dan 99% (α = 1%). Dari masing-masing tingkat kepercayaan tersebut secara berurutan didapatkan nilai tdf, 0,025 sebesar 2.09302, dan tdf, 0,005 sebesar 2.86093. Dari nilai tdf diatas, didapat nilai rentang X1 = -0.0151 dan X2 = 0.0251 untuk kala 1, X1 = -0.0175 dan X2 = 0.0275 untuk kala 2, dan X1 = -0.0069 dan X2 = 0.0169 tingkat kepercayaan 95%. Dengan nilai rentang yang didapat, maka semua data dapat diterima. Pada tingkat kepercayaan 99%, didapatkan nilai rentang X1 = -0.0226 dan X2 = 0.0326 untuk kala 1, X1 = -0.0258 dan X2 = 0.0358 untuk kala 2, X1 = -0.0113 dan X2 = 0.0213 untuk kala 3. Dengan nilai rentang yang didapat, maka semua data dapat diterima.
4.3 Analisa Regresi pada Setiap Pengukuran a. BPNRI-53
Gambar 4.2 Grafik regresi di titik BPNRI-53
39
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang curam dan pada kala berikutnya mengalami subsidence tidak curam dan hampir mengalami perbedaan yang tidak signifikan dari kala 1 dan kala 3 sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun dengan sangat tajam, hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua. b. BT01
Bulan
Gambar 4.3 Grafik regresi di titik BT01 Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3. Pada regresi polinomial ordo 2 juga menunjukkan bahwa titik tersebut mengalami penaikkan.
40
c. BT03
Bulan
Gambar 4.4 Grafik regresi di titik BT03 Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami Subsidence. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3 mengalami pernurunan yang relatif linier. Pada regresi polinomial ordo 2 juga menunjukkan bahwa titik tersebut mengalami penurunan. Dilihat dari kurva diatas, perubahan pada tiap kala relatif linier sehingga antara kurva regresi linier dan polinomial hampir berimpit.
d. BW13
Bulan
Gambar 4.5 Grafik regresi di titik BW13 41
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3. Terjadi penaikkan terbesar pada saat pengukuran kala kedua dari kala pertama. Sedangkan kala ketiga mengalami penaikkan tidak sebesar kala kedua. e. DG09
Bulan
Gambar 4.6 Grafik regresi di titik DG09 Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3. Terjadi penaikkan terbesar pada saat pengukuran kala kedua dari kala pertama. Sedangkan kala ketiga mengalami penaikkan tidak sebesar kala kedua.
42
f. DG11
Bulan
Gambar 4.7 Grafik regresi di titik DG11 Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3. Terjadi penaikkan terbesar pada saat pengukuran kala kedua dari kala pertama. Sedangkan kala ketiga mengalami penaikkan tidak sebesar kala kedua.
g. DG18
Gambar 4.8 Grafik regresi di titik DG18
43
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami Uplift. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3 mengalami penaikkan yang relatif linier. Pada regresi polinomial ordo 2 juga menunjukkan bahwa titik tersebut mengalami penurunan. Dilihat dari kurva diatas, perubahan pada tiap kala relatif linier sehingga antara kurva regresi linier dan polinomial hampir berimpit.
h. DG24
Gambar 4.9 Grafik regresi di titik DG24
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami Subsidence. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3 mengalami pernurunan yang relatif linier. Pada regresi polinomial ordo 2 juga menunjukkan bahwa titik tersebut mengalami penurunan. Dilihat dari kurva diatas, perubahan pada tiap kala relatif linier sehingga antara kurva regresi linier dan polinomial hampir berimpit.
44
i. DG34
Gambar 4.10 Grafik regresi di titik DG34
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami subsidence. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami subsidence yang curam dan pada kala berikutnya mengalami uplift yang cukup curam dan hampir mengalami perbedaan yang tidak signifikan dari kala 1 dan kala 3 sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva naik dengan sangat tajam, hal ini dikarenakan titik pada kala ketiga mengalami penaikkan dari titik kedua.
45
j. DGIJK
Bulan
Gambar 4.11 Grafik regresi di titik DGIJK
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami Subsidence. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3 mengalami pernurunan yang relatif linier. Pada regresi polinomial ordo 2 juga menunjukkan bahwa titik tersebut mengalami penurunan. Dilihat dari kurva diatas, perubahan pada tiap kala relatif linier sehingga antara kurva regresi linier dan polinomial hampir berimpit.
46
k. PLNG
Bulan
Gambar 4.12 Grafik regresi di titik PLNG
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang curam dan pada kala berikutnya mengalami subsidence tidak curam sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua.
47
l. PT06
Gambar 4.13 Grafik regresi di titik PT06
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami subsidence. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami subsidence yang curam dan pada kala berikutnya mengalami uplift tidak curam sehingga menyebabkan kurva regresi linier turun. Pada grafik polinomial, kurva naik hal ini dikarenakan titik pada kala ketiga mengalami penaikkan dari titik kedua.
48
m. PT11
Gambar 4.14 Grafik regresi di titik PT11 Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang curam dan pada kala berikutnya mengalami subsidence tidak curam sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua.
49
n. TTG 1304
Gambar 4.15 Grafik regresi di titik TTG1304
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang curam dan pada kala berikutnya mengalami subsidence tidak curam sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua.
50
o. TTG1305
Gambar 4.16 Grafik regresi di titik TTG1305
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang curam dan pada kala berikutnya mengalami subsidence tidak curam sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua.
51
p. TTG1307
Gambar 4.17 Grafik regresi di titik TTG1307
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami Subsidence. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3 mengalami pernurunan yang relatif linier. Pada regresi polinomial ordo 2 juga menunjukkan bahwa titik tersebut mengalami penurunan. Dilihat dari kurva diatas, perubahan pada tiap kala relatif linier sehingga antara kurva regresi linier dan polinomial hampir berimpit.
52
q. VK09
Gambar 4.18 Grafik regresi di titik VK09
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini disebabkan karena pada pengukuran kala ketiga, nilai tingginya sangat jauh berbeda dengan kala sebelumnya yang mana pada kala pertama mengalami subsidence yang kecil dan pada kala berikutnya mengalami uplift yang besar sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva naik hal ini dikarenakan titik pada kala ketiga mengalami penaikkan dari titik kedua.
53
r. VK13
Gambar 4.19 Grafik regresi di titik VK13
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplif dan subsidence. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang besar dan pada kala berikutnya mengalami subsidence kecil sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua..
54
s. VK14
Bulan
Gambar 4.20 Grafik regresi di titik VK14 Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplift. Hal ini bisa dilihat pada pengukuran mulai kala 1, kala 2 dan kala 3. Terjadi penaikkan terbesar pada saat pengukuran kala kedua dari kala pertama. Sedangkan kala ketiga mengalami penaikkan tidak sebesar kala kedua. t. VK23
Gambar 4.21 Grafik regresi di titik VK23 55
Dari diagram di atas bisa dilihat bahwa pada regresi linier, pengukuran tiap kala mengalami uplif dan subsidence. Hal ini disebabkan karena pada pengukuran kala pertama, nilai tingginya sangat jauh berbeda dengan kala berikutnya yang mana pada kala pertama mengalami uplift yang besar dan pada kala berikutnya mengalami subsidence kecil sehingga menyebabkan kurva regresi linier naik. Pada grafik polinomial, kurva menurun hal ini dikarenakan titik pada kala ketiga mengalami penurunan dari titik kedua.
Gambar 4.22 Hasil sebaran titik beserta keadaan tanaah selama pengamatan
Dari pemetaan diatas lokasi yang berdekatan dengan semburan lumpur dengan radius kurang dari 2 km mengalami subsidence dan uplift. Namun dari lokasi tersebut ada titik yang mengalami subsidence yaitu arah timur laut dari semburan / pada titik 1307 dan PT06. Ini disebabkan karena titik tersebut
56
merupaka kawasan dari jalur sesar watukosek. Menurut Mazzini (2007), Situasi LUSI pada busur belakang berjarak 10 km arah timur laut dari gunung Penanggungan. Penampang seismik dan pengamatan lapangan menunjukkan adanya sesar regional (yang dimaksud adalah Sesar Watukosek) melalui daerah Lusi. Dari penelitian ini ada beberapa perbedaan dengn hasil penelitian terdahulu. Hasil dari penelitian terdahulu bisa dilihat pada gambar berikut.
Gambar 4.23 Peta daerah yang mengalami uplift dan subsidence berdasarkan penelitian terdahulu
Dari gambar di atas bisa dilihat bahwa daerah yang mengalami subsidence adalah di sekitar semburan dan sebelah barat pusat semburan, sedangkan daerah yang mengalami uplift berada di sebelah timur laut pusat semburan. Sehingga dari penelitian ini terdapat beberapa perbedaan pada lokasi yang mengalami subsidence maupun uplift. Karena alat yang digunakan dalam penelitian berbeda, penelitian terdahulu menggunakan Insar, sedangkan penelitian ini menggunakan GPS.
57
Menurut Prawoko, Eko Yuli, dan Taufik (2008) melakukan penelitian tentang penurunan tanah disekitar Lumpur Lapindo dengan menggunakan GPS dengan metode analisa arah gerakan tanah dan memprediksi perubahan tinggi titik pada tiap kala. Hasil dari pengukuran GPS menunjukkan pergeseran vertikal yang terjadi rata-rata antara -0,1275 m sampai dengan 0,2574 m. Hutomo, Eko Yuli, (2010) melakukan studi tentang pemanfaatan CORS GPS milik ITS dalam hal ini di gedung Teknik Geomatika. Penelitian ini memanfaatkan CORS GPS sebagai alternatif stasiun referensi dalam pengukuran GPS. Nilai penurunan tanah terbesar adalah sebesar -0.018 m di titik BT03 dan nilai kenaikan tanah terbesar adalah sebesar 0.0115 m di titik 1304. Nilai korelasi dan determinasi yang terbesar adalah pada regresi polinomial ordo 3 yaitu dengan nilai korelasi rata-rata sebesar 0.893 dan nilai determinasi sebesar 0.810. Jika dibandingkan dengan penelitian ini dan penelitian sebelumnya, terjadi perubahan yang sangat signifikan, berikut tabel perubahan tiap titik. Tabel 4.6 Perbandingan dengan penelitian terdahulu no 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
posisi BPNRI-53 BT01 BT03 BW13 DG09 DG11 DG18 DG24 DG34 DGIJK PT06 PLNG PT11 TTG-1304 TTG-1305 TTG-1307 VK09 VK13 VK14 VK23
Maret 2014 - Mei 2011 (m) -0.479 -0.294 -0.468 -0.581 -0.603 -0.502 -0.442 -0.388 -0.397 -0.343 -0.497 -0.058 -
58
Perbedaan ini terjadi karena beberapa hal, antara lain: Hal-hal yang mempengaruhi penelitian ini juga berbeda. Penelitian ini menambahkan titik pengamatan dari penelitian sebelumnya, sehingga dalam tabel diatas data yang bisa dibandingkan hanya beberapa titik. Waktu dilakukannya pengamatan tidak sama dan terpaut jauh yaitu 3 tahun. Pada penelitian terdahulu dilakukan pada bulan Mei tahun 2014, sehingga sudah banyak terjadi perubahan.
59
“Halaman ini sengaja dikosongkan”
60
BAB V KESIMPULAN DAN SARAN
5.1 Kesimpulan Hasil dari pengamatan dalam penelitian ini mendapatkan beberapa hasil yang dapat disimpulkan, diantaranya: 1. Titik kontrol yang berjumlah 20 titik mengalami penurunan (subsidence) dan penaikkan (uplift) setelah dibuktikan dengan pengamatan GPS. 2. Titik yang mengalami penurunan ada 6 titik diantaranya adalah; BT03, DG24, DG34, DGIJK, PT06, dan TTG-1307 sedangkan yang mengalami penaikkan ada 14 titik adalah BPNRI-53, BT01, BW13, DG09, DG11, DG18, PLNG, PT11, TTG1304, TTG1305, VK09, VK13, VK14, dan VK23. 3. Penurunan terkecil adalah -0.060 m, yaitu pada titik PLNG terjadi pada kala kedua, serta terjadi penaikan tanah tertinggi sebesar 0.171 m yaitu pada titik VK23. Hasil yang kurang valid didapat pada VK14, yaitu mengalami penurunan hampir 50 cm. hal ini didapat karena terjadi kesalahan teknis pada receiver dan CORS sebagai base station. Sedangkan kala ketiga terjadi penurunan terkecil pada titik BT03 yaitu -0.057 m dan penaikkan terbesar terjadi pada titik VK14 dan DG18 sebesar 0.076 m.
5.2 Saran 1. Pengukuran menggunakan lebih banyak reveicer dan base station yang lebih dekat dari daerah pengukuran. 2. Pengecekan peralatan survey harus dilakukan secara berkala, agar tetap dapat digunakan dan mendapat hasil yang baik. 3. Penelitian ini harus dikorelasi dengan data pendukung yang lain sebagai pembanding karena nilai perambatan kesalahannya lebih besar..
61
“Halaman ini sengaja dikosongkan”
62
DAFTAR PUSTAKA
Abidin, H.Z. 2007. Penentuan Posisi Dengan GPS dan Aplikasinya. P.T. Pradnya Paramita, Jakarta. Second edition. ISBN 979-408-377-1. 268 pp Abidin, H.Z., Andreas,H., Kusuma,M.A., Sumintadireja,P., Gumilar,I. 2010. Ground Displacements around LUSI Mud Volcano Indonesia as Inferred from GPS Surveys. FIG Congress: Sydney, Australia Abidin, H.Z., Davies, R. J., Kusuma, M.A., Andreas, H., Deguchi, T. 2008. Subsidence and uplift of Sidoarjo (East Java) due to the eruption of the Lusi mud volcano (2006–present) Abidin, H.Z., Jones, A., Kahar, J. 2002. Survey dengan GPS. Jakarta: PT Pradnya Paramita Abidin, H.Z, Andreas.H, Gumilar.I, Sidiq.I.P, Gamal.M, Murdohardono.D, Supriyadi, Fukuda.Y. 2010. Studiying Land Subsidence in Semarang (Indonesia) Using Geodetics Methods. FIG Congress April 11-16 Andreas, H., Abidin, H. Z., Kusuma, M. A., Sumintadireja, P, Gumilar, I., 2010. Ground Displacement around LUSI Mud Volcano Indonesia as Inffered from GPS Surveys. Presented at the FIG Congress 2010 (Facing the Challenge-Building Capacity), Sydney Australia. Bakti, L.M., 2010. Kajian Sebaran Potensi Rob Kota Semarang dan Usulan Penanganannya.
Program
Magister
Teknik
Sipil,
Universitas
Diponegoro Semarang. Blewitt, Geoffrey. (1998), GPS Data Processing Methodology: From Theory To Applications, GPS for Geodesy, Hal. 231-270, Springer-Verlag, Berlin Bukhori, I., 2011. Model Estimasi Uplift dan Subsidence dari Hasil Ukuran GPS Menggunakan Metode Polinomial di Area Lumpur Sidoarjo. Program Studi Teknik Geomatika ITS), Kampus ITS Sukolilo, Surabaya 60111. Caspary, W. F. 1987. Concept of Network and Deformation Analysis. School of Surveying The University of New South Wales, Kensington, N.S.W, Australia.
63
Erif, L. O. M., Wirada, F. 2011. Tinjauan Kritis Terhadap Pemerintah dalam Penanganan Kasus Lumpur Lapindo. Sekolah Pascasarjana Universitas Gajah Mada. Yogyakarta. Fukushima, Y., Mori, J., Hashimoto, M., Kano, Y., 2009. Subsidence associated with the LUSI mud eruption, East Java, investigated by SAR interferometry. Marine and Petroleum Geology 26, 1740–1750. Gili, Josep, A., Corominas, Jordi., Rius, Joan., (2000), Using Global Positioning System Techniques in Landslide Monitoring, Engineering Geology, Vol. 55, Hal. 167-192. Goad, C.C., and L. Goodman (1974), A modified Hopfield tropospheric refraction correction model. In: American Geophysical Union Annual Fall Meeting,
12-17
December
1974,
San
Francisco,
California,
USA(abstract EOS Trans. AGU 55, 1106). Guntoro, 2007. Sidoarjo Mudflow Hypothesis in Geological-Expert. Discourse. LUSI MEDIA CENTER, Sidoarjo. Hartoko, A., Wirasatria,A., Helmi,M., dan Rochaddi,B., 2008. Aplikasi Teknologi Geomatika untuk Pemetaan Penurunan Tanah (Land Subsidence) di Pesisir Kota Semarang. Pusat Database Spasial Pusat Studi Geomatika Perikanan dan Kelautan, Fak. Perikanan dan llmu Kelautan – UNDIP, ILMU KELAUTAN. Vol. 13 (I): 13 -18. Hofmann-Wellenhof, Bernhard and Moritz, Helmut., 2006. Physical Geodesy second edition, Springer Wien Newyork. Kuang, S., 1996. Geodetic network analysis and optimal design: concepts and applications. Ann Arbor Press, Chelsea, Mich Kurniawan, A. 2011. Evaluasi PenurunanMuka Tanah Di Wilayah Kota Surabaya Dari Pengamatan Global Positioning System Dengan GAMIT/GLOBK. Program Studi Teknik Geomatika, ITS Kurniawan, D. 2008. Regresi Linier. Vienna, Austria: R Development Core Team Leonard, F. R., 2000. Penentuan Kecepatan dan Percepatan Penurunan Muka Tanah Menggunakan Data Beda Tinggi Departemen Geodesi ITB, Bandung.
64
Jaring Sipat Datar.
Mazzini, A., Svensen,H., Akhmanov,G.G., Aloisi,G., Planke,S., Sorenssen,A.M., Istiadi,B. 2007. Triggering and dynamic evolution of the LUSI mud volcano, Indonesia :Earth and Planetary Science Letters Prawoko, A,A. 2008 Pemantauan Land Subsidance untuk Daerah Sekitar Semburan Lumpur Sidoarjo dengan Menggunakan GPS, Program Studi Teknik Geomatika, ITS Satrio dan Pratikno, B. Studi Air Tanah pada Area Semburan Lumpur Lapindo Sidoarjo Menggunakan Isotop Alam. Pusat Aplikasi Teknologi Isotop dan Radiasi & BATAN Sophian, R. I., 2010. Penurunan Muka Air Tanah di Kota - Kota Besar Pesisir Pantai Utara Jawa Tengah (Studi Kasus : Kota Semarang). Bulletin of Scientific Contribution, 4 1-60 8 Nomor 1 Sudarsono, U dan Sudjarwo, I.B. 2008. Aspek Geologi Teknik Lumpur Sidoarjo, Jawa Timur. Bandung: Pusat Lingkungan Geologi Verhoef, H.M.E, dan de Heus H, M. 1994. On The Estimation of Polynomial Breakpoints in the Subsidance of the Groningen Gasfield. Belanda Widodo, A. 2007. Memahami Bencana Gunung Lumpur kasus Lumpur Sidoarjo. Surabaya: itspress Wijaya. D. S. A., 2013. A Review on Hazard Risk Reduction Systems and Reliability Estimate of the Dredging System of the Lusi Mud Volcano in Sidoarjo, East Java. International Journal of Civil & Environmenta l Engineering IJCEE -IJENS Vol: 13 No: 02. Whittaker B. N. & Reddish D. J. 1989, Subsidence Occurrence, Prediction and Control,. Developments in Geotechnical Engineering, 56, Elsevier.
65
“Halaman ini sengaja dikosongkan”
66
LAMPIRAN
1. Raw data koordinat
CORS NO
Tanggal Pengamatan
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
3/18/2014 9:26 3/18/2014 9:51 3/18/2014 14:52 3/18/2014 15:01 3/19/2014 9:29 3/19/2014 15:12 3/19/2014 16:05 3/20/2014 9:18 3/20/2014 9:39 3/20/2014 14:33 3/20/2014 15:13 3/24/2014 9:37 3/24/2014 9:44 3/24/2014 15:12 3/24/2014 15:19 3/25/2014 9:08 3/25/2014 14:31 3/26/2014 9:42 3/26/2014 15:09 3/26/2014 18:59
POSISI TTG‐1307 BT01 VK23 PT11 PLNG BT03 DGIJK BW13 DG18 TTG‐1305 DG34 DG11 PT06 BPNRI‐53 DG24 VK09 VK14 VK13 TTG‐1304 DG09
KALA 1 18‐20 Maret & 24‐26 Maret 2014 9194906.16 698128.921 47.951 Koordinat Standart Deviasi tinggi Y (m) X(m) Z(m) antena(m) N(m) E(m) U(m) 9171609.592 688853.560 32.399 1.185 0.012 0.038 0.039 9170380.973 687275.599 32.285 1.161 0.007 0.029 0.017 9169335.059 688400.679 32.771 1.14 0.014 0.041 0.028 9171031.650 685757.177 34.142 0.675 0.011 0.037 0.025 9168181.325 688091.594 31.408 1.214 0.014 0.044 0.042 9168700.731 687447.326 32.727 1.323 0.046 0.115 0.106 9168341.289 684425.374 33.903 1.032 0.06 0.136 0.126 9166787.789 687841.592 33.482 1.24 0.029 0.071 0.086 9165896.403 684993.558 35.934 1.12 0.009 0.011 0.026 9165632.400 687353.520 38.992 1.195 0.014 0.044 0.028 9166171.182 683191.065 36.447 1.041 0.01 0.03 0.019 9168770.289 690071.219 31.347 1.22 0.014 0.038 0.032 9170158.495 692979.909 30.497 1.217 0.026 0.051 0.05 9168497.123 691251.778 30.602 1.334 0.01 0.038 0.019 9169010.760 692360.482 31.123 1.128 0.063 0.089 0.122 9167321.956 690298.337 31.733 1.095 0.027 0.052 0.052 9165749.683 692869.647 36.890 1.175 0 0 0 9165694.08 688738.846 37.762 1.138 0.007 0.009 0.018 9162606.096 686566.847 43.258 1.168 0.006 0.025 0.012 9165308.179 689811.42 34.803 1.307 0.005 0.014 0.016
67
CORS NO
Tanggal Pengamatan
POSISI
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
4/21/2014 8:57 4/21/2014 9:25 4/21/2014 14:16 4/21/2014 14:47 4/22/2014 8:52 4/22/2014 14:18 4/22/2014 14:23 4/23/2014 8:40 4/23/2014 9:15 4/23/2014 14:50 4/24/2014 9:34 4/24/2014 13:57 4/24/2014 14:48 4/25/2014 7:31 4/25/2014 8:36 4/25/2014 13:00 4/25/2014 14:09 4/26/2014 8:00 4/26/2014 13:24 4/26/2014 18:29
BPNRI‐53 BT01 BT03 BW13 DG09 DG11 DG18 DG24 DG34 DGIJK PLNG PT06 PT11 TTG‐1304 TTG‐1305 TTG‐1307 VK09 VK13 VK14 VK23
KALA 2 21 APRIL ‐ 26 APRIL 2014 9194906.16 698128.921 47.951 Koordinat Standart Deviasi tinggi Y (m) X(m) Z(m) antena(m) N(m) E(m) U(m) 9168497.102 691251.551 30.685 1.185 0.006 0.016 0.011 9170380.942 687275.739 32.373 0.675 0.04 0.117 0.087 9168700.255 687447.343 32.696 1.145 0.019 0.035 0.042 9166787.838 687841.623 33.534 1.161 0.036 0.081 0.067 9165308.129 689811.409 34.810 1.22 0.094 0.168 0.198 9168770.551 690071.248 31.435 1.247 0.026 0.079 0.043 9165896.083 684993.663 35.997 1.402 0.048 0.085 0.097 9169010.962 692360.453 31.099 1.104 0.056 0.118 0.138 9166171.198 683191.130 36.425 1.324 0.027 0.083 0.049 9168341.123 684425.313 33.849 1.439 0.059 0.115 0.104 9168181.203 688091.667 31.348 1.082 0.004 0.018 0.011 9170158.541 692979.948 30.618 1.301 0.027 0.048 0.059 9171031.623 685757.094 34.202 1.1 0.105 0.204 0.238 9162606.096 686566.786 43.299 1.051 0.009 0.031 0.024 9165632.518 687353.247 39.024 1.254 0.039 0.08 0.069 9171609.568 688853.618 32.344 1.13 0.007 0.024 0.015 9167322.001 690298.498 31.729 1.112 0.005 0.024 0.012 9165694.077 688738.842 37.824 1.14 0.025 0.048 0.062 9165749.049 692869.567 36.411 1.093 0.007 0.021 0.015 9169334.688 688400.544 32.942 1.291 0.012 0.021 0.02
68
NO 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
KALA 3 26 MEI ‐ 29 MEI 2014 CORS 9194906.16 698128.921 Koordinat Tanggal POSISI Pengamatan Y (m) X(m) 5/26/2014 18:10 TTG‐1307 9171609.662 688853.599 5/26/2014 18:52 PT11 9171031.671 685757.001 5/26/2014 23:27 VK23 9169335.014 688400.652 5/27/2014 0:31 TTG‐1305 9165632.517 687353.509 5/27/2014 6:12 PLNG 9168181.209 688091.585 5/27/2014 11:36 BT01 9170381.008 687275.696 5/27/2014 21:41 BT03 9168700.867 687447.333 5/28/2014 3:12 VK14 9165749.571 692869.543 5/28/2014 7:34 DG11 9168770.556 690071.213 5/28/2014 9:03 VK09 9167321.993 690298.528 5/28/2014 12:57 DG24 9169011.047 692360.501 5/28/2014 14:32 VK13 9165694.084 688738.843 5/28/2014 20:07 BPNRI‐53 9168497.144 691251.721 5/29/2014 0:55 BW13 9166787.841 687841.622 5/29/2014 2:00 TTG‐1304 9162606.102 686566.79 5/29/2014 7:17 PT06 9170158.553 692979.926 5/29/2014 19:21 DG09 9165308.324 689811.332 5/30/2014 7:23 DG18 9165896.414 684993.592 5/30/2014 7:25 DG34 9166171.210 683191.125 5/30/2014 12:49 DGIJK 9168341.429 684425.414
69
47.951 Z(m) 32.292 34.189 32.852 39.049 31.306 32.400 32.605 36.504 31.464 31.784 31.069 37.802 30.633 33.535 43.262 30.569 34.823 36.090 36.441 33.781
tinggi antena(m) 1.190 0.703 1.230 1.173 1.297 1.181 1.416 1.181 1.230 1.143 1.243 1.206 1.342 1.139 1.096 1.347 1.126 1.201 1.236 1.156
Standart Deviasi N(m) E(m) U(m) 0.004 0.015 0.008 0.005 0.02 0.011 0.007 0.017 0.009 0.013 0.018 0.038 0.006 0.008 0.019 0.019 0.010 0.019 0.016 0.021 0.05 0.006 0.018 0.008 0.019 0.010 0.019 0.006 0.007 0.021 0.007 0.016 0.018 0.008 0.018 0.010 0.037 0.067 0.089 0.007 0.018 0.009 0.014 0.037 0.032 0.005 0.005 0.019 0.01 0.012 0.03 0.022 0.011 0.022 0.025 0.013 0.026 0.053 0.044 0.055
2. Tabel uji T student
70
LAMPIRAN GRAFIK
Perubahan tiap bulan titik BPNRI53 30.700 30.680 30.660 30.640 30.620 30.600 30.580 30.560 30.540 30.520 Maret
April
Mei
Perubahan tiap bulan titik BT01 32.380 32.360 32.340 32.320 32.300 32.280 32.260 32.240 Maret
April
Mei
Perubahan tiap bulan titik DGIJK 33.910 33.900 33.890 33.880 33.870 33.860 33.850 33.840 33.830 33.820 Maret
April
Mei
71
Perubahan tiap bulan titik PLNG 31.6 31.4 31.2 31 30.8 30.6 30.4 30.2 30 Maret
April
Mei
Perubahan tiap bulan titik BT03 32.730 32.725 32.720 32.715 32.710 32.705 32.700 32.695 32.690 32.685 32.680 Maret
April
Mei
Perubahan tiap bulan titik BW13 33.540 33.530 33.520 33.510 33.500 33.490 33.480 33.470 33.460 33.450 Maret
April
Mei
72
Perubahan tiap bulan titik DG09 34.825 34.82 34.815 34.81 34.805 34.8 34.795 34.79 Maret
April
Mei
Perubahan tiap bulan titik DG11 31.460 31.440 31.420 31.400 31.380 31.360 31.340 31.320 31.300 Maret
April
Mei
Perubahan tiap bulan titik DG24 31.130 31.120 31.110 31.100 31.090 31.080 31.070 31.060 31.050 31.040 Maret
April
Mei
73
Perubahan tiap bulan titik DG34 36.450 36.445 36.440 36.435 36.430 36.425 36.420 36.415 36.410 Maret
April
Mei
Perubahan tiap bulan titik PT06 31.600 31.400 31.200 31.000 30.800 30.600 30.400 30.200 30.000 Maret
April
Mei
Perubahan tiap bulan titik PT11 34.21 34.2 34.19 34.18 34.17 34.16 34.15 34.14 34.13 34.12 34.11 Maret
April
Mei
74
Perubahan tiap bulan titik TTG‐1304 43.31 43.3 43.29 43.28 43.27 43.26 43.25 43.24 43.23 Maret
April
Mei
Perubahan tiap bulan titik TTG‐1305 39.060 39.050 39.040 39.030 39.020 39.010 39.000 38.990 38.980 38.970 38.960 Maret
April
Mei
Perubahan tiap bulan titik TTG‐1307 32.480 32.460 32.440 32.420 32.400 32.380 32.360 32.340 32.320 32.300 32.280 Maret
April
Mei
75
Perubahan tiap bulan titik DG18 36.020 36.000 35.980 35.960 35.940 35.920 35.900 35.880 Maret
April
Mei
Perubahan tiap bulan titik VK09 31.790 31.780 31.770 31.760 31.750 31.740 31.730 31.720 31.710 31.700 Maret
April
Mei
Perubahan tiap bulan titik 13 37.83 37.82 37.81 37.8 37.79 37.78 37.77 37.76 37.75 37.74 37.73 Maret
April
Mei
76
Perubahan tiap bulan titik VK14 37.000 36.900 36.800 36.700 36.600 36.500 36.400 36.300 36.200 36.100 Maret
April
Mei
Perubahan tiap bulan titik DG18 33.000 32.950 32.900 32.850 32.800 32.750 32.700 32.650 Maret
April
Mei
77
LAMPIRAN FOTO
1. Fototitikpengamatan
TTG-1307
BW-13
VK-14
TTG-1304
PT-11
TTG-1305
78
PT-06
DG-18 VK-09
VK-23
DGIJK
79
DG-34
DG-11
DG-24
DG-09
VK-13
80
BIODATA PENULIS
Wisnu Pribadi merupakan anak pertama dari dua bersaudara, dilahirkan di Kota Ponorogo, 24 Agustus 1989. Penulis telah menempuh pendidikan formal di TK Arrohman 3 Kabupaten Ponorogo, SD Muhammadiyah Kabupaten Ponorogo, SMP Negeri 1 Kabupaten Ponorogo, SMA Negeri 1 Kabupaten Ponorogo dan memilih untuk melanjutkan kuliah S-1 dengan mengikuti jalur PMDK dan diterima di Jurusan Teknik Geomatika – FTSP ITS Surabaya dengan NRP 35 08 100 027 dan melanjutkan studi S-2 di Teknik Geomatika – FTSP melalui beasiswa fresh graduate dengan NRP 35 12 201 006. Selama menjadi mahasiswa penulis aktif di berbagai kegiatan organisasi di kampus ITS. Penulis memilih bidang keahlian geospasial untuk menyelesaikan studi S-1 sebelumnya di bidang penginderaan jauh dengan judul “Inventarisasi Lahan Pertanian Di Kabupaten Sidoarjo Menggunakan Citra Satelit Multitemporal”. Pada penelitian untuk studi S-2, penulis memilih bidang keahlian manajemen survey dan pemetaan dengan judul “Analisa Penurunan Tanah (Landsubsidence) Pada Daerah Semburan Lumpur Sidoarjo Dengan Menggunakan Data Satelit Global Positioning System (GPS)”. Dalam hidupnya, penulis berharap selalu bisa menjadi orang yang bermanfaat bagi orang lain dan lingkungan sekitar. Motto hidup penulis adalah kecerdasan seseorang bukan penentu utama suatu kesuksesan melainkan doa dan ikhtiar, keuletan dan sikap merupakan kunci utama. Setiap orang jika mempunyai modal tersebut pasti dapat meraih apa yang diinginkan serta apa yang dicita-citakan.