## ANALISIS PENGARUH PERUBAHAN PENGGUNAAN LAHAN TERHADAP DEBIT LIMPASAN PADA DAERAH ALIRAN SUNGAI BONDOYUDO KABUPATEN LUMAJANG DENGAN METODE RASIONAL
Akbar Kurniawan 1 , Fransedo Aminata 2
1,2,3 Departemen Teknik Geomatika, FTSLK-ITS, Kampus ITS Sukolilo, Surabaya, 60111, Indonesia e-mail: 1 akbar.geodesy@gmail.com, 2 fransedoaminata10@gmail.com

## Abstrak
Banjir  merupakan  peristiwa  meluapnya  air  sungai  dalam  jumlah  yang  besar  pada  saat  musim  hujan.  Wilayah  DAS Bondoyudo adalah salah satu wilayah yang rawan banjir di Kabupaten Lumajang. Penelitian ini menggunakan metode rasional untuk menghitung debit limpasan di setiap wilayah penelitian dengan rumus sebagai berikut, Q=0,278.C.I.A. Dari hasil penelitian ini, wilayah yang tergenang banjir berdasarkan metode rasional yaitu Kecamatan Rowokangkung, Jatiroto,  dan  Kedungjajang.  Pada  tahun  2002  wilayah  Kecamatan  Rowokangkung  dan  Kedungjajang  memiliki  selisih debit limpasan yang berlebih sebesar 3,65 m3/s dan 2,23 m3/s. Sementara pada tahun 2013 terjadi peningkatan selisih debit menjadi 3,72 m3/s dan 2,42 m3/s. Kecamatan Jatiroto juga berpotensi banjir pada tahun 2013 dengan selisih debit 0,35  m3/s.  Dan  pada  tahun  2018  selisih  debit  pada  masing-masing  Kecamatan  Rowokangkung,  Kedungjajang,  dan Jatiroto adalah 3,74 m3/s ; 2,54 m3/s; dan 0,39 m3/s. Besar pengaruh perubahan penggunaan lahan terhadap debit limpasan  dihitung  berdasarkan  rumus  persamaan  regresi  linier  ganda  dan  dihasilkan  dari  perhitungan  koefisien determinasi.  Nilai  koefisien  determinasi  dari  Kecamatan  Rowokangkung,  Kedungjajang,  dan  Jatiroto  adalah  92,5%  , 93,1% , dan 92,4%.
Kata Kunci: Banjir, Metode Rasional, Debit Limpasan, Regresi Linier Ganda.

## Abtract
Floods is an overflowing of a large amount of water on rainy seasons. Bondoyudo watershed is one of the most flood‐prone area in Lumajang district. This research use rational method for calculate the surface runoff in each research areas. The formula is : Q = 0,278 CIA. From this research, based on rational method, flooded areas are Rowokangkung sub‐district, Jatiroto sub‐district, and Kedungjajang sub‐district. In 2002, Rowokangkung sub‐ district and Kedungjajang sub‐district had excess runoff discharge as 3,65 m3/s and 2,23 m3/s. In 2010,  Rowokangkung, Kedungjajang, and Jatiroto sub‐district had excess runoff discharge as 3,72 m3/s, 2,42 m3/s, and 0,35 m3/s. In 2018, Rowokangkung, Kedungjajang, and Jatiroto sub‐district had excess runoff discharge as 3,74 m3/s, 2,54 m3/s, 0,39 m3/s. The value of the effect of land use change on runoff discharge is calculated  based on  multiple  linear regression and generated from the calculation of the coefficient of determination. The  value  of  determination  coefficient  for  Rowokangkung  Sub‐ distric, Kedungjajang Sub‐district, and Jatiroto Sub‐ district are 92.5%, 93.1%, and 92.4%.
Keywords: Flood, Rational Method, Runoff Discharge, Multiple Linear Regression

## PENDAHULUAN
Kabupaten  Lumajang  merupakan  salah  satu wilayah di Jawa Timur yang sering dilanda bencana banjir.  Hal  ini  karena  kondisi  hidrologi  wilayah Kabupaten  Lumajang  sendiri  yang  merupakan bagian  hilir  dari  6  Daerah  Aliran  Sungai  (DAS), antara  lain:  DAS  Bondoyudo,  DAS  Glidik,  DAS Mujur, DAS Rawaan, serta DAS Rejali. Banjir yang terjadi  rata-rata  disebabkan  oleh  meningkatnya debit  air  pada  saat  terjadi  hujan  secara  terus- menerus sehingga mengakibatkan air meluap dan terjadi  banjir  genangan  di  lingkungan  sekitar daerah aliran sungai. Selain disebabkan oleh curah hujan,  banjir  genangan  juga  disebabkan  oleh perubahan  penggunaan  lahan  di  daerah  aliran sungai. Kondisi ini dapat menyebabkan terjadinya gangguan  terhadap  sistem  hidrologi.  Perubahan kondisi  lahan  dari  waktu  ke  waktu  membuat ancaman terjadinya banjir semakin besar. Hal ini disebabkan oleh beberapa hal, antara lain:
(1) Daya  tampung  sungai  semakin  lama  semakin kecil akibat pendangkalan.
(2) Fluktuasi  debit  air  antara  musim  penghujan dengan musim kering semakin tinggi.
(3) Terjadi  konversi  lahan  pertanian  dan  daerah buffer  alami  ke  lahan  non  pertanian  dengan mengabaikan
konservasi sehingga
menyebabkan rusaknya daerah tangkapan air (cacthment  area).  Perubahan  penggunaan lahan dari lahan non-terbangun menjadi lahan terbangun  akan  menyebabkan  terjadinya peningkatan  aliran  permukaan  dan  semakin sedikitnya air tanah (Butler, 2011).
Bencana banjir genangan dapat diprediksi dengan  pendekatan  metode  rasional  (Suripin, 2003).  Metode  Rasional  merupakan  metode matematis  untuk  menghitung  debit  limpasan permukaan  di  wilayah  sekitar  aliran  sungai. Parameter-parameter  yang  diperlukan  dalam perhitungan  metode  rasional  ini  adalah  curah hujan, kemiringan lereng, jenis tanah, dan tutupan lahan.
Pada penelitian sebelumnya yang dilakukan  oleh  Ishmata  (2016)  yaitu  membahas tentang  aplikasi  dari  citra  tandem-X  untuk  studi potensi  genangan  dengan  metode  rasional.  Hasil dari  penelitian  tersebut  yaitu  wilayah  sub-DAS mana  yang  berpotensi  terjadi  banjir  genangan berdasarkan perhitungan metode rasional. Untuk
penelitian  kali  ini  penulis  menggunakan  citra satelit
Landsat
Multitemporal untuk mengidentifikasi perubahan penggunaan lahan di setiap  tahun  penelitian.  Analisis  yang  dilakukan yaitu  menghitung  debit  limpasan  di  setiap kecamatan dan menghitung selisih terhadap debit sungai.  Selain  menggunakan  metode  rasional, penelitian  kali  ini  juga  menggunakan  metode skoring  atau  pembobotan  untuk  mengetahui daerah mana saja yang rawan terhadap bencana banjir  di  DAS  Bondoyudo  Kabupaten  Lumajang. Hasil  peta  potensi  banjir  dengan  metode  skoring kemudian dibandingkan hasilnya dengan perhitungan  debit  limpasan  metode  rasional  di setiap wilayah kecamatan.

## METODE
Lokasi penelitian ini berada di wilayah Kabupaten Lumajang,  Provinsi  Jawa  Timur  yang  secara geografis  terletak  pada  koordinat  antara 112°50’00”  s/d  113°22’00”  BT  dan  7°52’00”  s/d 8°23’00”  LS  khususnya  wilayah  DAS  Bondoyudo. DAS Bondoyudo mencakup 12 wilayah Kecamatan yaitu Kecamatan Ranuyoso, Klakah, Kedungjajang, Randuagung, Jatiroto, Rowokangkung, Yosowilangun,  Tekung,  Lumajang,  Sukodono, Gucialit, dan Senduro. Lokasi penelitian bisa dilihat pada gambar 1 dibawah ini.

## Gambar 1. Lokasi Penelitian
Untuk menganalisis daerah potensi banjir genangan  dengan  metode  rasional  pada  DAS Bondoyudo diperlukan data-data antara lain Citra Landsat 7 Tahun 2002, Citra Landsat 8 Tahun 2013 dan  2018,  Data  Curah  Hujan  Tahun  2002-2018, Data  Jenis  Tanah,  Data  DEMNAS  Kabupaten Lumajang,  dan  Data  Debit  Sungai  Eksisting  DAS Bondoyudo.

## Diagram Alir Pengolahan Data Metode Rasional

## Gambar 2. Diagram Alir Pengolahan Data Metode Rasional
Proses  pengolahan  data  menggunakan  metode rasional:
1. Penentuan  Luas  Wilayah  Kecamatan  Data spasial  Kabupaten  Lumajang  berupa  batas kecamatan  yang  termasuk  dalam  DAS Bondoyudo. Kemudian dilakukan perhitungan luas area pada data spasial batas kecamatan, dimana  luas  wilayah  tiap  kecamatan  ini merupakan  salah  satu  parameter  dalam perhitungan debit limpasan menurut metode rasional.
2. Pengolahan Peta Tutupan Lahan Peta tutupan lahan  didapatkan  dari  hasil  klasifikasi terbimbing  citra  Landsat  7  Kabupaten Lumajang  Tahun  2002.  Dan  citra  Landsat  8 Kabupaten Lumajang tahun 2013 dan 2018.
3. Pengolahan Peta Kelerengan Peta Kelerengan diturunkan  dari  data  DEMNAS  Kabupaten Lumajang.  Dari  data  DEMNAS  tersebut kemudian dibuat peta kontur untuk selanjutanya  baru  diturunkan  menjadi  peta kelerengan.  Peta  kelerengan  diklasifikasikan sesuai karakteristik DAS menjadi empat kelas lereng, yaitu 0%-8%, 8%15%, 15%-25%, 25%- 45% dan >45%.
4. Perhitungan Intensitas Hujan Dalam perhitungan intensitas curah hujan, diperlukan beberapa tahap dalam perhitungannya.  Tahap  awal  dilakukan  untuk menghitung  curah  hujan  tiap  kecamatan. Selanjutnya adalah menghitung waktu konsentrasi  aliran  (Tc)  dengan  parameter ketinggian  yang  didapat  dari  data  DEMNAS. Dan  tahap  akhir  yang  dilakukan  adalah menghitung  intensitas  curah  hujan  dengan rumus Mononobe,
I = [((R24/24) x (24/Tc))] 2/3 (1)
5. Penentuan  Nilai  Koefisien  Limpasan  (C)  Nilai koefisien  limpasan  didapatkan  dari  nilai  hasil overlay peta kelerengan, peta tutupan lahan,
dan peta jenis tanah dengan rumus C=Ct + Cs + Cv.                             (2)
Harga C memiliki rentang antara 0 – 1. 6. Perhitungan  Debit  Limpasan  Perhitungan debit  maksimum  dilakukan  di  setiap  wilayah kecamatan  yang  termasuk  dalam  DAS Bondoyudo  dengan  menggunakan  Metode Rasional dengan 3 parameter utamanya yaitu nilai  koefisien  limpasan,  intensitas  curah hujan, dan luas area tiap kecamatan.
7. Perhitungan  Selisih  Debit  Analisis  statistik dilakukan  dengan  melakukan  perhitungan selisih  antara  data  debit  eksisting  dan  nilai debit  limpasan maksimum  hasil  perhitungan. Apabila  nilai  selisih  debit  bernilai  negative, maka  debit  limpasan  lebih  besar  daripada debit  sungai,  dan  dapat  dipastikan  daerah tersebut  berpotensi  banjir  genangan  ketika musim hujan.
8. Analisis  Pengaruh  Perubahan  Penggunaan Lahan  Terhadap  Debit  Limpasan  Analisis dilakukan  dengan  menghitung  luas  tutupan lahan  tahun  2002,  2013,  2018  dan dikorelasikan dengan debit sungai tahun 2002, 2013, 2018 dengan menggunakan persamaan regresi  linier  ganda. Dari  perhitungan  regresi linier ganda akan didapatkan nilai korelasi (R) dan  nilai  koefisien  determinasi  (R2).  Dimana koefisien  determinasi  merupakan  besaran pengaruh  antara  perubahan  penggunaan lahan terhadap debit limpasan.

## Diagram Alir Pengolahan Data Metode Skoring

## Gambar 3. Diagram Alir Pengolahan Data Metode Skoring
Analisis  Potensi  Banjir  menggunakan  metode skoring  (pembobotan)  didasarkan  pada  5 parameter penyebab bencana banjir yaitu tutupan lahan,  buffer  sungai,  kemiringan  lereng,  jenis tanah,  dan  curah  hujan.  Proses  ini  dilakukan dengan  cara  memberikan  nilai  pada  setiap parameter  penyebab  banjir.  Kemudian,  masing- masing parameter tersebut ditampalkan (overlay)
sehingga  dihasilkan  data  spasial  baru  berupa daerah  rawan  banjir.  Analisa  daerah  potensi bencana  didasarkan  pada  nilai  total  skor  dari masing-masing area. Area yang memiliki total skor atau nilai paling tinggi merupakan area yang rawan bencana  banjir.  Penentuan  bobot  dan  skor beracuan pada contoh tabel skoring yang terdapat pada  tinjauan  pustaka  namun  dengan  modifikasi untuk  penyesuaian  dengan  wilayah  Kabupaten Lumajang. Berikut ini adalah Tabel yang memuat skor  dan  bobot  dari  masing-masing  parameter penyebab banjir.

## HASIL DAN PEMBAHASAN
Gambar 4. Peta Tutupan Lahan Tahun 2002
Peta  Penggunaan  Lahan  Penggunaan  lahan merupakan salah satu parameter dalam perhitungan  potensi  genangan  menurut  metode rasional.  Hasil  klasifikasi  penggunaan  lahan Kabupaten  Lumajang  dibagi  menjadi  beberapa kelas penggunaan lahan antara lain permukiman, perkebunan, sawah, vegetasi, dan hutan. Gambar 4, 5, 6 serta Tabel 2, 3, 4 merupakan peta dan tabel perhitungan  luas  tiap  kelas  tutupan  lahan  tahun

## 2002, 2013, dan 2018.
Gambar 5. Peta Tutupan Lahan Tahun 2013
Gambar 6. Peta Tutupan Lahan Tahun 2018

## Peta Jaringan Sungai
Data  jaringan  sungai  diperoleh  dari  Dinas  PU Sumber Daya Air Lumajang. Dari hasil pengolahan dan  analisa  dapat  diketahui  wilayah  sungai  yang sering  meluap  dan  mengakibatkan  banjir genangan saat musim hujan yaitu Sungai Jatiroto yang melewati Kecamatan Jatiroto, Rowokangkung,  dan  Randuagung.  Sungai  yang kedua  yaitu  Sungai  Grobogan  yang  melewati Kecamatan  Kedungjajang  dan  Klakah.  Berikut merupakan peta sungai yang sering terjadi banjir di DAS Bondoyudo.

## Gambar 7. Peta Jaringan Sungai

## Peta Jenis Tanah
Jenis tanah merupakan salah satu faktor penyebab banjir genangan. Dimana jenis tanah berpengaruh terhadap proses infiltrasi air hujan ke dalam tanah. Semakin halus tekstur tanah maka tanah semakin resisten  dan  tidak  menyimpan  air  melainkan mengalirkannya  sehingga  menjadi  limpasan permukaan dan peka terhadap erosi.

## Gambar 8. Peta Jenis Tanah
Tanah  dengan  tekstur  berlapis  batu- batuan  memiliki  tekstur  tanah  yang  keras, sehingga  air  hujan  yang  jatuh  pada  wilayah tersebut  akan  langsung  menjadi  aliran  limpasan permukaan.  Akan  tetapi,  jenis  tanah  seperti  ini biasanya  ditemukan  di  daerah  dataran  tinggi sehingga  air  hujan  tidak  akan  terjebak  dan menimbulkan genangan.
Jenis tanah yang
berpotensi
menyebabkan  banjir  genangan  adalah  tanah dengan  tekstur  lanau  yang  memiliki  koefisien sebesar 0,16.  Tekstur tanah  lanau  sulit  ditembus oleh  air  karena  memiliki  ukuran  butiran  yang sangat  kecil  dan  halus.  Tanah  bertekstur  lanau banyak terdapat pada Kecamatan Rowokangkung, Kecamatan Jatiroto, dan Kecamatan Yosowilangun,  sehingga  kecamatan-kecamatan tersebut berpotensi terjadi banjir genangan.

## Peta Kelerengan
Kondisi  Topografi  atau  kelerengan  suatu  daerah mempengaruhi  kecepatan  dan  volume  limpasan permukaan  ketika  terjadi  hujan  di  wilayah tersebut. Semakin curam tingkat kelerengan maka kecepatan  aliran  permukaan  semakin  besar. Semakin datar tingkat kelerengan maka kecepatan aliran  permukaan  juga  semakin  kecil  (Ishmata, 2016).  Hal  ini  berarti  bahwa  air  hujan  akan semakin  lama  terjebak  pada  wilayah  datar tersebut dan kemungkinan terjadi banjir genangan semakin besar.
Kemiringan  lereng  pada  DAS  Bondoyudo Kabupaten  Lumajang  didominasi  oleh  wilayah yang  datar,  yaitu  pada  wilayah  kecamatan Rowokangkung,  Jatiroto,  Yosowilangun,  Tekung, Lumajang, Sukodono, dan kedungjajang, sehingga wilayah-wilayah  tersebut  memiliki  potensi  banjir
genangan yang lebih tinggi dibandingkan dengan wilayah lain.
Gambar 9. Peta Kelerengan

## Nilai Koefisien Limpasan
Koefisien  limpasan  adalah  persentase  jumlah  air yang dapat melimpas/mengalir melalui permukaan tanah dari keseluruhan air hujan yang jatuh  pada  suatu  daerah  (Suripin,  2003).  Nilai koefisien  limpasan  ditentukan  berdasarkan  3 parameter  yaitu  tutupan  lahan,  jenis  tanah,  dan kemiringan lereng. Dengan menggunakan rumus C =  Ct  +  Cs  +  Cv  dihasilkan nilai  koefisien  limpasan untuk tiap kecamatan seperti Tabel 7.
Dari  3  tahun  penelitian  tersebut, Kecamatan  Jatiroto  memiliki  koefisien  limpasan yang paling besar, diikuti Kecamatan Rowokangkung pada urutan kedua. Dimana pada tahun 2002 Kecamatan Jatiroto memiliki koefisien limpasan  0,530,  pada  tahun  2013  naik  menjadi 0,552,  dan  pada  tahun  2018  juga  naik  menjadi 0,565.
Sementara untuk Kecamatan
Rowokangkung  pada  tahun  2002,  2013,  2018 secara berturut-turut memiliki koefisien limpasan sebesar 0,533, 0,536, dan 0,531.
Dua  kecamatan  tersebut  berada  pada dataran  rendah  dan  jenis  tanah  dengan  tekstur lanau serta lempung yang memiliki nilai koefisien limpasan  yang  cukup  besar.  Sehingga  dua
kecamatan  tersebut  memiliki  potensi  banjir genangan yang besar juga.

## Intensitas Hujan
Intensitas  hujan  didefinisikan  sebagai  tren  curah hujan  per  satuan  waktu.  Dalam  hal  ini,  perlu diketahui  curah  hujan  yang  terjadi  dalam  satu tahun. Perhitungan dilakukan dalam kurun waktu 10  tahun  yaitu  antara  tahun  2009-2018.  Untuk mendapatkan  intensitas    hujan  selama    waktu konsentrasi, digunakan rumus Mononobe (Loebis, 1992).

## Debit Limpasan
Debit limpasan dihitung dari tiga parameter utama yaitu koefisien limpasan, intensitas hujan, dan luas wilayah  penelitian.  Hasil  perhitungan  debit maksimum  pada  setiap  wilayah  kecamatan disajikan pada tabel berikut:

## Debit Sungai
Debit Sungai eksisting merupakan volume air yang mengalir  per  satuan  waktu  pada  saat  musim hujan.  Debit  sungai  eksisting  bisa  menjadi  acuan kapasitas sungai dalam menampung air pada saat musim  hujan.  Berikut  merupakan  tabel  debit sungai di masing-masing kecamatan.

## Selisih Debit Limpasan dan Debit Sungai
Setelah  dilakukan  proses  perhitungan  debit limpasan maksimum di setiap wilayah kecamatan, maka selanjutnya dihitung selisih nilainya dengan data  debit  eksisting  tiap  sungai.  Nilai  debit eksisting  dibandingkan  dengan  nilai  debit limpasan yang lewat maka dapat diketahui daerah mana  saja  yang  aman  dan  daerah  mana  yang banjir  akibat  tidak  mampu  menampung  debit limpasan  yang  lewat  pada  saat  terjadi  hujan maksimum pada waktu tertentu. Untuk nilai selisih bertanda  negatif  mempunyai  arti  bahwa  daerah
tersebut berpotensi banjir karena debit limpasan yang  lebih  besar  daripada  debit  eksisting  sungai. Berikut  merupakan  hasil  dari  perhitungan  selisih antara debit limpasan dan debit sungai.

## Skoring (Pembobotan)
Potensi Banjir Skoring merupakan pemberian nilai terhadap  suatu  polygon  peta  untuk  memberikan tingkat  kedekatan,  keterkaitan  atau  beratnya dampak  tertentu  pada  suatu  fenomena  secara spasial  (Pratomo,  2009).  Peta  potensi  banjir dengan  metode  skoring  digunakan  sebagai
pembanding  hasil  perhitungan  banjir  genangan dengan  metode  rasional.  Tabel  12  memuat  skor dan  bobot  dari  masing-masing  parameter penyebab banjir.
Peta potensi rawan banjir didapatkan dari hasil overlay peta curah hujan, peta tutupan lahan, peta  kemiringan  lereng,  peta  jenis  tanah,  dan buffer  sungai.  Potensi  banjir  dalam  penelitian  ini dibagi  menjadi  3  kelas  tingkat  potensi  yaitu rendah,  sedang,  dan  tinggi.Berikut  merupakan peta potensi banjir dengan metode skoring tahun 2002, 2013, 2018.
Gambar 10. Perbandingan Hasil Metode Skoring dan Metode Rasional Tahun 2002
Gambar 11. Perbandingan Hasil Metode Skoring dan Metode Rasional Tahun 2013 Gambar 12. Perbandingan Hasil Metode Skoring dan Metode Rasional Tahun 2018

## Korelasi Perubahan Penggunaan Lahan Terhadap Debit Limpasan
Analisis  Korelasi  digunakan  untuk  mengetahui pengaruh  perubahan  luas  penggunaan  lahan terhadap jumlah debit limpasan. Dalam penelitian ini  digunakan  analisis  regresi  linier  berganda dengan  5  variabel  bebas  dan  1  variabel  terikat. Variabel bebas yang dimaksud yaitu kelas tutupan lahan  antaraa  lain  hutan  (X1),  kebun  (X2), permukiman  (X3), sawah  (X4),  dan vegetasi  (X5). Sedangkan  untuk  variabel  terikatnya  yaitu  debit limpasan hasil perhitungan (Y). Bentuk persamaan dari  regresi  linier  berganda  dengan  5  variabel bebas seperti berikut (Chapra, 1996):
Y = b 0 + b 1 .X 1 + b 2 .X 2 + b 3 .X 3 + b 4 .X 4 + b 5 .X 5
Nilai koefisien regresi dari masing-masing variabel  ditampilkan  pada  Tabel  16.  Dengan menggunakan  rumus  persamaan  regresi  linier berganda  maka  berdasarkan  nilai  koefisien  pada tabel  diatas  untuk  Kecamatan  Rowokangkung persamaannya adalah:
Y  =  44,895 – 0,692X1 – 0,102X2  +  0,082X3  + 0,110X4 – 0,098X5
Hal  ini  mengandung  pengertian  bahwa setiap ada pertambahan luas hutan (X1) sebesar 1 km2  maka  jumlah  debit  limpasan  (Y)  akan berkurang  sebesar  0,692  dengan  asumsi  bahwa variabel bebas yang lain dari model regresi adalah tetap.  Hal  yang  sama  juga  terjadi  apabila  luas kebun (X2)  bertambah 1 km2 maka jumlah debit limpasan berkurang sebesar 0,102. Untuk variabel permukiman  (X3)  apabila  terjadi  pertambahan luas 1 km2 maka jumlah debit limpasan juga ikut bertambah sebesar 0,082 karena bertanda positif. Untuk  variabel  sawah  (X4)  apabila  terjadi pertambahan  luas  1  km2  maka  jumlah  debit limpasan  bertambah  0,110.  Dan  untuk  variabel vegetasi  (X5)  apabila  terjadi  pertambahan  luas  1 km2 maka jumlah debit limpasan akan berkurang sebesar 0,098.
Selanjutnya dilakukan perhitungan koefisien  korelasi  (R)  dan  koefisien  determinasi (R2) untuk mengetahui hubungan korelasi antara perubahan  luas  penggunaan  lahan  dan  debit limpasan.  Nilai  korelasi  dihitung  berdasarkan rumus berikut:
Hasil  perhitungan  koefisien  korelasi  dari setiap  wilayah  kecamatan  penelitian  disajikan pada  Tabel  17.  Dilihat  dari  Tabel  17,  rata-rata setiap  kecamatan  memiliki  hubungan  korelasi antara  0,9  s/d  1  dimana  nilai  tersebut  masuk  ke dalam kategori korelasi kuat. Jadi dapat diartikan bahwa  perubahan  luas  ke-5  kelas  tutupan  lahan berpengaruh signifikan terhadap perubahan debit limpasan setiap tahunnya.
Untuk mengetahui besaran pengaruh dari perubahan  luas  tutupan  lahan  terhadap  jumlah debit  limpasan  dapat  dihitung  dengan  rumus koefisien  determinasi,  yang  disajikan  pada  tabel dibawah ini:
Nilai  koefisien  determinasi  menunjukkan besaran  pengaruh  dari  variabel  tutupan  lahan terhadap  jumlah  debit  limpasan.  Sebagai  contoh pada Kecamatan yang berpotensi banjir menurut metode rasional yaitu Kecamatan Rowokangkung dan  Jatiroto  memiliki  koefisien  determinasi sebesar 0,935 dan 0,924. Dengan dikalikan 100% maka  pengaruh  variabel  tutupan  lahan  pada Kecamatan Rowokangkung dan Jatiroto berpengaruh  93,5%  dan  92,4%  terhadap  jumlah debit limpasan. Sedangkan 6-7% dipengaruhi oleh faktor  lain  seperti  curah  hujan,  jenis  tanah,  dan kemiringan lereng.
Faktor  lain  yang  mempengaruhi  debit limpasan  yang  pertama  yaitu  curah  hujan. Menurut  Sudaryono  (2002)  hujan  berpengaruh terhadap  limpasan  yang  meliputi  tipe,  lama, intensitas dan sebaran hujan sangat menentukan
limpasan permukaan yang terjadi di suatu daerah aliran  sungai  (DAS).  Jumlah/volume  debit limpasan berkaitan dengan intensitas dan lamanya hujan  turun  yang  terjadi  pada  DAS  yang bersangkutan.
Faktor  jenis  tanah  juga  berkontribusi terhadap peningkatan debit limpasan. Jenis tanah yang  terdapat  pada  Kecamatan  Rowokangkung dan  Jatiroto  merupakan  tanah  dengan  tekstur lanau.  Dimana  tanah  bertekstur  lanau  memiliki koefisien yang cukup tinggi yaitu 0,16. Tanah lanau memiliki  tekstur  yang  sangat  halus  sehingga memiliki  daya  infiltrasi  yang  rendah, maka  tanah semakin  resisten  dan  tidak  meyimpan  air melainkan  mengalirkannya  sehingga  menjadi limpasan  air  dan  peka  terhadap  erosi  (Arsyad, 2010).
Bentuk topografi suatu daerah
(kelerengan)  juga berpengaruh pada volume  dan debit limpasan permukaan. Kecamatan Rowokangkung  dan  Jatiroto  berada  pada  tingkat kemiringan  lereng  yang  datar  (08%).  Dengan kondisi seperti ini maka air hujan yang turun lebih banyak terjebak dan akan menjadi genangan.

## KESIMPULAN
Kesimpulan dari hasil penelitian ini adalah:
1. Wilayah  yang  berpotensi  banjir  berdasarkan perhitungan  metode  rasional  pada  tahun 2002  adalah  Kecamatan  Rowokangkung  dan Kecamatan  Kedungjajang.  Sedangkan  pada tahun  2013  wilayah  yang  berpotensi  banjir bertambah  satu  dari  tahun  2002  yaitu Kecamatan  Jatiroto.  Dan  pada  tahun  2018 wilayah  yang  berpotensi  banjir  yaitu Kecamatan Rowokangkung, Kecamatan Kedungjajang, dan Kecamatan Jatiroto.
2. Pada  tahun  2002  wilayah  Kecamatan Rowokangkung dan Kecamatan Kedungjajang memiliki  kelebihan  debit  limpasan  dengan selisih  masing-masing  3,65  m3/s  dan  2,23 m3/s.  Sementara  pada  tahun  2013  terjadi peningkatan selisih debit pada dua kecamatan tersebut menjadi 3,72 m3/s  dan
2,42 m3/s. Kecamatan Jatiroto juga berpotensi  banjir  pada  tahun  2013  dengan selisih debit 0,35 m3/s. Dan pada tahun 2018 selisih debit pada masing-masing Kecamatan Rowokangkung,  Kecamatan  Kedungjajang,
dan  Kecamatan  Jatiroto  adalah  3,74  m3/s  ; 2,54 m3/s; dan 0,39 m3/s.
3. Besar  pengaruh  perubahan  tutupan  lahan terhadap  peningkatan  debit  limpasan  untuk Kecamatan  Rowokangkung  adalah  93,5%. Kecamatan Kedungjajang sebesar 93,1%. Dan untuk  Kecamatan  Jatiroto  memiliki  memiliki pengaruh sebesar 92,4%.

## UCAPAN TERIMA KASIH
Penulis mengucapkan terima kasih kepada Bapak Akbar  Kurniawan  ST,  MT.  selaku  dosen pembimbing  dalam  penelitian  ini.  Terima  kasih atas  kesempatan,  bimbingan,  serta  dukungan hingga dapat terselesaikannya penelitian ini. Dan ucapan  terima  kasih  kepada  UPT  Dinas  PU Sumber Daya Air Kabupaten Lumajang yang telah memberikan dukungan berupa data curah hujan dan  debit  sungai  sehingga  penelitian  ini  dapat berjalan dengan lancar.

## DAFTAR PUSTAKA
Arsyad, S. 2010. Konserversi Tanah dan Air. Bogor: IPB
Press.
Butler, D. dan Davies, J. W. 2011. Urban Drainage, Third Edition. New York: Taylor and Francis Group.
Chapra,  Steven.  1996.  Metode  Numerik  Jilid  1 (Diterjemahkan  oleh  I  Nyoman  Susila).  Jakarta: Erlangga.
Ishmata, A. 2016. Penggunaan Citra TANDEM-X Untuk Studi  Potensi  Genangan  Dengan  Prakiraan  Debit Sungai  Menggunakan  Metode  Rasional.  Skripsi. Teknik  Geomatika,  Institut  Teknologi  Sepuluh Nopember
Loebis, Joerson. 1992. Banjir Rencana Untuk Bangunan Air. Departemen Pekerjaan Umum.
Pratomo.  2009.  Analisis  Kerentanan  Banjir  di  Daerah Aliran  Sungai  Sengkarang  Kabupaten  Pekalongan Provinsi  Jawa  Tengah  dengan  Bantuan  Sistem Informasi  Geografis.  Fakultas  Geografi  Universitas Muhammadiyah. Surakarta
Primayuda  A,  2006.  Pemetaan  Daerah  Bahaya  dan Resiko  Banjir  Menggunakan  Sistem  Informasi Geografis: studi kasus Kabupaten Trenggalek, Jawa Timur.  Tugas  Akhir.  Bogor:  Fakultas  Pertanian, Institut Pertanian Bogor.
Sudaryono.  2002.  Pengelolaan  Daerah  Aliran  Sungai (DAS) Terpadu.Bogor: Angkasa.
Suripin.  2003.  Sistem  Drainase  Perkotaan  yang Berkelanjutan. Yogyakarta : Andi.