ANALISIS INDIKATOR KUANTITATIF FUNGSI HIDROLOGI AKIBAT ALIH GUNA LAHAN DAS CILIWUNG HULU
APRIAN PURNOTO
DEPARTEMEN GEOFISIKA DAN METEOROLOGI FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT PERTANIAN BOGOR 2007
ABSTRAKSI APRIAN PURNOTO. Analisis Indikator Kuantitatif Fungsi Hidrologi Akibat Alih Guna Lahan DAS Ciliwung Hulu. Dibimbing oleh Bambang Dwi Dasanto. Daerah Aliran Sungai (DAS) Ciliwung merupakan salah satu DAS yang penting untuk mendapat perhatian, karena DAS ini mempunyai topografi yang cukup beragam dengan daerah hulu yang berbukit – bukit dan daerah hilir yang relatif datar yang berpeluang besar menyebabkab bencana banjir. Oleh karena itu pengamanan DAS ini perlu mendapat perhatian mulai dari kawasan hulu sampai hilir. Penelitian ini bertujuan untuk mengevaluasi fungsi hidrologi DAS di kawasan hulu dengan pendekatan kuantitatif selama periode 1993 – 2005 dan usaha untuk mengurangi debit banjir dengan menerapkan skenario alih guna lahan menggunakan model GenRiver. Model simulasi hidrologi GenRiver telah digunakan untuk mempelajari perubahan aliran sungai sebagai akibat adanya alih guna lahan, dan selanjutnya dipakai sebagai dasar untuk mempelajari beberapa skenario alih guna lahan yang mungkin terjadi di masa yang akan datang. GenRiver adalah sebuah model yang berbasis pada proses neraca air pada skala plot dan dikembangkan menjadi skala sub-DAS. Aplikasi simulasi model GenRiver untuk daerah Ciliwung hulu menggunakan skenario pertama alih guna lahan yaitu merubah 9,1% tegalan dan 4,3% perkebunan menjadi lahan hutan sehingga luas hutan menjadi 35%. Skenario kedua merubah 17,8% perkebunan, dan 7,3% sawah menjadi lahan hutan, serta 9,1% tegalan menjadi hutan dan 1% menjadi pemukiman, sehingga luas hutan menjadi 55,8%. Hasil analisis menunjukkan bahwa dengan curah hujan normal dan variasi curah hujan ± 25% komposisi lahan 55,8% hutan, perkebunan 13,8%, pemukiman 26,8%, tegalan 1,5%, dan sawah 2,1% secara nyata dapat menurunkan debit maksimum tetapi kurang nyata pada peningkatan debit minimum. Dari analisis sensitivitas parameter, peubah curah hujan memberikan perubahan luaran yang paling besar disusul peubah penggunaan lahan hutan dan perkebunan.
ANALISIS INDIKATOR KUANTITATIF FUNGSI HIDROLOGI AKIBAT ALIH GUNA LAHAN DAS CILIWUNG HULU
APRIAN PURNOTO
Skripsi sebagai salah satu syarat untuk memperoleh gelar Sarjana Sains pada Fakultas Matematika dan Ilmu Pengetahuan Alam Institut Pertanian Bogor
DEPARTEMEN GEOFISIKA DAN METEOROLOGI FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT PERTANIAN BOGOR 2007
JUDUL : Analisis Indikator Kuantitatif Fungsi Hidrologi Akibat Alih Guna Lahan DAS Ciliwung Hulu NAMA : Aprian Purnoto NRP
: G24102040
Menyetujui, Pembimbing
Drs. Bambang Dwi Dasanto, M.Si NIP. 132 014 045
Mengetahui, Dekan Fakultas Matematika dan Ilmu Pengetahuan Alam Insitut Pertanian Bogor
Prof.Dr.Ir.Yonny Koesmaryono, MS NIP. 131 473 999
Tanggal lulus :
RIWAYAT HIDUP Penulis dilahirkan di Purbalingga pada tanggal 22 April 1983 dari pasangan Bapak Supono dan Ibu Warsiyah, merupakan anak ke-5 dari 5 bersaudara. Pada tahun 1990 – 1996 penulis menempuh pendidikan dasar di Madrasah Ibtidaiyah (MI) Kembaran Wetan, dilanjutkan dengan pendidikan lanjutan tingkat pertama di SLTP N 1 Kaligondang. Pada tahun 1999 penulis melanjutkan sekolah menengah umum di SMU N 1 Purbalingga sampai tahun 2002. Pada tahun yang sama penulis lulus dan diterima sebagai mahasiswa di Departemen Geofisika dan Meteorologi IPB melalui jalur Seleksi Penerimaan Mahasiswa Baru (SPMB). Selama kuliah penulis aktif di beberapa kegiatan keorganisasian. Tahun 2004 penulis menjadi panitia bagian publikasi dan dokumentasi pada kegiatan penyambutan mahasiswa baru di Departemen. Pada tahun yang sama penulis juga masuk dalam kepanitiaan ”PESTA SAINS 2004” yang dilaksanakan oleh BEM FMIPA sebagai seksi Publikasi dan Dokumentasi (Pubdok). Tahun 2005 penulis masuk sebagai anggota Departemen Keilmuan dan Keprofesian Himpunan Mahasiswa Geofisika dan Meteorologi (HIMAGRETO). Pada tahun yang sama penulis dipercaya sebagai ketua pelaksana pada kegiatan ”Talk Show PL dan Skripsi” yang diselenggarakan oleh HIMAGRETO. Selain itu penulis juga aktif berpartisipasi pada seminar dan workshop baik didalam maupun diluar Departemen. Tahun 2005 penulis melaksanakan Praktik Kerja Lapang (PKL) di Balai Penelitian Agroklimat dan Hidrologi, Bogor. Pada akhir tahun 2006 penulis dipercaya untuk membantu pengolahan data spasial pada proyek TISDA BPPT Jakarta.
KATA PENGANTAR Puji syukur Alhamdulillah senantiasa penulis panjatkan kapada Allah SWT yang telah memberikan nikmat dan kesehatan yang tiada hentinya kepada umat-Nya. Shalawat dan salam semoga selalu tercurah kepada junjungan kita Nabi Besar Muhammad SAW beserta keluarga, sahabat, dan para pengikutnya. Atas rakhmat Allah Yang Maha Kuasa, penulis dapat melaksanakan skripsi dengan topik “ANALISIS INDIKATOR KUANTITATIF FUNGSI HIDROLOGI AKIBAT ALIH GUNA LAHAN DAS CILIWUNG HULU” yang dilanjutkan dengan penyusunan laporan sebagai salah satu syarat untuk memperoleh gelar sarjana sains pada Departemen Geofisika dan Meteorologi, Fakultas Matematika dan Ilmu Pengetahuan Alam, Institut Pertanian Bogor. Penulis juga mengucapkan terima kasih yang sebesar – besarnya kepada : 1. 2. 3. 4. 5. 6. 7. 8. 9.
Ayah, Ibu, dan Kakak, yang tiada henti – hentinya memberikan dorongan dan kekuatan kepada penulis untuk dapat menyelesaikan skripsi. Bapak Drs. Bambang Dwi Dasanto. M.Si selaku dosen pembimbing yang senantiasa memberikan masukan dan arahan kepada penulis. Bapak Prof. Dr. Hidayat Pawitan, M.Sc dan Bapak I Putu Santikayasa, S.Si, M.Sc selaku dosen penguji yang telah menguji dan memberikan masukan – masukan kepada penulis. Staff TU dan jajarannya yang telah memudahkan penulis dalam mengurus administrasi selama penulis belajar di Departemen Geofisika dan Meteorologi. Pak Pono yang tiada lelah menjaga perpustakaan GFM sehingga penulis dapat meminjam buku kapanpun penulis butuhkan. Pak Udin dan Pak Khairun yang telah membantu penulis dalam mempersiapkan sidang. All crew GFM `39 “I love u so much” (Basyar, Fio, Ani, Kiki, Ode, Noni, Misna, An-an, Hesti, Vivi, Gian, Joko, Lupi, Lina, Nana, Anton, Ipit, Dwi, Linda, Sapta, Samba, Deni, Tado, Mian, Ana, Sasat, Mamat, Rudi, Eko Tas, Nida, Dwinita, Zainul, dan Away). Teman – teman kosan di ”Wisma Sri”: Away, Rudi, Dwi, thanks banget untuk persahabatannya, jujur kalian telah memberikan kenangan – kenangan yang luar biasa, The unforgetable moments!! Adik – adik kelas GFM 40, 41, dan 42 terima kasih telah memberikan warna yang berbeda selama kuliah.
Penulis menyadari masih banyak kekurangan dan keterbatasan pada tulisan ini sehingga kritik dan saran senantiasa penulis nantikan untuk perbakan di masa yang akan datang. Akhirnya semoga tulisan ini memberikan banyak manfaat kepada pembacanya.
Bogor, Februari 2007 Aprian Purnoto G24102040
DAFTAR ISI
hal DAFTAR TABEL ………………………………………………………..viii DAFTAR GAMBAR ......…………………………………………………....viii DAFTAR LAMPIRAN..........................................................................................ix I. PENDAHULUAN ............................................................................1 1.1 Latar Belakang ............................................................................1 1.2 Tujuan .........................................................................................1 II. TINJAUAN PUSTAKA............................................................................1 2.1 Definisi DAS ............................................................................1 2.2 Definsis Penggunaan Lahan ..................................................2 2.3 Jenis Penggunaan Lahan ..................................................4 2.4 Alih Guna Lahan ...............................................................4 2.5 Kriteria dan Indikator Fungsi Hidrologi .....................................4 2.6 Kejadian Banjir di Jakarta ..................................................5 2.7 Deskripsi Singkat Model GenRiver .....................................5 2.8 Model Aplikasi ............................................................................6 III. METODOLOGI ............................................................................7 3.1 Waktu dan Tempat ...............................................................7 3.2 Bahan dan Alat ............................................................................7 3.3 Metodologi 3.3.1 Studi Literatur ...............................................................7 3.3.2 Analisis Alih Guna Lahan .....................................7 3.3.3 Menentukan Kriteria dan Indikator Fungsi Hidrologi ....7 3.3.3.1 Transmisi Air ..................................................8 3.3.3.2 Penyangga pada Puncak Kejadian Hujan ...........8 3.3.3.3 Koefisien Rejim Sungai .....................................8 3.3.4 Simulasi Model GenRiver .....................................8 3.3.5 Verifikasi dan kalibrasi Model ....................................11 3.3.6 Analisis Sensitivitas Model ....................................12 3.3.7 Uji Keabsahan Model .................................................12 3.3.8 Skenario Tutupan Lahan dan Hujan .......................12 IV. KONDISI UMUM ...........................................................................13 4.1 Kondisi Umum DAS Ciliwung Hulu ....................................13 4.2 Iklim ........................................................................................13 4.3 Tanah ........................................................................................13 V. HASIL DAN PEMBAHASAN ..............................................................14 5.1 Analisis Perubahan Penggunaan Lahan ....................................14 5.2 Analisis Curah Hujan ..............................................................15 5.3 Analisis Evapotranspirasi .................................................16 5.4 Analisis Debit Sungai ..............................................................16 5.5 Analisis Indikator Penyangga DAS ....................................17 5.6 Koefisien Rejim Sungai ..............................................................18 5.7 Simulasi Model GenRiver untuk Prediksi Debit .......................19 5.8 Analisis Sensitivitas Tutupan Lahan ....................................20 5.9 Analisis Sensitivitas Curah Hujan ....................................21 6.0 Hasil Prediksi Debit ..............................................................22 VI. KASIMPULAN DAN SARAN ..............................................................23 6.1 Kesimpulan ...........................................................................23 6.2 Saran ........................................................................................23
DAFTAR TABEL 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14.
hal Akibat yang terjadi dari proses urbanisasi terhadap respon hidrologis .................3 Jenis penggunaan lahan di DAS Ciliwung hulu ..................................................4 Morfometri DAS Ciliwung hulu ..........................................................................10 Skenario perubahan penggunaan lahan .............................................................13 Luas penggunaan/penutupan lahan DAS Ciliwung hulu ...................................14 Curah hujan bulanan rata-rata DAS Ciliwung hulu (1993-2005) ......................15 Perbandingan luaran model pada beberapa nilai parameter penggunaan lahan dengan luaran model dan nilai awal ...................................20 Perbandingan luaran model pada beberapa nilai peubah curah hujan dan nilai awal .....................................................................................................21 Debit hasil luaran model pada berbagai skenario lahan dan hujan ......................22 Analisis sifat curah hujan bulanan wilayah DAS Ciliwung hulu ......................28 Contoh perhitungan indikator fungsi hidrologi ................................................33 Perubahan dominan tipe dan luas penggunaan lahan di sub DAS Ciliwung hulu tahun 1994 – 2004 .............................................................34 Data debit harian bendung katulampa 1994 ................................................35 ................................................36 Data debit harian bendung katulampa 2004
DAFTAR GAMBAR hal 1. Daur hidrologi DAS .........................................................................................2 2. Lima faktor yang mempengaruhi partisi air hujan menjadi komponen debit sungai dan evapotranspiras............................................................................3 3. Hasil simulasi debit Sumberjaya Lampung ………………………………..7 4. Hasil simulasi debit di Mae Chaem Thailand …………………………….….7 5. Diagram konsep model ……………………………………………….............…9 6. Diagram alur proses hidrologi ............................................................................9 7. Diagram alur model dalam Stella .............................................................11 8. Diagram sub model neraca air ..........................................................................11 9. Diagram sub model aliran sungai .............................................................11 10. Peta lokasi kajian .......................................................................................14 11. Grafik curah hujan bulanan 3 stasiun di wilayah Ciliwung hulu ......................15 12. Grafik curah hujan wilayah bulanan rata – rata ................................................16 13. Grafik nilai evapotranspirasi ..........................................................................16 14. Grafik time series hubungan curah hujan dan debit harian ......................17 15. Grafik hubungan antara curah hujan dengan distribusi debit bersarkan probabilitas............................................................................................................17 16a.Hubungan indikator penyangga terhadap TWY ................................................17 16b.Hubungan indikator penyangga puncak kejadian hujan terhadap TWY..............18 16c.Hubungan indikator penyangga relative terhadap total debit terhadap TWY......18 17. Grafik nilai KRS tahun 1993 – 2005 .............................................................18 18a.Hasil simulasi GenRiver tahun 1994 ……………………………………….19 18b Hasil simulasi GenRiver tahun 2004 ……………………………………….19 19 Plot hubungan debit model skenario terhadap debit model rata-rata .........22
DAFTAR LAMPIRAN hal Lampiran 1 .....................................................................................................26 Lampiran 2 .....................................................................................................27 Lampiran 3 .....................................................................................................28 Lampiran 4 .....................................................................................................29 Lampiran 5 .....................................................................................................30 Lampiran 6 .....................................................................................................31 Lampiran 7 .....................................................................................................32 Lampiran 8 .....................................................................................................33 Lampiran 9 .....................................................................................................34 Lampiran10.....................................................................................................35 Lampiran11.....................................................................................................36
PENDAHULUAN 1.1
Latar Belakang
Daerah Aliran Sungai (DAS) Ciliwung merupakan salah satu DAS yang penting untuk mendapat perhatian, karena DAS ini mempunyai topografi yang cukup beragam dengan daerah hulu yang berbukit – bukit dan daerah hilir yang relatif datar, sehingga pengamanan DAS ini perlu mendapat perhatian. Hal ini juga ditunjang dengan Perda No. 2 Tahun 2003 tentang RTRW Propinsi Jawa Barat 2010, telah ditetapkan kebijakan untuk “meningkatkan fungsi dan kualitas kawasan lindung di Jawa Barat,” termasuk kawasan lindung di Kawasan Bodebek dan Bopunjur. Berdasarkan kebijakan tersebut di atas, upaya pengelolaan sumber daya air di DAS Ciliwung Cisadane untuk mengatasi krisis air Jakarta adalah melalui penataan situ, waduk, dan sungai sebagai sarana dan prasarana konservasi, penyedia air baku, dan pengendali banjir, serta melalui konservasi lahan yang diprioritaskan pada kawasan lindung baik di dalam kawasan hutan maupun di luar kawasan hutan melalui rehabilitasi lahan kritis, pengendalian pemanfaatan lahan dan pengendalian kualitas air. Pengelolaan DAS seringkali dihubungkan dengan tingkat penutupan lahan oleh hutan, dengan asumsi bahwa ‘reforestasi’ atau ‘reboisasi’ dapat mengembalikan dampak negative dari terjadinya deforestasi (penggundulan hutan). Pada saat ini telah tersedia model simulasi yang dapat dipakai untuk mempelajari dinamika pori makro tanah yang berhubungan dengan sifat hujan menurut skala waktu dan ruang. Model tersebut disusun berdasarkan hasil pengukuran yang intensif dari berbagai sub DAS dan dapat digunakan untuk memprediksi pengaruh alih guna lahan terhadap indikator fungsi hidrologi DAS. Dengan demikian, model tersebut dapat digunakan untuk ekstrapolasi berbagai skenario sistem penggunaan lahan di masa yang akan datang. Untuk mengkaji fungsi hidrologi suatu DAS dari waktu ke waktu dibutuhkan data histori yang relatif panjang agar kesimpulan yang diperoleh dapat mewakili. DAS Ciliwung merupakan salah satu DAS dengan pengelolaan yang terpadu dan intensif sehingga data – data daerah ini
cukup tersedia. Alasan tersebut yang mendasari pengkajian fungsi hidrologi di DAS Ciliwung hulu. 1.2
Tujuan Penelitian
Tujuan dari penelitian ini adalah untuk mengetahui pengaruh alih guna lahan terhadap fungsi hidrologi tahun 1993 sampai dengan 2005 melalui pendekatan kuantitatif dan memprediksi debit apabila terjadi alih guna lahan.
II. 2.1
TINJAUAN PUSTAKA Definisi DAS
Daerah Aliran Sungai (DAS) menurut Manan (1978) didefinisikan sebagai suatu kawasan yang dibatasi oleh pemisah topografis yang menampung, menyimpan, dan mengalirkan curah hujan yang jatuh di atasnya, ke sungai utama yang bermuara ke danau atau lautan. Pemisah topografi adalah punggung bukit. Dari definisi di atas, dapat dikemukakan bahwa DAS merupakan ekosistem, dimana unsur organisme dan lingkungan biofisik serta unsur kimia berinteraksi secara dinamis dan di dalamnya terdapat keseimbangan inflow dan outflow dari material dan energi. Selain itu pengelolaan DAS dapat disebutkan merupakan suatu bentuk pengembangan wilayah yang menempatkan DAS sebagai suatu unit pengelolaan sumber daya alam (SDA) yang secara umum untuk mencapai tujuan peningkatan produksi pertanian dan kehutanan yang optimum dan berkelanjutan (lestari) dengan upaya menekan kerusakan seminimum mungkin agar distribusi aliran air sungai yang berasal dari DAS dapat merata sepanjang tahun. Dalam pendefinisian DAS pemahaman akan konsep daur hidrologi sangat diperlukan terutama untuk melihat masukan berupa curah hujan yang selanjutnya didistribusikan melalui beberapa cara seperti diperlihatkan pada Gambar 1. Konsep daur hidrologi DAS menjelaskan bahwa air hujan langsung sampai ke permukaan tanah untuk kemudian terbagi menjadi air limpasan, evaporasi dan air infiltrasi, yang kemudian akan mengalir ke sungai sebagai debit aliran. Selama berlangsungnya siklus hidrologi, yakni sepanjang air dari permukaan laut ke atmosfer kemudian ke permukaan tanah dan kembali lagi ke laut yang terus menerus bersirkulasi, penguapan, presipitasi dan pengaliran ke luar.
Air menguap ke udara dari permukaan tanah dan laut, berubah menjadi awan sesudah melalui beberapa proses dan kemudian jatuh sebagai hujan atau salju ke permukaan laut, sungai, atau daratan. Sebelum tiba ke permukaan bumi, tidak semua bagian hujan yang jatuh ke permukaan bumi mencapai permukaan tanah. Sebagian akan tertahan oleh tumbuh – tumbuhan dimana sebagian akan menguap (intersepsi) dan sebagian lagi akan jatuh atau mengalir melalui dahan – dahan menuju ke permukaan tanah (troughfall dan steamflow). Sebagian air hujan yang tiba ke permukaan tanah akan masuk ke dalam tanah (infiltrasi). Bagian lain yang merupakan kelebihan akan mengisi lekukan – lekukan permukaan tanah, kemudian mengalir ke daerah – daerah rendah, masuk ke sungai, dan akhirnya ke laut.dalam perjalanannya ke laut sebagian akan menguap kembali ke udara. Sebagian air yang masuk ke dalam tanah sebelum menjadi air bawah tanah keluar kembali segera ke sungai sebagai aliran bawah permukaan (interflow), tetapi sebagian besar akan tersimpan sebagai air bawah tanah (groundwater) yang akan keluar sedikit demi sedikit dalam jangka waktu yang lama ke sungai sebagai aliran bawah tanah (groundwater flow).
Gambar 1. Daur Hidrologi DAS (Pawitan, 2004) 2.2
Definisi Penggunaan Lahan dan Tutupan Lahan
Penggunaan lahan (land use) adalah setiap bentuk campur tangan (intervensi) manusia terhadap lahan dalam rangka memenuhi kebutuhan hidupnya baik material maupun spiritual (Vink, 1975). Penggunaan lahan menurut Sitorus (2001) secara luas dapat dikelompokkan kedalam 2 kategori, yaitu: 1.
Penggunaan lahan pedesaan (rural land use) yang secara luas meliputi kegiatan pertanian, kehutanan, konservasi satwa liar serta
pengembangan tempat rekreasi. 2.
dan
pengelolaan
Penggunaan lahan perkotaan dan industri (urban and industrial land use) termasuk kota, komplek industri, jalan raya, serta fasilitas pertambangan.
Adapun pembagian lainnya dikelompokkan kedalam dua kelompok besar, yaitu pertanian dan non pertanian. Penggunaan lahan pertanian meliputi berbagai macam penggunaan lahan berdasarkan atas penyediaan air dan komoditas yang diusahakan. Sedangkan penggunaan lahan non pertanian meliputi pemukiman desa dan kota, industri, rekreasi, dan pertambangan. Dalam Undang-Undang Republik Indonesia nomor 24 tahun 1992 tentang Penataan Ruang tertulis: pemanfaatan ruang meliputi kawasan perdesaan, kawasan perkotaan, kawasan lindung serta kawasan budidaya. Kawasan lindung adalah kawasan yang ditetapkan dengan fungsi utama melindungi kelestarian lingkungan hidup yang mencakup sumberdaya alam dan sumberdaya buatan. Sedangkan kawasan budidaya merupakan kawasan yang ditetapkan dengan fungsi utama untuk dibudidayakan atas dasar kondisi dan potensi sumberdaya alam, sumberdaya manusia, dan sumberdaya buatan. Pemanfaatan kawasan budidaya untuk memenuhi kebutuhan hidup manusia diusahakan tanpa mengganggu dan merusak ekosistem kawasan lindung. Namun sebaliknya, seiring dengan meningkatnya kebutuhan manusia akan pangan dan perumahan seringkali terjadi konversi lahan. Banyak kawasan lindung beralih fungsi menjadi kawasan budidaya, akibatnya fungsi kawasan lindung menjadi terganggu. Sebagaimana yang terjadi pada tragedi “Banjir bandang di Bukit Lawang Bohorok”, Sumatera Utara, yang dipicu oleh kerusakan ekosistem hutan lindung. Diperkirakan tingkat kerusakan Taman Nasional Gunung Leuser (TNGL) saat ini adalah 170.000 Ha (22 %) dari total luas TNGL yang 788.000 Ha akibat penebangan liar dan perambahan hutan (Walhi, 2003). Tutupan lahan oleh pohon (tutupan pohon) dengan segala bentuknya dapat mempengaruhi aliran air. Tutupan pohon tersebut dapat berupa hutan alami, atau sebagai permudaan alam (natural regeneration), pohon yang dibudidayakan,
pohon sebagai tanaman pagar, atau pohon monokultur (misalnya hutan)
Gambar 2. Lima faktor yang mempengaruhi partisi air hujan menjadi komponen debit sungai dan evapotranspirasi (Van Noordwijk et al., 2004).
permukaan dan berbagai lapisan dalam profil tanah. Struktur tanah juga dipengaruhi oleh aktivitas biota yang sumber energinya tergantung kepada bahan organik (seresah di permukaan, eksudasi organik oleh akar, dan akar-akar yang mati). Ketersediaan makanan bagi biota (terutama cacing tanah), penting untuk mengantisipasi adanya proses peluruhan dan penyumbatan pori makro tanah. •
Faktor – faktor yang mempengaruhi jumlah serapan air oleh pohon adalah fenologi pohon, distribusi akar dan respon fisiologi pohon terhadap cekaman parsial air tersedia. Serapan air oleh pohon diantara kejadian hujan akan mempengaruhi jumlah air yang dapat disimpan dari kejadian hujan berikutnya, sehingga selanjutnya akan mempengaruhi proses infiltrasi dan aliran permukaan. Serapan air pada musim kemarau, khususnya dari lapisan tanah bawah akan mempengaruhi jumlah air tersedia untuk ‘aliran lambat’ (slow flow).
Pengaruh tutupan pohon terhadap aliran air berdasarkan Van Noordwijk et al., (2004) adalah dalam bentuk berikut: •
Intersepsi air hujan Banyaknya air yang diintersepsi dan dievaporasi tergantung pada indeks luas daun (LAI), karakteristik permukaan daun, dan karakteristik hujan. Intersepsi merupakan komponen penting jika jumlah curah hujan rendah, tetapi dapat diabaikan jika curah hujan tinggi. Apabila curah hujan tinggi, peran intersepsi pohon penting dalam kaitannya dengan pengurangan banjir.
•
Daya pukul air hujan Vegetasi melindungi permukaan tanah dari pukulan langsung tetesan air hujan yang dapat menghancurkan agregat tanah yang akan menyumbat pori tanah makro sehingga menghambat infiltrasi air tanah, akibatnya limpasan permukaan akan meningkat.
•
Infiltrasi air
Serapan air
•
Drainase lansekap Besarnya drainase suatu lansekap (bentang lahan) dipengaruhi oleh beberapa faktor antara lain kekasaran permukaan tanah, relief permukaan tanah yang memungkinkan air tinggal di permukaan tanah lebih lama sehingga mendorong terjadinya infiltrasi, tipe saluran yang terbentuk akibat aliran permukaan yang dapat memicu terjadinya ‘aliran cepat air tanah’ (quick flow). Menurut Kibler (1982) dalam Rogers (1998) akibat yang terjadi karena proses urbanisasi terhadap hidrologi adalah seperti tertera pada Tabel 1.
Proses infiltrasi tergantung pada struktur tanah pada lapisan Tabel 1. Akibat yang terjadi dari proses urbanisasi terhadap respon hidrologis Pengaruh urbanisasi Respon hidrologis Perubahan vegetasi
Penurunan evapotranspirasi dan intersepsi dan peningkatan sedimentasi sungai
Konstruksi dasar perumahan dan inftrastruktur
Pengembangan pembangunan dan perdagangan Konstruksi drainase dan pengembangan saluran air
Penurunan infiltrasi dan penurunan muka air tanah; pernurunan aliran debit dan penurunan aliran dasar selama musim kemarau Meningkatkan volume aliran permukaan Memicu peluapan banjir di aliran sungai
Sumber: Kibler (1992) dalam Rogers (1998) 2.3
Jenis Penggunaan Lahan
Harimurti (1999) memberikan definisi dan batasan yang jelas mengenai tipe – tipe penggunaan lahan yang ada di No 1 2
3 4 5
6
7 8
2.4
kawasan DAS Ciliwung hulu seperti yang dijelaskan pada Tabel 2 berikut.
Tabel 2. Jenis Penggunaan Lahan di DAS Ciliwung Hulu Nama Lahan Definisi Hutan Lebat wilayah yang ditutupi oleh vegetasi pepohonan, baik alami maupun yang dikelola, dengan tajuk yang rimbun dan besar/lebat. Hutan Semak/belukar hutan yang telah dirambah/dibuka, merupakan area transisi dari hutan lebat menjadi kebun atau lahan pertanian, bisa berupa hutan dengan semak/belukar dengan tajuk yang relatif kurang rimbun. Kebun Campuran daerah yang ditumbuhi vegetasi tahunan satu jenis maupun campuran baik dengan pola acak, maupun teratur sebagai pembatas tegalan Pemukiman kombinasi antara jalan, bangunan, pekarangan, dan bangunan itu sendiri. Sawah daerah pertanian yang ditanami padi sebagai tanaman utama dengan rotasi tertentu yang biasanya diairi sejak saat penanaman hingga beberapa hari sebelum panen. Tegalan daerah yang umumnya ditanami tanaman semusim, namun pada sebagian lahan tidak ditanami, dengan vegetasi yang umum dijumpai seperti padi gogo, singkong, jagung, kentang, kedelai, dan kacang tanah. Lahan Terbuka daerah yang tidak ditemukan vegetasi berkayu, umumnya hanya jenis rerumputan maupun penggunaan lain akibat aktivitas manusia. Kebun Teh daerah yang digunakan sebagai perkebunan the baik yang diusahakan pemerintah maupun swasta.
Alih guna lahan
Alih guna lahan adalah bertambahnya suatu penggunaan lahan dari satu sisi penggunaan ke penggunaan yang lainnya diikuti dengan berkurangnya tipe penggunaan lahan yang lain dari suatu waktu ke waktu berikutnya, atau berubahnya fungsi suatu lahan pada kurun waktu yang berbeda (Wahyunto et al., 2001). Penelitian Janudianto (2004) menjelaskan alih guna lahan di Sub DAS Ciliwung Hulu didominasi oleh kecenderungan perubahan lahan pertanian (sawah) menjadi lahan pemukiman dan
perubahan hutan menjadi lahan perkebunan (kebun teh). 2.5
Kriteria dan Indikator Fungsi Hidrologi
Kriteria dan indikator kuantitatif diperlukan dalam mempelajari fungsi hidrologi DAS. Kriteria dan indikator yang ditetapkan berdasarkan pemahaman kuantitatif hujan yang terbagi menjadi evapotranspirasi, aliran sungai dan perubahan penutupan serta pola penggunaan lahan sesuai dengan karakteristik lokal. Fluktuasi debit sungai dan curah hujan dijadikan parameter utama untuk menilai indikator
penyangga (buffering indicator) akibat alih guna lahan. Indikator penyangga (buffering indicator) cenderung berkorelasi negatif dengan total debit sungai sehingga peningkatan debit akan menurunkan kapasitas menyangga dari sungai (Farida et al., 2004). Indikator penyangga menunjukkan tingkat penurunan yang relatif rendah pada kondisi puncak kejadian hujan (buffering peak events). Peningkatan total debit tidak selalu diikuti dengan peningkatan debit terendah (bulanan) akibat adanya variabilitas hujan antar tahun (inter-annual). 2.6
Kejadian Banjir di Jakarta
Berdasarkan dokumentasi yang tersedia (Bappenas, 2007), Kota Jakarta dilanda banjir pada tahun 1621, 1654 dan 1918, kemudian pada tahun 1976, 1996, 2002 dan 2007. Banjir di Jakarta yang terjadi pada tahun 1996 selain menggenangi hampir seluruh penjuru kota juga menjadi tragedi nasional yang menjadi perhatian dunia. Banjir besar ini dipercaya sebagai banjir lima tahunan yang akan berulang setiap lima tahun. Pada awal 2002 banjir melanda Jakarta dan sekitarnya dan pada awal 2007 banjir kembali melanda wilayah Jakarta dan sekitarnya dengan cakupan wilayah genangan yang lebih luas. Namun menurut BMG (dikutip dari Tempointeraktif.com, Maret 2007), banjir 2007 bukanlah yang terbesar karena bukan hanya dihitung dari curah hujan sesaat, namun dari jumlah curah hujan 29 hari terakhir Data curah hujan kumulatif 5, 10, 20 dan 29 hari terakhir di DKI Jakarta pada tahun – tahun banjir menunjukkan Februari 2006 berturut – turut 221,4mm, 285, 7mm, 341,7mm, 442,1mm. Januari 2002 berturut – turut 232,9mm, 361,7mm, 572,7mm, 668mm. Bulan Februari 2007 berturut – turut 327mm, 401,5mm, 427mm, 427mm. Sebelum banjir besar 2007, banjir terbesar yang pernah melanda DKI Jakarta adalah kejadian banjir Januari 2002 dan awal Februari 2002, dan merupakan banjir terburuk yang pernah melanda DKI Jakarta pada kurun waktu ke belakang (www.tumoutou.net/pps702_9145/nana_mul yana.pdf/diakses April, 2006). Berdasarkan pengamatan data curah hujan yang terekam di beberapa stasiun pada saat terjadi banjir tahun 2002 terlihat bahwa curah hujan harian yang turun di Halim PK dan Ciledug merupakan periode ulang 2 tahunan, di daerah Depok,
Citeko, Tanjung Priok, dan Darmaga periode ulang 5 tahun, dan di kantor pusat BMG merupakan priode ulang 10 tahun, dengan demikian waktu frekuensi terjadinya banjir seperti tahun 2002 kisaran waktunya 2 – 5 tahun, sehingga perlu mendapat perhatian yang sangat serius. Akan tetapi apabila dilihat dari curah hujan kumulatif 2 hari (Duration Dept frekuensi /DDF) curah hujan yang tergolong tinggi dengan kisaran 130 – 295 mm. Akibat curah hujan yang turun selama awal Januari, menyebabkan kondisi tanah jenuh air, sehingga sangat sedikit air yang diinfiltrasikan. Pada tanggal 30 Januari terjadi pengaruh pasang air laut yang tertinggi di pantai utara Jakarta, sehingga curah hujan yang tinggi di bagian hulu DAS Ciliwung bersamaan dengan terjadinya pasang tertinggi sehingga banjir pada akhir Januari merupakan gabungan antara kondisi pasang surut dan pola curah hujan tinggi yang terjadi di bagian hulu dan hilir DAS Ciliwung serta Jakarta dan sekitarnya. Berdasarkan pola induk yang telah dibuat tahun 1973 dan kemudian disempurnakan tahun 1997 setelah ada banjir besar yang melanda tahun 1996, nampak telah terjadi kenaikan debit rencana pada semua badan sungai yang ada di DKI Jakarta. Master plan Cengkareng Drain telah dinaikkan dari 370 m3/det menjadi 620 m3/det, sementara sungai Ciliwung telah dinaikkan dari 370 m3/det menjadi 570 m3/det. Perubahan pola induk ini untuk mengantisipasi kenaikan debit sungai akibat perubahan tata guna lahan, khususnya kurangnya daerah resapan dan terlalu dominannya pemukiman akibatnya pesatnya pertumbuhan di kawasan Jabodetabek. 2.7
Deskripsi singkat model GenRiver
Model Aliran Sungai Generik merupakan model yang (GenRiver) dikembangkan berdasarkan proses hidrologi (process based model). Simulasi model GenRiver menggunakan Stella sebagai software yang dihubungkan dengan file microsoft excel. Input utama dari model ini adalah curah hujan, tingkat penutupan lahan dan sifat fisik tanah dengan keluaran utama berupa aliran sungai dan neraca air untuk skala DAS Bagian utama dari GenRiver meliputi neraca air pada skala plot (patch level water balance) berdasarkan curah hujan dan modifikasi sifat fisik tanah dan
penutupan lahan. Plot – plot ini memiliki kontribusi terhadap aliran sungai melalui aliran permukaan pada saat terjadinya hujan (surface quick flow), aliran air dalam tanah yang terjadi setelah hujan (soil quick flow) dan aliran dasar (base flow) yang berasal dari pelepasan air tanah secara bertahap menuju sungai (gradual release of groundwater).
menghasilkan aliran permukaan yang dibatasi oleh infiltrasi (infiltration limited runoff’), sedangkan pada kondisi kedua aliran permukaan yang terjadi merupakan aliran jenuh permukaan (saturation overland flow).
Curah hujan harian. Curah hujan untuk skala sub-DAS dapat diambil dari data empiris atau menggunakan data bangkitan dari pembangkit data acak (random generator) yang mempertimbangkan pola temporal (seperti model rantai Markov) atau model yang mempertimbangkan korelasi ruang (spatial correlation) dari hujan pada waktu tertentu.
Evapotranspirasi Total evapotranspirasi yang digunakan model ini mengikuti pada evapotranspirasi potensial Penman – Monteith dengan faktor koreksi yang dipengaruhi oleh: (a) air yang terintersepsi oleh kanopi, (b) kondisi tutupan lahan yang terkait dengan sensitivitas setiap jenis penutupan lahan terhadap kekeringan, (c) faktor pembobot pada evapotranspirasi potensial harian yang mengikuti fenologi dan pola tanam, (d) relatif potensial evapotranspirasi (bulanan) untuk setiap tipe penutupan lahan.
Intensitas hujan dan waktu untuk infiltrasi. Intensitas hujan dihitung dari rata - rata data empiris intensitas hujan (mm/jam) dengan mempertimbangkan koefisien variasi dari kumpulan data tersebut. Lamanya hujan menentukan waktu yang tersedia untuk proses infiltrasi. Namun parameter ini dapat dimodifikasi dengan mempertimbangkan intersepsi oleh kanopi dan lamanya penetesan air dari kanopi (dripping phase) dengan penetapan awal (default) 30 menit.
Redistribusi air tanah. Selama kejadian hujan, tanah dapat mencapai kondisijenuh air, namun sehari setelah hujan kondisi akan kembali pada kapasitas lapang (kondisi air tanah setelah 24 jam dari kejadian hujan lebat). Perbedaan antara kondisi jenuh dan kapasitas lapang dipengaruhi oleh: (a) transpirasi, (b) adanya aliran air ke zona bawah, (c) adanya aliran air ke sungai sebagai aliran cepat air tanah (soil quick flow) apabila air yang ada melebihi kapasitas lapang
Pelepasan air tanah menuju sungai (melalui aliran dasar) Jarak (routing distance). Jarak titik pengamatan atau outlet DAS ditentukan dari titik pusat masing – masing subDAS. Waktu tempuh (routing time) dari masing–masing sub-DAS dapat diturunkan dari data jarak dan asumsi rata-rata kecepatan aliran air.
Komponen utama model GenRiver dan proses-proses yang terlibat sebagai berikut:
Intersepsi. Kapasitas penyimpanan air terintersepsi merupakan fungsi linier dari luas area daun dan ranting dari berbagai tipe penutupan lahan. Evaporasi dari air yang terintersepsi (interceptionevaporation) mempunyai prioritas sesuai dengan kebutuhan transpirasi tanaman.
Infiltrasi dan aliran permukaan. Proses infiltrasi dihitung berdasarkan nilai minimum dari : (a) kapasitas infiltrasi harian dan waktu yang tersedia untuk infiltrasi (ditentukan oleh intensitas hujan dan kapasitas penyimpanan lapisan permukaan tanah), (b) jumlah air yang dapat disimpan oleh tanah pada kondisi jenuh dan jumlah air yang dapat memasuki zona air tanah pada rentang waktu satu hari. Apabila kondisi pertama yang terjadi maka model akan
Keluaran dari model ini berupa debit sungai harian dan kumulatif neraca air tahunan. Pengolahan lebih lanjut dari output model ini dapat digunakan sebagai indikator dalam mempelajari fungsi DAS. 2.8
Aplikasi Model
Model GenRiver telah digunakan untuk dua lokasi yaitu DAS Sumbarjaya Lampung dan Mae Chaem Thailand Utara. Pada kedua DAS ini telah disimulasikan untuk debit prediksi yang diperoleh dari input curah hujan dan penutupan lahan. Namun
tidak dicantumkan nilai koefisien determinasinya hanya dari pola grafik yang teratur. Di DAS Sumberjaya dengan total luas 404 km2 menggunakan tiga skenario, yaitu seluruh DAS merupakan lahan hutan, seluruh DAS berupa lahan alang – alang, dan kondisi penutupan lahan sebenarnya. Sedangkan untuk DAS Mae Chaem dengan luas 40.000 km2 menggunakan tiga skenario, yaitu kondisi penutupan lahan sebenarnya, seluruhnya berupa hutan hijau sepanjang tahun, dan hutan yang berganti daun sepanjang tahun. Berikut beberapa gambar hasil simulasi GenRiver di DAS Sumberjaya Lampung dan Mae Chaem Thailand (Van Noordwijk et al., 2004).
2.
3. 4. 5.
Dan data penggunaan lahan tahun 2004 dari citra Ikonos yang dipetakan oleh Prasatya (2006) Data debit harian tahun 1993 2005 Bendung Katulampa, Ciliwung Hulu dari Balai Pengelolaan Sumber Daya Air Ciliwung – Cisadane Bogor. Data curah hujan harian dari stasiun Katulampa, Citeko, dan Gunung Mas tahun 1993 - 2005. Data jenis tanah Ciliwung Hulu Data morfometri sungai (subDAS, panjang sungai utama).
Alat yang digunakan meliputi: Seperangkat komputer dengan software Arcview 3.3, dan MS Office 2003(Excel dan Word), Minitab 14, Stella Vr. 511. 3.3
Metodologi Penelitian ini meliputi kegiatan:
3.3.1
Gambar 3. Hasil simulasi debit Sumberjaya Lampung
Tahap persiapan diawali dengan studi literatur dan pengumpulan data-data berupa Peta Penggunaan Lahan Sub DAS Ciliwung Hulu tahun 1994 dan 2004, data debit, data tanah, dan data penunjang lainnya. Langkah berikutnya adalah menyeragamkan atau kalibrasi data, khususnya untuk data penggunaan lahan sehingga memungkinkan proses analisis spasial. 3.3.2
Gambar 4. Hasil simulasi debit di Mae Chaem Thailand
III. 3.1
METODOLOGI
Waktu dan Tempat
Penelitian dilaksanakan mulai bulan September sampai dengan Desember 2006 di Laboratorium Klimatologi Departemen Geofisika dan Meteorologi. 3.2
Bahan dan Alat Bahan yang diguanakan penelitian ini meliputi: 1.
dalam
Peta penggunaan lahan DAS Ciliwung Hulu tahun 1994 dan 2004, skala 1:150.000 hasil penelitian Janudianto (2004).
Studi Literatur
Analisis Alih guna lahan
Operasi tumpang susun (overlay) dilakukan menggunakan data digital peta penggunaan/penutupan lahan dengan bantuan ArcView 3.3. Operasi tumpang susun dilakukan antara peta penggunaan/penutupan lahan tahun 1994 dan 2004 yang bertujuan untuk melihat arah dan pola perubahan penggunaan/penutupan lahan. Ekstraksi data atribut hasil dari operasi tumpang susun ini digunakan sebagai data dalam teknik analisis selanjutnya. 3.3.3
Menentukan kriteria indikator fungsi hidrologi
dan
Kriteria dan indikator kuantitatif diperlukan dalam mempelajari fungsi hidrologi DAS. Fluktuasi debit sungai dan curah hujan dijadikan parameter utama untuk menilai indikator penyangga (buffering indicator) akibat alih guna lahan. Kriteria ini dapat dihubungkan langsung dengan
pengertian kuantitatif bagaimana hujan atau presipitasi (P) terurai menjadi aliran sungai (Q) dan evapotranspirasi (E) pada suatu sistem neraca air. Hubungan antara faktorfaktor tersebut dapat membantu kita dalam memahami logika dan tarik ulur antara perubahan transmisi air, daya sangga kejadian puncak hujan dan fungsi DAS dalam menyalurkan air secara perlahan. Melalui pemanfaatan data empiris curah hujan dan aliran sungai dan/ atau hasil simulasi model, maka dapat dikembangkan beberapa indikator kuantitatif untuk tiga kriteria utama fungsi DAS. 3.3.3.1 Transmisi air Merupakan hasil air per unit curah hujan tiap suatu luasan DAS. Sebagai indikator adalah Total debit sungai per unit hujan (TWY) yang dirumuskan sebagai berikut:
TWY = ∑ Q /( A * ∑ P)
atau curah hujan di atas rata – rata selama bulan – bulan basah atau musim hujan. BPE= 1− Max(daily_ Q − Qmean) /(A* Max(daily_ P − Pmean))
Satuan debit yang digunakan untuk analisis indikator adalah dalam mm/hari. 3.3.3.3 Koefisien Rejim Sungai (KRS) Besarnya fluktuasi debit aliran sungai yang terjadi dapat diketahui dengan menggunakan koefisien rejim sungai. Koefisien rejim sungai merupakan perbandingan antara debit harian rata – rata maksimum dan debit harian rata – rata minimum. KRS bisa digunakan untuk mengevaluasi kondisi suatu DAS. Semakin kecil koefisien ini berarti kondisi hidrologi dari suatu wilayah DAS semakin baik Batasan yang diberikan untuk menilai indikator ini adalah sebagai berikut (Asdak, 1995):
Dengan: Q = debit sungai (mm/hari) P = presipitasi (mm/hari) A = luasan DAS (km2) Satuan unit parameter tersebut selanjutnya digunakan untuk perhitungan indikator penyangga. 3.3.3.2 Penyangga pada puncak kejadian hujan. Merupakan indikator kejadian banjir relatif terhadap kejadian hujan pada suatu luasan DAS. Sebagai indikator dari kriteria ini adalah: indicator (BI)/Indikator a. Buffering penyangga Merupakan koefisien yang menyatakan kemampuan suatu DAS untuk menyangga hasil air (debit) tiap satuan waktu curah hujan pada suatu luasan DAS.
BI = ( PabAvg − (QabAvg / A)) / PabAvg = 1 − QabAvg /( A * PabAvg )
b.
dimana : PabAvg = ∑ max(P-Pmean,0) QabAvg = ∑ max(Q-Qmean,0) buffering indicator Relative (RBI)/Indikator penyangga relatif terhadap total debit
RBI = 1 − (QabAvg / Qtot ) /(PabAvg / Ptot ) c.
Buffering peak event (BPE)/Indikator penyangga puncak kejadian hujan Koefsien yang menyatakan kemampuan suatu DAS dalam menyangga total hasil air (debit) pada saat hujan maksimum
• KRS < 50 ; baik • 50 ≤ KRS < 120 ; sedang • KRS ≥ 120 ; buruk 3.3.4
Simulasi model GenRiver
Simulasi model GenRiver menggunakan Stella sebagai software yang dihubungkan dengan file microsoft excel. Input utama dari model ini adalah curah hujan, tingkat penutupan lahan dan sifat fisik tanah dengan keluaran utama berupa aliran sungai dan neraca air untuk skala DAS (Gambar 6). Komponen utama model GenRiver dan proses-proses yang terlibat sebagai berikut : curah hujan harian, intensitas hujan dan waktu untuk infiltrasi, intersepsi, infiltrasi dan aliran permukaan, evapotranspirasi, redistribusi air tanah, pelepasan air tanah menuju sungai (melalui aliran dasar), jarak (routing distance).
Input: CH, tutupan lahan,tanah,
Model
Keluaran
Data pengukuran
model Kesalahan minimum
tidak
kalibrasi
ya Aplikasi model
Gambar 5. Diagram konsep model (Sutoyo, 1999) CURAH HUJAN TOTAL evapotranspirasi
Curah hujan langsung
Intersepsi Air
Air hilang terintersepsi
Air Lolos
Aliran batang
Curah Hujan Bersih
Air larian
Permukaan Tanah
Evaporasi tanah
Kelembaban Tanah
Transpirasi vegetasi
Debit
Gambar 6. Diagram alur proses hidrologi (Asdak, 1995)
Model GenRiver dihubungkan dengan sebuah file Excel GenRiver.xls. File ini berisi parameter masukan yang terhubung pada GenRiver.STM. Parameter utama masukan untuk model GenRiver pada file excelnya adalah sebagai berikut: • Data hujan dan debit Data hujan dalam satuan mm/hari sedangkan untuk data debit dalam m3/detik. Data hujan yang digunakan adalah data dari tiga stasiun yang dianggap mewakili wilayah DAS Ciliwung hulu yaitu stasiun Katulampa, stasiun Citeko, dan stasiun Gunung Mas dari tahun 1994 sampai dengan 2004. Sedangkan data debit diambil dari Stasiun Pengamatan Aliran Sungai (SPAS) Katulampa pada periode 1993 – 2005. Curah hujan wilayah dihitung dengan metode Thiessen, yaitu:
•
ET0 = w * Rn + (1 − w) f (u)(ea − ed ) dengan w faktor pembobot merupakan fungsi suhu dan ketinggian. ea tekanan uap jenuh; ed tekanan uap aktual ⎡ 17 . 27 T dew ⎤ e a = e 0 ( T dew ) = 0 . 6108 exp ⎢ ⎥ ⎣ T dew + 237 . 3 ⎦
n
Pg = ∑ Wi ⋅ Pi dengn i =1
Wi = Ai / ∑ Ai
•
Dimana: Pg : curah hujan wilayah (mm) Pi : curah hujan pada stasiun ke (mm) Ai : luas polygon ke-i (Ha) Wi : bobot stasiun ke-i i :jumlah stasiun pengamatan • Data tutupan lahan Terdiri atas jenis tutupan lahan, parameter – parameter fisik tanah, dan evapotranspirasi. Data tutupan lahan yang dianalisis yaitu data tahun 1994 dan 2004 (di lampiran) hasil penelitian Janudianto (2004) dan Prasatya (2006) . Untuk paremeter – parameter fisik tanah mengacu pada database yang ada pada
Sub DAS Ciliwung Hulu
file GenRiver.xls yang disesuaikan dengan kondisi wilayah Ciliwung hulu. Nilai evapotranspirasi bulanan diperoleh dengan software ETo (Irsis) versi 1.0 keluaran Leuven University. Penentuan evapotranspirasi menggunakan metode Penman FAO yang dirumuskan sebagai berikut:
Input data yang digunakan untuk analisis evapotranspirasi model Penman FAO adalah: -Temperatur : Trata – rata (oC), T max (oC), Tmin (oC) - Kelembaban: RH rata – rata (%) - Radiasi surya - Kecepatan angin: Urata – rata (km/jam) Informasi sub DAS DAS dibagi menjadi beberapa sub-DAS (dalam km2) dan jarak suatu titik terhadap pelepasannya (routing distance) tiap – tiap sub-DAS. Data morfometri DAS Ciliwung hulu diperoleh dari interpretasi peta Landsat tahun 2000 hasil penelitian Irianto (2000) yang didasarkan atas nilai bilangan kurva. Nilai bilangan kurva sendiri telah ditetapkan berdasarkan klasifikasi kelompok hidrologi tanah, klasifikasi komplek penutup tanah dan kondisi kandungan air dalam tanah sebelumnya. Data tersebut selanjutnya dianggap berlaku pada saat penelitian.
Tabel 3. Morfometri DAS Ciliwung Hulu Panjang No. Sub Sub DAS Luas sungai utama ha m 1 Ciliwung Hulu 1 4908 14750 2 Cisarua 2237 16500 3 Ciliwung Hulu 2 240 2300 4 Cibogo 1376 14000 5 Ciliwung Hulu 3 106 1050 6 Cisukabirus 1688 15300 7 Ciliwung Hulu 4 131 1300 8 Ciesek 12557 14800 9 Ciseuseupan 1190 12100 10 Ciliwung Hulu 5 595 4200
Panjang sungai utama km 14.75 16.5 2.3 14 1.05 15.3 1.3 14.8 12.1 4.2
1.
Model Sektor
initialization
Patch water balance
Stream network
Output conversion
Reservoir dynamic
Additional output
Gambar 7. Daigram alur model dalam stella 2.
Sub model Neraca Air
Jumlah curah hujan tiap – tiap jenis penutupan lahan pada setiap sub-DAS dihitung per unit area dan masing – masing fraksi area. Curah hujan yang jatuh pada sebuah DAS, setelah diuapkan sisanya akan mengalir ke sungai, biasa disebut hasil air (water yield). Neraca air sebuah DAS yang berhutan dapat digambarkan dengan persamaan matematika sebagai berikut:
3.
Sub model aliran sungai
Sub model aliran sungai merupakan sub model yang mengumpulkan dan menyusun proses – proses hidrologi dari sub model lain menjadi debit aliran sungai. Pada sub model ini, komponen – komponen tersebut membentuk debit aliran sungai untuk periode harian serta akumulasi harian selama setahun.
Pg = (T + Ic + If + Es + w) + Q + ∆S ± L + U Jika (T + Ic + If + Es + w) sama dengan total veapotranspirasi (Et) dan dianggap tidak ada kebocoran (L) dan aliran sungai bawah tanah (U), maka persamaan neraca air sebuah DAS dapat disederhanakan sebagai berikut (Manan, 1978 dan Ward 1975 dalam Bruijnzeel, 1982): Pg = Et + Q + ∆S Gambar 9. Diagram sub model aliran sungai (Van Noordwijk et al., 2004 dalam modul GenRiver) 3.3.5
Gambar 8. Diagram sub model neraca air (Van Noordwijk et al., 2004 dalam modul GenRiver)
Verifikasi dan Kalibrasi Model
Verifikasi model dalam penelitian ini dilakukan pada data tahun 1994 dan 2004. Proses verifikasi ini melibatkan kegiatan kalibrasi model untuk mendapatkan nilai koefisien determinasi tertinggi dan kesalahan minimum terrendah sehingga hasil model yang berupa data nilai debit mendekati data debit harian hasil pengukuran di lapangan. Proses kalibrasi dilakukan dengan memasukkan beberapa parameter secara coba
ulang (trial error). Dari proses ini didapatkan nilai koefisien determinasi yang menyatakan bahwa secara statistik hasil keluaran model telah menggambarkan kebenaran dengan nilai R-square berkisar antara 0 – 1. R-square yang semakin besar, yaitu mendekati 1 menunjukkan bahwa model yang dirumuskan untuk menjelaskan keragaman data sangat baik (Sutoyo, 1999). 3.3.6
Analisis Sensitivitas Model
Tujuan utama analisis ini pada proses permodelan adalah untuk menentukan peubah keputusan mana yang cukup penting untuk ditelaah lebih lanjut pada aplikasi model. Peubah keputusan ini dapat berupa parameter rancang bangun atau input peubah keputusan. Analisis ini mampu menghilangkan faktor yang kurang penting sehingga pemusatan studi lebih dapat ditekankan pada peubah keputusan kunci serta menaikkan efisiensi dari proses pengambilan keputusan (Suwarto, 2006). Menurut Syaifullah (2004) model dikatakan sensitive terhadap perubahan nilai parameter apabila perubahan nilai parameter sebesar X % diikuti oleh perubahan nilai keluaran model minimal dalam persentase yang sama. Model dikatakan tidak sensitive apabila persentase perubahan luaran model terlalu kecil dibandingkan persentase perubahan nilai parameter tersebut. Belum ada referensi nilai persentase tertentu dalam membandingkan perubahan luaran model GenRiver terhadap perubahan parameternya. 3.3.7
Uji Keabsahan Model
Tolok ukur uji keabsahan model didasarkan pada: 1.
2.
Penampilan hubungan antara debit dugaan dan debit aktual secara grafik sehingga dapat ditentukan nilai mutlak (maksimum dan minimum) data yang diperoleh. model (R2) Nilai koefisien (Wardhani, 2002) dengan persamaan:
R = 2
1−
[∑ (Y − y ) ] 2
i
∑Y
i
i
2
Dengan: Yi : debit aktual ke-i yi : debit model ke-i 3.3.8
Skenario Tutupan Lahan dan Curah Hujan
Penggunaan model dilakukan dengan menerapkan beberapa kemungkinan skenario untuk mengetahui penggunaan lahan yang paling optimal dalam menekan fluktuasi debit, maka skenario yang digunakan adalah perubahan luas penggunaan lahan. Untuk skenario curah hujan digunakan nilai curah hujan apabila naik 15 20% dari rata – rata dan turun 15 - 20% dari kondisi rata – rata normalnya untuk masing – masing skenario perubahan penggunaan lahan. Beberapa skenario yang dilakukan seluruhnya mengacu pada kondisi awal yaitu tahun 2004. Skenario 1 merubah tegalan dan perkebunan masing – masing 9,1% dan 4,3% menjadi lahan hutan sedangkan lainnya tetap. Sehingga hutan pada skenario 1 menjadi 35% Sedangkan skenario kedua merubah 17,8% lahan perkebunan dan 7,3% lahan sawah menjadi lahan hutan, serta merubah 9,1% tegalan menjadi hutan dan 1% menjadi pemukiman. Sehingga pada skenario 2 proporsi hutan menjadi 55,8%. Dengan asumsi bahwa kriteria debit harian maksimum yang normal ditetapkan berdasarkan batas debit harian maksimum normal dari bendungan katulampa, yaitu debit yang tidak melebihi 244 m3/detik. Debit yang melebihi angka tersebut sudah termasuk kategori status siaga I yaitu debit yang berada antara 244 m3/detik dan 411 m3/detik. Besarnya debit diantara nilai tersebut sudah memiliki potensi untuk menimbulkan banjir di Jakarta. Dalam skenario ini input data iklim dianggap tetap. Berdasarkan hal tersebut, maka kriteria penggunaan lahan yang optimal dalam menekan fluktuasi debit dalam penelitian ini adalah tipe penggunaan lahan yang mampu memberikan nilai KRS yang rendah. Kriteria lain yang digunakan dalam penentuan lahan yang paling optimal adalah debit harian maksimum yang tidak melebihi batas debit normal di bendung katulampa.
Komposisi Hutan Perkebunan Pemukiman Tegalan
Tabel 4. Skenario perubahan penggunaan lahan *Awal (%) Skenario 1 (%) Skenario 2 (%)
21.6 31.6 25.8 11.6 9.4
Sawah * awal: mengacu pada kondisi tahun 2004
IV.
KONDISI UMUM WILAYAH KAJIAN
4.1
Kondisi Umum DAS Ciliwung Hulu
Secara astronomis Sungai Ciliwung berada pada letak lintang dan bujur 6°05` 6°50` LS dan 106°40` - 107°00` BT. Sungai ini bermula (hulu) di Gunung Mandalawangi. (Talaga) dan bermuara (hilir) di Teluk Jakarta. Wilayah DAS dengan luas sekitar 322 km2 ini dibatasi oleh DAS Cisadane di sebelah barat dan DAS Citarum di sebelah Timur. Sungai ini mengalir dari arah Selatan ke Utara dengan bentuk melebar di bagian hulu dan menyempit di bagian hilir. Sungai ini mengalir melalui daerah – daerah yang termasuk wilayah administrasi: a) Kabupaten Bogor khususnya kecamatan Cisarua, Ciawi, Kedunghalang, Cibinong dan Cimanggis; b) Kotamadya Bogor; c) Kota Administratif Depok; dan d) wilayah DKI Jakarta. Bagian hulu merupakan pegunungan dan berada pada batas ketinggian 300 m sampai 3000 dpl. Dengan luas 146 km2 bagian hulu DAS ini meliputi kecamatan Cisarua, Ciawi, dan Kedunghalang yang dibatasi oleh bendungan Katulampa sebagai outletnya, serta dikelilingi oleh G. Gede, G. Pengrango, G. Hambalang, dan Megamendung. Bagian DAS hulu ini terdiri dari sepuluh anak sungai yaitu: Citamiang, Cimegamendung, Cilember, Ciesek, Cisukabirus, dan Ciseuseupan. 4.2
Iklim
Iklim di daerah penelitian tergolong ke dalam iklim tropika. Suhu berkisar antara 23-24 °C dengan kelembaban nisbi antara 7382 %. Radiasi surya minimum terjadi pada bulan Januari (27,36 %) dan maksimum pada bulan September (81,85 %). Rata-rata penguapan minimum sebesar 2,08 mm terjadi pada bulan Januari sedangkan rata-rata penguapan maksimum sebesar 3,56 mm pada bulan Oktober.
35
55.8
27.3
13.8
25.8
26.8
2.5
1.5
9.4
2.1
Menurut Klasifikasi Iklim Schmidt Ferguson dalam Handoko (1994), iklim Sub DAS Ciliwung Hulu adalah termasuk ke dalam Zona Agroklimat A yang berarti daerah sangat basah dengan vegetasi hutan hujan tropika. Klasifikasi ini ditentukan berdasar dari jumlah Bulan Basah (hujan bulanan jangka panjang >100 mm) dan Bulan Kering (hujan bulanan jangka panjang <60 mm). Klasifikasi iklim A karena daerah ini mempunyai bulan – bulan basah berturut – turut sepanjang tahun. Hal ini sesuai dengan kondisi daerah hulu yang selalu tertutup awan 4.3
Tanah
Tanah – tanah yang terbentuk umumnya berasal dari bahan induk abu volkan dan batuan piroklastik. Berdasarkan Peta Tanah Semi Detil Tahun 1992 skala 1:50.000 yang dikeluarkan oleh Pusat Penelitian Tanah dan Agroklimat, jenis tanah yang terdapat di daerah penelitian meliputi order Andisol, Ultisol, Inceptisol, dan Entisol yang masing-masing sebesar 38.9 %, 11 %, 48 %, dan 2,1%. DAS Ciliwung Bagian Hulu dibangun oleh formasi geologi vulkanik yaitu komplek utama Gunung Salak dan komplek Gunung Pangrango. Deskripsi Litologi Kawasan DAS Ciliwung Bagian Hulu adalah tufa glas lhitnik kristal, tufa fumice dan batu pasiran tufa, sedangkan kondisi fisiografi daerah kawasan DAS Ciliwung Bagian Hulu merupakan daerah pegunungan dan berbukit. Elevasi umumnya diatas 150 m dpl dan terdiri atas daerah lungur volkan tua dan muda. Bahan induk tanah yang terdapat di DAS Ciliwung Bagian Hulu adalah berupa tufa volkanik dan derivatifnya merupakan bahan dasar pembentuk tanah jenis tanah Latosol Coklat Kemerahan adalah jenis tanah yang dominan. Adanya pencampuran bahan vulkanik tua dan yang lebih muda memungkinkan terbentuknya jenis-jenis tanah lain yang berasosiasi dengan Latosol antara lain adalah tanah Andosol dan Regosol
Gambar 10. Peta lokasi kajian
V.
HASIL DAN PEMBAHASAN
5.1
Analisis Perubahan Penggunaan Lahan
Pola penggunaan/penutupan lahan di daerah penelitian hasil pengamatan tahun 1994 sampai dengan tahun 2005 masingmasing digambarkan pada peta-peta yang disajikan pada Lampiran. Berdasarkan petapeta tersebut, daerah penelitian memiliki luas total 14.617 hektar dengan 8 tipe penggunaan/penutupan lahan, yaitu hutan lebat, hutan semak/belukar, kebun campuran, kebun teh, lahan terbuka, pemukiman, sawah, dan tegalan/ladang. Luas masing-masing tipe penggunaan/penutupan lahan tersebut disajikan pada Tabel 5 Perubahan penggunaan/penutupan lahan dari tahun 1994 – 2004 dapat diamati
melalui proses tumpang susun (overlay) peta pada ArcView. Data lengkap perubahan tipe dan luas penggunaan lahan dapat dilihat pada lampiran 9. Berdasarkan analisis tersebut lahan yang paling banyak terkonversi adalah kebun teh dimana pada tahun 2004 kebun teh telah terkonversi menjadi hutan semak belukar, kebun campuran, lahan terbuka, pemukiman, sawah, dan tegalan dengan persentase perubahan tertinggi menjadi daerah tegalan. Menyusul lahan sawah yang terkonversi menjadi kebun campuran, kebun teh, daerah pemukiman, sawah, dan tegalan, dengan perubahan tertinggi menjadi daerah tegalan. Luas penggunaan/penutupan lahan lainnya, seperti hutan semak/belukar, kebun campuran, lahan terbuka, dan pemukiman relatif tidak terlalu besar perubahannya.
Tabel 5. Luas penggunaan/penutupan lahan DAS Ciliwung Hulu Laju ∆ 1994 2004 landuse perubahan/th Nama lahan luas (Ha) % Luas (Ha) % (%) (%) hutan lebat
3571.7
21.1
3403.43
20.1
-1
-0.1%
hutan semak
576.9
4.4
256.8
1.5
-2.9
-0.29%
kebun campuran
1843.5
10.8
1873.4
11
+0.2
+0.02%
kebun teh
4248.2
25.1
3478.3
20.6
-4.5
+0.45%
lahan terbuka
50
0.2
14.8
0.1
-0.1
-0.01%
pemukiman
3315.5
19.6
4379.8
25.8
+6.2
+0.62%
sawah tegalan
2889.4 421.6
17.1 2.5
1581.2 1966.6
9.4 11.6
-7.7 +9.1
-0.77% +0.91%
5.2
Analisis Curah Ciliwung Hulu
Hujan
DAS
DAS Ciliwung terletak di Pulau Jawa dengan curah hujan yang dipengaruhi oleh monsun Asia dan Australia. Pola grafik curah hujan bulanan di daerah ini menyerupai huruf ‘V’, dengan nilai curah hujan minimum terjadi pada bulan Juni sampai dengan Septembar, sedangkan curah hujan maksimum terjadi pada bulan Oktober sapmai dengan April. Faktor iklim lain yang dapat mempengaruhi pola curah hujan di Indonesia adalah siklus Hadley dan siklus Walker, kedua siklus tersebut terjadi karena perubahan tekanan pada belahan bumi utara dan selatan serta pada bagian barat dan timur.
hujan ynag dianggap dapat mewakili pada tiga ketinggian yang berbeda di DAS Ciliwung Hulu selama periode 1994 – 2005. Data harian tersebut kemudian dirata – ratakan menjadi data bulanan untuk mengetahui pola grafik curah hujan dari setiap setasiun. Ketiga stasiun hujan tersebut terletak mulai dari kawasan puncak Cisarua dengan ketinggian 1160 m dpl hingga daerah outlet katulampa dengan ketinggian 480 m dpl. Data selengkapnya dapat dilihat pada Tabel 6 dan Tabel 10 (Lampiran 3) untuk analisis sifat hujan wilayah. Pola grafik curah hujan bulanan tiap stasiun hujan pada daerah hulu DAS Ciliwung dapat dilihat pada Gambar11.
Data curah hujan harian yang diamati pada penelitian didapat dari 3 stasiun Tabel 6. Curah hujan bulanan rata – rata DAS Ciliwung Hulu (1993 – 2005) Stasiun Elevasi jan feb mar apr mei jun ags sep okt nov katulampa 480m dpl 490 438 411 369 300 229 220 240 410 437 Gn Mas 1160m dpl 626 615 354 322 223 180 124 165 265 393 Citeko 920m dpl 576 470 312 290 187 121 116 114 211 299
des 368 410 378
700 Curah Hujan (mm)
600 500 400 300 200 100 0 jan
feb
mar
apr
mei
jun
jul
ags
sep
okt
nov
des
Bulan Katulampa
Gunung Mas
Citeko
Gambar 11. Grafik curah hujan bulanan tiga stasiun di wilayah Ciliwung Hulu Dari Tabel 6 dan Gambar 9 dapat dilihat bahwa curah hujan daerah Gunung Mas lebih besar daripada daerah Citeko yang berada dibawahnya. Namun demikian daerah Katulampa cenderung memiliki curah hujan yang paling tinggi diantara keduanya. Menurut Badan Meteorologi dan Geofisika (BMG) bahwa bulan yang memiliki nilai curah hujan bulanan < 150 mm disebut bulan kering, sedangkan bulan dengan curah hujan ≥ 150 mm disebut bulan basah.
Pada daerah hulu Ciliwung cenderung lebih banyak mengalami bulan basah sepanjang tahun yang disebabkan oleh pengaruh dari lokasi wilayah yang berada pada dataran sedang hingga tinggi. Bulan – bulan basah ini berturut – turut terjdi mulai September hingga Mei. Hal tersebut terjadi karena selain faktor ketinggian pada bulan – bulan tersebut tekanan udara di daratan Asia lebih tinggi dibandingkan di daratan Australia sehingga angin muson barat yang berasal dari benua Asia bergerak menuju benua Australia.
Kejadian sebaliknya bulan kering yang lebih banyak terjadi di daerah Citeko dan Gunung Mas. Tekanan udara di daratan Asia pada periode Juni - September mulai melemah sehingga angin muson timur bergerak dari Australia menuju Asia. Massa udara yang dibawa angin tersebut umumnya tidak mengandung uap air bahkan memiliki sifat yang kering sehingga tidak potensial untuk mendatangkan hujan. Hujan yang terjadi pada periode tersebut hanya dipengaruhi oleh sumber – sumber lokal dan awan yang dibentuk karena keadaan orografik dan juga proses konveksi. Pada periode tersebut sebagian besar daerah di Pulau Jawa mengalami musim kemarau. 600
Curah Hujan (mm)
500 400 300 200 100 0 jan
feb mar
apr
mei
jun
jul
ags sep
okt
nov des
Bulan
Gambar 12. Grafik curah hujan wilayah bulanan rata – rata DAS Ciliwung Hulu 5.3
Analisis Evapotranspirasi
Hasil perhitungan metode Penman FAO keluaran software Irsis diketahui bahwa nilai evapotranspirasi bulanan dan tahunan untuk ketiga stasiun yang dianalisis tidak nenunjukkan perbedaan yang nyata. Untuk stasiun Citeko nilai evapotranspirasinya 117 mm/bulan dan 1410 mm/tahun, stasiun Gunung Mas 113 mm/bulan dan 1366/tahun, sedangkan stasiun Katulampa 118 mm/bulan dan 1417 mm/tahun. Perbedaan topografi yang relatif rendah untuk daerah Katulampa
menyebabkan daerah ini mrempunyai nilai evapotranspirasi yang relatif lebih tinggi karena penerimaan radiasi surya yang lebih banyak dengan jumlah curah hujan yang tinggi pula. Vegetasi di daerah ini pun sudah relatif berkurang sehingga penguapan air dari permukaan tanah akan lebih banyak menyokong dari proses evapotranspirasi. 160
Evapotranspirasi (mm/bulan)
Angin tersebut mengangkut massa udara yang mengandung uap air yang berasal dari proses penguapan di atas Samudera Hindia dan kemudian membentuk awan potensial yang merata dan menutupi sebagian besar pulau Jawa dan Sumatera. Awan tersebut sangat potensial untuk menghasilkan hujan, sehingga pada periode tersebut di Pulau Jawa dan sekitarnya mengalami musim penghujan.
140 120 100 80 60 40 20 0 1
2
3
4
5
6
7
8
9
10
11
12
Bulan Citeko
Gn Mas
Katulampa
Gambar 13. Grafik nilai evapotranspirasi 5.4
Analisis Debit Sungai Ciliwung Hulu (1993 – 2005)
Peningkatan debit sungai Ciliwung hulu terjadi pada musim hujan dengan nilai maksimum 195 m3/det. Penurunan debit sungai pada musim kemarau mencapai 0,02 m3/det yang terjadi pada minggu ke-2 Oktober tahun 1997. Penurunan debit minimum yang sangat drastis pada tahun 1997/1998 diduga kuat akibat adanya El Nino pada periode tersebut yang menyebabkan debit sungai menyusut pada musim kering. Selanjutnya analisis debit sungai dibagi menjadi dua periode pengamatan. Periode pertama mulai tahun 1993 sampai dengan 1996 sedangkan periode kedua dari tahun 1997 sampai dengan 2005. Pembagian ini didasarkan pada nilai koefisien rejim sungai yang mulai naik signifikan setelah tahun 1996. Untuk melihat kecenderungan grafik kenaikan debit selama periode tersebut maka data debit dipilah menurut kesamaan probabilitas atau kemungkinanan peluang kejadian suatu nilai melebihi suatu nilai tertentu yang dalam analsis ini digunakan peluang terlampaui di atas 75%. Berikut hasil dari analisis data debit sungai di outlet Katulampa.
250
0 20
200
150
60 80
100
CH (mm/hari)
Debit (m3/dt)
40
100 50 120 0
140 1
251
501
751
1001
1251
1501
1751
2001
2251
2501
2751
3001
3251
Hari (Julian date) CH Debit
Gambar 14. Grafik time series hubungan curah hujan dan debit harian meningkat tajam hingga 6,2%. Sehingga daerah hulu yang seharusnya menjadi daerah resapan air (catchment area) perlahan – lahan berubah fungsi menjadi daerah padat hunian akibat maraknya pembukaan lahan.
250
D eb it ( m 3/d t )
200 150 100
5.5
Analisis DAS
50
0
20
40
60
80
100
120
140
Hujan (mm/hari) per 1997-2005
per 1993-1996
Gambar 15. Hubungan antara curah hujan dengan distribusi debit periode 1993 - 2005. Curah hujan dan debit harian telah dipilah berdasarkan kesamaan exeedance probability (peluang kejadian suatu nilai melebihi suatu nilai tertentu). Secara umum peningkatan debit seiring dengan peningkatan curah hujan Terjadi peningkatan debit harian yang cukup signifikan antara periode tahun pertama (1993 – 1996) dan periode kedua (1997 – 2005). Pada periode tahun kedua rata – rata debit hariannya lebih tinggi dibandingkan periode pertama dengan puncaknya mencapai dua kali debit puncak periode tahun pertama. Debit harian rata – rata periode pertama sebesar 12,32 m3/det sedangkan periode kedua sebesar 18,41 m3/det. Hal ini dapat dilihat dari pengaruh ekosistem bagian hulu yang lahannya terus terkonversi untuk lahan hutan dan kebun sebagai penyangga air hujan pada musim hujan dan penyimpanan air pada saat musim kering. Tercatat untuk klasifikasi hutan baik hutan lebat maupun hutan semak telah berkurang 3,9%, kebun teh telah berkurang 4,5%, sdangkan pemukiman
Penyangga
Analisis lain yang dilakukan dengan data empiris curah hujan dan debit sungai Ciliwung hulu adalah aplikasi perhitungan kuantitatif dari beberapa indikator. Hasil perhitungan kuantitatif dari beberapa indikator fungsi hidrologi DAS disajikan pada Gambar 16. Indikator penyangga (buffering indicator) cenderung berkorelasi negatif dengan total debit sungai sehingga peningkatan debit akan menurunkan kapasitas menyangga dari sungai. Indikator penyangga menunjukkan tingkat penurunan yang relatif rendah pada kondisi puncak kejadian hujan (buffering peak events). 0.999
0.9985
0.998 BI
0
Indikator
0.9975
0.997 y = -0.6534x + 0.9999 R2 = 0.902 0.9965 0
0.001
0.002
0.003
0.004
0.005
TWY
Gambar 16a. Hubungan indikator penyangga terhadap TWY
tanam terhadap debit sungai terjadi karena vegetasi dapat menghalangi jalannya limpasan langsung permukaan dan memperbesar jumlah air yang tertahan di atas permukaan tanah sehingga akan menurunkan laju debit sungai.
1.000000 0.999000 0.998000
BPE
0.997000 0.996000 0.995000
y = -1.1931x + 1.0009 R2 = 0.4478
0.994000 0.993000 0.992000
5.6 0
0.001
0.002
0.003
0.004
0.005
TWY
Gambar 16b. Hubungan indikator penyangga puncak kejadian hujan terhadap TWY 0.6 0.5
RBI
0.4 0.3 0.2 y = 16.306x + 0.2553 R2 = 0.0543
0.1 0 0
0.001
0.002
0.003
0.004
0.005
TWY
Gambar 16c. Hubungan indikator penyangga relative terhadap total debit terhadap TWY Dari Gambar 16 terlihat bahwa indikator penyangga DAS berkorelasi negatif dan cenderung menurun terhadap nilai transmisi air pada saat bulan – bulan hujan. Berbeda pada indikator penyangga relatif terhadap total debit yang tidak memperhitungkan luasan area cenderung tidak terpengaruh terhadap perubahan nilai transmisi air (TWY) yang ditunjukkan dengan pola grafik yang tidak teratur. Hal ini sesuai dengan penelitian Van Noordwidjk et al., (2004) di DAS Sumber Jaya Lampung yang menunjukkan indikator penyangga DAS semakin menurun apabila korelasinya negatif terhadap transmisi air. Kenyataan ini menunjukkan bahwa kemampuan menyangga DAS semakin berkurang terhadap produksi limpasan pada musim hujan setiap tahunnya. Penurunan kemampuan menyangga ini dapat disebabkan oleh menurunnya daya dukung daerah disekitar aliran sungai yang telah berubah menjadi daerah pertanian dan pemukiman pada kurun waktu 1993 – 2005. Kemampuan suatu penggunaan lahan dalam menahan curah hujan dan mengurangi terjadinya debit sungai serta menyerapkan air ke dalam tanah dipengaruhi oleh kondisi permukaan tanah terutama vegetasi (Asdak 1995, Chapman et al., 2003) dan sifat tanah (Purwowidodo, 1999) dari lahan tersebut. Asdak (1995) menyatakan bahwa pengaruh vegetasi dan cara bercocok
Koefisien Rejim Sungai (KRS)
Kriteria lain yang digunakan untuk menilai fungsi hidrologi suatu DAS adalah dengan melihat nilai koefisien rejim sungai tiap tahunnya. Koefisien rejim sungai merupakan perbandingan antara debit harian rata – rata maksimum dan debit harian rata – rata minimum. Kecendrungan kenaikan nilai KRS menunjukkan bahwa fungsi hidrologi DAS semakin menurun demikian sebaliknya. semakin kecil koefisien ini berarti kondisi hidrologi dari suatu wilayah DAS semakin baik Batasan yang diberikan untuk menilai indikator ini adalah sebagai berikut (Asdak, 1995): KRS < 50 baik; 50 ≤ KRS < 120 sedang; KRS ≥ 120 buruk. 2500
2000 1500 KRS
0.991000
y = 25.503x 0.8151 R2 = 0.1267
1000 500
0 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 Tahun
Gambar 17. Grafik nilai KRS tahun 1993 – 2005 Dari analisis data didapatkan bahwa KRS mulai mengalami kenaikan yang cukup signifkan mulai tahun 1996 sampai dengan tahun 1999 dengan niai KRS tertinggi dicapai pada tahun 1998 sebesar 2040 hasil perbandingan dari debit maksimum 61,2 m3/det dan debit minimumnya 0,03 m3/det. Selisih yang sangat besar antara debit maksimum dan debit mimimum pada saat musim hujan dan musim kering mengindikasikan bahwa pada saat musim hujan debit sungai dapat mendatangkan banjir tetpai pada saat musim kemarau debit sungai sangat kecil hingga menyebabkan kekeringan disekitar daerah aliran sungai. Sepanjang rentang data pengamatan nilai KRS menunjukkan tahun 1993 -1995 koefisien baik karena nilai berkisar kurang dari 50, tahun 1996 – 2001 menunjukkan nilai KRS yang sangat buruk dengan kisaran 182 hingga 353,7, dan pada tahun 2002 – 2005
menunjukkan nilai KRS baik hingga sedang dengan nilai 36,6 sampai dengan 103,3. Simulasi Model GenRiver untuk Evaluasi dan Prediksi Debit Untuk mempelajari hubungan curah hujan, debit sungai dan alih guna lahan dilakukan simulasi model GenRiver menggunakan data-data daerah Ciliwung hulu. Untuk itu dilakukan simulasi model dengan komposisi 25,5% hutan pada awal simulasi dengan penurunan hingga 21,6% pada akhir simulasi selama periode 13 tahun.
Peningkatan luasan pemukiman dari 19,6% hingga 25,8% dengan penurunan luas areal pertanian dari 17,1% hingga 9,4%. Perbandingan debit dari data empiris (data pengukuran) dengan hasil simulasi model GenRiver untuk tahun ke -1(1994) dan ke11(2004) disajikan pada Gambar 18. Tahun ke – 1 mewakili kondisi awal simulasi (25,5% areal hutan) dan tahun ke-11 mewakili kondisi akhir simulasi dengan 21,6% areal hutan
35
Debit Aktual
30
Debit Prediksi
Debit (mm/hari)
25 20 15 10 5 0 1
18 35 52 69
86 103 120 137 154 171 188 205 222 239 256 273 290 307 324 341 358 Hari
Gambar 18a. Hasil simulasi GenRiver pada tahun 1994 60 55
Debit Aktual
Debit Prediksi
50 Debit (mm/hari)
45 40 35 30 25 20 15 10 5 0 1
18 35 52 69 86 103 120 137 154 171 188 205 222 239 256 273 290 307 324 341 358 Hari
Gambar 18b. Hasil simulasi GenRiver pada tahun 2004 Perbandingan hasil simulasi dengan data pengukuran tidak dapat dilakukan dengan melihat kedekatan setiap titik hasil simulasi dengan data pengukuran. Hasil tersebut secara umum berarti simulasi model dapat menghasilkan pola debit yang sama dengan data pengukuran walaupun masih belum bisa mendekati beberapa titik puncak dan aliran dasar. Dari uji keabsahan model simulasi pertama, dilihat dari kesesuaian pola bahwa debit prediksi mendekati pola debit pengukuran dengan nilai koefisien
determinasi model sebesar 0,55 yang berarti peluang data terwakili sebesar 55% sedangkan nilai rmsenya 1,7. Hal ini karena ragam (variance) dari debit prediksi yang lebih besar dari data pengukuran. Sedangkan untuk simulasi kedua, pola debit prediksi juga mendekati pola debit aktualnya dengan koefisien determinasi model (R2) mencapai 0,71 yang berarti peluang data terwakili sebesar 71% dan nilai rmsenya sebesar 2,7. Debit sungai pada tahun
ke-11 relatif lebih tinggi dibandingkan debit pada tahun ke-1(Gambar 18). Namun demikian peningkatan debit maksimum pada tahun ke-11 cukup signifikan dari tahun ke-1. Debit sungai tahun ke-11 lebih berfluktuasi jika dibandingkan dengan tahun ke-1 yang cenderung stabil setelah memasuki musim kering yaitu hari ke-201 hingga akhir tahun 1994.
meskipun waktu terjadinya hampir dapat bersamaan. Debit maksimum aktual dan prediksi tercapai pada tengah bulan Januari. 5.8
Analisis Sensitivitas Tutupan Lahan.
Parameter
Analisis sensitivitas dilakukan baik terhadap parameter tipe penggunaan lahan maupun terhadap peubah curah hujan. Parameter tipe penggunaan lahan yang dipilih adalah perubahan tipe tutupan lahan yaitu hutan, perkebunan, pemukiman, dan sawah. Hal ini dimaksudkan untuk melihat sensitivitas model terhadap perubahan tutupan lahan di suatu DAS akibat dari perubahan penggunaan lahan yang taerhadap nilai – nilai luarannya.
Pada simulasi tahun pertama debit maksimum tercapai pada besaran 44,40 m3/dt untuk data aktual sedangkan untuk debit prediksi mencapai 49,16 m3/dt. Debit minimum aktual 1.93 m3/dt dan prediksinya 1,13 m3/dt Perbedaan ini tidak terlalu jauh yang disebabkan oleh parameterisasi yang masih harus diuji lebih lanjut karena secara grafik untuk pola debit harian model mendekati pola debit harian aktualnya. Debit maksimum ini tercapai pada saat memasuki puncak musim hujan atau pada bulan Februari. Selanjutnya untuk simulasi tahun ke-11, debit maksimum aktual tercapai pada besaran 82,46 m3/dt sedangkan untuk debit prediksinya 99,36 m3/dt. Debit minimum aktual 1,20 dan prediksinya 2,35. Pada simulasi tahun ke-11 ini penyimpangan data aktual dan prediksinya juga tidak terlalu besar karena parameternya mendekati ketelitian
Langkah yang dilakukan dalam analisis sensitivitas parameter penggunaan lahan adalah menambah 5 – 10% dan mengurangi 5 – 10% dari masing – masing parameter. Penambahan dan pengurangan lahan 5% akan mengurangi dan menambah 1,25% pada setiap penggunaan lahan lainnya sedangkan penambahan dan pengurangan lahan 10% akan mengurangi dan menambah 2,5% pada setiap penggunaan lahan lainnya. Model dijalankan dengan data parameter penggunaan lahan yang telah diubah.
Tabel 7. Perbandingan luaran model pada beberapa nilai paramater penggunaan lahan dengan luaran model dan nilai awal. luaran model pada perubahan hutan (%) +5 +10
0 var luaran Total debit
∆ (mm)
(mm)
1432
1349
0
(mm)
Total debit
%
mm
(mm)
1432
1494
0
(mm)
mm
mm
(mm)
1432
1380
(mm) -80
%
mm
mm
%
97
6.8
1608
176
12.3
-5
-10 ∆
%
mm
62 4.3 1572 140 9.8 1375 luaran model pada perubahan perkebunan (%) +5 +10 ∆
(mm)
∆
mm
∆ %
-10 ∆
mm
∆ (mm)
var luaran Total debit
∆
% mm -83 -5.8 1249 -183 12.8 1529 luaran model pada perubahan pemukiman (%) +5 +10
var luaran
-5
∆
mm
%
mm
mm
%
-57
-4
1321
-111
-7.8
-5
∆
-10 ∆
∆
%
mm
mm
%
mm
mm
%
mm
mm
%
-5.6
1274
-158
-11
1504
72
5
1585
153
10.7
luaran model pada perubahan sawah (%) +5 +10
0 var luaran Total debit
∆ (mm)
(mm)
1432
1462
0
Total debit
∆
(mm)
%
mm
mm
∆ (mm)
(mm)
1432
1458
mm
26
%
mm
mm
%
-29
-2
1383
-49
-3.4
-5
-10 ∆
%
mm
mm
%
mm
mm
1.8
1469
37
2.6
1408
-24
Hasil analisis sensitivitas pada Tabel 7 terlihat bahwa variasi perubahan penggunaan/tutupan lahan selalu diikuti dengan perubahan luaran model yang berupa total debit selama setahun. Dari lima parameter yang dirubah untuk menganalisa sensitivitas model GenRiver perubahan hutan dan perkebunan yang paling sesuai dengan persentase perubahan parameter. Hal ini karena hutan dan perkebunan mempunyai kemampuan untuk menahan laju air hujan ketika jatuh ke permukaan tanah melalui kanopi pada bagian atas dan serasah pada bagian bawah. Dengan demikian air hujan tidak seluruhnya secara langsung akan membentuk limpasan akan tetapi sebagian tertahan oleh tajuk untuk kemudian diuapkan dan sebagian lagi terserap oleh akar – akar pohon besar dan terinfiltrasi pada saat curah hujan menurun. Kemampuan masing –masing lahan di sub DAS Ciliwung Hulu dalam menghasilkan debit terutama terkait dengan kemampuan masing – masing lahan tersebut dalam menghasilkan runoff. Ini disebabkan runoff merupakan komponen pembentuk debit yang memberikan sumbangan aliran air
∆
mm
∆
(mm)
-10 ∆
%
30 2.1 1472 42 2.9 1403 luaran model pada perubahan tegalan (%) +5 +10
var luaran
-5
∆ % 1.7
mm
mm
%
1392
-40
-2.8
paling besar dibandingkan dengan sumbangan air dari komponen pembentuk debit lainnya yaitu interflow dan baseflow. 5.9
Analisis Sensitivitas terhadap Curah hujan
Model
Pertimbangan ini dilakukan karena parameter curah hujan sangat berperan dalam proses kesetimbangan air. Curah hujan akan langsung berpengaruh pada besarnya aliran. Analisis sensitivitas dilakukan dengan melihat pengaruh perubahan luaran model apabila peubah ini divariasikan. Langkah yang dilakukan adalah dengan menambah 5 – 10%, dan mengurangi 5 – 10% dari peubah curah hujan sementara peubah yang lainya tetap. Model dijalankan dengan data yang telah dirubah. Luaran yang dibandingkan adalah total debit selama setahun. Data masukan yang akan digunakan adalah data simulasi tahun 2004 dengan parameter tanah dan penggunaan lahan yang telah dikalibrasi. Data curah hujan rata – rata normal digunakan untuk menjalankan model luarannya merupakan nilai awal.
Tabel 8. Perbandingan luaran model pada beberapa nilai peubah curah hujan dan nilai awal. luaran model pada perubahan hujan (%) 0 +5 +10 -5 -10 var luaran
∆ (mm)
Total debit
1432
(mm) 1528
(mm) 96
∆ %
mm
mm
%
mm
6.7
1618
186
13
1332
Berdasarkan hasil skenario di atas, kontribusi curah hujan terhadap perubahan luaran model cukup nyata. Kenaikan 5% hujan menyebabkan perubahan total debit
∆ mm % 100 - 7
∆ mm 1246
mm 186
% 13
sebesar 6,7%, kenaikan hujan 10% menaikkan 13% total debit. Sedangkan untuk penurunan hujan 5% telah menurunkan total debit hingga 4,9% dan penurunan hujan 10%
menurunkan total debit hingga 10%. Berdasarkan hasil tersebut dapat dikatakan model GenRiver lebih sensitive terhadap perubahan curah hujan apabila dibandingkan dengan perubahan penggunaan/tutupan lahan. Pengaruh curah hujan yang secara langsung terhadap kesetimbangan air di permukaan tanah menyebabkan peubah ini sangat sensitive mempengaruhi fluktuasi debit. Hasil Prediksi Debit menggunakan Model GenRiver Untuk melihat pengaruh alih guna lahan dan variasi curah hujan terhadap debit harian di DAS Ciliwung hulu dilakukan simulasi GenRiver. Penggunaan model dilakukan dengan menerapkan beberapa kemungkinan skenario untuk mengetahui penggunaan lahan yang paling optimal dalam menekan fluktuasi debit, skenario yang akan digunakan adalah perubahan luas penggunaan
lahan seperti yang tercantum pada Tabel 4 dengan skenario variasi curah hujan rata – rata normal, + 25%, dan -25%. Skenario pertama menggunakan perubahan lahan tegalan dan perkebunan masing – masing 9,1% dan 4,3% menjadi lahan hutan, sementara penggunaan lahan lainnya tetap. Dengan penambahan ini maka lahan hutan menjadi 35%. Sedangkan skenario kedua mengalihgunakan 8,6% lahan perkebunan, 12,2% pemukiman, 9,1% tegalan, dan 4,3% sawah menjadi lahan hutan. Dengan pengalihugunaan ini lahan hutan menjadi 55,8% dari area DAS. Secara konsep penambahan lahan hutan dapat menurunkan debit maksimum dan debit harian karena intersepsi oleh tajuk dan penyerapan oleh akar pohon – pohon besar. Berikut hasil skenario selengkapnya pada Tabel 9 sedangkan untuk gambar grafiknya ada di Lampiran 4.
Tabel 9. Debit hasil luaran model pada berbagai skenario lahan dan hujan (satuan dalam m3/dt) Skenario 1
Var. Luaran
Skenario 2
CH rata-rata
CH + 25%
CH -25%
CH rata-rata
CH +25%
CH -25%
Debit rata2
13.25
16.48
9.73
8.89
11.1
7.13
Debit max
92.43
113.84
69.32
61.44
76.73
49.31
Debit min
2.37
2.96
0.56
2.90
3.63
2.17
KRS
39
38.46
123.79
21.19
21.14
22.72
Penggunaan berbagai skenario memberikan perubahan yang nyata terhadap KRS dan debit maksimum. Namun perubahan tersebut kurang nyata terhadap debit minimum, sehingga nilai KRS lebih
banyak dipengaruhi oleh debit maksimum. Berikut grafik yang menggambarkan kecenderungan pada berbagai skenario yang telah dilakukan. 40 Discharge skenario (mm/day)
Dari hasil simulasi debit model GenRiver skenario 1 maka dapat diketahui besaran volume debit maksimum selama satu hari sebesar 7.985.952 m3 pada kondisi curah hujan normal dengan debit minimum 204.768 m3, 9.835.776 m3 dan debit minimum 255.744 m3 dengan penambahan hujan 25%, dan 5.989.248 m3 dan debit minimum 48.384 m3 dengan pengurangan hujan 25%. Sedangkan pada skenario 2 diperoleh volume debit maksimum pada kondisi hujan normal turun menjadi 6.447.168 m3 dan debit minimumnya 217.728, 7.285.248 m3 dan debit minimum 243.648 pada penambahan hujan 25%, dan 4.923.072 m3 dan debit minimumnya 185.760 m3 pada pengurangan hujan 25%.
35 30 25 20 15 10 5 0 0
5
10
15
20
25
Discharge aktual rata-rata (mm /day) rata-rata aktual SK2CHrata-rata
SK1CHrata-rata SK2CH-25%
SK1CH-25% SK2CH+25%
SK1CH+25%
Gambar 19. Plot hubungan debit model skenario terhadap debit model kondisi ratarata aktual. Skenario 1 dengan penambahan hujan 25% ternyata nilainya melebihi pada kondisi aktual normal. Sedangkan skenario lainnya mempunyai kecenderungan nilai
debit berada di bawah debit pada kondisi normal aktualnya. Sehingga dapat dipastikan bahwa skenario 1 dengan hujan naik 25% akan meningkatkan rata – rata debit normal aktualnya. Perubahan debit sebagai hasil perubahan luasan lahan serta penambahan dan pengurangan curah hujan pada masing – masing skenario dipengaruhi oleh kemampuan masing – masing lahan dalam menghasilkan debit. Tipe penggunaan lahan dan kondisinya mempengaruhi besarnya runoff melalui pengaruhnya terhadap laju infiltrasi tanah. Dedaunan dan serasah menjaga infiltrasi potensial tanah dengan mencegah penutupan permukaan tanah akibat tetesan air hujan. Sebagian tetesan air hujan tertahan pada permukaan daun, meningkatkan kemungkinan air hujan tersebut terevaporasi kembali ke atmosfer. Sebagian kelembaban yang diintersepsi akan sudah mengering selama mengalir dari tanaman ke bawah menuju tanah sehingga tidak memberikan kontribusi pada masa awal terjadinya runoff.
VI. 6.1
KESIMPULAN DAN SARAN Kesimpulan
1.
Perubahan lahan dari tahun 1994 – 2004 didominasi oleh tegalan yang meningkat luasannya hingga 9,1%, kedua areal sawah yang menurun hingga 7,7%, dan ketiga pemukiman yang meningkat hingga 6,2% dari total luas lahan.
2.
Hasil perhitungan dan analisis yang dilakukan dalam penelitian ini menunjukkan penurunan indikator penyangga (buffering indicator). Penurunan indikator penyangga dari tahun 1993 – 2005 menurun dari 0,998883 sampai 0,996705 cenderung diikuti dengan peningkatan nilai transmisi air dari 0,001538 sampai dengan 0,004708.
3.
Nilai koefisien rejim sungai (KRS) memiliki kecenderungsn naik selama periode tahun 1993 – 2005 yang mengindikasikan DAS Ciliwung Hulu semakin kritis.
4.
Model GenRiver cukup peka terhadap variasi perubahan penggunaan lahan khususnya hutan dan perkebunan dan variasi hujan.
6.2
Saran 1.
Hasil simulasi model GenRiver untuk tahun pertama (1994) masih harus dilakukan parameterisasi lebih lanjut sehingga didapatkan tidak hanya model yang mendekati pola aktual tetapi secara statistik juga memiliki nilai korelasi yang tinggi.
2.
Nilai debit hasil simulasi untuk setiap periode cenderung selalu di atas nilai aktualnya (over estimate) sehingga untuk nilai debit skenario dapat juga tidak tepat.
DAFTAR PUSTAKA Asdak,
Bapeda
Chay. 1995. Hidrologi dan Pengelolaan Daerah Aliran Sungai. Gadjah Mada University. Yogyakarta. Prop. Jawa Barat. Kebiajakan Pengelolaan Sumber Daya Air di SWS Ciliwung – Cisadane untuk mengatasi Krisis Air Jakarta. http://www.air.bappenas.go.id/doc /pdf/makalah/Makalah%20Ciliwu ng.pdf (Diakses Oktober, 2006)
BAPPENAS. 2007. Laporan perkiraan kerusakan dan kerugian pasca bencana banjir awal feburari 2007 di wilayah JABODETABEK. Kementrian Negara Perencanaan Pembanngunan Nasional. Jakarta. Bruijnzeel, L.A. 1982. Hydrologycal and biogeochemical aspects of manmade forests in South-Central Java, Indonesia. Nuffic Project ITC/GUA/VU. Chapman, Dana C, dan Dorald W. Lake Jr. 2003. Section 4: Computing runoff. New York standards and specifications for erosion and sediment control. New York. Handoko. 1994. Klimatologi Dasar. Edisi ke2. Jakarta: PT Dunia Pustaka Jaya. Harimurti. 1999. Interpretasi Visual Foto Udara Digital pada Layar Monitor. Skripsi S1. Bogor: Jurusan Tanah, Fakultas Pertanian, Institut Pertanian Bogor. Tidak dipublikasikan. Irianto, S. 2001. Kajian Hidrologi Daerah Aliran Sungai Ciliwung menggunakan Model HEC-1. Tesis S2 Program Pascasarjana. IPB Bogor. Tidak dipublikasikan. Janudianto. 2003. Analisis Perubahan Penggunaan/Penutupan Lahan dan Pengaruhnya terhadap Debit Maksimum – Minimum di Sub DAS Ciliwung Hulu. Skripsi. Jurusan Tanah, Fakultas Pertanian, Institut Pertanian Bogor. Bogor. Tidak dipublikasikan. Koesmaryono, Y, Boer, R, Pawitan, H, Yusmin dan Las, I. 1999.
Pendekatan IPTEK dalam mengatasi penyimpangan iklim. Di dalam Prosiding Antisipasi El – Nino. PERHIMPI, Bogor. Manan, S. 1978. Kaidah dan pengertian dasar manajemen daerah aliran sungai. Didalam Manan, S. 1998. Hutan rimbawan dan masyarakat. Penerbit IPB Press, Bogor. Nijssen B, Lettenmaier DP, Liang X, Wetzel SW, Wood EF. 1997. Streamflow simulation for continental-scale river basins. Water Resources Research 33: 711-724. Pawitan, H. 2004. Metode Analisis Hidrologi. Di dalam Pelatihan Dosen tentang Teknologi Informasi untuk Manajemen Sumber Daya Alam. Bogor, 9 – 21 Agustus Prasatya, RD. 2006. Kajian Spasial Sebaran Vegetasi menggunakan Citra Ikonos dan Sistem Informasi Geografis, Studi Kasus di Sub DAS Ciliwung Hulu. Skripsi. Jurusan Manajemen Hutan, Fakultas Kehutanan IPB. Bogor. Tidak dipublikasikan Purwowidodo. 1999. Konservasi Tanah di Kawasan Hutan. Laboratorium Pengaruh Hutan. Jurusan Manajemen Hutan. Fakultas Kehutanan IPB. Bogor. Tidak dipublikasikan. Rogers, P. 1998. Hidrology and water quality in: changes in land use and land cover. Meyer, WB and BL Turner II, (Editors) A Global Perspective. The Press Syndicate of The University of Cambridge. Cambridge: pp 231-258 Sitorus, SRP. 2003. Pengembangan sumber lahan berkelanjutan. daya Laboratorium Perencanaan Pengembangan Sumber daya Lahan. Jurusan Tanah, Fakultas Pertanian, Institut Pertanian Bogor. Bogor. Sutoyo. 1999. Pendugaan Debit Sungai Hujan dengan Berdasarkan menggunakan Model Tangki di DAS Cidanau Serang. Skripsi. Jurusan Teknik Pertanian, Fakultas Tegnologi Pertanian,
Institut Pertanian Bogor. Bogor. Tidak dipublikasikan. Suwarto. 2006. Sistem dan model. Makalah pada Pelatihan Penyusunan Produk Domestik Regional Bruto (Pdrb) Hijau dan Perencanaan Kehutanan Berbasis Pentaan Ruang. Badan Planologi Kehutunan, Departemen Kehutanan. Bogor 4 – 9 Juni 2006. Syaifullah, MD. 2004. Kalibrasi dan Analisis Sensitivitas Model SEWAB di Sub DAS Sengguruh Brantas. Tesis. Sekolah Pasca Sarjana, Institut Pertanian Bogor. Bogor. Tidak dipublikasikan. Trisnadi, D. 2006. Optimasi Penggunaan Lahan dengan menggunakan Model Simulasi Hidrologi Daerah Aliran Sungai (Studi Kasus di Sub DAS Ciliwung Hulu, Jawa Barat). Skripsi. Jurusan Manajemen Hutan, Fakultas Kehutanan, Institut Pertanian Bogor. Bogor. Tidak dipublikasikan. Van Noordwijk M, Farida, A. 2004. Analisis debit sungai akibat alih guna lahan dan aplikasi model genriver pada DAS Way Besai, SumberJaya. Jurnal AGRIVITA (26) (1): 39 – 47. Visa,
J,
Komarudin, R, Ratag, M, Prabotosari, SK, Siswanto, B, Halimurrahman. 2006. Perubahan Klimatologis curah hujan di beberapa wilayah Indonesia ( Ciamis, Kerawang, Cilacap, Banyuwangi dan Maros ). Pusat Pemanfaatan Sains Atmosfer dan Iklim, LAPAN. Bandung. Dalam http://www.bdg.lapan.go.id/index. php?nama=reinstra&opt=detail&i d=37 (Diakses November, 2006)
Vink, APA. 1975. Land use in advancing agriculture. Springer Verlaag. New York. 394p. Wahana
Lingkungan Hidup Indonesia Banjir Bandang [WALHI]. Bahorok Akibat Rusaknya Hutan. www.walhi.or.id/Indonesia/siaran %20pers/siarpers_bjrbdg_akbt_rs k_HUTAN_04112003.htm [29 Desember 2003].
Wardhani, AN. 2002. Pendugaan Kualitas Air Sungai di Daerah Aliran Sungai Cidanau, Banten. Skripsi. Jurusan Teknik Pertanian, Fakultas Teknologi Pertanian IPB. Bogor. Tidak dipublikasikan. Wahyunto, Abidin, M Z., Priyono, A, dan 2001. Studi Sunaryanto. perubahan penggunaan lahan DAS Citarik, Jawa Barat dan DAS Garang, Jawa Timur. Makalah Seminar Nasional Multifungsi Lahan Sawah. Asean Secretariate Maff Japan dan Puslitbang Tanah dan Agroklimat. Bogor. Mulyana, N. Analisis karakteristik banjir di DKI Jakarta dan alternatif penanggulangannya. Makalah Pribadi dalam Pengantar ke Falsafah Sains. Sekolah Pasca Sarjana, Institut Pertanian Bogor. (www.tumoutou.net/pps702_9145 /nana_mulyana.pdf (diakses April, 2006) http://www.tempointeraktif.com/hg/nasional/ 2007/02/09/brk,2007020992838,i d.html (diakses Maret 2007)
Lampiran 1
Lampiran 2
Lampiran 3 Tabel 10. Analisis sifat curah hujan bulanan wilayah DAS Ciliwung hulu 1 2 3 4 5 6 7 8 9 658 N 383 BN 187 BN 208 BN 268 N 171 N 142 N 289 AN 160 476 BN 603 AN 279 BN 407 AN 126 BN 83 BN 109 BN 121 BN 200 515 N 543 N 349 N 337 N 141 BN 199 N 94 BN 208 AN 236 609 N 488 N 401 N 343 N 196 BN 220 AN 204 AN 223 AN 213 723 AN 394 BN 272 BN 303 N 156 BN 47 BN 107 BN 79 BN 161 534 N 604 AN 399 N 304 N 246 N 211 N 143 N 95 BN 241 496 N 497 N 481 AN 343 N 257 N 159 BN 172 N 229 AN 210 558 N 424 BN 294 BN 303 N 377 AN 192 N 118 BN 94 BN 17 448 BN 436 N 563 AN 414 AN 330 AN 361 AN 246 AN 239 AN 155 596 N 545 N 345 N 289 N 253 N 191 N 143 N 64 BN 107 516 N 440 N 312 BN 289 N 221 N 121 N 219 AN 134 N 106 732 AN 553 N 415 N 467 AN 506 AN 109 BN 133 BN 25 BN 408 676 AN 730 AN 637 AN 302 N 167 BN 406 AN 218 AN 216 AN 252 580 511 380 332 250 190 158 155 190 rata2 94.0 99.2 124.3 66.5 106.1 100.8 49.5 82.5 93.0 St Dev 434 323 282 212 161 134 132 161 85% 493 115% 667 587 436 381 287 218 181 178 218 Ket: 85% ≤ CH < 115% = CH Normal (N) CH < 85% = CH Di bawah Normal (BN) CH > 115% = CH Di atas normal (AN)
Tahun 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005
B N AN N N AN N BN BN BN BN AN AN
10 229 251 465 264 272 374 261 227 313 471 202 290 329 304 86.0 258 349
BN BN AN N N AN N BN N AN BN N N
11 195 334 339 460 378 457 336 511 337 417 312 782 400 404 138.9 344 465
BN BN BN N N N BN A BN N BN A N
12 495 515 531 470 355 448 348 475 128 240 277 501 235 386 131.2 328 444
AN AN AN AN N AN N AN BN BN BN AN BN
0.0
60
0.0
50
10.0
50
10.0
40
20.0
40
20.0
30
30.0
30
30.0
20
40.0
20
40.0
10
50.0
10
50.0
0
60.0
0
1
51
101
151
201
251
301
D e b it (m m /h a ri)
60
351
60.0 1
51
101
151
201
Hari
251
301
351
Hari
CH (mm/hari)
CH
Debit (mm/hari)
Debit luaran model dengan hujan rata-rata 45
Debit
Debit luaran model dengan hujan rata-rata
0.00
40
0.00
40
10.00
35
10.00
25
30.00
20 15
40.00
10
20.00 20
30.00 40.00
10
50.00
5 0
50.00
60.00 1
51
101
151
201
251
301
C H (m m /h a ri)
20.00
D e b it ( m m /h a ri)
30
30
C H (m m /h ari)
D eb it (m m /h ari)
C H (m m /h a ri)
Skenario 2
CH (m m /hari)
Deb it (m m /h ari)
Lampiran 4 Skenario 1
0
351
60.00 1
51
101
151
201
Hari
251
301
351
Hari
CH (mm/hari)
Debit (mm/hari)
CH
Debit luaran model dengan pangurangan hujan 25% 80 70
Debit
Debit luaran model dengan pengurangan hujan 25%
0.00
60
0.00
10.00
50
10.00
40
20.00
30
30.00
20
40.00
50.00
10
50.00
60.00
0
30.00
40 30
40.00
20 10 0 1
51
101
151
201
251
301
351
Hari CH (mm/hari)
C H ( m m /h a ri)
50
D e b it ( m m /h a ri)
20.00
C H (m m /h ari)
Deb it (m m /h ari)
60
60.00 1
51
101
151
201
251
301
351
Hari Debit (mm/hari)
Debit luaran model dengan penambahan hujan 25%
CH
Debit
Debit luaran model dengan penambahan hujan 25%
Lampiran 5 Model GenRiverVersion 1.0 1. Tampilan Model pada Excel File Excel berisi parameter – parameter input yang akan terhubung dengan model GenRiver versi 1.0 yang berjalan pada Stella versi 5 atau lebih ataupun versi regular.
•
Sheets Excel “Rainfall Database” dan “Debit Database” berisi data curah hujan dan debit dari DAS Ciliwung hulu dari tahun 1993 - 2004 (sepanjang tahun simulasi). Sheets Excel “Landcover Database Data” berisi properti dari 11 tipe penutupan lahan yang berhubungan dengan pola evapotranspirasi bulanan serta struktur tanah. Jumlah tipe penutupan lahan ini dapat disesuikan pada kondisi DAS setempat. Dalam hal ini DAS Ciliwung hanya menggunakan 5 tipe penutupan lahan sehingga tipe yang lain dikosongkan.
•
Sheets Excel “Subcatchinfo” berisi sub – sub DAS yang me-routing jarak dari tiap – tiap sub DAS menuju outlet dan bagian – bagian area dengan tipe tanah yang berbeda.
Lapmiran 6
•
2.
Sheets Excel “Soilproperties” menyediakan pendugaan nilai kapasitas lapang, titik layu permanent, dan kejenuhan air pada 10 kelompok tanah yang diturunkan dari database tanah. Tampilan Model pada Stella • Input Parameter Model GenRiver mempunyai ekstension .STM yang dapat dijalankan pada Stella versi 5 atau lebih mempunyai komponen – komponen berikut: bagian initialisasi, sektor dinamik, (berhubungan dengan kesetimbangan air), aliran sungai, dan jalur operasional reservoir didalam jaringan sungai, serta dua sektor yang menjaga jalur seluruh luaran parameter – parameter yang dibutuhkan.
•
Pengguna (User) dapat mendefinisikan panjang simulasi dan awal serta akhir periode pengukuran. Modifikasi input parameter dilakukan dengan menggeser slider ke kanan atau ke kiri sehingga nilainya berubah. Run dan Output Pengaturan input selesai menuju pada eksekusi model. Tampilan akan menuju pada halaman “Run & Output”.
Lampiran 7
Apabila hasil luaran model masih jauh dari yang diharapkan, kalibrasi dapat dilakukan dengan kembali pada halaman “Input Parameter” dan atur ulangi pada bagian input tanah dan intensitas hujan sampai hasil model mendekati nilai aktualnya. Selanjutnya output dari model dapat disajikan dalam bentuk grafik dan angka yang dapat dieksport ke file Excel. Satuan unit output dalam mm/hari.
Lampiran 8 Tabel 11. Contoh Perhitungan Indikator Fungsi Hidrologi (Th 1994) kejadian Q (m3/dt) Qabavg P abavg Q (mm/hari) P (mm/hari) 1 8.78 5.20 24.0 12.25 90.0 2 7.74 4.58 26.0 9.71 73.0 3 7.94 4.70 16.0 9.17 71.0 4 7.35 4.35 29.0 9.17 70.0 5 5.80 3.43 13.0 9.00 60.0 6 5.72 3.38 40.0 8.82 60.0 7 9.44 5.59 25.0 8.46 60.0 8 9.44 5.59 0.0 7.46 58.0 9 6.96 4.12 0.0 7.46 58.0 10 6.24 3.69 6.5 7.46 57.5 11 9.44 5.59 21.0 7.16 56.0 12 8.57 5.07 33.0 7.16 55.0 13 10.10 5.98 41.0 7.16 51.0 14 7.97 4.72 6.5 6.86 51.0 15 7.15 4.23 0.0 6.39 50.0 16 7.54 4.46 24.0 6.27 50.0 17 9.67 5.72 18.0 6.27 49.0 18 7.54 4.46 23.0 5.98 41.0 19 7.35 4.35 28.0 5.98 40.0 20 10.10 5.98 58.0 5.72 40.0 21 20.70 12.25 57.5 5.72 40.0 22 15.20 9.00 8.0 5.72 40.0 23 11.60 6.86 11.0 5.59 39.0 -
177 178 179 180
2.09 2.29 2.38 2.38
1.24 1.36 1.41 1.41
PabAvg = ∑ max(P-Pmean,0) = 2593.5 QabAvg = ∑ max(Q-Qmean,0) = 423.03 • Total debit sungai per unit hujan (TWY) TWY = ∑ Q /( A * ∑ P) •
= 0.00153812 Buffering Indicator (BI)
BI = ( PabAvg − (QabAvg / A)) / PabAvg = 1 − QabAvg /( A * PabAvg ) = 0.9988828
8.0 0.0 0.0 0.0 •
1.36 1.36 1.24 1.24
0.0 0.0 0.0 0.0
Relative Buffering Indicator (RBI)
RBI = 1 − (QabAvg / Qtot ) /(PabAvg / Ptot ) = 0.27365952 •
Buffering Peak Event (BPE)
BPE= 1− Max(daily_ Q − Qmean) /(A* Max(daily_ P − Pmean))
= 0.99906774
Lampiran 9 Tabel 12. Perubahan dominan tipe dan luas penggunaan lahan di sub DAS 1994 – 2004 Tipe penggunaan lahan Luas (Ha) 1994 2004 Sawah Tegalan 667,97 Sawah Permukiman 502,08 Kebun teh Tegalan 473,03 Sawah Kebun campuran 294,88 Kebun teh Permukiman 275,76 Kebun campuran Sawah 247,54 Kebun campuran Tegalan 202,89 Kebun teh Kebun campuran 170,17 Hutan semak/belukar Tegalan 143,79 Hutan lebat Kebun teh 135,50 Kebun campuran Permukiman 120,76 Tegalan Sawah 107,67 Sawah Kebun teh 68,20 Hutan semak/belukar Kebun teh 58,08 Tegalan Kebun campuran 53,70 Hutan semak/belukar Kebun campuran 46,99 Lahan terbuka Kebun campuran 43,69 Kebun teh Sawah 43,46 Kebun campuran Kebun teh 40,74 Tegalan Permukiman 28,20 Tegalan Kebun teh 26,39 Kebun teh Hutan semak/belukar 21,08 Kebun teh Lahan terbuka 9,80 Hutan lebat Hutan semak/belukar 8,74 Hutan semak/belukar Sawah 7,94 Hutan semak/belukar Permukiman 6,39 Hutan lebat Permukiman 5,26 Kebun campuran Lahan terbuka 1,90 Lahan terbuka Permukiman 0,42 Lahan terbuka Tegalan 0,33 Permukiman Kebun campuran 0,01 Hutan lebat Kebun campuran 0,00 Kebun teh Hutan lebat 0,00 Total 3813,36
Ciliwung Hulu tahun
% 17,52 13,17 12,40 7,73 7,23 6,49 5,32 4,46 3,77 3,55 3,17 2,82 1,79 1,52 1,41 1,23 1,15 1,14 1,07 0,74 0,69 0,55 0,26 0,23 0,21 0,17 0,14 0,05 0,01 0,01 0,00 0,00 0,00 100,00
Lampiran 10 Tabel 13. Data debit harian bendung katulampa 1994 (dalam m3/dt) Tanggal
1994 Jan
Feb
Maret
April
Mei
Juni
Juli
Agust
Sept
Okt
Nov
Des
1
8.78
10.8
7.74
5.72
4.61
3.39
2.29
1.91
1.58
1.58
1.83
4.92
2
7.74
10.6
7.35
4.76
4.32
3.27
2.19
1.83
1.58
1.58
1.83
5.55
3
7.94
12.1
8.36
5.72
4.04
3.15
2.19
1.91
1.58
1.58
4.32
4.18
4
7.35
9.67
7.94
6.24
5.55
3.15
2.38
2
1.66
1.58
2.92
3.65
5
5.8
15.5
6.24
7.15
5.89
3.03
2.29
2
1.66
1.58
2.38
4.47
6
5.72
9.67
5.55
6.78
5.89
3.03
2.29
2.09
1.58
1.5
2.29
5.55
7
9.44
8.57
6.24
5.23
6.06
3.27
2.29
2
1.74
1.58
2.38
5.89
8
9.44
10.6
6.78
4.92
5.89
3.78
2.19
1.91
1.58
1.83
2.09
5.55
9
6.96
8.36
7.74
5.39
6.59
2.92
2.09
1.83
2
2.81
2
6.41
10
6.24
7.94
7.15
4.61
9.44
2.81
2.19
1.74
1.74
2.19
1.91
10.6
11
9.44
7.54
6.78
4.18
7.54
2.7
2.29
1.74
1.83
2
1.91
6.06
12
8.57
8.57
6.24
4.04
6.24
2.59
2.19
1.66
1.74
1.91
1.91
4.76
13
10.1
8.57
6.41
6.96
7.54
2.59
2.19
1.66
2.09
3.78
2.29
5.07
14
7.97
8.57
6.78
5.39
6.06
2.7
2.19
1.66
1.83
2.49
2
5.39
15
7.15
9
5.72
8.15
5.72
2.7
2.09
1.66
1.74
2.09
3.52
7.35
16
7.54
8.57
5.23
6.06
5.39
2.82
2
1.66
1.66
3.39
2.59
5.72
17
9.67
7.15
6.59
5.55
4.61
2.81
2.09
1.66
1.66
2.49
3.15
6.41
18
7.54
6.59
7.74
5.39
4.32
2.59
2.19
1.66
1.83
2.7
3.39
5.07
19
7.35
7.35
8.36
5.39
4.18
2.49
2.09
1.66
1.74
2.29
3.39
4.18
20
10.1
6.06
7.74
5.89
3.91
2.49
2
1.66
1.74
2.09
3.03
3.91
21
20.7
5.72
5.55
6.06
3.91
2.38
2
1.58
1.74
2
3.03
4.32
22
15.2
6.78
5.72
5.72
3.65
2.38
2
1.66
2.49
1.91
3.15
4.04
23
11.6
5.55
6.24
7.94
3.52
2.29
2
1.58
2.49
1.91
3.03
3.65
24
14.3
5.35
9
8.15
3.52
2.29
2
1.58
1.91
1.83
3.03
3.52
25
14.9
5.23
6.59
7.15
3.39
2.09
2
1.58
1.83
1.91
3.15
3.39
26
16.4
7.74
7.35
6.41
3.39
2.09
1.91
1.58
1.83
1.83
4.92
3.52
27
12.1
12.6
5.55
5.89
3.27
2.29
2
1.5
1.83
2
3.15
3.52
28
15.5
7.54
5.07
5.55
3.03
2.38
1.91
1.58
1.66
1.91
6.78
3.52
29
12.6
5.72
5.07
3.03
2.38
1.91
1.66
1.74
1.83
9.44
3.52
30
12.6
5.5
5.07
3.03
2.29
1.91
1.58
1.58
1.74
6.96
3.03
31
12.1
5.07
2
1.58
Debit Rata2
10.29
8.51
6.65
2.11
1.72
2
Aliran km (I/det)
65.1
53.9
42.1
37.2
30.8
17.1
13.3
10.9
11.3
13
20.6
30.7
Tinggi aliran (mm)
174
130
113
96.5
82.5
44.4
35.7
29.2
29.3
34.8
53.5
82.3
Meter kubik (10**6)
27.6
20.6
17.8
15.3
13
7.02
5.65
4.61
4.63
5.5
8.45
13
Debit maksimum
10.6
3.27 5.88
4.86
2.70
1.74 1.79
2.05
3.78 3.26
4.85
20.7
15.5
9
8.15
9.44
3.78
2.38
2.09
2.49
3.78
9.44
Aliran km (I/det)
131
98
57
52
60
24
15
13
16
24
60
67
Tinggi aliran (mm)
350
237
153
134
160
62
40
35
41
64
155
180
Meter kubik (10**6)
56
38
24
21
25
10
6
6
6
10
24
28
2
Debit bulanan 2
318.8 238.3
206 176.5 150.8 81.14 65.35 53.36 53.66 63.65 97.77 150.5
Aliran km (I/det)
2018
1509
1305
1116
955
513
412
338
339
403
618
952
Tinggi aliran (mm)
5394
3640
3503
2895
2558
1332
1107
905
879
1079
1605
2551
Meter kubik (10**6)
856
577
552
459
403
211
175
143
139
171
254
403
Debit minimum
5.72
5.23
5.07
4.04
3.03
2.09
1.91
1.5
1.58
1.5
1.83
3.03
Rasio Q mak/Qmin
3.6
3.0
1.8
2.0
3.1
1.8
1.2
1.4
1.6
2.5
5.2
3.5
Lampiran 11 Tabel 14. Data debit harian bendung katulampa 2004 (dalam m 3/dt) Tanggal
2004 Jan
Feb
Maret
April
Mei
Juni
Juli
Agust
Sept
Okt
Nov
Des
1
4.6
10.5
26.7
4.22
13.6
5.27
11.1
2
4.6
9.18
27.6
4.18
15.7
5.02
5.27
1.85
2
1.54
12.1
11.3
1.7
2.15
1.69
13.4
3
8.96
8.85
22
4.98
9.41
5.26
10.1
4.02
1.39
2.46
1.39
12.8
9.48
4
5.92
9.73
19.6
4.22
6.91
5
5.04
11.7
20.4
4.38
10.9
5.52
3.38
1.38
2.17
1.28
13.5
13.3
5.28
3.95
1.83
2.16
1.98
13.2
6
5.75
9.55
16.8
4.99
15.3
9.82
6.51
4.37
1.55
2.01
1.55
10.8
7
5.78
9.75
13.3
16.2
7.49
8.27
7.85
4.59
1.54
1.7
2.14
9.8
12.6
8
7.27
12
9
6.09
12.4
11.8
5.11
5.63
6.36
4.4
1.39
1.39
2.01
9.47
9.28
9.55
4.61
5.28
5.56
3.99
1.38
1.28
1.85
9.16
10
6.04
12.5
12.3
17.5
4.2
4.42
5.03
3.58
1.68
1.18
2.45
7.95
14.1
11
4.65
12
5.49
11.1
17
3.98
5.23
4.61
3.56
2.44
2.12
3.72
7.9
12.9
15.6
13.6
5.71
7.25
4.6
4.16
2.02
1.41
2.37
6.12
13
16.4
10.2
14.6
25.6
5.29
5.6
5.49
3.59
3.5
1.28
3.11
6.04
19.9
14
13.7
18.9
30.1
7.5
6.27
7.51
4.95
1.28
1.18
2.65
5.05
14.4
15
12.2
16.3
19.1
7.34
7.29
9.08
4.81
1.08
1.18
4.32
9.28
19.8
16
11.1
15.8
14
8.77
4.9
7.95
5.5
1.08
1.37
19.4
10.7
17.5
17
11.4
15.8
27
7.94
4.01
9.69
4.63
1.37
1.38
9.58
9.51
16.5
18
10.7
7.67
13.4
5.37
3.78
10.1
3.8
1.38
2.73
8.25
12.3
12.3
19
11
5.35
11.5
4.82
3.97
10.4
4.96
1.38
1.34
7.61
11.4
9.55
20
10.1
4.42
10.4
4.8
3.58
9.19
4.61
1.28
1.53
9.98
12.4
9.16
21
10.4
9.85
10.1
5
10.4
9.44
4.4
1.18
1.39
9.48
9.58
11.6
22
11.7
7.69
9.78
4.41
8.31
4.39
4.79
1.18
1.38
13.3
7.96
11.5
23
13.4
18.2
9.47
3.99
9.11
3.59
3.61
1.18
2.13
10.8
10.9
9.82
24
18.9
21.7
7.17
2.38
9.15
4.16
3.36
1.18
1.56
12
8.62
10.1
25
23
20.1
5.34
2.62
8.25
2.24
3.16
1.82
1.39
10.5
10
8.85
26
21.4
25.4
3.03
2.18
7.91
2.16
1.9
1.85
1.28
12
10.7
11
27
22.6
19.8
4.13
2.16
5.62
2.76
2.3
1.3
2.27
8.96
21.4
11.1
28
21.4
27.3
3.98
3.9
7.26
2.93
2.17
1.98
2.47
10.6
16.8
10.1
29
20.5
4.17
7.45
7.08
3.73
2.16
2.3
1.57
10.1
14.3
10.1
30
19.2
4.78
8.78
5.84
4.76
2.01
2.47
1.69
12.7
10.2
13.4
31
20
5
2
1.87
4.03
1.64
Rata2 2
11.71 13.63 14.00
5.29 5.09
7.29
5.88
11.8
11.1
1.71
6.81
10.78 12.61
Aliran km (I/det)
50.7
59
60.6
22
31.6
25.5
17.5
7.1
7.38
29.5
46.7
Tinggi aliran (mm)
136
143
162
57.1
84.6
66
46.8
19
19.1
79
121
146
Meter kubik (10**6)
31.4
33
37.5
13.2
19.5
15.2
10.8
4.39
4.42
18.2
27.9
33.1 19.9
Debit maksimum
54.6
23
27.3
30.1
8.78
15.7
10.4
11.1
3.5
2.73
19.4
21.4
Aliran km (I/det)
100
118
130
38
68
45
48
15
12
84
93
86
Tinggi aliran (mm)
267
286
348
98
182
117
129
41
31
225
240
230
Meter kubik (10**6)
62
66
81
23
42
27
30
9
7
52
55
52
2
Debit bulanan
363.1 381.7 433.9 152.8
226
Aliran km (I/det)
1572
1652
1879
660
980
765
543
220
221
915
1401
1693
Tinggi aliran (mm)
4216
4004
5022
1713
2623
1980
1451
589
573
2449
3630
4526
Meter kubik (10**6)
973
924
1163
396
605
456
335
136
133
564
837
1026
Debit minimum
4.6
4.42
3.03
2.16
3.58
2.16
1.9
1.08
1.18
1.28
5.05
8.85
2
176.4 125.1 50.81 51.15 211.1 323.3
391
