## Pemodelan Hidrodinamis 2D Skenario Penggenangan Banjir Luapan Krueng Tamiang
Rifa Maulina 1 , Ella Meilianda 2* , Ahmad Reza Kasury 3
1,2,3 Jurusan Teknik Sipil, Universitas Syiah Kuala, Banda Aceh 23111 Indonesia Email: 1 rifa_m@mhs.unsyiah.ac.id , 2 ella_meilianda@unsyiah.ac.id , 3 ahmadrezakasury@unsyiah.ac.id *Corresponding Author

## Abstract
Flooding is a natural disaster that often occurs in Indonesia. Generally, flood events can be predicted if there is high rainfall in an area and it lasts for several days. However, floods can also occur suddenly if the area is in a damaged  ecosystem.  Currently,  flood  events  in  Indonesia  are  classified  as  extreme.  One  of  the  areas  prone  to flooding is the Krueng Tamiang watershed. This watershed is the largest watershed in Aceh Province with an area of 4,826.56 Km 2 with a main river length of 217 km which has a low soil infiltration capacity in the rainy season, causing high surface flow that can cause flooding. From these problems, this study aims to determine the potential for flood overflow using 2D hydrodynamic modeling for various scenarios of flood discharge return periods in the Krueng Tamiang watershed and review the validation of the flood inundation model against flood location data from field surveys in the Krueng Tamiang watershed based on historical floods using Geo-HECRAS software, and determine flood conditions (depth and extent of flood inundation) at return periods of 2, 5, 10, 25, 50, and 100 years. The modeling results obtained the flood inundation area for each return period as follows L2th = 43.373 km2, L5th = 69.389 km2, L10th = 92.451 km2, L25th = 138.092 km2, L50th = 145.939 km2, and L100th = 179.836 km2. The depth of flood inundation ranges from 0 to more than 2 m. The validation results of flood inundation with flood distribution points resulted in information that 10 points were in the flood inundation area and 2 points were not in the flood inundation area.
Keywords: flood, flood inundation model, Geo-HECRAS.

## Abstrak
Banjir merupakan bencana alam yang sering terjadi di Indonesia. Umumnya kejadian banjir dapat diprediksi jika terjadi curah hujan yang tinggi di suatu kawasan dan berlangsung selama beberapa hari. Akan tetapi, banjir juga dapat terjadi secara tiba-tiba jika kawasan tersebut berada dalam ekosistem yang rusak. Saat ini kejadian banjir di  Indonesia  tergolong  ekstrim.  Salah  satu  daerah  yang  rawan  terhadap  banjir  adalah  Daerah  aliran  Sungai (DAS) Krueng Tamiang. DAS ini merupakan DAS terbesar di Provinsi Aceh dengan luas 4.826,56 Km 2 dengan panjang  sungai  utama  217  km  yang  memiliki  kapasitas  infiltrasi  tanah  yang  rendah  pada  musim  penghujan, sehingga menyebabkan tingginya aliran permukaan yang dapat menimbulkan banjir. Dari permasalahan tersebut, maka  penelitian  ini  bertujuan  untuk  mengetahui  potensi  luapan  banjir  dengan  menggunakan  pemodelan hidrodinamis 2D untuk berbagai skenario periode ulang debit banjir pada DAS Krueng Tamiang dan meninjau validasi  model  genangan  banjir  terhadap  data  lokasi  banjir  dari  survey  lapangan  di  DAS  Krueng  Tamiang berdasarkan banjir historis menggunakan Software Geo-HECRAS, serta mengetahui kondisi banjir (kedalaman dan luasan genangan banjir) pada kala ulang 2, 5, 10, 25, 50, dan 100 tahun. Dari hasil pemodelan didapatkan luas genangan banjir untuk setiap kala ulang sebagai berikut L 2th = 43,373km 2 , L 5th = 69,389 km 2 , L 10th = 92,451 km 2 , L 25th = 138,092 km 2 , L 50th = 145,939 km 2 , dan L 100th = 179,836 km 2 . Kedalaman genangan banjir berkisar 0 sampai dengan lebih dari 2 m. Hasil validasi genangan banjir dengan titik sebaran banjir menghasilkan informasi 10 titik berada pada area genangan banjir serta 2 titik tidak berada pada area genangan banjir.
Kata kunci: banjir, pemodelan genangan banjir, Geo-HECRAS.

## I. Pendahuluan
Banjir merupakan bencana alam yang sering terjadi  di  Indonesia.  Umumnya  kejadian  banjir dapat  diprediksi  jika  terjadi  curah  hujan  yang tinggi di suatu kawasan dan berlangsung selama beberapa  hari.  Akan  tetapi,  banjir  juga  dapat terjadi  secara  tiba-tiba  jika  kawasan  tersebut berada  dalam  ekosistem  yang  rusak.  Dalam
kurun  waktu  10  tahun  terakhir  menunjukkan bahwa  besaran  banjir  baik  intensitas,  frekuensi, durasi,  maupun  dampak  yang  ditimbulkan memiliki  tren  yang  meningkat    (Puslitbang Sumber  Daya  Air,  2014).  Hal  ini  dipengaruhi oleh  aktivitas  manusia  yang  menyebabkan perubahan  iklim  dan  perubahan  tutupan  lahan akibat  urbanisasi  dan  penambahan  jumlah penduduk yang signifikan. Peristiwa banjir yang berulang  sangat  merugikan  banyak  pihak
terutama  masyarakat  yang  tinggal  di  wilayah rawan banjir dan juga masyarakat yang melintasi daerah banjir yang merendam bedan jalan (BPBD Aceh Tamiang,2022).
Aceh Tamiang menjadi salah satu daerah di Provinsi  Aceh  yang  berpotensi  terjadinya bencana  banjir.  Setiap  tahunnya  Aceh  Tamiang selalu  dilanda  bencana  banjir  yang  dipicu  oleh curah  hujan  yang  sangat  tinggi.  Hal  ini mengakibatkan  volume  air  di  sungai  menjadi meluap  dan  menggenangi  daerah  sekitar  aliran sungai. Banjir yang melanda Aceh Tamiang pada November 2022 merupakan banjir terparah sejak banjir bandang di tahun 2006. Tingginya debit air luapan  sungai  bukan  hanya  merendam  ribuan bangunan, areal pertanian, perkebunan dan budi daya air, tetapi juga memutus akses jalan nasional Banda  Aceh  –  Medan  yang  mengakibatkan terhambatnya  lalu  lintas  kedua  arah  dan menimbulkan  kerugian  mencapai  200  miliar (BPBD Aceh Tamiang, 2022).
Daerah Aliran Sungai (DAS)
KruengTamiang  merupakan  salah  satu  DAS terbesar  yang  melintasi  empat  kabupaten  di Provinsi Aceh yaitu Aceh Tamiang, Aceh Timur, Gayo  Lues,  dan  sebagian  Aceh  Tenggara.  DAS Tamiang  memiliki  luas  ±  4.826,56  Km 2 atau mencapai  92%  dari  luas  total  Wilayah  Sungai (WS)  Tamiang-Langsa.  Panjang  Sungai  utama Krueng Tamiang sampai ke muara adalah sekitar 217  Km  dengan  jaringan  sungai  yang  relatif kompleks  dan  banyak  pertemuan  cabang  anak sungai  yang masih  bersifat  alami. Selanjutnya, jalur  sungai  hilir  dan  tengah  mempunyai bentangan  yang  sempit  dengan  topografi rata,  sedangkan  bagian  hulu  berada  pada belahan bukit dan gunung dengan bentangan lebar.
DAS  Krueng  Tamiang  memiliki  kapasitas infiltrasi  tanah  yang  rendah  pada  musim penghujan,  sehingga  menyebabkan  tingginya aliran  permukaan  yang  dapat  menimbulkan banjir (Azizah,2020). Banjir yang terjadi di DAS Krueng Tamiang mencapai dua sampai tiga kali dalam  satu  tahun  yang  diantaranya  banjir bandang (Menhut RI 2009; Nugroho 2017; Hadi 2018). Banjir dalam satu dekade terakhir terjadi lebih  dari  30  kali  dan  pada  tahun  2018,  telah mengalami  delapan  kejadian  banjir  (BNPB 2018).  Dengan  kondisi  dampak  banjir  tersebut dan  sumber  permasalahan  banjir  di  Aceh Tamiang,  diperlukan  strategi  dan  kebijakan
pengendalian banjir baik secara fisik ataupun non fisik.
Oleh  karena  itu,  perlu  dilakukan  pemetaan genangan banjir untuk perencanaan penggunaan lahan,  sistem  peringatan  dini,  desain  tanggap darurat, dan tindakan pengurangan risiko banjir. Pemetaan  genangan  banjir  pada  wilayah Tamiang  menggunakan  pemodelan  dinamis hidraulika  2D  dengan software Geo-HECRAS untuk  mendapatkan  titik-titik  luapan  banjir  dan luasan daerah penggenangannya.

## II. Landasan Teori

## Analisis Distribusi Curah Hujan
Perhitungan curah hujan wilayah berdasarkan data curah hujan untuk stasiun hujan yang  ada  pada  DAS  Krueng  Tamiang  dengan menggunakan Poligon Thiessen .

## Analisis Frekuensi Hujan Rencana
Analisis  frekuensi  hujan  rencana  dapat dihitung menggunakan 4 metode distribusi:
Metode distribusi Gumbel:
R T = Ṝ + K T S (1)
K T = Y±Yn S𝑛 (2)
Y T = −𝑙𝑛  [ −𝑙𝑛 (𝑇 − 1) 𝑇 ] (3)
Metode distribusi Normal:
R T = Ṝ + K T S (4)
Metode distribusi Log Normall:
Log R T = Log Ṝ + K T S (5)
Cs = 3Cv +Cv 3 (6)
Ck = Cv 8 + 6Cv 6 + 15Cv 4 + 16Cv 2 + 3      (7)
Metode distribusi Log Pearson III:
Log R T = Log R + K T S (8)

## Pengujian Kesesuaian Distribusi Frekuensi
Pengujian  kesesuaian  frekuensi  untuk mengetahui  kebenaran  dari  data  yang  diperoleh dengan uji Smirnov - Kolmogorov .

## Intensitas Hujan
Intensitas hujan adalah kedalaman air hujan persatuan  waktu  (Yulius,2014).  Perhitungan intensitas  hujan  dengan data  curah  hujan  harian maka digunakan rumus mononobe (Harto, 1993):
𝐼 = ( 𝑅24 24 ) ( 24 𝑡𝑐 ) 2/3 (9)
Waktu konsentrasi digunakan rumus Kirpich :
𝑇𝑐 = 0,06628 𝐿 0,77
𝑆 0,385
(10)

## Distribusi Hujan Jam-Jaman
Analisis distribusi curah hujan rencana jam- jaman  setiap  periode  ulang  untuk  memperoleh curah  hujan  efektif  yang  digunakan  dalam analisis  debit  banjir  rencana.  Penelitian  ini menggunakan  pola  distribusi  hujan Alternating Block Method (ABM).

## Analisis Debit Banjir Rencana
Analisis debit banjir rencana untuk menghitung debit/hidrogra banjir rencana berdasarkan  curah  hujan  rencana  setiap  periode ulang dengan metode hidrograph Satuan Sintetik (HSS) Soil Conservation Service (SCS).

## HSS SCS
Waktu puncak:
𝑡𝑝 = 0,6  × 𝑇𝑐 (11)
Debit puncak:
𝑄𝑝 = 0,208 𝐴 𝑃𝑟 (12)
Waktu  yang  diperlukan  untuk  mencapai aliran puncak:
𝑃𝑟 = tr 2 + 𝑡𝑝 (13)

## III. Metode Penelitian
Penelitian  ini  dilakukan  pada  DAS  Krueng Tamiang  dengan  luas  4826,54  km 2 atau  sama dengan 92% dari luas total wilayah Sungai (WS) Tamiang-Langsa dengan sungai utama sepanjang 217 Km. DAS Tamiang secara geografis terletak antara 3˚44’13,04”- 4˚28’5,61” LU –
97˚17’8,84”  -  98˚14’9,55”  BT.  Secara Administrasi  Wilayah  DAS  Tamiang  termasuk dalam  5  Kabupaten  di  Provinsi  Aceh  yaitu
Kabupaten  Aceh  Tamiang,  Aceh  Timur,  Aceh Tengah,  Gayo  Lues,  dan  Aceh  Tenggara  dan Kabupaten  Langkat  Provinsi  Sumatera  Utara. Peta lokasi penelitian dapat dilihat pada Gambar 1.
Gambar 1. Peta Lokasi Penelitian
Data yang digunakan dalam penelitian ini adalah data  sekunder  berupa  data  curah  hujan,  data  jenis tanah,  data  tutupan  lahan,  data  DEM,  dan  data  titik sebaran  banjir.  Setelah  data  diperoleh  dari  beberapa sumber kemudian dilakukan analisis untuk memperoleh  hasil  penelitian.  Beberapa  hal  yang menjadi  proses  dari  penelitian  ini  adalah  sebagai berikut:
1. Analisis karakteristik DAS
2. Analisis Hidrologi
3. Penentuan Koefisien Manning
4. Pemodelan Banjir
5. Pembuatan peta genangan Banjir
6. Validasi model

## IV. Hasil dan Pembahasan
Analisis Karakteristik DAS 1. Batas DAS Batas  DAS  Krueng  Tamiang  didapatkan melalui proses delineasi DAS dengan menggunakan software ArcGIS 10.5 dengan data input berupa data DEM. Setelah dilakukan proses delineasi,  diperoleh  luas  DAS  Krueng  Tamiang sebesar  4.826,56  km 2 .    Delineasi  DAS  Krueng Tamiang dapat dilihat pada Gambar 2.
Gambar 2. DAS Krueng Tamiang
2. Kemiringan Sungai
Kemiringan sungai Krueng Tamiang diperoleh dari pengerjaan menggunakan software ArcGIS  10.5  pada  menu tools  slope pada  menu 3d analyst tools . Nilai kemiringan sungai Krueng Tamiang diperoleh sebesar 0.0145.

## Analisis Hidrologi

## Distribusi Curah Hujan
Data  curah  hujan  harian  maksimum  yang digunakan dalam analisis ini adalah curah hujan tahun  2009-2021  dari  stasiun  Krueng  Langsa, Pulo Tiga,  Kp.  Selamat,  dan  Gumpang.  Berikut merupakan hasil distribusi hujan harian maksimum  DAS  Krueng  Tamiang  dapat  dilihat pada Tabel 1
Tabel  1 .  Distribusi  Hujan  Harian  Maksimum DAS Krueng Tamiang

## Analisis Frekuensi Curah Hujan
Distribusi  frekuensi  yang  digunakan  untuk menentukan besarnya curah hujan rencana, yaitu distribusi  Normal,  Log  Normal,  Gumbel,  dan Log  Pearson  III.  Parameter  statistik  diperlukan untuk  menentukan  jenis  distribusi  yang  cocok dengan  distribusi  curah  hujan  rencana  DAS Krueng Tamiang. Perhitungan parameter statistik distribusi  curah  hujan,  persyaratan  pemilihan jenis sebaran distribusi, dan curah hujan rencana dapat  dilihat  pada  Tabel  2,3,  dan  4.  Kemudian dilakukan pengujian kesesuaian distribusi menggunakan uji Smirnov-Kolmogorov .
Tabel  2 .  Perhitungan  Parameter  Statistik Distribusi Curah Hujan
Parameter Normal dan Gumbel Log Normal dan Log Pearson III Hujan rata-rata (mm) 99,38 1,98 Standar Deviasi (Sd) 34,52 0,11 Koef. Skewness (Cs) 2,96 2,27 Koef. Kurtosis (Ck) 9,70 9,70 Koef. Variasi (Cv) 0,34 0,058
Tabel  3. Persyaratan  Pemilihan  Jenis  Sebaran Distribusi
Tabel 4 . Curah Hujan Rencana
No Tr (Tahun) Pr (%) k Log Rt Rt (mm) 1 2 50 -0,338 1,942 87,440 2 5 20 0,560 2,045  110,920 3 10 10 1,277 2,127  134,091 4 25 4 2,246 2,239  173,308 5 50 2 2,989 2,324  211,023 6 100 1 3,739 2,411  257,384 Luas Poligon (km²) Curah Hujan (mm) Kr. Langsa Pulo Tiga Kp. Selamat Gum pang Thiessen (mm) 1262,91 82,.40 1985,03 758,20 No Th 1 2009 91,5 91,5 91,5 85,1 90.4 2 2010 88,9 88,9 88,9 82,7 87.9 3 2011 81,8 81,8 81,8 76,1 80.9 4 2012 70,8 70,8 70,8 65,9 70.0 5 2013 59,4 92 86,4 85,6 80.1 6 2014 84,6 98,3 92,6 91,45 91.3 7 2015 73,4 112,9 107,9 100,92 98.6 8 2016 95 148 107 100,9 109.8 9 2017 40 227 130 38,1 108.5 10 2018 43 118 120 40,8 87.0 11 2019 43 192,1 119 27,1 97.1 12 2020 105 114,1 380,2 30.3 207.9 13 2021 122,7 88,5 73,5 28,7 81.8
No Jenis Distribusi Syarat Perhitungan  Keterangan 1 Normal Cs ≈ 0 Cs = 2,96 Tidak Memenuhi Ck ≈ 3 Ck = 9,70 2 Log Normal Cs = 3Cv +Cv 3 0,17 ≠ 2,27 Tidak Memenuhi Ck = Cv 8
+ 6Cv 6 + 15Cv 4 + 16Cv 2 + 3
3,05 ≠ -9,70 3 Gumbel Cs = 1,1396 Cs = 2,96 Tidak Memenuhi Ck = 5,4002 Ck = 9,70 4 Log Pearson III Selain dari nilai diatas Memenuhi

## Distribusi Hujan Jam-jaman
Perhitungan  debit  banjir  rencana  dalam memperkirakan  besaran  debit  diperlukan  curah hujan  harian  maksimum  jam-jaman  dan  hujan efektif.  Grafik  Hyetograph  ABM  dapat  dilihat pada Gambar 3.
Gambar 3. Grafik Hyetograph ABM.

## Analisis Debit Banjir Rencana
Debit banjir diperoleh dari perkalian antara hujan  efektif  dan  unit  hidrograf  satuan  (UHS). Hidrograf banjir dapat dilihat pada Gambar 4.
Gambar 4 . Hidrograf Banjir

## Pemodelan Banjir
Batasan daerah atau mesh untuk pemodelan dimulai  dari  kecamatan  Rantau  sampai  ke  hilir Krueng Tamiang di Kecamatan Bendahara yang dapat  dilihat  pada  Gambar  4.5.  Pemodelan dilakukan  pada  daerah  model  seperti  pada Gambar  5  dan  dipilih  bentuk adaptive  mesh dengan Boundary element spacing nya digunakan  sebesar  30  m  dan  nilai  kekasaran Manning saluran sebesar 0,06. Selanjutnya garis boundary lines di bagian hulu di input data debit pada pilihan flow hydrograph sesuai periode kala ulang  dan  bagian  hilir  dianggap  pada  saat keadaan muka air laut pasang maka dipilih stage hydrograph dengan  nilai  0,4  m.  batasan  daerah pemodelan dapat dilihat pada Gambar 5.
Gambar 5. Batasan Daerah Pemodelan

## Peta Genangan Banjir
Hasil  sebaran  genangan banjir  yang diproses menggunakan software Geo-HECRAS  adalah pola  aliran  dan  genangan  banjir  pada  wilayah DAS  Krueng  Tamiang.  Peta  yang  dihasilkan berisi informasi sebaran aliran air selama periode banjir,  kedalaman  banjir,  dan  wilayah  yang tergenang oleh banjir. Peta ini dibuat berdasarkan nilai  debit  untuk periode  ulang  2, 5, 10,  25, 50, dan  100  tahun.  Luas  genangan  untuk  setiap periode ulang disajikan pada Tabel 5.
Tabel  5 .  Luas  Genangan  Banjir  Berdasarkan Periode Ulang Periode Ulang Luas Genangan (km2) Q (m3/dt) 2 43,373 985.26 5 69,389 1730.154 10 92,451 2549.389 25 138,092 4048.108 50 145,939 5601.817 100 179,836 7554.197
Untuk  melihat  perbandingan  nilai  debit dengan  luasan  genangan banjir,  disajikan  grafik hubungan luas genangan banjir dan debit. Grafik dapat dilihat pada Gambar 6.
Gambar  6 .  Grafik  Luas  Genangan  Banjir  dan Debit
Dari Gambar 6 dapat dilihat hubungan antara luas genangan  banjir  dan  debit,  semakinn  besar periode  ulang,  semakin  besar  pula  debit  yang dihasilkan  serta  semakin  luas  genangan  banjir yang akan terjadi.
Hasil peta genangan dan kedalaman banjir DAS Krueng  Tamiang  untuk  setiap  periode  ulang dapat dilihat pada Gambar 7 sampai Gambar 12.
Gambar 7. Peta Genangan dan Kedalaman Banjir DAS Krueng Tamiang Peiode Ulang 2 Tahun
Gambar 8. Peta Genangan dan Kedalaman Banjir DAS Krueng Tamiang Peiode Ulang 5 Tahun
Gambar 9. Peta Genangan dan Kedalaman Banjir DAS Krueng Tamiang Peiode Ulang 10 Tahun
Gambar  10.  Peta  Genangan  dan  Kedalaman Banjir  DAS  Krueng  Tamiang  Peiode  Ulang  25 Tahun
Gambar  11.  Peta  Genangan  dan  Kedalaman Banjir  DAS  Krueng  Tamiang  Peiode  Ulang  50
Tahun
Gambar  12.  Peta  Genangan  dan  Kedalaman Banjir DAS Krueng Tamiang Peiode Ulang 100 Tahun

## Validasi Model
Validasi  adalah  sebuah  proses  evaluasi terhadap  model  guna  mendapatkan  gambaran tentang tingkat ketidakpastian yang dimiliki oleh suatu  model.  Validasi  pada  penelitian  ini dilakukan  dengan  mencocokkan  hasil  model genangan banjir dengan data titik sebaran banjir yang diperoleh dari Dinas Pengairan sejak tahun 2003-2021.  Hasil overlay peta  genangan  banjir periode  ulang  25  tahun  dengan  titik  sebaran banjir terdapat 12 titik sebaran banjir di wilayah DAS  Krueng  Tamiang.  Dimana  10  titik  berada pada area genangan banjir dan 2 titik tidak berada di  area  genangan  banjir.  Hasil  validasi  dapat dilihat pada Gambar 13.
Gambar 13 . Peta Validasi Titik Banjir Periode Ulang 25 Tahun DAS Krueng Tamiang IV. Kesimpulan dan Saran

## A. Kesimpulan
Berdasarkan  hasil  analisis  didapatkan kesimpulan sebagai berikut:
1. Dari hasil pemodelan diperoleh luas genangan banjir  untuk    setiap  periode  ulang  sebagai berikut: L 2th = 43,373km 2 , L 5th = 69,389 km 2 , L 10th = 92,451 km 2 , L 25th = 138,092 km 2 , L 50th = 145,939 km 2 , dan L 100th = 179,836 km 2 serta diperoleh juga ketinggian banjir untuk setiap periode  ulang  adalah  0  m  sampai  dengan  > 2m.
2. Validasi  dilakukan  secara  kualitatif,  yaitu dengan mengindentifikasi kecocokan wilayah genangan  banjir  yang  dihasilkan  oleh  model dengan  lokasi  genangan  banjir  berdasarkan data survey lapangan.
B. Saran
Adapun saran yang diberikan oleh penulis terkait penelitian yang telah dilakukan, sebagai berikut:
1. Perlu  dilakukan  penyempurnaan  data  yang akan  digunakan  dalam  pemodelan,  seperti data  DEM  dan  data  curah  hujan  untuk meningkatkan akurasi hasil pemodelan.
2. diharapkan hasil pemodelan hidrodinamis 2D ini  dapat  diintegrasikan  dalam  sistem peringatan  dini  banir  untuk  memberikan informasi yang lebih tepat waktu dan akurat
3. Diharapkan penelitian tugas akhir ini dapat menjadi referensi bagi penelitian yang mempunyai topik yang serupa.

## Daftar Pustaka
[1] Azizah,  C.,  Pawitan,  H.,  Dasanto,  B.D., Ridwansyah,  I.,  Taufik,  M.  (2020).  Sifat  Fisik Tanah  dan  Hubungannya  dengan  Kapasitas Infiltrasi DAS Tamiang. Jurnal Tanah dan Iklim Vol. 43 No. 2, Desember 2019: 167-173
[2] BNPB (Badan Nasional Penanggulangan Bencana).  2018.  Data  dan  Informasi  Bencana Indonesia.
[3] Chow, V.T., Maidment, D.R., and Mays., L.W. Applied  Hydrology.  Singapore:  McGraw-Hill Book Company, 1998.
[4] Harto,  Sri.  (1993). Analisis  Hidrologi .  Jakarta: PT Gramedia Pustaka Utama.
[5] Menhut RI. 2009. Keputusan Menteri Kehutanan Republik  Indonesia  Nomor:  Sk.  328/Menhut- Ii/2009  Tentang  Penetapan  Daerah  Aliran Sungai (DAS) Prioritas Dalam Rangka Rencana Pembangunan Jangka Menengah (RPJM) Tahun 2010-2014.
[6] Puslitbang  Sumber  Daya  Air.  2014.  Penelitian Prakiraan  dan  Pengendalian  Kekeringan  Serta Pengembangan Peta Risiko Banjir dan
Kekeringan. Kementerian Pekerjaan Umum
[7] Soewarno.  (1995). Hidrologi,  Aplikasi  Metode Statistik Untuk Analisa Data . Bandung: Nova.
[8] Triatmodjo,  B.  (2008). Hidrologi  Terapan . Yogyakarta: Beta Offset.
[9] Yulius,  E. (2014).  Analisa  Curah  Hujan  Dalam Membuat  Kurva Intensity  Duration Frequency (IDF) Pada DAS Bekasi. In Jurnal Bentang (Vol. 2, Issue 1).