Ϯϯ

## ANALISIS KERAPATAN JARINGAN STASIUN CURAH HUJAN PADA WILAYAH SUNGAI (WS) AESESA DI PULAU FLORES
Yerison Dimu Ratu (Jeri.dimu@yahoo.com) 1)
Denik Sri Krisnayanti 2)
I Made Udiana 3)

## ABSTRACT
The rain that falls on the surface of the earth is influenced by several factors that are not evenly spread to a region of the river. It can be determined by the proper placement of rain station, good location, number and pattern of spread. However, the placement of the rain stations are generally based only on immediate needs, so do not pay attention to the development of water resources as a whole.
This research is specifically done on Aesesa with area of 8.202,41 km2. Methods to be used in this study is the method of Kagan.
Analysis of the WS Aesesa obtained daily rainfall correlation coefficient 0,184 and 0,293 monthly. Existing density of rainfall stations in the WS Aesesa 1.025,0 km 2 /Station. For daily rainfall, averaging 5% error obtained an area 195,295 km 2 /stations and 42 stations, alignment errors of 10% obtained an area of 585.886 km 2 /stations and 14 stations. While 5% for monthly rainfall, obtained  an area of 182,276 km 2 / stations and 45 stations, 10% error obtained an area 546,827 km 2 /stations and 15 stations. The number of stations is not the only factor affecting the level of rainfall forecast accuracy, but the pattern of spread also plays a role in determining the level of accuracy.
Keywords : Number of Rainfall Station, Density, Kagan method

## ABSTRAK
Hujan yang jatuh dipermukaan bumi dipengaruhi oleh beberapa faktor sehingga penyebarannya tidak merata untuk suatu wilayah sungai. Hal tersebut dapat diketahui dengan penempatan stasiun hujan yang tepat, baik lokasi, jumlah dan pola penyebarannya. Namun penempatan stasiun hujan pada umumnya hanya didasarkan pada kebutuhan sesaat, sehingga belum memperhatikan pengembangan sumber daya air secara menyeluruh.
Penelitian ini dilakukan pada WS Aesesa dengan luas 8,202.41 km 2 . Metode yang akan digunakan dalam penelitian ini adalah metode kagan.
Analisis pada WS Aesesa didapatkan koefisien korelasi hujan harian 0,184 dan bulanan 0,293. Kerapatan stasiun hujan existing pada WS Aesesa 1.025,0 km 2 / Stasiun. Untuk hujan harian, kesalahan perataan 5% didapat luasan 195,295 km 2 /stasiun dan 42 stasiun, kesalahan perataan 10% didapat luasan 585,886 km 2 /stasiun dan 14 stasiun. Sedangkan 5% untuk hujan bulanan, didapatkan luasan 182,276 km 2 /stasiun dengan 45 stasiun, kesalahan 10% didapat luasan 546,827 km 2 /stasiun dan 15 stasiunn. Jumlah stasiun bukan satu-satunya faktor
yang berpengaruh terhadap tingkat ketelitian perkiraan hujan, pola penyebaran juga berperan dalam menentukan tingkat ketelitian hitungan Kata Kunci : Jumlah Stasiun Curah Hujan, Kerapatan, Metode Kagan

## 1. PENDAHULUAN
Secara umum analisis hidrologi merupakan satu bagian analisis awal dalam perancangan bangunan-bangunan hidraulik. Data hidrologi adalah kumpulan keterangan atau fakta mengenai fenomena hidrologi (Soewarno 1995, dalam Joksan). Salah satu data hidrologi yang penting dalam analisis hidrologi adalah data curah hujan. Data curah hujan didapat dari pengukuran pada stasiun hujan. Karena intensitas, penyebaran, serta kedalaman hujan berbeda dan tidak merata disetiap wilayah, maka pola penempatan dan penyebaran stasiun pencatatan curah hujan harus tepat sehingga diharapkan dapat memberikan data yang mewakili lokasi dimana stasiun tersebut berada.

## 2. TINJAUAN PUSTAKA

## A. Pengertian Siklus Hidrologi
Secara keseluruhan jumlah air di planet bumi relatif tetap dari masa ke masa. Air di bumi mengalami suatu siklus melalui serangkaian peristiwa yang berlangsung terus-menerus, dimana kita tidak tahu kapan dan dari mana berawalnya dan kapan pula akan berakhir. Serangkaian  peristiwa tersebut disebut  atau dinamakan siklus hidrologi (Suripin, 2003: 20).

## B. Curah Hujan

## 1. Cara Mengukur Curah Hujan
Besarnya curah hujan diukur dengan menggunakan alat penakar curah hujan. Alat penakar curah hujan dibedakan menjadi dua grup, yakni Alat penakar hujan manual dan Alat penakar hujan otomatis, (Suripin, 2003 : 24).

## 2. Cara Perhitungan Curah Hujan Daerah
Curah hujan yang diperlukan untuk menyusun suatu rancangan pemanfaatan air adalah curah hujan rata-rata di daerah yang bersangkutan, bukan hanya pada satu
titik tertentu. Curah hujan ini disebut curah hujan wilayah atau daerah dan dinyatakan dalam mm.
Cara-cara perhitungan curah hujan daerah dari pengamatan curah hujan di beberapa pos stasiun hujan adalah sebagai berikut :
a) Cara Rata rata Aljabar
b) Cara Poligon Thiessen
c) Metode Isohyet
C. Jaringan Pengukuran Hujan
Dalam merencanakan jaringan stasiun hujan  terdapat dua hal penting yang harus diperhatikan, yaitu:
1. Berapa jumlah stasiun hujan, dinyatakan dalam Km 2 / stasiun. Kerapatan jaringan dinyatakan dalam satu stasiun tiap luas tertentu, misalnya 1 stasiun setiap 200 km 2 .
2. Pola penempatan stasiun dalam wilayah sungai atau dimana stasiun-
stasiun tersebut akan dipasang.
Hal tersebut diperlukan, karena dalam jaringan stasiun hujan, perbedaan jumlah stasiun yang digunakan dalam memperkirakan besar hujan yang terjadi dalam suatu WS memberi perbedaan dalam besaran hujan yang didapat. Selain itu pola penyebaran  stasiun hujan dalam WS yang bersangkutan juga ternyata mempunyai pengaruh yang nyata terhadap ketelitian hitungan hujan rata-rata WS.
Pada dasarnya cara Kagan mempergunakan analisis statistik dan mengaitkan kerapatan jaringan pengukur hujan dengan kesalahan interpolasi dan kesalahan perataan (interpolation error and averaging error). Persamaan-persamaan yang dipergunakan :
r (d) =r (0) e –d /d(0) ……………………………. (1)
N N d A r Cv Z ) 0 ( ) 0 ( 1 23 , 0 1 + − = …………...………(2) ( ) A/N /d 0,52r /3 r 1 Cv Z (0) (0) (0) 3 + − = …...……..…(3)
A/N 1,07 l = ……………………....…..(4)
Dimana :
r(d) =  koefisien korelasi untuk jarak d km
r(0) =  koefisien korelasi untuk jarak sangat dekat
d =  jarak antar stasiun, dalam km
d(0)  =  radius korelasi, yaitu     jarak    antar stasiun dimana korelasi berkurang dengan faktor e.
Z1  =  kesalahan perataan, dalam %
Cv  =  koefisien variasi
A  =  luas Daerah Aliran Sungai (DAS), dalam km2
N  =  jumlah stasiun
Z3  =  kesalahan interpolasi, dalam %
l =  jarak antar  stasiun
Dari hubungan antara jarak antar stasiun dan koefisien korelasi (r), dapat dapat digambarkan grafik lengkung eksponensial, seperti yang nampak pada gambar 1 sebagai berikut.
Gambar 1.  Korelasi antar stasiun hujan pada Suatu WS Sumber : Sri Harto Br, 2000
Berdasarkan persamaan (2) dan (3), dapat diperoleh grafik hubungan antara jumlah stasiun dan ketelitian yang diperoleh baik untuk hujan harian maupun hujan bulanan sebagai berikut.

## Gambar 2. Kesalahan sebagai fungsi jumlah Stasiun pada suatu WS
Sumber : Sri Harto Br, 2000
Koefisien variasi (variation coefficient) adalah nilai perbandingan antara deviasi standar dengan nilai rata-rata hitung dari suatu distribusi. Koefisien variasi dapat dihitung dengan rumus sebagai berikut. (Chay Asdak, 2004 : 296)
x S Cv = ……………………………….(5)
Dengan :
Cv =  koefisien variasi S =  deviasi standar
x =  rata-rata hitung

## 3. METODE PENELITIAN
Analisis yang biasanya digunakan dalam pengelolaan data hidrologi adalah analisis statistik.
Secara garis besar langkah-langkah yang ditempuh sebagai berikut ini :
1. Menentukan stasiun – stasiun hujan yang aktif pada WS Aesesa
2. Menentukan curah hujan maksimum harian dan curah hujan maksimum bulanan pada stasiun – stasiun hujan yang aktif pada WS Aesesa
3. Menghitung jarak antar stasiun hujan
4. Menghitung korelasi antar stasiun curah hujan, baik untuk hujan harian maupun hujan bulanan, sesuai dengan yang diperlukan.
5. Hubungan yang diperoleh di atas digambarkan dalam sebuah grafik lengkung eksponensial, dari grafik ini dapat diperoleh besaran d (0) dengan menggunakan nilai rata-rata d dan r (d) .
6. Dengan besaran tersebut, maka kesalahan perataan dan kesalahan interpolasi dapat dihitung dengan persamaan (2) dan (3), setelah tinggi ketelitian ditetapkan.
7. Setelah jumlah stasiun ditetapkan untuk Wilayah Sungai tersebut maka penetapan stasiun hujan dapat dilakukan dengan menggunakan persamaan (4) dan menggambarkan jaring-jaring segitiga sama sisi dengan panjang sisi sama dengan satu

## 4. HASIL DAN PEMBAHASAN
Wilayah Sungai (WS) Aesesa terletak Pulau Flores. Luas Wilayah Sungainya adalah 8.202,4 km 2 yang meliputi empat Kabupaten, yaitu Kabupaten Ngada, Kabupaten Nagekeo, Kabupaten Ende dan Kabupaten Sikka.
WS Aesesa tergolong daerah yang beriklim tropis yang dipengaruhi musim barat yang membawa hujan dari laut Jawa dan laut Flores serta angin musim panas dari Benua Australia. Kerapatan stasiun hujan yang ada sekarang (eksisting) pada WS Aesesa rata-rata adalah 1.025,00 km 2 /stasiun.
Berdasarkan  studi literatur, diketahui bahwa pada WS Aesesa terdapat 22 Stasiun Curah Hujan MRG (Manual Rain fall Gauge), 7 diantaranya tidak berfungsi karena mengalami kerusakan. Jadi stasiun yang masih aktif berjumlah 15 buah. Dari 15 buah stasiun yang aktif tersebut, 7 diantaranya baru di bangun.
Dengan demikian ketujuh stasiun tersebut dianggap tidak aktif (data yang dipakai adalah data dari stasiun dengan data di atas 10 tahun). Stasiun curah hujan yang masih aktif akan dipakai dalam analisis selanjutnya yaitu sebanyak 8 stasiun.
Sedangkan jarak  antar  stasiun  curah hujan di ukur  dari satu stasiun ke stasiun yang lain dan sebaliknya dengan menggunakan Software Map Source, sehingga didapat jarak semua stasiun dalam WS Aesesa, seperti pada tabel 1.
Tabel 1. Jarak Antar Stasiun Hujan Pada WS Aesesa
Analisis korelasi diperlukan guna menentukan besarnya korelasi stasiun hujan yang satu dengan stasiun hujan yang lain dan untuk memperkirakan besarnya kesalahan yang terjadi baik untuk curah hujan harian maupun bulanan. Perhitungan analisis korelasi antar stasiun untuk curah hujan harian dan curah hujan bulanan ditunjukan pada tabel 2, sebagai berikut :
Tabel 2. Jarak Rerata  Dan Korelasi masing-masing Stasiun Hujan Pada WS Aesesa
Hubungan antara jarak antar stasiun dan koefisien korelasi untuk hujan harian maupun hujan bulanan untuk semua stasiun pada WS Aesesa seperti pada Gambar 3 berikut:
Bajawa Waepana Riung Mbay Nanganio Maumere Ledalero Welamosa 1 Bajawa 0,00 18,92 36,12 37,84 60,20 138,26 130,72 87,72 2 Waepana 18,92 0,00 24,08 19,78 47,30 126,42 118,68 74,82 3 Riung 36,12 24,08 0,00 27,52 52,89 132,87 127,28 77,83 4 Mbay 37,84 19,78 27,52 0,00 30,96 108,36 101,48 57,19 5 Nanganio 60,20 47,30 52,89 30,96 0,00 78,26 71,81 27,95 6 Maumere 138,26 126,42 132,87 108,36 78,26 0,00 7,31 50,74 7 Ledalero 130,72 118,68 127,28 101,48 71,81 7,31 0,00 47,30 8 Welamosa 87,72 74,82 77,83 57,19 27,95 50,74 47,30 0,00 91,746 86,369 60,507 No Stasiun Jarak Terhadap Stasiun (Km) Jarak rerata 72,826 61,429 68,370 54,733 52,767
1 2 3 4 5 6 7
8 0,237 0,528 Nanganio 52,767 0,616 0,590 Mbay 54,733 86,369
60,507 68,593
0,167 Ledalero Welamosa 0,155 Riung 68,370 0,110 0,392 0,424 0,367 0,255 No. Stasiun Hujan
Jarak Rerata Terhadap Waepana 61,429 r hari r bulan Stasiun Lainnya Bajawa 72,826 0,315 0,483 0,257 0,404 0,179 Maumere 91,746 0,281 0,404
Gambar 3. Grafik Korelasi Antar Stasiun Hujan Pada WS Aesesa
Berdasarkan Gambar 3 di atas, dengan persamaan  1 : ( ) ( ) ( )       − = 0 d d e o r d r maka akan diperoleh  besaran-besaran sebagai berikut:
Hujan harian
r (0) =  0,453
d (0) = 50,251 Hujan bulanan r (0) =  0,529 d (0) = 57,803
Perhitungan curah hujan rata-rata menggunakan cara garis isohyet dengan interval tinggi hujan harian 50 mm. Perhitungan koefisien variasi hujan harian seperti pada Tabel 3 dan 4 di bawah ini.
ƌ;ĚͿсϬ͘ϰϱϯĞ ͲϬ͘ϬϭϵϵĚ ZϸсϬϯϰ r(d) = 0529e -0.0173d R² = 0.086 Ϭ͘Ϭ Ϭ͘Ϯ Ϭ͘ϰ Ϭ͘ϲ Ϭ͘ϴ ϱϬ͘Ϭ ϲϬ͘Ϭ ϳϬ͘Ϭ ϴϬ͘Ϭ ϵϬ͘Ϭ ;< ŽĞ ĨŝƐ ŝĞ Ŷ <Ž ƌĞ ůĂ ƐŝͿ Jarak (Km)

## ϭͿ DĂŚĂƐŝƐǁĂdĞŬŶŝŬ^ŝƉŝůhŶŝǀĞƌƐŝƚĂƐEƵƐĂĞŶĚĂŶĂ Ϯ͕ϯͿ ŽƐĞŶdĞŬŶŝŬ^ŝƉŝůhŶŝǀĞƌƐŝƚĂƐEƵƐĂĞŶĚĂŶĂ
ϯϭ
Tabel 3  Perhitungan Koefisien Variasi Curah Hujan Harian (1) (10) (12) 1 100 - 150 -2 4.00 2 150 - 200 -1 1.00 3 200 - 250 1 1.00 4 250 - 300 2 4.00 No. Curah Hujan Luas Luas Titik Tengah Frekuensi Xi.fi (2) (3) (4) (5) (6) (mm) (km) (%) Xi fi │Xi -    │ fi.│Xi-   │ Ci 1534.23 18.70 125 1534.23 191778.8 75 115067.3 2178.39 2178.39 -3068.46 6136.92 2456.34 29.95 175 2456.34 429859.5 25 61408.5 -2456.34 2456.34 559198.8 75 2178.39 26.56 225 2178.39 490137.8 25 54459.75 720.49 Jumlah 8202.41 100.00 800 8202.41 383444.3 2033.45 24.79 275 2033.45 18905.45 152508.8 4066.9 8133.8 (7) fi.Ci 2 Ci.fi Ci 2 (8) (9) (11) (13) __ X
__ X

## Tabel 4 Perhitungan Koefisien Variasi Curah Hujan Bulanan
(1) (10) (12) 1 375 - 525 -2 4.00 2 525 - 675 -1 1.00 3 675 - 825 1 1.00 4 825 - 975 2 4.00 No. Curah Hujan Luas Luas Titik Tengah Frekuensi (2) (3) (4) (5) (6) (7) = 5*6 fi.Ci 2 (mm) (km) (%) Xi fi Xi.fi │Xi -  │ fi.│Xi-  │ Ci Ci.fi Ci 2 (8) (9) (11) (13) 1812.33 22.10 450 1812.33 815548.5 225 407774.3 -3624.66 7249.32 1778.69 21.68 600 1778.69 1067214 75 133401.8 -1778.69 1778.69 2644.35 32.24 750 2644.35 1983263 75.0 198326.3 2644.35 2644.35 1967.04 23.98 900 1967.04 1770336 225.0 442584 3934.08 7868.16 Jumlah 8202.41 100.00 2700 8202.410 1182086 19540.52 1175.08 X X
ϯϮ
Menghitung Curah Hujan Rata-Rata Curah hujan rata-rata (      )   = 4 0 . 800 =  200. mm

## 1. Menghitung Deviasi Rata-Rata
Deviasi rata-rata =                             = =   46.74 mm/ hari 2. Menghitung Standar Deviasi Rata-Rata ( S )
Standar deviasi rata-rata  (S) =   75.782 mm/ hari
3. Menghitung Koefisien Variasi Koefisien Variasi ( Cv )  = __ S X = 75.782 200.0 =   0.379 Jadi koefisien variasi curah hujan harian adalah 0.379 Menghitung Curah Hujan Rata-Rata Curah hujan rata-rata (      )   = 4 0 . 800 =  200. mm
1. Menghitung Deviasi Rata-Rata 2
50
18905.45 - 720.49
8202.410 8202.410 __ X ∑ ∑ = = − n i i i n i i f X X f 1 1 383444.3 8202.41 2 1 1 1 1 2 . .             − ∑ ∑ ∑ ∑ = = = = n i i n i i i n i i n i i i f C f f C f i __ X ∑ ∑ = = − n i i i n i i f X X f 1 1 383444.3 8202.41
ϯϯ
Deviasi rata-rata =                             =   46.74 mm/ hari
2. Menghitung Standar Deviasi Rata-Rata ( S ) Standar deviasi rata-rata  (S) =   75.782 mm/ hari
3. Menghitung Koefisien Variasi Koefisien Variasi ( Cv )  = __ S X = 75.782 200.0 =   0.379
Jadi koefisien variasi curah hujan harian adalah 0.379
Perhitungan curah hujan rata-rata menggunakan cara garis isohyet dengan interval tinggi hujan bulanan 150 mm. Perhitungan koefisien variasi hujan bulanan seperti pada tabel 4
1. Menghitung Curah Hujan Rata-Rata
Curah hujan rata-rata (      )   = 4 0 . 2700 =  675 mm/bulan 2. Menghitung Deviasi Rata-Rata Deviasi rata-rata = = 1182086 8202.41 =  144 __ X ∑ ∑ = = − n i i i n i i f X X f 1 1 2 50 18905.45 - 720.49
8202.410 8202.410 2 1 1 1 1 2 . .             − ∑ ∑ ∑ ∑ = = = = n i i n i i i n i i n i i i f C f f C f i

## 3. Menghitung Standar Deviasi Rata-Rata ( S )
Standar deviasi rata-rata  (S)  = = =   230.521 mm/ bulan
4. Menghitung Koefisien Variasi Koefisien Variasi ( Cv )  = __ S X = 230.521 675.0 =   0.342
Jadi koefisien variasi curah hujan bulanan adalah 0.342
Dengan  persamaan (2) dan (3) dihitung besarnya kesalahan perataan (Z1) dan kesalahan interpolasi (Z3) baik untuk curah hujan harian maupun bulanan yang hasilnya dapat dilihat pada gambar 4 dan 5 sebagai berikut :
Gambar 4. Kesalahan perataan dan interpolasi untuk curah hujan harian pada WS Aesesa 0 10 20 30 40 50 0 10 20 30 40 50 60 70 80 90 <Ğ ƐĂ ůĂ ŚĂ Ŷ ;й Ϳ :ƵŵůĂŚ^ƚĂƐŝƵŶ;ŶͿ keralahan perataan kesalahan interpolasi 2 1 1 1 1 2 . .             − ∑ ∑ ∑ ∑ = = = = n i i n i i i n i i n i i i f C f f C f i 2 150 8202.410 19540.52 - 1175.08 8202.410
Gambar 5. Kesalahan perataan dan interpolasi untuk curah hujan Bulanan pada WS Aesesa
Berdasar gambar tersebut, atau persamaan (2) dapat ditetapkan jumlah stasiun hujan yang diperlukan apabila besar kesalahan yang diperlukan dapat ditetapkan. Untuk kesalahan 5% dan 10%, jaringan stasiun hujan harian memerlukan masing-masing  42 dan 14 stasiun hujan. Apabila didasarkan pada hujan bulanan diperoleh masing-masing 45dan 15 stasiun. Hasil analisis data dengan Cara Kagan menunjukkan bahwa jumlah stasiun yang ada dalam WS Aesesa sekarang belum memenuhi atau masih lebih kecil dibandingkan dengan jumlah stasiun yang dituntut dengan Cara Kagan, Sehingga perlu penetapan jumlah stasiun hujan. Untuk kesalahan 5 % pada hujan harian memiliki kerapatan 195,295 km2/stasiun dan kesalahan 10 % memiliki kerapatan 585,886 km2/stasiun, sedangkan kesalahan 5 % pada hujan bulanan memiliki kerapatan 182,276 km2/stasiun, dan kesalahan 10 % untuk hujan bulanan memiliki kerapatan 546,827 km2/stasiun.
Ϭ ϭϬ ϮϬ ϯϬ ϰϬ Ϭ ϭϬ ϮϬ ϯϬ ϰϬ ϱϬ ϲϬ ϳϬ ϴϬ ϵϬ <Ğ ƐĂ ůĂ ŚĂ Ŷ ;й Ϳ :ƵŵůĂŚ^ƚĂƐŝƵŶ;ŶͿ <ĞƐĂůĂŚĂŶ ƉĞƌĂƚĂĂŶ <ĞƐĂůĂŚĂŶ /ŶƚĞƌƉŽůĂƐŝ
Gambar 6. Peta WS Aesesa untuk Curah Hujan Harian dengan Kesalahan 5%

## 5. KESIMPULAN
Berdasarkan hasil analisis dan pengolahan data untuk evaluasi kerapatan jaringan stasiun hujan Wilayah Sungai (WS) Aesesa di Pulau Flores didapatkan bahwa :
1. Koefisien korelasi rata-rata hujan harian adalah 0.184 dan hujan bulanan
0.293. Nilai  koefisien korelasi ini menunjukan korelasi yang lemah antar dua variabel.
2. Besarnya kesalahan perataan untuk 8 buah stasiun yang sudah ada pada
WS Aesesa adalah 12.925 % untuk curah hujan harian dan 10.108 % untuk curah hujan bulanan.
3. Hasil analisis dengan Cara Kagan memberikan besarnya kesalahan untuk hujan harian dan bulanan. Untuk curah hujan harian dengan kesalahan 5%, didapat jumlah stasiun (n) = 42, dengan jarak antar stasiun (l) = 14.953 km dan untuk kesalahan 10%, didapat jumlah stasiun (n) = 14, dengan jarak antar stasiun (l) = 25.899 km. Sedangkan untuk hujan bulanan kesalahan 5%, didapat  jumlah stasiun (n) = 45, dengan jarak antar stasiun (l) = 14.446 km dan kesalahan 10% didapat jumlah stasiun (n) = 15 dengan jarak antar stasiun (l) = 25.021 km.
4. Rekomendasi stasiun baru pada daerah sabana untuk tingkat kesalahan 10
% yakni 6 stasiun untuk curah hujan harian dan 7 stasiun untuk curah hujan bulanan.
Gambar 7. Peta WS Aesesa untuk Curah Hujan Harian dengan Kesalahan 5%

## DAFTAR PUSTAKA
Anonim.  2009.   Profil  Balai  Wilayah Sungai  Nusa  Tenggara  II, BWS  NT  II, Nusa Tenggara Timur.
Chay  Asdak.  2004.  Hidrologi dan Pengelolaan Daerah Aliran Sungai, Gajah Mada University Press,
Yogyakarta
Soemarto CD. 1987. Hidrologi Teknik, Erlangga, Jakarta
Soewarno. 1995. Hidrologi Jilid 1, Penerbit Nova , Bandung
Sri Harto Br. 1993. Analisis Hidrologi, PT. Gramedia Pustaka Utama, Jakarta
Sri Harto Br. 2000.  Hidrologi, Teori Masalah Penyelesaian,  Nafiri Offset,
Yogyakarta
Suyono Sosrodarsono, Kensaku Takeda. 2003. Hidrologi Untuk Pengairan, PT. Pradnya Paramita, Jakarta.
Suripin. 2004. Sistem Drainase Perkotaan Yang Berkelanjutan, Andi, ogyakarta.