RENCANA RESTORASI LAHAN RAWA

9 RENCANA RESTORASI LAHAN RAWA Wetland Restoration Plan Hairul Basri Program Studi Ilmu Tanah Fakultas Pertanian Unsyiah, Universitas Syiah Kuala, Darussalam Banda Aceh, Indonesia 23111, [email protected]

Ahmad Reza Kasuri Program Studi Teknik Sipil Fakultas Teknik Unsyiah, Universitas Syiah Kuala, Darussalam Banda Aceh, Indonesia 23111, [email protected] I. PENDAHULUAN A. Latar Belakang

L

ahan basah (wetlands) adalah salah satu ekosistem yang paling penting di bumi karena kondisi hidrologi yang unik dan perannya sebagai zona peralihan antara sistem daratan dan perairan (Mitsch dan Gosselink, 2011). Lahan basah, sebagai zona peralihan antara tanah dan air, memberikan perlindungan alami terhadap banjir ekstrim, sebagai simpanan air tawar, penyimpanan karbon jangka panjang serta memiliki keanekaragaman hayati yang tinggi. Namun, banyak lahan basah, khususnya dataran banjir sungai, delta dan muara, telah terdegradasi akibat aktivitas manusia (Verhoeven et al., 2006). Lahan rawa juga termasuk ke dalam lahan basah. Keberadaan rawa di kawasan pesisir Selatan Pulau Sumatera sangat dipengaruhi oleh faktor daratan dan Samudera Hindia. Faktor-faktor yang mempengaruhi pembentukan rawa di pesisir ini diantaranya vegetasi di pesisir, iklim, debit sungai, pergerakan sedimen pantai dan litologi. Pada beberapa kawasan pesisir yang memiliki potensi kejadian banjir yang dominan dan vegetasi yang mendukung, berpotensi membentuk rawa gambut. Gambut terjadi akibat akumulasi bahan organik yang melimpah yang melebihi proses pelapukannya dalam kondisi tergenang (an-aerob) yang menghambat terjadinya dekomposisi aerob sisa tumbuhan, kemudian menumpuk menjadi gambut. Dalam konteks hidrologi, rawa gambut memiliki peran yang sangat penting sebagai pengatur hidrologi. Gambut memiliki porositas yang tinggi sehingga mempunyai daya menyerap air yang sangat besar hingga 850% dari berat keringnya (Dreissen and Suhardjo, 1975). Kawasan gambut dan terutama kubah gambut berfungsi penambat (reservoir) air tawar dan melepasnya kembali secara perlahan‐lahan ke sungai. Pada musim kemarau kandungan air yang ada di bawah permukaan gambut akan terlepas secara perlahan, namun dengan debit air yang masih cukup besar dan pada saat musim hujan kawasan rawa gambut akan terendam/banjir kembali (Barkah dan Sidiq, 2009).

309

310

| SCIENTIFIC STUDIES FOR THE REHABILITATION AND MANAGEMENT OF THE TRIPA PEAT-SWAMP FOREST

Namun, fungsi hidrologi lahan rawa tersebut terus menurun seiring dengan semakin meningkatnya tuntutan ekonomi akan lahan yang memacu alih fungsi lahan lahan rawa menjadi lahan budidaya. Hutan Rawa Bergambut di Kawasan Rawa Tripa (Tripa Peat Swamp Forest = TPSF) terdapat di Provinsi Aceh atau tepatnya di Kecamatan Darul Makmur (Kabupaten Nagan Raya) dan Kecamatan Babahrot (Kabupaten Aceh Barat Daya).Ditinjau dari segi pengelolaan wilayah sungai, kawasan Rawa Tripa terletak di hilir Wilayah Sungai (WS)Tripa–Batee. Kondisi hidro-klimatologi Rawa Tripa sangat dipengaruhi oleh kondisi sungai-sungai yang mengalir melalui rawa dan Samudera Hindia. Imbuhan air kedalam sistem rawa sangat dipengaruhi oleh over bank flow banjir di sungai-sungai tersebut dan imbuhan akibat air hujan. Rawa gambut memiliki peran penting sebagai pengatur hidrologi karena berfungsi sebagai daerah penangkap air pada saat banjir dan kemudian melepaskannya secara perlahan pada saat musim kering. Endapan sedimen di areal TPSF umumnya terletak pada daerah rendah yang relatif datar dengan aliran sungai-sungai melalui rawa Tripa memiliki pola meander. Mata air sungaisungai ini sebagian besar bersumber dari daerah bukan gambut yang berada jauh di hulu. Sedangkan sumber air lainnya berasal dari kubah gambut yang mengalir melalui saluran-saluran alam atau buatan menuju sungai utama. Pola aliran alam yang ada di TPSF merupakan pola radial. Sedangkan saluran buatan yang dibangun untuk kepentingan perkebunan sawit dan memiliki kemampuan drainase yang besar. Di samping areal TPSF merupakan ekosistem air yang sangat penting sebagai pengatur hidrologi, Rawa Tripa juga mempunyai keanekaragaman hayati yang relatif tinggi. Terdapat aneka biotik dan non-biotik yang merupakan plasma nuftah endemik kawasan ini. Beberapa hewan dan tumbuhan yang ada di kawasan ini telah lama dimanfaatkan secara lestari oleh penduduk setempat sebagai sumber pangan dan sumber ekonomi. Alih fungsi lahan Gambut Rawa Tripa menjadi lahan budidaya pertanian untuk perkebunan kelapa sawit di satu sisi memiliki manfaat ekonomi kawasan, baik yang dilakukan oleh beberapa perkebunan yang beroperasi dan masyarakat di kawasan ini. Namun di sisi lain, aktivitas perkebunan tersebut menyebabkan fungsi lahan gambut sebagai pengatur hidrologi, sumber keanekaragaman hayati dan penyimpan karbon telah terganggu. Aktivitas perkebunan di kawasan rawa Tripa dengan membangun parit atau kanal cenderung tidak sesuai dengan kaidah konservasi karena memiliki kemampuan drainase yang besar. Hal ini telah mempercepat proses menurunnya secara berlebihan kandungan air di areal hutan rawa gambut. Penurunan air permukaan menyebabkan lahan gambut mengalami kekeringan dan pada kondisi kekeringan yang ekstrim, gambut akan sulit menyerap air kembali karena memiliki sifat kering tak balik (irresversible drying).Drainase berlebihan menyebabkan gambut akan mengalami subsidence (penurunan permukaan tanah), kondisi ini disebabkan oleh berkurangnya kandungan air dan mempercepat proses mineralisasi gambut akibat oksidasi. Selain itu, oksidasi tersebut meningkatkan emisi gas rumah kaca ke dalam atmosfir. Selanjutnya, pengeringan lahan gambut juga menyebabkan fungsi gambut sebagai penyimpan air menjadi terganggu. Oleh karena itu, areal hutan gambut rawa Tripa perlu direstorasi agar dapat berfungsi kembali sebagaimana yang diharapkan oleh semua pihak. Restorasi merupakan salah satu upaya untuk menyelamatkan dan memperbaiki kondisi ekosistem yang telah rusak (Ward, and Elliot, 1995). Restorasi dapat didefinisikan sebagai proses yang intens dalam membantu pemulihan dan pengelolaan integritas

LAPORAN UTAMA

RENCANA RESTORASI LAHAN RAWA

| 311

ekologi suatu ekosistem yang rusak termasuk berbagai variabel keragaman hayati penting, struktur dan proses-proses ekologi konteks sejarah dan kewilayahan, dan kelestarian praktik-praktik budaya (Clewell et al., 2005, Perrow & Davy, 2002). Tujuan utama restorasi lahan rawa adalah mengembalikan struktur dan fungsi ekosistem rawa awal. Namun, tujuan utama ini tidak dapat dilakukan dengan mudah dan dalam waktu singkat. Karena realita di lapangan, pemanfaatan lahan rawa sebagai areal perkebunan kelapa sawit telah berlangsung lama dan cenderung tidak ramah lingkungan dan menyebabkan terdegradasinya areal hutan gambut rawa Tripa (TPSF). Untuk itu perlu dilakukan kajian yang mendasar dan komprehensif terhadap dinamika ekologis rawa gambut di areal TPSF untuk mendukung upaya rencana restorasi areal TPSF. B. Tujuan Kajian rencana restorasi lahan rawa ini memiliki beberapa tujuan yaitu : 1. Analisis hidrologi sungai di TPSF ditujukan untuk memberikan gambaran tentang kondisi sungai-sungai yang melintasi TPSF, informasi tentang berbagai komponen hidrologi, pemetaan kondisi hidrologi, hubungan hujan-aliran, banjir rencana serta lama genangan di saat banjir,hubungan komponen hidrologi dengan emisi karbon. 2. Analisa hidraulika sungai di TPSF ditujukan untuk mengetahui karakteristik hidrolika sungai, parameter hidrolis sungai kondisi normal dan banjir tahunan, karakteristik sedimen, hubungan antar berbagai parameter hidrolika sungai, pembentukan meander sungai dan perkembangan konfigurasi dasar sungai. 3. Analisa kualitas air di TPSF ditujukan untuk memberikan informasi tentang kondisi kualitas air di sungai-sungai yang melintasi Rawa Tripa;kualitas air di saluran kebun sawit dan kualitas air di lahan dan genangan di Rawa Tripa; kondisi kualitas air berdasarkan baku mutu air. 4. Manajemen restorasi lahan rawa di TPSF ditujukan untuk memberikan informasi kondisi antropogenik; rencana restorasi rawa Tripa melalui rekayasa aliran sungai pada manajemen tata air rawa, upaya konservasi DAS dan lokasi utama yang perlu dikonservasi; manajemen air rawa dengan sistem drainase terkendali; 5. Monitoring dan evaluasi TPSF ditujukan untuk memberikan informasi pentingnya pemantauan kondisi aktual Rawa Tripa pasca kajian dan mengetahui realisasi rencana implementasi beberapa rekomendasi penting di masa yang akan datang untuk memperbaiki dan menjaga agar ekosistem TPSF menjadi lebih baik, lestari, berkelanjutan serta bermanfaat. C. Ruang Lingkup Studi Ruang lingkup kajian rencana restorasi lahan Rawa Tripa ini adalah sebagai berikut: 1. Analisis Hidrologi DAS di TPSF meliputi pembahasan karakteristik DAS, menganalis morfologi wilayah, analisis curah hujan wilayah; Analisis hujan-aliran, analisis banjir, pola hujan badai berbagai kala ulang, distribusi hujan jam-jam;Analisis curah hujan rancangan dengan periode kala ulang tertentu menggunakan berbagai distribusi seperti Gumbel Tipe I, Normal, Log Normal 3 dan Log – Pearson Tipe III; Analisa banjir rencana menggunakan metoda HSS Gama dan Nakayasu, kemungkinan banjir maksimum (PMF) dan analisis komponen hidrologi terhadap emisi karbon. 2. Analisa Hidraulika Sungai di TPSF mencakup karakteristik hidrolika sungai; Parameter hidrolis sungai kondisi normal dan banjir tahunan; Karakteristik sedimen; hubungan antar parameter hidrolika sungai seperti hubungan kecepatan geser dan debit,

Project Implementation Unit - Studi Ekosistem Rawa Tripa Universitas Syiah Kuala

312


hubungan kecepatan kritis dengan debit, hubungan bilangan Froude dengan kecepatan; pembentukan meander sungai; perkembangan konfigurasi dasar sungai. 3. Analisa Kualitas Air di TPSF melalui pengujian kualitas air di sungai-sungai Rawa Tripa, Seumayam, Batee Seunaam, Tripa;pengujian kualitas air di saluran perkebunan kelapa sawit, pengujian kualitas air di lahan dan genangan di areal Rawa Tripa; Pembahasan tentang kondisi kualitas air berdasarkan baku mutu air. 4. Manajemen Restorasi Lahan Rawa di TPSF meliputi kajian antropogenik; Rencana restorasi sistem tata air rawa melalui rekayasa aliran sungai pada manajemen tata air dengan skema recharge dan discharge aktualrawa Tripa pada musim kemarau, upaya menaikkan elevasi muka air hingga mencapai lapisan gambut pada saat tertentu melalui kontruksi bendung sederhana, skema recharge dan discharge rencana rawa Tripa pada musim kemarau setelah bendung dibuat, tipikal konstruksi bendung yang direkomendasikan untuk merekayasa elevasi aliran sungai di Rawa Tripa; Konservasi DAS membahas karakteristik tanah, kedalaman bahan gambut di lokasi rawa Tripa dan lokasi utama yang perlu dikonservasi;Uraian restorasi lahan gambut melalui manajemen air dengan sistem drainase terkendali. 5. Monitoring dan Evaluasi Kondisi Rawa Tripa membahas pentingnya aktivitas monitoring dan evaluasi kondisi Rawa Tripa pasca penelitian dan mengetahui realisasi rencana implementasi beberapa rekomendasi penting di masa yang akan datang. Hal ini dilakukan agar hasil kajian menjadi bermanfaat dan ekosistem TPSF menjadi lebih baik dan berkelanjutan. Output dari kajian restorasi lahan rawa adalah Rancangan Teknis untuk restorasi lahan rawa di areal TPSF. II. METODOLOGI PENELITIAN A. Tempat dan Waktu Kegiatan studi ini dilakukan di areal Hutan Rawa Gambut Tripa (TPSF) Provinsi Aceh seluas lebih kurang 60.657,29 hektar yang mencakup wilayah Kecamatan Darul Makmur Kabupaten Nagan Raya seluas 36.394 ha atau sekitar 60 persen, dan wilayah Kecamatan Babah Rot Kabupaten Aceh Barat Daya seluas 24.263 ha atau sekitar 40 persen dari luas areal. Luas kawasan rawa gambut Tripa tersebut termasuk dalam Kawasan Ekosistem Leuser (KEL). Secara geografis kawasan ini terletak pada 030 44’-030 56’ LU dan 960 23’ 960 46’ BT. Areal yang termasuk ke dalam wilayah studi adalah bagian hilir Wilayah Sungai (WS) Tripa-Batee. Kawasan WS Tripa Batee memiliki beberapa DAS dengan tiga sungai utama yaitu :DAS Krueng Tripa, DAS Krueng Seumayam dan DAS Krueng Batee. Kegiatan studi ini dilaksanakan mulai Mei sampai dengan Agustus 2013 dengan jumlah hari kalender sekitar 120 hari. B. Bahan dan Alat Analisa kondisi hidrologi untuk menyusun program restorasi lahan di kawasan gambut rawa Tripa membutuhkan ketersediaan data sekunder dan data pengamatan langsung. Data sekunder yang dibutuhkan bersumber dari stasiun pengumpulan data tetap. Sedangkan kajian terdahulu yang dijadikan rujukan bersumber dari instansi teknis terkait. Data dan peralatan yang dipergunakan pada kajian ini disajikan pada Tabel 1 dan Tabel 2. .

LAPORAN UTAMA


| 313

Tabel 1. Kebutuhan Data No

Jenis Data

1

Data Hujan

2

Data Debit

Sumber Data

Tahun Data

PT. Socfindo BPP Tadu BPP Darul Makmur Data Stasiun Beutong Data Alu Ie Mirah Kanwil PU/ADU

PT Socfindo BWS I - Sumatera BWS I - Sumatera

1996 – 2003, 2006 - 2012 2010 - 2012 2010 – 2012 1998, 2012-2013 2012-2013 Sesuai ketersedian data lapangan Sesuai ketersedian data lapangan 2005 – Jan 2013 2005 – 2012 2009

BAPPEDA Prov. Aceh Pengukuran Langsung Pengukuran Langsung Pengukuran Langsung Pengukuran Langsung

2007 2013 2013 2013 2013

BWS I - Sumatera 3 4 5 6 7 8 9 10

Data iklim Data kejadian banjir WS Tripa - Batee Pola Pengelolaan Wilayah Sungai Tripa - Batee Study ketersediaan Air Provinsi NAD Hidrometri Sungai Kualitas Air Sungai Hidrometri Kanal Rawa Kualitas Air Rawa

Tabel 2. Peralatan Survai yang Digunakan No

Peralatan

Unit

Vol

1 buah sewa

8 hari

1 unit sewa 2x18 botol air kemasan 1,5ltr

15 hari

1

Perahu untuk hidrometri/ Boats for hydrometric.

2

Water Quality Checker (WQC)

3

Air kemasan pembilas sensor WQC

4

Current meter

1 buah sewa

5

Echosounder

1 buah sewa

8 hari

6 7

Measure tape Tali Tambang

2 unit ls

2 unit 1 rol

8

Sedimen grap

1 unit sewa

15 hari

9

GPS

36 botol 15 hari

2 unit sewa

30 hari

10

Kontainer sampel + Analisis sedimen di laboratorium

unit

18 unit

11

Kontainer sampel + Analisis kualitas air di laboratorium

unit

18 unit

12

Mobil

unit

15 hari

13

Motor roda 2 Kontainer untuk pengangkutan sampel kualitas air (sungai dan rawa)

unit

13 hari

2 unit

2 unit

14

C. Metodologi Sesuai dengan tujuan dan ruang lingkup pekerjaan maka metodologi yang digunakan dalam studi ini secara umum adalah metode deskriptif kualitatif dan kuantitatif yaitu melalui kegiatan survai lapangan dan analisis laboratorium. Pada kegiatan analisa hidrologi, data utama yang dijadikan rujukan adalah data sekunder berupa data iklim dan karakteristik DAS. Pengumpulan data primer ditujukan untuk memverifikasi hasil analisa dan kebutuhan analisa yang lebih detail untuk hidraulika. Secara garis besar, kegiatan yang dilakukan oleh Tim Kajian Rencana Restorasi Lahan Rawa adalah sebagai berikut:


314


Kegiatan Persiapan Kegiatan ini bertujuan untuk mengidentifikasi dan merumuskan masalah secara lebih spesifik berdasarkan TOR dan studi terdahulu yang telah dilakukan, menyusun hipotesa dan menyusun rencana kerja. Kegiatan persiapan meliputi : (1) penyusunan rencana kerja, (2) studi Literatur dan (3) pengumpulan data sekunder. Penyusunan Rencana Kerja Penyusunan rencana kerja dilakukan berdasarkan penjabaran dari scope of work dan activities yang diharapkan di dalam TOR Wetland-Restoration Plan For In The TPSF Area. Berdasarkan TOR, cakupan analisa yang diharapkan meliputi dampak kondisi hidrologi terhadap tingkat cadangan karbon, emisi karbon, dan kondisi hidrologi terhadap tata guna lahan dan fungsi lahan. Studi Literatur. Studi literatur dilakukan dengan mengkaji laporan terdahulu terhadap kegiatan yang telah pernah dilakukan di Rawa Tripa. Selain itu, sebagai pengayaan referensi penyusunan teori, dilakukan kajian tentang studi lain yang berkaitan dengan rawa. Literature yang dijadikan rujukan merupakan literature yang memuat teori dan aplikasi mengenai pengembangan rawa, petunjuk teknis kegiatan pembangunan rawa, hidraulika sungai, kondisi hidrologi dan hidraulika rawa, antropogenik rawa, dan teori lainnya. Pengumpulan Data Sekunder. Pengumpulan data sekunder dilakukan di institusi yang terletak di Kota Banda Aceh. Pengumpulan data sekunder dilakukan dengan mengidentifikasi lokasi sumber-sumber data yang terkait langsung dengan pengembangan dan pembangunan di Rawa Tripa secara langsung dan tidak. Berdasarkan kebutuhan hidrologi, data yang dibutuhkan untuk menganalisa kondisi rawa Tripa meliputi data hidroklimatologi, karakteristik DAS, karakteristik sungai, kualitas air, morfologi rawa, sistem drainase, prasarana dasar pengairan. Data yang digunakan bersumber dari instansi teknis dan perkebunan yang terdapat di rawa Tripa dan DAS-DAS yang menjadi inflow aliran air ke rawa. Tahap Pra Survai Tahap pra-survai yaitu melakukan survai pendahuluan ke lokasi studi yang diikuti oleh seluruh Tim Ahli dan Ketua/Koordinator Tim untuk melakukan orientasi lapangan agar mendapatkan kepastian tentang rencana titik-titik pengamatan lapangan dan prediksi jumlah sampel dan/atau responden saat pelaksanaan survai utama. Kegiatan Pra Survai ini dilaksanakan sejak 27 April - 30 April 2013. Hasil kegiatan prasurvai ini telah didiskusikan kembali pada tanggal 3 Mei 2013 untuk perencanaan kegiatan selanjutnya terutama persiapan untuk survai utama. Survai utama dilaksanakan secara terpadu pada Minggu ke dua bulan Mei 2013. Survai lapangan dilaksanakan untuk mengumpulkan data fisik hidraulika sungai secara langsung yang akan dijadikan rujukan dalam perencanaan restorasi.Titik survai yang diidentifikasi pada tahapan pra survai disajikan pada Tabel 3.

LAPORAN UTAMA


| 315

Tabel 3. Titik Survai Hidrometri dan Kualitas Air No

Kode Titik

Koordinat

Deskripsi

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31

T-00 T-01 T-01a T-02 T-02a T-02b T-03 T-03a T-03b T-04 T-05 T-07 T-10 T-11 T-12 T-12a T-12b T-13 T-14 T-16 T-17 T-18 T-20 T-21 T-22 H-01 H-02 H-03 H-04 H-05 H-06

3°49'22.87"U ; 96°40'32.31"T 3°48'39.89"U ; 96°41'53.04"T 3°45'13.74"U ; 96°42'11.85"T 3°46'47.97"U ; 96°40'11.55"T 3°47'8.83"U ; 96°38'44.82"T 3°46'17.82"U ; 96°37'35.13"T 3°45'27.06"U ; 96°39'43.92"T 3°45'4.60"U ; 96°39'24.41"T 3°46'4.08"U ; 96°38'1.30"T 3°49'22.62"U ; 96°37'35.35"T 3°48'10.84"U ; 96°38'0.31"T 3°50'21.13"U ; 96°38'18.69"T 3°51'37.92"U ; 96°37'7.53"T 3°51'47.32"U ; 96°34'32.63"T 3°49'56.26"U ; 96°34'16.05"T 3°48'29.45"U ; 96°34'31.20"T 3°49'16.83"U ; 96°32'15.89"T 3°45'21.20"U ; 96°36'24.05"T 3°45'6.90"U ; 96°35'54.38"T 3°46'14.62"U ; 96°34'28.30"T 3°48'34.06"U ; 96°30'13.90"T 3°54'5.23"U ; 96°26'41.81"T 3°52'20.53"U ; 96°27'27.24"T 3°56'15.79"U ; 96°28'22.17"T 3°51'36.21"U ; 96°29'1.95"T 3°52'55.50"U ; 96°42'11.42"T 3°53'15.54"U ; 96°35'57.81"T 3°47'51.86"U ; 96°36'44.74"T 3°45'2.44"U ; 96°39'35.68"T 3°55'7.75"U ; 96°25'24.58"T 3°55'53.89"U ; 96°26'6.10"T

Saluran Buatan Kolam alam Kolam alam Saluran Buatan Saluran buatan Saluran buatan Saluran alam Kolam alam Kolam alam Kolam alam Saluran buatan Saluran buatan Saluran buatan Saluran buatan Saluran buatan Saluran buatan Saluran buatan Kolam alam Saluran buatan Saluran buatan Saluran alam Saluran buatan Saluran buatan Saluran buatan Saluran buatan Sungai Batee Sungai Seumanyam Sungai Seunaam Sungai Seunaam Sungai Tripa Sungai Tripa

Tahap Survai utama (a) Analisa Hidrologi Sungai Kegiatan ini menghasilkan gambaran tentang kondisi hidrologi di sungai-sungai yang menjadi melintasi Rawa Tripa (hulu - hilir), pemetaan kondisi hidrologi, pengaruh perubahan lahan terhadap kondisi hidrologi, tata cara pemantauan kondisi hidroklimatologi dan debit aliran. Kegiatan analisa hidrologi meliputi : pengumpulan data hujan dan iklim, analisa debit normal dan debit banjir dengan berbagai kala ulang, analisa imbangan air di rawa Tripa, dan analisa imbangan air di rawa Tripa. Pengumpulan data hujan dan iklim dilakukan untuk data yang tidak terdapat di Banda Aceh. Analisa debit normal dilakukan dengan menggunakan data perencanaan bangunan pelindung tebing sungai yang terdapat di instansi terkait. Data ini kemudian dikompilasikan dengan data AWLR yang ada dan data hasil pengukuran pertama hidrometri sesaat. Hasil analisa berupa kondisi hidrologi dan hidraulika sungai dikalibrasi melalui pengukuran hidrometri sesaat tahap ke dua. Analisa imbangan air dilakukan untuk melihat kondisi imbangan air di Rawa Tripa berdasarkan kondisi hidroklimatologi rawa. Kajian ini dilakukan dengan dua asumsi.


316


Asumsi pertama adalah sungai yang terdapat di sekitar rawa member pengaruh besar terhadap imbangan air di rawa. Asumsi kedua adalah keseimbangan air di rawa bersifat independen dan tidak terkait dengan sungai yang ada. Asumsi ini dilakukan untuk menyederhanakan model yang akan dibangun mengingat keterbatasan waktu yang diberikan. Pengukuran parameter fisik kualitas air rawa dan pengambilan sampel untuk mengetahui tingkat keasaman air akan menjadi indikator untuk melihat kondisi limpasan sungai. Secara teori, pada air permukaan rawa yang memiliki tingkat keasaman lebih tinggi dari air sungai dan air rawa dalam tanah mengindikasikan bahwa imbuhan air dari hujan lebih dominan. Sedangkan pada kondisi air permukaan rawa lebih rendah tingkat keasamannya dari air rawa dalam tanah mengindikasikan bahwa terjadi aliran air sungai yang dominan. (b) Analisa Hidraulika Sungai Rawa Tripa memiliki 3 sungai yang melintasi dan berhubungan langsung dengan sistem rawa. Untuk mengetahui sebesarap besar pengaruh inflow air sungai ke dalam sistem rawa, perlu dilakukan analisa kondisi hidraulika sungai. Pada analisa hidarulika sungai, parameter yang dijadikan rujukan adalah data aliran dengan periode minimal 10 tahun, tampang melintang dan memanjang sungai, kualitas air sungai dan pergerakan sedimen sungai. Data ini akan dipakai untuk menganalisa kondisi hidraulika saluran dan kemungkinan perubahan morfologi sungai. Pada kegiatan ini, tidak terdapat data dari stasiun pengamatan tetap. Untuk memenuhi kebutuhan minimal data analisa, dilakukan survai hidrometri sesaat. Tahapan pada kegiatan survai ini meliputi :Penentuan titik survai, Pengukuran hidrometri sesaat, Jenis sampel yang akan diambil dan Jenis parameter. Penentuan titik survai. Idealnya minimal titik survai pada analisa limpasan sungai ke rawa adalah 2, yaitu di hulu dan di hilir limpasan. Survai lapangan hidrologi akan melibatkan 2 komponen, yaitu hidrometri dan kualitas air. Survai ini dilakukan untuk kebutuhan studi secara umum dan bukan untuk mendapatkan data detail. Pengambilan titik sampel dilakukan pada setiap sungai yang melewati kawasan rawa Tripa. Terdapat 3 sungai utama yang menjadi lokasi pengamatan hidrometeri. Sedangkan untuk pengamatan di saluran, dipilih lokasi saluran induk drainase, saluran alam dan kolam alam yang masih ada. Titik survai hidrometri dan kualitas air dapat dilihat pada Tabel 3. Pengukuran hidrometri sesaat. Kegiatan dilakukan untuk mendapatkan data hidraulika sungai minimal. Kegiatan pengukuran hidrometri sesaat dilakukan sebanyak 2 kali. Pengukuran awal merupakan pengambilan data awal, dan pengukuran kedua merupakan proses kalibrasi. Item kegiatan yang dilakukan pada ke dua pengukuran sama. Jenis sampel yang diambil. Pada setiap titik pengamatan arus sungai, dilakukan pengambilan 1 buah sampel sedimen dasar sungai. Jadi untuk 2 kali pengukuran pada 3 sungai untuk 3 pias sehingga sampel sebanyak 18 buah sampel sedimen dasar sungai. Jumlah sampel sangat ditentukan oleh keseragaman morfologi dan tingkat kerapatan drainase. Sampel kualitas air sungai dan kualitas air rawa masing-masing sebanyak 18 sampel.

LAPORAN UTAMA


| 317

Jenis parameter (pengukuran insitu dan eksitu). Pengukuran langsung yang dilakukan di lapangan meliputi : Kecepatan aliran, pH air, Turbidity, Conductivity, Suhu dan Oksigen terlarut. Kegiatan Laboratorium Kegiatan ini meliputi pengujian sampel dari lapangan di laboratorium. Analisis Laboratorium meliputi : Uji butiran, sediment content, parameter fisik dan kimia air. Penentuan baku mutu air mengacu kepada Pasal 1 PP No. 82 tahun 2001. Rencana Restorasi Sistem Tata Air Rawa Kajian ini menguraikan rencana restorasi sistem tata air di areal hutan gambut rawa Tripa (TPSF) melalui rekayasa aliran sungai untuk manajemen tata air, termasuk upaya rehabilitasi dan manajemen tata air melalui pengaturan drainase untuk mengatasi drainase yang berlebihan (over drainage). Monitoring dan Evaluasi Monitoring dan evaluasi TPSF ditujukan untuk memberikan informasi pentingnya pemantauan kondisi aktual hidrologi Rawa Tripa pasca kajian dan mengetahui realisasi rencana implementasi beberapa rekomendasi penting di masa yang akan datang untuk memperbaiki dan menjaga agar ekosistem TPSF menjadi lebih baik, lestari, berkelanjutan serta bermanfaat. D. Tinjauan Teoritis Untuk analisis hidrologi syarat data yang digunakan harus; konsisten dan homogen, independent, representatif, continue serta panjang data yang cukup. Sebelum digunakan, data harus disaring atau diperiksa, secara manual dan statistik, untuk melihat pemenuhannya terhadap syarat di atas. Untuk menetapkan curah hujan maksimum boleh jadi (HMB atau PMP) dan banjir maksimum boleh jadi (BMB atau PMF), agar hasil analisisnya akurat diperlukan data pengamatan jangka panjang lebih dari 30 tahun. Data hujan yang digunakan adalah hujan harian, hujan harian maksimum tahunan. Hujan harian maksimum tahunan digunakan untuk estimasi hujan rencana. Pengujian data dilakukan untuk menghindari penggunaan data tidak valid yang disebabkan oleh : karena salah pencatatan, data dengan nilai sama terjadi berurutan, data kosong dan data pencilan (out layer). Pemeriksaan dan Penyaringan Data Pengolahan data hujan untuk analisa hujan rancangan (design rainfall) dan analisa hidrologi memerlukan analisa data yang hilang (missing data), perbaikan data dengan analisis kepanggahan. Kemudian penyaringan data (screenning data) dilakukan untuk mengetahui mutu data yang dipakai dalam analisis, yang meliputi: (1) Penyaringan data secara manual, (2) Penyaringan data hujan secara statistik dan (3) Pengujian data hujan hasil penyaringan. Penyaringan data secara manual dilakukan dengan menyaring data hujan dari setiap stasiun hujan, yaitu dengan metode partial series dengan menetapkan suatu batas bawah tertentu (threshold) dengan pertimbangan-pertimbangan tertentu, selanjutnya semua besaran hujan yang lebih besar daripada batas bawah tersebut diambil dan dijadikan bagian seri data untuk kemudian dianalisis dengan metode-metode yang telah ditentukan.


318


Penyaringan data hujan secara statistik yang dimaksud adalah analisa sebagai berikut: (1) Analisis ketidakadaan trend pada data hujan (homogenitas), (2) Analisis kestabilan data (stationary), (3) Analisis ketidakadaan persistensi data, data bersifat acak (randomness) atau tingkat independensi data dan (4) Analisis data hujan yang jarang terjadi (outlier). Pengujian data hujan hasil penyaringan dilakukan untuk melihat tingkat kevalidan data. Metode pengujian yang digunakan adalah uji Wald-Wolflowtz, Mann Whitney, dan Grubbs & Beck dengan derajat kepercayaan 95%. Prosedur penyaringan data hujan harian maksimum tahunan dapat dilihat pada Gambar 1. Periksa simbol pencatatan data

Periksa panjang pencatatan data

Tidak lolos

Tidak digunakan

Tidak lolos

Tidak digunakan

Lolos dan meragukan Periksa HHMT < 20 mm

Pengumpulan data Hujan Harian Maksimum Tahunan (HHMT)

Lolos dan meragukan Periksa HHMT  H

bulannya Lolos dan meragukan Periksa HHMT >400 mm H

bulannya

Tidak lolos Tidak lolos

Lolos dan meragukan Periksa HHMT > 400 mmH1..H2 Lolos dan meragukan

START

Pengujian statustik

Tidak lolos Tidak lolos

Tidak digunakan Tidak digunakan Tidak digunakan

Tidak digunakan

Lolos dan meragukan

Pengujian spasial

Tidak lolos

Lolos dan meragukan

Pengujian R max, R100  Sn END

Tidak lolos

Tidak digunakan Tidak digunakan

Pos yang lolos dan HHMT yang lolos

Gambar 1. Prosedur Penyaringan Data Hujan Harian Maksimum Tahunan Ketersediaan Data Debit Data yang diperlukan untuk analisis ketersediaan air adalah data debit harian atau bulanan dengan periode pencatatan yang cukup panjang, minimal 10 tahun. Data yang dipakai merupakan hasil rekaman pos duga air di lokasi kajian atau dekat sebelah hulu atau hilirnya. Bila data terlalu pendek atau tidak tersedia, debit bulanan dapat disimulasi berdasarkan data hujan dan data evapotranspirasi potensial pada daerah studi dengan bantuan model matematika hubungan hujan-limpasan. Tingkat keandalan ketersediaan air di dalam rawa disimulasikan berdasarkan konsep tampungan. Besarnya volume tampungan rawa netto yang terjadi ditentukan secara simulasi neraca air sebagai fungsi inflow dan outflow dalam interval waktu tertentu (misalnya bulanan atau tahunan).

LAPORAN UTAMA


| 319

Analisis Distribusi Hujan Distribusi hujan per jam ditetapkan dengan cara pengamatan langsung terhadap data pencatatan hujan setiap jam pada stasiun yang paling berpengaruh pada DAS. Distribusi tersebut diperoleh dengan pengelompokkan tinggi hujan ke dalam range dengan tinggi tertentu. Dari data yang telah disusun dalam range tinggi hujan tersebut dipilih distribusi tinggi hujan rancangan dengan berdasarkan hasil analisis frekuensi kemunculan tertinggi pada distribusi hujan jam-jaman tertentu. Selanjutnya prosentase hujan tiap jam diperoleh dengan menbandingkan tinggi pada setiap jam terhadap tinggi hujan total pada distribusi hujan yang ditetapkan. Namun demikian jika tersedia data hujan otomatis (Automatic Rainfall Recorder, ARR), maka pola distribusi hujan jam-jaman dapat dibuat dengan menggunakan metode kurva massa (mass curve) untuk tiap kejadian hujan lebat dengan mengabaikan waktu kejadian. Distribusi hujan per jam dengan interval tertentu perlu diketahui untuk menghitung hidrograf banjir rancangan dengan cara hidrograf satuan (unit hidrograf). Salah satu persamaan dari agihan hujan jam-jaman dirumuskan oleh Dr. Ishiguro yang dikenal dengan rumus Mononobe (Brook et al., 2003), yakni :

Ro 

R 24 t ............................................................................ 1) 2/3

t .................................................................... 2) Rt  Ro   T RT= t . Rt – (t - 1) R(t-1) ......................................................... 3)

Dimana, Rt = curah hujan rerata sampai jam ke-T (mm) Ro = curah hujan harian rerata (mm) T = waktu mulai hujan sampai ke-t (jam) R24 = curah hujan netto dalam 24 jam (mm) t = waktu konsentrasi (jam) RT = curah hujan sampai jam ke-T (mm) Waktu konsentrasi dapat dihitung dengan rumus Kinematik sebagai berikut:  L0.6  n 0.6  tc  0.93 0.4 0.3   i  S  ..................................................................... 4) Dimana, tc = waktu konsentrasi (jam) L = panjang over flow n = Manning overland roughness A = luas DAS S = kemiringan lahan Analisis Frekuensi Analisis frekuensi adalah analisis untuk memperkirakan harga besaran hidrologi (variate) yang masa ulangnya panjang, atau digunakan untuk peramalan dalam arti menentukan peluang terjadinya suatu peristiwa bagi tujuan perencanaan di masa datang. Variate terbesar yang didapatkan dari pengamatan hujan dan banjir, biasanya tidak ada sebesar atau lebih dari pada variate yang besarnya diperkirakan sebelumnya. Karena itu perlu Project Implementation Unit - Studi Ekosistem Rawa Tripa Universitas Syiah Kuala

320


dibuat suatu ekstrapolasi secara tepat, hanya mungkin jika persamaan matematis dari lengkungnya diketahui berdasarkan hasil perhitungan. Distribusi Gumbel. Distribusi ini menghasilkan estimasi paling besar diantara distribusi peluang yang lain. Persamaan umum estimasi banjir/hujan rencana periode T tahun (Sri Harto, 2000) adalah: X T  X  K G  S x ................................................................. 5)

KG 

YT  Yn 1 YT   ln(  ln(1  )) ......................... 6) dan Sn T

Persamaan di atas dapat disubstitusikan menjadi:

YT  Yn  Sx Sn ......................................................... 7) Yn 1 Sx jika: a  S dan b  X  a n XT  X 

Maka persamaan umum distribusi Gumbel menjadi : XT  b 

1 YT .................................................................... 8) a

Diamana, XT = curah hujan rancangan dengan periode ulang T tahun = curah hujan rerata X KG = faktor frekuensi distribusi Gumbel Sx = standar deviasi T = periode ulang T tahun YT = reduced variate untuk periode ulang T tahun Yn = reduced mean untuk jumlah n data Sn = reduced standart deviation untuk jumlah n data Distribusi Log Pearson Type III. Metode yang dianjurkan dalam pemakaian distribusi Log Pearson Type III adalah dengan mengkonversikan datanya menjadi bentuk logaritmis. Persamaan dari distribusi ini adalah(Sri Harto, 2000): log X T  log X  G  Slog x ..................................................... 9)

Dimana : XT = curah hujan rancangan dengan periode ulang T tahun log X = rerata dari logaritma curah hujan G = faktor frekuensi yang tergantung pada nilai koefisien Skewness Cs Slog X = standart deviasi dari logaritma curah hujan Koefisien skewness (kepencengan) dari logaritma curah hujan dihitung dari : n(log X i  log X )3 Cs  3 (n  1)(n  2)S logX ................................................. 10)

LAPORAN UTAMA


| 321

Distribusi Normal. Persamaan umum distribusi Normal (Sri Harto, 2000) adalah:

X T  X  K N  S x ................................................................ 11) Koefisien kekerapan normal KNdirumuskan seperti berikut:

KN  w 

2,515517  0,802853w  0,010328w 2 1  1,432788w  0,189269w 2  0,001308w 3 ......... 12)

Variabel w untuk kemungkinan terlampaui p = 1/T, dapat dihitung dengan:  1 w  ln  2 p

  untuk p < 0,5 dan nilai KN bertanda positif ................ 13)  

 1   w  ln  2  untuk p > 0,5 dan nilai KN bertanda negatif ........ 14)  (1  p) 

Hujan dengan Periode Ulang Curah hujan rancangan adalah curah hujan terbesar yang mungkin terjadi dalam suatu daerah dengan kala ulang tertentu, yang dipakai sebagai dasar perhitungan perencanaan dimensi suatu bangunan. Analisis probabilitas dilakukan untuk memperoleh curah hujan rancangan dengan periode ulang 2, 5, 10, 25, 50, 100, 200, 500 dan 1000 tahunan. Hujan rancangan ini diestimasikan dari distribusi frekuensi yang terpilih. Pemilihan kala ulang ditentukan oleh pertimbangan-pertimbangan hidro ekonomis yang didasarkan terutama pada: (1) Besarnya kerugian yang diderita apabila bangunan-bangunan rusak akibat banjir dan sering tidaknya perusakan itu terjadi, (2) Umur ekonomis bangunan, (3) Biaya pembangunan. Kriteria pemilihan kala ulang curah hujan rancangan, adalah: (1) Faktor-faktor penetuan curah hujan rancangan, (2) Pentingnya bangunan dan derajat resiko yang masih dipikul, (3) Ekonomi dengan memperhatikan usia guna. Dalam menentukan curah hujan rancangan dikategorikan sebagai berikut: (a) Bila kehancuran dari bangunan tersebut mengakibatkan penderitaan penduduk dengan korban jiwa, maka inflow hujan rancangan dibuat semaksimal mungkin, (b) Bila kehancuran dari bangunan tersebut hanya menimbulkan sedikit kerusakan pada penduduk tanpa menimbulkan korban jiwa, maka inflow hujan rancangan bisa berkurang dari ketentuan. Uji Kecocokan Distribusi Frekuensi Dalam pemilihan distribusi yang sesuai biasanya didasarkan pada uji kecocokan, uji kecocokan yang sering digunakan adalah Chi-Square dan Smirnov Kolmogorov. Dalam studi ini ditambah dengan uji Least Square dan Standard Error. Semua uji menggunakan derajat kepercayaan 95%. Uji Chi-Square. Metode ini menganggap pengamatan membentuk variable acak dan dilakukan secara statistik dengan mengikuti asimtot distribusi chi square dengan derajat kebebasan k-p-1, p adalah jumlah parameter yang diesitimasi dari data. Uji statistik ini berdasarkan pada bobot jumlah kuadrat perbedaan antara pengamatan dan teoritisnya


322


yang dibagi dalam kelompok kelas. Persamaan uji kecocokan ini (Sri Harto, 2000)sebagai berikut. k

h2   i 1

Dimana: h k Oi Ei

Oi  Ei 2 Ei

.......................................................... 15)

= parameter chi-kuadrat terhitung = jumlah sub kelompok = jumlah nilai pengamatan pada sub kelompok ke I = jumlah nilai teoritis pada sub kelompok ke i

Jika hasilnya besar menunjukkan bahwa distribusi yang dipilih tidak cocok, uji ini akan memberikan hasil yang baik jika mempunyai panjang pencatatan (n) yang besar. Kottegoda (1996) menyarankan sebaiknya n  50 tahun dan jumlah kelas interval  5. Uji Smirnov – Kolmogorov. Uji kecocokan ini adalah uji kecocokan non parametrik karena tidak mengikuti distribusi tertentu. Uji ini menghitung besarnya jarak maximum secara vertical antara pengamatan dan teotitisnya dari distribusi sampelnya. Perbedaan jarak maksimum untuk Smirnov – Kolmogorov (Sri Harto, 2000) tertera pada persamaan berikut.

Dn  max P(x)  P0 (x) ...................................................... 16) Dimana : Dn = Jarak vertikal maximum antara pengamatan dan teoritisnya P(x) = Probability dari sample data P0(x) = Probability dari teoritisnya Distribusi dikatakan cocok jika nilai Dn lebih kecil dari D kritisnya pada derajat kepercayaan yang diinginkan. Uji Least Square. Metode kuadrat terkecil mencocokkan distribusi yang berbeda dengan sample data dengan cara menghitung jumlah kuadrad terkecil nilai variate x hasil perhitungan dengan nilai pengamatan. Kesulitan yang dihadapi pada saat interpolasi hasil perhitungan yang mempunyai periode ulang sama dengan periode ulang pada variate x hasil pengamatan. Jenis distribusi yang terpilih adalah distribusi yang mempunyai nilai terkecil dan diformulasikan (Sri Harto, 2000) sebagai berikut.  X  YiT  LS   iT   n  m j 

1/ 2

............................................................ 17)

Dimana : LS = nilai Least Square XiT = nilai pengamatan pada kejadian i periode ulang T YiT = nilai perhitungan (teoritisnya) pada kejadian i periode ulang T n = jumlah kejadian mj = jumlah parameter yang diestimasi

LAPORAN UTAMA


| 323

Probable Maximum Precipitation (PMP) Hujan maksimum boleh jadi atau probable maximum precipitation (PMP) diartikan sebagai besaran hujan terbesar dengan durasi tertentu yang secara hidrologis dan meteorologis dimungkinkan bagi suatu daerah pengaliran dalam suatu waktu dalam tahun, tanpa adanya kelonggaran yang dibuat untuk trend klimatologis jangka panjang. Rumus yang digunakan dan kemudian dimodifikasi oleh Herschfield dengan menggunakan data curah hujan maksimum harian, didasarkan atas persamaan frekuensi umum (Subramanya, 2009), yaitu: X m  X n  K m  Sn ............................................................. 18) Dimana : Xm

Xn Km Sn

= curah hujan rancangan dengan periode ulang T = rerata (mean) curah hujan yang dikoreksi = faktor frekuensi = standar deviasi data yang dikoreksi

Langkah-langkah perhitungan PMP adalah sebagai berikut: (1) Perhitungan parameter statistik. Rerata X n dan standar deviasi Sn hujan maksimum harian dapat dihitung dari data, kemudian dihitung X nm dan Sn-m dengan tidak menyertakan data curah hujan maksimum terbesar; (2) Penyesuaian X n dan Sn untuk Maximum Observed Event. Faktor penyesuaian X n (fx1) dan Sn(fs1) untuk maximum observed rainfall didapat dari kurva penyesuaian Hershfield, (3) Penyesuaian X n dan Sn untuk ukuran sampel. Faktor penyesuaian X n (fx1) dan Sn(fs1) untuk panjang pengamatan didapat dari kurva penyesuaian Hershfield (Gambar 2); (4). Faktor frekuensi KSm . Faktor frekuensi Km didapatkan dari kurva empiris (Gambar 3); (5) Penyesuaian untuk waktu pengamatan yang tertentuPengamatan hujan pada umumnya dilakukan setiap hari. Karena data hujan dihitung berdasarkan suatu waktu pengamatan yang tertentu, nilai PMP yang diperoleh harus dikalikan dengan suatu faktor penyesuaian fo; (6) Penyesuaian untuk luas daerah. Konversi dari hujan pengamatan di suatu titik ke hujan rata-rata DAS dilakukan dengan faktor reduksiluas ; (7) Perhitungan PMP (Ayob dan Wahab ,2003).Penyesuaian terhadap rerata curah hujan maksimum harian adalah: '

X n  fx1  fx 2  X n ............................................................ 19) Penyesuaian terhadap standar deviasi curah hujan maksimum harian adalah:

Sn'  fx1  fx 2  Sn .............................................................. 20) PMP yang belum disesuaikan adalah: '

X m  X n  K m  Sn' ............................................................. 21) PMP yang telah disesuaikan adalah:

PMP  fo  X m .................................................................... 22) PMP titik disesuaikan untuk luas DAS adalah: PMP'    PMP .................................................................. 23) Hidrograf Satuan Untuk memperkirakan besarnya debit banjir yang terjadi pada suatu aliran sungai sangat dipengaruhi oleh ketersediaan data, apakah data debit (AWLR) atau data hujan. Jika dalam pelaksanaan pekerjaan nantinya data yang tersedia hanya berupa data hujan,


324


untuk memperkirakan debit banjir akan menggunakan metode empiris, metode unit hidrograf.

Gambar 2. Faktor Penyesuaian rerata (a) dan standar deviasi (b) terhadap Pengamatan Maksimum

Gambar 3. Hubungan Km dan rerata Hujan Maksimum Tahunan Dimana : Cp Qp tp t A tL

= koefisien debit puncak = ordinat maksimum dari unit hidrograf = waktu dimana Qp terjadi (jam) = durasi dari hujan lebih = luas Daerah Aliran Sungai (km2) = waktu tenggang = Tp/5,5 untuk unit hidrograf standard.

Hidrograf Satuan Sintetik Nakayasu. Hidrograf dapat diartikan sebagai penyajian grafis antara salah satu unsur aliran dengan waktu. Hidrograf ini menunjukkan tanggapan menyeluruh (integral responce) DAS terhadap masukan tertentu sesuai dengan sifat dan prilaku DAS yang bersangkutan. Besarnya debit puncak hidrograf akibat hujan satuan (Sri Harto, 2000) adalah:

LAPORAN UTAMA


Qp 

| 325

1,2 A  Ro 3,6  ( 0,3 Tp  T0,3 ) .................................................. 26)

Tp  Tg  0,8 Tr .................................................................. 27) T0,3 adalah waktu yang diperlukan untuk penurunan debit, dari debit puncak sampai menjadi 30 % dari debit puncak.

T0,3    Tg Tg  0,2  L0,7

........................................................................ 28) untuk L  15 km

Tg  0,4  0,058  L

untuk L  15 km

Bagian lengkung naik (rising limb) hidrograf mempunyai persamaan: t  Qa  Qp    Tp 

2,4

............................................................................. 29)

0 < t < Tp Bagian lengkung turun (decreasing limb) hidrograf dirumuskan : Tp< t < (Tp + T0,3) : Q d  Q p  0.30 (Tp + T0,3) < t < (Tp + T0,3 + 1,5T0,3) : t  (Tp + T0,3 + 1,5T0,3) :

t Tp T0.3

............................ 30)

Q d  Q p  0.30

Q d  Q p  0.30

t Tp 1.50T0.3 2T0.3

t Tp  0.50T0.3 1.5T0.3

........... 31)

........ 32)

Dimana : Qp= debit puncak banjir, m3/det Ro= besarnya hujan satuan, Ro = 1 mm Tp= tenggang waktu dari permulaan hujan sampai puncak banjir, jam Tg= waktu konsentrasi, jam Tr= waktu lamanya hujan, Tr = 1 jam L= panjang alur sungai utama, km A= luas daerah aliran sungai, km2 = koefisien sesuai kondisi daerah aliran sungai,  = 1,5 ~ 3,5 Qa=debit limpasan sebelum mencapai debit puncak, m3/dt Qd=debit limpasan setelah melewati debit puncak, m3/dt Hidrograf Satuan Sintetik Gama I. Beberapa karakteristik parameter DAS yang diperlukan adalah: (1) Luas Daerah Aliran Sungai (A), (2) Panjang sungai utama (L), (3) Selisih ketinggian hulu dan hilir sungai (H), (4) Kemiringan sungai (S), (5) Kerapatan jaringan kuras yaitu jumlah panjang semua sungai dan anak sungainya per satuan luas DAS (D), (6) Faktor sumber yaitu jumlah panjang semua sungai tingkat satu dibandingkan dengan panjang sungai semua tingkat (SF), (7) Frekuensi sumber yaitu jumlah pangsa sungai tingkat satu dibandingkan dengan jumlah seluruh pangsa sungai (SN), (8) Faktor lebar yaitu perbandingan antara lebar DAS yang diukur di titik sungai yang berjarak 0,75 panjang sungai dari stasiun hidrometri, dengan lebar DAS yang diukur di titik di sungai yang berjarak 0,25 panjang sungai (WF), (9) Luas DAS bagian hulu yaitu luas relatif DAS Project Implementation Unit - Studi Ekosistem Rawa Tripa Universitas Syiah Kuala

326


di bagian hulu garis yang ditarik melalui titik di sungai yang terdekat dengan titik berat DAS, tegak lurus garis hubung titik tersebut dengan stasiun hidrometri (RUA), (10) Faktor ` yaitu waktu yang diukur dari saat hidrograf mulai naik sampai saat terjadinya debit puncak, dinyatakan dengan persamaan (Sri Harto, 2000): 3

 L  TR  0,43    1,0665 SIM  1,2775 ...................... 33) 100 SF 

Debit puncak hidrograf dinyatakan dengan persamaan : Q p  0,1836 A 0,5886  JN 0,2381  TR 0,4008 ......................... 34)

Waktu dasar hidrograf yaitu waktu yang diukur dari saat hidrograf mulai naik sampai berakhirnya limpasan langsung, atau debit sama dengan nol, dinyatakan dengan persamaan:

TB  27,4132 TR 0,1457  S 0,0986  SN 0,7344  RUA 0,2574 ......... 35) Koefisien tampungan yang menunjukkan kemampuan DAS dalam fungsinya menampung air, dinyatakan dengan persamaan:

K  0,5617 A 0,1798  S 0,1446  SF 1,0897  D 0,0452 .................. 36) III. HASIL DAN PEMBAHASAN A. Karakteristik DAS Wilayah sungai (WS) Tripa – Batee memiliki 3 sungai utama, yaitu sungai Krueng Tripa, sungai Krueng Batee dan sungai Krueng Seumayam. WS Tripa – Batee mempunyai hulu di koridor tengah Aceh, meliputi Kabupaten Aceh Tengah dan Gayo Lues. Profil memotong WS Tripa- Batee disajikan pada Gambar 4. Morfologi wilayah ini berupa dataran landai di bagian Selatan, perbukitan bergelombang di bagian tengah, dan perbukitan terjal di bagian Utara (Gambar 4). Batuan-batuan lunak hasil sedimentasi material sungai dan pantai membentuk dataran landai di bagian Selatan, sedangkan batuan-batuan keras yang merupakan kerangka Pulau Sumatera terdiri dari batuan beku dan batuan metamorf membentuk perbukitan terjal di bagian Utara. Kawasan Rawa Tripa yang menjadi wilayah studi terletak di hilir WS. Secara litologi, wilayah ini terbentuk dari sedimentasi sungai dan pembusukan bahan organik. Proses yang berlangsung dalam waktu yang lama, mengakibatkan terjadinya rawa gambut yang luas. Berdasarkan data topografi parameter hidrologis DAS Krueng Seumayam sebagai berikut : - Catchment area (pada site AWLR) = 110,60 km2 - Elevasi hulu sungai = + 1.660,00 - Elevasi hilir sungai (pada site bendung) = + 29,00 - Slope rata-rata sungai = 0,000041 - Panjang sungai utama = 64.17 km Kemudian ditinjau dari metode kuantitatif lain dalam jaringan sungai suatu DAS adalah penentuan kepadatan aliran (drainage density) yang dinyatakan dalam rumus (Asdak, 1995).

LAPORAN UTAMA


Dd 

| 327

L ……………………………………………………………………….37) A

Dimana : Dd = Kepadatan aliran (km/km2) L = Panjang sungai total (km) A = Luas DAS (km2) Jika nilai kepadatan aliran lebih kecil dari 1 mile/mile2 (0,62 km/km2), DAS akan mengalami penggenangan, sedangkan jika nilai kepadatan aliran lebih besar dari 5 mile/mile2 (3,10 km/km2), DAS sering mengalami kekeringan. Dalam artian lain semakin besar angka kerapatan maka makin memperpendek waktu konsentrasi, sehingga memperbesar laju aliran permukaan. DAS Krueng Seumayam mempunyai kerapatan drainase (Dd) sebesar 1,32 km/km2 sehingga dapat dikategorikan sebagai DAS yang mempunyai kemungkinan penggenangan sedang. Secara umum, penggenangan di DAS Seumayam terjadi pada areal-areal tertentu yang mempunyai cekungan dan di hilir. Penggenangan di hilir ini kemudian disebut sebagai Kawasan Rawa Tripa.

Gambar 4. Profil Memotong WS Tripa- Batee B. Analisa Hujan Areal Data curah hujan yang dipakai untuk perhitungan debit banjir adalah hujan yang terjadi pada daerah aliran sungai pada waktu yang sama. Curah hujan yang diperlukan untuk penyusunan suatu rancangan pemanfaatan air dan rancangan pengendalian banjir adalah curah hujan rata-rata di seluruh daerah yang bersangkutan, bukan curah hujan pada suatu titik tertentu. Curah hujan ini disebut curah hujan area dan dinyatakan dalam mm (Sri Harto, 2000). Pada studi ini, metode yang digunakan untuk menentukan hujan areal adalah dengan metode polygon Thiessen. Metode perhitungan berdasarkan rata-rata timbang (weighted average). Metode ini memberikan proporsi luasan daerah pengaruh stasiun hujan untuk mengakomodasi ketidakseragaman jarak. Daerah pengaruh dibentuk dengan menggambarkan garis-garis sumbu tegak lurus terhadap garis penghubung antara dua stasiun hujan terdekat. Metode ini didasarkan pada asumsi bahwa variasi hujan antara stasiun hujan yang satu dengan lainnya adalah linear dan stasiun hujannya


328


dianggap dapat mewakili kawasan terdekat. Peta pengaruh hujan areal disajikan pada Gambar 5. Pada WS Tripa – Batee, terdapat beberapa stasiun hujan. Stasiun tersebut memberikan pengaruh yang besar terhadap kondisi hujan aliran di WS Tripa – Batee. Diantara stasiun tersebut, terdapat 5 stasiun yang member informasi mengenai kondisi hujan-aliran di WS tripa Batee. Lokasi dari masing-masing stasiun hujan disajikan pada Lampiran Peta Hujan Areal. Berdasarkan hasil analisa, distribusi hujan areal yang terjadi di WS Tripa – Batee disajikan pada Tabel 4. Berdasarkan Tabel 4 terlihat bahwa hujan – aliran di WS Tripa – Batee sangat dipengaruhi oleh hujan yang terjadi di hulu WS. Pada kondisi tertentu, wilayah hilir sungai-sungai di WS ini mengalami banjir walau tidak terjadi hujan di hilir. Akan tetapi, bila hujan terjadi di hulu dan tengah, maka dapat mengakibatkan banjir di sungai. Tabel 4. Hujan Areal WS Tripa – Batee dan Rawa Tripa Pengaruh Terhadap WS

Pengaruh Terhadap Kawasan Rawa Tripa

No

Nama Stasiun

1

Sta BMKG Cut Nyak Dhien

2,82%

0,00%

2

Sta Socfindo

46,86%

100,00%

3

BPP Beutong

2,95%

0,00%

4

101b (Sta Ceumpa)

41,17%

0,00%

5

97b (BPP Kota Baro)

6,21%

0,00%

Total

100,00%

100,00%

Gambar 5. Peta Pengaruh Hujan Areal di Wilayah Sungai Tripa – Batee

LAPORAN UTAMA


| 329

C. Hujan – Aliran Konsep Daerah Aliran Sungai (DAS) merupakan dasar dari semua perencanaan hidrologi. Mengingat DAS yang besar pada dasarnya tersusun dari DAS-DAS kecil, dan DAS kecil ini juga tersusun dari DAS-DAS yang lebih kecil lagi. Secara umum DAS dapat didefinisikan sebagai suatu wilayah yang dibatasi oleh batas alam seperti punggung bukit-bukit atau gunung, maupun batas buatan seperti jalan atau tanggul dimana air hujan yang turun di wilayah tersebut terdistribusi berupa aliran pengeluarannya. Batas aliran tidak ditetapkan berdasarkan air bawah tanah karena permukaan air tanah selalu berubah sesuai dengan musim dan tingkat kegiatan pemakaian. Aliran permukaan yang terjadi di kawasan studi (Rawa Tripa) secara alamiah terjadi akibat hujan yang turun di daerah tangkapan. Pada studi ini, aliran permukaan yang terjadi dianalisa berdasarkan kejadian hujan yang menyebabkan aliran air di dalam sungai. Berdasarkan kondisi hujan areal dan debit hasil pengukuran AWLR di Krueng Seumayam, hubungan antara hujan aliran di Rawa Tripa disajikan pada Gambar 6. Banjir di kawasan Rawa Tripa terjadi pada bulan Februari, Maret, Agustus, Oktober, November dan Desember. Analisa Banjir Kajian banjir yang terjadi di kawasan Rawa Tripa bersumber dari banjir akibat limpasan sungai dan banjir yang terjadi akibat hujan setempat. Banjir yang terjadi akibat limpasan sungai di dominasi oleh kejadian banjir dari sungai Krueng Seumayam dan sungai-sungai kecil yang melintasi kawasan Rawa Tripa. Beberapa sungai kecil yang memberikan konstribusi dominan terhadap kejadian banjir di Rawa Tripa seperti Alue Ie Mirah, Suak Palembang, Alue Beringen, Gunoeng Samarinda. Berdasarkan kondisi ini, maka stasiun utama yang dijadikan rujukan untuk menganalisa banjir di Rawa Tripa adalah stasiun hujan PT. Socfindo. Pola Hujan Badai di Aceh Penentuan pola hujan badai diperlukan untuk melihat trend kejadian hujan badai di suatu kawasan. Hal ini disebabkan oleh terbatasnya ketersediaan data dalam waktu yang cukup panjang di Provinsi Aceh. Pola hujan badai untuk kawasan Selat Malaka pada dasarnya berbeda dengan kondisi badai di pesisir Sumatera yang berada di sisi Samudera Hindia. Curah hujan badai rancangan untuk Selat Malaka dan koridor Tengah Aceh disajikan pada Tabel 5. Selanjutnya, pola hujan badai rancangan di kawasan pesisir Barat Selatan Sumatera yang dipengaruhi oleh Samudera Hindia disajikan pada Tabel 6. Tabel 5. Pola Perulangan Hujan Rencana berbagai Kala Ulang (mm) di Selat Malaka No. Pos

Nama Pos

Periode Ulang (Tahun) 2

5

10

20

25

50

100

200

500

1000

94

131

157

182

191

216

242

267

302

328

115

Bireuen Lhokseumaw e

103

142

176

212

225

264

305

348

407

453

116b

Telong

70

91

106

121

126

142

158

174

195

212

116c

Lampahan

63

84

96

106

109

118

127

135

144

152

116d

Baleg

74

95

110

123

127

140

152

165

181

193

117

Lhoksukon

110

155

187

218

228

259

290

322

363

395

117a

Panton Labu

102

130

149

166

171

187

203

218

238

253

114

Sumber : Dinas Sumber Daya Air NAD, 2003


330


Gambar 6. Hubungan Hujan – aliran di Kawasan Rawa Tripa Tabel 6. Desain Pola Hujan Badai (%) Pesisir Barat Selatan Sumatera Waktu (jam ke) No

Durasi (jam)

1 2 3 4 5 6 7 8 9 10 11

3 4 5 6 7 8 9 10 11 12 13

persen (%) 1 2 35 47 28 40 18 40 15 35 16 22 10 17 14 16 8 12 8 10 6 9 6 14

3 18 26 27 25 32 23 20 18 12 11 14

4

5

6

7

8

9

10

11

12

6 13 15 16 25 16 16 14 12 15

2 6 8 11 13 13 18 14 17

4 4 8 9 11 12 16 22

2 4 6 9 11 10 9

2 4 6 7 7 5

2 4 4 5 3

3 2 4 2

2 3 1.5

3 0.5

Hujan rancangan untuk kawasan Rawa Tripa diasumsikan selama 6 jam. Asumsi ini dilakukan berdasarkan kondisi hujan yang terjadi pada saat dilakukan survai lapangan. Distribusi hujan dianalisa dengan menggunakan metode yang dikembangkan oleh

LAPORAN UTAMA


| 331

Mononobe. Pola distribusi hujan jam-jaman di kawasan Rawa Tripa disajikan pada Tabel 7 dan Gambar 6. Berdasarkan pola distribusi tersebut, hujan di kawasan rawa Tripa mencapai puncak pada jam hujan ke – 2. Tabel 7. Distribusi Hujan di Rawa Tripa JAM

DISTRIBUSI

1

15,00%

2

35,00%

3

25,00%

4

15,00%

5

6,00%

6

4,00%

Jumlah

100,00%

1. 2. 3. 4. 5. 6.

Gambar 6. Pola Distribusi Hujan di Rawa Tripa Analisis Curah Hujan Rencana Curah hujan rancangan diperlukan sebagai data masukan pada analisis debit banjir rancangan. Untuk itu perlu dilakukan analisis curah hujan rancangan. Metode yang digunakan untuk melakukan analisis curah hujan rancangan dengan periode kala ulang tertentu adalah sebagai berikut : (1) Distribusi Gumbel Tipe I, (2) Distribusi Log – Normal Dua Parameter, (3) Distribusi Log – Pearson Tipe III. Hasil analisa curah hujan rencana dengan menggunakan ke tiga metode tersebut (Gambar 7) menunjukkan bahwa pola distribusi hujan di kawasan Rawa Tripa memiliki trend yang serupa. Akan tetapi, perbedaan akan semakin besar pada pengujian pola distribusi untuk kala ulang diatas 10 tahun atau lebih lama. Untuk itu perlu dilakukan uji kecocokan distribusi. Hujan Probalitas Maksimum (PMP) Hujan maksimum boleh jadi atau probable maximum precipitation (PMP) diartikan sebagai besaran hujan terbesar dengan durasi tertentu yang secara hidrologis dan dimungkinkan bagi suatu daerah pengaliran dalam suatu waktu dalam tahun, tanpa adanya kelonggaran yang dibuat untuk trend klimatologis jangka panjang. Rumus yang


332


digunakan dan kemudian dimodifikasi oleh Herschfield. Perhitungan PMP untuk kawasan Rawa Tripa disajikan pada Tabel 8. Tabel 8. Perhitungan PMP untuk Kawasan Rawa Tripa No. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

Parameter Sx Sx-m Xn Xn-m (Xn-m)/Xn (Sx-m)/Sx Grafik Km Faktor Adjusment: Xn Sn Xn Terkoreksi Sn Terkoreksi Hujan Terpusat Faktor Reduksi Luas DAS 2 110,60 km Fixed Time Internal Hujan PMP Terkoreksi Hujan PMP diambil dengan koefisien Sehingga hujan PMP adalah

= = = = = = = = = = = = = = = = =

Nilai 28,339 28,339 125,125 126,826 1,014 1,000 11,57

Unit mm mm

108,00 119,40 135,132 33,836 526,534

% % mm mm mm

85,36 101,00 445,020

% % mm

100 445,020

% mm

Keterangan

GRAFIK GRAFIK GRAFIK

GRAFIK

Gambar 7. Grafik Distribusi Curah Hujan di Rawa Tripa Berdasarkan pengujian distribusi untuk metode sebaran data yang ada dan PMP, maka dilakukan pendugaan curah hujan rancangan. Pengujian kecocokan distribusi sebaran dilakukan dengan menggunakan metode Smirnov – Kormonogorov. Pengujian uji kecocokan distribusi disajikan pada Tabel 9. Berdasarkan Tabel 9, distribusi yang digunakan pada analisa hujan di Rawa Tripa adalah distribusi dengan metode Log Pearson Tipe III . Hasil perhitungan curah hujan rancangan untuk kawasan Rawa Tripa disajikan pada Tabel 10.

LAPORAN UTAMA


| 333

Tabel 9. Pengujian kecocokan distribusi dengan metode Smirnov-Kosmonogorov pada hujan rencana Rawa Tripa

Data

Prob Empiris

Prob.Teoritis (%)

EJ. Log Log Pearson Gumbel Pearson Normal (1) (2) (3) (4) (5) (6) 195 3,45 3,6367 2,631 1,982 2,673 175 6,90 7,8823 6,385 6,090 5,962 173 10,34 8,5067 6,937 6,729 6,560 165 13,79 11,5055 9,210 9,285 9,024 155 17,24 16,6372 14,223 14,926 14,035 150 20,69 19,9106 17,228 18,098 17,057 148 24,14 21,3701 18,459 19,368 18,295 135 27,59 33,1890 32,005 34,153 31,873 133 31,03 35,3842 34,390 36,511 34,286 132 34,48 36,5197 35,596 37,690 35,506 131 37,93 37,6800 36,811 38,869 36,735 129 41,38 40,0741 39,269 41,227 39,223 127 44,83 42,5629 41,766 43,586 41,749 121 48,28 50,5400 49,500 50,854 49,573 116 51,72 57,6408 57,263 58,477 57,941 111 55,17 64,9389 65,454 66,100 66,786 107 58,62 70,7389 72,277 72,198 74,154 107 62,07 70,7389 72,277 72,198 74,154 106 65,52 72,1629 74,022 73,723 76,039 106 68,97 72,1629 74,022 73,723 76,039 106 72,41 72,1629 74,022 73,723 76,039 104 75,86 74,9609 77,564 76,772 79,862 104 79,31 74,9609 77,564 76,772 79,862 101 82,76 78,9930 82,173 81,077 82,635 99 86,21 81,5400 84,861 83,517 84,479 98 89,66 82,7637 86,226 84,737 85,414 97 93,10 83,9510 87,604 85,957 86,360 86 96,55 94,1237 96,947 95,879 94,209 Delta Maksimum Empiris dan Teoritis Delta Maksimum Kritis Empiris dan Teoritis dengan α 0,05 = (%)

Delta Probability EJ. Gumbel (7) 0,188 0,986 1,838 2,288 0,604 0,779 2,768 5,603 4,350 2,037 0,251 1,305 2,265 2,264 5,917 9,767 12,118 8,670 6,646 3,197 0,251 0,901 4,349 3,766 4,667 6,891 9,152 2,428

Log Pearson (8) 0,817 0,512 3,408 4,583 3,019 3,461 5,679 4,418 3,355 1,113 1,120 2,110 3,062 1,224 5,539 10,282 13,656 10,208 8,505 5,057 1,609 1,702 1,747 0,586 1,346 3,430 5,499 0,395

12,1182

13,6562

(9) 1,466 0,806 3,615 4,508 2,316 2,591 4,770 6,567 5,477 3,207 0,938 0,152 1,242 2,578 6,753 10,927 13,578 10,129 8,206 4,758 1,309 0,910 2,538 1,681 2,690 4,918 7,146 0,672

Log Normal (10) 0,775 0,935 3,785 4,769 3,207 3,632 5,843 4,286 3,251 1,023 1,196 2,157 3,079 1,297 6,217 11,614 15,533 12,085 10,521 7,073 3,625 4,000 0,552 0,124 1,728 4,241 6,744 2,342

13,5777

15,5330

Pearson

25,2044

Analisa Banjir Rencana Analisa banjir rencana dapat dilakukan dengan mengkaji pola hujan aliran atau menganalisa dengan menggunakan data AWLR yang tersedia. Pada kondisi kajian hidrologi di Rawa Tripa, ketersediaan data AWLR sangat terbatas. Data yang tersedia adalah data AWLR Krueng Seumayam dan AWLR Krueng Tripa. AWLR Krueng Seumayam terletak di Desa Alue Ie Mirah, tepatnya pada 3:58’12.9” dan 96:39’06.9”. AWLR ini mempunyai rentang data sejak bulan Januari 2013. Sedangkan AWLR lainnya terletak di sungai Krueng Tripa yaitu di kawasan Gunong Khong. AWLR ini mempunyai periode data terbatas. Oleh karena kondisi tersebut, kajian hidrologi di Rawa Tripa melakukan analisa empirik. Hasil perhitungan secara empiris ini akan dikalibrasi dengan hujan aliran. Metode yang dipakai untuk menganalisa banjir yang terjadi secara empiris adalah HSS Gama I dan metode Nakayasu (Sri Harto, 2000). Perhitungan banjir rencana dengan metode Gama I disajikan pada Tabel 11.


334


Berdasarkan hasil analisa, terlihat bahwa jam puncak banjir dengan analisa menggunakan metode GAMA-1 dalah 3,4658 jam. Hidrograf banjir di kawasan Rawa Tripa dengan menggunakan metode GAMA-1 disajikan pada Gambar 8.

Gambar 8. Hidrograf Satuan Sintetik metode GAMA – I untuk banjir di Kawasan Rawa Tripa

Tabel 10. Hujan Rencana Rawa Tripa METODE NO

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

KALA ULANG

1,1 2 5 10 15 20 25 50 75 100 150 200 500 1000 PMP

GUMBEL TYPE I

NORMAL

LOG NORMAL 3

LOG PEARSON 3

54,35 120,47 145,51 162,09 171,45 178,00 183,04 198,59 207,62 214,01 223,01 229,38 249,66 264,99

53,99 125,12 148,97 161,45 167,67 171,75 174,75 183,34 187,95 191,06 195,27 198,13 206,70 212,71

71,42 122,31 146,48 160,98 168,74 174,02 178,01 189,96 196,70 201,39 207,90 212,46 226,68 237,21

79,38 119,88 145,22 162,49 172,42 179,46 184,95 202,17 212,45 219,85 230,47 238,14 263,45 283,57 445,02

LAPORAN UTAMA


| 335

Tabel 11. Karakteristik DAS dan Parameter Banjir Metoda Satuan GAMA-1 No

Karakteristik DAS

1

Jumlah pangsa sungai tingkat 1 (N1)

=

25,00

buah

Sumber data Peta RBI

2

Jumlah pangsa sungai semua tingkat (Nt)

=

49,00

buah

Peta RBI

3

Panjang pangsa sungai tingkat 1(L1)

=

93,49

km

Peta RBI

4

Panjang pangsa sungai semua tingkat (Lt)

=

169,07

km

Peta RBI

5

Jumlah pertemuan sungai (JN)

=

22,00

buah

Peta RBI

6

Nilai

Luas DTA (A)

=

Unit

122,41

2

Peta RBI

2

km

7

Luas DTA hulu (AU)

=

101,02

km

Peta RBI

8

Panjang sungai utama (L)

=

57,04

km

Peta RBI

9

0.75 L

=

42,78

km

Peta RBI

10

0.25 L

=

14,26

km

Peta RBI

11

Kemiringan sungai rata-rata (S)

=

0,0067

12

Faktor sumber ( SF)

=

0,55

L1 / Lt

13

Frekuensi sumber (SN)

=

0,51

N1 / Nt

14

Kerapatan jaringan kuras (D)

=

1,38

15

B . 0.75 L (Wu)

=

16

B . 0.25 L (Wl)

17

Peta RBI

km/km 2

Lt / A

8,38

km

Peta RBI

=

4,30

km

Peta RBI

Faktor lebar (WF)

=

1,95

Wu / Wl

18

Perbandingan DTA hulu dan DTA (RUA)

=

0,83

Au / A

19

SIM = RUA * WF

=

1,61

Parameter banjir dengan metode GAMA – I: 20 21 22 23 24 25 26

3

TR = 0,43 (L / 100. SF) + 1,0665 SIM + 1,2775 QP = 0,1836 A

0,5886

TB = 27,4132 TR K = 0,5617 A

JN

0,1457

0.1798

0,2381

TR

-0,0956

S

-0.1446

S

-0.4008

SN

0,7344

-1.0897

SF

RUA

0,2574

0.0452

D

=

3,4658

=

3,945

=

31,259

=

5,323

=

10,433

=

14,270

–t / k

Qt = Qp . e -6 2 Phi indeks = 10,4903 – 3,859.10 . A + 1,6985 . -13 4 10 (A/SN) 0,6444 0,943 QB = 0,4751. A .D

Metode lain yang digunakan pada analisa banjir rencana adalah dengan menggunakan metode Nakayasu. Metode ini digunakan dengan anggapan bahwa metode ini dikembangkan di kawasan dengan kondisi yang serupa seperti di Indonesia (kawasan pulau). Perhitungan hidrograf satuan dengan menggunakan metode Nakayasu disajikan pada Tabel 12. Selanjutnya hidrograf satuan Nakayasu disajikan pada Gambar 9.


336


Tabel 12. Karakteristik DAS untuk Perhitungan hidrograf Satuan Nakayasu No. Karakteristik DAS 1 2 3 4 5 6 7 8

Panjang sungai ( Km ) Luas DAS ( Km2 ) tg = 0,4 + 0,058L tr = ( 0,5 – 1,0 ) tg Tp = tg + 0,8 tr T0,3 = 2tg Qp = A Ro/3,6/(0,3 Tp + T0,3) Konstanta Untuk Nilai T0,3

Nilai = = = = = = = =

57,04 122,41 3,71 1,85 5,19 7,42 3,79 2,00

Gambar 9. Hidrograf Satuan Sintetik Metode Nakayasu untuk Banjir di Kawasan Rawa Tripa Hasil analisa secara empirik untuk kejadian banjir tahunan di Kawasan Rawa Tripa kemudian diverifikasi berdasarkan kondisi survai. Berdasarkan hasil wawancara dengan masyarakat di desa Kuala Seumayam, bila banjir terjadi di Alue Ie Mirah (lokasi AWLR Krueng Seumayam), desa mereka akan mengalami banjir dalam 2,0 – 4,0 jam. Kondisi tersebut tergantung dari besarnya banjir yang terjadi. Selain itu, debit banjir yang terjadi mendekati debit banjir puncak (berdasarkan data grafik AWLR Krueng Seumayam).

LAPORAN UTAMA


| 337

Berdasarkan hasil perhitungan secara empiris dan kondisi lapangan, metode yang lebih mendekati dengan kondisi yang terjadi di lapangan adalah Metode Gama-1. Rekapitulasi debit banjir di kawasan Rawa Tripa berdasarkan metode Gama-1 disajikan pada Gambar 10 dan Tabel 13.

Gambar 10. Rekapitulasi hidrograf periode ulang banjir tahunan di Kawasan Rawa Tripa


338


Tabel 13. Rekapitulasi Perhitungan Hidrograf Banjir dengan Berbagai Kala ulang Waktu (jam) 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45

2 Tahun

5 Tahun

10 Tahun

15 Tahun

14,27 29,18 78,89 160,39 242,30 271,62 265,47 237,89 206,78 173,81 146,49 123,84 105,08 89,52 76,63 65,95 57,10 49,77 43,69 38,65 34,47 31,01 28,15 25,77 23,80 22,17 20,82 19,69 18,77 18,00 17,36 16,83 16,39 16,03 15,73 15,48 15,27 15,10 14,96 14,84 14,74 14,66 14,59 14,54 14,49 14,45

14,27 34,18 100,56 209,40 318,79 357,94 349,73 312,90 271,36 227,33 190,84 160,60 135,53 114,77 97,55 83,29 71,47 61,67 53,55 46,83 41,25 36,63 32,80 29,63 27,00 24,82 23,01 21,51 20,27 19,24 18,39 17,69 17,10 16,62 16,21 15,88 15,61 15,38 15,19 15,03 14,90 14,79 14,70 14,63 14,57 14,52

14,27 37,63 115,52 243,22 371,57 417,50 407,87 364,66 315,91 264,25 221,44 185,96 156,55 132,18 111,99 95,25 81,38 69,89 60,36 52,47 45,93 40,50 36,01 32,29 29,20 26,64 24,53 22,77 21,31 20,11 19,11 18,28 17,59 17,02 16,55 16,16 15,84 15,57 15,35 15,16 15,01 14,88 14,78 14,69 14,62 14,56

14,27 39,63 124,17 262,78 402,09 451,95 441,49 394,59 341,68 285,61 239,13 200,62 168,71 142,26 120,34 102,17 87,12 74,64 64,30 55,73 48,63 42,75 37,87 33,83 30,48 27,70 25,40 23,49 21,91 20,61 19,52 18,62 17,88 17,26 16,75 16,32 15,97 15,68 15,44 15,24 15,07 14,94 14,82 14,73 14,65 14,58

Kala Ulang 20 25 Tahun Tahun 14,27 41,05 130,32 276,70 423,81 476,47 465,42 415,89 360,02 300,80 251,73 211,06 177,36 149,42 126,28 107,09 91,20 78,02 67,10 58,05 50,55 44,34 39,19 34,92 31,39 28,45 26,02 24,01 22,34 20,96 19,81 18,87 18,08 17,43 16,89 16,44 16,07 15,76 15,50 15,29 15,12 14,97 14,85 14,75 14,67 14,60

14,27 42,16 135,12 287,55 440,76 495,59 484,09 432,50 374,32 312,66 261,55 219,20 184,10 155,01 130,91 110,93 94,38 80,66 69,29 59,86 52,06 45,58 40,22 35,78 32,09 29,04 26,51 24,41 22,68 21,24 20,04 19,06 18,24 17,56 16,99 16,53 16,14 15,82 15,55 15,33 15,15 15,00 14,88 14,77 14,69 14,62

50 tahun

75 Tahun

100 Tahun

14,27 45,65 150,26 321,80 494,20 555,90 542,96 484,91 419,44 350,05 292,54 244,88 205,38 172,65 145,53 123,05 104,42 88,98 76,18 65,58 56,79 49,51 43,47 38,47 34,33 30,89 28,04 25,69 23,73 22,11 20,77 19,65 18,73 17,97 17,33 16,81 16,38 16,01 15,72 15,47 15,26 15,09 14,95 14,84 14,74 14,66

14,27 47,75 159,34 342,32 526,22 592,03 578,23 516,31 446,47 372,45 311,10 260,27 218,13 183,22 154,28 130,30 110,43 93,96 80,31 69,00 59,63 51,86 45,42 40,09 35,66 32,00 28,96 26,45 24,36 22,63 21,20 20,01 19,03 18,21 17,54 16,98 16,52 16,13 15,81 15,55 15,33 15,15 15,00 14,87 14,77 14,68

14,27 49,26 165,89 357,13 549,34 618,13 603,70 538,99 465,99 388,63 324,51 271,38 227,34 190,85 160,61 135,54 114,77 97,56 83,30 71,47 61,68 53,56 46,83 41,25 36,63 32,80 29,63 27,00 24,82 23,01 21,51 20,27 19,25 18,39 17,69 17,10 16,62 16,21 15,88 15,61 15,38 15,19 15,03 14,90 14,79 14,70

Untuk memperhitungkan nilai ekstrim banjir yang terjadi di Rawa Tripa, dilakukan analisa kemungkinan banjir maksimum (PMF). Analisa secara empiris dilakukan dengan menggunakan metode Gama-1. Hasil analisa hidrograf banjir PMF dengan metode Gama-1 disajikan pada Gambar 11.

LAPORAN UTAMA


| 339

Banjir yang terjadi di Kawasan Rawa Tripa selain berasal dari sungai Krueng Seumayam, sungai Krueng Tripa dan Krueng Batee. Rekapitulasi kejadian banjir di sungai Krueng Semayam, Krueng Tripa dan Krueng Batee disajikan pada Tabel 14.

Gambar 11. Hidrograf Kemungkinan Banjir Maksimum (PMF) dengan Metode Gama-1 di Rawa Tripa Tabel 14. Rekapitulasi Debit Banjir Rawa Tripa dari Krueng Seumayam, Krueng Tripa dan Krueng Batee No.

Kala Ulang (Tahun)

Debit banjir Krueng 3 Seumayam (m /det)

Debit banjir Krueng Tripa 3 (m /det)

Debit banjir Krueng 3 Batee (m /det)

1 2 3 4 5 6 7 8 9 10

2 5 10 15 20 25 50 75 100 PMF

271,62 357,94 417,50 451,95 476,47 495,59 555,90 592,03 618,13 1.427,03

293,52 740,63 866,01 933,74 983,85 1.022,93 1.141,92 1.220,07 1.276,18 2.926,80

272,57 459,74 531,15 576,97 606,61 624,77 697,09 740,41 771,70 1.755,84

Untuk menentukan frekuensi kejadian banjir berdasarkan tampang sungai, dilakukan pengukuran lapangan tampang basah sungai. Pengukuran lapangan ini dilakukan untuk melihat kemampuan sungai terhadap debit banjir yang terjadi. Kondisi aliran sungai pada saat banjir dibahas pada analisa hidraulika sungai.


340


D. Analisa Hidraulika Sungai Kajian kondisi hidrometri WS Tripa – Batee, secara garis besar merupakan kajian terhadap hidrometri di sungai dengan sedimentasi/alluvium. Pengetahuan hidrolika fluvial ini di tekankan pada transportasi sedimen sebagai dasar dalam perencanaan bangunan pengendalian sungai. Faktor-faktor yang menentukan transportasi sedimen adalah sifat-sifat aliran air (flow characteristic); sifat-sifat sedimen (sediment characteristic); dan pengaruh timbal balik (interaction). Transportasi sedimen pada sungai sangat dipengaruhi oleh debit yang masuk ke sungai tersebut, oleh sebab itu diperlukan suatu pengukuran hidrometri yang meliputi pengukuran debit sesaat dan pengambilan sampel sedimen. Sungai adalah jalur aliran air di atas permukaan bumi yang selain mengalirkan air, juga mengangkut sedimen terkandung dalam air sungai tersebut. Besarnya material terangkut dari suatu sungai dapat ditentukan dengan pengukuran pengangkutan sedimen pada titik kontrol dan semi empiris yang ada. Beberapa faktor yang mempengaruhi laju proses sedimentasi di daerah pengaliran sungai diuraikan sebagai berikut: Cakupan Areal Daerah. Pengaliran.Kapasitas sedimen yang dihanyutkan oleh suatu sungai biasanya berbanding lurus dengan luas daerah pengaliran. Kondisi Geologi Daerah Pengaliran. Kondisi geologi sangat mempengaruhi intensitas proses-proses degradasi dan erosi pada batuan yang selanjutnya akan mempengaruhi intensitas sedimentasi pada sungai yang bersangkutan. Kondisi Topografi. Elevasi suatu daerah pengaliran, kondisi perbukitan maupun pegunungan, tingkat kemiringan akan sangat mempengaruhi intensitas degradasi dari batuan yang terdapat di daerah pengaliran. Kondisi Meteorologi. Karakteristik hujan yang jatuh pada suatu daerah pengaliran sangat mempengaruhi intensitas degredasi dan erosi dari pada batuan yang membentuk daerah pengaliran khususnya pada daerah pegunungan, dimana air hujan yang mengalir di atas permukaan tanah dengan mudah mengikis lapisan atasnya serta menghanyutkannya masuk ke dalam alur sungai. Karakteristik Hidraulika Sungai. Karakteristik yang dimaksud adalah debit sungai, kecepatan aliran, konfigurasi alur sungai, bentuk penamnpang lintang sungai, kemiringan dan kekasaran batuan pembentuk alur sungai. Vegetasi Daerah Pengaliran. Vegetasi di daerah pengaliran akan sangat membantu menurunkan intensitas proses degredasi maupun erosi yang terjadi. Kegiatan Manusia. Kegiatan manusia di atas daerah pengaliran seperti daerah perkampungan, pertanian, pembangunan jalan dan pengusahaan hutan juga akan mempengaruhi kapasitas sedimentasi pada suatu sungai. Angkutan material (sedimen) yang mengalir pada suatu sungai terdiri dari sedimen yang berasal dari muatan dasar (bed load) dan muatan melayang (suspended load). Muatan dasar (bed load), adalah sedimen dasar yang bergerak dalam aliran sungai dengan cara bergulir, meluncur dan meloncat-loncat di atas permukaan sungai. Muatan melayang (suspended load), adalah sedimen yang terdiri dari butiran halus dengan ukuran lebih kecil dari 0.1 mm dan senantiasa melayang. Besarnya volume sedimen dasar dan sedimen melayang yang terangkut dihitung dengan persamaan Jorgen Fredsoe dan Rolf Deigaard (Fredso and Deigaard, 1992).

LAPORAN UTAMA


| 341

Pengukuran Hidrometri Sesaat Pengukuran debit sesaat dilakukan untuk mendapatkan debit yang terjadi persatuan waktu sehingga diperoleh kecepatan air yang mengalirkan sedimen dari hulu ke hilir. Pengukuran debit sesaat dilakukan untuk mendapatkan angka koefisien yang nantinya dapat digunakan sebagai koefisien routing aliran dari stasiun hidrometri yang ada . Debit (discharge) atau besarnya aliran sungai (stream flow) adalah volume aliran yang mengalir melalui suatu penampang melintang sungai persatuan waktu. Biasanya dinyatakan dalam satuan meter kubik per detik (m3/det) atau liter per detik (l/det). Menurut Asdak (1995), salah satu teknik pengukuran debit sesaat di lapangan adalah dengan cara mengukur kecepatan aliran dan menentukan luas penampang melintang sungai. Pengukuran ini dilakukan dengan mengukur kecepatan air dan penampang melintang sungai yang diteliti. Pengukuran debit yang dilaksanakan di suatu pos duga air tujuannya adalah untuk membuat lengkung debit dari pos duga air yang bersangkutan. Lengkung debit dapat merupakan hubungan yang sederhana antara tinggi muka air dan debit, dapat pula merupakan hubungan yang komplek apabila debit di samping fungsi dari tinggi muka air juga merupakan fungsi dari kemiringan muka air, tingkat perubahan muka air dan fungsi dari faktor lainnya. Pada dasarnya pengukuran debit adalah pengukuran luas penampang basah, kecepatan aliran dan tinggi muka air. Keadaan aliran yang terdapat pada dasar sungai berbeda-beda, kadang-kadang ada kalanya sungai tersebut berarus cepat dan ada kalanya berarus lambat, keadaan aliran tersebut sangat ditentukan oleh besarnya kecepatan aliran yang terjadi pada sungai tersebut, dimana secara formulasi dapat dirumuskan sebagai bilangan Froude (Fr) sebagai berikut (Chow, 1959).

Fr 

V ……………………………………………………………… 38) gxh

keterangan : V = kecepatan aliran (m/dt) g = kecepatan gravitasi (m/dt2) h = kedalaman sungai (m) Bilangan Froude ini mempunyai batasannya, dimana jika batasan bilangan Froude < 1, maka aliran yang terjadi adalah sub kritis, jika yang terjadi bilangan Froude = 1, maka aliran yang terjadi dikatakan kritis, jika bilangan Froude > 1, maka aliran yang terjadi dikatakan superkritis, semakin kritis suatu aliran maka akan semakin besar dan kecepatan aliran (V) yang timbul begitu juga sebaliknya. Keadaan super kritis jarang terjadi hal ini terjadi pada saat-saat debit banjir atau sungai dalam keadaan tidak stabil dan aliran yang terjadi sangat deras, dalam keadaan normal biasanya aliran yang terjadi sub kritis atau kritis, dengan batasan bilangan Froude lebih kecil atau sama dengan 1. Hasil pengukuran parameter hidrolis sungai langsung dilapangan, diperoleh data: kedalaman sungai (h), lebar sungai (b) sehingga dapat diketahui kondisi aliran pada daerah tersebut. Hasil pengukuran tampang basah sungai dan debit normal serta hubungannya dengan debit banjir disajikan pada Tabel 15. Tabel 15 memberikan informasi bahwa seluruh sungai yang melintasi Rawa Tripa akan mengalami peristiwa melimpahnya air keluar dari sungai. Dengan kata lain, Rawa Tripa adalah derah yang selau digenangi banjir tahunan. Selain itu, pada kondisi normal, sungai ini mengalami proses pergerakan sedimen sungai yang dominan.


342


Hasil perhitungan parameter hidrolis sungai dapat dilihat pada Tabel 16 dan 17. Berdasarkan data tersebut, terlihat bahwa kondisi aliran di sungai utama dalam rawa Tripa berada pada kondisi sub kritis pada kondisi normal dan banjir. Tabel 15. Luas Tampang Basah, Debit Normal dan Debit Banjir Tahunan Normal

Lokasi

Asungai (m2)

Banjir Tahunan (m3/det)

Banjir 2 Tahunan (m3/det)

PMF (m3/det)

71,64

290,67

274,30

293,53

2926,80

1

Tripa

137,76

V (m/det ) 0,52

2

Tripa

144,60

0,55

79,53

267,80

274,30

3

Seumayam

77,10

0,43

33,15

90,87

256,05

4

Seumayam

70,70

0,46

32,53

89,30

256,05

271,62

1427,03

5

Seumayam

82,67

0,33

27,28

107,80

256,05

271,62

1427,03

6

Batee

98,80

0,58

57,30

286,70

264,60

272,57

1755,84

7

Batee

91,70

0,51

46,77

140,40

264,60

272,57

1755,84

No

Abasah (m2)

Q (m3/det)

293,53

271,62

Keterangan

Normal Melimpah pada banjir tahunan Melimpah pada banjir tahunan Melimpah pada banjir tahunan Melimpah pada banjir tahunan

2926,80

1427,03

Normal Melimpah pada banjir tahunan

Tabel 16. Parameter Hidrolis Sungai Kondisi Normal

1 2 3 4 5

Tripa

Berat Jenis (Sg) 2,3

Tripa

2,8

0,0009

3,10

1,55

74,40

0,55

0,100

sub kritis

Seumayam

2,62

0,0027

2,10

1,05

44,00

0,43

0,095

sub kritis

Seumayam

2,71

0,0021

2,50

1,25

63,20

0,46

0,093

sub kritis

Seumayam

2,66

0,0003

2,80

1,40

50,40

0,33

0,063

sub kritis

6

Batee

2,43

0,0028

1,80

0,90

74,70

0,58

0,138

sub kritis

Batee

2,65

0,0008

1,70

0,85

80,70

0,51

0,125

sub kritis

No Lokasi

7

Parameter Hidrolis Ket.

I (%)

H (m)

R (m)

b (m)

V (m/dt)

Fr

0,0021

4,20

2,10

77,90

0,52

0,081

sub kritis

Tabel 17. Parameter Hidrolis Sungai Kondisi Banjir Tahunan No

Lokasi

Berat Jenis (Sg)

Parameter Hidrolis I (%)

H (m)

Ket.

R (m)

b (m)

V (m/dt)

Fr

1

Tripa

2,3

0,0021

6,60

3,30

77,90

4,617

0,574

sub kritis

2

Tripa

2,8

0,0009

5,00

2,50

74,40

2,441

0,349

sub kritis

3

Seumayam

2,62

0,0027

4,40

2,20

44,00

3,995

0,608

sub kritis

4

Seumayam

2,71

0,0021

4,50

2,25

63,20

3,594

0,541

sub kritis

5

Seumayam

2,66

0,0003

6,68

3,34

50,40

1,700

0,210

sub kritis

6

Batee

2,43

0,0028

5,58

2,79

74,70

4,768

0,644

sub kritis

7

Batee

2,65

0,0008

4,78

2,39

80,70

2,299

0,336

sub kritis

LAPORAN UTAMA


| 343

Karakteristik Sedimen Sedimentasi dapat didefinisikan sebagai pengangkutan, melayangnya atau mengendapnya material fragmental oleh air. Sedimentasi merupakan akibat adanya erosi dan menimbulkan banyak dampak sebagai berikut: (1) Di sungai; pengendapan di sungai mengakibatkan naiknya dasar sungai, meningginya muka air dan mengakibatkan banjir. Hal ini dapat pula menyebabkan aliran terganggu dan mencari alur baru, 2) Di saluran;jika saluran irigasi dialiri oleh air yang banyak sedimen maka akan terjadi pengendapan sedimen di dasar saluran sehingga akan diperlukan biaya yang cukup besar untuk pengerukan tersebut, (3) Di daerah sepanjang sungai; daerah yang dilindungi oleh tanggul akan aman, selama tanggulnya selalu dipertinggi sesuai dengan kenaikan dasar sungai, dan permukaan airnya akan mempengaruhi drainase daerah sekitarnya. Lama kelamaan drainase dengan cara gravitasi tidak memungkinkan lagi. Untuk transportasi sedimen di sungai terdiri dari 2 (dua) cara yaitu sedimen melayang (suspended load) dan sedimen dasar (bed load). Pada suspended load butir bergerak di atas dasar secara melayang. Berat butir terus menerus dipengaruhi oleh gerak turbulensi air (Brownlie, 1981). Suspended load ini dapat di ukur tetapi menimbulkan kesukaran dalam perhitungan. Pada bed load butir bergerak di dasar secara menggelinding (rolling); menggeser (sliding); meloncat (jumping). Bed load ini dapat dihitung dengan rumus semi empirik, tetapi pada pengukuran menimbulkan kesulitan. Beban sedimen dasar (bedload) biasanya terdiri dari lempung, lanau atau pasir dengan ukuran butiran (0,06-2) mm dan kerikil dengan ukuran (2-20) mm. Ukuran, bentuk dan berat jenis sedimen berbeda-beda tergantung dari tempat sedimen tersebut berada. Keadaan besar butiran/gradasi bahan dan sifat pergerakannya berhubungan dengan sifat aliran. Pergerakan sedimen dalam air tergantung kepada kecepatan dan arah aliran, turbulensi aliran, butir dan gradasi bahan sedimen, berat jenis sedimen, bentuk butir sedimen, suhu dan sebagainya. Perhitungan Sedimen Dasar (bed load) Pada dasarnya butiran-butiran tidak akan bergerak jika kecepatan aliran tidak begitu kuat membuat butiran-butiran tersebut bergerak. Keadaan di mana butiran pada saat diam dan akan bergerak dinamakan kecepatan aliran kritis. Partikel yang terbawa pada suatu daerah dapat merupakan partikel yang datang dari daerah lain atau partikel yang diam (mengendap) pada daerah itu. Pada kondisi keseimbangan dapat dikatakan bahwa; partikel yang bergerak dengan partikel yang diam mempunyai jumlah yang sama. Perhitungan analisis sedimen dasar (bed load) dilakukan dengan menggunakan 2 metode. Metode yang digunakan yaitu MPM, dan metode Frijlink, sehingga dapat dilihat perbandingan besarnya sedimen menurut masing-masing metode. Kedua metode ini di sarankan untuk di uji lebih lanjut sesuai dengan kondisi di lapangan (terutama pada Krueng Batee). Metode Frijlink . Frijlink pada tahun 1952 mengeluarkan rumus untuk menghitung sedimen bed load dengan memperhatikan pengaruh konfigurasi dasar sungai (Leopold et al 1957 dalam Chang, 1986), dapat dirumuskan sebagai:

Tb  ( S  9) …………………………………… 39) dm gRI Project Implementation Unit - Studi Ekosistem Rawa Tripa Universitas Syiah Kuala

344


Dimana: dm = Diameter butiran rata-rata (d50) R = Jari – jari hidrolis sungai (m) I = Kemiringan sungai (%), T’b = Berat sedimen (padat) dalam air tiap satuan lebar tiap satuan waktu (m3/m’/dt), (S-9) = Grafik yang menghubungkan parameter intensitas aliran (    ) ripple dan intensitas bed load,

 = Ripple yang merupakan hubungan dari  

dm dan dengan menggunakan grafik h.I

(S-10) untuk parameter Niger dan Benue, maka akan didapatkan Ripple Faktor (  ) pada sungai. Langkah perhitungan jumlah sedimen pada pias Kr. Seumayam dengan menggunakan metode Frijlink dalah sebagai berikut: γs = 2,666 t/m3; γw = 1

Q w  s Q

Q 

t/m3; R= H;  s   1 Q

   R I = 0,047 (s - w ) dm + 0,25 

 w     g 

1/ 3

(Tb)2/3

 1  1x1x0,484x2,600x0,010 = 0,047x1,666x0.0791x10-3+0,25x    9,81  Tb= 0,0035 t/m.det

1/ 3

x Tb2/3

Hasil perhitungan sedimentasi dengan metode Frijlink untuk Sungai Krueng Seumayam dapat dilihat pada Tabel 18. Metode Meyer-Peter dan Mouller (MPM). Meyer Peter & Mueller mengemukakan rumus untuk menghitung total sedimen dalam bentuk bed load, dengan menggunakan prinsip kehilangan tenaga akibat bentuk dasar sungai dan akibat gesekan dengan butir di dasar sungai (Leopold et al 1957 dalam Chang, 1986), dapat dirumuskan sebagai berikut.

 Qs   Ks  wx  x   Q   K's 

3/ 2

 w  2/3 xhxI  0,047 x(s  w)dm  0,25 x  xTb …….. 40)  g 

Dimana: w = Berat jenis air (1000 kg/m3 atau 1 ton/m3)

Qs = Faktor koreksi yang berhubungan dengan tampang saluran biasanya besarnya = 1 Q Ks   = Ripple Faktor K's dm = Diameter butiran median (d50) s = Berat jenis sedimen (gravity spesifik) T’b = Berat sedimen (padat) dalam air tiap satuan lebar tiap satuan waktu (t/m’/dt)

LAPORAN UTAMA


| 345

Tabel 18. Perhitungan Sedimentasi dengan Metode Frijlink untuk Sungai Krueng Seumayam Titik Pengamatan

Parameter

1

2

3

4

5

6

7

D50

0.0791

0.0441

0.1124

0.3625

0.0310

0.0310

0.2551

m

0.484

0.481

0.385

0.593

0.387

0.351

0.534

gw

1

1

1

1

1

1

1

Qs/Q

1

1

1

1

1

1

1

R=H

2.6000

2.1000

5.4000

5.5020

4.6820

5.5820

4.7820

gw*(Qs/Q)*m*R*I

0.0027

0.0026

0.0011

0.0013

0.0140

0.0149

0.0128

s

2.666

3.017

2.736

2.825

2.656

2.656

2.178

.047*(s-gw)*d

0.000006 0.000004 0.000009 0.000031 0.000002 0.000002 0.000014

.25*(gw/g)^(1/3)

0.1168

0.1168

0.1168

0.1168

0.1168

0.1168

0.1168

0.003

0.003

0.001

0.001

0.014

0.015

0.013

(Tb)^(2/3)

0.023

0.023

0.009

0.011

0.120

0.127

0.109

Tb

0.0035

0.0034

0.0009

0.0011

0.0417

0.0455

0.0360

Sumber : Hasil Perhitungan

Perhitungan sedimen bed load dipakai untuk pias Kr. Seumayam sebagai berikut: (1) Lebar sungai (B) = 44.35 m, (2) Kedalaman rata-rata sungai (h) = 2.10 m, (3) Jari – jari hidrolis sungai (R) = 1.30 m, (4) Kemiringan memanjang sungai (I) = 0.0392, (5) Diameter butiran; d35 = 0.118 mm, d50 = 0.132, mm, d90 = 0.184 mm, (6) s  2.666 kg/m3.



s  w .2666  1000 =

w

1000

 1,666

Langkah perhitungan untuk pias 1.2 adalah sebagai berikut: 1. Kontrol H (Tinggi) terhadap B (Lebar) sungai

H x 100% B

=

2,10 x100% 44.35

= 0.047 % < 0.050 % ; maka R = H = 2,10 meter 2. Kecepatan geser butiran I = 0,0392 U* =

gRI = 9,81x2,10 x0,0392 = 0.1626 m/det

3. Kecepatan geser butiran kritis Untuk dmedian = d50 = 0,079 mm ;



 s   w 2666  1000  1,666 = w 1000

U*cr 2 = 0,0034 ; dari grafik S3 Δgd U*cr =

0,009 x1,666 x9,81x0,079  0,0034 m/det U*cr< U*


(butiran Bergerak)

346


4. Menghitung ripple factor (  )

U = 33.5778 RI UR 33.5778 x 2,60 = = 45332583 Re  υ 10 6 C

Dari grafik Sd ; R/k = 2

Dari grafik Sd ; R/k

k = 1.300

1 12R 12x 2,600 = 18 Log = 24,844 m 2 /det k 1,300 1 12R 12 x 2,600 2 Cd90  18 Log = 18 Log = 40,319 /det m d90 0,18 x10 3

C  18 Log

 C     Cd 90 

3/ 2

 24,844  =    40,319 

3/ 2

= 0,484

5. Menghitung volume timbunan sedimen berdasarkan Frijlink

 

0,07912  d m  2660  1000    = 0,1011 = x  1000 RI    0,484 x2,60 x0,010  dari grafik S9 ;

Tb =5 d m μgRI Tb = 5 x 0,07912x10-3 x 0,484 x9,81x2,600 x0,010 = 0,172 t/m.det Dengan cara yang sama perhitungan dilakukan pada sampel pias yang lain, sehingga diperoleh Estimasi Debit Bed Load pada sungai-sungai di kawasan studi seperti terlihat pada Tabel 19. Selanjutnya, rekapitulasi perhitungan jumlah sedimentasi pada kondisi normal (m3/det) dengan menggunakan metode Frijlink dan MPM dapat dilihat pada Tabel 20. Suspended Load Suspended load adalah bagian dari jumlah sedimen yang bergerak tanpa kontak dengan dasar saluran dan gerakannya dipengaruhi oleh gerak turbulensi air. Secara umum, ratarata konsentrasi sedimen melayang seluruh penampang pengukuran (Leopold et al. 1957 dalam Chang, 1986) dapat dihitung dengan menggunakan persamaan berikut ini. C=

 Cq q p qp

…………………………………………………..……41)

Keterangan : C = konsentrasi rata-rata, (kg/m3); qp = debit di sub bagian penampang, (m3/dt); Cq = konsentrasi rata-rata di vertikal pada sub bagian penampang, (kg/m3);

LAPORAN UTAMA


| 347

Tabel 19. Perhitungan Sedimentasi Metode MPM sampel

1

2

3

4

5

6

7

H

2.6000

2.1000

5.4000

5.5020

4.6820

5.5820

4.7820 88.45

B

71.35

71.35

44.00

44.00

88.45

88.45

<0.05 R=H

0.036

0.029

0.123

0.125

0.053

0.063

0.054

U*

0.1626

0.1611

0.1042

0.1125

0.3712

0.3821

0.3537

D50

0.07912

0.04406

0.11242

0.36251

0.03101

0.03101

0.25514

Grafik S3

0.009

0.200

0.008

0.032

0.037

0.037

0.039

s

2.666

3.017

2.736

2.825

2.656

2.656

2.178

U*cr

0.0034

0.0132

0.0039

0.0144

0.0043

0.0043

0.0107

ket

bergerak

bergerak

bergerak

bergerak

bergerak

bergerak

bergerak

U

1.74

1.67

1.26

1.37

4.39

4.65

4.20

C

33.5778

32.4036

37.9277

38.0462

37.0364

38.1378

37.1671

Re

4533253.814

3500322.711

6816714.761

7515615.201

20551189.56

25973157.06

20070474.55

R/k

2

2

2

2

2

2

2

k

1.300

1.050

2.700

2.751

2.341

2.791

12R/k

24.00

24.00

24.00

24.00

24.00

24.00

24.00

24.8438

24.8438

24.8438

24.8438

24.8438

24.8438

24.8438

c D90

2.391

0.180

0.142

0.160

0.732

0.142

0.112

0.458

12R/D90

173.7470

177.3610

405.6815

90.2184

395.4306

595.8317

125.2489 37.7599

C

40.3185

40.4794

46.9473

35.1953

46.7473

49.9522

m

0.484

0.481

0.385

0.593

0.387

0.351

0.534

y

0.1011

0.0699

0.4577

0.8654

0.0094

0.0098

0.0442

5

4.8

4.8

4.5

4.8

4.8

4.8

0.172

0.082

0.301

1.140

0.078

0.081

0.757

Grafik S9 Tb Sumber

: Hasil Perhitungan

Tabel 20. Rekapitulasi Perhitungan Jumlah Sedimentasi pada Kondisi Normal (m3/det) No

Lokasi pengukuran

frijlink

MPM

Rerata

1 2 3 4 5 6 7

Tripa Tripa Seumayam Seumayam Seumayam Batee Batee

0.1716 0.0822 0.3010 1.1399 0.0776 0.0806 0.7568

0.0035 0.0034 0.0009 0.0011 0.0417 0.0455 0.0360

0.0876 0.0428 0.1509 0.5705 0.0596 0.0630 0.3964

Debit Sedimen Melayang Pada periode waktu-waktu tertentu debit muatan sedimen melayang dapat didefinisikan sebagai hasil perkalian konsentrasi dan debitnya yang dapat dirumuskan sebagai berikut (Blench 1961 dalam Gregory, 1977). Qsi = kCQi …………………………………………………….………… 42) Dimana : Qsi = debit sedimen melayang; K = faktor konversi; C = konsentrasi sedimen melayang; Qi = Debit.


348


Umumnya untuk perhitungan debit sedimen melayang pengukuran, ditulis sebagai berikut: Qs = 0,0864 x C x Q ……………………………………… 43) Dimana : Qs=debit sedimen melayang (ton/hari); C=konsentrasi sedimen melayang (mg/l atau g/m3); Q=debit (m3/dt). Apabila didalam periode 1 hari dilaksanakan pengukuran debit, maka debit rata-rata pada hari itu dapat dihitung dengan persamaan sebagai berikut:

Q

1 n  Qi ti …………………………………….....… 44) 24 1

Dimana : = debit rata-rata harian, (m3/dt); = debit yang diukur pada saat t i , (m3/dt);

Q Qi

t i = interval waktu pengukuran debit, (jam); n = banyaknya waktu pengukuran debit. Konsentrasi sedimen rata-rata harian dapat dihitung dengan persamaan sebagai berikut:

C Dimana : C Ci

1 n  Ci ti ………………………………….......…. 45) 24 1 = konsentrasi sedimen melayang rerata harian, (mg/l); = konsentrasi sedimen pada saat t i , (mg/l);

t i

= interval waktu pengukuran debit, (jam);

n

= banyaknya waktu pengukuran debit;

Oleh karena Qi dan Ci kedua-duanya tidak tetap selama periode waktu 24 jam, maka besarnya debit sedimen rata-rata hariannya dihitung dengan persamaan sebagai berikut(Blench, 1961 dalam Gregory , 1977): n

Qs =

 1

0,0864 Ci Qi ti …………………………..… 46) 24

Keterangan : Qs = debit sedimen melayang rata-rata harian, (ton/hari); Ci = konsentrasi sedimen pada saat t i , (mg/l); Qi

=

t i

=

debit yang diukur pada saat t i , (m3/dt); interval waktu pengukuran debit, (jam);

n

=

banyaknya waktu pengukuran debit.

Hasil akhir dari persamaan di atas tidak sama dengan nilai debit sedimen rata-rata harian yang dihitung dengan persamaan berikut (Blench, 1961 dalam Gregory, 1977). Selanjutnya, hasil perhitungan suspended load dapat dilihat pada Tabel 21. Qs = 0,0864 C Qw …………………………………………. 47)

LAPORAN UTAMA


| 349

Keterangan : Qs = debit sedimen melayang rata-rata harian, (ton/hari); C = konsentrasi rata-rata harian, (mg/l); Qw = debit rata-rata, (m3/dt).

Tabel 21. Perhitungan Suspended Load (m3/det) No.

Lokasi Pengukuran

1 Tripa 2 Tripa 3 Seumayam 4 Seumayam 5 Seumayam 6 Batee 7 Batee Sumber : Hasil Perhitungan

Turb (mg/ltr)

Ts (m3/det)

Ts (m3/hari)

257.40 242.50 236.00 230.70 271.10 253.50 204.30

0.927 0.873 0.850 0.831 0.976 0.913 0.735

80,061.70 75,427.20 73,405.44 71,756.93 84,322.94 78,848.64 63,545.47

Sedimen Total Beban sedimen total merupakan gabungan dari bedload dan suspended load sehingga dapat dituliskan (Blench 1961 dalam Gregory, 1977): Tt= Tb + Ts ………………………………………………………….. 48) Dimana: Tb = jumlah sedimen bed load (m3/dt) Ts = jumlah sedimen suspended load (m3/dt) Tt = jumlah sedimen bed load dan sedimen suspended load (m3/dt) Jumlah sedimen yang terangkut pada setiap lokasi adalah kumulasi antara jumlah sedimen suspended dan bed load. Besarnya jumlah sedimen pada setiap pias pengamatan dapat dilihat pada Tabel 22. Tabel 22. Sedimen Total (m3/det) No

Lokasi Pengukuran

Bed Load

Suspended Load

Jumlah

1 2 3 4

Tripa Tripa Seumayam

3.741 1.829 3.977

0.927 0.873 0.850

4.668 2.702 4.826

Seumayam Seumayam Batee

15.032 3.158 3.338

0.831 0.976 0.913

15.862 4.134 4.251

5 6

Parameter Hidrolika Sungai Parameter hidraulis sungai di WS Tripa – Batee ditentukan berdasarkan parameter yang didapat pada pengukuran/survey hidrometri sesaat. Parameter ini meliputi kecepatan geser dan debit dan faktor pembentuk meandering sungai.


350


Hubungan Kecepatan Geser dan Debit Secara alamiah, debit sungai tidaklah tetap, melainkan selalu berubah dengan kecepatan Hubungan Kecepatan Kritis dengan Debit gesernya. Intensitas transpor sedimen juga akan berubah dengan perubahan debitnya. Grafik hubungan kecepatan geser dengan debit dapat dilihat pada Gambar 12.

2500

Debit

2000 1500 1000 500 0 0.0034

0.0132

0.0039

0.0144

0.0043 0.0043

0.0107

Kecepatan Kritis Gambar 12. Grafik Hubungan Kecepatan Kritis dengan Debit Dari hasil analisa data dan perhitungan yang telah didapatkan, maka dapat juga diperoleh hubungan antara beberapa parameter yang dijelaskan dalam bentuk grafik yaitu : (1) grafik hubungan antara debit aliran (Q) dengan kecepatan aliran (V), (2) grafik hubungan antara debit aliran (Q)dengan sedimen bedload (Tb) dan (3) grafik hubungan antara kecepatan aliran (V) dengan bilangan Froude (Fr). Grafik hubungan Kecepatan dengan Debit (V-Q) disajikan pada Gambar 13. Dari grafik QV biasanya terlihat bahwa hubungan antara debit aliran dan kecepatan aliran berbanding lurus, yaitu semakin besar kecepatan aliran (V) maka debit (Q) yang dihasilkan akan semakin besar juga begitu juga sebaliknya. Selanjutnya, hubungan bilangan Froude dengan kecepatan dapat dilihat pada Gambar 14. Pada grafik V-Fr biasanya terlihat bahwa hubungan antara kecepatan aliran (V) dengan bilangan Froude (Fr) berbanding lurus, yaitu semakin besar kecepatan aliran (V) maka bilangan Froude (Fr) yang dihasilkan akan semakin besar, begitu juga sebaliknya. GRAFIK HUBUNGAN KECEPATAN DENGAN DEBIT

DEBIT (m^3/det)

2500

0. 6364x

y = 115.79e R 2 = 0.9751

2000 1500 1000 500 0 0.00

1.00

2.00

3.00

4.00

5.00

KECEPATAN (m/det)

Gambar 13. Grafik Hubungan Kecepatan dengan Debit

LAPORAN UTAMA


| 351

Gambar 14. Hubungan Bilangan Froude dengan Kecepatan Pembentukan Meander Sungai Pada analisa debit sesaat ada analisa tentang gerak turbulensi air yang dapat menjelaskan pembagian kecepatan dan nilai rata-rata kecepatan. Dalam keadaan ini keadaan batas (boundary) yang sangat menentukan adalah dasar/dinding saluran. Dalam hidraulika sedimen diketahui bahwa boundary tersebut terdiri atas sedimen granuler sehingga terjadi pengaruh timbal balik antara dasar dan aliran yang tampak dari perubahan kekasaran dasar. Gerak longitudinal dari sungai seiring dengan gerakan yang membujur, dan ini dinyatakan sebagai pola penting bagi saluran yang dapat memberi pengaruh secara langsung. Leopotd et al (1964 dalam Gregory, 1977) menginformasikan tiga macam bentuk saluran yaitu meandering, jalinan (braided), dan berliku-liku (sinuous). Bentuk saluran ini dapat diidentifikasi dengan melihat ciri dari pola penting yakni perbandingan antara panjangnya saluran dengan panjangnya lembah. Jika perbandingan ini sama dengan atau lebih besar dari 1,5 maka sungai ini disebut sungai dengan pola aliran meander, jika lebih kecil dari 1,5 maka disebut sebagai sungai dengan pola aliran berliku. Untuk debit yang terjadi pada sungai yang berliku-liku dengan slope yang lebih kecil berpotensi menyebabkan pengendapan yang lebih besar, sedangkan pada slope besar akan terjadi transpor sedimen yang tinggi. Kondisi slope sungai di Rawa Tripa cenderung kecil. Kondisi ini menyebabkan proses pengendapan yang dominan, terutama di sungai Tripa, Seumayam dan Batee. Kondisi pengendapan ini semakin berpengaruh pada pias sungai yang berdekatan dengan muara. Selanjutnya, konsep morfologi sungai dapat dilihat pada Gambar 15. Kondisi yang berbeda terjadi bila kita meninjau pada zona peralihan pegunungan ke kawasan yang lebih datar di DAS-DAS yang ada di Rawa Tripa. Pada zona ini, slope sungai yang besar menyebabkan sedimen transport yang terjadi lebih dominan disebabkan oleh erosi di badan sungai. Pias sungai yang berada di hilir dan lembah DAS-DAS Tripa – Batee mempunyai badan sungai yang lebih lebar jika dibandingkan dengan sungaisungai yang terdapat di hulu. Hal ini disebabkan oleh pengangkutan sedimen yang tinggi, sungai akan mengalami pelebaran yang besar dan kedalaman menjadi sangat dangkal. Pembentukan meander akibat dari proses pengikisan dan pengendapan oleh arus sungai. Meander biasanya terbentuk di daerah lembah yang memiliki kelerengan yang tidak terlalu tajam. Proses pengerjaan sungai ini dimulai dari arah hulu dengan arus yang deras kemudian ketika memasuki bagian tengah, relief lembah semakin datar dan


352


kecepatan aliran sungai mulai menurun. Aliran sungai ini membentuk alur di lekukan dalam dan lekukan luar sungai. Lekukan luar memiliki kecepatan aliran yang lebih besar dibanding lekukan dalam, sehingga di bagian lekukan dalam sering terjadi erosi.

Tebing sungai Erosi tebing

sungai

Pengendapan

Gambar 15. Konsep Morfologi Sungai Perkembangan Konfigurasi Dasar Sungai Pada kondisi awal diketahui diperoleh informasi bahwa air bersih/tawar, dasar rata, bed material granuler/non kohesif dan butir material seragam. Klasifikasi tahapan kejadian konfigurasi dasar sungai diuraikan pada Tabel 23. Tabel 23. Klasifikasi Tahapan Kejadian Konfigurasi Dasar Sungai No 1

2

Flow Regime/Tahapan

Keterangan

Lower Flow Regime (< 0.4-0.1) - Tahap 1: plane bed(dasar rata)

Gerak butir masing-masing secara menggelinding, menggeser atau meloncat, at random pada sembarang tempat. Frekuensi gerak butir bertambah. Untuk bahan yang sangat halus dapat terjadi saltasi (saltation) – awan (clouds) – suspended load. Dasar belum mengalami deformasi.

- Tahap 2 : Ripple (dasar bergelombang)

Jika 0 bertambah, di beberapa tempat secara random dan simultan terjadi gundukan kecil. Lama-lama dasar bergelombang teratur dengan amplitudo relatif kecil terhadap panjang gelombang Hr<> d; Makin lama makin tak teratur dan berdimensi 3.

- Tahap 3: Dunes

Sisi sebelah hulu lebih landai dan sisi sebelah hilir lebih curam. Bentuk kurang teratur dan a-simentrik. Kekasaran bertambah, Hr membesar, Lr memanjang. Terjadi erosi di sisi sebelah hulu dan endapan di bagian bawah sisi sebelah hilir, terlihat dari Gambar 16.

Upper Flow Regime (>0.4-1)

Terjadi di Hilir WS Tripa – Bateeterutama pada deskripsi geologi aluvium.

- Tahap 4: Sheet Flow

Bila 0 bertambah besar maka dasar yang bergelombang menjadi rata lagi dengan Tb sangat besar. Kekasaran sangat berkurang. Dapat dianggap k mendekati d. Tahap ini disebut transisi. Fr mendekati 1.

- Tahap 5 : Anti dunes

Jika kecepatan bertambah besar lagi, hingga Fr > 1 maka aliran menjadi tak stabil. Interaksi antara gelombang air dan dasar menghasilkan anti dunes dengan bentuk gelombang pasir agak simetrik. Endapan terjadi pada sisi hulu dan erosi perlu sisi hilir sehingga seolah-olah gelombang pasir merambat ke arah hulu. Pada muka air dapat terjadi standing wave atau breaking wave.

- Tahap 6: Chute and Pool

Bentuk pool dapat dilihat pada Gambar 17.

LAPORAN UTAMA


| 353

Gambar 16. Bentuk Dunes

Gambar 17. Bentuk Pool E. Analisa Kualitas Air Pengelolaan kualitas air dan pengendalian pencemaran air diperlukan untuk perlindungan pemanfaatan air dan ekosistem sumber daya air. Kegiatan pemantauan kualitas air untuk mendukung pengambilan kebijakan pengelolaan Rawa Tripa di Kabupaten Nagan Raya dan Kabupaten Aceh Barat Daya. Berdasarkan pasal 4 ayat 3 PP Nomor 81 tahun 2001 tentang Pengelolaan Kualitas Air dan Pengendalian Pencemaran Air berbunyi “ Upaya pengelolaan kualitas air sebagaimana dimaksud dalam ayat (1) dilakukan pada : (1) sumber air yang terdapat di dalam hutan lindung; (1) mata air yang terdapat di luar hutan lindung; dan (3) akuifer air tanah dalam”.


354


Pengaturan yang terkait dengan kualitas air meliputi zoning, regulasi, peraturanperaturan spesifik tentang air dan tanah, pengendalian, perizinan, larangan dan lisensi. Pengaturan ini menyangkut hal-hal yang berkaitan dengan pencemaran air baik dari polusi sumber titik (point source pollution) maupun polusi bukan sumber titik (non-point source pollution. Pengertian sumber titik adalah titik-titik dimana limbah cair dikeluarkan, misalnya limbah cair yang dialirkan melalui pipa-pipa dan dibuang ke sungai. Sedangkan polusi bukan sumber titik adalah polusi yang dihasilkan dari suatu kawasan tertentu. Zat pencemar atau kontaminan dapat dikelompokkan dalam beberapa kategori. Zat pencemar yang ditinjau disini adalah zat pencemar yang berbentuk cair atau dapat larut dalam air, yang dapat dibagi menjadi : (1) Kontaminan anorganik;(2) Kontaminan organik; (3) Material radioaktif, dan (4) Mikroorganisme. Selanjutnya, penentuan baku mutu pemanfaatan air mengacu kepada Pasal 1 PP No. 82 tahun 2001. Berdasarkan hasil pengujian di lapangan, kualitas air di sungai-sungai utama Rawa Tripa berada dalam kondisi baik sampai tercemar ringan. Pada umumnya parameter pH untuk titik pemantauan ini masih memenuhi baku mutu kualitas air sungai dengan PPRI Nomor 82 Tahun 2001. Terjadi peningkatan jumlah padatan tersuspensi (TSS) yang signifikan pada kondisi sesaat sebelum sampai dengan sesaat setelah banjir. Kondisi ini seiring dengan peningkatan jumlah erosi lahan di hulu yang masuk ke badan sungai. Kondisi sebaliknya terjadi pada jumlah padatan terlarut (TDS). Jumlah padatan terlarut di Sungai Krueng Tripa lebih besar dari sungai-sungai lain di kawasan ini. Hal ini disebabkan pemanfaatan lahan di hulu DAS yang yang sangat dominan untuk pertanian semusim atau pengelolaan tanah yang tidak sesuai dengan kaidah konservasi lahan. Untuk mendukung program pemantauan lanjutan secara lengkap, dilakukan pengambilan sampel yang di uji di Laboratorium Kesehatan Provinsi Aceh. Hasil pengujian selengkapnya disajikan pada Tabel 24, Lmpiran 1, Lampiran 2, dan Lampiran 3. Tabel 24. Hasil Rerata Pengujian Sampel Air Kondisi Normal Parameter Suhu PH TSS TDS Pb Cu Minyak Lemak NH3-N DO COD BOD

Tripa

Nama Sungai Seumayam

Batee

27.8 5,80 28 255 0.001 0.0156 0.23 6.9 7.6 4.4

27.7 6.2 44 130 0.011 0.004 0.18 0.32 6.5 7.3 4.7

27.8 6.43 51 149 0.0012 0.0042 2 0.325 6.9 7.3 4.5

Satuan o

C

mg/l mg/l mg/l mg/l ug/l mg/l mg/l mg/l mg/l

F. Rencana Restorasi Lahan Rawa di TPSF Kajian Antropogenik Penjelasan sebelumnya memberikan informasi bahwa rawa tripa merupakan hilir dari 3 (tiga) DAS utama yakni DAS Krueng Tripa, Krueng Batee dan Krueng Seumayam. Selanjutnya, beberapa DAS lain yang berukuran lebih kecil seperti DAS Krueng Seunaam

LAPORAN UTAMA


| 355

dan DAS Alu Ie Mirah juga melintasi Rawa Tripa. Kawasan Rawa Tripa yang secara administratif berlokasi di Kecamatan Darul Makmur, Kabupaten Nagan Raya dan Kecamatan Bahbarot, Kabupaten Abdya. Desa-desa yang berada di kawasan Rawa Tripa secara umum berlokasi di hilir sungaisungai yang melintasi Rawa Tripa. Hasil dari wawancara dengan penduduk memberikan informasi bahwa beberapa desa yang berlokasi dihilir berdekatan dengan bantaran sungai sering mengalami banjir. Frekuensi banjir yang terjadi setiap tahun pada curah hujan cukup tinggi yang menyebabkan badan sungai tidak mampu lagi mengalirkan limpasan permukaan. Dari konteks hidrologi, kawasan Rawa Tripa dikatagorikan sebagai daerah banjir. Hal ini disebabkan oleh lokasinya yang berada di hilir beberapa sungai yang melintasi Rawa Tripa. Hasil wawancara dan observasi lapangan memberikan informasi bahwa sejumlah desa di Kecamatan Darul Makmur seperti Desa Kuta Trieng dan Desa Seunaam sering mengalami banjir dan tinggi air mencapai 1-2 meter. Frekuensi banjir yang terjadi lebih dari satu kali dalam setahun. Banjir tersebut menyebabkan warga mengungsi ke tempat yang lebih tinggi. Selanjutnya, banjir yang melanda Kecamatan Darul Makmur dan sebagian Kecamatan Bahbarot akibat meluapnya Krueng Ie Mirah menyebabkan penduduk mengungsi. Hasil pemetaan banjir seperti yang dijelaskan sebelumnya memberikan informasi bahwa periode (lamanya) banjir 1-5 hari. Informasi yang diperoleh dari penduduk bahwa lamanya banjir tersebut cenderung lebih panjang di bandingkan dengan lamanya banjir 20 (dua puluh) tahun yang lalu, yang pada umumnya lamanya banjir 1-2 hari. Dengan kata lain, banjir telah surut pada hari ke dua. Jika terjadi banjir, sumber air untuk minum penduduk adalah menggunakan air gallon (aqua) dan air sumur yang ada di lokasi pengungsian. Lama genangan di Kawasan Rawa Tripa disajikan pada Gambar 18. Dari konteks sanitasi lingkungan, sebagian besar penduduk tidak memiliki WC/kakus dan sumur di dalam rumah. Penduduk menggunakan sungai sebagai tempat untuk MCK (mandi, cuci dan kakus). Hal ini menyebabkan kondisi sungai terkontaminasi dengan limbah rumah tangga. Namun, belum diperoleh informasi bahwa kondisi ini menyebabkan terjadinya wabah penyakit seperti disentri dan penyakit kulit yang dialami penduduk. Dalam periode 5 (lima) tahun terakhir, pada musim kemarau sumur penduduk di kawasan Rawa Tripa seperti Desa Alue Bili Kecamatan Bahbarot yang memiliki topografi yang lebih tinggi mengalami penurunan muka air sumur 1-2 meter. Sementara itu, sumur penduduk yang berada di desa yang memiliki topografi yang lebih rendah mengalami penurunan muka air sumur sekitar 20 cm. Penurunan muka sumur pada musim kemarau tersebut belum menjadi permasalahan serius bagi penduduk pada saat ini. Dengan kata lain, air sumur masih dapat diandalkan untuk memenuhi kebutuhan air untuk MCK. Mata pencaharian penduduk pada umumnya adalah petani dan nelayan. Tanaman utama yang dibudidayakan oleh petani antara lain kelapa sawit, kakao, padi, palawija, pisang dan jagung. Namun, fenomena yang terjadi lima tahun terakhir menunjukkan bahwa tanaman kelapa sawit menjadi tanaman yang paling digemari oleh penduduk. Kelapa sawit tersebut pada umumnya dibudidayakan di bantaran sungai-sungai di Kawasan Rawa Tripa. Hasil oberservasi lapangan menunjukkan bahwa bantaran sungai bagian hilir Krueng Tripa, Krueng Batee dan Krueng Seumayam yang digunakan sebagai


356


lahan budidaya kelapa sawit telah mengalami erosi tebing (stream bank erosion) yang menyebabkan tanaman kelapa sawit jatuh ke badan sungai.

Gambar 18. Peta Lama Genangan Banjir di Kawasan Rawa Tripa G. Rencana Restorasi Sistem Tata Air Rawa Manajemen Air Rawa Untuk merestorasi lahan gambut yang sudah rusak akibat kekeringan dan kebakaran, manajemen air rawa merupakan strategi kunci untuk mengembalikan fungsi‐fungsi ekologi dan hidrologi. Selanjutnya, restorasi lahan gambut Tripa melalui manajemen air rawa tersebut harus mempertimbangkan hal-hal sebagai berikut : (1) Mendesain dan memodifikasi kanal dan bangunan air agar kelebihan air dapat dibuang ke dalam keadaan air sungai/laut; (2) Mengatur keseimbangan air di lahan sesuai dengan kebutuhan; dan (3) Penempatan dan pembangunan bangunan air yang berfungsi untuk mengatur keseimbangan air di lahan sesuai dengan kebutuhan . Manajemen air yang diterapkan untuk pengelolaan air di lahan gambut adalah sistem drainase terkendali. Pengontrolan tinggi muka air (water table control) diperlukan untuk menjaga agar supaya tidak terjadi over drainage dan tinggi air tetap konstan sesuai dengan kebutuhan tanaman kelapa sawit. Tinggi muka air yang harus dipertahankan untuk kelapa sawit pada lahan gambut 60 cm (Hasil Kajian 5). Prinsip utama dalam penyekatan parit atau saluran di lahan gambut dapat di lihat pada Gambar 19. Melalui penyekatan ini diharapkan muka air tanah di lahan gambut meningkat dan gambut tidak mengalami kekeringan. Barkah dan Sidiq (2009)

LAPORAN UTAMA


| 357

menyatakan bahwa ukuran bangunan pengendali terutama lebar saluran tergantung komoditas yang diusahakan, untuk tanaman padi memerlukan kondisi lahan tetap tergenang sehingga relatif sempit agar aliran muka air relatif terkendali, dan untuk tanaman perkebunan yang memerlukan kedalaman muka air tanah relatif dalam sehingga perlu dikendalikan sesuai dengan kedalaman zone perakarannya. Ilustrasi manajemen air dengan sistem terkendali dapat dilihat pada Gambar 20.

Sumber : Barkah dan Sidiq, 2009

Gambar 19. Prinsip Utama dalam Penyekatan Parit atau Saluran di Lahan Gambut

Bangunan kontrol Tanggul Sungai

Elevasi muka air

Dasar saluran Elevasi lahan Gambar 20. Manajemen Air dengan Sistem Drainase Terkendali Rekayasa Aliran Sungai pada Manajemen Air di Rawa Tripa Uraian sebelumnya menjelaskan bahwa manajemen air di Rawa Tipa pada prinsipnya ditujukan untuk mengatur keseimbangan air. Melalui pengaturan keseimbangan air ini diharapkan tidak terjadi over drainage dan tinggi air tetap konstan sesuai dengan kebutuhan tanaman kelapa sawit. Tinggi jagaan air maksimum di saluran tidak melebihi 60 cm. Selanjutnya, pengaturan keseimbangan air ini bukan hanya memerlukan modifikasi desain kanal untuk pengendalian drainase, tetapi juga memerlukan suatu


358


rekayasa aliran air yang dapat mengatur pemasukan air (recharge water) dan pengeluaran air (discharge water) secara seimbang dan terkendali pada lahan rawa gambut. Dengan kata lain, pada saat air di lahan rawa gambut dalam kondisi berlebihan di musim penghujan maka air harus dikeluarkan (discharge). Di sisi lain, pada musim kemarau dimana air telah berkurang drastis di lahan rawa gambut dan tidak dapat memenuhi kebutuhan air minimal yang dapat dipertahankan setinggi 60 cm, maka diperlukan pemasukan air (recharge water) dari sungai yang dapat dialirkan ke areal lahan rawa gambut tersebut. Sumber air alternatif untuk menjaga keseimbangan di lahan rawa gambut Tripa adalah sungai yang melintasi lahan rawa gambut. Sungai Krueng Seumayam merupakan sungai yang memiliki pengaruh besar terhadap keseimbangan air di lahan rawa gambut. Oleh karena itu, rekayasa aliran Sungai Krueng Seumayam menjadi alternatif untuk pembuatan bendung. Rekayasa aliran sungai pada manajemen air di Rawa Tripa yang berkaitan dengan deskripsi pola aliran dan tipikal konstruksi bangunan air (bendung) serta penentuan elevasi mercu bendung untuk menaikkan elevasi muka air air sungai agar dapat dialirkan ke lahan rawa gambut diuraikan berikut ini. Deskripsi Pola Aliran. Sistem tata air rawa merupakan bentuk sistem yang terintegrasi antara banjir, topografi dan geologi permukaan. Pada kondisi rawa Tripa, pembentukan rawa akibat banjir sangat dipengaruhi oleh struktur geologi pembentukan dataran banjir (aluvial) yang terlapisi oleh struktur organik (gambut). Satuan endapan aluvium sungai dan pantai di sepadan antara Krueng Tripa sampai Krueng Batee, terdiri dari lempung, pasir, dan kerakal. Deskripsi sayatan geologi di Rawa Tripa dideskripsikan pada Gambar 21.

Gambut

Sungai

Aluvial Batua nindu

k Gambar 21. Sayatan Geologi Permukaan Rawa Tripa Gambar 21 menunjukkan bahwa formasi gambut selalu berada di atas formasi aluvial. Pada titik-titik tertentu, ketinggian air di sungai mencapai lapisan gambut. Pada titik yang lain, aliran sungai berada pada formasi aluvial. Pada kondisi alamiah, terjadi aliran secara lateral dan diagonal di dalam gambut. Proses aliran ini sangat mudah diamati pada musim kemarau. Pada puncak musim kemarau, kemungkinan elevasi aliran air sungai berada di bawah formasi gambut sangat besar. Untuk menaikkan ketinggian air sungai sehingga mencapai lapisan gambut, perlu dilakukan suatu rekayasa terhadap elevasi sungai. Rekayasa aliran sungai bertujuan untuk menaikkan elevasi muka air hingga mencapai lapisan gambut pada saat tertentu. Skematisasi aliran terhadap sistem air rawa pada musim kemarau disajikan pada Gambar 22.

LAPORAN UTAMA


| 359

Rawa Tripa

salura Saluran alami

n alami Saluran Sungai Krueng Seumayam drainase

Samudera Hindia

Sungai Krueng Batee

Saluran

Rawa Tripa drainase Saluran alami

Sungai Krueng Tripa

Gambar 22. Skema Recharge Water dan Discharge Water Rawa Tripa Aktual pada Musim Kemarau Rekayasa aliran tersebut harus tetap memperhatikan kondisi banjir dan dampak terhadap lingkungan. Pada rawa yang dimanfaatkan untuk kegiatan pertanian dan perkebunan, rekayasa aliran sungai tidak dilakukan. Pada kegiatan yang bersifat merehabilitasi lahan dan mengendalikan over drainage, rekayasa aliran berfungsi untuk menjaga muka air tanah dan muka air di saluran transpor utama. Secara alamiah, suplai air ke dalam Rawa Tripa diperoleh melalui hujan dan aliran saluran alamiah. Sedangkan aliran air yang keluar dari sistem rawa terjadi akibat adanya saluran drainase lahan. Saluran drainase lahan yang terdapat di Rawa Tripa terdiri atas saluran drainase alami dan saluran drainase buatan. Permasalahan akan terjadi bila air yang dikeluarkan dari Rawa Tripa lebih besar daripada air yang masuk kedalam sistem rawa. Tipikal Konstruksi. Konstruksi yang direkomendasikan sebagai bangunan untuk merekayasa elevasi aliran sungai di Rawa Tripa adalah bangunan pelimpah sederhana berbentuk bendung. Jenis bangunan yang direkomendasikan berupa bendung tetap atau bendung gerak. Kedua tipe bangunan ini nantinya dapat dipilih salah satu setelah dilakukan perencanaan teknis detail serta dampak bangunan terhadap lingkungan. Pemilihan lokasi bendung untuk keperluan pengendalian recharge ke dalam Rawa Tripa adalah di Sungai Krueng Seumayam. Berdasarkan data topografi parameter hidrologis DAS Krueng Seumayam dapat dilihat pada Tabel 25.


| SCIENTIFIC STUDIES FOR THE REHABILITATION AND MANAGEMENT OF

360

THE TRIPA PEAT-SWAMP FOREST

Tabel 25. Parameter Hidrologis DAS Krueng Seumayam No 1 2 3 4 5

Parameter Hidrologis

Keterangan

Catchment area (pada site AWLR) Elevasi hulu sungai Elevasi hilir sungai (pada site bendung) Slope rata-rata sungai Panjang sungai utama

= 110,60 km = + 1.660,00 m = + 13,00 m = 0,000041 = 64.17 km

2

Selain kondisi tersebut di atas, sungai Krueng Seumayam dipilih sebagai lokasi bendung karena : (1) merupakan sungai utama yang membelah Rawa Tripa; (2) memiliki tampang sungai yang lebih stabil bila dibandingkan dengan sungai Krueng Tripa dan sungai Krueng Batee; (3) memiliki pengaruh banjir lebih dominan bila dibandingkan dengan sungai Krueng Tripa dan sungai Krueng Batee terhadap kejadian banjir di Rawa Tripa; (4) areal sepadan sungai merupakan kawasan perkebunan yang memiliki ketebalan gambut berkisar antara 0,50 m sampai >2,00 m; (5) merupakan outlet utama bagi 64% areal gambut yang drain sehingga akan memudahkan untuk dilakukan pengawasan dan pengendalian; (6) memiliki stasiun pemantauan debit dan hujan. Penentuan Elevasi Mercu Bendung. Elevasi mercu bending ditentukan berdasarkan: (1) elevasi rawa tertinggi, dan (2) faktor tinggi tekan untuk pembilasan. Uraian masingmasing faktor yang diperhitungkan dapat dilihat pada Tabel 26 dan 27. Dari dua metode perhitungan ditetapkan bahwa elevasi mercu bendung adalah yang tertinggi, yaitu pada elevasi +11,50 m. Selanjutnya, penentuan panjang mercu bendung, lebar lubang dan pilar pembilas, panjang mercu bendung efektif tinggi muka air di udik bendung, nilai jari-jari mercu bendung dan resume perhitungan hidraulik bendung dapat dilihat pada Tabel 28. Tabel 26. Penentuan Elevasi Mercu Bendung berdasarkan Elevasi Rawa Tertinggi No

Faktor yang diperhitungkan

Satuan (m)

1 2 3 4 5 6

Elevasi rawa yang akan diairi Tinggi air di rawa Kehilangan tinggi tekanan dari rawa ke saluran tersier Kehilangan tinggi tekanan dari saluran tersier ke saluran sekunder Kehilangan tinggi tekanan dari saluran sekunder ke saluran induk Kehilangan tinggi tekanan dari slope saluran induk ke sedimen trap Kehilangan tinggi tekanan akibat bangunan ukur Kehilangan tinggi tekanan dari sedimen trap ke intake Kehilangan tinggi tekanan pada intake Kehilangan tinggi tekan akibat eksploitasi Ketinggian elevasi mercu bendung

10,00 0,10 0,10 0,10 0,10 0,15

7 8 9 10

0,40 0,25 0,20 0,10 11,50

LAPORAN UTAMA


| 361

Tabel 27. Penentuan Elevasi Bendung berdasarkan Faktor Tinggi Tekan untuk Pembilasan No

Faktor yang diperhitungkan

1

Penangkap sedimen

2

Cara perhitungan

Keterangan - Panjang = 50,0 m; - Panjang saluran pengantar = 20,0 m; - Pemiringan permukaan sedimen di penangkap sedimen = 0,00015; - Elevasi hilir penangkap sedimen + 8,50 m; - Elevasi muka air penangkap sedimen bagian hilir + 11,0 m. - Elevasi muka air dikantong sedimen bagian udik: + 11,0 + (50x0,00015) = +11,005 m. - Elevasi muka air di udik saluran pengantar/tepat di hilir intake bendung: + 11,0 + (50+20)x0,00015 = +11,0105 m. - Kehilangan tinggi tekan di intake diambil = 0,20 m. - Elevasi muka air di udik intake : + 11,0105+0,20 = +11,2105 m. - Kehilangan tinggi tekan akibat eksploitasi diambil = 0,10 m.

Ketinggian elevasi mercu bendung

+11,2105 + 0,10 = +11,3105 m

Tabel 28. Penentuan Hidraulika Bendung No

Hidraulika Bendung

1

Penentuan Panjang Mercu Bendung

2

Penentuan Lebar Lubang Dan Pilar Pembilas

3

Perhitungan Panjang Mercu Bendung Efektif

4

Perhitungan Tinggi Muka Air di Udik Bendung

5

Koefisien debit pelimpah (C)

Keterangan - Panjang mercu bendung ditentukan 120% lebar sungai rata-rata yang lebarnya bervariasi antara 35 – 40 m, sehingga dipilih lebar sungai = 37,50 m. - Jadi panjang mercu bendung: 1,2 x 37,50 m = 45,0 m. - Lebar bangunan pembilas diambil 10% dari lebar sungai = 0,1x37,5 = 3,75 m. - Pembilas dibuat dua buah masing-masing lebarnya 1,90 m. - Lebar pilar pembilas ditetapkan masing-masing 1,50 m. - Panjang efektif mercu bendung dihitung dengan rumus: Be = Bb – 2 (n kp + ka) He, dengan : - Be = panjang efektif mercu bendung (m); - Bb = panjang bruto mercu bendung (45 m); - k = jumlah pilar pembilas (k = 2); - kp = koefisien kontraksi pilas (kp = 0,01); - ka = koefisien kontraksi pangkal bendung (ka = 0,10); - He = tinggi energi (m). - Perhitungannya menjadi: - Be = Bb – 2 (n kp + ka) - He = 45 – 2 (0,01 + 0,1) - He = 45 – 0,24 He - Elevasi muka air banjir di udik bendung dapat diketahui dengan menghitung tinggi energi dengan persamaan berikut ini. 3/2 - Qd = C.Be.He , dengan - Qd = debit banjir rencana - C = koefisien pelimpah debit. - Be = panjang efektif mercu bendung, m. - He = tinggi energi, m. - C = 3,97 (He/Hd)0,12, dimana He = Ha, dari persamaan tersebut diperoleh nilai C = 2,19; - Tinggi energi dihitung dengan:


362


3/2

-He

=

Qd C.Be

= 2,73 m (perhitungan ‘trial and error’ )

- Be = 45 – 0,24 (2.73) = 44.35 m.

6 7

Penentuan Nilai Jari-Jari Mercu Bendung Resume Perhitungan Hidraulik Bendung

- Tinggi tekanan (desain head), Ha dihitung dengan: 2 - Hd = He – V /2 g dan dianggap V2/2g = 0, maka Ha = He = 2.73 m. - Jadi tinggi muka air di udik bendung Ha = 2.73 m - Elevasi muka air banjir = Elevasi Mercu + Ha = 11,50 + 2.73 = +14,23 m. - Jari-jari mercu bendung ditentukan berdasarkan ha = Ha = 2,73 m; - q = 439,5/44,35 = 9,8 m3/det/m’ diperoleh R = 2,25 m. -

Lebar pembilas 2 x 1,90 = 3.80 m. Lebar pilar pembilas 2 x 1,50 = 3,00 m. Tinggi muka air diudik bending =73 m. Elevasi muka air banjir = + 14.23 m. Tinggi pembendungan = 4.0 m. Kemiringan tubuh bending = 1:1

Pemilihan tipe peredam energi. Jenis sungai merupakan sungai aluvial dengan angkutan sedimen dominan dari fraksi pasir dan kerikil. Dengan sungai yang demikian maka bangunan energi yang dipilih adalah tipe lantai datar dengan ujung berkotak-kotak (MDO). Setelah dilakukan rekayasa terhadap aliran di Seumayam, maka sistem tata air di Rawa Tripa pada puncak musim kemarau mengalami perubahan. Perubahan tersebut adalah terjadinya recharge kedalam sistem rawa dan pengendalian atau penutupan drainase keluar dari rawa. Skema di Rawa Tripa akibat adanya pengaturan recharge dan discharge disajikan pada Gambar 23 yang menunjukkan bahwa seluruh saluran drainase dilakukan pengurangan volume drain. Tindakan ini dilakukan dengan merekayasa saluran yang terdapat di dalam areal perkebunan secara sistem jaringan, pembuatan bangunan pelengkap, rekayasa pelindung tebing atau lainnya. Tujuan dari rekayasa di saluran adalah untuk mempertahankan water table dari atas permukaan tanah tidak kurang dari 60 cm. Konservasi DAS Daerah Aliran Sungai (DAS) memiliki keterkaitan secara menyeluruh antara bagian hulu, tengah dan hilir DAS. Kawasan Rawa Tripa yang merupakan bagian hilir dari DAS-DAS yang melintasinya merupakan daerah dimana tempat terkonsentrasi air sebelum mencapai laut. Kawasan Rawa Tripa didominasi oleh tanah Histosol seluas 34.551 Ha (56,96 %), Entisol seluas 19.458 Ha (32,08 %) dan Inceptisol seluas 6.648 Ha ( 10,96 %), secara rinci disajikan pada Tabel 28. Tabel 28 memberikan informasi bahwa jenis tanah histosol merupakan jenis tanah yang mendominasi di kawasan Rawa Tripa jika dibandingkan dengan dari jenis tanah lainnya. Pembentukan jenis tanah histosol atau dikenal dengan peat soil akibat dari proses pelapukan bahan organik (flora dan fauna) pada kondisi an-aerob yang pernah hidup di Kawasan Tripa pada masa lalu. Selanjutnya, pembentukan tanah Entisol dan Inceptisol pada kawasan Rawa Tripa yang disebut juga sebagai jenis tanah alluvial terbentuk dari bahan alluvium (sedimen) yang mengendap akibat banjir. Secara perlahan, proses endapan sedimen pada kawasan Rawa Tripa membentuk tanah alluvial. Hal ini dapat dilihat dari lapisan tanah alluvial tersebut masih bercampur dengan tanah histosol pada lapisan dasar tanahnya. Sehingga, hal ini mengindikasikan bahwa ratusan tahun yang

LAPORAN UTAMA


| 363

lalu seluruh jenis tanah di kawasan Rawa Tripa pada awalnya adalah termasuk ke dalam jenis tanah Histosol.

Sungai Krueng Batee

salur

Bendung

Samudera Hindia

Saluran Hujan dan saluran alami

RawaTripa drainase

an alami Saluran Sungai Krueng Seumayam drainase

Saluran

Rawa Tripa drainase Hujan dan saluran alami

salur

an alami

Sungai KruengSaluran Tripa

drainase Gambar 23. Skema Recharge dan Discharge Rencana Rawa Tripa pada musim kemarau setelah dibangun bendung Tabel 28. Karakteristik Tanah di Kawasan Rawa Tripa No 1

2

3

Sistem Klasifikasi Tanah Soil Taxanomi USDA Nasional Ordo Great Group Indonesia 2013 Histosol Typic haplo fibrist Organosol fibrik Typic haplo hemist Organosol hemik Typic haplo saprist Organosol saprik Sub total Entisol Typic tropopluvent Alluvial eutrik Typic tropaquent Alluvial gleik Sub total Inceptisol Typic distropepts Alluvial distrik Total

Luas (Ha) 4.931 3.679 25.941 34.551 15.698 3.760 19.458 6.648 60.657

(%) 8,13 6,07 42,77 56,96 25,88 6,20 32,08 10,96 100,0

Sumber : Laporan kajian 5 (2013)

Perubahan dari penggunaan lahan dari rawa menjadi lahan budidaya yang memerlukan drainase (pembuangan air gambut) menyebabkan terjadi perubahan kondisi an-aerob menjadi aerob yang memacu terjadinya pelapukan bahan organik pada tanah gambut. Proses pelapukan bahan organik ini menyebabkan hilangnya bahan organik dan menyebabkan terjadinya subsidance (penurunan muka tanah gambut) dan berkurangnya kemampuan gambut untuk menyerap air. Hasil dari observasi lapangan tentang berbagai ke dalaman tanah di Rawa Tripa disajikan pada Tabel 29.


364


Tabel 29. Kedalaman Bahan Gambut di Lokasi Rawa Tripa No 1 a b c 2

Kedalaman Gambut (cm) < 200 200-300 >300 Tanah mineral Total

Luas Areal (Ha) 2.844, 46 19.411,40 12.296,22 26.105,20 60.657,29

(%) 4,69 32,00 20,27 43,04 100,00

Sumber : Laporan kajian 5 (2013)

Gambut yang memiliki kedalaman > 3 m harus menjadi daerah konservasi, sementara gambut yang < 3 m dapat dipertimbangkan menjadi areal budidaya. Demikian juga daerah bantaran sungai (sisi kiri dan kanan) sungai juga termasuk areal konservasi. Arahan alokasi ruang di areal TPSF dapat dilihat pada Gambar 24. AK

AK

Budidaya

Areal Konservasi (AK)

Budidaya

AK

AK H1

Kubah gambut H2 <3m

Tanah gambut> 3 m

<3m Sungai

Sungai Bantaran Sungai

Tanah mineral

Bantaran Sungai

Gambar 24. Arahan Alokasi Ruang di Lahan TPSF Upaya mitigasi banjir dilakukan merelokasi perumahan penduduk yang yang berada di bantaran sungai bagian hilir di Kawasan Rawa Tripa atau menyesuaikan konstruksi bangunan rumah seperti rumah panggung yang lebih tinggi dari tinggi ancaman banjir yang terjadi. Selanjutnya, lokasi yang perlu dikonservasi khususnya untuk areal sempadan sungai dan sempandan pantai dapat dilihat pada Gambar 25. Analisis Komponen Hidrologi terhadap Emisi Karbon Emisi CO2 yang terjadi di Kawasan Gambut Rawa Tripa dipengaruhi oleh beberapa faktor yaitu: kedalaman muka air gambut, kadar air gambut, temperatur gambut dan electrical conductivity gambut. Secara umum, hubungan berbagai komponen hidrologi menunjukkan bahwa hujan (rainfall) meningkatkan water table dan volumetric water contents (VWC) tanah gambut di areal TPSF, demikian sebaliknya pada kondisi tidak ada hujan maka water table dan VWC tanah gambut mengalami penurunan. Namun, hubungan antara water table dengan CO2 flux masih belum dapat dijelaskan dengan baik karena memiliki nilai koefisien korelasi yang tidak significant (Kajian 9, 2013). Hal ini perlu dikaji lebih dalam karena peneliti lain berpendapat bahwa penurunan water table mempunyai korelasi positif dengan laju emisi karbon (Hooijer et al., 2006)

LAPORAN UTAMA


| 365

Gambar 25. Lokasi Rawa Tripa yang harus Dikonservasi Dari konteks hidrologi, faktor utama yang menyebabkan tingginya emisi CO2 yang terjadi adalah kedalaman drainase karena berkaitan langsung dengan kedalaman air tanah (water table). Gambar 26 memperlihatkan adanya korelasi antara jumlah emisi CO2 (ton/ha/tahun) dengan tinggi muka air tanah gambut pada 3 wilayah geografi berbeda. Dari gambar tersebut terlihat bahwa, semakin dalam air tanah gambut diturunkan dari permukaan, maka semakin besar emisi CO2 yang ditimbulkan. Kondisi demikian dampaknya semakin besar di negara tropis dibandingkan di wilayah boreal dan temperate (Hooijer et al., 2006). Perkebunan kelapa sawit di areal TPSF membangun dua jenis saluran yaitu saluran kolektor dan saluran transport. Fungsi saluran kolektor adalah untuk pengeringan lahan sedangkan saluran transport menerima air dari saluran kolektor untuk dibuang atau dialirkan ke sungai. Saluran kolektor dibuat pada sisi kiri dan kanan jalan blok dengan lebar 2,0 m dan dalam 1,5 m berdinding tanah (non lining) dengan kemiringan talud relatif terjal, yaitu 1H:5V. Kemiringan dasar saluran relatif kecil mendekati datar sehingga aliran air terjadi diakibatkan oleh adanya beda potensial muka air. Saluran transport dibuat sama dengan saluran kolektor hanya lebarnya yang lebih lebar yaitu sekitar 5.0 m (Kajian 6, 2013). Secara umum dapat dikatakan bahwa ke dua jenis saluran drainase tersebut memiliki kemampuan drainase yang besar jika tidak dikendalikan. Hooijer et al. (2006) menyatakaan bahwa pada setiap penambahan kedalaman drainase satu centimeter berpotensi meningkatkan emisi sebesar 0,91 ton CO2/ha/tahun (Gambar 26). Sehubungan dengan hal tersebut, upaya restorasi untuk pengendalian drainase dan rekayasa aliran air untuk mempertahankan tinggi jagaan muka air konstan sangat diperlukan.


366


Gambar 26. Hubungan antara emisi CO2 dengan Kedalaman Air Tanah Gambut Monitoring dan Evaluasi Monitoring dan evaluasi kondisi lahan rawa gambut Tripa perlu dilakukan dengan tujuan antara lain: (1) memonitor apakah rekomendasi kajian ini ditindaklanjuti oleh pihakpihak terkait, (2) mengevaluasi secara lebih mendalam tentang kondisi kekinian areal TPSF pasca kajian setelah rekomendasi rencana restorasi dilaksanakan termasuk kajian dampak positif dan negatifnya terhadap lingkungan di TPSF. Tabel 30. Jadwal Rencana Kegiatan dan Monev Restorasi Areal TPSF

No A 1 2 3 4 5 6 7 B 1 2 3 C 1 2

Aktvitas Restorasi Lahan

Jadwal Kegiatan/Monev (Tahun ke- ) 1 2 3 4 5

Perencanaan : Pengaturan Alokasi Ruang Detil di TPSF DED Konservasl DAS Kajian terhadap water yield DED pengendalian drainase (canal blocking) DED pembangunan bendung AMDAL Konservasi DAS AMDAL Pengendalian drainase dan pembangunan bendung Pelaksanaan: Konservasi DAS Pembangunan pengendali drainase (canal blocking) Pembangunan bendung Monitoring dan Evaluasi: Monev terhadap implementasi UPL dan UKL Pembangunan bendung dan pengendalian drainase Monev terhadap implementasi UPL dan UKL Pembangunan bendung dan pengendalian drainase

LAPORAN UTAMA


| 367

Beberapa rekomendasi untuk restorasi areal TPSF seperti meredesain kembali alokasi ruang di TPSF, sistem tata air untuk pengendalian drainase, rekayasa aliran melalui pembuatan bendung, upaya konservasi DAS serta dampak postif dan negatif terhadap lingkungan muncul akibat pelaksanaaan kegiatan tersebut penting ditindaklanjuti serta dimonitor secara konprehensif. Demikian juga dengan keterlibatan masyarakat dan perusahaaan perkebunan di TPSF dalam upaya restorasi areal TPSF perlu dimonitor dan dievaluasi secara periodik. Bagaimana pun juga, restorasi areal TPSF tidak hanya mempertimbangkan aspek lingkungan, tetapi juga aspek ekonomi, khususnya masyarakat yang berada di sekitar TPSF. Selanjutnya, rekomendasi kegiatan monitoring dan evaluasi tersebut dapat dilihat pada Tabel 30.

IV. KESIMPULAN DAN REKOMENDASI A. Kesimpulan 1. Kawasan Rawa Tripa pada dasarnya merupakan kawasan banjir yang terletak di hilir dari WS Woyla – Tripa. Banjir yang terjadi di kawasan ini merupakan banjir tahunan dengan lama genangan berkisar antara 1 sampai 5 hari; 2. Terdapat 3 sungai utama pada Kawasan Rawa Tripa. Sungai terbesar di kawasan ini adalah sungai Krueng Tripa. Sedangkan 2 sungai besar lainnya adalah sungai Krueng Batee dan sungai Krueng Seumayam. Sungai Krueng Seumayam ini terletak tepat di tengah-tengah Rawa Tripa; 3. Dari sungai-sungai yang terdapat di Rawa Tripa, sungai Krueng Seumayam memiliki pengaruh yang sangat dominan terhadap banjir di Rawa Tripa. Limpasan banjir dari sungai ini sangat mempengaruhi sistem imbangan air rawa dan kualitas air rawa. Perubahan terhadap kualitas air rawa memberikan pengaruh yang sangat besar pada tingkat keasaman lahan di Kawasan Rawa Tripa; 4. Limpasan banjir yang terjadi dari Krueng Tripa mempengaruhi bahagian Barat dari kawasan Rawa Tripa. Kawasan yang dipengaruhi meliputi areal perkebunan PT. GSM, perkebunan campuran milik rakyat dan kawasan sepadan sungai. Banjir dari sungai Tripa membawa sedimen dalam jumlah besar dan memberi suplai tanah mineral ke dalam sistem rawa Tripa; 5. Banjir yang terjadi dari sungai Krueng Batee mempengaruhi sistem rawa yang terletak di sebelah Timur kawasan. Banjir dari Krueng Batee memberi pengaruh sangat dominan pada kawasan gambut di PT. Cemerlang Abadi dan PT. Patriot GSA. Banjir dari sungai Krueng Batee mempengaruhi secara langsung kondisi kolam-kolam rawa yang terletak di sepadan pantai dari kawasan Timur Rawa Tripa; 6. Berdasarkan kondisi imbangan air, alih fungsi penggunaan lahan di rawa Tripa memberi dampak yang sangat ekstrim. Perubahan penggunaan lahan dari kawasan hutan rawa menjadi perkebunan sawit menyebabkan terjadinya pengeringan air dari Rawa Tripa ke sungai-sungai yang ada di kawasan rawa. Proses ini menyebabkan terjadinya penurunan muka tanah akibat berkurangnya air di rawa yang mengakibatkan meningkatnya emisi karbon yang berasal dari TPSF. 7. Rencana restorasi lahan rawa di TPSF dapat diimplementasikan melalui manajemen air yang terkendali, rekayasa aliran sungai pada manajemen tata air, dan konservasi DAS. Manajemen air terkendali dilakukan dengan medesain kembali pola aliran, ukuran dan bentuk kanal. 8. Untuk menjaga keseimbangan air di lahan rawa TPSF dapat dilakukan dengan rekayasa aliran air Sungai Krueng Seumayam dengan pembuatan bendung sederhana yang dapat mengendalikan recharge water dan discharge water di areal TPSF. Project Implementation Unit - Studi Ekosistem Rawa Tripa Universitas Syiah Kuala

368


9. Konservasi DAS di aeral TPSF dapat dilakukan dengan mengatur alokasi ruang gambut sesuai dengan yang direkomendasikan, perbaikan sistem pembuang air rawa, perlindungan sepadan sungai dan pantai. B. Rekomendasi 1. Rencana restorasi lahan di TPSF perlu ditindaklanjuti melalui penyiapan dokumen untuk pengaturan alokasi ruang TPSF skala detil, DED pengendalian drainase, DED pembangunan bendung dan DED konservasi DAS, pelaksanaan kegiatan restorasi serta pelaksanaan monitoring dan evaluasi kegiatan restorasi. 2. Perlu segera dilakukan kajian mengenai hasil air (water yield) terhadap pola penggunaan lahan di TPSF untuk mendukung program monitoring dan evaluasi kegiatan restorasi. 3. Pembukaan lahan baru di TPSF perlu dihentikan dan rehabilitasi lahan yang telah rusak di luar konsesi perkebunan kelapa sawit seperti bantaran sungai dan sempadan pantai serta areal terlantar lainnya perlu segera dilakukan. 4. Kerjasama semua pihak yang melibatkan pihak pemerintah (pemerintah pusat, pemerintah Aceh, pemerintah Kabupaten Nagan Raya dan Aceh Barat Daya), perusahaan kelapa sawit, masyarakat di sekitar areal TPSF, perguruan tinggi serta pihak terkait lainnya sangat diperlukan agar area Hutan Rawa Gambut Tripa (TPSF) dapat bermanfaat secara adil dan berkelanjutan. DAFTAR PUSTAKA Asdak, C (1995). Hidrologi dan pengelolaan DAS, Gadjah Mada University Press, Yokyakarta. Ayob, K dan Wahab, A.K.A (2003). Hydrological Analysis of A Drained Peat Basin Using Time Series Correlation and Cross-correlation Fungtions,Jurnal Teknologi p63-74, Universiti Teknologi Malaysia, Malaysia. Barkah, B.S dan Sidiq, M (1999). Panduan Penyekatan Parit/Kanal dan Pengelolaannya Bersama Masyarakat Di areal Hutan Rawa Gambut MRPP Kabupaten Musi Banyuasin, Dinas Kehutanan Kabupaten Musi Banyuasin, Sumatera Selatan. Blench (1961 dalam Gregory, K. J, 1977). River Channel Change. A Wiley Interscience, Publication, New York, USA. Brooks, K.N., Ffolliott. P.F., Gregersen, H.M., and DeBano, L.F (2003). Hydrology and the management of watersheds. Iowa State. Blackwell Publishing Company, USA. Brownlie, W. R. (1981). Prediction of flow depth and sediment discharge in open channels, Report No. KH-R-43A, W. M. Keck Laboratory of Hydraulics and Water Resources, California Institute of Technology, Pasadena, California, USA. Chow, V.T (1959). Open Channel Hydraulic. McGraw-Hill, International Edition, New York, USA. Chow, V.T., Maidment, D.R., and Mays, L.W (1988). Applied hydrology. McGraw-Hill, International Edition, New York, USA. Clewell, A., Rieger, J., Munro, J (2005). Guidelines for Developing and Managing Ecological Restoration Projects. 2nd Edition. Society for Ecological Restoration International. Driessen, P.M and Suhardjo, H (1976). On The Defective Grain Formation of Sawah Rice on Peat, Bulletin 3. Soil Research Insitute, Bogor. Fredsoe, J. and Deigaard, R (1992). Mechanics of Coastal Sediment Transport. Advanced Series on Ocean Engineerjng - Volume 3. World Scientific Publishing Co. Pe. Ltd, UK office, London.

LAPORAN UTAMA


| 369

Hooijer, A., Silvius, M., Wösten, H. and Page, S (2006). PEAT-CO2, Assessment of CO2 emissions from drained peatlands in SE Asia. Delft Hydraulics report Q3943. Leopod, et al. (1957 dalam Chang H.H. 1986). River Channel Change. A Justment of Equilibirium. J. Hydraulic Engineering. Vol. 112 (1), p 43-45. Leopod, et al. (1964 dalam Gregory, K. J. 1977). River Channel Change. A Wiley Interscience, Publication, New York, USA. Natural Resources Conservation Service. (1986). Conservation engineering division. Technical Release 55 . Urban Hidrology for Small Watershed, US Departement of Agriculture, Washington. Mitsch, W.J., and Gosselink, J. G. (2011). Wetlands. Ecological Studies, Vol. 190, John Wiley & Sons. Perrow, M.R., and Davy, A.J. 2002. Handbook of Ecological Restoration. Volume 1. Principles of restoration. Cambridge: Cambridge University Press. Sri Harto (2000). Hidrologi ; Teori Masalah Penyelesaian, Nafiri, Yokyakarta. Subramanya, K. (2009), Engineering Hydrology, Mc Graw Hill, Singapore. Verhoeven, J.T.A., Beltman, B., Bobbink, R., and Whigham, D.F. (2006) Wetlands and Natural Resource Management, Springer-Verlag Berlin Heidelberg. Ward, A.D., and Elliot, W.I. (1995). Environmental hydrology. Lewis Publishers, CRC Press Inc., New York.


370


Lampiran 1. Hasil Pengujian Kualitas Air di Sungai-Sungai Rawa Tripa NO

PARAMETER

A 1 2 3 4 5 6 B 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Catatan :

FISIKA Bau Rasa Warna Kekeruhan/turbidity Suhu Zat Padatan Terlarut KIMIA Air Raksa (Hg) Arsen (As) Besi (Fe) Flourida (F) Kadmium (Cd) Kesadahan (CaCO3) Klorida (Cl ) Kromium (Cr total) Mangan (Mn) Nitrat (NO3) Nitrit (NO2) pH Seng (Zn) Sianida (Cn) Sulfat (SO4) Timbal (Pb) Selenium (Se) Zat Organik

SATUAN

BAKU MUTU

METODE ANALISA

TCU NTU :C mg/l

Tidak Berbau Tidak Berasa 50 25

mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l

HASIL ANALISA (KODE LOKASI) SM 1 A

SM 2 A

SS 01

STR 4 A

STR 2 A

SB A

1500

Organoleptik Organoleptik Spektofotometri Termometri Termometri Gravimeter

Tidak Berbau Tidak Berasa 13,12 41,2 27,3 25,6

0,625 2,65 27,3 42,8

5 5,57 27,6 34,3

10,16 19,02 27,2 82

15,31 15,32 27,5 66,7

2,187 4,65 27,4 32,3

0,001 0,05 1 1,5 0,005 500 600 0,05 0,5 10 1 6,5 – 9 15 0,1 400 0,05 0,01 10

SSA SSA SSA Photometri SSA Trimetri Argentometri SSA SSA Brucin Spektrofotometri Elektrometri SSA Photometri Turbidimetri SSA SSA Trimetri

TT TT 0,0708 1,03 0,0009 44 5,68 0,001 0,0741 0,061 0,021 6,56 0,0011 0,006 6,544 0,0014 TT 0,217

TT TT 0,01 1,23 0,0003 52 5,68 0,0016 0,1174 0,003 0,008 7,06 0,0016 0,004 2,279 0,0003 TT 0,821

TT TT 0,0502 1,53 0,0001 26 9,23 0,0008 0,1242 0,006 0,023 5,96 0,0012 0,011 0,172 0,0006 TT 0,236

TT TT 0,065 0,74 0,0003 76 8,52 0,0022 0,1068 0,168 0,03 6,85 0,0031 0,043 19,78 0,0004 TT 0,304

TT TT 0,0687 0,49 0,0012 76 11,36 0,0025 0,1121 0,213 0,062 6,72 0,0051 0,03 6,709 0,0008 TT 0,623

TT TT 0,0099 0,65 0,0006 44 6,39 0,0012 0,0991 0,012 0,004 6,8 0,0009 0,007 0,441 TT TT 0,405

Suhu udara ±3 :C

SM = Sungai Krueng Seumayam, SB = Sungai Krueng Batee, SS = Sungai Krueng Seunaam, STR = Sungai Krueng Tripa

LAPORAN UTAMA


| 371

Lampiran 2. Hasil Pengujian Kualitas Air di Saluran Kebun Sawit Rawa Tripa NO

PARAMETER

SATUAN

BAKU MUTU

METODE ANALISA

A FISIKA 1 Bau Tidak Berbau Organoleptik 2 Rasa Tidak Berasa Organoleptik 3 Warna TCU 50 Spektofotometri 4 Kekeruhan/turbidity NTU 25 Termometri 5 Suhu :C Suhu udara ±3 :C Termometri 6 Zat Padatan Terlarut mg/l 1500 Gravimeter B KIMIA 7 Air Raksa (Hg) mg/l 0,001 SSA 8 Arsen (As) mg/l 0,05 SSA 9 Besi (Fe) mg/l 1 SSA 10 Flourida (F) mg/l 1,5 Photometri 11 Kadmium (Cd) mg/l 0,005 SSA 12 Kesadahan (CaCO3) mg/l 500 Trimetri 13 Klorida (Cl ) mg/l 600 Argentometri 14 Kromium (Cr total) mg/l 0,05 SSA 15 Mangan (Mn) mg/l 0,5 SSA 16 Nitrat (NO3) mg/l 10 Brucin 17 Nitrit (NO2) mg/l 1 Spektrofotometri 18 pH 6,5 – 9 Elektrometri 19 Seng (Zn) mg/l 15 SSA 20 Sianida (Cn) mg/l 0,1 Photometri 21 Sulfat (SO4) mg/l 400 Turbidimetri 22 Timbal (Pb) mg/l 0,05 SSA 23 Selenium (Se) mg/l 0,01 SSA 24 Zat Organik mg/l 10 Trimetri Catatan : RO 2, RO 3 dan RO 6 berada di antara sungai Krueng Tripa dan sungai Krueng Seumayam RO 7, RO 12 dan RO 13 berada di antara sungai Krueng Seumayam dan Krueng Batee


RO 2 A

RO 2 C

HASIL ANALISA (KODE LOKASI) RO 3 RO 6 A RO 7 A

RO 12 A

RO 13

Tidak Berbau Tidak Berasa 28,906 3,18 27,9 83,6

35,469 2,8 28,1 24,9

135,14 3,34 27,7 33,6

21,875 3,73 27,7 33,6

6,094 2,75 27,9 24,9

11,72 7,33 27,3 43,2

22,34 8,23 27,5 61,7

TT TT 0,0988 0,77 0,0011 110 9,23 0,0023 0,1397 2,332 0,072 6,46 0,0014 0,004 14,14 0,0005 TT 0,724

TT TT 0,1371 0,69 0,0017 60 7,1 0,0039 0,1654 1,141 0,018 5,93 0,0032 0,003 6,299 0,0004 TT 0,186

TT TT 0,1986 0,39 0,0027 100 7,1 0,004 0,1243 0,855 0,024 5,63 0,0031 0,007 18,43 0,0011 TT 0,245

TT TT 0,0877 0,83 0,0017 26 9,94 0,0038 0,1134 0,523 0,005 5,8 0,0023 0,005 12,55 0,0099 TT 0,512

TT TT 0,021 1,95 0,0016 44 8,52 0,0041 0,1261 0,003 0,001 6,08 0,0008 0,003 7,525 0,0002 TT 0,167

TT TT 0,1104 1 0,0012 50 7,81 0,0037 0,1267 0,614 0,118 6,23 0,0051 0,006 14,26 0,0003 TT 1,186

TT TT 0,1518 2,28 0,0009 56 7,81 0,0035 0,002 1,404 1,966 6,39 0,015 0,012 4,093 0,0006 TT 10,97

372


Lampiran 3. Hasil Pengujian Kualitas Air di lahan dan Genangan di lahan Rawa Tripa NO

PARAMETER

A 1 2 3 4 5 6 B 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24

FISIKA Bau Rasa Warna Kekeruhan/turbidity Suhu Zat Padatan Terlarut KIMIA Air Raksa (Hg) Arsen (As) Besi (Fe) Flourida (F) Kadmium (Cd) Kesadahan (CaCO3) Klorida (Cl ) Kromium (Cr total) Mangan (Mn) Nitrat (NO3) Nitrit (NO2) pH Seng (Zn) Sianida (Cn) Sulfat (SO4) Timbal (Pb) Selenium (Se) Zat Organik

SATUAN

BAKU MUTU

METODE ANALISA

TCU NTU :C mg/l

Tidak Berbau Tidak Berasa 50 25 Suhu udara ±3 :C 1500

mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l mg/l

0,001 0,05 1 1,5 0,005 500 600 0,05 0,5 10 1 6,5 – 9 15 0,1 400 0,05 0,01 10

HASIL ANALISA (KODE LOKASI) 2A 1A

4A

3A

5A

Organoleptik Organoleptik Spektofotometri Termometri Termometri Gravimeter

Tidak Berbau Tidak Berasa 233,44 136 27,4 20,7

418,59 263 27,7 32,4

468,75 323 28 34,1

173,75 170,3 28,1 21,2

182,5 61,8 27,6 33,9

SSA SSA SSA Photometri SSA Trimetri Argentometri SSA SSA Brucin Spektrofotometri Elektrometri SSA Photometri Turbidimetri SSA SSA Trimetri

TT TT 0,0686 2,49 0,0013 560 71 0,0023 0,0964 3,832 0,114 4,86 0,0116 0,122 89,39 0,0004 TT 1,149

TT TT 0,2793 2,88 0,0018 540 8,52 0,0022 0,1158 7,255 0,589 5,29 0,0111 0,211 162,4 0,0014 TT 5,892

TT TT 0,3312 3,36 0,0016 560 78,1 0,0016 0,0676 12,82 0,812 5,27 0,0137 0,255 208,5 0,0006 TT 8,123

TT TT 0,0491 0,95 0,0013 400 78,1 0,0011 0,122 2,961 0,278 5,05 0,012 0,087 95,88 0,0001 TT 2,781

TT TT 0,0964 1,28 0,0008 720 71 0,002 0,0028 2,675 0,187 4,36 0,0179 0,063 87,79 0,0001 TT 1,873

LAPORAN UTAMA

RENCANA RESTORASI LAHAN RAWA

Recommend Documents