## KESESUAIAN HIDROGRAF SATUAN SINTETIK TERHADAP HIDROGRAF SATUAN TERUKUR(STUDI KASUS SUB SUB DAERAH ALIRAN SUNGAI PEDINDANG BAGIAN TENGAH)

## Gustama
Mahasiswa Jurusan Teknik Sipil Univesitas Bangka Belitung Email : tama.tekniksipil@yahoo.com

## Fadillah Sabri
Staf Pengajar Jurusan Teknik Sipil Universitas Bangka Belitung Email : sabrifadillah@yahoo.com

## Donny Fransiskus Manalu
Staf Pengajar Jurusan Teknik Sipil Universitas Bangka Belitung Email : donny_fm@yahoo.com

## ABSTRACT
Flooding can be caused by the inability of the river to accommodate the flow  of water that passes  through  the  river  so  that  the  water  overflows.  A  widely  used  method  for  analyzing river  flow  for  flood  forecasts  is  hydrograph  unit.  The  hydrograph  unit  is  a  direct  runoff hydrograph that can be created when there are AWLR record data, debit measurements and rainfall data. Synthetic Unit Hydrograph (SUH) is a unit hydrograph derived based on river data  in  the  same  watershed  or  nearby  watershed  but  has  the  same  characteristics,  ie  HSS Gama I, HSS Nakayasu, Limasan HSS, HSS Snyder and HSS SCS. Of the two hydrographs, there  will  be  suitability  of  the  hydrograph  form  that  is  going  to  be  made. Sub  territory  of  Pedindang    River  Basin  has  four  flood  incidents,  namely,  date  23-24 February  2016;  March  2-3,  2016;  March  3-4,  2016;  and  date  5-6  March  2016.  In  the analysis  of  each  flood  event,  the  peak  discharge  of  synthetic  unit  hydrograph  is  very different  from  the  peak  discharge  of  the  measured  unit  hydrograph.  The  average  peak discharge  of  synthetic  unit  hydrograph  occurs  in  the  range  of  2  or  3  hours,  while  the measured unit hydrograph of Pedindang River occurs in the range of 7 or 8 hours. In four flood  events  it  is  stated  that,  HSS  Gama  I  approaches  RMSE  value  (validation  <10%)  to HST form of Pedindang River with value: RMSE incidence I (23,601%); RMSE incidence II (16.315%); RMSE incidence III (50,400%); RMSE incidence IV (22.322%). With this result, it is stated that there is no synthetic unit hydrograph model that has compatibility with the measured unit hydrograph of Pedindang River.
Key Word: flood hydrograph, synthetic unit hydrograph, RMSE

## PENDAHULUAN
Daerah Aliran Sungai (DAS) adalah suatu  wilayah  daratan  yang  secara topografik dibatasi oleh punggung- punggung  gunung  yang  menampung  dan menyimpan  air  hujan  untuk  kemudian menyalurkan ke laut melalui sungai utama (Asdak,  2014).  Daerah  Aliran  Sungai (DAS)  dapat  digunakan  sebagai  batas kajian  suatu  penelitian,  mengingat  DAS merupakan suatu sistem yang kompleks. Ilmu  yang  mempelajari  air  dalam  segala bentuk  (cairan,  gas,  padat)  pada,  dalam dan  di  atas  permukaan  tanah.  Termasuk didalamnya  adalah  penyebaran,  daur  dan perilakunya, sifat-sifat fisika dan kimianya.  Serta  hubungannya  dengan unsur-unsur  hidup  dalam  air  itu  sendiri disebut  dengan  hidrologi  (Asdak,  2014). Apabila  dilihat  dari  pernyataan  ini, hidrologi  adalah  ilmu    yang  tepat mengkaji masalah-masalah yang berkaitan dengan  air,  salah  satunya  masalah  banjir. Oleh  karena  hal  ini,  debit  aliran  sungai merupakan hal  yang sangat  penting untuk dikaji  dalam  prakiraan  banjir.  Metode yang banyak digunakan untuk
menganalisis debit sungai untuk prakiraan banjir adalah hidrograf satuan. Hidrograf  satuan  adalah  hidrograf limpasan langsung ( direct runoff hydrograf )  yang  dihasilkan  oleh  hujan efektif yang terjadi merata di seluruh DAS dengan intensitas tetap dalam suatu satuan tertentu,  yaitu  untuk  satuan  hujan  1 mm/jam.  Hirograf  satuan  dapat  dibuat apabila tersedia data seperti  data rekaman AWLR ( Automatic Water Level Recorder ),  data  pengukuran  debit  juga data hujan.
Hidrograf  Satuan  Sintetik  (HSS) adalah  hidrograf  satuan  yang  diturunkan berdasarkan  data  sungai  pada  DAS  yang sama  atau  DAS  terdekat  tetapi  memiliki karakteristik yang sama. Dalam penelitian dilakukan  analisis  pemodelan  hidrograf satuan  sintetik  yang  dibuat  dengan menggunakan  karakteristik  DAS  untuk mencari  kesesuiaan  model  hidrograf satuan  manakah  yang  sesuai  dengan hidrograf satuan terukur  di  Sub Sub DAS Pedindang, model hidrograf satuan sintetik dibatasi dengan hanya menggunakan  HSS  Snyder,  Nakayasu, Gama I, Limantara, SCS.
Sungai  Pedindang  merupakan  Sub Sub  DAS  Baturusa.  Sungai  Pedindang yang  dibatasi  oleh  tinjauan  lokasi  yaitu pada outlet yang  terletak  pada  stasiun pengukuran  tinggi  muka  air  mempunyai luas  DAS  35,14  km 2 dengan  panjang sungai  15,93  km  yang  membentang  dari Desa  Mangkol  sampai  Desa  Parit  Lalang yang  terdapat  alat  AWLR.  Sungai  inilah yang  menjadi  bahan  untuk  melakukan pengamatan  karena  pada  kawasan  ini tidak  tersedia  pengukuran-pengukuran langsung mengenai hidrograf banjirnya.

## LANDASAN TEORI A. Siklus Hidrologi
Siklus  hidrologi  merupakan  proses kontinyu dimana air bergerak dari bumi ke atmosfer  dan  kemudian  kembali  kebumi lagi  (Chow,  1998  dalam  Triatmodjo, (2014)).

## B. Daerah Aliran Sungai (DAS)
Daerah  Aliran  Sungai  (DAS)  adalah wilayah  daratan  yang  secara  topografik
dibatasi oleh punggung-punggung gunung yang  menampung  dan  menyimpan  air hujan  kemudian  menyalurkannya  ke  laut melalui sungai utama.

## C. Hujan
Presipitasi adalah turunya air dari atmosfer  ke  permukaan  bumi,  yang  bisa berupa  hujan,  hujan  salju,  kabut,  embun dan hujan es (Triatmodjo, 2014). Mengingat bahwa di daerah Indonesia dan lokasi  penelitian  adalah  daerah  iklim tropis,  maka  presipitasi  dalam  konteks daerah tropis sama dengan curah hujan.

## D. Infiltrasi
Infiltrasi adalah proses aliran air (umumnya  berasal  dari  curah  hujan) masuk kedalam tanah. Proses mengalirnya air hujan ke dalam tanah disebabkan oleh tarikan  gaya  gravitasi  dan  gaya  kapiler tanah (Asdak, 2014).
1. Pengukuran Infiltrasi
Ada tiga cara untuk menentukan besarnya infiltrasi  menurut  Knapp,  (1978) dalam Asdak, (2014):
a. Simulasi di laboratorium.
b. Menggunakan alat infiltrometer.
c. Teknik  pemisahan  hidrograf  aliran dari data aliran air hujan.
2. Indeks Infiltrasi (Φ indeks )
Indeks infiltrasi adalah laju rerata kehilangan air karena infiltrasi, sedemikian  sehingga  volume  air  hujan yang  lebih  dari  laju  tersebut  adalah  sama dengan  aliran  permukaan  (Triatmodjo, 2014),  dapat  dihitung  dengan  persamaan berikut:
Tr Q P Tr F Indeks     ......................(1)
dengan:
F =   infiltrasi total (mm) P =   hujan total (mm)
Q =   aliran total (mm) Tr =   waktu terjadinya hujan (jam)

## E. Hidrometri
1. Pengukuran Kecepatan Aliran Air Pengukuran kecepatan arus secara langsung dapat dilakukan dengan menggunakan  pelampung,  yaitu  dengan mengukur  selang  waktu  yang  diperlukan untuk  oleh  pelampung  untuk  menempuh suatu  jarak  tertentu,  dan  dapat  dihitung dengan persamaan berikut:
t L V  ..................................................(2)
dengan:
V =   kecepatan (m/det)
L =  jarak (m) t =   waktu (detik)
2. Hitungan Debit Aliran
Soewarno (2013), memberikan teori untuk menghitung  debit.  Debit  sungai  dapat dihitung  dengan  cara  mengukur  luas penampang basah dan kecepatan alirannya, dengan persamaan berikut:
    v a Q . .........................................(3)
p V k v .  .............................................(4)
dengan:
Q =   debit sungai total (m 3 /det)
a =   luas bagian penampang basah v =   kecepatan  aliran  rata-rata  pada bagian penampang basah (m/det)
k =   faktor  koreksi  (0,85  menurut
Triatmodjo, 2014)
V p
=   kecepatan  lintasan  pelampung (m 3 /det)
3. Persamaan Lengkung Debit
Soewarno  (2013)  menyatakan  kurva lengkung  debit  dapat  dihitung  dengan metode  logaritmik,  dengan  persamaan umum berikut:
  n H H K Q 0   ..................................(5) dengan: Q = debit (m 3 /d) H = tinggi muka air (m) H 0 = tinggi muka air pada saat aliran nol
(m ) (saat Q = 0)
K,n = bilangan konstanta Semua data pengukuran
debit
digambarkan pada kertas grafik logaritmik.  Pilih  harga  debit  Q 1 ,  Q 2 dan Q 3 dari  kurva  lengkung  debitnya,  dan harga  Q 2 ditentukan  dari  persamaan  rata- rata geometris sebagai berikut:
  5 , 0 3 1 2 Q Q Q   .................................(6)
Penentuan  harga  H 0 dapat  ditentukan dengan  tiga  cara,  yaitu  coba-ulang  ( tial and  eror ),  aritmatik  dan  grafis.  Untuk penelitian  ini  harga  H 0 ditentukan  dengan metode  aritmatik  dengan  persamaan sebagai berikut:
2 3 1 2 2 3 1 0 2 H H H
H H H H     .......................(7)
Untuk  mencari K dan  n  dapat  dibantu dengan dua persamaan dibawah ini:
 
       0 log x n k m y ..............(8)             0 log 2 x n k x xy ....(9)
dengan:
   y =  jumlah nilai log Q (m 3 /d)    x =  jumlah nilai log (H-H 0 ) (m)  
 2 x =  jumlah nilai kuadrat dari X
 (xy) =  jumlah nilai X dikalikan Y
m
=  jumlah pengamatan

## F. Hidrograf
Hidrograf adalah kurva  yang memberikan hubungan  antara  parameter  aliran  dan waktu  (Triatmodjo,  2014).  Parameter tersebut  bisa  berupa  kedalaman  aliran (elevasi)  atau  debit  aliran,  sehingga terdapat  dua  macam  hidrograf  yaitu hidrograf muka air dan hidrograf debit. 1. Hidrograf Satuan Metode Collins
Hidrograf satuan ditarifkan sebagai hidrograf limpasan langsung ( direct runoff hydrograf )  yang  dihasilkan  oleh  hujan efektif yang terjadi merata di seluruh DAS dengan  intensitas  tetap  dalam  satu  satuan tertentu.
2. Hidrograf Satuan Sintetik
Hidrograf satuan sintetik adalah hidrograf satuan  yang  diturunkan  berdasarkan  data sungai  pada  DAS  yang  sama  atau  DAS terdekat tetapi memiliki karakteristik yang sama  (Kamiana,  2012).  Indarto  (2016) menyatakan  idealnya  hidrograf  satuan diturunkan  dari  data  hujan  dan  AWLR, tetapi  pada  kebanyakan  DAS  jarang  data yang lengkap.
a) Hidrograf Satuan Sintetik Nakayasu Nakayasu (1950) telah menyelidiki hidrograf satuan di Jepang dan memberikan seperangkat persamaan untuk membentuk suatu hidrograf satuan sebagai berikut:
         
3 , 0 3 , 0 Re 6 , 3 1 T Tp A Q p ...................(10)
r g p T t T 8 , 0   ...................................(11)
L t g 058 , 0 4 , 0   L > 15 km  …......(12)
7 , 0 21 , 0 L t g  ... L < 15 km  ..............(13)
g t T   3 , 0 ......................................(14)
g g r t t T sampai 5 , 0  .......................(15)
dengan:
Q p = Debit puncak banjir
A = Luas DAS (km 2 )
R e = Curah hujan efektif (1 mm)
T p = Waktu  dari  permulaan  banjir sampai puncak hidrograf (jam)
T 0,3 = Waktu  dari  puncak  banjir  sampai
0,3 kali debit puncak (jam)
t g = Waktu konsentrasi (jam)
T r = Satuan  waktu  dari  curah  hujan
(jam)
α = Koefisien karakteristik DAS
biasanya diambil 2
L = Panjang sungai utama (km)
Bentuk  hidrograf  satuan  diberikan  oleh
persamaan berikut:
4 , 2          p p t T t Q
Q ................................(16)
Pada kurva ( T p < t < T p +T 0,3 ) 3 , 0 / ) ( 3 , 0 T T t p r p Q Q    ........................(17)
Pada  kurva  turun  ( T p +T 0,3 < t  <  T p +T 0,3 + 1,5 T 0,3 )       3 , 0 3 , 0 5 , 1 / 5 , 0 ) ( 3 , 0 T T T t p r p Q Q     ...........(18)
Pada kurva turun ( t > T p +T 0,3 + 1,5 T 0,3 )       3 , 0 3 , 0 2 / 5 , 0 ) ( 3 , 0 T T T t p r p Q Q     ...........(19)
Sumber: Bambang Triadmodjo, 2014
Gambar 1 Hidrograf satuan sintetik
Nakayasu b) Hidrograf Satuan Sintetik Snyder Snyder (1938) mendapatkan dan mengembangkan hidrograf satuan DAS di Amerika  Serikat  yang  berukuran  30 sampai 30.000 km 2 dengan
menghubungkan  unsur-unsur  hidrograf satuan  dengan  karanteristik  DAS  akibat hujan 1 cm. Persamaan berikut:
  3 , 0 c t p L L C t  ...................................(20)
p p p t A C Q  .........................................(21) 8 3 p t T   ..........................................(22)
5 , 5 p D t t  .............................................(23)
Apabila durasi  hujan efektif t p tidak sama dengan durasi standar t D , maka:   D r p pR t t t t    25 , 0 ........................(24)
pR p p pR t t Q Q  .....................................(25)
dengan:
t D = Durasi  standar  dari  hujan  efektif (jam) t r
= Durasi hujan efektif (jam)
t p = Waktu dari titik berat durasi hujan efeftif t D ke puncak hidrograf satuan (jam)
t pR = Waktu dari titik berat durasi hujan t r ke puncak hidrograf satuan (jam)
T = Waktu  dasar  hidrograf  satuan (jam) Q p = Debit puncak untuk durasi t D Q pR = Debit puncak untuk durasi t r L = Panjang  sungai  utama  terhadap titik kontrol yang ditinjau (km)
R0= 1mm Tg 0,8 Tr tp t 1,5 t 0,3 0,3 Qp 0,3 Qp 0,3² Qp Tr t (jam)
L c = Jarak  antara  titik  kontrol  ke  titik
yang terdekat dengan titik
berat DAS (km)
A = Luas DAS (km 2 )
C t = Koefisien yang tergantung
kemiringan DAS (0,75 sampai 3)
C p = Koefisien  yang  tergantung  pada karakteristik DAS (0,9 sampai 1,4) Dengan menggunakan rumus-rumus tersebut  di  atas  dapat  digambarkan hidrograf  satuan.  Untuk  memudahkan penggambaran,  beberapa  rumus  diberikan sebagai berikut:
08 , 1 08 , 1 50 23 , 0 pR Q A W  .................................(26) 08 , 1 08 , 1 75 13 , 0 pR Q A
W  ..................................(27)
Dengan W 50 dan W 75 adalah  lebar  unit hidrograf  pada  debit  50%  dan  75%  dari debit puncak, yang dinyatakan dalam jam.
c) Hidrograf Satuan Gama I
Hidrograf satuan sintetis Gama I
dikembangkan  oleh  Sri  Harto  (1993, 2000)  berdasar  perilaku  hidrologis  30 DAS di Pulau Jawa. Meskipun diturunkan dari  DAS  di  Pulau  Jawa,  ternyata hidrograf  satuan  sintetis  Gama  I  juga berfungsi  baik  untuk  berbagai  daerah  lain di  Indonesia.  HSS  Gama  I  terdiri  dari empat  variabel  pokok,  yaitu  waktu  naik ( time  of  rise  -  TR ),  debit  puncak  ( Q p ), waktu  dasar  ( TB )  dan  sisi  resesi  yang ditentukan oleh nilai koefisien tampungan ( K ) yang mengikuti persamaan berikut:
K t p t e Q
Q /   ......................................(28)
dengan:
Q t = Debit pada jam ke t ( m 3 /d )
Q p = Debit puncak ( m 3 /d )
t = Waktu  dari  saat  terjadinya  debit
pucak ( jam )
K = Koefisien tampungan ( jam )
Untuk  menghitung  Hidrograf  Gama  I dapat  digunakan  persamaan-persamaan berikut:
Waktu puncak HSS Gama I ( TR )
2775 , 1 0665 , 1 3 100 43 , 0          SIM SF L TR (29)
1) Debit puncak banjir ( QP ) 2381 , 0 4008 , 0 5886 , 0 1836 , 0 JN TR A
QP   (30) 2) Waktu dasar ( TB ) 2574 , 0 7344 , 0 0986 , 0 1457 , 0 4132 , 27 RUA SN S TR TB   (31)
3) Koefisien resesi ( K ) 0452 , 0 0897 , 1 1446 , 0 1798 , 0 5617 , 0 D SF S A K    (32)
4) Aliran dasar ( QB ) 8430 , 0 6444 , 0 4715 , 0 D A QB  (33)
dengan:
A = Luas DAS ( km 2 )
L = Panjang sungai utama ( km ) S
= Kemiringan dasar sungai
SF = Faktor sumber, perbandingan antara jumlah panjang pangsa  sungai  tingkat  satu  dengan jumlah panjang sungai tingkat SN = Frekuensi  sumber,  perbandingan antara jumlah pangsa sungai  tingkat satu  dengan  jumlah  pangsa  sungai semua tingkat
WF = Faktor  lebar,  perbandingan  antara lebar  DAS  yang  diukur  di  titik sungai  berjarak 0,75 L dengan lebar DAS  yang  diukur  di  sungai  yang berjarak 0,25 L dari stasiun hidrometri (Gambar 2.16).
JN = Jumlah pertemuan sungai
SIM = Faktor  simetri,  hasil  kali  antara faktor  lebar  ( WF )  dengan  luas  DAS sebelah hulu ( RUA )
RUA = Luas DAS sebelah
hulu,
perbandingan antara luas DAS yang diukur  di  hulu  garis  yang  ditarik tegak  lurus  garis  hubungan  antara stasiun hidrometri dengan titik yang paling dekat dengan titik berat DAS, melalui titik tersebut (Gambar 2.16).
D = Kerapatan  jaringan  kuras,  jumlah panjang  sungai  semua  tingkat  tiap satuan luas.
d) Hidrograf Satuan Sintetik Limantara
Hidrograf  satuan  sintetik  limantara,  yang asalnya  dari  Indonesia,  ditemukan  oleh Lily  Montarcih  Limantara  pada  tahun 2006.  Lokasi  penelitian  di  sebagian Indonesia  dianggap  mewakili  :  Jawa  (6 DAS, 67 Sub DAS), Bali (2 DAS, 13 Sub DAS), Lombok (1 DAS, 5 Sub DAS) dan Kalimantan  Timur  (1  DAS,  9  Sub  DAS). Untuk persamaan- persamaan yang digunakan  Hidrograf  Satuan  Sintetik
Limantara berikut:
1) Debit puncak 168 , 0 . 131 , 0 . 356 , 0 . 497 , 0 . 451 , 0 . 042 , 0 n S Lc L A Qp   ...(34)
dengan Qp = debit  puncak  banjir  hidrograf satuan ( m 3 /det ) A = luas DAS ( km 2 )
L = panjang sungai utama ( km )
Lc = kemiringan sungai utama
n = koefisien kekasaran DAS
2) Persamaan kurva naik     107 , 1
/ . p n T t Qp Q  .......................(35) dengan
Q n = debit  pada  persamaan  kurva  naik ( m 3 /det )
Qp = debit  puncak  banjir  hidrograf satuan ( m 3 /det ) t = waktu hidrograf ( jam )
Tp = waktu  naik  hidrograf  atau  waktu mencapai puncak hidrograf ( jam )

## 3) Persamaan kurva turun
) ( 175 , 0 10 . t Tp t Qp Q
  ...........................(36) dengan
Q t = debit  pada  persamaan  kurva  turun ( m3/det )
Qp = debit  puncak  hidrograf  satuan ( m3/det )
Tp = waktu  naik  hidrograf  atau  waktu mencapai puncak hidrograf ( jam )
t = waktu hidrograf ( jam )
4) Perkiraan waktu puncak banjir (Tp)
Untuk  memperkirakan  waktu  puncak banjir  (Tp)  bisa  dipakai  rumus  seperti pada Nakayasu berikut:
tr tg Tp 8 , 0   ..................................(37) dengan
Tp = tengangan  waktu  ( time  lag )  dari permulaan  hujan  sampai  puncak banjir ( jam )
tg = waktu  konsentrasi  hujan  ( jam ), untuk menentukannya dapat digunakan persamaan berikut : Jika, L ≥ 15 km, maka L tg 058 , 0 40 , 0   ..(38)
L < 15 km, maka 7 , 0 21 , 0 L tg  ...... (39)
g g r t sampai 0,5t T  ..............(40) e) Hidrograf Satuan Sintetik SCS SCS menggunakan hidrograf tak
berdimensi  yang  dikembangkan  dari
analisis  sejumlah  besar  hidrograf  satuan dari  data  lapangan  dengan  berbagai ukuran  DAS  dan  lokasi  berbeda.  Ordinat hidrograf  satuan  untuk  periode  waktu berbeda dapat diperoleh dari tabel berikut, dengan nilai (Gupta, 1989):
p p T A Q 208 , 0  ....................................(41)
g r p t t
T   2 .........................................(42)
385 , 0 2 1000 87 , 0
         S L t c ...........................(43) c g t t   6 , 0 ........................................(44)
c r t t   133 , 0 .....................................(45) dengan
Q p = debit puncak (m 3 /det)
A = luas DAS (km 2 ) T p = waktu puncak (jam)
S = kemiringan sungai utama t g
= lag time (jam) t r
= durasi hujan efektif (jam)
Tabel 1 Hidrograf satuan metode SCS t/p r Q/Q p t/p r Q/Q p t/p r Q/Q p
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,015 0,075 0,16 0,28 0,43 0,60 0,77 0,89 0,97 1,0 1,1 1,2 1,3 1,4 1,5 1,6 1,8 2,0 2,2 1,0 0,98 0,92 0,84 0,75 0,66 0,56 0,42 0,32 0,24 2,4 2,6 2,8 3,0 3.5 4,0 4,5 5,0 0,18 0,13 0,098 0,075 0,036 0,018 0,009 0,004 0

## Sumber : Triadmodjo, 2014

## G. Kalibrasi Model Hidrograf Debit
RMSE ( root mean square erorr ) merupakan    metode  alternatif  untuk
mengevaluasi  teknik  peramalan  yang digunakan untuk mengukur tingkat akurasi  hasil  perkiraaan  suatu  model. Parameter model yang dipilih bila dengan uji  statistik  diperoleh  nilai  RMSE  ( root mean square erorr ) yang kurang dari 10% (Soewarno,  2013).  RMSE  dapat  dihitung dengan rumus:
         n n Qo Qe RMSE 1
2 .................(46) dengan : RMSE =  Root Mean Square Error Qe = Debit estimasi (m 3 /det)
Qo = Debit observasi (m 3 /det) n = Jumlah pengamatan

## H. Korelasi
Korelasi  adalah  istilah  statistik  yang menyatakan derajat hubungan linier antara dua  variabel  atau  lebih,  yang  ditemukan oleh  Karl  Pearson  pada  awal  1900  dan dikenal  dengan Korelasi Pearson  Product Moment (PPM), (Susanti, 2014).
Untuk  data  kualitatif  dalam  mengukur kuatnya  hubungan  disebut Contigency Coefficient (koefisien  bersyarat,  Cc).
Adapun  batas  ketentuan  Cc  sebagai berikut:
a. Jika  nilai  perbandingan  Cc  dengan batas < 0,50 maka hubungan lemah
b. Jika nilai Cc antara 0,50 dan 0,75 maka hunbungan sedang/ cukup
c. Jika nilai Cc antara 0,75 dan 0,90 maka hubungan kuat
d. Jika  nilai  Cc  anata  0,90  dan  1  maka hubungan sangat kuat
e. Jika  nilai  Cc  =  1  maka  hubungan sempurna

## METODE PENELITIAN

## A. Lokasi Penelitian
Secara geografis Sub Sub DAS Pedindang berada pada 106 o 00’BT sampai 106 o 10’BT  dan  2 o 05’LS  sampai  2 o 15’LS.
Sub  DAS  Pedindang  mempunyai  panjang 15,93 km yang membentang dan memotong  dua  daerah  administrasi,  yaitu hulu  sungai  berada  di  Desa  Mangkol Kabupaten  Bangka  Tengah  dan  hilir sungai  berada  di  Stasiun  Pengukuran
Tinggi Muka Air (AWLR) Kota
Pangkalpinang.  Sub  Sub  DAS  Pedindang mempunyai  luas  DAS  yaitu  dengan  luas 35,14 km 2 .
Sumber : Data diolah, 2018
Gambar 2 Bentuk Sub Sub DAS

## Pedindang
Mata  air  Sungai  Pedindang  berada  di Gunung  Mangkol  (397,5  mDPL).  Dari Gunung  Mangkol,  Sungai  Pedindang mengalir  menuju  timur  laut  melewati Kolong  Pedindang  di  Desa  Parit  Lalang. Kolong  Pedindang  adalah  kolong  yang dipergunakan  sebagai  salah  satu  sumber air  baku  untuk  Kota  Pangkalpinang. Sungai  Pedindang  bermuara  di  Kolong Pasar Ikan yang berada di Kota
Pangkalpinang  sebelum  bertemu  dengan Sungai Rangkui.

## B. Alat dan Bahan
Alat  yang  digunakan  dalam  penelitian  ini berupa  seperangkat  komputer  dengan software, AutoCAD, meteran, GPS, Stopwatch, Microsoft Excel  2013 , ArcGIS 10.1 .  Sedangkan  data  yang  diperlukan dalam penelitian ini meliputi:
1) Data  karakteristik  Sub  Sub  DAS Pedindang.
2) Data  hujan  jam-jaman  Stasiun  BMKG dan  Stasiun  Cuaca  Univesitas  Bangka Belitung periode 1 Februari 2016 s.d. 7
Maret 2016. 3) Data  tinggi  muka  air  jam-jaman Kementerian SDA Wil. Sungai Sumatera VIII periode 1 Februari 2016
s.d. 7 Maret 2016.
4) Data  primer  Sungai  Pedindang  yaitu, kecepatan, luas penampang, dan debit.

## C. Diagraf Alir
Diagram  alir  ( flowchart )  penelitian  yang digunakan  dalam  penelitian  ini  tersaji seperti pada Gambar 3.
Gambar 3 Diagram alir ( flowchart )

## penelitian

## ANALISIS DAN PEMBAHASAN A. Analisis Hidrometri

## 1. Pengukuran Tampang Sungai

## Sumber : Hasil analisis, 2017
Gambar 4 Bentuk Sub Sub DAS Pedindang
Pengukuran tampang sungai ini dilakukan bersamaan dengan mengukur tinggi  muka air.  Dari  pengukuran  tinggi  muka  air tersebut akan mendapatkan luas
penampang  basah  yang  akan  digunakan untuk menghitung debit aliran.
2. Perhitungan Debit Aliran Air
Perhitungan debit aliran air dapat dihitung setelah  memperoleh  data  tinggi  muka  air, luas tampang basah serta kecepatan aliran air.  Debit  dihitung  berdasarkan  data tersebut  yang  diukur  1  hari  sebanyak  3
kali  yaitu  waktu  pagi,  siang  dan  sore selama 30 hari.
Tabel 2 Peritungan kecepatan aliran air
Sumber : Hasil analisis, 2017
Tabel 3 Perhitungan debit aliran air
Sumber : Hasil analisis, 2017
Adapun contoh rekapitulasi data pengukuran debit aliran air terlampir pada Lampiran-5.

## B. Persamaan Lengkung Debit
Dari perhitungan debit diperoleh sebanyak
90 data pengukuran debit. Tahap perhitungan  dalam  pembuatan  lengkung debit sebagai berikut:
1. Menentukan harga H 0 Menetukan harga H 0 dengan
menggunakan  cara  aritmatik,  berdasarkan persamaan (6) dan (7) sebagai berikut: H 1 = 0,6347 m dan Q 1 = 1,8966 m 3 /det H 3 = 1,0245 m dan Q 3 = 4,4046 m 3 /det Berdasarkan persamaan (6), maka:
Q 2 = 0,5 4,4046) (1,8966  = 2,8903 m 3 /det dan  untuk  mencari  H 2 dengan  cara interpolasi linier sebagai
mendapatkan nilai H 2 = 0,7891 m. Setelah  didapatkan  nilai  H 2 ,  maka berdasarkan  persamaan  (7)  didapatkan  H 0 sebesar :
Datum ±0.00 MSL - 6.0 - 3.00 - 2.00 - 1.00 ± 0.00 + 1.00 + 2.00 ? - 4.486 - 0.664 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11 P12 P15 P1 P16 Jarak (m) Elevasi (m) STA ± 0.00 - 2.581 - 4.486 - 0.091 + 6.0
H 0 =           7891 0 2 0245 1 6347 0 7891 0 0245 1 6347 0 2 , , , , , ,      H 0 = 0,3399 m
Hasil  perhitungan  secara  keseluruhan telampir  pada  Lampiran-3,  dan  diperoleh hasil sebagai berikut:
Σ(x)  : -41,5047 Σ(y)  : 30,8991 Σ(xy) : -13,5374 Σ(x 2 ) : 19,8763
m : 90
2. Menetukan nilai K dan n Setelah mendapatkan nilai tersebut, langkah  selanjutnya  menentukan  nilai K dan n , berdasarkan persamaan (8) dan (9) maka diperoleh sebagai berikut:
30,8991 – 90 log K + 41,5047n  = 0 .. (47)
-13,5374 + 41,5047log K – 19,8763n = 0
(48)
Dengan  persamaan  (47)  dan  persamaan (48),  untuk  memperoleh  nilai  K  dan  n digunakan metode eliminasi dan substitusi dan mendapatkan nilai K = 6,1597 dan n = 0,9676. Dengan nilai H 0 , K dan n tersebut maka diperoleh persamaan lengkung debit Sungai  Pedindang  untuk  tanggal  2 Agustus  2017  sampai  dengan  31  Agustus 2017  menurut  persamaan  (5)  adalah Q  =
6,1597 (H – 0,3399) 0,9676 .

## C. Data Hujan
Data  Data  hujan  berupa  data  curah  hujan jam-jaman yang dilakukan analisis korelasi  antara  data  hujan  per  3  jam-an pada  Stasiun  Iklim  Univesitas  Bangka Belitung  terhadap  Stasiun  Meteorologi Kelas I Depati Amir Kota Pangkalpinang. Analisis  ini  dilakukan  untuk  menetukan seberapa kuat hubungan antara kedua data hujan tersebut.

## Sumber : Hasil analisis, 2017
Gambar 5 Grafik korelasi data hujan antar
stasiun
Dari  grafik  curah  hujan  per  3  jam-an kemudian dianalisis menggunkaan
metode  korelasi.  Dengan  metode  korelasi menggunakan  microsoft  excel  diperoleh hasil sebagai berikut.
Tabel 4 Nilai korelasi data hujan antar
stasiun Data Hujan UBB Data Hujan BMKG Data Hujan UBB 1 Data Hujan BMKG 0,837965 1

## Sumber : Hasil analisis, 2018
Dari  Tabel  4  diperoleh  nilai  korelasi antara data hujan Stasiun Iklim Univesitas Bangka  Belitung  terhadap  data  hujan Stasiun  Meteorologi  Klas  I  Depati  Amir Kota  Pangkalpinang  sebesar  0,837965, dapat  ditarik  kesimpulan  bahwa,    dari hasil tersebut  menunjukan data hujan dari kedua stasiun tersebut mempunyai korelasi (hubungan) yang kuat.

## D. Penentuan Kejadian Banjir
Dari  data  debit  dan  hujan  Sungai Pedindang  didapatkan  empat  kejadian banjir  yaitu  tanggal  23-24  Februari  2016, 2-3  Maret  2016,  3-4  Maret  2016  dan  5-6
0 20 40 60 80 100 120 140 160 180 200 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 C u r ah H u jan ( m m ) Tanggal Grafik Korelasi Data Hujan Data Hujan Univesitas Bangka Belitung Data Hujan BMKG Kota Pangkalpinang
Maret  2017.  Adapun  data  penentuan kejadian  banjir  terlampir  pada  Lampiran 3.

## Sumber : Hasil analisis, 2017
Gambar 6 kejadian banjir tanggal 23-24 Februari 2016
Sumber : Hasil analisis, 2017
Gambar 7 kejadian banjir tanggal 2-3 Maret 2016
Sumber : Hasil analisis, 2017
Gambar 8 kejadian banjir tanggal 3-4 Maret 2016
Sumber : Hasil analisis, 2017

## E. Hidrograf Collins
Penentuan  kejadian  banjir  telah  diperoleh dan  analisis  Φ indeks telah  diperoleh hidrograf  limpasan  langsung  serta  hujan efektif,  selanjutnya  dilakukan  analisis hidrograf  limpasan  langsung  dan  hujan efektif  dengan  menggunakan  metode Collins. Untuk bentuk

## Sumber : Hasil analisis, 2017
Gambar 10 Hidrograf Collins tanggal 23- 24 Februari 2016
Sumber : Hasil analisis, 2017
Gambar 11 Hidrograf Collins tanggal 2-3 Maret 2016

## Sumber : Hasil analisis, 2017
Gambar 12 Hidrograf Collins tanggal 3-4
Maret 2016
0,0 1,0 2,0 3,0 4,0 5,0 6,0 7,0 8,0 9,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 D eb it (m 3 /d et ) Waktu (jam) Kejadian Banjir I Tanggal 23-24 Maret 2016 Debit Banjir 0,0 1,0 2,0 3,0 4,0 5,0 6,0 7,0 8,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 D eb it (m 3 /d et ) Waktu (jam) Kejadian Banjir II Tanggal 2-3 Maret 2016 Debit Banjir 0,0 2,0 4,0 6,0 8,0 10,0 12,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 D eb it (m 3 /d et ) Waktu (jam) Kejadian Banjir III Tanggal 3-4 Maret 2016 Debit Banjir 0,0 1,0 2,0 3,0 4,0 5,0 6,0 7,0 8,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
D eb it (m 3 /d et) Waktu (jam) Kejadian Banjir IV Tanggal 5-6 Maret 2016 Debit Banjir 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 D eb it ( m 3 /d et ) Waktu (jam) Hidrograf Collins Tanggal 23-24 Februari 2016 Debit Banjir 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 D eb it (m 3 /d et ) Waktu (jam) Hidrograf Collins Tanggal 2-3 Maret 2016 Debit Banjir 0,0 0,2 0,4 0,6 0,8 1,0 1,2 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 D eb it (m 3 /d et ) Waktu (jam) Metode Collins Kejadiaan Banjir 3 - 4 Maret 2016 Debit Banjir

## Sumber : Hasil analisis, 2017
Gambar 13 Hidrograf Collins tanggal 5-6 Maret 2016

## F. Hidrograf Satuan Sintetik
Prosedur  pembuatan  hidrograf  satuan sintetik  yang  digunakan  adalah  HSS Gama  I,  HSS  Nakayasu,  HSS  Snyder, HSS Limantara dan HSS SCS.
1. HSS Gama I

## Sumber : Hasil analisis, 2017
Gambar 14 Grafik HSS Gama I Sub Sub DAS Pedindang

## 2. HSS Nakayasu
Sumber : Hasil analisis, 2017
Gambar 15 Grafik HSS Nakayasu Sub Sub DAS Pedindang

## 3. HSS Snyder

## Sumber : Hasil analisis, 2017
Gambar 16 Grafik HSS Snyder Sub Sub DAS Pedindang

## 4. HSS Limantara
Sumber : Hasil analisis, 2017
Gambar 17 Grafik HSS Limantara Sub Sub DAS Pedindang

## 5. HSS SCS

## Sumber : Hasil analisis, 2017
Gambar 18 Grafik HSS SCS Sub Sub

## DAS Pedindang
G. Kesesuaian Hidrograf Satuan

## Sintetik Dengan Terukur
Dari  analisis  hidrograf  satuan  Collins diperoleh  bentuk  hidrograf  satuan  terukur Sungai  Pedindang.  Hidrograf  satuan terukur  yang  diperoleh  dengan  metode hidrograf  collins  sebanyak  4  kejadian yaitu:  kejadian  banjir  23-24  Februari 2016,  kejadian  banjir  2-3  Maret  2016,
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 D eb it (m 3 /d et ) Waktu (jam) Metode Collins Kejadian Banjir Tanggal 5 -6 Maret 2016 Debit Banjir 0,0 0,5 1,0 1,5 2,0 2,5 3,0 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 D eb it (m 3/ de t) Waktu (jam) Hidrograf Satuan Sintetik Gama I Sub Sub DAS Pedindang Q Awa l Q Terkoreksi 0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 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 D eb it (m 3 /d et ) Waktu (jam) Hidrograf Satuan Sintetik Nakayasu Sub Sub DAS Pedindang Q Awal Q Terkoreksi 0,00 0,20 0,40 0,60 0,80 1,00 1,20 1,40 1,60 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 D eb it (m 3/ de t) Waktu (jam) Hidrograf Satuan Sintetik Snyder Sub Sub DAS Pedindang Q Awal Q Koreksi 0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 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 D eb it (m 3 /d et ) Waktu (jam) Hidrograf Satuan Sintetik Limantara Sub Sub DAS Pedindang HSS Limatara Awal HSS Limantara Terkoreksi 0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 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 D eb it ( m 3 /d et ) Waktu (jam) Hidrograf Satuan Sintetik Limantara Sub Sub DAS Pedindang HSS Limatara Awal HSS Limantara Terkoreksi
kejadian  banjir  3-4  Maret  2016  dan kejadian  banjir  5-6  Maret  2016.  Untuk bentuk  hidrograf satuan  sintetik  Sub Sub Daerah  Aliran  Sungai  diperoleh  dari analisis  hidrograf  satuan  sintetik  dengan menggunakan  5  prosedur  perhitungan yaitu:  HSS  Gama  I,  HSS  Nakayasu,  HSS Snyder, HSS Limantara dan HSS SCS.

## 1. Kejadian  Banjir  23-24  Februari
2016

## Sumber : Hasil analisis, 2017
Gambar 19 Grafik kesesuaian HSS Sub Sub DAS Pedindang Terhadap HST Sungai Pedindang 23-24 Februari 2016

## Sumber : Hasil analisis, 2017
Gambar 20 Grafik kesesuaian HSS Sub Sub DAS Pedindang Terhadap HST Sungai Pedindang 2-3 Maret 2016

## Sumber : Hasil analisis, 2017
Gambar 21 Grafik kesesuaian HSS Sub Sub DAS Pedindang Terhadap HST Sungai Pedindang 3-4 Maret 2016
Sumber : Hasil analisis, 2017
Gambar 22 Grafik kesesuaian HSS Sub Sub DAS Pedindang Terhadap HST Sungai Pedindang 5-6 Maret 2016

## H. Validasi  Kesesuaian  Hidrograf Metode RMSE
Metode  RMSE  membandingkan  debit estimasi  dengan  debit  observasi,  dimana nilai  RMSE  kurang  dari  (<10%),  maka model tersebut paling sesuai.
1. Kejadian  Banjir  23-24  Februari 2016
Tabel 5 Hasil nilai RMSE kejadian banjir I
Sumber: Hasil analisis, 2018
2. Kejadian Banjir 2-3 Maret 2016
Tabel 6 Hasil nilai RMSE kejadian banjir II

## Sumber: Hasil analisis, 2018
3. Kejadian Banjir 3-4 Maret 2016
Tabel 7 Hasil nilai RMSE kejadian banjir III
0,0 2,0 4,0 6,0 8,0 10,0 12,0 14,0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 D eb it (m 3 /s ) Waktu (jam) HS Terukur S. Pedindang HS Gama I S. Pedindang HSS Nakayasu HSS Limantara S. Pedindang HSS Snyder S. Pedindang HSS SCS S. Pedindang 0,00 2,00 4,00 6,00 8,00 10,00 12,00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 D eb it (m 3 /s) Waktu (jam) HS Terukur S. Pedindang HSS Gama I S. Pedindang HSS Nakayasu S. Pedindang HSS Limantara S. Pedindang HSS Snyder S. Pedindang HSS SCS S. Pedindang 0,00 2,00 4,00 6,00 8,00 10,00 12,00 14,00 16,00 18,00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 D eb it (m 3/ s) Waktu (jam) HS Terukur S. Pedindang HSS Gamma 1 S. Pedindang HSS Nakayasu S. Pedindang HSS Limantara S. Pedindang HSS SCS S. Pedindang HSS SCS S. Pedindang
0,00 2,00 4,00 6,00 8,00 10,00 12,00 14,00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 D eb it (m 3 /s ) Waktu (jam) HS Terukur S. Pedindang HSS Gamma 1 S. Pedindang HSS Nakayasu S. Pedindang HSS Limantara S. Pedindang HSS Snyder S. Pedindang HSS SCS S. Pedindang
Sumber: Hasil analisis, 2018

## 4. Kejadian Banjir 5-6 Maret 2016
Tabel 8 Hasil nilai RMSE kejadian banjir
IV
Sumber: Hasil analisis, 2018

## KESIMPULAN DAN SARAN
A. Kesimpulan
Berdasarkan hasil analisis dan pembahasan mengenai kesesuaian hidrograf satuan sintetik terhadap hidrograf  satuan  terukur  Sub  Sub  Daerah Aliran  Sungai  Pedindang  Bagian  Tengah maka  dapat  disimpulkan  beberapa  hal sebagai berikut:
1. Kesesuaian  hidrograf  satuan  terukur terhadap  hidrograf  satuan  sintetik Sungai Pedindang sangat berbeda satu dengan  lainnya,  hal  ini  dapat  dilihat pada hasil analisis. Pada hasil analisis untuk  setiap  kejadian  banjir,  debit puncak  hidrograf  satuan  sintetik Sungai  Pedindang  sangat  berbeda dengan debit puncak hidrograf satuan terukurnya.  Rata-rata  debit  puncak hidrograf  satuan  sintetik  terjadi  pada kisaran  2  atau  3  jam,  sedangkan hidrograf  satuan  terukur  Sungai Pedindang terjadi pada kisaran 7 atau 8 jam.
2. Dari  keempat  kejadian  banjir  Sungai Pedindang    tidak  diperoleh  nilai
RMSE yang memenuhi syarat kurang dari  10%  (<  10%),  namun  yang mendekati syarat hanya satu hidrograf satuan sintetik dengan nilai berikut: a. Tanggal  23-24  Februari  2016 sebesar  23,601%  adalah  Hidrograf Satuan Sintetik Gama I
b. Tanggal  2-3  Maret  2016  sebesar 16,315%  adalah  Hidrograf  Satuan Sintetik Gama I
c. Tanggal  3-4  Maret  2016  sebesar
50,400%  adalah  Hidrograf  Satuan Sintetik Gama I
d. Tanggal 5-6 Maret sebesar
22,322%  adalah  Hidrograf  Satuan Sintetik Gama I 3. Dengan  diperolehnya  nilai  RMSE (validasi) digunakan untuk menyatakan  model  hidrograf  satuan sintetik  yang  paling  sesuai  dengan hidrograf  satuan  terukur  Sungai Pedindang.  Maka,  dari  keempat kejadian  banjir  Sungai  Pedindang, tidak  ada  model  hidrograf  satuan sintetik  mempunyai  nilai  RMSE  < 10%  dan  menyatakan  bahwa  seluruh model  hidrograf  satuan  sintetik  tidak berkesesuaian terhadap hidrograf satuan terukur Sungai Pedindang.

## B. Saran
1. Dalam  penelitian  ini  pada  analisis kesesuaian  bentuk  hidrograf  satuannya sangat  berbeda  jauh.  Dari  hasil  analisis, debit  puncak  hirograf  satuan  sintetik berada  pada  kisaran  waktu  2  atau  3  jam. Sedangkan  pada  hidrograf  satuan  terukur, debit puncak hirograf berada pada kisaran 7  atau  8  jam.  Hasil  penelitian  ini dipengaruhi  oleh  penentuan  kejadian
banjirnya. Pada penentuan kejadian banjir, data  yang  digunakan  adalah  data  tinggi muka  air  dari  alat  AWLR  pada  tanggal- tanggal  yang  bukan  merupakan  puncak banjir  pada  Sungai  Pedindang.  Hasil analisis  juga  dipengaruhi  data  hujan  jam- jaman  yang  digunakan  terlalu  jauh  dari lokasi  penelitian,  dimana  data  hujan  yang digunakan  adalah  data  hujan  dari  Stasiun Cuaca  Universitas  Bangka  Belitung. Untuk  penelitian  selanjutnya,  dalam analisis hidrograf satuan terukur sebaiknya digunakan lebih banyak kejadian  banjir  sebagai  acuan  serta  data hujan  yang  digunakan  adalah  data  hujan terdekat pada lokasi penelitian.

## DAFTAR PUSTAKA
Asdak,  Chay.  2014.  Hidrologi  dan Pengelolaan  Daerah  Aliran  Sungai. Cetakan
Keenam. Yogyakarta:
Gadjah Mada University Press.
Kamiana, I Made. 2012. Teknik Perhitungan Debit Rencana Bangunan  Air.  Cetakan  Kedua. Yogyakarta: Graha Ilmu Soewarno. 2013. Hidrometri Dan Aplikasi Teknosabo Dalam Pengelolaan Sumber Daya Air. Cetakan Pertama. Yogyakarta: Graha Ilmu. Soewarno. 2015. Analisisn Data Hidrologi Menggunakan Metode Statitiska dan Stokastik. Cetakan Pertama. Yogyakarta: Graha Ilmu Sri  Harto.  2000.  Hidrologi:  Teori, Masalah, Penyelesaian. Cetakan Pertama. Yogyakarta: Beta Offset.
Susanti, Meilia Nur Indah. 2014. Statistika Deskriptif  &  Induktif.  Cetakan Kedua. Yogyakarta: Graha Ilmu.
Triatmodjo,  Bambang.  2014.  Hidrologi Terapan. Cetakan Keempat. Yogyakarta: Beta Offset.